Abstract
Lipid biosynthesis is a complex process, which is regulated by multiple factors including lncRNA. However, the role of lncRNA in chicken abdominal fat accumulation is still unclear. In this research, we collected liver tissues from six high abdominal fat rate Sanhuang broilers and six low abdominal fat rate Sanhuang broilers to perform lncRNA sequencing and small RNA sequencing. A total of 2,265 lncRNAs, 245 miRNAs, and 5,315 mRNAs were differently expressed. Among of them, 1,136 differently expressed genes were enriched in the metabolic process. A total of 36 differently expressed genes, which were considered as differently expressed lncRNAs’ targets, were enriched in the metabolic process. In addition, we also found out that eight differently expressed miRNAs could target 19 differently expressed genes. FNIP2 and PEX5L were shared in a cis-regulatory network and a differently expressed miRNA target relationship network. LncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as a potential candidate that may participate in lipid synthesis. Experimentally, the objective reality of lncRNA-FNIP2/miR-24-3p/FNIP2 axis was clarified and the regulation effect of lncRNA-FNIP2/miR-24-3p/FNIP2 axis on synthesis was validated. In brief, our study reveals a potential novel regulatory mechanism that lncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as being involved in lipid synthesis during chicken adipogenesis in liver.
Background
Adipose tissue is a signal hub that regulates systemic metabolism through paracrine and endocrine signals. The liver is an important organ for the synthesis and accumulation of fat and has a regulatory role in lipid metabolism (). During the past few decades, intensive genetic selection breeding in commercial broilers has gained satisfactory achievements in increased growth rates and egg production, but also led to undesirable fat deposition (; ). Excessive fat deposition in broilers has thus become a major problem in today’s broiler industry, with the reason that it causes a reduction on the carcass yield and feed efficiency, and loss of consumers’ favor for the meat as well (). It has been proven that fat deposition is regulated by a variety of signal pathways and gene networks, among which, non-coding RNAs are the main regulators (; ; ). A better understanding of the mechanism underlying adipogenesis in chicken will bring benefits to the sustainable development of poultry industry. Whole-transcriptome sequencing has enabled the identification and characterization of a variety of putative relevant non-coding RNAs (; ).
Long non-coding RNAs (lncRNAs) are a class of non-coding RNA with a length longer than 200 nucleotides (; ). MicroRNAs (miRNAs) are a class of small non-coding RNA of approximately 22 nucleotides in length (; ). It is well established that lncRNA can act as a miRNA sponge through competing endogenous RNA (ceRNA) activity, thereby regulating the gene expression of miRNA (; ). Fat acid could be synthesized by hepatic de novo lipogenesis (DNL) from carbohydrates (especially glucose) and the liver is the main site for fatty acid metabolism and triglyceride synthesis in vertebrates (; ; ; ; ), but the regulatory mechanism of adipogenesis in chicken liver is still unclear. Investigating the RNA interaction will shed light on the gene regulatory network in the progression of chicken adipogenesis. To achieve this goal, we carried out an extensive profiling of transcriptome in the liver tissues from high-fat and low-fat chickens by RNA-Seq, screening out the differently expressed lncRNA (DEL), miRNA (DES), and mRNA (DEM).
In this study, we obtained 2,265 DELs, 245 DESs, and 5,315 mRNAs, many of which were involved in lipid metabolism. In addition, we identified multiple candidate genes and regulatory networks, including the lncRNA-FNIP2/miR-24-3p/FNIP2 axis, which was considered as being involved in lipid biosynthesis. In total, our data not only gained novel insights into regulatory mechanism in chicken adipogenesis but also contributed to a better understanding of epigenetics.
Materials and Methods
Ethics Statement
All animals used in this research were sourced from Wens. All animal experiments were conducted with the guarantee that the animals suffer minimal pain. The animal experiments were approved by the Animal Care Committee of South China Agricultural University (Approval ID: SCAU#2017015; September 13, 2017).
Samples for RNA Sequencing
After birth, 12 broilers (Sanghuang chicken) were placed in a cage with free feeding. At the 100th day, the 12 broilers were euthanized by cervical dislocation. The abdominal fat tissues were isolated and weighed. According to their abdominal fat ratio, they were divided into two groups, high-fat group and low-fat group. The abdominal fat ratio in the low-fat group ranged from 2.81 to 3.89%. The abdominal fat ratio in the high-fat group ranged from 8.04 to 9.75%. The liver tissues from the two groups were collected for RNA sequencing.
Total RNA Extraction, cDNA Library Construction, and RNA Sequencing
Total RNA was extracted by using the Trizol Reagent (Invitrogen, Carlsbad, CA, United States) following its manufacturer’s protocol. RNA integrity was detected in agarose gel electrophoresis and RNA concentration and purity were determined by the NanoDrop 2100 (Thermo Fisher Scientific, Fremont, CA, United States). Ribo-ZeroTM rRNA Removal Kit (Epicentre, Madison, Wisconsin, United States) was used to remove the ribosomal RNA (rRNA). After purification, the RNA was fragmented. Subsequently, cDNA first strand was synthesized by using the TruSeq® Stranded kit (Illumina, San Diego, CA, United States) according its protocol. DNase I (Thermo Fisher Scientific, Fremont, CA, United States) and RNase H (Thermo Fisher Scientific, Fremont, CA, United States) were used in the synthesis of the second strand of cDNA. The poly A tail and adapter were attached into the double-strand cDNA. After amplification and purification, cDNA library was obtained for RNA sequencing. BGI was responsible for RNA sequencing. The raw data obtained from sequencing were filtered to obtain clean data. The clean data was mapped to the reference genome (GRCg6a) by using HISAT (). StringTie () was used in the transcript assembly. For novel transcripts, they were divided into mRNA and lncRNA according to their coding ability, which were predicted by the CPC (), txCdsPredict (), CNC (), and pfam database (). CPC score < 0, txCdsPredict score < 500, CNC score < 0, and out of pfam database were considered as the evaluation standard for lncRNA. When at least three of the four methods were consistent, the transcript could be identified as an mRNA or lncRNA. The clean reads were aligned to reference sequences by using Bowtie2 (). The expression of genes or transcripts was calculated by RSEM (). The filtering conditions for the significant differently expressed transcripts were set as |Fold Change| ≥ 2 and Q-value ≤ 0.001. Cluster analysis for differently expressed transcripts was performed by pheatmap 1.0.8. The raw data obtained from RNA sequencing was submitted in the SRA database (accession link1).
Function Annotation of mRNA
The mRNA (including novel and known mRNA) was annotated in the NR and KEGG databases. Blast 2.2.23 () and Diamond 0.8.31 () were used in the mRNA NR and KEGG annotation. Blast2GO 2.5.0 () and NR annotation were used in the Gene Ontology (GO) annotation. The parameters used in the annotation were all default parameters.
LncRNA Target Gene Prediction and Cis-regulated Networks Construction
Cis-regulation or trans regulation on target genes were considered as the ways for lncRNA function. For cis regulation, lncRNA function was associated with its adjacent coding gene which could be regarded as its target gene. If lncRNA was located within the upstream 10 kb or downstream 20 kb from its adjacent mRNA, it would be judged as a cis regulation. If the lncRNA was out of the range and MEF (between lncRNA and mRNA) < -30 kcal/mol, it would be judged as a trans regulation. The Spearman correlation ≥ 0.6 and Pearson correlation ≥ 0.6 between lncRNA and mRNA were taken as the criteria for their target relationship. The intersection of differentially expressed mRNA and predicted target mRNA was selected to construct the cis-regulated networks between DELs and potential targets. Cytoscape 3.7.2 was used to construct the cis-regulated networks.
Small RNA Library Construction and Sequencing
The RNA sized 18–30 nt was isolated from total RNA by agarose gel electrophoresis, then their 3′ region with 5-adenylated and 3-blocked single strand DNA adapter were ligated. The RNA 5′ region was ligated with another adapter. The RT primer with Unique Molecular Index (UMI), hybridized with the 3′ adapter, was used to synthesize the first strand of cDNA. After amplification, the PCR productions sized 100–120 bp were isolated for quality examination. The qualified library was used in small RNA sequencing. The raw data obtained from small RNA sequencing were submitted on the SRA database (accession link2). After removing the none-insertion fragment sequences, long-insertion fragments sequences, low-quality sequences, polyA sequences, and small-fragment sequences, the clean data was mapped to the reference genome (GRCg6a) by using the AASRA software (). All unique RNA was sorted into a different kind RNA by annotating in the small RNA database. The priority for annotation was MiRbase > pirnabank > snoRNA > Rfam. Novel miRNAs were predicted by using miRDeep23 (). DEGseq4 () was used to analyze the differentially expressed miRNA. Significantly differentially expressed miRNAs with |Fold Change| ≥ 2 and Q-value ≤ 0.001 were considered.
Mirna Target Gene Prediction and DES-DEM Target Network Construction
Target gene prediction was performed by using miRanda (accession link5) () and RNAhybrid (accession link6) (). For the identification of potential target genes, RNAhybrid ≤ -30 kcal/mol, MEF miRanda MEF ≤ -45 kcal/mol, and miRanda score ≥ 300 were considered as the filter criteria. These selected potential targets were compared with the DEMs in the above RNA sequencing and differently expressed potential targets were selected for DES-DEM target network construction. Cytoscape 3.7.2 was used for network construction.
Reverse Transcription Reaction and Quantified Real Time PCR (qRT-PCR)
Reverse transcription reaction was performed by using the HiScript® II Q RT SuperMix for qPCR (+gDNA wiper) (Vazyme, Nanjing, Jiangsu, China) and miRNA 1st Strand cDNA Synthesis Kit (by stem-loop) (Vazyme, Nanjing, Jiangsu, China) following their manufacturers’ protocol, respectively. Quantified real time PCR (qRT-PCR) was carried out in the ABI QuantStudio 5 instrument (Thermo Fisher Scientific, Fremont, NY, United States) by using the ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, Jiangsu, China). U6 and GAPDH were respectively considered as the reference genes in qRT-PCR for miRNA, mRNA, and lncRNA expressions. The primers used in miRNA qRT-PCR were designed and synthesized by RioBio (Guangzhou, Guangdong, China). The primers used in lncRNA and mRNA qRT-PCR were designed by the Premier Primer 5.0 software and synthesized by TSINGKE Biotech (Guangzhou, Guangdong, China). The primers we designed were listed in Supplementary Table 1.
Oligonucleotides Synthesis and Plasmid Construction
The oligonucleotides used in this study (including miR-24-3p mimic and miR-24-3p inhibitor) were synthesized by RioBio (Guangzhou, Guangdong, China) and they were listed in Supplementary Table 1. The wild type sequences and mutated type sequences of lncRNA-FNIP2 and FNIP2 3′UTR with miR-24-3p binding sites were synthesized and cloned into pmiR-GLO vector by TSINGKE Biotech (Guangzhou, Guangdong, China). The mutated type sequences were transformed from wild type sequences with the binding sites which mutated from “CUGAGCU” to “GCACAUC” or from “CUGAGCU” to “GCACAUC.”
Cell Culture and Cell Transfection
Immortalized chicken preadipocytes 1 (ICP1) was provided by the Hui Li research team from Northeast Agricultural University. ICP1 was cultured in DMEM/F12 medium (Gibco, Grand Island, NY, United States) with 10% fetal bovine serum (Gibco, Grand Island, NY, United States) and 0.1%/0.2% penicillin/streptomycin (Gibco, Grand Island, NY, United States) in an incubator with 5% CO2 at 37°C. Lipofectamine 3000 reagent (Invitrogen, Carlsbad, CA, United States) was used to perform cell transfection by following its manufacturer’s protocol. The transfection dose of DNA was 1 μg/well or 0.1 μg/well for a 12-well plate or a 96-well plate, respectively. For oligonucleotides, the final transfection concentration was 50 μM.
Dual-Luciferase Reporter Assay
Wild sequences (WT) and mutational sequences (MT) of lnc-FNIP2 and FNIP2 3′UTR were respectively cloned into the pmiR-GLO vector. The transfection treatment combinations in the dual-luciferase reporter assay were respectively “WT + mimic,” “MT + mimic,” “WT + mimic NC,” and “MT + mimic NC.” Dual-luciferase reporter assay was performed by using a dual-luciferase reporter assay kit (Vazyme, Nanjing, Jiangsu, China) by following its manufacturer’s protocol. The firefly luciferase and Renilla luminescence activities were measured in a multi-function microplate reader (Biotek, Winooski, VT, United States).
Oil Red O Staining
The ICP1 cells were induced to differentiate by using the culture medium with 15% fetal serum and 0.2% oleic acid (Sigma, St Louis, CA, United States) after 12-h transfection. After 48-h induction, oil red O staining was performed by using an Oil Red Staining kit (Solarbio, Beijing, China) according to its protocol. The stained cells were captured in an electric microscope (Nikon, Tokyo, Japan). The oil red dye was extracted from the cells by using isopropanol solution. At last, the dye was quantified by a microreader (Biotek, Winooski, VT, United States).
Western Blot
The ice-cold radio immunoprecipitation (RIPA) lysis buffer (Beyotime, Shanghai, China) with 1 mM phenylmethyl sulfonyl fluoride (Beyotime, Shanghai, China) was used in the total protein extracted from liver tissues or ICP1 cells. 10% SDS-page gel was used to separate proteins and the separated proteins were transferred to the polyvinylidene fluoride (PVDF) membrane (Bio-Rad, Hercules, CA, United States). After 30-min of blocking, the membrane was incubated with anti-PPARγ (1:1,000; bs-0530R, BIOSS), anti-FNIP2 (1:500; bs-13194R, BIOSS), or anti-GAPDH (1:2,000; bsm-33033M, BIOSS) at 4°C for 12 h. Anti-mouse secondary antibody (1:10,000; 7076P2, CST) or anti-rabbit secondary antibody (1:10,000; 7074P2, CST) was used to incubate the membranes. ECL Peroxidase Color Development Kit (Vazyme, Nanjing, Jiangsu, China) was used in chromogenic reaction by following the manufacturer’s protocol. The protein bands visualization was performed in the Odyssey instrument (Li-cor, Lincoln, NE, United States).
Data Analysis
The data was represented as mean ± SEM. The statistical significance of differences between groups were evaluated by the Student’s t-test. ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001; ****P < 0.0001; ns, no significance.
Results
Blast Analysis of Transcriptome Sequencing
According to abdominal fat ratio, twelve 100-day-old Sanhuang broilers were divided in high-fat and low-fat groups. The liver tissues from the 12 broilers were collected for whole transcriptome sequencing. We obtained 12.86 Gb of raw data for each sample on average (Supplementary Table 2) and they were submitted to the NCBI database (accession ID: PRJNA684949). After filtrating the raw reads, 12.54 Gb of total clean reads were obtained. The average GC content reached 47.05% and the base proportion reaching Q30 standard was not less than 92.41%. The clean reads were mapped to the reference genome (GRCg6a), average mapping ratio was 89.57%, and unique mapping ratio was not less than 83.08% (Supplementary Table 3).
In order to distinguish mRNA and lncRNA, the coding capacity of transcripts was predicted by three softwares and the pfam database. In this sequencing, a total of 52,306 transcripts were detected and 25,760 of them were novel transcripts, including 9,861 novel lncRNAs (Supplementary Figure 1a), 12,978 novel mRNAs (Supplementary Figure 1b), 5,106 known lncRNAs, and 25,433 known mRNAs. Here, we showed the distribution of RNA length (including the novel and the known). The main length distribution of lncRNA is from 500 to 1,500 bp, while the mRNA is from 1,000 to 3,000 bp (Supplementary Figure 1c). Bowtie2 was used to align clean reads to the reference sequences, and RESM was then used to calculate gene and transcript expression. The number statistics of isoforms (lncRNA and mRNA) for 12 samples was listed in Supplementary Table 4. In order to visualize the number of genes in each sample in different FPKM intervals, we counted the number of genes in three conditions of FPKM (FPKM ≤ 1, 1 < FPKM < 10, FPKM ≥ 10). The number of transcripts which FPKM ≥ 10 was not less than 3,206 (Supplementary Figure 1d).
DEG Between Liver Tissues From High-Fat to Low-Fat Chickens
To identify the potential candidate genes related to adipogenesis, the expression levels of lncRNA and mRNA were examined in liver tissues from high-fat and low-fat chickens. Q-value ≤ 0.001 and | fold change| ≥ 2 were set as the standard for differential expression. A total of 2,265 lncRNAs (Supplementary Table 5) and 5,315 mRNAs (Supplementary Table 7) were differently expressed between liver tissues in high-fat and low-fat chickens. The differently expressed lncRNAs (DELs) and differently expressed mRNAs (DEMs) were clustered (Figures 1A,B) by the pheatmap software, according to their fold change. The number of DELs (Supplementary Table 5), differently expressed genes (DEGs) (Supplementary Table 6), and DEMs (Supplementary Table 7) between the different groups are shown in Figure 2A. In these 2,265 DELs, including 1,821 novel lncRNAs and 444 known lncRNAs, 1,350 lncRNAs were upregulated and 915 lncRNAs were downregulated. Among the 5,315 DEMs, 3,138 and 2,177 mRNAs were upregulated and downregulated, respectively. Besides, we made statistics on the differentially expressed genes at the isoform and gene levels, respectively (Figure 2B). A total of 1,864 DEGs were upregulated and 1,275 DEGs were downregulated.
FIGURE 1
FIGURE 2
The top 10 abundant DELs and top 10 abundant DEMs were listed in Tables 1 and 2, respectively. All the top 10 abundant DELs were novel lncRNAs and only two of them were upregulated. LTCONS_00020831 (log2FC = 1.21, Q-value = 0), the most abundant DEL in this sequencing, was upregulated in liver tissues from high-fat chickens. Eight of the top 10 abundant DEMs were downregulated, and seven of the top 10 abundant DEMs were known transcripts. For DEMs, ENSGALT00000002892 (VTG2, vitellogenin 2, log2FC = –1.73992, Q-value = 0) was highly expressed in the liver of low-fat chickens and showed the highest abundance in DEMs.
TABLE 1
| LncRNA ID | Length | L Reads number | H Reads number | log2FC(H/L) | H/L | q-Value |
| LTCONS_00020831 | 1,186 | 79,374.24 | 168,343.7 | 1.213923 | Up | 0 |
| LTCONS_00020832 | 1,351 | 230,312.4 | 5,371.91 | –5.29276 | Down | 0 |
| LTCONS_00055244 | 43,793 | 113,243 | 49,591 | –1.06202 | Down | 0 |
| LTCONS_00029013 | 512 | 158,164.8 | 4,411.57 | –5.03474 | Down | 0 |
| LTCONS_00040763 | 2,932 | 73,239.73 | 30,631.24 | –1.12837 | Down | 0 |
| LTCONS_00046846 | 24,766 | 77,193.92 | 3,376.17 | –4.38577 | Down | 0 |
| LTCONS_00000249 | 10,649 | 50,129.79 | 19,380 | –1.24184 | Down | 0 |
| LTCONS_00025768 | 41,916 | 47,501.63 | 20,827.96 | –1.0602 | Down | 0 |
| LTCONS_00010452 | 68,068 | 21,954.88 | 40,581.95 | 1.015552 | Up | 0 |
| LTCONS_00000621 | 12,396 | 45,945.24 | 12,643.59 | –1.73225 | Down | 0 |
The top 10 abundant DELs between the high-fat and low-fat groups.
TABLE 2
| mRNA ID | Length | L Reads number | H Reads number | log2FC (H/L) | H/L | q-Value |
| ENSGALT00000002892 | 5,565 | 2,165,799 | 592,843.8 | –1.73992 | Down | 0 |
| ENSGALT00000056823 | 5,951 | 959,658 | 234,419 | –1.90418 | Down | 0 |
| MTCONS_00058597 | 6,098 | 282,406.2 | 119,824.5 | –1.10759 | Down | 0 |
| ENSGALT00000083486 | 1,927 | 371,304 | 112,513 | –1.59325 | Down | 0 |
| ENSGALT00000027531 | 4,885 | 24,398.34 | 84,566.61 | 1.922561 | Up | 0 |
| ENSGALT00000019639 | 1,603 | 168,144 | 70,319 | –1.12846 | Down | 0 |
| ENSGALT00000068577 | 805 | 26,347 | 59,925.02 | 1.314775 | Up | 0 |
| MTCONS_00056041 | 5,126 | 545,124.6 | 59,193.44 | –3.07382 | Down | 0 |
| ENSGALT00000024496 | 2,764 | 159,362 | 58,913 | –1.30639 | Down | 0 |
| MTCONS_00043232 | 2,325 | 132,570.9 | 58,860.85 | –1.04213 | Down | 0 |
The top 10 abundant DEMs between the high-fat and low-fat groups.
Functional Enrichment Analysis of DEMs Involved in Lipid Metabolism
In order to explore the functions of DEMs involved in chicken lipid metabolism, GO enrichment analysis was performed in this study (Supplementary Table 8). Here, we only analyze the DEGs corresponding to DEMs. As shown in Figure 3A, 1,136 DEGs were enriched in the metabolic process (650 upregulated and 486 downregulated DEGs), which may be associated with chicken abdominal fat deposition. In addition, 65 upregulated and 61 downregulated DEGs were enriched in cell proliferation. To determine in which pathways the DEGs between high-fat and low-fat chickens are more concentrated, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed. Q-value < 1 was set as the cut off for the significantly enriched pathways. The top 20 enriched KEGG pathways were shown in Figure 3B. Our result suggested that the DEGs were significantly enriched in Cell adhesion molecules (CAMs), Graft-versus-host disease, Type I diabetes mellitus, and so on. Here, we showed all the DEGs enriched pathways (Supplementary Table 9). A total of 494 DEGs were significantly enriched in the Metabolic pathways, which may be involved with lipid metabolism during chicken abdominal fat accumulation.
FIGURE 3
Target Gene Prediction of DEL
The function of lncRNA is mainly achieved by acting on the target genes in cis mode. The basic principle of cis target gene prediction is that the function of lncRNA is related to the protein coding genes adjacent to its coordinates, so the mRNA adjacent to lncRNA is screened out as its target gene. Here, we totally obtained 9,242 target pairs between lncRNA and their adjacent mRNA (Supplementary Table 10). A total of 2,634 of these target pairs were localized at the 10K range upstream of the mRNAs, and the target pair number of the distance between the upstream lncRNA and the mRNA ranging from 1 to 1,000 nt was 796. On the other hand, 4,003 target pairs were localized at the 20K range downstream of the mRNAs, and the distance of 513 pairs were within 1,000 nt from downstream lncRNA to mRNA. In addition, there was an overlap between 1,820 lncRNA and 2,075 mRNA, and 2,512 target pairs were formed (Figure 4A and Supplementary Table 11). There were 628 lncRNAs that overlapped with 774 mRNAs, forming 907 target pairs. In addition, 591 lncRNAs anti-overlapped with 724 mRNAs forming 833 target pairs. Spearman-correlation ≥ 0.6 and Pearson-correlation ≥ 0.6 between lncRNA and mRNA were considered as the cut off for DELs’ target genes, and we obtained 778 DELs’ target genes. The 778 target genes were compared with DEMs, 184 mRNAs of which were differently expressed (Figure 4B).
FIGURE 4
Functional Enrichment Analysis of Differently Expressed DELs’ Target Genes
To evaluate the biological function of these differently expressed DELs’ target genes, we carried out the functional enrichment analysis. GO function enrichment analysis showed that 36 DEMs were enriched in the metabolic process (Figure 5A). The 36 DEMs might participate in lipid metabolism and chicken abdominal fat deposition, regulated by lncRNAs to achieve this biological process. In addition, there were 29, 14, and seven DEMs, respectively, enriched in the regulation of biological process, developmental process, and cell proliferation, which may be involved in abdominal fat adipocytes proliferation and differentiation. KEGG enrichment analysis was also performed to determine which pathways were more concentrated in which target genes were significantly different between groups. Q-value < 0.05 was considered as the criteria for the significant differences and the DEMs could be enriched in 142 pathways (Supplementary Table 12). Here, we showed the top 20 significantly enriched pathways in Figure 5B. The most enriched pathways are the Osteoclast differentiation, ECM-receptor interaction, Platelet activation, Graft-versus-host disease, and Type I diabetes mellitus. The 23 DEMs, enriched in the Metabolic pathways, were considered to associate with the liver regulation of abdominal fat deposition (Supplementary Table 12).
FIGURE 5
DEL-DEG Interaction Network Construction
To further evaluate how DELs come into play in regulating their targets, we carried out the DEL-DEG interaction network construction. Here, we found that the 184 DEMs, which we obtained above, were transcribed from 68 genes and we obtained 99 target pairs between the 68 DEGs and 97 DELs (Supplementary Table 13). We constructed two target interaction networks between the upregulated DELs and upregulated DEGs or downregulated DELs and downregulated DEGs, respectively (Figure 6). A total of 45 target pairs were formed by 45 upregulated DELs and 32 upregulated DEGs (Figure 6A). In this network, lncRNA-GBE1 (LTCONS_00002483, log2FC = 1.07), lncRNA-PEX5L (LTCONS_00058214, log2FC = 1.34), lncRNA-PARD3 (LTCONS_00021300, log2FC = 1.22), lncRNA-NTNG1 (LTCONS_00056507, log2FC = 1.14), and lncRNA-FNIP2 (LTCONS_00043262, log2FC = 1.20) were the top five abundant lncRNAs, which could respectively target to GBE1 (log2FC = 2.60), PEX5L (log2FC = 1.20), PARD3 (log2FC = 1.67), NTNG1 (log2FC = 1.44), and FNIP2 (log2FC = 1.01), indicating their potential cis regulatory relationship. On the other hand, 54 target pairs were formed by 52 downregulated DELs and 36 downregulated DEGs (Figure 6B). In the downregulated DEL-downregulated DEG interaction network, the top five abundant lncRNAs were lncRNA-SPIA3 (LTCONS_00050642, log2FC = –3.71), lncRNA-SLC38A2 (LTCONS_00005023, log2FC = –1.47), lncRNA-IGF1 (LTCONS_00001244, log2FC = –1.63), lncRNA-JPH2 (LTCONS_00060413, log2FC = –1.25), and lncRNA-SOX7 (LTCONS_00037430, log2FC = –2.29). In addition, in our previous sequencing, LPIN1 was downregulated in the abdominal fat tissue from the high-fat Sanhuang chicken, which, in this sequence, its downregulation was also found. This network revealed a novel potential regulatory mechanism that a low expression of lncRNA-LPIN1 (LTCONS_00039704) mediates the down-regulation of LPIN1 expression and, thus, affects lipid metabolism.
FIGURE 6
MiRNA Expression Profile
In order to further clarify the potential regulatory mechanism of abdominal fat deposition in an epigenetic perspective, the RNA, used in lncRNA sequencing, was used in a small RNA sequencing on the BG1SEQ-500 platform. In total, we obtained an average of 24.24 Gb of raw data for each sample and the raw data were submitted to the SRA database (accession ID: PRJNA686699). After filtering, an average of 21.73 Gb of clean data for each sample was obtained and the percentage of clean tag reached 89.63% (Supplementary Table 14). The clean tag was mapped to the known small RNA databases, including miRBase, Rfam, siRNA, piRNA, and snoRNA, and the mapped rate ranged from 83.62 to 94.25% (Supplementary Table 14). In this sequencing, we detected 1,100 microRNAs (miRNAs), including 545 known miRNAs and 555 novel miRNAs. | Fold change| ≥ 2 and q-Value ≤ 0.001 were considered as a significant difference of miRNA expression. In total, 245 DESs were found, including 137 novel miRNAs and 108 known miRNAs (Figure 7A). Among them, 82 miRNAs were upregulated and 163 miRNAs were downregulated (Figure 7B). We listed the top 20 abundant differently expressed miRNAs in Table 3. There were nine and 11 miRNAs that were upregulated and downregulated, respectively. The most abundant DES was miR-92a_1 (FC = -1.78953), which was downregulated in the liver of chickens with a high abdominal fat rate.
FIGURE 7
TABLE 3
| miRNA id | Expression (L) | Expression (H) | log2Ratio (H/L) | Up/down | q-Value |
| miR-92a_1 | 35,074 | 7,418 | –1.78953 | DOWN | 0 |
| miR-191-5p | 27,280.33 | 9,062.833 | –1.13806 | DOWN | 0 |
| miR-100-5p_1 | 10,997.17 | 24,209 | 1.590181 | UP | 0 |
| miR-126-3p | 22,141 | 2,196.833 | –2.88145 | DOWN | 0 |
| miR-26c_1 | 412.5 | 21,795.17 | 6.17524 | UP | 0 |
| miR-126-3p_1 | 2,961 | 17,333.33 | 3.001163 | UP | 0 |
| miR-92a_2 | 1,384.167 | 17,077 | 4.076735 | UP | 0 |
| novel_mir182 | 7,226.5 | 303.3333 | –4.12255 | DOWN | 0 |
| miR-429 | 2,103 | 4,070.5 | 1.404527 | UP | 0 |
| miR-125b_1 | 2,115.667 | 3,120.167 | 1.012281 | UP | 0 |
| miR-10a-5p | 4,031 | 995 | –1.5666 | DOWN | 0 |
| miR-24-3p | 3,729.667 | 1,236.5 | –1.14101 | DOWN | 0 |
| novel_mir4 | 3,346.167 | 1,189.333 | –1.04059 | DOWN | 0 |
| let-7b | 3,158.333 | 1,033.833 | –1.15939 | DOWN | 0 |
| miR-146a_1 | 3,093.333 | 716.6667 | –1.65802 | DOWN | 0 |
| miR-140-5p_1 | 1,253.667 | 2,049.5 | 1.160888 | UP | 0 |
| miR-146a-5p_1 | 1,103.667 | 2,170.167 | 1.427271 | UP | 0 |
| miR-429_1 | 3,218.167 | 42.16667 | –5.80222 | DOWN | 0 |
| miR-10a_2 | 638.8333 | 2,542.5 | 2.444506 | UP | 0 |
| miR-103_1 | 1,747.5 | 602 | –1.08569 | DOWN | 0 |
The top 20 abundant differently expressed miRNAs.
The Target Prediction of Differently Expressed miRNAs
To evaluate the regulatory mechanism of these differently expressed miRNAs, target prediction was performed by RNAhybrid and miRanda. The filter parameters of screening potential target mRNA were as follows: RNAhybrid MEF ≤ -30 kcal/mol, miRanda MEF ≤ -45 kcal/mol, and miRanda score ≥ 300. In total, 1,407 target pairs were formed by 32 known miRNAs and 1,068 mRNAs. These potential mRNAs were compared to the DEGs we obtained in the above whole transcriptome sequencing, and we found 23 target pairs between 19 DEGs and 8 DESs, including 13 negative target pairs and 10 positive target pairs. Based on the circumstances of regulation on mRNA caused by miRNA, we performed the target relationship network construction (Figure 8). In this network, we could find that miR-24-3p was at a critical node position, interacted with six DEGs, including three negative interaction and three positive interaction. miR-17-3p_2 had four positive targets, three of them were negatively regulated by miR-17_1.
FIGURE 8
The regulatory networks between DELs and DEGs were compared with the regulatory network between DESs and DEGs, in which FNIP2 and PEX5L was found out to be common to both networks. The low expression of miR-15a may be responsible for the upregulation of PEX5L which may also be cis regulated by lncRNA-PEX5L. Similar in PEX5L, the upregulation of FNIP2 may be caused by miR-24-3p low expression and lncRNA-FNIP2 high expression. Interestingly, it was found that lncRNA-FNIP2 has the potential to target and bind miR-24-3p (Supplementary Figure 2), suggesting that a miR-24-3p low expression may be caused by lncRNA-FNIP2. The target relationship between lncRNA-FNIP2 and miR-24-3p means that lncRNA-FNIP2 not only affects the expression of FNIP2 through cis regulation, but also by the lncRNA-FNIP2/miR-24-3p/FNIP2 axis.
FNIP2 Promotes Preadipocyte Lipid Synthesis
In order to validate the expression difference of FNIP2 between high-fat and low-fat groups, qRT-PCR was performed and the result showed that the FNIP2 RNA level in the high-fat group was higher than that in the low-fat group (Figure 9A). In addition, the same trend was observed in the FNIP2 protein level (Figure 9B). The results indicate that FNIP2 was high-expressed in the liver of high-fat individuals. To verify the potential role of FNIP2 in adipogenesis, we constructed the overexpression plasmid and synthesized the specific siRNA of FNIP2. The transfection efficiency of the plasmid and siRNA in ICP1 cells was validated both in the mRNA and protein levels (Figures 9C,D). ICP1 cells were induced to differentiate and the lipid droplet formation was detected by Oil Red O Staining. It could be found that FNIP2 was able to induce more lipid droplet formation, while the knockdown of FNIP2 could suppress lipid droplet formation in the ICP1 cells (Figures 9E,F). FNIP2 could upregulate the mRNA expression of some lipid biosynthesis-related genes (including PPARγ, CEBP/α, CEBP/β, LPL, and ADIPOR1), while the knockdown of FNIP2 downregulated their mRNA expression (Figures 9G,H). Besides, PPARγ protein level was increased by FNIP2, indicating that PPARγ signal could be intensified by FNIP2 (Figure 9I). In short, these results suggest that the upregulated FNIP2 in high-fat individuals could facilitate preadipocytes lipid synthesis.
FIGURE 9
miR-24-3p Suppresses Lipid Synthesis by Targeting to FNIP2
The expression difference of miR-24-3p was detected by qRT-PCR and the expression trend was consistent with the sequencing result (Figure 10A). Based on the target prediction above between miR-24-3p and FNIP2 (Figure 10B), we carried out a dual-luciferase reporter assay to verify the binding between miR-24-3p and FNIP2. It was showed that the co-transfection between wild type FNIP2 3′UTR and miR-24-3p could reduce the luciferase activity (Figure 10C), indicating a target relationship between FNIP2 and miR-24-3p. qRT-PCR and the Western Blot results showed that miR-24-3p does have a posttranscriptional regulation on FNIP2, which miR-24-3p downregulated FNIP2 not only at the mRNA level but also at the protein level (Figures 10D,E). Considering their target relationship, we performed another co-transfection to verify whether the promotion of FNIP2 on lipid synthesis can be modulated by miR-24-3p. We found that the upregulation of PPARγ, CEBP/α, CEBP/β, LPL, and ADIPOR1 caused from FNIP2 could be restored to normal levels by miR-24-3p (Figure 10F), being in accord with our expectation. In the Oil Red O Staining, it was found that miR-24-3p expression could decrease lipid droplet formation (Figures 10G,H). In addition, miR-24-3p was able to downregulate PPARγ, CEBP/α, CEBP/β, LPL, and ADIPOR1, while miR-24-3p inhibitor upregulated their expression (Figures 10I,J). With the expression of miR-24-3p, PPARγ protein level was reduced, which manifested an inhibition of the miR-24-3p on the PPARγ signal (Figure 10K). These results suggest that miR-24-3p may cause lipid synthesis inhibition by targeting and downregulating FNIP2.
FIGURE 10
LncRNA-FNIP2 Accelerates Lipid Synthesis by Releasing FNIP2 From miR-24-3p
The expression difference of lncRNA-FNIP2 between the high-fat and low-fat groups was validated by qRT-PCR (Figure 11A). The above results indicated that lncRNA-FNIP2 may affect FNIP2 through ceRNA mechanism (by miR-24-3p) during adipogenesis, and lncRNA-FNIP2 was predicted to adsorb miR-24-3p at positions 952–958 (Figure 11B). To verify the potential target relationship between lncRNA-FNIP2 and miR-24-3p, we performed a dual-luciferase reporter assay. It could be found that the luciferase activity in the co-transfection group of wild type lncRNA-FNIP2 and miR-24-3p was significantly lower than the other groups (Figure 11C), revealing their target relationship. In restored validation, we found that the mRNA expression of PPARγ, CEBP/α, CEBP/β, LPL, and ADIPOR1, which have been reduced by miR-24-3p, could be restored or reversed by lncRNA-FNIP2 overexpression (Figure 11D). Based on the situation that both lncRNA-FNIP2 and FNIP2 have a target relationship with miR-24-3p, we tried to verify whether lncRNA-FNIP2 could release FNIP2 from miR-24-3p. It could be found that lncRNA-FNIP2 was be able to restore the FNIP2 mRNA level from miR-24-3p, while the mutation of miR-24-3p binding site on lncRNA-FNIP2 lost this modulation function (Figure 11E). Subsequently, lncRNA-FNIP2 also showed a promoting effect on FNIP2 both at the mRNA and protein levels (Figures 11F,G). Similarly, Oil Red O staining indicated an effect of lncRNA-FNIP2 on accelerating lipid droplet formation (Figures 11H,I). Additionally, lncRNA-FNIP2 increased the mRNA levels of PPARγ, CEBP/α, CEBP/β, LPL, and ADIPOR1, while lncRNA-FNIP2 knockdown decreased them (Figures 11J,K). LncRNA-FNIP2 could facilitate PPARγ protein expression, suggesting a promoting effect of lncRNA-FNIP2 on the PPARγ signal. These results reveal a positive role of lncRNA-FNIP2 on lipid synthesis, mediated by lncRNA-FNIP2/miR-24-3p/FNIP2 axis during the adipogenesis process.
FIGURE 11
Discussion
Lipid synthesis is a crucial process in animal fat accumulation and epigenetics has been characterized in vertebrate lipid metabolism and lipogenesis (; ; ). Recently, lncRNAs have been identified to be associated with adipogenesis (; ). However, the regulatory role of lncRNA on chicken abdominal fat deposition is still unclear. In this study, we tried to gain an insight of the regulatory mechanism on chicken abdominal fat deposition by performing a high throughput sequencing analysis at different levels (including lncRNA, mRNA, and miRNA). A total of 2,265 DELs, 5,315 DEMs, and 245 DESs were obtained from our sequencing analysis. We constructed two interaction networks based on the cis regulation relationship between 68 DEGs and 97 DELs. Besides, the interaction network between 19 DEGs and eight DESs was constructed. Based on these interaction networks, the potential target axis lncRNA-FNIP2/miR-24-3p/FNIP2 was hypothesized as a candidate regulatory mechanism in chicken adipogenesis.
Cross-talk between adjacent genes is a common phenomenon, involving multiple regulatory mechanisms and cis-regulatory signals (). Recent researches suggested that the activity of transcription or DNA elements in lncRNA locus may be the way for lncRNA to regulate adjacent genes in cis (). The potential of lncRNA during chicken abdominal fat adipocytes differentiation has been characterized (,). In this research, we constructed the cis-regulated networks between DEL and DEG. The potential cis-regulation relationships of five up-regulated lncRNAs and five down-regulated lncRNAs with their adjacent genes were respectively speculated, including lncRNA-GBE1/GBE1, lncRNA-PEX5L/PEX5L, lncRNA-PARD3/PARD3, lncRNA-NTNG1/NTNG1, lncRNA-FNIP2/FNIP2, lncRNA-SPIA3/SPIA3, lncRNA-SLC38A2/SLC38A2, lncRNA-IGF1/IGF1, lncRNA-JPH2/JPH2, and lncRNA-SOX7/SOX7. GBE1 (glycogen-branching enzyme 1) is required in glycogen accumulation for maintaining glucose metabolism balance (; ; ). Our previous study found that GBE1 was upregulated in the liver from the fast-growing WRR chickens (). Recently, GBE1 was suggested to be associated with backfat thickness traits (), implying that lncRNA-GBE1 may participate in the glycogen metabolic process by regulating GBE1 in cis. Parental high-fat diet increased offspring obesity and type 2 diabetes mellitus risks through the regulation of SLC28A2 (solute carrier family 38 member 2) expression (; ), which indicated that lncRNA-SLC28A2 may be associated with chicken lipid synthesis.
In our previous research, LPIN1, downregulated by miR-429 in high-fat group abdominal fat tissues (the same individuals used in this study), inhibited lipid droplet formation. For epigenetics, miRNA play a pivotal role in the chicken lipid metabolic process (; ; ). In this study, we also performed a small RNA sequencing and obtained 245 DESs. Compared with the sequencing data in abdominal fat tissues, relatively few miRNAs were detected in the liver and we did not observe miR-429-3p, indicating that it may be expressed ectopic in abdominal fat. Here, we hypothesized a potential mechanism that LPIN1 expression abnormity may be caused by the transcriptional activity reduction of the upstream lncRNA-LPIN1 gene locus. We constructed a target network between DEMs and DESs, which was formed by 13 negative target pairs and 10 positive target pairs. Two genes, including FNIP2 and PEX5L, are shared between the DESs-DEMs and DELs-DEMs networks. The cytosolic receptor, PEX5L (peroxisome biogenesis factor 5 long isoform), could recognize peroxisomal targeting signals 1 or 2 and peroxisomal import of matrix proteins would be initiated, which are involved in several metabolic processes (including lipid biosynthesis, fatty acid α-, and β-oxidation) (; ). In our study, PEX5L was upregulated by the synergistic action of low expressed miR-15a and high expressed lncRNA-PEX5L, implying their potential on chicken liver lipid biosynthesis. Similar to miR-429-3p, miR-15a and miR-24-3p did not show a different expression trend in chicken abdominal fat tissues, which also characterized their tissue specificity. The reduction of miR-24-3p may be the reason for the increase of the fat mass and leptin levels during the treatment improvement of HD (Huntington’s disease) patients (). In addition, miR-24-3p was also downregulated during the bovine preadipocytes differentiation process, which was predicted to target the FASN (fatty acid synthase) gene (). Our research also revealed that miR-24-3p low expression in liver might cause more lipid formation.
The novel lncRNA-FNIP2 was predicted to sponge miR-24-3p and dual-luciferase reporter assay indeed identified their target relationship, as well as FNIP2. Here, we verified their role on lipid synthesis and it was found that lncRNA-FNIP2 and FNIP2 were able to facilitate ICP1 lipid synthesis, while miR-24-3p inhibited lipid synthesis. Nutrient-rich conditions activate mTORC1 to trigger downstream anabolic reactions (). FNIP2 (folliculin interacting proteins 2) could form a complex with FLCN. The complex would directly contact the Rag GTPases to stimulate GTP hydrolysis to GDP-bound state, which would promote mTORC1 activation (), indicating that FNIP2 may contribute to the mTORC1 signal to accelerate lipid synthesis, but further verification is needed to clarify the underlying mechanism. In this study, we performed a restored experiment and it showed that lncRNA-FNIP2 did release FNIP2 from miR-24-3p, suggesting the objective reality of the lncRNA-FNIP2/miR-24-3p/FNIP2 axis. Besides, the promotion of the lncRNA-FNIP2/miR-24-3p/FNIP2 axis on the PPARγ signal was preliminary demonstrated, suggesting a positive effect of the lncRNA-FNIP2/miR-24-3p/FNIP2 axis on lipid metabolism. In short, we screened the lncRNA-FNIP2/miR-24-3p/FNIP2 axis that may be related to abdominal fat deposition through the whole transcriptome analysis, and verified the real existence of the lncRNA-FNIP2/miR-24-3p/FNIP2 axis and its promoting effect on lipid synthesis.
Statements
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/, PRJNA684949 and https://www.ncbi.nlm.nih.gov/, PRJNA686699.
Ethics statement
All animals used in this research were sourced from Wens. All animal experiments were conducted with the guarantee that the animals suffer minimal pain. The animal experiments were approved by the Animal Care Committee of South China Agricultural University (Approval ID: SCAU#2017015; September 13, 2017).
Author contributions
LG was responsible for the research designing, data analysis, and manuscript writing. XC participated partly in the experiments and data analysis. WH and ZL reviewed and modified the manuscript. KL, MY, SZ, ML, HL, WL, and XZ were responsible for the animal experiments and sample collection. QN participated in the research designing. QL carried out the design of the whole research and guided the research progress. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Guangdong Basic and Applied Basic Research Fund project (2019B1515210005), the Natural Science Foundation of Guangdong Province (2018A030310206), the Guangdong Province Modern Agricultural Industry Technology System Project (2020KJ128), and the Science and Technology Program of Guangzhou (201803020016).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.653798/full#supplementary-material
Supplementary Figure 1The coding capacity of novel transcripts, distribution of RNA length, and isoforms expression in each sample.
Supplementary Figure 2The target relationship prediction of miR-24-3p.
Supplementary Table 1Fragment sequence information.
Supplementary Table 2The information of sequencing data.
Supplementary Table 3Genome alignment statistical of clean reads.
Supplementary Table 4The number statistics of lncRNA and mRNA isoforms.
Supplementary Table 5The differently expressed lncRNA between high-fat and low-fat groups.
Supplementary Table 6The differently expressed genes between high-fat and low-fat groups.
Supplementary Table 7The differently expressed mRNAs between high-fat and low-fat groups.
Supplementary Table 8The GO term DEMs enriched in.
Supplementary Table 9The KEGG pathway DEMs enriched in.
Supplementary Table 10The statistics of DELs cis regulation.
Supplementary Table 11Classification statistics of DEL overlap with mRNA.
Supplementary Table 12The KEGG pathway differently expressed DELs’ target genes enriched in.
Supplementary Table 13Target pairs between DELs and DEGs.
Supplementary Table 14Small RNA sequencing data statistics.
Footnotes
1.^https://www.ncbi.nlm.nih.gov/bioproject/PRJNA684949
2.^https://www.ncbi.nlm.nih.gov/bioproject/PRJNA686699
3.^https://www.mdc-berlin.de/8551903/en/
4.^http://www.bioconductor.org/packages/release/bioc/html/DEGseq.html
References
1
AganzoM.MontojoM. T.LopezD. L. H. M.Martinez-DescalsA.Ricote-VilaM.SanzR.et al (2018). Customized dietary intervention avoids unintentional weight loss and modulates circulating miRNAs footprint in Huntington’s disease.Mol. Nutr. Food Res.62:e1800619. 10.1002/mnfr.201800619
2
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool.J. Mol. Biol.215403–410. 10.1016/S0022-2836(05)80360-2
3
Alvarez-DominguezJ. R.BaiZ.XuD.YuanB.LoK. A.YoonM. J.et al (2015). De novo reconstruction of adipose tissue transcriptomes reveals long non-coding RNA regulators of brown adipocyte development.Cell Metab.21764–776. 10.1016/j.cmet.2015.04.003
4
BaoY.KishnaniP.WuJ. Y.ChenY. T. (1996). Hepatic and neuromuscular forms of glycogen storage disease type IV caused by mutations in the same glycogen-branching enzyme gene.J. Clin. Invest.97941–948. 10.1172/JCI118517
5
BuchfinkB.XieC.HusonD. H. (2015). Fast and sensitive protein alignment using DIAMOND.Nat. Methods1259–60. 10.1038/nmeth.3176
6
ChakravarthyM. V.LodhiI. J.YinL.MalapakaR. R.XuH. E.TurkJ.et al (2009). Identification of a physiologically relevant endogenous ligand for PPARalpha in liver.Cell138476–488. 10.1016/j.cell.2009.05.036
7
ChakravarthyM. V.PanZ.ZhuY.TordjmanK.SchneiderJ. G.ColemanT.et al (2005). “New” hepatic fat activates PPARalpha to maintain glucose, lipid, and cholesterol homeostasis.Cell Metab.1309–322. 10.1016/j.cmet.2005.04.002
8
ChangB. H.LiL.PaulA.TaniguchiS.NannegariV.HeirdW. C.et al (2006). Protection against fatty liver but normal adipogenesis in mice lacking adipose differentiation-related protein.Mol. Cell. Biol.261063–1076. 10.1128/MCB.26.3.1063-1076.2006
9
ChenY.ZhaoY.JinW.LiY.ZhangY.MaX.et al (2019). MicroRNAs and their regulatory networks in Chinese Gushi chicken abdominal adipose tissue during postnatal late development.BMC Genomics20:778. 10.1186/s12864-019-6094-2
10
ChongT.XieY.WeiY. (2017). AASRA: an anchor alignment-based small RNA annotation pipeline.BioRxiv [Preprint]. 10.1101/132928
11
ClaireD. H.PaulW.ShenX.JiaX.ZhangR.SunL.et al (2013). Identification and characterization of genes that control fat deposition in chickens.J. Anim. Sci. Biotechnol.4:43. 10.1186/2049-1891-4-43
12
Claycombe-LarsonK. G.BundyA. N.RoemmichJ. N. (2020). Paternal high-fat diet and exercise regulate sperm miRNA and histone methylation to modify placental inflammation, nutrient transporter mRNA expression and fetal weight in a sex-dependent manner.J. Nutr. Biochem.81:1024. 10.1016/j.jnutbio.2020.108373
13
ConesaA.GotzS.Garcia-GomezJ. M.TerolJ.TalonM.RoblesM. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research.Bioinformatics213674–3676. 10.1093/bioinformatics/bti610
14
EngreitzJ. M.HainesJ. E.PerezE. M.MunsonG.ChenJ.KaneM.et al (2016). Local regulation of gene expression by lncRNA promoters, transcription and splicing.Nature539452–455. 10.1038/nature20149
15
FinnR. D.CoggillP.EberhardtR. Y.EddyS. R.MistryJ.MitchellA. L.et al (2016). The Pfam protein families database: towards a more sustainable future.Nucleic Acids Res.44D279–D285. 10.1093/nar/gkv1344
16
Friedl NderM. R.ChenW.AdamidiC.MaaskolaJ.EinspanierR.KnespelS.et al (2008). Discovering microRNAs from deep sequencing data using miRDeep.Nat. Biotechnol.26407–415.
17
FroeseD. S.MichaeliA.McCorvieT. J.KrojerT.SasiM.MelaevE.et al (2015). Structural basis of glycogen branching enzyme deficiency and pharmacologic rescue by rational peptide design.Hum. Mol. Genet.245667–5676. 10.1093/hmg/ddv280
18
FuX.DongB.TianY.LefebvreP.MengZ.WangX.et al (2015). MicroRNA-26a regulates insulin sensitivity and metabolism of glucose and lipids.J. Clin. Invest.1252497–2509. 10.1172/JCI75438
19
GaoX.LinS. H.RenF.LiJ. T.ChenJ. J.YaoC. B.et al (2016). Acetate functions as an epigenetic metabolite to promote lipid synthesis under hypoxia.Nat. Commun.7:11960. 10.1038/ncomms11960
20
HahnO.GronkeS.StubbsT. M.FiczG.HendrichO.KruegerF.et al (2017). Dietary restriction protects from age-associated DNA methylation and induces epigenetic reprogramming of lipid metabolism.Genome Biol.18:56. 10.1186/s13059-017-1187-1
21
HasanS.PlattaH. W.ErdmannR. (2013). Import of proteins into the peroxisomal matrix.Front. Physiol.4:261. 10.3389/fphys.2013.00261
22
HootmanK. C.TrezziJ. P.KraemerL.BurwellL. S.DongX.GuertinK. A.et al (2017). Erythritol is a pentose-phosphate pathway metabolite and associated with adiposity gain in young adults.Proc. Natl. Acad. Sci. U.S.A.114E4233–E4240. 10.1073/pnas.1620079114
23
HudginsL. C.HellersteinM.SeidmanC.NeeseR.DiakunJ.HirschJ. (1996). Human fatty acid synthesis is stimulated by a eucaloric low fat, high carbohydrate diet.J. Clin. Invest.972081–2091. 10.1172/JCI118645
24
Jensen-UrstadA. P.SemenkovichC. F. (2012). Fatty acid synthase and liver triglyceride metabolism: housekeeper or messenger?Biochim. Biophys. Acta1821747–753. 10.1016/j.bbalip.2011.09.017
25
JohnB.EnrightA. J.AravinA.TuschlT.SanderC.MarksD. S. (2004). Human MicroRNA targets.PLoS Biol.2:e363. 10.1371/journal.pbio.0020363
26
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements.Nat. Methods12357–360. 10.1038/nmeth.3317
27
KongL.ZhangY.YeZ. Q.LiuX. Q.ZhaoS. Q.WeiL.et al (2007). CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine.Nucleic Acids Res.35W345–W349. 10.1093/nar/gkm391
28
KoppF.MendellJ. T. (2018). Functional classification and experimental dissection of long noncoding RNAs.Cell172393–407. 10.1016/j.cell.2018.01.011
29
KrolJ.LoedigeI.FilipowiczW. (2010). The widespread regulation of microRNA biogenesis, function and decay.Nat. Rev. Genet.11597–610. 10.1038/nrg2843
30
KroutD.RoemmichJ. N.BundyA.GarciaR. A.YanL.Claycombe-LarsonK. J. (2018). Paternal exercise protects mouse offspring from high-fat-diet-induced type 2 diabetes risk by increasing skeletal muscle insulin signaling.J. Nutr. Biochem.5735–44. 10.1016/j.jnutbio.2018.03.013
31
KrugerJ.RehmsmeierM. (2006). RNAhybrid: microRNA target prediction easy, fast and flexible.Nucleic Acids Res.34W451–W454. 10.1093/nar/gkl243
32
LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2.Nat. Methods9357–359. 10.1038/nmeth.1923
33
LiB.DeweyC. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome.BMC Bioinformatics12:323. 10.1186/1471-2105-12-323
34
LiM.SunX.CaiH.SunY.PlathM.LiC.et al (2016). Long non-coding RNA ADNCR suppresses adipogenic differentiation by targeting miR-204.Biochim. Biophys. Acta1859871–882. 10.1016/j.bbagrm.2016.05.003
35
LiangM. J.WangZ. P.XuL.LengL.WangS. Z.LuanP.et al (2015). Estimating the genetic parameters for liver fat traits in broiler lines divergently selected for abdominal fat.Genet. Mol. Res.149646–9654. 10.4238/2015.August.14.27
36
LuT. X.RothenbergM. E. (2018). MicroRNA.J. Allergy Clin. Immunol.1411202–1207. 10.1016/j.jaci.2017.08.034
37
MaC.HanM.HeinrichB.FuQ.ZhangQ.SandhuM.et al (2018). Gut microbiome-mediated bile acid metabolism regulates liver cancer via NKT cells.Science360:eaan5931. 10.1126/science.aan5931
38
MaH.ZhangS.ZhangK.ZhanH.PengX.XieS.et al (2019). Identifying selection signatures for backfat thickness in yorkshire pigs highlights new regions affecting fat metabolism.Genes10:254. 10.3390/genes10040254
39
PerteaM.PerteaG. M.AntonescuC. M.ChangT. C.MendellJ. T.SalzbergS. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads.Nat. Biotechnol.33290–295. 10.1038/nbt.3122
40
QuinnJ. J.ChangH. Y. (2016). Unique features of long non-coding RNA biogenesis and function.Nat. Rev. Genet.1747–62. 10.1038/nrg.2015.10
41
RansohoffJ. D.WeiY.KhavariP. A. (2018). The functions and unique features of long intergenic non-coding RNA.Nat. Rev. Mol. Cell Biol.19143–157. 10.1038/nrm.2017.104
42
ResnykC. W.CarreW.WangX.PorterT. E.SimonJ.Le Bihan-DuvalE.et al (2017). Transcriptional analysis of abdominal fat in chickens divergently selected on bodyweight at two ages reveals novel mechanisms controlling adiposity: validating visceral adipose tissue as a dynamic endocrine and metabolic organ.BMC Genomics18:626. 10.1186/s12864-017-4035-5
43
ShenK.RogalaK. B.ChouH. T.HuangR. K.YuZ.SabatiniD. M. (2019). Cryo-EM structure of the human FLCN-FNIP2-Rag-ragulator complex.Cell1791319.e8–1329.e8. 10.1016/j.cell.2019.10.036
44
SunG.LiF.MaX.SunJ.JiangR.TianY.et al (2019). gga-miRNA-18b-3p inhibits intramuscular adipocytes differentiation in chicken by targeting the ACOT13 gene.Cells8:556. 10.3390/cells8060556
45
SunL.LuoH.BuD.ZhaoG.YuK.ZhangC.et al (2013). Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts.Nucleic Acids Res.41:e166. 10.1093/nar/gkt646
46
TayY.RinnJ.PandolfiP. P. (2014). The multilayered complexity of ceRNA crosstalk and competition.Nature505344–352. 10.1038/nature12986
47
ThonV. J.KhalilM.CannonJ. F. (1993). Isolation of human glycogen branching enzyme cDNAs by screening complementation in yeast.J. Biol. Chem.2687509–7513.
48
TsunZ. Y.Bar-PeledL.ChantranupongL.ZoncuR.WangT.KimC.et al (2013). The folliculin tumor suppressor is a GAP for the RagC/D GTPases that signal amino acid levels to mTORC1.Mol. Cell52495–505. 10.1016/j.molcel.2013.09.016
49
WandersR. J. (2014). Metabolic functions of peroxisomes in health and disease.Biochimie9836–44. 10.1016/j.biochi.2013.08.022
50
WangG.KimW. K.ClineM. A.GilbertE. R. (2017). Factors affecting adipose tissue development in chickens: a review.Poult. Sci.963687–3699. 10.3382/ps/pex184
51
WangL.FengZ.WangX.WangX.ZhangX. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data.Bioinformatics26136–138. 10.1093/bioinformatics/btp612
52
WangW.DuZ. Q.ChengB.WangY.YaoJ.LiY.et al (2015). Expression profiling of preadipocyte microRNAs by deep sequencing on chicken lines divergently selected for abdominal fatness.PLoS One10:e0117843. 10.1371/journal.pone.0117843
53
XiaoT.LiuL.LiH.SunY.LuoH.LiT.et al (2015). Long noncoding RNA ADINR regulates adipogenesis by transcriptionally activating C/EBPalpha.Stem Cell Rep.5856–865. 10.1016/j.stemcr.2015.09.007
54
YuX.FangX.GaoM.MiJ.ZhangX.XiaL.et al (2020). Isolation and identification of bovine preadipocytes and screening of MicroRNAs associated with adipogenesis.Animals10:818. 10.3390/ani10050818
55
ZhangG.LiS.LuJ.GeY.WangQ.MaG.et al (2018). LncRNA MT1JP functions as a ceRNA in regulating FBXW7 through competitively binding to miR-92a-3p in gastric cancer.Mol. Cancer17:87. 10.1186/s12943-018-0829-6
56
ZhangM.HanY.ZhaiY.MaX.AnX.ZhangS.et al (2020a). Integrative analysis of circRNAs, miRNAs, and mRNAs profiles to reveal ceRNAs networks in chicken intramuscular and abdominal adipogenesis.BMC Genomics21:594. 10.1186/s12864-020-07000-3
57
ZhangM.MaX.ZhaiY.ZhangD.SuiL.LiW.et al (2020b). Comprehensive transcriptome analysis of lncRNAs reveals the role of lncAD in chicken intramuscular and abdominal adipogenesis.J. Agric. Food Chem.683678–3688. 10.1021/acs.jafc.9b07405
58
ZhangT.ZhangX.HanK.ZhangG.WangJ.XieK.et al (2017a). Analysis of long noncoding RNA and mRNA using RNA sequencing during the differentiation of intramuscular preadipocytes in chicken.PLoS One12:e0172389. 10.1371/journal.pone.0172389
59
ZhangT.ZhangX.HanK.ZhangG.WangJ.XieK.et al (2017b). Genome-wide analysis of lncRNA and mRNA expression during differentiation of abdominal preadipocytes in the chicken.G37953–966. 10.1534/g3.116.037069
60
ZhengQ.LinJ.HuangJ.ZhangH.ZhangR.ZhangX.et al (2017). Reconstitution of UCP1 using CRISPR/Cas9 in the white adipose tissue of pigs decreases fat deposition and improves thermogenic capacity.Proc. Natl. Acad. Sci. U.S.A.114E9474–E9482. 10.1073/pnas.1707853114
Summary
Keywords
whole transcriptome analysis, lncRNA, miRNA, mRNA, lipid synthesis, adipogenesis
Citation
Guo L, Chao X, Huang W, Li Z, Luan K, Ye M, Zhang S, Liu M, Li H, Luo W, Nie Q, Zhang X and Luo Q (2021) Whole Transcriptome Analysis Reveals a Potential Regulatory Mechanism of LncRNA-FNIP2/miR-24-3p/FNIP2 Axis in Chicken Adipogenesis. Front. Cell Dev. Biol. 9:653798. doi: 10.3389/fcell.2021.653798
Received
15 January 2021
Accepted
07 April 2021
Published
24 June 2021
Volume
9 - 2021
Edited by
Xiajun Li, ShanghaiTech University, China
Reviewed by
Guang Hu, National Institute of Environmental Health Sciences, United States; Yanjie Chao, Institut Pasteur of Shanghai (CAS), China
Updates
Copyright
© 2021 Guo, Chao, Huang, Li, Luan, Ye, Zhang, Liu, Li, Luo, Nie, Zhang and Luo.
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: Qingbin Luo, qbluo@scau.edu.cn
†These authors have contributed equally to this work
This article was submitted to Developmental Epigenetics, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
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.