Yongdong Peng
Mingzhu Zhu
Yunyi Gong
Changfa Wang*- Liaocheng Research Institute of Donkey High-Efficiency Breeding and Ecological Feeding, Agricultural Science and Engineering School, Liaocheng University, Liaocheng, China
Many studies have shown that long non-coding RNAs (lncRNAs) play key regulatory roles in various biological processes. However, the importance and molecular regulatory mechanisms of lncRNAs in donkey intramuscular fat deposition remain to be further investigated. In this study, we used published transcriptomic data from the longissimus dorsi muscle of Guangling donkeys to identify lncRNAs and obtained 196 novel lncRNAs. Compared with the coding genes, the novel lncRNAs and the known lncRNAs exhibited some typical features, such as shorter transcript length and smaller exons. A total of 272 coding genes and 52 lncRNAs were differentially expressed between the longissimus dorsi muscles of the low-fat and high-fat groups. The differentially expressed genes were found to be involved in various biological processes related to lipid metabolism. The potential target genes of differentially expressed lncRNAs were predicted by cis and trans. Functional analysis of lncRNA targets showed that some lncRNAs may act on potential target genes involved in lipid metabolism processes and regulate lipid deposition in the longissimus dorsi muscle. This study provides valuable information for further investigation of the molecular mechanisms of lipid deposition traits in donkeys, which may improve meat traits and facilitate the selection process of donkeys in future breeding.
1 Introduction
Long non-coding RNA (lncRNA) is a class of non-coding RNAs with a length greater than 200 nucleotides. Increasing evidence suggests that lncRNA plays important roles in various biological processes, such as embryonic development (1, 2), gene expression regulation (3, 4), reprogramming (5, 6), and genomic imprinting (7, 8). Additionally, many lncRNAs have been found to be involved in regulating lipid metabolism. For example, lncRNA IMFNCR promotes chicken myoblast differentiation by sequestering miR-128-3p and miR-27b-3p (9), lncRNA 332,443 inhibits preadipocyte differentiation by targeting Runx1, p38 MAPK, and ERK1/2-MAPK signaling pathways (10), and LncLSTR forms a molecular complex with TDP-43 to regulate the expression level of Cyp8b1, thereby affecting the FXR regulatory pathway, leading to increased apoC2 levels and influencing triglyceride levels. LNCLSR directly binds to TDP-43 to inhibit Cyp8b1 expression and subsequently regulate triglyceride levels (11).
Guangling donkeys are distributed in Guangling County, Shanxi Province, China, and they are a local dominant breed that is carefully reared by local people using traditional production practices (12). Guangling donkeys have a stout physique and full muscles. Guangling donkeys used for meat production have a high intramuscular fat (IMF) content; however, the underlying molecular mechanisms underlying the IMF variation among donkey species are not fully understood.
IMF, also known as marbling, is an important indicator of the lean meat-to-fat ratio, which directly affects the tenderness, juiciness, and flavor of the meat (13). IMF content is one of the most important indicators used to evaluate meat quality (14). IMF content is a polygenic trait that is regulated by many genes affecting adipogenesis and lipid metabolism (15). At present, the underlying molecular variations affecting IMF content among donkey breeds are unclear.
This study used published transcriptome data from the longissimus dorsi muscle of a Guangling donkey to identify lncRNA and conducted differential expression studies of coding genes and lncRNA, constructing an expression regulation network, which lays the foundation for further analyzing the molecular mechanisms of lipid deposition traits in donkeys.
2 Materials and methods
2.1 Data sources
A total of 30 Guangling donkeys were raised on a commercial donkey farm in Fanshi County, Xinzhou City, Shanxi Province, China, and 6 donkeys with IMF differences and similar ages were selected (age: 2–3 years old, weight: 232–245 kg; female) for use in this study. All Guangling donkeys were reared under the same natural conditions of uncontrolled room temperature and light with unrestricted access to food and water. The longissimus dorsi samples at the 13th rib were aseptically and quickly obtained within 30 min of harvest. The collected samples were stored in liquid nitrogen for immediate storage, and long-term storage was carried out at −80°C. According to the China National Standard GB5009.6-2016 “Determination of Fat in Foods in National Food Safety Standard,” the IMF content was determined by the Soxhlet method. A Soxhlet extraction apparatus was used to remove fat and dry the ground meat samples for fat extraction. Petroleum ether was used as a solvent. This was recycled and dried for 8 h. Then, it was weighed to obtain the weight of the bottle containing fat. The IMF content was calculated by a formula. The three longissimus dorsi samples with the highest IMF contents and the other three with the lowest IMF contents were selected for transcriptome analysis (12). A total of six RNA-seq datasets were obtained from a previously published study and downloaded from the NCBI’s GEO database (PRJNA658642) (12). The donkey gene annotations were downloaded from https://ftp.ensembl.org/pub/release-110/gtf/equus_asinus. Moreover, the Non-Redundant Protein Sequence (NR) Database was downloaded from ftp://ftp.ncbi.nih.gov/blast/db/. The uniref90 database was downloaded from https://ftp.uniprot.org/pub/databases/uniprot/uniref/uniref90.
2.2 RNA-seq reads mapping and transcriptome assembly
The quality of sequencing reads was evaluated by FastQC command. The raw reads were filtered and trimmed by Trimmomatic (version 0.39) with default parameters (16). The clean reads were then mapped to the donkey reference genome [Ensembl: ASM1607732v2 (GCA_016077325.2)] by HISAT2 v2.2.1 with the default parameters (17–20). StringTie (version 2.2.1) was used to assemble the mapped reads with default parameters (17). Then, the merge tool of StringTie was used to merge the six assembled transcript files (GTF format) of the two groups into a non-redundant transcriptome. In addition, by using the assembled GTF file, StringTie software was used to estimate the expression levels of genes and transcripts in all samples for subsequent studies with the parameters “-e” and “-B” (17, 18, 21).
2.3 LincRNAs identification pipeline
The pipeline for lincRNA (long intergenic non-coding RNA) identification was as follows (Figure 1): (1) retained those transcripts with “u” category categorized by using gffcompare, which indicated intergenic transcripts (20, 21). (2) According to the merged GTF file, the transcripts with single exons and less than 200 bp in length were removed (18, 21). (3) The CPC2, CNCI, PLEK and LGC were used to assess the protein-coding potential of complete transcript sequences, and the transcripts that cannot encode proteins based on protein-coding potential were retained (22). (4) The HMMER was used to identify the transcripts translated in all six possible frames with homologs that were concluded in any of the known protein family domains in the Pfam database, and transcripts that matched to the Pfam hit (E-value < 1e-5) were excluded (18–21, 23). (5) BLASTX program (24) was used to filter out any transcripts that have similarities to known proteins in the NCBI NR and UniRef90 databases (E-value < 1e-5) (20, 21). (6) Reserve transcripts with FPKM values greater than 0 in at least one sample (18–20).
Figure 1. Pipeline for the identification of putative lincRNAs in this study. The frames in the direction of the arrow show the filtering process and the number of screened transcripts. “u”: Unknown, intergenic transcript; CPC2: coding potential calculator 2; nr: non-redundant protein sequence database; FPKM: fragments per kilobase of transcript per million mapped reads.
2.4 Comparisons between lncRNAs and protein-coding transcripts
We selected the transcripts annotated as “protein-coding” in the gene annotation file, and the obtained lncRNAs were screened with “known” and “novel” by “blastn” command. The transcript length, exon length, and exon number of lncRNAs were compared with those of protein-coding transcripts (20).
2.5 Analysis of differentially expressed genes and differentially expressed lncRNAs
DESeq2 tool was used to perform differential expression analysis of protein-coding genes and lncRNAs between the high (H) fat group and low (L) fat group (25). |log2 fold change | ≥ 1 and adjusted p-value (padj) < 0.05 were used to screen differentially expressed genes (DEGs) and lncRNAs (26).
2.6 Prediction of potential target genes
We predicted the molecular functions of protein-coding genes regulated by RNA in cis and trans. First, the neighboring protein-coding genes near DELs (<100 kb) were identified based on cis-prediction principles using Bedtools (18–20, 27, 28). For the trans-regulation of DELs, we calculated the Pearson’s correlation coefficient (r) between DELs and protein-coding genes. We selected protein-coding genes with a Pearson’s correlation coefficient |r| ≥ 0.95, p-value ≤ 0.01 as potential target genes (PTGs) of DELs (21, 29).
2.7 Functional enrichment analysis
Gene Ontology (GO) enrichment analysis was performed by clusterProfiler (30, 31). KOBAS v3.01 was used for the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis (32). A p-value of less than 0.05 was considered statistically significant (33).
3 Results
3.1 Summary of RNA-seq data mapping and transcripts assembly in longissimus dorsi muscles
RNA-seq data involving two groups of Guangling donkeys were obtained from a previously published study (12). The six longissimus dorsi muscles (three in each group) with the lowest and highest IMF contents were named the L group and H group, respectively. The individual samples in the groups were named L1, L2, L3, H1, H2, and H3. The clean reads were mapped to the donkey reference genome using HISAT2. Approximately 86.02–90.51% of clean reads from each library were mapped to the donkey reference genome, and 79.76–84.78% of the reads were uniquely mapped to the genome. Then, the transcriptome was assembled for each library by StringTie, and all transcripts were synthesized into non-redundant transcripts using StringTie-Merge. After merging non-redundant transcripts, approximately 1.07% (944 of 87,886) of the transcripts were intergenic transcripts. The 196 putative lincRNAs were obtained according to the illustration shown in Figure 1 (Supplementary Table S1).
3.2 Comparison of coding genes and lncRNA features
Previous studies showed that there are many differences between protein-coding transcripts and lncRNAs (18–20, 27, 34–36). According to the assembled transcriptome, the characteristics of lncRNA and protein-coding transcripts were compared. A total of 51,112 protein-coding transcripts, corresponding to 20,553 protein-coding genes annotated in donkeys, were acquired. In addition, the donkey annotation file contains 7,615 known lncRNA transcripts that correspond to 4,709 lncRNA genes (Supplementary Table S1).
The average transcript length of the protein-coding transcripts (2,919 bp) was longer than the novel lncRNA transcripts (1,217 bp) and the known lncRNA transcripts (2,394 bp) (Figure 2A). In terms of average exon length, the novel lncRNA is 448 bp in length, shorter than the known lncRNA gene (766 bp) but longer than the protein-coding transcripts (271 bp) (Figure 2B). In addition, we found that the average number of exons in the protein-coding transcripts is 10.7, which is significantly higher than that of the novel lncRNA (2.7) and the known lncRNA (3.1) (Figure 2C).
Figure 2. Comparison of the characteristics of protein-coding genes and lncRNA genes. (A) Comparison of transcript length, (B) comparison of exon length, and (C) comparison of exon number.
3.3 Differential expression analysis of coding genes and lncRNAs
To explore their potential biological functions, we performed differential expression analysis of coding genes and lncRNAs. A total of 272 DEGs were obtained by comparing the high-fat group with the low-fat group of the longissimus dorsi muscle samples. Among them, 147 genes were upregulated in the high-fat group, and 125 genes were downregulated (Figure 3A). Additionally, 52 differentially expressed lncRNAs (DELs) were identified, with 25 lncRNAs upregulated and 27 lncRNAs downregulated in the high-fat group compared to the low-fat group (Figure 3B).
Figure 3. Volcano plots in analyzing differentially expressed genes (DEGs) (A) and differentially expressed lncRNAs (DELs) (B) in longissimus dorsi muscles with high and low intramuscular fat contents. The green plot represents upregulated expression in the high group; the blue plot represents downregulated expression in the high group; and the gray plot represents no significance.
3.4 Functional analysis of differentially expressed genes
To investigate the function of DEGs, we performed GO and KEGG analyses, respectively. GO analysis of DEGs showed that integrin-mediated signaling pathway, cellular response to lipopolysaccharide, trachea development, cellular response to molecule of bacterial origin, and positive regulation of myeloid leukocyte-mediated immunity were the most abundant terms in the biological process category. In terms of cellular component category, actin cytoskeleton, phagocytic cup, cell leading edge, TORC2 complex, and basolateral plasma membrane were the top five terms, while rRNA binding, pattern recognition receptor activity, motor activity, protein kinase binding, and ion channel binding were most prevalent in the molecular function (Supplementary Table S2). Some GO terms were significantly associated with lipid metabolisms, such as response to lipids, cellular response to lipids, glycerophospholipid metabolic process, glycerolipid metabolic process, unsaturated fatty acid biosynthetic process, and phospholipid dephosphorylation. KEGG analysis indicated that DEGs were significantly enriched in 62 KEGG pathways, of which several pathways were related to lipid metabolism, such as the Sphingolipid signaling pathway and MAPK signaling pathway. In addition, some other pathways are related to lipid metabolism, namely fatty acid metabolism, glycerophospholipid metabolism, biosynthesis of unsaturated fatty acids, PI3K-Akt signaling pathway, fatty acid degradation, ether lipid metabolism, cholesterol metabolism, PPAR signaling pathway, TGF-beta signaling pathway, and Wnt signaling pathway (Supplementary Table S3).
3.5 Prediction and functional analysis of lncRNA target genes
Many studies have indicated that lncRNA may regulate adjacent genes in a cis manner (18–20, 29, 37–41). For PTGs regulated by lncRNAs in cis (<100 kb), we identified a total of 323 PTGs, corresponding to 52 DELs (Supplementary Table S4). To explore the function of putative lncRNAs, GO and KEGG analyses were performed on expressed protein-coding genes transcribed near lncRNA (<100 kb). The results indicated that 102 of 323 PTGs were significantly involved in 143 biological processes (Supplementary Table S5). Of them, 59 PTGs significantly participated in 32 pathways (Supplementary Table S6), including 2 pathways related to lipid metabolism, such as biosynthesis of unsaturated fatty acids and fatty acid metabolism.
LncRNAs can not only regulate the expression of neighboring protein-coding genes through a cis mechanism but also regulate the expression of genes located on other chromosomes via a trans mechanism (42). In this study, we carried out the trans analysis to find the PTGs that were significantly correlated (|r| ≥ 0.95, p ≤ 0.01) to the DELs. In total, 3,366 PTGs were highly correlated with 52 DELs. Among these genes, 132 PTGs were differentially expressed in groups as DEPTGs, suggesting that most of the lncRNAs regulated gene expression through trans regulation. GO enrichment analysis showed that 3,366 PTGs were enriched in 835 biological processes and 132 DEPTGs were enriched in 143 biological processes (Supplementary Tables S7, S8). In cases of biological process, some GO terms were significantly associated with lipid metabolism, such as negative regulation of the lipid catabolic process, regulation of the lipid metabolic process, negative regulation of the lipid metabolic process, regulation of the lipid biosynthetic process, cellular lipid catabolic process, negative regulation of the lipid biosynthetic process, regulation of fatty acid transport and regulation of the lipid catabolic process (Figures 4A,C). In addition, 3,366 PTGs and 132 DEPTGs were enriched in 74 pathways and 7 pathways, respectively (Supplementary Tables S9, S10). KEGG pathways were involved in the MAPK signaling pathway, PI3K-Akt signaling pathway, glycerolipid metabolism, ether lipid metabolism, fat digestion, and absorption (Figures 4B,D). The results indicated that DELs had an important role in regulating their PTGs that regulate lipid metabolism in longissimus dorsi muscles.
Figure 4. Gene Ontology and pathway analysis of PTGs of DELs. (A) Biological processes analysis of PTGs of DELs, (B) pathway analysis of PTGs of DELs, (C) biological processes analysis of DEPTGs of DELs, and (D) pathway analysis of DEPTGs of DELs. PTGs, potential target genes; DELs, differentially expressed lncRNAs; DEPTGs, differentially expressed potential target genes.
4 Discussion
In the present study, transcriptome sequencing of longissimus dorsi tissues with different IMF contents from Guangling donkeys was used to investigate genes and lncRNA related to lipid metabolism in the longissimus dorsi muscle. A total of 20,553 protein-coding genes, 7,615 known lncRNAs, and 196 novel lincRNAs were obtained. We compared the known and novel lincRNAs with the donkey protein-coding genes and found that known and novel lincRNAs have shorter transcript lengths, longer exon lengths, and fewer exon numbers compared with protein-coding transcripts, which is consistent with some previous studies (18–20, 27, 34–36, 43–46).
Differential expression analysis identified 272 DEGs, with 147 genes upregulated and 125 genes downregulated in the high-fat group compared to the low-fat group. Some of these genes may play a key role in lipid metabolism. As expected, GO analysis revealed the involvement of a significant number of DEGs in lipid metabolism-related biological processes, including response to lipid, cellular response to lipid, glycerophospholipid metabolic process, glycerolipid metabolic process, unsaturated fatty acid biosynthetic process, and phospholipid dephosphorylation. KEGG analysis showed significant enrichment of 62 KEGG pathways in the DEGs. Enriched pathways included the sphingolipid signaling pathway, MAPK signaling pathway, fatty acid metabolism, glycerophospholipid metabolism, biosynthesis of unsaturated fatty acids, PI3K-Akt signaling pathway, fatty acid degradation, ether lipid metabolism, cholesterol metabolism, PPAR signaling pathway, TGF-beta signaling pathway, and Wnt signaling pathway, which are involved in lipid metabolism.
Previous studies indicated that lncRNAs can regulate gene expression in a cis-acting manner (18, 20, 29, 35, 41). In the present study, a total of 196 novel lincRNAs and 7,615 known lncRNAs, were identified. Moreover, 52 DELs were detected, 25 of which were upregulated and 27 downregulated in the high-fat group compared with the low-fat group. To predict the function of these lncRNAs, protein-coding genes transcribed near lncRNAs (<100 kb) were screened. A total of 323 PTGs were identified, corresponding to 52 DELs. Exploration of lncRNA function through GO and KEGG analysis of PTGs revealed that 102 of 323 PTGs were significantly involved in 143 biological processes, and 59 PTGs significantly participated in 7 pathways, including 2 pathways related to lipid metabolism, such as biosynthesis of unsaturated fatty acids and fatty acid metabolism. In addition, four other pathways are related to lipid metabolism, namely the Wnt signaling pathway, glycerolipid metabolism, PPAR signaling pathway, and glycerophospholipid metabolism. These pathways can affect lipid metabolism to some extent.
An intriguing observation is that the same DEL ENSEAST00005078264 acts on both FADS1 and FADS2, which are the key enzymes that catalyze adenylation of flavin mononucleotide (FMN) to form flavin adenine dinucleotide (FAD) coenzyme (47, 48). These two PTGs are involved in the biosynthesis of unsaturated fatty acids and fatty acid metabolism. The DEL ENSEAST00005078264 may be involved in lipid metabolism via biosynthesis of unsaturated fatty acids and fatty acid metabolism signaling pathway. Another DEL ENSEAST00005069204 was found to target AGPAT5, which converts lysophosphatidic acid to phosphatidic acid, the second step in de novo phospholipid biosynthesis (49). These lncRNAs may regulate target genes involved in lipid metabolism pathways through cis-acting mechanisms.
The analysis of PTGs by trans-acting lncRNAs was conducted, and enriched analysis was also performed. Among the 3,366 PTGs, 1,352 were significantly involved in 835 biological processes. Pathway analysis showed that 710 PTGs were significantly involved in 74 pathways, including MAPK signaling pathway, PI3K-Akt signaling pathway, glycerolipid metabolism, ether lipid metabolism, fat digestion, and absorption, all of which are related to lipid metabolism. Other pathways related to lipid metabolisms, such as glycerophospholipid metabolism, fatty acid metabolism, biosynthesis of unsaturated fatty acids, sphingolipid signaling pathway, cholesterol metabolism, adipocytokine signaling pathway, Wnt signaling pathway, fatty acid elongation, regulation of lipolysis in adipocytes, sphingolipid metabolism, TGF-beta signaling pathway, fatty acid biosynthesis, and fatty acid degradation, were also identified. For 132 DEPTG, 56 were significantly involved in 143 biological processes. Pathway analysis showed that 27 DEPTGs were significantly involved in 7 pathways, including arachidonic acid metabolism. Other pathways involved in lipid metabolism have also been identified, including biosynthesis of unsaturated fatty acids, fat digestion and absorption, ether lipid metabolism, fatty acid metabolism, MAPK signaling pathway, PPAR signaling pathway, glycerophospholipid metabolism, and PI3K-Akt signaling pathway.
It was found that SCD, targeted by four DELs (ENSEAST00005042127, ENSEAST00005051768, ENSEAST00005052324, and ENSEAST00005072263), is a key gene that regulates lipid metabolism. Studies have shown that SCD plays an important role in lipid biosynthesis (50–52). In addition, two DELs, ENSEAST00005052438 and ENSEAST00005076477, were all found to act on PLA2G3, which is involved in lipid metabolism and catalyzes the calcium-dependent hydrolysis of the sn-2 acyl bond of phospholipids to release arachidonic acid and lysophospholipids (53–56). It is noteworthy that THRSP was targeted by four DELs (ENSEAST00005042127, ENSEAST00005048792, ENSEAST00005051768, and ENSEAST00005052324). Previous studies showed that THRSP is important for the biosynthesis of triglycerides with medium-length fatty acid chains and plays a role in the regulation of lipogenesis, which may modulate lipogenesis by interacting with MID1IP1 and preventing its interaction with ACACA, may function as a transcriptional coactivator and may modulate the transcription factor activity of THRB (57, 58). These lncRNAs may regulate target genes involved in lipid metabolism pathways through trans-acting mechanisms.
SCD plays an important role in the regulation of lipid deposition, and it can catalyze the conversion of saturated fatty acids (SFAs) to monounsaturated fatty acids (MUFAs) (59). It was found that SCD expression was positively correlated with SCD protein levels, and pigs with higher SCD protein levels had higher levels of IMF, indicating a conserved trend between transcriptional and translational levels (60). Smith et al. (61) found that SCD expression was closely related to marbled adipocyte differentiation and that grain diets increased SCD expression, resulting in higher intracellular levels of MUFA (especially oleic acid), and thus increased adipocyte differentiation. In addition, increased FA levels in IMF tissues were found to be associated with tissue-specific activation of SCD expression under the influence of a low-protein protein diet, and the low-protein protein diet significantly increased the expression and activity of SCD proteins in muscle, but not in subcutaneous adipose tissue (AT) (62). These results suggest that the high specific expression of SCD in IMF cells may be related to meat quality and that reduced protein intake may affect IMF content by modulating SCD activity. Notably, SCD was significantly up-regulated in the high-fat group compared with the low-fat group in this study. The results show that SCD may play an important function in donkey IMF deposition. The thyroid hormone-sensitive protein (THRSP; Spot14; S14) is a nuclear protein that is abundantly expressed in lipogenic tissues such as in the liver, mammary gland, AT and lipogenic breast cancers (63–66). A previous study indicated that the polymorphisms and genotype distribution of THRSP were closely related to the potential for fat production in pig breeds (67). It was found that miR-195 may inhibit lipid accumulation in adipocytes by regulating THRSP (68). THRSP is regulated by insulin both in vivo in human AT and in vitro in adipocytes, and its expression is downregulated by insulin resistance. As THRSP silencing decreases mitochondrial respiration and fatty acid oxidation, its downregulation in human AT could contribute to mitochondrial dysfunction. Furthermore, disturbed sphingolipid metabolism could contribute to metabolic dysfunction in obese AT (69). The current RNA-seq results revealed that THRSP had higher expression in the high-fat group than in the low-fat group, which suggested that THRSP may be a key gene involved in IMF deposition in donkeys. Taken together, some DEGs, including SCD and THRSP, may be key candidate genes for donkey meat quality improvement.
It is noteworthy that SCD was targeted by four DELs (ENSEAST00005042127, ENSEAST00005051768, ENSEAST00005052324, and ENSEAST00005072263). Moreover, four DELs (ENSEAST00005042127, ENSEAST00005048792, ENSEAST00005051768, and ENSEAST00005052324) were found to target THRSP. These lncRNAs may influence donkey IMF deposition by regulating the expression of their target genes. The relationship of these lncRNAs with their target genes and with fat deposition in donkey muscle needs to be further investigated in the future using molecular biology, gene editing, etc. These lncRNAs may improve meat quality and facilitate the selection process of donkeys in future breeding.
5 Conclusion
In the study, we identified 196 putative lncRNAs and analyzed the characteristics of lncRNAs compared with protein-coding genes in the longissimus dorsi muscles of Guangling donkeys. We observed numerous DELs and protein-coding genes in longissimus dorsi muscles with different IMF contents. Some DEGs were found to be involved in various biological processes related to lipid metabolism. Functional enrichment analysis of PTGs by DELs revealed that some lncRNAs (such as ENSEAST00005042127, ENSEAST00005051768, ENSEAST00005052324, ENSEAST00005072263, ENSEAST00005052438, ENSEAST00005076477, ENSEAST00005042127, ENSEAST00005048792, ENSEAST00005051768, and ENSEAST00005052324) may act on PTGs (such as SCD, PLA2G3, and THRSP), participate in lipid metabolism processes, and regulate IMF deposition in the longissimus dorsi muscle. This study provides valuable resources for future analyses of lipid deposition traits and may contribute to the improvement of donkey meat quality and the selection process in donkey breeding.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal studies were approved by Animal Ethics Committee of Liaocheng University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
YP: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing. MZ: Data curation, Investigation, Methodology, Visualization, Writing – review & editing. YG: Data curation, Formal analysis, Methodology, Software, Writing – review & editing. CW: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Shandong Province Modern Agricultural Industrial Technology System Project (grant number SDAIT-27), the Construction of Molecular ID Card for Donkey and Camel Species, the Livestock and Poultry Seed Industry Project of the Ministry of Agriculture and Rural Affairs (grant number 19211162), the Key R&D Project of Shandong Province: Innovation and Demonstration of Key Technologies for the Integrated Development of the Donge Black Donkey Industry (grant number 2021TZXD012), and the Liaocheng University Scientific Research Fund (grant number 318052339).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1410109/full#supplementary-material
Footnotes
References
1. Zhang, K, Huang, K, Luo, Y, and Li, S. Identification and functional analysis of long non-coding RNAs in mouse cleavage stage embryonic development based on single cell transcriptome data. BMC Genomics. (2014) 15:845. doi: 10.1186/1471-2164-15-845
2. Li, H, Wang, Z, Zhao, B, Zhang, H, Fan, D, Ma, H, et al. Sperm-borne lncRNA loc100847420 improves development of early bovine embryos. Anim Reprod Sci. (2023) 257:107333. doi: 10.1016/j.anireprosci.2023.107333
3. Ransohoff, JD, Wei, Y, and Khavari, PA. The functions and unique features of long intergenic non-coding RNA. Nat Rev Mol Cell Biol. (2018) 19:143–57. doi: 10.1038/nrm.2017.104
4. Natarajan, P, Shrinivas, K, and Chakraborty, AK. A model for cis-regulation of transcriptional condensates and gene expression by proximal LncRNAs. Biophys J. (2023) 122:2757–72. doi: 10.1016/j.bpj.2023.05.032
5. Loewer, S, Cabili, MN, Guttman, M, Loh, YH, Thomas, K, Park, IH, et al. Large intergenic non-coding RNA-RoR modulates reprogramming of human induced pluripotent stem cells. Nat Genet. (2010) 42:1113–7. doi: 10.1038/ng.710
6. He, XY, Fan, X, Qu, L, Wang, X, Jiang, L, Sang, LJ, et al. LncRNA modulates hippo-yap signaling to reprogram iron metabolism. Nat Commun. (2023) 14:2253. doi: 10.1038/s41467-023-37871-5
7. Prickett, AR, and Oakey, RJ. A survey of tissue-specific genomic imprinting in mammals. Mol Gen Genomics. (2012) 287:621–30. doi: 10.1007/s00438-012-0708-6
8. Di Michele, F, Chillon, I, and Feil, R. Imprinted long non-coding RNAs in mammalian development and disease. Int J Mol Sci. (2023) 24:13647. doi: 10.3390/ijms241713647
9. Zhang, M, Li, F, Sun, JW, Li, DH, Li, WT, Jiang, RR, et al. LncRNA IMFNCR promotes intramuscular adipocyte differentiation by sponging miR-128-3p and miR-27b-3p. Front Genet. (2019) 10:42. doi: 10.3389/fgene.2019.00042
10. Xiao, F, Tang, CY, Tang, HN, Wu, HX, Hu, N, Li, L, et al. Long non-coding RNA 332443 inhibits Preadipocyte differentiation by targeting Runx1 and P38-Mapk and Erk1/2-Mapk signaling pathways. Front Cell Dev Biol. (2021) 9:663959. doi: 10.3389/fcell.2021.663959
11. Li, P, Ruan, X, Yang, L, Kiesewetter, K, Zhao, Y, Luo, H, et al. A liver-enriched Long non-coding RNA, lncLSTR, regulates systemic lipid metabolism in mice. Cell Metab. (2015) 21:455–67. doi: 10.1016/j.cmet.2015.02.004
12. Li, W, Qiu, L, Guan, J, Sun, Y, Zhao, J, and Du, M. Comparative transcriptome analysis of longissimus Dorsi tissues with different intramuscular fat contents from Guangling donkeys. BMC Genomics. (2022) 23:644. doi: 10.1186/s12864-022-08857-2
13. van Laack, RL, Stevens, SG, and Stalder, KJ. The influence of ultimate PH and intramuscular fat content on pork tenderness and tenderization. J Anim Sci. (2001) 79:392–7. doi: 10.2527/2001.792392x
14. Jeremiah, LE, Dugan, ME, Aalhus, JL, and Gibson, LL. Assessment of the chemical and cooking properties of the major beef muscles and muscle groups. Meat Sci. (2003) 65:985–92. doi: 10.1016/S0309-1740(02)00308-X
15. Michal, JJ, Zhang, ZW, Gaskins, CT, and Jiang, Z. The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu X Limousin F2 crosses. Anim Genet. (2006) 37:400–2. doi: 10.1111/j.1365-2052.2006.01464.x
16. Bolger, AM, Lohse, M, and Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. (2014) 30:2114–20. Epub 2014/04/04. doi: 10.1093/bioinformatics/btu170
17. Pertea, M, Kim, D, Pertea, GM, Leek, JT, and Salzberg, SL. Transcript-level expression analysis of RNA-Seq experiments with HISAT, Stringtie and Ballgown. Nat Protoc. (2016) 11:1650–67. doi: 10.1038/nprot.2016.095
18. Chen, L, Shi, G, Chen, G, Li, J, Li, M, Zou, C, et al. Transcriptome analysis suggests the roles of Long intergenic non-coding RNAs in the growth performance of weaned piglets. Front Genet. (2019) 10:196. doi: 10.3389/fgene.2019.00196
19. Chen, G, Cheng, X, Shi, G, Zou, C, Chen, L, Li, J, et al. Transcriptome analysis reveals the effect of Long intergenic noncoding RNAs on pig muscle growth and fat deposition. Biomed Res Int. (2019) 2019:2951427. doi: 10.1155/2019/2951427
20. Shi, G, Chen, L, Chen, G, Zou, C, Li, J, Li, M, et al. Identification and functional prediction of Long intergenic non-coding RNAs related to subcutaneous adipose development in pigs. Front Genet. (2019) 10:160. doi: 10.3389/fgene.2019.00160
21. Zhao, W, Li, Z, Liu, Q, Xie, S, Li, M, Wang, Y, et al. Analysis of Long intergenic non-coding RNAs transcriptomic profiling in skeletal muscle growth during porcine embryonic development. Sci Rep. (2021) 11:15240. doi: 10.1038/s41598-021-94014-w
22. Kang, YJ, Yang, DC, Kong, L, Hou, M, Meng, YQ, Wei, L, et al. Cpc2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res. (2017) 45:W12–6. doi: 10.1093/nar/gkx428
23. Prakash, A, Jeffryes, M, Bateman, A, and Finn, RD. The Hmmer web server for protein sequence similarity search. Curr Protoc Bioinformatics. (2017) 60:3 15 1–3 23. doi: 10.1002/cpbi.40
24. Pirooznia, M, Perkins, EJ, and Deng, Y. Batch blast extractor: an automated Blastx parser application. BMC Genomics. (2008) 9:S10. doi: 10.1186/1471-2164-9-S2-S10
25. Love, MI, Huber, W, and Anders, S. Moderated estimation of fold change and dispersion for RNA-Seq data with Deseq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
26. Benjamini, Y, Drai, D, Elmer, G, Kafkafi, N, and Golani, I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res. (2001) 125:279–84. doi: 10.1016/S0166-4328(01)00297-2
27. Zou, C, Li, L, Cheng, X, Li, C, Fu, Y, Fang, C, et al. Identification and functional analysis of Long intergenic non-coding RNAs underlying intramuscular fat content in pigs. Front Genet. (2018) 9:102. doi: 10.3389/fgene.2018.00102
28. Quinlan, AR, and Hall, IM. Bedtools: a flexible suite of utilities for comparing genomic features. Bioinformatics. (2010) 26:841–2. doi: 10.1093/bioinformatics/btq033
29. Cabili, MN, Trapnell, C, Goff, L, Koziol, M, Tazon-Vega, B, Regev, A, et al. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. (2011) 25:1915–27. doi: 10.1101/gad.17446611
30. Yu, G, Wang, LG, Han, Y, and He, QY. Clusterprofiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
31. Wu, T, Hu, E, Xu, S, Chen, M, Guo, P, Dai, Z, et al. Clusterprofiler 4.0: a universal enrichment tool for interpreting Omics data. Innovation (Camb). (2021) 2:100141. doi: 10.1016/j.xinn.2021.100141
32. Xie, C, Mao, X, Huang, J, Ding, Y, Wu, J, Dong, S, et al. Kobas 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. (2011) 39:W316–22. doi: 10.1093/nar/gkr483
33. Benjamini, Y, and Cohen, R. Weighted false discovery rate controlling procedures for clinical trials. Biostatistics. (2017) 18:91–104. doi: 10.1093/biostatistics/kxw030
34. Zou, C, Li, J, Luo, W, Li, L, Hu, A, Fu, Y, et al. Transcriptome analysis reveals Long intergenic non-coding RNAs involved in skeletal muscle growth and development in pig. Sci Rep. (2017) 7:8704. doi: 10.1038/s41598-017-07998-9
35. Zou, C, Li, S, Deng, L, Guan, Y, Chen, D, Yuan, X, et al. Transcriptome analysis reveals Long intergenic noncoding RNAs contributed to growth and meat quality differences between Yorkshire and Wannanhua pig. Genes (Basel). (2017) 8:203. doi: 10.3390/genes8080203
36. Gong, Y, He, J, Li, B, Xiao, Y, Zeng, Q, Xu, K, et al. Integrated analysis of lncRNA and mRNA in subcutaneous adipose tissue of Ningxiang pig. Biology (Basel). (2021) 10:726. doi: 10.3390/biology10080726
37. Gil, N, and Ulitsky, I. Regulation of gene expression by cis-acting long non-coding RNAs. Nat Rev Genet. (2020) 21:102–17. doi: 10.1038/s41576-019-0184-5
38. Zhou, H, Simion, V, Pierce, JB, Haemmig, S, Chen, AF, and Feinberg, MW. Lncrna-Map3k4 regulates vascular inflammation through the P38 Mapk Signaling pathway and cis-modulation of Map3k4. FASEB J. (2021) 35:e21133. doi: 10.1096/fj.202001654RR
39. Szafranski, P, Gambin, T, Karolak, JA, Popek, E, and Stankiewicz, P. Lung-specific distant enhancer cis regulates expression of Foxf1 and lncRNA FENDRR. Hum Mutat. (2021) 42:694–8. doi: 10.1002/humu.24198
40. Degani, N, Lubelsky, Y, Perry, RB, Ainbinder, E, and Ulitsky, I. Highly conserved and cis-acting Lncrnas produced from paralogous regions in the Center of Hoxa and Hoxb Clusters in the endoderm lineage. PLoS Genet. (2021) 17:e1009681. doi: 10.1371/journal.pgen.1009681
41. Hao, Z, Fan, C, Cheng, T, Su, Y, Wei, Q, and Li, G. Genome-wide identification, characterization and evolutionary analysis of Long intergenic noncoding Rnas in cucumber. PLoS One. (2015) 10:e0121800. doi: 10.1371/journal.pone.0121800
42. Feng, H, Liu, T, Yousuf, S, Zhang, X, Huang, W, Li, A, et al. Identification and analysis of lncRNA, miRNA and mRNA related to subcutaneous and intramuscular fat in Laiwu pigs. Front Endocrinol. (2022) 13:1081460. doi: 10.3389/fendo.2022.1081460
43. Yu, H, Zhao, X, and Li, Q. Genome-wide identification and characterization of long intergenic noncoding RNAs and their potential association with larval development in the Pacific oyster. Sci Rep. (2016) 6:20796. doi: 10.1038/srep20796
44. Wu, Y, Cheng, T, Liu, C, Liu, D, Zhang, Q, Long, R, et al. Systematic identification and characterization of long non-coding RNAs in the silkworm, Bombyx Mori. PLoS One. (2016) 11:e0147147. doi: 10.1371/journal.pone.0147147
45. Wang, J, Fu, L, Koganti, PP, Wang, L, Hand, JM, Ma, H, et al. Identification and functional prediction of large intergenic noncoding RNAs (linCRNAs) in rainbow trout (Oncorhynchus Mykiss). Mar Biotechnol (NY). (2016) 18:271–82. doi: 10.1007/s10126-016-9689-5
46. Zhang, YC, Liao, JY, Li, ZY, Yu, Y, Zhang, JP, Li, QF, et al. Genome-wide screening and functional analysis identify a large number of Long noncoding RNAs involved in the sexual reproduction of Rice. Genome Biol. (2014) 15:512. doi: 10.1186/s13059-014-0512-1
47. Brizio, C, Galluccio, M, Wait, R, Torchetti, EM, Bafunno, V, Accardi, R, et al. Over-expression in Escherichia coli and characterization of two recombinant isoforms of human fad synthetase. Biochem Biophys Res Commun. (2006) 344:1008–16. doi: 10.1016/j.bbrc.2006.04.003
48. Olsen, RKJ, Konarikova, E, Giancaspero, TA, Mosegaard, S, Boczonadi, V, Matakovic, L, et al. Riboflavin-responsive and -non-responsive mutations in fad synthase cause multiple acyl-CoA dehydrogenase and combined respiratory-chain deficiency. Am J Hum Genet. (2016) 98:1130–45. doi: 10.1016/j.ajhg.2016.04.006
49. Prasad, SS, Garg, A, and Agarwal, AK. Enzymatic activities of the human Agpat isoform 3 and isoform 5: localization of Agpat5 to mitochondria. J Lipid Res. (2011) 52:451–62. doi: 10.1194/jlr.M007575
50. Wang, J, Yu, L, Schmidt, RE, Su, C, Huang, X, Gould, K, et al. Characterization of HSCD5, a novel human Stearoyl-CoA desaturase unique to primates. Biochem Biophys Res Commun. (2005) 332:735–42. doi: 10.1016/j.bbrc.2005.05.013
51. Goren, MA, and Fox, BG. Wheat germ cell-free translation, purification, and assembly of a functional human Stearoyl-CoA desaturase complex. Protein Expr Purif. (2008) 62:171–8. doi: 10.1016/j.pep.2008.08.002
52. Zhang, S, Yang, Y, and Shi, Y. Characterization of human SCD2, an oligomeric desaturase with improved stability and enzyme activity by cross-linking in intact cells. Biochem J. (2005) 388:135–42. doi: 10.1042/BJ20041554
53. Murakami, M, Masuda, S, Shimbara, S, Bezzine, S, Lazdunski, M, Lambeau, G, et al. Cellular Arachidonate-releasing function of novel classes of secretory phospholipase A2s (groups iii and xii). J Biol Chem. (2003) 278:10657–67. doi: 10.1074/jbc.M211325200
54. Sato, H, Kato, R, Isogai, Y, Saka, G, Ohtsuki, M, Taketomi, Y, et al. Analyses of group iii secreted phospholipase A2 transgenic mice reveal potential participation of this enzyme in plasma lipoprotein modification, macrophage foam cell formation, and atherosclerosis. J Biol Chem. (2008) 283:33483–97. doi: 10.1074/jbc.M804628200
55. Murakami, M, Masuda, S, Shimbara, S, Ishikawa, Y, Ishii, T, and Kudo, I. Cellular distribution, post-translational modification, and tumorigenic potential of human group iii secreted phospholipase a(2). J Biol Chem. (2005) 280:24987–98. doi: 10.1074/jbc.M502088200
56. Murase, R, Taketomi, Y, Miki, Y, Nishito, Y, Saito, M, Fukami, K, et al. Group iii phospholipase a(2) promotes colitis and colorectal cancer. Sci Rep. (2017) 7:12261. doi: 10.1038/s41598-017-12434-z
57. Chou, WY, Cheng, YS, Ho, CL, Liu, ST, Liu, PY, Kuo, CC, et al. Human spot 14 protein interacts physically and functionally with the thyroid receptor. Biochem Biophys Res Commun. (2007) 357:133–8. doi: 10.1016/j.bbrc.2007.03.103
58. Chou, WY, Ho, CL, Tseng, ML, Liu, ST, Yen, LC, and Huang, SM. Human spot 14 protein is a P53-dependent transcriptional coactivator via the recruitment of thyroid receptor and ZAC1. Int J Biochem Cell Biol. (2008) 40:1826–34. doi: 10.1016/j.biocel.2008.01.014
59. Dumas, SN, Guo, CA, Kim, JK, Friedline, RH, and Ntambi, JM. Interleukin-6 derived from cutaneous deficiency of Stearoyl-CoA desaturase- 1 may mediate metabolic organ crosstalk among skin, adipose tissue and liver. Biochem Biophys Res Commun. (2019) 508:87–91. doi: 10.1016/j.bbrc.2018.11.083
60. Zappaterra, M, Deserti, M, Mazza, R, Braglia, S, Zambonelli, P, and Davoli, R. A gene and Protein expression study on four porcine genes related to intramuscular fat deposition. Meat Sci. (2016) 121:27–32. doi: 10.1016/j.meatsci.2016.05.007
61. Smith, SB, Kawachi, H, Choi, CB, Choi, CW, Wu, G, and Sawyer, JE. Cellular regulation of bovine intramuscular adipose tissue development and composition. J Anim Sci. (2009) 87:E72–82. doi: 10.2527/jas.2008-1340
62. Doran, O, Moule, SK, Teye, GA, Whittington, FM, Hallett, KG, and Wood, JD. A reduced protein diet induces Stearoyl-CoA desaturase protein expression in pig muscle but not in subcutaneous adipose tissue: relationship with intramuscular lipid formation. Br J Nutr. (2006) 95:609–17. doi: 10.1079/BJN20051526
63. Freake, HC, and Moon, YK. Hormonal and nutritional regulation of Lipogenic enzyme mRNA levels in rat primary white and Brown adipocytes. J Nutr Sci Vitaminol (Tokyo). (2003) 49:40–6. doi: 10.3177/jnsv.49.40
64. Freake, HC, and Oppenheimer, JH. Stimulation of S14 mRNA and lipogenesis in Brown fat by hypothyroidism, cold exposure, and cafeteria feeding: evidence supporting a general role for S14 in lipogenesis and lipogenesis in the maintenance of thermogenesis. Proc Natl Acad Sci USA. (1987) 84:3070–4. doi: 10.1073/pnas.84.9.3070
65. Jump, DB. Rapid induction of rat liver S14 gene transcription by thyroid hormone. J Biol Chem. (1989) 264:4698–703. doi: 10.1016/S0021-9258(18)83799-1
66. Jump, DB, Narayan, P, Towle, H, and Oppenheimer, JH. Rapid effects of triiodothyronine on hepatic gene expression. Hybridization analysis of tissue-specific triiodothyronine regulation of Mrnas14. J Biol Chem. (1984) 259:2789–97. doi: 10.1016/S0021-9258(17)43215-7
67. Wang, X, Cheng, J, Qin, W, Chen, H, Chen, G, Shang, X, et al. Polymorphisms in 5’ proximal regulating region of THRSP gene are associated with fat production in pigs. 3 Biotech. (2020) 10:267. doi: 10.1007/s13205-020-02266-6
68. Wang, S, Pan, C, Ma, X, Yang, C, Tang, L, Huang, J, et al. Identification and functional verification reveals that MiR-195 inhibiting THRSP to affect fat deposition in Xinyang Buffalo. Front Genet. (2021) 12:736441. doi: 10.3389/fgene.2021.736441
Keywords: Guangling donkey, lncRNA, longissimus dorsi muscle, lipid deposition, regulate
Citation: Peng Y, Zhu M, Gong Y and Wang C (2024) Identification and functional prediction of lncRNAs associated with intramuscular lipid deposition in Guangling donkeys. Front. Vet. Sci. 11:1410109. doi: 10.3389/fvets.2024.1410109
Edited by:
Zhihong Liu, Inner Mongolia Agricultural University, ChinaReviewed by:
Gan Shangquan, Xinjiang Academy of Agricultural and Reclamation Sciences (XAARS), ChinaJianning He, Qingdao Agricultural University, China
Copyright © 2024 Peng, Zhu, Gong and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yongdong Peng, cHlvbmdkb25nMTIzQDEyNi5jb20=; Changfa Wang, d2FuZ2NoYW5nZmFAbGN1LmVkdS5jbg==