Changsheng He
Yong Wang1,2
Jiangjiang Zhu- 1Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization of Education Ministry, Southwest Minzu University, Chengdu, China
- 2Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Exploitation of Sichuan Province, Southwest Minzu University, Chengdu, China
- 3College of Animal and Veterinary Science, Southwest Minzu University, Chengdu, China
- 4College of Food Science and Technology, Southwest Minzu University, Chengdu, China
Goats are popular in China because of their superior meat quality, delicate flesh, and unique flavor. Long noncoding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation of gene expression. However, the effects of lncRNAs on adipocyte differentiation in goat has not been fully elucidated yet. In this investigation, we performed RNA-Seq analysis of intramuscular and subcutaneous adipocytes from Jianzhou Daer goat before and after differentiation, including both intramuscular preadipocytes (IMPA) vs. intramuscular adipocytes (IMA) and subcutaneous preadipocytes (SPA) vs. subcutaneous adipocytes (SA). A total of 289.49 G clean reads and 12,519 lncRNAs were obtained from 20 samples. In total, 3,733 differentially expressed RNAs (182 lncRNAs and 3,551 mRNAs) were identified by pairwise comparison. There were 135 differentially expressed lncRNAs (DELs) specific to intramuscular adipocytes, 39 DELs specific to subcutaneous adipocytes, and 8 DELs common to both adipocytes in these 182 DELs. Some well-known and novel pathways associated with preadipocyte differentiation were identified: fat acid metabolism, TGF-beta signaling pathway and PI3K-Akt signaling pathway. By integrating miRNA-seq data from another study, we also identified hub miRNAs in both types of fat cells. Our analysis revealed the unique and common lncRNA-miRNA-mRNA networks of two kinds of adipocytes. Several lncRNAs that regulate potentially goat preadipocyte differentiation were identified, such as XR_001918 647.1, XR_001917728.1, XR_001297263.2 and LNC_004191. Furthermore, our findings from the present study may contribute to a better understanding of the molecular mechanisms underlying in goat meat quality and provide a theoretical basis for further goat molecular breeding.
Introduction
Fat is generally divided into subcutaneous fat, intermuscular fat and intramuscular fat. The amount of intermuscular fat is small and its composition is similar to subcutaneous fat, so it is generally considered that fat is mainly divided into subcutaneous fat and intramuscular fat (Timón et al., 2001). The distribution of fat in meat animals and the fat content of each part have an important impact on meat production and quality. Subcutaneous fat determines the carcass lean rate to a certain extent. The thicker the fat means the lower the lean meat rate. Intramuscular fat has an important influence on the flavor and tenderness of goat meat. Within a certain range, the quality of meat gradually improves with the increase of intramuscular fat content. One of the current research directions for meat is to study how to regulate the content of intramuscular fat and subcutaneous fat to improve meat quality of goats.
Adipogenesis is a complex process regulated by various transcription factors, non-coding RNA and signal pathways (Mattar et al., 2018). RNA sequencing (RNA-seq) is a revolutionary tool to identify differentially expressed genes (DEGs) regulating various biological processes. It enables us to discover new genes and therefore to describe unannotated transcriptional activity by identifying numerous noncoding transcripts (Wang et al., 2016). Long non-coding RNAs (lncRNAs) are non-coding RNAs with a length greater than 200 nucleotides and low protein coding ability. MicroRNAs (miRNAs) are a type of endogenous non-coding RNA with a length of approximately 22 nucleotides. Previous studies have shown that lncRNAs and mRNAs containing the same miRNA binding site can regulate mutual expression levels by competitively binding miRNAs. The recent explosion in knowledge demonstrating the importance of miRNAs and lncRNAs in the regulation of multiple major biological processes, which impacts the development, differentiation, and metabolism have brought these neglected molecular players (Huang et al., 2011; Wang and Chang, 2011; Ren et al., 2016). Some studies have shown that the adipogenic differentiation ability of intramuscular adipocytes is significantly lower than that of subcutaneous adipocytes. The expression of related genes is low, indicating that there are specific regulatory mechanisms for intramuscular fat and subcutaneous fat in animals (Zhou et al., 2010; Sun et al., 2013). In mammals, the differentiation of preadipocytes has been well studied, especially in bovine and porcine (Chen et al., 2016; Dong et al., 2016; Sun et al., 2016; Yonekura et al., 2016; Zappaterra et al., 2016). Recently, goat meat is gradually welcomed by consumers because of its high protein content, low fat and cholesterol content. Moreover, the consumers demand for goat meat and quality requirements continue to rise (Ivanovic et al., 2016; Teixeira et al., 2019). However, there are few studies on the molecular mechanism of the difference in the deposition of different fat tissues in goats.
Herein, we provided a comprehensive transcriptome profile on intramuscular and subcutaneous adipocytes in before and after differentiation of Jianzhou Daer goat. The results revealed the expression patterns of mRNAs and lncRNAs, which were important in the developmental stages of two different adipocytes. An integrated analysis of differentially expressed lncRNAs (DELs) and mRNAs (DEGs) was performed, and the lncRNA and its target miRNA/mRNA that could potentially regulate the differentiation of intramuscular and subcutaneous preadipocytes were screened through bioinformatics analysis. At the same time, a lncRNA-miRNA-mRNA interaction network was constructed based on miRNAs known to play a role in adipogenesis, which were identified by another study (manuscript in preparation). The lncRNA that may be combined with miRNA and mRNA was identified and verified by RT-qPCR technique. The results from the present investigation provided a theoretical basis for in-depth analysis of the regulation mechanism of goat fat cell differentiation and improvement of meat quality in goats.
Materials and Methods
Experimental Animals and Sample Collection
All experimental procedures were reviewed and approved by the Institutional Animal Care and Use Committee, Southwest Minzu University. Also, all the experiments complied with the requirements of the directory of the Ethical Treatment of Experimental Animals of China.
The 7-day-old male Jianzhou Daer goat (Capra hircus) (n = 5) was purchased from Sichuan Jianyang Dageda Aminal Husbandry Co., Ltd (Sichuan, China), being euthanized by bloodletting. The longissimus dorsi and subcutaneous fat were excised from the goats and minced. The isolation and culture of goat intramuscular and subcutaneous preadipocytes was performed as described assay by (Xu et al., 2018a). The intramuscular and subcutaneous adipocytes were isolated by using twice the volume of Type I collagenase (Sigma, St. Louis, MO, United States) in PBS (Hyclone, Logan, UT, United States). When the cell confluence reaches 80%, we discarded the medium (DMEM/F12 (Hyclone, Logan, UT, United States) with 10% fetal bovine serum (FBS) (Hyclone, Logan, UT, United States)) and added appropriate amount of trypsin (Hyclone, Logan, UT, United States), and the cells were digested at 37°C for 1 min. And then, the cells were harvested in a centrifuge tube, followed by a centrifuge at 800 rpm/min for 3 min to obtain cell pellet. The collected cells were re-suspended for subsequent experiments. The cells of passage 3, were used for experimental treatment. The cells were coaxed into differentiating using oleic acid induction solution (DMEM/F12 (Hyclone, Logan, UT, United States) containing 10% FBS (Hyclone, Logan, UT, United States) and 50 μm/L oleic acid (Sigma, St. Louis, MO, United States), and harvested after induction for 0 and 3 days (Xiong et al., 2018; Song et al., 2020; Xiong et al., 2021; Li et al., 2022). These cell samples were named intramuscular preadipocytes (IMPA), subcutaneous preadipocytes (SPA), intramuscular adipocytes (IMA) and subcutaneous adipocytes (SA), respectively, and five biological replicates were set for each group.
Total RNA Extraction and Sequencing
A total of 20 cell samples were successfully collected. Total RNA was extracted using Trizol reagent (Takara, Dalian, China). RNA quality was determined using NanoPhotometer spectrophotometer and Agilent 2,100 bioanalyzer, which analyzes the integrity of the RNA （RIN）. The lowest RIN accepted for RNA analyses is 6.8. The lncRNA library is constructed using a chain-specific library. The method for synthesizing the first strand of cDNA by reverse transcription is the same as the normal method of NEB library construction. The difference is that when the second strand is synthesized, dTTP in dNTPs is replaced by dUTP. After that, cDNA end repair, A-tailing, ligation of sequencing adapters, and length screening were also performed, and then the second strand of cDNA containing U was degraded by USER enzyme, and then PCR amplification was performed to obtain a library. Finally, twenty libraries were sequenced at Novogene Co. Ltd. (Beijing, China) on Illumina HiSeq Sequencing System. The RNA-Seq dataset supporting the conclusions of this article is available in the Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information (NCBI) under accession numbers GSE186988.
Quality Control and Transcript Assembly
After sequencing, the raw data were stored in fastq format (de Sena Brandine and Smith, 2019). We removed reads containing connectors, poly-N and low quality from Raw Reads to get clean reads. The Q20, Q30 and GC content of the clean data were calculated (Bai et al., 2015). The obtained clean, high-quality data were used for further analyses. The clean, paired-end reads were aligned to the goat genome sequence assembly using HISAT2 (v2.0.4) (Kim et al., 2015), and the transcripts were assembled using both Scripture (beta2) (Guttman et al., 2010) and Cufflinks (v2.1.1) (Trapnell et al., 2010).
Identification of lncRNAs and Their Expression Analysis
To ensure the quality of the obtained lncRNAs, three criteria were used to identify the desired lncRNAs in the transcriptome assemblies: 1) transcripts with length >200 bp and exon number ≥2 were selected; 2) Cuffcompare (v2.1.1) was used to calculate the read coverage of every transcript, and transcripts with an FPKM value (Cuffquant) of more than 0.05 were removed; and 3) the coding potential of the transcripts was predicted using the coding potential calculator (CPC <1) (Kong et al., 2007), Coding Noncoding Index (CNCI, v2) (Sun et al., 2013), and Pfam (v1.3) (Finn et al., 2016) protein domain families to further remove coding genes. Subsequent steps were performed based on the intersection of the results obtained using the twenty databases.
The FPKMs values of lncRNAs and mRNAs in IMPA, IMA, SPA and SA libraries were calculated by StringTie (v1.3.1), Ballgown (Pertea et al., 2016) and Cuffdiff (v2.1.1) (Trapnell C et al., 2010). The transcripts in the 20 libraries were analyzed for differential expression, and the p value was used to screen the DELs and DEGs between two different adipocyte samples. When p < 0.05, the lncRNAs and mRNAs are considered to be differentially expressed.
Gene Functional Annotation
The trans function describes the co-expression relationship between lncRNAs and mRNAs. Pearson correlation coefficient (R > 0.95 or R < 0.95) was calculated by custom scripts. David was used to cluster the target genes among 20 samples for functional enrichment analysis of lncRNA target genes (Huang et al., 2009). Adjusted p-value < 0.05 was set as significant threshold. Go (gene ontology) can enrich and analyze the target genes of DELs. Goseq (Release2.12) (Young et al., 2010) is used to enrich and analyze the target genes of DELs. Adjusted p < 0.05 is considered to be the significant enrichment of GEGs. KEGG is a database (http://www.genome.jp/kegg/) for understanding the advanced functions and utilities of biological systems such as cells, organisms, and ecosystems (Kanehisa and Goto, 2000; Kanehisa et al., 2008; Kanehisa, 2019; Kanehisa et al., 2021). We used KOBAS (v2.0) (Mao et al., 2005) software to detect the pathway enrichment analysis of DELs target genes in the KEGG pathway.
Differentially Expressed miRNA and Construction of lncRNA-miRNA-mRNA Networks
We have previously established small RNA libraries by RNA-Seq before and after differentiation of intramuscular and subcutaneous adipocytes (data not yet published). The raw reads were filtered to remove low quality reads and reads with connectors to obtain clean reads. Novoaligen software was used to match the clean reads with the miRBase database to identify known miRNAs, and MirDeep software (Xie et al., 2012) was used to predict novel miRNAs were finally quantified using Novoaligen and Samtools software and normalised using Reads Per Million (RPM). DEGSeq software (Wang et al., 2010) was used for DEMs analysis, and P adjust values were used to indicate the significance of differentially expressed genes. The edgeR software (Robinson et al., 2010) was used to calculate DEMs, in two adipocytes. The screening criterion for significant differences is p-value <0.05. miRWalk, miRanda, RNAhybrid, and Targetscan were used for screening DEMs target DELs and DEGs. Then, the intersections of co-expressed lncRNAs, miRNAs, and mRNAs in the four programs were used to construct a lncRNA–miRNA–mRNA network.
DELs and DEGs were combined with the target lncRNAs and mRNAs of DEMs respectively to obtain the miRNA-mRNA relationship pair and miRNA-lnRNA relationship pair. Finally, the lncRNA–miRNA–mRNA network was visualized using Cytoscape v3.7 software (Shannon et al., 2003).
Validation of Gene Expression of by RT-qPCR Technique
Primers were designed using Primer-BLAST on the NCBI website (Table 1). 1 μg of RNA was synthesized into the first cDNA strand using the RevertAid First Strand cDNA Synthesis Kit (Thermo, Waltham, MA, United States), in accordance with the user manual. Ubiquitously expressed transcript gene (UXT) was used as a housekeeping gene (Xu et al., 2018b). The qPCR reaction procedure was composed of four steps, including pre-degeneration (95°C, 3 min), degeneration (95°C, 10 s), annealing (60°C, 10 s) and extension (72°C, 15 s), of which degeneration, annealing and extension were running for 40 cycles. Melt curve stage was 65–95°C in 0.5°C increments for 5s. Quantification of selected gene expression was performed using the comparative threshold cycle (2−ΔΔCT) method (Livak and Schmittgen, 2001). The experiment was repeated for three times.
TABLE 1. Primer sequences of lncRNA amplified by qRT-PCR.
Results
Transcript Sequencing and Assembly
Three days after differentiation of intramuscular and subcutaneous adipogenesis, lipid droplets could be observed with Oil Red O staining (Figures 1A,B) and the relative expression level of the adipocyte differentiation marker gene Pref-1 was significantly lower than that of preadipocytes (Figure 1C), which could indicate that differentiated intramuscular and subcutaneous adipocyte models were successfully established. We established twenty cDNA libraries that represented two different adipocytes: intramuscular preadipocytes (IMPA, 1–5) and intramuscular adipocytes (IMA, 1–5) from the longissimus dorsi of 7-day-old Jianzhou Daer goat, subcutaneous preadipocytes (SPA, 1–5) and subcutaneous adipocytes (SA, 1–5) from the subcutaneous fat of 7-day-old Jianzhou Daer goat. The RNA sequencing obtained a total of 289.49 Gb of data, with each stage averaging 14.4745 Gb of data. The Q30 results in each sample were >91%, and the GC percentage was less than 52%, as listed in Table 2. More than 47.76% of the clean reads were perfectly mapped to the goat reference genome (assembly ARS1, https://www.ncbi.nlm.nih.gov/genome/?term = goat), and 44.56–87.33% uniquely mapped reads were obtained from the total mapped reads from the twenty samples (Table 3).
FIGURE 1. Adipocytes differentiated for 3 days (A) Oil-red O staining analysis in goat intramuscular adipocyte adipogenesis. “*”means p < 0.05 and “**” means p < 0.01, same below. (B) Oil-red O staining analysis in goat subcutaneous adipocyte adipogenesis. (C) The mRNA level of Pref-1 in goat intramuscular and subcutaneous adipocyte adipogenesis.
TABLE 2. Output statistics of the sequencing reads for each sample.
TABLE 3. Summary of the clean reads alignment to the goat reference genome.
Finally, 12519 assumed non-coding transcripts were retained by CNCI, CPC2 and PFAM softwares (Figure 2A), including 16.9% lncRNAs, 76.1% intronic lncRNAs, and 7% anti-sense lncRNAs. In addition, 43767 mRNAs were identified. The expression level of each RNA was standardized to fragments per kilobase of exon model per million mapped reads (FPKM), and it was found that a small part of them was not expressed or was expressed at a relatively low level. The expression levels of a large number of RNAs were mainly between 0.1 and 2 of log10 (FPKM+1), and a small number of genes had very high expression level (log10 |FPKM+ 1| > 2) (Figure 2B). LncRNAs are lesser expressed than mRNAs (Figure 2C), the extremely significant differences of genes were helpful to explore the molecular mechanism of fat deposition in different tissues.
FIGURE 2. Identification and characterization of lncRNAs in goat, Capra hircus. (A) Venn diagrams of coding potential analysis by using stringent criteria. Three tools (CPC2, CNCI and PFAM) were employed to analyze the coding potential of lncRNAs. Those simultaneously shared by three analytical tools were designated as candidate lncRNAs and used in subsequent analysis. (B) Volcanic maps of differentially expressed RNA. The expression levels of mRNAs and lncRNAs were indicated by log10 (FPKM +1). (C) LncRNA, mRNA expression violin graph. The horizontal coordinate is the molecular type, the vertical coordinate is log10 (FPKM+1) and the width of each violin plot reflects the number of points at that expression level. (D) Length distribution of mRNAs and lncRNAs, unit of the length is bp. (E) Distribution of exon number in the mRNAs and lncRNAs. (F) Open reading frame (ORF) length distribution for mRNAs and lncRNAs.
Characteristics of lncRNAs and mRNAs in Intramuscular and Subcutaneous Adipocytes of Goats Before and After Differentiation
Since most lncRNAs are produced by RNA polymerase II transcription, they have structural features similar to those of mRNAs (Gao et al., 2017). Therefore, in order to observe the characteristics of lncRNAs and the difference between lncRNAs and mRNAs, we compared the length of lncRNAs and mRNAs, including the number of exons, and the open reading frame (ORF). The results showed that most of lncRNAs tended to be shorter in length, and they contained less exons than mRNAs (Figure 2D). Also, the length of ORFs in the lncRNAs was shorter than those of the mRNAs (Figures 2E,F).
Differentially Expressed lncRNAs and mRNAs in Goat Different Adipocytes
The correlation coefficient can represent the degree of similarity between samples. We found that the data between SPA and SA, and between IMPA and IMA were highly correlated in expression (Figure 3A). We found 143 lncRNAs and 33258 mRNAs differentially expressed between IMA and IMPA (73 lncRNAs and 1,378 mRNAs were up-regulated, 70 lncRNAs and 1880 mRNAs were down-regulated), 47 lncRNAs and 1,659 mRNAs differentially expressed between SA and SPA (26 lncRNAs and 825 mRNAs were up-regulated, 21 lncRNAs and 834 mRNAs were down-regulated) (Figure 3B) (Supplements 1, 2, 5, 6). The PCA score plots for each of mRNA and lncRNA showed a high degree of similarity in the same group of samples (Figure 3C). To further explore the differences in lncRNAs expression between subcutaneous and intramuscular adipocytes at different developmental stages, we performed a clustered heat map on the differentially expressed genes. The results of the cluster analysis showed that distinct lncRNAs expression patterns are associated with the differentiation of both subcutaneous and intramuscular adipocytes in goats (Figure 3D).
FIGURE 3. Differentially expressed coding genes and lncRNAs in goat adipocytes between intramuscular and subcutaneous of differentiation (Fold Change ≥2 and P-adjust value ≤ 0.05). (A) Representative of the degree of similarity between samples. The correlation coefficient is represented by color; deeper color represents a stronger correlation. (B) Clustering analysis of differentially expressed lncRNAs and mRNAs in IMA vs. IMPA and SA vs. SPA. Red means highly expressed genes and blue means low expressed genes. The color is from red to blue, which means log10(FPKM+1) from large to small. (C) PCA plots of differentially expressed lncRNAs and mRNAs in IMA vs. IMPA and SA vs. SPA. Scatter plots are drawn using principal component 1 and principal component 2 as the x- and y-axes.(D) Volcano plot of differentially expressed lncRNAs and mRNAs in IMA vs. IMPA and SA vs. SPA. Green and red represent downregulated and upregulated expression, respectively.
Function Prediction of lncRNAs and Corresponding Genes During Intramuscular and Subcutaneous Preadipocyte Differentiation
LncRNAs not only can regulate the expression of neighboring protein-coding genes through a cis mechanism (Bu et al., 2012; Han et al., 2012), but also regulate the expression of genes located on other chromosomes via a trans mechanism (Yang et al., 2016; Cai et al., 2017). In this study, co-expression analysis was used to predict the potential target genes of DELs via a trans mechanism. The absolute value of Pearson’s correlation coefficient was set as greater than 0.95 (Supplement 7). There were 39012 lncRNA-mRNA relationship pairs. Among them, LNC_004374, LNC_007272, LNC_010037.
LNC_010037, LNC_004066, LNC_010143 and other lncRNAs can target genes related to adipocyte differentiation or lipid metabolism, CD36 (Lewis, 2006), FABP3 (Nakamura et al., 2013), FGF11 (Lee et al., 2019), FOXO6 (Lee et al., 2019), SMAD1 (Shin et al., 2016), TGFB2 (Takahashi et al., 2019), FGFR2 (Fischer et al., 2017). The results are shown in Table 4.
TABLE 4. Expression of some lncRNAs and their associated genes affecting adipogenesis.
We drew Venn diagrams (Figure 4A) to visualize the number of differential lncRNAs that are common versus unique to the two adipocyte differential differentiation processes. 135 DELs were specific to intramuscular adipocytes, 39 were specific to subcutaneous adipocytes, and 8 DELs were common to both adipocytes. Through the use of Goseq, we compare the GO classifications of DELs using mRNAs showing a correlated expression. (adjusted p < 0.05) (Young et al., 2010). The DELs target genes are enriched in biological process (BP), cellular component (CC) and molecular function (MF) (Table 5). Unique DELs target genes of goat intramuscular adipocytes before and after differentiation are mainly enriched in binding (Figure 4B), and the unique DELs target genes of subcutaneous adipocytes are mainly enriched in glucose metabolism and chemokine (Figure 4C), while the shared DELs target genes are mainly enriched in binding (Figure 4D).
FIGURE 4. GO and KEGG enrichment analyses of differentially expressed lncRNAs in goat different adipocytes. (A) Venny maps of differentially expressed lncRNAs during intramuscular and subcutaneous preadipocyte differentiation. (B) Categories of biological processes, cellular components and molecular functions of the target genes of differentially expressed lncRNAs (intramuscular adipocytes). (C) Categories of the biological processes and molecular functions of the target genes of differentially expressed lncRNAs (subcutaneous adipocytes). (D) Categories of molecular functions of the target genes of differentially expressed lncRNAs (both subcutaneous and intramuscular adipocytes). (E) Scatter plot of the top 20 pathways enriched for differentially expressed lncRNAs in intramuscular adipocytes before and after differentiation. (F) Scatter plot of the top 20 pathways enriched for differentially expressed lncRNAs in subcutaneous adipocytes before and after differentiation. (G) Scatter plot of the top 15 pathways enriched for differentially expressed lncRNAs in both intramuscular and subcutaneous adipocytes before and after differentiation. The abscissa represents the richness factor, and the ordinate represents the enriched pathway terms. The Q-value represents the corrected p value.
TABLE 5. GO enrichment analysis of the target genes of IncRNAs in intramuscular and subcutaneous adipocytes.
The KEGG annotation results of the differentially expressed target genes of lncRNAs are classified according to the pathway types of the KEGG database. As shown in Figure 4E, the unique DELs target genes of intramuscular adipocytes before and after differentiation are most significantly enriched in hypertrophic cardiomyopathy (HCM), dilated cardiomyopathy, focal adhesion, arrhythmogenic right ventricular cardiomyopathy (ARVC), ECM-receptor interaction, proteoglycans in cancer and cardiac muscle contraction (adjusted p < 0.05). In addition, the analyzed genes also was enriched in fat acid metabolism and TGF-beta signaling pathway, which are related to adipocyte formation (Supplement 10). The unique DELs target genes of subcutaneous adipocytes before and after differentiation are most significantly enriched in biosynthesis of amino acids, glycolysis/Gluconeogenesis, carbon metabolism, valine, leucine and isoleucine degradation, histidine metabolism, fatty acid degradation, melanoma, fatty acid metabolism (adjusted p < 0.05), among which fatty acid degradation, fatty acid metabolism, TGF-beta signaling pathway, PI3K-Akt signaling pathway have a significant impact on the adipogenesis (Figure 4F). The shared DELs target genes are most significantly enriched in biosynthesis of amino acids, carbon metabolism, glycolysis/gluconeogenesis and alzheimer’s disease (adjusted p < 0.05) (Figure 4G).
Construction of Regulation Networks for Adipocytes in Goats
As to miRNAs analysis, 210 known miRNAs and 95 novel miRNAs were found to have significantly different expression between intramuscular preadipocytes and adipocytes (Supplement 3), 175 known miRNAs and 67 novel miRNAs were found to have significantly different expression between subcutaneous preadipocytes and adipocytes (Supplement 4). Among them, we obtained 146 differentially expressed miRNAs (DEMs) unique to intramuscular adipocytes, 83 (DEMs) unique to subcutaneous adipocytes, and 159 (DEMs) shared by two adipocytes (Figure 5A). After targeted binding prediction, 30,6571 mRNAs (Supplement 9) and 32,029 lncRNAs.
FIGURE 5. Venny maps of differentially expressed RNAs. Venny maps showing specific miRNAs (A) and mRNAs (B), respectively, in intramuscular adipocytes or subcutaneous adipocytes.
32,029 lncRNAs (Supplement 8) were found to have potential targeted binding relationships with DEMs. Combined with the sequenced intramuscular adipocyte-specific 2,724 DEGs and 135 DELs, 6,572 miRNA-mRNA relationship pairs and 1,240 miRNA-lncRNA relationship pairs were obtained (Figure 5B). The 1,125 DEGs and 39 DELs unique to subcutaneous adipocytes were combined to obtain 3,287 miRNA-mRNA relationship pairs and 307 miRNA-lncRNA relationship pairs (Figure 5B). At the same time, the intersection with the shared 534 DEGs and 8 DELs yielded 1,681 miRNA-mRNA relationship pairs and 52 miRNA-lncRNA relationship pairs. The lncRNA-miRNA-mRNA interaction network was constructed by Cytoscape software, and the top 20 miRNAs were selected after calculating the degree value of each factor (Table 6). Among them, we found miR-20, miR-194, miR-335, miR-363, miR-200, miR-199, and miR-302 related to fat formation, thereby, we constructed a ceRNA network view. As shown in Figure 6A, the unique lncRNA-miRNA-mRNA network of intramuscular adipocytes includes 30 lncRNAs, six miRNAs and 426 mRNAs, and the lncRNA-miRNA-mRNA network unique to subcutaneous adipocytes includes 8 lncRNAs, four miRNAs and 182 mRNAs (Figure 6B). Whereas, two types of adipocyte-common lncRNA-miRNA-mRNA networks include 1 lncRNA, 1 miRNAs and 26 mRNAs (Figure 6C).
TABLE 6. Node degree of mRNA-miRNA-lncRNA regualatory network.
FIGURE 6. lncRNA–miRNA–mRNA network. (A) Regulatory network of intramuscular adipocytes before and after differentiation. Nodes with Closeness Centrality >6.4815E-4 are discarded. Closeness Centrality is the reciprocal of the sum of the shortest paths from a point to all other points, and the larger it is, the shorter the path from that point to all other points and the closer it is to the center in space. (B) Regulatory network of subcutaneous adipocytes before and after differentiation. (C) There is a common network of intramuscular and subcutaneous adipocytes before and after differentiation.
Verification of lncRNA Expression Profiles Using qRT-PCR
To verify the reliability of RNA-seq results, seven candidate lncRNAs were randomly selected from the DELs obtained from the screening, and UXT was used as an internal reference for qRT-PCR analysis. The results showed that XR_001917557.1, and LNC_004191 were significantly down-regulated in expression during the differentiation of intramuscular adipose tissue. Whereas, the XR_001918647.1, XR_001917728.1 were significantly up-regulated in expression during the differentiation of intramuscular adipose tissue. XR_001295810.1, XR_001917637.1 and XR_001297263.2 were significantly down-regulated expression during the differentiation of subcutaneous adipose tissue, and the LNC_004191 was significantly up-regulated expression during the differentiation of subcutaneous adipose tissue, These results were consistent with the trend of RNA-seq, indicating the credibility of the RNA-seq results (Figure 7).
FIGURE 7. Validation of lncRNAs in intramuscular and subcutaneous adipocytes using RT-qPCR technique. Data were analyzed by the 2−ΔΔCt method where the UXT was used as a reference gene. Each column represents the mean ± SE. White bars represent read from RNA-Seq, and blue bars represent the results from qRT-PCR analysis.
Discussion
The distribution and deposition of adipose tissue in different parts of the body are the key factors affecting carcass quality and meat flavor. Subcutaneous fat mainly affects carcass quality. Intramuscular fat (IMF) is the material basis of marbling, and an important factor affecting meat flavor. A large number of studies have shown that IMF is directly involved in the formation of meat tenderness, juiciness and flavor (Van Laack et al., 2001; Suzuki et al., 2005). Goat is an indispensable animal in China’s agricultural production, and the molecular regulation mechanism of its lipid deposition has not been fully elucidated yet.
LncRNA is a kind of noncoding RNA longer than 200 nt, which has attracted substantial attention in the last few years. Studies have shown that lncRNAs regulate metabolic tissue development and function, including adipogenesis, hepatic lipid metabolism, islet function, and energy balance (Alvarez-Dominguez et al., 2015; Chen et al., 2015; Luan et al., 2015; You et al., 2015; Zhao and Lin, 2015). Despite the fact that many studies have indicated the importance of lncRNAs in different tissues, little is known about their biological function in goat fat deposition, especially in the differentiation of goat intramuscular and subcutaneous preadipocytes. To the best of our knowledge, our study is the first to screen for lncRNAs and mRNAs regulating goat preadipocyte differentiation by sequencing and annotating the transcriptome of intramuscular and subcutaneous preadipocytes. A total of 1,118,110,544 reads were successfully mapped to the goat reference genome assembly. We identified 12,519 lncRNAs. The average sequence length of lncRNAs was shorter than that of mRNAs, and the number of exons was less than that of mRNAs, with the ORF length being shorter than that of mRNAs. Our results indicated that the predicted lncRNAs were shorter with fewer exons than mRNAs, which are in agreement with the those reproted in previous studies (Trapnell et al., 2010; Ran et al., 2016; Wang et al., 2016). The Pearson correlation (R2) of each sample is greater than 0.8, which indicated that our experiment was reliable and the sample selection was reasonable.
LncRNA functions by regulating mRNA. At present, the mechanism of interaction between lncRNA and mRNA is not clear. We predict the biological function of lncRNAs through its co-expression with protein coding genes. Consequently, we found that many target genes of DELs were also differentially expressed in goat intramuscular and subcutaneous preadipocytes. This suggested that lncRNAs may function through complementary target genes, which can play critical roles in the differentiation of goat intramuscular and subcutaneous preadipocytes. For example, SMAD1 is a target gene of the differentially expressed lncRNAs LNC_009792, LNC_007731, LNC_000706, LNC_008467, LNC_006192 and LNC_004878, and it has been reported to regulate the differentiation of preadipocytes (Shin et al., 2016). These findings suggest that these lncRNAs might be involved in the differentiation of intramuscular preadipocytes by affecting the expression of SMAD1. However, there is no DEGs in the differentiation process of subcutaneous preadipocytes in our selected fat development related genes, so we speculate that the network regulating intramuscular adipogenesis is more complex than subcutaneous fat. The higher number of DELs, DEMs and DEGs in IMF than that of subcutaneous fat also supported this, which is consistent with the fact that intramuscular preadipocytes have stronger ability to deposit fat than that of subcutaneous preadipocytes (Chen et al., 2010).
To explore the similarities and differences of different adipocytes, DELs target genes (IMPA vs. IMA and SPA vs. SA) were subjected to GO and KEGG pathway enrichment analyses. We found that few common term was found between the IMPA vs. IMA and SPA vs. SA comparisons. Several pathways involved in preadipocyte differentiation were previously identified, including the TGF-β signaling pathway (IMF and subcutaneous fat) (Stewart et al., 2010), PI3K/AKT signaling pathway (subcutaneous fat) (Dong et al., 2016), and arachidonic acid metabolism (subcutaneous fat) (Nikolopoulou et al., 2014). For example, the LncRNA GAS5 inhibits lipogenesis in 3T3-L1 cells through the miR-21a-5p/PTEN signaling pathway (Liu et al., 2018). FDNCR1 affects porcine lipogenesis by competitively binding miR-204 to regulate the TGF-β pathway (Zhang et al., 2020). However, for some pathways identified here, their involvement in the goat preadipocyte differentiation process is being reported for the first time. Interestingly, in the pathway analysis, we found that the components of two pathways, fatty acid metabolism (IMF and subcutaneous fat) and fatty acid degradation (subcutaneous fat), which have been reported to be involved in lipid metabolism, and were enriched in the entire process of differentiation of intramuscular and subcutaneous preadipocytes (Wang et al., 2006; Jump, 2011). The common enrichment pathways during differentiation of both adipocytes involve amino acid metabolism, gluconeogenesis and carbohydrate metabolism, suggesting that IMF and subcutaneous fat are largely different in differentiation pathways and lipid metabolism pathways, while there are similarities in communication with other metabolic pathways. In addition to the specific pathways of the two adipocytes, there are also differences in the common components of the two pathways. Just as IMF is enriched in 19 target genes of TGF-β signaling pathway, 8 target genes of fatty acid metabolism, and subcutaneous fat is enriched in 8 and 7 targat genes. Here, we hypothesize that there are differences in the pathways regulating intramuscular and subcutaneous adipose differentiation, and that there are differences in the downstream target genes in the same pathways, which indirectly demonstrates that gene expression is tissue-specific in goats.
To date, many genes have been reported to regulate the differentiation of preadipocytes. However, few studies have been conducted on the roles of lncRNAs in intramuscular and subcutaneous preadipocytes differentiation. The molecular and cellular mechanisms regulating goat preadipocytes differentiation are thus still poorly understood. Here, we constructed a miRNA-lncRNA-mRNA interaction network, and calculated the degree (the number of times each factor interacts with other factors) of each factor through Cytoscape (CentiScaPe (Scardoni et al., 2009)). We selected the top 20 miRNAs in the dgree, and found that miR-20 (Wang et al., 2008), miR-194 (Jeong et al., 2014), miR-335 (Fernandez-Hernando et al., 2011), miR-363 (Chen et al., 2014), miR-200 (Kennell et al., 2008), miR-199 (Shi et al., 2014), and miR-302 (Kim et al., 2014) were related to fat development, and visualized them in Cytoscape. We identified a number of highly connected lncRNAs and mRNAs in the three modules, including the two kinds of adipocytes are unique and common. For example, XM_005693834.3 (ACLY) (Han et al., 2016) and XM_013975359.2 (ANGPT2) (Bae et al., 2020) are unique to IMF, and XM_018056618.1 (MEDAG) (Li et al., 2021) and XM_005678357.3(PLPP3) (Bae et al., 2016) are unique to subcutaneous fat, whereas, XM_018040030.1 (CAPN10) (Patel and Lane, 1999) and XM_018050348.1 (TGFβ1) (Stewart et al., 2010) are shared by the two kinds of preadipocytes in differentiation processes. Seven DELs from the IMF vs. subcutaneous fat comparison were validated by qRT-PCR technique and the results were in excellent generally agreement with the RNA-seq findings. This suggests that our RNA-seq findings are reliable. Triangles miRNA-mRNA-lncRNA should be represented by a gene expression negatively correlated between miRNA-mRNA and miRNA-lncRNA and a positively correlated between lncRNA and mRNA as discussed in a pure sponge prediction module (Paci et al., 2014). Therefore, we constructed several examples of a pure sponge module based on qPCR-validated lncRNAs (Figure 8). Some notable phenomena such as: it is expected that both XR_001917637.1 and XR_001297263.2 should be upregulated in subcutaneous adipocytes to promote the expression of PLPP3 and MEDAG, yet both lncRNAs are actually downregulated in subcutaneous adipocytes. This phenomenon suggests that the triangle of miRNA-mRNA-lncRNA is a complex regulatory network, with negatively correlated gene expression between miRNA-lncRNAs or positively correlated gene expression between lncRNA-mRNAs contributing to this difference in final outcome. Our next experimental plan is to investigate the interactions between mRNA, lncRNA, and miRNA.
FIGURE 8. Example of a pure sponge prediction module. (A) A pure sponge prediction module extracted from intramuscular adipocytes network. (B) A pure sponge prediction module extracted from intramuscular adipocytes network.
In conclusion, we first generated the expression profiles of lncRNAs from intramuscular and subcutaneous adipocytes of Jianzhou Daer goat (IMA vs. IMPA and SA vs. SPA) based on RNA-Seq technique. We found that the number of lncRNAs regulating IMF differentiation was more than that of subcutaneous adipocytes. Our results suggest that those lncRNAs might play important roles in adipocyte differentiation. Collectively, this study takes the first step toward understanding the molecular mechanisms underlying variations in goat adipogenesis. Also, our results provided a theoretical basis for molecular breeding to improve the meat quality in goats.
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 below: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE186988 GEO accession: GSE186988.
Ethics Statement
The animal study was reviewed and approved by the Animal Experimental Ethical Inspection of Southwest University for Nationalities. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
CH and YL conceived the research; CH performed the exp -eriments; YW and JZ shared regents; CH and YL analysed the data; CH. JC and YL wrote the paper. All authors read and approved the final manuscript.
Funding
This work were supported by National Natural Sciences Foundation of China (No. 3207 2723), supported by Sichuan Science and Technology Program(2022JDTD0030), and Fundamental Research Funds for the Central Universities, Southwest Minzu University (2021057).
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.
Acknowledgments
The authors acknowledge assistance from the goat lipid metabolism research lab.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.900179/full#supplementary-material
References
Alvarez-Dominguez J. R., Bai Z., Xu D., Yuan B., Lo K. A., Yoon M. J., et al. (2015). De Novo Reconstruction of Adipose Tissue Transcriptomes Reveals Long Non-coding RNA Regulators of Brown Adipocyte Development. Cell Metab. 21, 764–776. doi:10.1016/j.cmet.2015.04.003
Bae H., Hong K. Y., Lee C.-k., Jang C., Lee S.-J., Choe K., et al. (2020). Angiopoietin-2-integrin α5β1 Signaling Enhances Vascular Fatty Acid Transport and Prevents Ectopic Lipid-Induced Insulin Resistance. Nat. Commun. 11, 2980. doi:10.1038/s41467-020-16795-4
Bae Y. J., Kim S.-E., Hong S. Y., Park T., Lee S. G., Choi M.-S., et al. (2016). Time-course Microarray Analysis for Identifying Candidate Genes Involved in Obesity-Associated Pathological Changes in the Mouse Colon. Genes Nutr. 11, 30–41. doi:10.1186/s12263-016-0547-x
Bai B., Wu J., Sheng W.-T., Zhou B., Zhou L.-J., Zhuang W., et al. (2015). Comparative Analysis of Anther Transcriptome Profiles of Two Different Rice Male Sterile Lines Genotypes under Cold Stress. Ijms 16, 11398–11416. doi:10.3390/ijms160511398
Bu Q., Hu Z., Chen F., Zhu R., Deng Y., Shao X., et al. (2012). Transcriptome Analysis of Long Non-coding RNAs of the Nucleus Accumbens in Cocaine-Conditioned Mice. J. Neurochem. 123, 790–799. doi:10.1111/jnc.12006
Cai B., Li Z., Ma M., Wang Z., Han P., Abdalla B. A., et al. (2017). LncRNA-Six1 Encodes a Micropeptide to Activate Six1 in Cis and Is Involved in Cell Proliferation and Muscle Growth. Front. Physiol. 8, 230. doi:10.3389/fphys.2017.00230
Chen J., Cui X., Shi C., Chen L., Yang L., Pang L., et al. (2015). Differential lncRNA Expression Profiles in Brown and White Adipose Tissues. Mol. Genet. Genomics 290, 699–707. doi:10.1007/s00438-014-0954-x
Chen J., Dodson M. V., Jiang Z. (2010). Cellular and Molecular Comparison of Redifferentiation of Intramuscular- and Visceral-Adipocyte Derived Progeny Cells. Int. J. Biol. Sci. 6, 80–88. doi:10.7150/ijbs.6.80
Chen L., Cui J., Hou J., Long J., Li C., Liu L. (2014). A Novel Negative Regulator of Adipogenesis: microRNA-363. Stem Cells 32, 510–520. doi:10.1002/stem.1549
Chen X., Zhou B., Luo Y., Huang Z., Jia G., Liu G., et al. (2016). Tissue Distribution of Porcine FTO and its Effect on Porcine Intramuscular Preadipocytes Proliferation and Differentiation. PloS one 11, e0151056. doi:10.1371/journal.pone.0151056
de Sena Brandine G., Smith A. D. (2019). Falco: High-Speed FastQC Emulation for Quality Control of Sequencing Data. F1000Res 8, 1874. doi:10.12688/f1000research.21142.2
Dong X., Tang S., Zhang W., Gao W., Chen Y. (2016). GPR39 Activates Proliferation and Differentiation of Porcine Intramuscular Preadipocytes through Targeting the PI3K/AKT Cell Signaling Pathway. J. Recept. Signal Transduct. 36, 130–138. doi:10.3109/10799893.2015.1056308
Fernández-Hernando C., Suárez Y., Rayner K. J., Moore K. J. (2011). MicroRNAs in Lipid Metabolism. Curr. Opin. Lipidol. 22, 86–92. doi:10.1097/MOL.0b013e3283428d9d
Finn R. D., Coggill P., Eberhardt R. Y., Eddy S. R., Mistry J., Mitchell A. L., et al. (2016). The Pfam Protein Families Database: towards a More Sustainable Future. Nucleic Acids Res. 44, D279–D285. doi:10.1093/nar/gkv1344
Fischer C., Seki T., Lim S., Nakamura M., Andersson P., Yang Y., et al. (2017). A miR-327-FGF10-FGFR2-Mediated Autocrine Signaling Mechanism Controls White Fat Browning. Nat. Commun. 8, 2079. doi:10.1038/s41467-017-02158-z
Gao X., Ye J., Yang C., Zhang K., Li X., Luo L., et al. (2017). Screening and Evaluating of Long Noncoding RNAs in the Puberty of Goats. BMC genomics 18, 164. doi:10.1186/s12864-017-3578-9
Guttman M., Garber M., Levin J. Z., Donaghey J., Robinson J., Adiconis X., et al. (2010). Ab Initio reconstruction of Cell Type-specific Transcriptomes in Mouse Reveals the Conserved Multi-Exonic Structure of lincRNAs. Nat. Biotechnol. 28, 503–510. doi:10.1038/nbt.1633
Han J., Li L., Wang D., Ma H. (2016). (−)-Hydroxycitric Acid Reduced Fat Deposition via Regulating Lipid Metabolism-Related Gene Expression in Broiler Chickens. Lipids Health Dis. 15, 37. doi:10.1186/s12944-016-0208-5
Han L., Zhang K., Shi Z., Zhang J., Zhu J., Zhu S., et al. (2012). LncRNA Profile of Glioblastoma Reveals the Potential Role of lncRNAs in Contributing to Glioblastoma Pathogenesis FIle of Glioblastoma Reveals the Potential Role of lncRNAs in Contributing to Glioblastoma Pathogenesis. Int. J. Oncol. 40, 2004–2012. doi:10.3892/ijo.2012.1413
Huang D. W., Sherman B. T., Lempicki R. A. (2009). Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 37, 1–13. doi:10.1093/nar/gkn923
Huang Y., Shen X. J., Zou Q., Wang S. P., Tang S. M., Zhang G. Z. (2011). Biological Functions of microRNAs: a Review. J. Physiol. Biochem. 67, 129–139. doi:10.1007/s13105-010-0050-6
Ivanovic S., Pavlovic I., Pisinov B. (2016). The Quality of Goat Meat and It's Impact on Human Health. Biotec Anim. Husb. 32, 111–122. doi:10.2298/bah1602111i
Jeong B.-C., Kang I.-H., Hwang Y.-C., Kim S.-H., Koh J.-T. (2014). MicroRNA-194 Reciprocally Stimulates Osteogenesis and Inhibits Adipogenesis via Regulating COUP-TFII Expression. Cell Death Dis. 5–e1532. doi:10.1038/cddis.2014.485
Jump D. B. (2011). Fatty Acid Regulation of Hepatic Lipid Metabolism. Curr. Opin. Clin. Nutr. Metabolic Care 14, 115–120. doi:10.1097/MCO.0b013e328342991c
Kanehisa M., Araki M., Goto S., Hattori M., Hirakawa M., Itoh M., et al. (2008). KEGG for Linking Genomes to Life and the Environment. Nucleic Acids Res. 36, D480–D484. doi:10.1093/nar/gkm882
Kanehisa M., Furumichi M., Sato Y., Ishiguro-Watanabe M., Tanabe M. (2021). KEGG: Integrating Viruses and Cellular Organisms. Nucleic Acids Res. 49, D545–D551. doi:10.1093/nar/gkaa970
Kanehisa M., Goto S. (2000). KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Kanehisa M. (2019). Toward Understanding the Origin and Evolution of Cellular Organisms. Protein Sci. 28, 1947–1951. doi:10.1002/pro.3715
Kennell J. A., Gerin I., MacDougald O. A., Cadigan K. M. (2008). The microRNA miR-8 Is a Conserved Negative Regulator of Wnt Signaling. Proc. Natl. Acad. Sci. U.S.A. 105, 15417–15422. doi:10.1073/pnas.0807763105
Kim D., Langmead B., Salzberg S. L. (2015). HISAT: a Fast Spliced Aligner with Low Memory Requirements. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kim J. Y., Shin K. K., Lee A. L., Kim Y. S., Park H. J., Park Y. K., et al. (2014). MicroRNA-302 Induces Proliferation and Inhibits Oxidant-Induced Cell Death in Human Adipose Tissue-Derived Mesenchymal Stem Cells. Cell Death Dis. 5, e1385. doi:10.1038/cddis.2014.344
Kong L., Zhang Y., Ye Z.-Q., Liu X.-Q., Zhao S.-Q., Wei L., et al. (2007). CPC: Assess the Protein-Coding Potential of Transcripts Using Sequence Features and Support Vector Machine. Nucleic Acids Res. 35, W345–W349. doi:10.1093/nar/gkm391
Lee K. W., Jeong J. Y., An Y. J., Lee J. H., Yim H. S. (2019). FGF11 Influences 3T3‐L1 Preadipocyte Differentiation by Modulating the Expression of PPARγ Regulators. FEBS Open Bio 9, 769–780. doi:10.1002/2211-5463.12619
Lewis G. F. (2006). Lipid Metabolism. Curr. Opin. Lipidol. 17, 205–208. doi:10.1097/01.mol.0000217906.36057.2b
Li X., Zhang H., Wang Y., Li Y., He C., Zhu J., et al. (2022). RNA-seq Analysis Reveals the Positive Role of KLF5 in the Differentiation of Subcutaneous Adipocyte in Goats. Gene 808, 145969. doi:10.1016/j.gene.2021.145969
Li Z., Li C., Wu Q., Tu Y., Wang C., Yu X., et al. (2021). MEDAG Enhances Breast Cancer Progression and Reduces Epirubicin Sensitivity through the AKT/AMPK/mTOR Pathway. Cell Death Dis. 12, 97. doi:10.1038/s41419-020-03340-w
Liu H., Li H., Jin L., Li G., Hu S., Ning C., et al. (2018). Long Noncoding RNA GAS5 Suppresses 3T3-L1 Cells Adipogenesis Through miR-21a-5p/PTEN Signal Pathway. DNA Cell Biol. 37, 767–777. doi:10.1089/dna.2018.4264
Livak K. J., Schmittgen T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Luan A., Paik K. J., Li J., Zielins E. R., Atashroo D. A., Spencley A., et al. (2015). RNA Sequencing for Identification of Differentially Expressed Noncoding Transcripts during Adipogenic Differentiation of Adipose-Derived Stromal Cells. Plastic Reconstr. Surg. 136, 752–763. doi:10.1097/PRS.0000000000001582
Mao X., Cai T., Olyarchuk J. G., Wei L. (2005). Automated Genome Annotation and Pathway Identification Using the KEGG Orthology (KO) as a Controlled Vocabulary. Bioinformatics 21, 3787–3793. doi:10.1093/bioinformatics/bti430
Mattar P., Bravo-Sagua R., Tobar N., Fuentes C., Troncoso R., Breitwieser G., et al. (2018). Autophagy Mediates Calcium-Sensing Receptor-Induced TNFα Production in Human Preadipocytes. Biochimica Biophysica Acta (BBA) - Mol. Basis Dis. 1864, 3585–3594. doi:10.1016/j.bbadis.2018.08.020
Nakamura Y., Sato T., Shiimura Y., Miura Y., Kojima M. (2013). FABP3 and Brown Adipocyte-Characteristic Mitochondrial Fatty Acid Oxidation Enzymes Are Induced in Beige Cells in a Different Pathway from UCP1. Biochem. Biophysical Res. Commun. 441, 42–46. doi:10.1016/j.bbrc.2013.10.014
Nikolopoulou E., Papacleovoulou G., Jean-Alphonse F., Grimaldi G., Parker M. G., Hanyaloglu A. C., et al. (2014). Arachidonic Acid-dependent Gene Regulation during Preadipocyte Differentiation Controls Adipocyte Potential. J. Lipid Res. 55, 2479–2490. doi:10.1194/jlr.M049551
Paci P., Colombo T., Farina L. (2014). Computational Analysis Identifies a Sponge Interaction Network between Long Non-coding RNAs and Messenger RNAs in Human Breast Cancer. BMC Syst. Biol. 8, 83. doi:10.1186/1752-0509-8-83
Patel Y. M., Lane M. D. (1999). Role of Calpain in Adipocyte Differentiation. Proc. Natl. Acad. Sci. U.S.A. 96, 1279–1284. doi:10.1073/pnas.96.4.1279
Pertea M., Kim D., Pertea G. M., Leek J. T., Salzberg S. L. (2016). Transcript-level Expression Analysis of RNA-Seq Experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi:10.1038/nprot.2016.095
Ran M., Chen B., Li Z., Wu M., Liu X., He C., et al. (2016). Systematic Identification of Long Noncoding RNAs in Immature and Mature Porcine Testes1. Biol. Reprod. 94, 77. doi:10.1095/biolreprod.115.136911
Ren H., Wang G., Chen L., Jiang J., Liu L., Li N., et al. (2016). Genome-wide Analysis of Long Non-coding RNAs at Early Stage of Skin Pigmentation in Goats (Capra hircus). BMC genomics 17, 67. doi:10.1186/s12864-016-2365-3
Robinson M. D., McCarthy D. J., Smyth G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616
Scardoni G., Petterlini M., Laudanna C. (2009). Analyzing Biological Network Parameters with CentiScaPe. Bioinformatics 25, 2857–2859. doi:10.1093/bioinformatics/btp517
Shannon P., Markiel A., Ozier O., Baliga N. S., Wang J. T., Ramage D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Shi X.-E., Li Y.-F., Jia L., Ji H.-L., Song Z.-Y., Cheng J., et al. (2014). MicroRNA-199a-5p Affects Porcine Preadipocyte Proliferation and Differentiation. Ijms 15, 8526–8538. doi:10.3390/ijms15058526
Shin S., Seong J. K., Bae Y. S. (2016). Ahnak Stimulates BMP2-Mediated Adipocyte Differentiation through Smad1 Activation. Obesity 24, 398–407. doi:10.1002/oby.21367
Song W., Zhong C., Yuan Y., Zhu Q., Wang Y., Yin H., et al. (2020). Peroxisome Proliferator-Activated Receptor-Coactivator 1-beta (PGC-1β) Modulates the Expression of Genes Involved in Adipogenesis during Preadipocyte Differentiation in Chicken. Gene 741, 144516. doi:10.1016/j.gene.2020.144516
Stewart A., Guan H., Yang K. (2010). BMP-3 Promotes Mesenchymal Stem Cell Proliferation through the TGF-Beta/activin Signaling Pathway. J. Cell Physiol. 223, 658–666. doi:10.1002/jcp.22064
Sun W. X., Dodson M. V., Jiang Z. H., Yu S. G., Chu W. W., Chen J. (2016). Myostatin Inhibits Porcine Intramuscular Preadipocyte Differentiation In Vitro. Domest. Anim. Endocrinol. 55, 25–31. doi:10.1016/j.domaniend.2015.10.005
Sun W. X., Wang H. H., Jiang B. C., Zhao Y. Y., Xie Z. R., Xiong K., et al. (2013). Global Comparison of Gene Expression between Subcutaneous and Intramuscular Adipose Tissue of Mature Erhualian Pig. Genet. Mol. Res. 12, 5085–5101. doi:10.4238/2013.october.29.3
Suzuki K., Irie M., Kadowaki H., Shibata T., Kumagai M., Nishida A. (2005). Genetic Parameter Estimates of Meat Quality Traits in Duroc Pigs Selected for Average Daily Gain, Longissimus Muscle Area, Backfat Thickness, and Intramuscular Fat Content. J. Anim. Sci. 83, 2058–2065. doi:10.2527/2005.8392058x
Takahashi H., Alves C. R. R., Stanford K. I., Middelbeek R. J. W., Nigro P., Ryan R. E., et al. (2019). TGF-β2 Is an Exercise-Induced Adipokine that Regulates Glucose and Fatty Acid Metabolism. Nat. Metab. 1, 291–303. doi:10.1038/s42255-018-0030-7
Teixeira A., Silva S., Rodrigues S. (2019). Advances in Sheep and Goat Meat Products Research. Adv. Food Nutr. Res. 87, 305–370. doi:10.1016/bs.afnr.2018.09.002
Timón M. L., Ventanas J., Carrapiso A. I., Jurado A., Garcı́a C. (2001). Subcutaneous and Intermuscular Fat Characterisation of Dry-Cured Iberian Hams. Meat Sci. 58, 85–91. doi:10.1016/s0309-1740(00)00136-4
Trapnell C., Williams B. A., Pertea G., Mortazavi A., Kwan G., van Baren M. J., et al. (2010). Transcript Assembly and Quantification by RNA-Seq Reveals Unannotated Transcripts and Isoform Switching during Cell Differentiation. Nat. Biotechnol. 28, 511–515. doi:10.1038/nbt.1621
Van Laack R. L., Stevens S. G., Stalder K. J. (2001). The Influence of Ultimate pH and Intramuscular Fat Content on Pork Tenderness and Tenderization. J. Anim. Sci. 79, 392–397. doi:10.2527/2001.792392x
Wang K. C., Chang H. Y. (2011). Molecular Mechanisms of Long Noncoding Rnas. Mol. Cell 43, 904–914. doi:10.1016/j.molcel.2011.08.018
Wang L., Feng Z., Wang X., Wang X., Zhang X. (2010). DEGseq: an R Package for Identifying Differentially Expressed Genes from RNA-Seq Data. Bioinformatics 26, 136–138. doi:10.1093/bioinformatics/btp612
Wang Q., Li Y. C., Wang J., Kong J., Qi Y., Quigg R. J., et al. (2008). miR-17-92 Cluster Accelerates Adipocyte Differentiation by Negatively Regulating Tumor-Suppressor Rb2/p130. Proc. Natl. Acad. Sci. U.S.A. 105, 2889–2894. doi:10.1073/pnas.0800178105
Wang Y., Botolin D., Xu J., Christian B., Mitchell E., Jayaprakasam B., et al. (2006). Regulation of Hepatic Fatty Acid Elongase and Desaturase Expression in Diabetes and Obesity. J. Lipid Res. 47, 2028–2041. doi:10.1194/jlr.m600177-jlr200
Wang Y., Xue S., Liu X., Liu H., Hu T., Qiu X., et al. (2016). Analyses of Long Non-coding RNA and mRNA Profiling Using RNA Sequencing during the Pre-implantation Phases in Pig Endometrium. Sci. Rep. 6, 20238. doi:10.1038/srep20238
Xie F., Xiao P., Chen D., Xu L., Zhang B. (2012). miRDeepFinder: a miRNA Analysis Tool for Deep Sequencing of Plant Small RNAs. Plant Mol. Biol. 80, 75–84. doi:10.1007/s11103-012-9885-2
Xiong Y., Wang Y., Xu Q., Li A., Yue Y., Ma Y., et al. (2021). LKB1 Regulates Goat Intramuscular Adipogenesis through Focal Adhesion Pathway. Front. Physiol. 12, 755598. doi:10.3389/fphys.2021.755598
Xiong Y., Xu Q., Lin S., Wang Y., Lin Y., Zhu J. (2018). Knockdown of LXRα Inhibits Goat Intramuscular Preadipocyte Differentiation. Ijms 19, 3037. doi:10.3390/ijms19103037
Xu Q., Lin S., Wang Y., Zhu J., Lin Y. (2018a). Fibroblast Growth Factor 10 (FGF10) Promotes the Adipogenesis of Intramuscular Preadipocytes in Goat. Mol. Biol. Rep. 45, 1881–1888. doi:10.1007/s11033-018-4334-1
Xu Q., Lin S., Zhu J. J., Wang Y., Lin Y. Q. (2018b). The Expression Stability Analysis of Reference in the Process of Goat Intramuscular Preadipocytes Differentiation in Goat. Acta Veterinaria Zootechnica Sinica 495, 907–918. doi:10.11843/j.issn.0366-6964
Yang L., Yi K., Wang H., Zhao Y., Xi M. (2016). Comprehensive Analysis of lncRNAs Microarray Profile and mRNA-lncRNA Co-expression in Oncogenic HPV-Positive Cervical Cancer Cell Lines. Oncotarget 7, 49917–49929. doi:10.18632/oncotarget.10232
Yonekura S., Hirota S., Miyazaki H., Tokutake Y. (2016). Subcellular Localization and Polymorphism of Bovine FABP4 in Bovine Intramuscular Adipocytes. Anim. Biotechnol. 27, 96–103. doi:10.1080/10495398.2015.1102148
You L. H., Zhu L. J., Yang L., Shi C. M., Pang L. X., Zhang J., et al. (2015). Transcriptome Analysis Reveals the Potential Contribution of Long Noncoding RNAs to Brown Adipocyte Differentiation. Mol. Genet. Genomics 290, 1659–1671. doi:10.1007/s00438-015-1026-6
Young M. D., Wakefield M. J., Smyth G. K., Oshlack A. (2010). Gene Ontology Analysis for RNA-Seq: Accounting for Selection Bias. Genome Biol. 11, R14. doi:10.1186/gb-2010-11-2-r14
Zappaterra M., Deserti M., Mazza R., Braglia S., Zambonelli P., Davoli R. (2016). A Gene and Protein Expression Study on Four Porcine Genes Related to Intramuscular Fat Deposition. Meat Sci. 121, 27–32. doi:10.1016/j.meatsci.2016.05.007
Zhang Z., Meng Y., Gao F., Xiao Y., Zheng Y., Wang H.-Q., et al. (2020). TGF-β1-Mediated FDNCR1 Regulates Porcine Preadipocyte Differentiation via the TGF-β Signaling Pathway. Animals 10, 1399. doi:10.3390/ani10081399
Zhao X.-Y., Lin J. D. (2015). Long Noncoding RNAs: A New Regulatory Code in Metabolic Control. Trends Biochem. Sci. 40, 586–596. doi:10.1016/j.tibs.2015.08.002
Keywords: RNA-Seq, intramuscular adipocyte, subcutaneous adipocyte, goat, lncRNA, mRNA
Citation: He C, Wang Y, Zhu J, Li Y, Chen J and Lin Y (2022) Integrative Analysis of lncRNA-miRNA-mRNA Regulatory Network Reveals the Key lncRNAs Implicated Potentially in the Differentiation of Adipocyte in Goats. Front. Physiol. 13:900179. doi: 10.3389/fphys.2022.900179
Received: 20 March 2022; Accepted: 20 April 2022;
Published: 05 May 2022.
Edited by:
Dmitri Samovski, Washington University in St. Louis, United StatesReviewed by:
Zhuanjian Li, Henan Agricultural University, ChinaLiang Guo, Shanghai University of Sport, China
Copyright © 2022 He, Wang, Zhu, Li, Chen and Lin. 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: Yaqiu Lin, linyq1999@163.com