https://orcid.org/0000-0001-8442-8861
Figures
Abstract
Circular RNAs (circRNAs) play important roles in adipogenesis. However, studies on circRNA expression profiles associated with the development of abdominal adipose tissue are lacking in chickens. In this study, 12 cDNA libraries were constructed from the abdominal adipose tissue of Chinese domestic Gushi chickens at 6, 14, 22, and 30 weeks. A total of 1,766 circRNAs were identified by Illumina HiSeq 2500 sequencing. These circRNAs were primarily distributed on chr1 through chr10 and sex chromosomes, and 84.95% of the circRNAs were from gene exons. Bioinformatic analysis showed that each circRNA has 35 miRNA binding sites on average, and 62.71% have internal ribosome entry site (IRES) elements. Meanwhile, these circRNAs were primarily concentrated in TPM < 0.1 and TPM > 60, and their numbers accounted for 18.90% and 80.51%, respectively, exhibiting specific expression patterns in chicken abdominal adipose tissue. In addition, 275 differentially expressed (DE) circRNAs were identified by comparison analysis. Functional enrichment analysis showed that the parental genes of DE circRNAs were primarily involved in biological processes and pathways related to lipid metabolism, such as regulation of fat cell differentiation, fatty acid homeostasis, and triglyceride homeostasis, as well as fatty acid biosynthesis, fatty acid metabolism, and glycerolipid metabolism. Furthermore, ceRNA regulatory networks related to abdominal adipose development were constructed. The results of this study indicated that circRNAs can regulate lipid metabolism, adipocyte proliferation and differentiation, and cell junctions during abdominal adipose tissue development in chickens through complex ceRNA networks between circRNAs, miRNAs, genes, and pathways. The results of this study may help to expand the number of known circRNAs in abdominal adipose tissue and provide a valuable resource for further research on the function of circRNAs in chicken abdominal adipose tissue.
Citation: Jin W, Zhao Y, Zhai B, Li Y, Fan S, Yuan P, et al. (2021) Characteristics and expression profiles of circRNAs during abdominal adipose tissue development in Chinese Gushi chickens. PLoS ONE 16(4): e0249288. https://doi.org/10.1371/journal.pone.0249288
Editor: Michael H. Kogut, USDA-Agricultural Research Service, UNITED STATES
Received: October 24, 2020; Accepted: March 15, 2021; Published: April 15, 2021
Copyright: © 2021 Jin et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The raw data from this study have been stored in NCBI Sequence Read Archive (accession numbers PRJNA528858 and PRJNA551368).
Funding: This study was funded by a grant from the National Natural Science Foundation of China (32072692 and 31572356), the Program for Innovation Research Team of Ministry of Education (IRT16R23) and the Scientific Studio of Zhongyuan Scholars (No. 30601985).
Competing interests: The authors declare that they have no competing interests.
Introduction
Chickens are an economically important agricultural animal [1]. Excessive fat accumulation in the abdominal cavity not only reduces feed utilization efficiency but also affects meat quality in poultry production [2]. Therefore, demonstrating the molecular mechanism governing chicken abdominal adipose tissue development is of great significance to chicken genetic breeding and production. CircRNAs are a type of nonphosphorylated, single-stranded, covalently closed noncoding RNA that is widely present in various tissues and cell types of different species [3]. Many studies have shown that circRNAs may be involved in the regulation of various biological processes and affect cell physiology through various molecular mechanisms [4]. In particular, some studies have shown that circRNAs are also involved in the regulation of adipogenesis. For example, circFUT10 has been reported to promote adipocyte proliferation and inhibit adipocyte differentiation through sponge adsorption of let-7 [5]. Therefore, circRNAs can be employed as new targets to elucidate the molecular mechanism underlying chicken abdominal fat development.
CircRNAs are produced by a special reverse splicing process that covalently closes the 5’ end of the upstream exon and the 3’ end of the downstream exon of pre-mRNA, which primarily exist in the cytoplasm [6]. It is known that more than 20% of the expressed genes in the genome can produce circRNAs, which present notable abundance, stability, and diversity of expression profiles in various types of tissues and cells [3]. CircRNAs can be used as miRNA- or RNA-binding protein decoys to regulate gene expression or protein translation [4] and thus play pivotal roles in various physiological activities and diseases. Therefore, circRNAs have recently become a heavily researched topic in the field of RNA. At present, research on circRNAs related to fat deposition is primarily focused on identifying novel circRNAs from adipose tissue or adipocytes in different anatomic sites. For example, Sun (2020) detected the circRNA expression profiles of preadipocytes and adipocytes in human visceral adipose tissue by microarray technology. Compared with preadipocytes, 2,215 circRNAs were significantly upregulated, and 1,865 were downregulated, in adipocytes [7]. In the Chinese Erhualian pig, Liu (2019) identified many DE circRNAs in adipose tissue of intramuscular, subcutaneous, peritoneal and mesenteric adipose tissues [8] and between subcutaneous preadipocytes and adipocytes [9]. In yaks (Bos grunniens), a previous study found 211 DE circRNAs from adipose tissue in yaks of two age groups [10], and another study identified 7,203 circRNAs and 136 DE circRNAs during yak adipocyte differentiation [11]. In addition, Huang (2019) identified 5,141 circRNAs and 252 DE circRNAs from the adipose tissue of adult and young Chinese buffaloes [12]. These studies indicated that many circRNAs were involved in regulating adipose tissue development or adipose deposition in animals. Although these results helped to elucidate the molecular mechanisms underlying adipose tissue development or fat deposition, current research has primarily concentrated on humans and a few other species. In addition, only a few studies have described the mechanism of specific circRNAs in adipogenesis [5, 13, 14]. At present, the function of circRNAs in the development or deposition of animal adipose tissue has not been fully elucidated. The expression of circRNAs is known to have obvious tissue or cell specificity [15]. Therefore, the identification of circRNAs related to adipose tissue development in economically important animals is an important step to better study circRNA function in the formation of adipose tissue-related traits.
In view of the importance of chickens in poultry production, the roles played circRNAs in the formation of economic traits have been studied increasingly frequently by researchers. Previous studies have identified many novel circRNAs. For instance, Shen (2019, 2020) quantified circRNAs from theca cells [16] and granulosa cells [17] during follicular development in chickens, counting 14,502 and 1,1642, respectively. Zhang (2019) studied the differential expression profiles of circRNAs in chicken hepatocytes in GHR antisense transcript overexpression groups and control groups and identified 4,772 circRNAs and 92 DE circRNAs [18]. Ouyang (2018) identified 13,377 circRNAs and 462 DE circRNAs from the leg muscles of female Xinghua chickens of different embryo ages and at one day post-hatching [19]. In addition, several studies have demonstrated the mechanism of action of specific circRNAs in chicken skeletal muscle satellite cells or myoblasts [20–23]. These studies have proven that many circRNAs participate in the formation and regulation of economically important traits in chickens, which helps to elucidate the molecular mechanism governing the formation and regulation of economically important traits in chicken. To the best of our knowledge, however, there have been no reports characterizing circRNAs in relation to adipose tissue development in chickens.
Gushi chickens are native Chinese chickens that originate from Gushi County, Henan Province, China. These chickens have been bred on a large scale due to their excellent meat quality. To elucidate the role played by circRNAs in the development of abdominal fat, we constructed the circRNA transcriptome profile of 6-, 14-, 22- and 30-week-old Gushi chickens and the regulatory network of circRNAs related to abdominal fat development in Gushi chickens. The results of this study may not only help to elucidate the posttranscriptional regulation of abdominal fat deposition in chickens but also provide basic data for resource protection and the mining of useful characteristics in Gushi chickens.
Materials and methods
Ethics statement
All animal experiments were performed according to the program approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University (Permit Number: 17–0118).
Experimental animals and total RNA extraction
In this study, domestic Gushi chickens were used as experimental animals. The feeding method was as described by Fu (2018) [24]. Before reaching 14 weeks of age, chickens were fed a diet containing 18.5% crude protein and 12.35 MJ/kg energy, and after reaching the age of 14 weeks, they were fed a diet containing 15.6% crude protein and 12.75 MJ/kg energy. At 6, 14, 22, and 30 weeks of age, three individuals were randomly selected at each developmental stage. After dissection, abdominal adipose tissues were collected and stored at -80°C. Total RNA was extracted by using TRIzol reagent (Takara, Dalian, China) according to the instructions of the manufacturer. A 1% agarose gel electrophoresis and Agilent Bioanalyzer 2100 system (Agilent, Santa Clara, CA, USA) were used to evaluate the quality of the total RNA.
Library construction and sequence analysis
The initial total RNA content of each sample used for library construction was 5 μg. rRNA was removed from total RNA by an Epicentre Ribo-ZeroTM kit (Illumina, San Diego, USA). A short fragment of 250–300 bp was prepared as a template for the synthesis of double-stranded cDNA. Subsequently, the A tail was added, the sequencing adapter was connected, and a chain-specific cDNA library was obtained by PCR amplification and enrichment. The library was diluted to 1 ng/μl, and an Agilent 2100 (Agilent, Santa Clara, CA, USA) was used to detect the insert size of the library to meet the expected size. Quantitative PCR (qPCR) was used to quantify the library to ensure that the effective concentration of the library was > 2 nm. Finally, the qualified library was sequenced on the Illumina HiSeq 2500 platform (Novogene, Beijing, China). The quality of sequencing data was evaluated rigorously. The adaptor sequences and low-quality sequences in raw data were eliminated. The Q20, Q30 and GC contents of clean data were calculated. The obtained high-quality reads were used for subsequent analysis. The paired-end clean reads were compared with the reference genome by Bowtie software (Releases 1.2.3) [25].
Identification of circRNAs
CircRNAs were identified using find_circ [26] and CIRI2 [27] software. Circos software (v0.69–9) [28] was used to construct circos diagrams. The reverse shear method was used to extract the connection points that were not aligned with reads. The prediction analysis of circRNAs for each software was performed according to the system default parameters. Only those sequences containing more than 2 independent junction-spanning reads that met the Breathnach-Chambon rule (GU/AG rule) were identified as potential circRNAs. To reduce false positives, the shared circRNAs in the results predicted by the two software programs were identified as circRNAs based on their position in the chromosomes.
Identification of DE circRNAs
The read count of each circRNA in different libraries was calculated according to back-spliced reads per million mapped reads, and the transcripts per million clean tags (TPM) were used to correct it. The normalized expression value was calculated by the following formula: normalized expression level = (read count × 1,000,000)/libsize (libsize is the sum of circRNA read count). Based on the read count, DE circRNAs were analyzed by DESeq2 software (v3.12) [29] in six comparison groups, including W14 vs. W6, W22 vs. W6, W30 vs. W6, W14 vs. W22, W14 vs. W30, and W22 vs. W30. p < 0.05 was considered to be the threshold criterion for determining the DE circRNAs.
Potential function analysis of circRNAs
First, a functional enrichment analysis was performed on the parental genes of DE circRNAs. The GoSeq R package (v3.12) [30] was used for enrichment analysis of Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was conducted using KOBAS software (v3.0) [31]. p < 0.05 was used as the threshold for evaluating GO terms or KEGG pathways for significant enrichment. Second, miRNA target sites were predicted by miRanda software (Releases 3.3a), and the analysis of the interaction between circRNAs and miRNAs was completed. In addition, IRESfinder software (v1.1.0) [32] was utilized to predict and identify the internal ribosome entry site (IRES) on the circRNA sequences to determine the circRNAs with the potential to translate polypeptides or proteins.
Interactive network analysis of circRNA-miRNA-gene pairs
The mRNA and miRNA transcriptome profiles of Gushi chicken abdominal adipose tissue at 6, 14, 22, and 30 weeks were also constructed using the same tissue samples. Based on these sequencing data, miRanda software (Releases 3.3a) was used to predict the binding sites of DE miRNAs on circRNA and miRNA target genes. By Pearson’s correlation analysis, 5,455 DE circRNA-miRNA-gene pairs were identified according to the expression levels of circRNAs, miRNAs and genes in abdominal adipose tissue of Gushi chickens. The 431 target genes in these circRNA-miRNA-gene pairs were annotated with GO and KEGG using DAVID online tools (https://david.ncifcrf.gov/). In addition, the circRNA-miRNA-gene pairs were selected from the pathways associated with lipid metabolism, cell proliferation and differentiation, and cell junctions, and the interaction networks were constructed by Cytoscape software (v3.2).
Reliability verification of circRNAs
First, reverse transcription PCR analysis and Sanger sequencing were used to verify the circRNA sequences from RNA-seq analysis. Total RNA was extracted from abdominal adipose tissue samples of Gushi chickens at different weeks. After the RNA quality was qualified, it was mixed in equal amounts for cDNA synthesis. Ten circRNAs, namely, gga_circ_0000348, gga_circ_0003969, gga_circ_0005065, gga_circ_0001623, gga_circ_0000833, gga_circ_0004577, gga_circ_0003686, gga_circ_0003828, gga_circ_0003244 and gga_circ_0002520, were randomly selected and verified by PCR amplification using divergent primers (S1 Table). Subsequently, the PCR product fragment size was detected by gel electrophoresis; the sequence was determined by Sanger sequencing; finally, DNAMAN software (v6.25) was used to compare the sequence obtained by Sanger sequencing with RNA-seq data and the chicken reference genome.
Second, the resistance of the above nine circRNAs to RNase R (Geneseed, Guangzhou, China) digestion was detected by quantitative real-time PCR (qRT-PCR). Before cDNA synthesis, total RNA was treated with RNase R (RNR-07250, Epicentre). Subsequently, using an RT-PCR kit (Takara, Dalian, China), RNA treated with RNase R and RNA without RNase R treatment were separately synthesized into cDNA. All qRT-PCR analyses were completed in accordance with the SYBR Green kit (Takara, Dalian, China) product instructions. All reactions were repeated three times. The 2−ΔΔct method [33] was used to calculate the relative expression of circRNAs, and the linear gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a control.
Plasmid construction and dual luciferase reporter assay
The miR-215-5p mimics and mimics NC were synthesized by GenePharma (Shanghai, China). For the construction of the dual-luciferase reporter vector, wild-type and mutated sequences in the 3′UTR of NCOA3 and the partial region of gga_circ_0002520 (named WT-NCOA3, Mut-NCOA3, gga_circ_2520-WT and gga_circ_2520-Mut, respectively), which contain the binding sites of miR-215-5p, were synthesized and inserted into psiCHECK2 vectors (Promega, Madison, WI, USA) according to instructions using XhoI and NotI restriction sites.
After planting DF-1 cells evenly in 96-well plates and growing to 70%, four groups (wild-type psiCHECK2 plasmids and miR-215-5p mimics as the treatment, mutated psiCHECK2 plasmids and miR-215-5p mimics, wild-type psiCHECK2 plasmids and miR-215-5p mimics NC, mutated type psiCHECK2 plasmids and miR-215-5p mimics NC) were set and cotransfected. To confirm the target relationship between gga_circ_2520 and miR-215-5p, another dual-luciferase reporter assay method was employed. Three comparison groups (WT-NCOA3; WT-NCOA3 and miR-215-5p mimics; WT-NCOA3, miR-215-5p mimics and gga_circ_2520-WT) were set and cotransfected. After 48 h, a dual luciferase assay system kit (Promega, Madison, WI, USA) and fluorescence/multidetection microplate reader (BioTek, Winooski, VT, USA) were employed to detect luminescent signals. Firefly luciferase activities were normalized to Renilla luminescence in each well. Detailed protocols were described in the manufacturer’s instructions.
Results and discussion
Overview of library sequencing
A total of 12 cDNA libraries were constructed using abdominal fat tissue obtained from Gushi chickens aged 6, 14, 22, and 30 weeks. Raw reads obtained from each library ranged from 93,729,158 to 116,176,278. After removing adapter reads and low-quality reads, the amount of clean base data obtained by each library ranged from 13.49 G to 16.73 G, accounting for more than 95%. In addition, reference genome mapping of the sequence was also completed. The results showed that the percent range of total mapped reads/clean reads in different libraries was 78.27–82.8%, of which the percent range of multiple mapped reads/clean reads was 2.21–8.29%, and the percent range of uniquely mapped reads/clean reads was 74.02–80.08%. Moreover, 64.46–70.7% reads were aligned to the protein-coding region, and 19.01–26.09% of the reads were compared to the intron region. The sequence output and quality evaluation of each library are shown in Table 1 and S1 Fig.
CircRNAs expressed in chicken abdominal adipose tissues
A total of 1,766 novel circRNAs were identified from the 12 cDNA libraries in Gushi chicken abdominal adipose tissue (S1 File). In total, 1,637, 1,711, 1,735 and 1,663 circRNAs were expressed in the abdominal adipose tissue of Gushi chicken at 6, 14, 22, and 30 weeks of age, respectively. These circRNAs were primarily located in exons of genes with an average ratio of 84.95% followed by intergenic regions (9.62%) and introns (5.43%). As shown in S2 Fig, the full length of these circRNAs ranged from 185 nt to 99,591 nt, and their splicing length was 43–1842 nt. The full length of 543 circRNAs was more than 10,000 nt, and the length range of 1,071 circRNAs was 1,000–9,999 nt. Moreover, these circRNAs were distributed on chr1 through chr28 and sex chromosomes, especially in chr1 (19.03%), chr2 (14.04%), chr3 (13.93%), chr4 (7.19%) chr5 (5.83%), chrZ (4.70%), chr6 (4.13%), chr7 (3.96%), chr8 (3.40%), chr10 (2.66%), and chr9 (1.98%) (S2 Fig).
Previously, several studies have reported circRNAs expressed in adipose tissue, preadipocytes and adipocytes in different anatomical parts of humans and several other species, such as pigs, yaks, and buffalo [7–12]. These studies indicate that a large quantity of circRNAs participate in the regulation of fat tissue development or adipogenesis. The results of this study are consistent with the findings of previous studies, suggesting that the development or deposition of chicken abdominal fat is also regulated by many circRNAs. Ten circRNAs were randomly selected and verified by reverse transcription PCR analysis, Sanger sequencing and RNase R digestion tests. These results proved that the identification of circRNAs in this study was accurate and reliable (Figs 1A, 1B, 9A and 9B). Therefore, the results of this study provide a basis for further research on the molecular mechanism for the regulation of chicken abdominal adipogenesis by circRNAs.
(A) Electrophoretogram of the PCR products. M: marker 500. (B) The head-to-tail junctions of 9 circRNAs were confirmed by Sanger sequencing. The black line indicates the position of reverse splicing. The area marked with red (or black) is the end (or start) sequence of the circRNAs (junction). (C) qRT-PCR was used to detect the resistance of circRNA to RNase R digestion. GAPDH was used as a linearity control. Data are expressed as means ± SEM (n = 3).
Expression profiles of circRNAs during abdominal adipose development
The expression of circRNAs is known to be tissue- and developmental stage-specific [4, 34]. The TPM value distribution of circRNA in this study showed that the circRNA expression levels of all samples had two peaks, indicating that the expression patterns of circRNAs in the four stages of Gushi chicken abdominal fat development were highly consistent (Fig 2A). These circRNAs were mainly concentrated in TPM < 0.1 and TPM > 60, and their numbers accounted for 18.90% and 80.51%, respectively. This finding indicated that circRNAs had specific expression patterns in chicken abdominal adipose tissue. This specific expression pattern is similar to that of circRNAs in porcine intramuscular, subcutaneous, retroperitoneal, and mesenteric adipose tissues [8].
(A) Density distribution of TPM values of circRNAs in each library. (B) Box plot of TPM values of circRNAs in 12 libraries. (C) Two subclusters of the circRNA expression profile in Gushi chicken abdominal adipose tissue. K-means clustering was used to represent the circRNA dynamic profile. The gray line is a line chart representing the relative expression level of a given circRNA for the four developmental stages. The blue line is a line chart representing the relative mean expression level of all circRNAs for the four developmental stages. The red line is for reference; circRNAs above the red line were upregulated, and those below the red line were downregulated. (D) Hierarchical clustering of the DE circRNAs in 12 libraries. (E) K-means cluster map of the DE circRNAs.
However, the expression profiles of circRNAs were different in different developmental stages of abdominal adipose tissue in Gushi chickens. In particular, there were many circRNAs with notable changes in expression in each developmental stage (Fig 2B). Meanwhile, we found that 1,518 of the identified circRNAs were expressed in all four developmental stages. The expression profiles of these circRNAs exhibited clear temporal characteristics and were divided into eight different subclusters (Fig 2C and S3 Fig). We also found that one, three, six, and two circRNAs were expressed specifically in abdominal adipose tissue at 6, 14, 22, and 30 weeks of age, respectively (S1 File). In addition, 275 DE circRNAs were identified from six combinations (S2 File). Accordingly, the DE circRNAs in W14 vs. W6, W22 vs. W6, W30 vs. W6, W22 vs. W14, W30 vs. W14, W30 vs. W22 combinations were 50, 80, 87, 61, 67, 57, respectively. These DE circRNAs were clustered together repeatedly within the group, indicating that the differences within the group were less than the differences between the groups and proving that the data were reliable (Fig 2D). The expression patterns of these DE circRNAs were divided into four different modes and had clear temporal characteristics at different developmental stages of abdominal fat tissues (Fig 2E). The distribution of DE circRNAs showed that there was no common DE circRNA between all six comparison combinations or any five comparison combinations (Fig 3). Only gga_circ_0001526 and gga_circ_0002515 were DE in four comparison combinations, including W22 vs. W6, W30 vs. W6, W22 vs. W14, and W30 vs. W14. These findings indicated that the expression profile of circRNAs in chicken abdominal adipose tissue had clear developmental stage specificity.
Potential function of circRNAs in chicken abdominal adipose tissues
A large number of circRNAs are derived from protein-coding genes [35–37]. In this study, we observed that 1,579 circRNAs were derived from 1,076 known genes, and the parental genes of some circRNAs were closely related to lipid metabolism or adipogenesis. For example, the host gene ASXL2 of gga_circ_0005102 can increase the activity of liver X receptor alpha and play a key role in lipid metabolism [38, 39], and the host gene cytochrome P450 family 39 subfamily A member 1 (CYP39A1) for gga_circ_0005166 can participate in lipoprotein transport and cholesterol metabolism [40]. Moreover, the diacylglycerol O-acyltransferase 2 (DGAT2) gene, which can encode three circRNAs, gga_circ_0002513, gga_circ_0002515 and gga_circ_0002517, plays a vital role in the synthesis of triacylglycerol [41]. To further elucidate the potential roles played by circRNAs in chicken abdominal fat deposition, functional enrichment analysis of the parental genes of DE circRNAs was performed. The results showed that the functions of DE circRNAs were significantly different in different developmental stages of abdominal adipose tissues. The biological functions of DE circRNAs in abdominal adipose tissue were mainly related to the regulation of fat cell differentiation between 6 and 14 weeks of age, and they functioned primarily in biological processes, such as fatty acid homeostasis, long-chain fatty-acyl-CoA metabolic process, fat pad development, and triglyceride homeostasis, after 14 weeks of age (Fig 4A). The regulatory effects of DE circRNAs were determined to be consistent with the physiological characteristics of abdominal adipose tissue development in Gushi chickens.
(A) The objective was to compare the different biological processes related to lipid metabolism in W14 and W06, W22 and W14, and W30 and W22 combinations. The y-axis represents the GO terms that were significantly enriched in different comparisons, and the x-axis represents the enrichment factors in the corresponding GO terms. (B) The first 20 KEGG pathways were enriched in the W22 vs. W6 combination. (C) The first 20 KEGG pathways were enriched in the W30 vs. W22 combination.
Moreover, the results of KEGG pathway analysis showed that some pathways related to lipid metabolism, such as pyruvate metabolism, fatty acid metabolism, fatty acid biosynthesis, the insulin signaling pathway and the Wnt signaling pathway, were enriched in the W22 vs. W6 combination (Fig 4B). The phosphatidylinositol signaling system and glycerolipid metabolism were significantly enriched in the W30 vs. W22 combination (Fig 4C). In particular, we found that the enrichment pathways propionic acid metabolism, fatty acid biosynthesis and pyruvate metabolism formed an interactive regulatory network through the DE circRNAs and their parental genes, as well as miRNAs and their targets (Fig 5). Acetyl-CoA carboxylase alpha (ACACA) and acetyl-CoA carboxylase beta (ACACB) were significantly enriched in all three pathways, and they were the parental genes for gga_circ_0001526 and gga_circ_0001148, respectively. MiR-125b-5p, miR-1559-5p, miR-1736-5p and miR-187-3p had target sites with ACACA. Acyl-CoA synthetase short-chain family member 3 (ACSS3) was the parental gene of gga_circ_0002744 and gga_circ_0002746, and miR-193a-3p, miR-193b-3p, miR-194 and miR-30a-3p had target sites. These results indicate that circRNAs may participate in the regulation of abdominal fat development or deposition through their parental genes and corresponding pathways in chickens.
The genes within the propanoate metabolism, fatty acid biosynthesis, and pyruvate metabolism pathways are shown in green, yellow, and blue, respectively. The DE mRNAs are the parental gene of DE circRNAs. There are target interactions between DE mRNAs and miRNAs.
MiRNA target sites on circRNAs
It is well-known that the most significant function of circRNAs is to serve as miRNA sponges and regulate the expression of target genes by inhibiting miRNA activity [42]. Therefore, we predicted the binding sites of known chicken miRNAs on the identified circRNA sequences (S3 File). The results showed that there were 61,988 miRNA binding sites in 1,766 circRNAs with an average of 35 miRNA binding sites on each circRNA and a maximum of 185 miRNA binding sites. These sites could be bound by 1,235 known miRNAs. Thus, a single circRNA can bind to a considerable number of miRNAs, and a miRNA can also target many circRNAs simultaneously. These results are consistent with previous studies in human visceral preadipocytes and adipocytes [7] and yak adipocytes [11]. This finding suggests that circRNAs, as miRNA sponges, have a complex regulatory relationship in abdominal fat tissue development or adipogenesis in chickens.
Literature mining determined that 35 miRNAs of the 1,235 miRNAs with binding sites on circRNAs, such as let-7b, miR-101-3p, miR-103-3p, and miR-10b-5p, were related to lipid metabolism or adipogenesis [43–68]. These miRNAs were also expressed during the development of abdominal fat in Gushi chickens [69]. To characterize the complex interaction between circRNAs and miRNAs, we constructed a circRNA-miRNA interaction network using these 35 miRNAs (Fig 6). In the network, some circRNAs contained multiple target sites for more than four different miRNAs, such as gga_circ_0000582, gga_circ_0001623, gga_circ_0008796, and gga_circ_0008809. Among these miRNAs decoyed by circRNAs, miR-20a can regulate adipocyte differentiation by targeting lysine-specific demethylase 6b [43], miR-144-3p promotes fat formation by releasing C/EBPα from Klf3 and CtBP2 [44], and miR-148a can modulate adipocyte differentiation of mesenchymal stem cells [45]. In addition, miR-204 can target 11 different circRNAs, such as gga_circ_0007413, gga_circ_0007670, gga_circ_0007776, and gga_circ_0008100, in the network, and it has been shown to be involved in the regulation of bovine adipocyte proliferation, differentiation, and apoptosis [46]. These studies indicate that circRNAs may play endogenous competitive regulatory roles through interaction with multiple miRNAs during the development of chicken abdominal fat tissue.
circRNAs and miRNAs are represented by ellipses and triangles, respectively.
IRES elements on circRNAs
The internal ribosome entry site (IRES) is an important regulatory element for RNA translation without relying on the 5′ cap structure, and it can mediate the assembly of ribosomal larger and smaller subunits. Studies have shown that circRNAs contain IRES, and some circRNAs can also encode proteins or functional peptides in vivo and in vitro [70–72]. Therefore, we predicted the potential IRES elements in the identified circRNA sequences. There were IRES elements on 1,108 circRNAs, accounting for 62.71% of the predicted circRNA number (Fig 7). The top 10 circRNAs in the predicted scores were gga_circ_0007953, gga_circ_0005736, gga_circ_0005910, gga_circ_0001863, gga_circ_0001990, gga_circ_0005598, gga_circ_0006417, gga_circ_0007094, gga_circ_0004846, and gga_circ_0002156. These results indicated that most circRNAs in chicken abdominal adipose tissue may have the potential to encode peptides or proteins.
The x-axis indicates the interval of the predicted score, and the y-axis represents the percentage of circRNAs within the corresponding predicted score interval for all circRNAs with IRES. CircRNAs with predicted scores of > 0.5 were considered to have IRES elements. A total of 1,108 circRNAs with IRES were identified.
CircRNA-miRNA-gene networks during abdominal adipose development
Abdominal fat tissue development is a complex process participating in many biological processes and cellular events, which are regulated by multiple signaling pathways. It is known that circRNAs regulate the function of miRNAs by a ceRNA mechanism [73], thereby impacting the posttranscriptional regulation of miRNAs in downstream genes [74] and regulating biological processes by affecting signaling pathways [4]. For this reason, we analyzed the interactions among genes, miRNAs, circRNAs, and pathways during the development of Gushi chicken abdominal fat. Based on the results of KEGG enrichment analyses using the whole transcriptome sequencing data of Gushi chicken abdominal fat at 6, 14, 22, and 30 weeks of age, we identified 53, 83, and 56 circRNA-miRNA-gene pairs associated with lipid metabolism, cell proliferation and differentiation, and cell junction pathways, respectively, and constructed three ceRNA networks (Fig 8, S4 and S5 Figs). The lipid metabolism-related ceRNA network contains 32 circRNAs, 19 miRNAs, and 13 genes and involves 13 pathways, including glycerolipid metabolism, glycerophospholipid metabolism, fatty acid elongation, sphingolipid metabolism, fatty acid degradation, alpha-linolenic acid metabolism, adipocytokine signaling pathway, fatty acid metabolism, the PPAR signaling pathway, fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, ether lipid metabolism and insulin signaling pathway (Fig 8).
The ellipse, triangle and box nodes represent miRNAs, circRNAs, mRNAs, respectively. Pink indicates upregulation, and blue indicates downregulation.
In these networks, some miRNAs, such as gga-miR-130b-3p, gga-miR-1744-3p, gga-miR-200a-3p, gga-miR-1a-3p, gga-miR-34c-5p, gga-miR-449a, gga-miR-454-3p, and gga-miR-9-5p, targeted multiple circRNAs, and some genes acted as components of multiple pathways and were targeted by multiple miRNAs. Meanwhile, some circRNAs, such as gga_circ_0000348, gga_circ_0000525, gga_circ_0007512, and gga_circ_0008796, interacted with more than two miRNAs. As key nodes, these factors constituted a complex ceRNA network and synergistically regulated lipid metabolism and multiple cellular events during the development of abdominal fat tissue in Gushi chickens. In addition, the ceRNA gga_circ_0002520-miR-215-5p-NCOA3 in the lipid metabolism-related circRNA-miRNA-gene interaction network was selected and verified by luciferase reporter assay. The results showed that there were target sites of miR-215-5p in the gga_circ_0002520 sequence and NCOA3 mRNA 3´-UTR (Fig 9), which demonstrated that the results of our circRNA-miRNA-gene interaction analysis were reliable. These results indicated that circRNAs play critical roles in the regulation of chicken abdominal fat tissue development or fat deposition through the complex networks formed by circRNAs, miRNAs and genes and corresponding pathways.
(A) Convergent and divergent primers amplify circRNAs in cDNA. The head-to-tail junctions of gga_circ_0002520 were confirmed by Sanger sequencing. The black line indicates the position of reverse splicing. (B) qRT-PCR was used to detect the resistance of gga_circ_0002520 to RNase R digestion. β-actin was used as a linearity control. Data are expressed as means ± SEM (n = 3). (C) Luminescence was measured after cotransfecting the wild-type or mutant sequence of NCOA3-3′UTR with miR-215-5p mimics (or mimics NC) in DF-1 cells. (D) Luminescence was measured after cotransfecting the wild-type or mutant linear sequence of gga_circ_0002520 with miR-215-5p mimics (or mimics NC) in DF-1 cells. (E) A double luciferase assay was performed in DF-1 cells 48 h after miR-215-5p mimic transfection with NCOA3-3′UTR-WT and gga_circ_0002520-WT.
Conclusions
The present study identified many circRNAs and DE circRNAs related to abdominal adipose tissue development and assessed the characteristics of circRNAs and their expression profiles in abdominal adipose tissue during late development in Gushi chickens. We found that circRNAs had specific expression patterns and obvious temporal characteristics in chicken abdominal adipose tissue, and most circRNAs had the potential to encode proteins or functional peptides. In addition, we also found that the parental genes of DE circRNAs were primarily enriched in some pathways related to lipid metabolism, such as pyruvate metabolism, fatty acid metabolism, and fatty acid biosynthesis, and these circRNAs play critical roles in the regulation of chicken abdominal adipose tissue development by complex ceRNA networks. Our results provide novel insight and a valuable resource helping to elucidate the regulatory roles played by circRNAs in chicken abdominal adipose tissue. Further research is required to clarify the function of the identified circRNAs in chicken abdominal adipose tissue development or deposition.
Supporting information
S1 Fig. Quality evaluation and processing of sequencing data.
(A) The sequencing error rate of each base position. (B) The content of each base during the sequencing process. (C) The proportion of different types of raw data. (D) Classification of mapped reads in 12 libraries.
https://doi.org/10.1371/journal.pone.0249288.s001
(TIF)
S2 Fig. Characteristics of the identified circRNAs.
(A) Length distribution. (B) The distribution of circRNAs in different chromosomes.
https://doi.org/10.1371/journal.pone.0249288.s002
(TIF)
S3 Fig. Dynamic expression profile of the identified circRNAs.
https://doi.org/10.1371/journal.pone.0249288.s003
(TIF)
S4 Fig. CircRNA-miRNA-gene interaction network contains five pathways related to cell proliferation and differentiation.
The ellipse, triangle and box nodes represent DE miRNAs, DE circRNAs, and DE mRNAs, respectively. Pink indicates upregulation, and blue indicates downregulation. The five pathways are the Wnt signaling pathway, FoxO signaling pathway, p53 signaling pathway, TGF-beta signaling pathway and MAPK signaling pathway.
https://doi.org/10.1371/journal.pone.0249288.s004
(TIF)
S5 Fig. CircRNA-miRNA-gene interaction network containing five cell junction-related pathways.
The ellipse, triangle and box nodes represent miRNAs, circRNAs, mRNAs, respectively. Pink indicates upregulation, and blue indicates downregulation. The five pathways are cell adhesion molecules, focal adhesion, adherens junctions, tight junctions and gap junctions.
https://doi.org/10.1371/journal.pone.0249288.s005
(TIF)
S1 Table. Primers used for the validation of circRNAs.
https://doi.org/10.1371/journal.pone.0249288.s006
(DOCX)
S1 File. Details of the circRNAs identified in this study.
https://doi.org/10.1371/journal.pone.0249288.s007
(XLSX)
S2 File. Details of the significant differentially expressed circRNAs identified in this study.
https://doi.org/10.1371/journal.pone.0249288.s008
(XLSX)
S3 File. Binding site analysis of miRNAs on the circRNAs.
https://doi.org/10.1371/journal.pone.0249288.s009
(XLSX)
References
- 1. Burt DW. Emergence of the chicken as a model organism: implications for agriculture and biology. Poult Sci. 2007; 86(7):1460–71. pmid:17575197
- 2. Hood RL. The cellular basis for growth of the abdominal fat pad in broiler-type chickens. Poult Sci. 1982; 61(1):117–21. pmid:7088775
- 3. Bose R, Ain R. Regulation of Transcription by Circular RNAs. Adv Exp Med Biol. 2018; 1087:81–94. pmid:30259359
- 4. Yu CY, Kuo HC. The emerging roles and functions of circular RNAs and their generation. J Biomed Sci. 2019; 26(1):29. pmid:31027496
- 5. Jiang R, Li H, Yang J, Shen X, Song C, Yang Z, et al. circRNA Profiling Reveals an Abundant circFUT10 that Promotes Adipocyte Proliferation and Inhibits Adipocyte Differentiation via Sponging let-7. Mol Ther Nucleic Acids. 2020; 20:491–501. pmid:32305019
- 6. Enuka Y, Lauriola M, Feldman ME, Sas-Chen A, Ulitsky I, Yarden Y. Circular RNAs are long-lived and display only minimal early alterations in response to a growth factor. Nucleic Acids Res. 2016; 44(3):1370–83. pmid:26657629
- 7. Sun W, Sun X, Chu W, Yu S, Dong F, Xu G. CircRNA expression profiles in human visceral preadipocytes and adipocytes. Mol Med Rep. 2020; 21(2):815–21. pmid:31974620
- 8. Liu X, Wei S, Deng S, Li D, Liu K, Shan B, et al. Genome-wide identification and comparison of mRNAs, lncRNAs and circRNAs in porcine intramuscular, subcutaneous, retroperitoneal and mesenteric adipose tissues. Anim Genet. 2019; 50(3):228–41. pmid:30982992
- 9. Liu X, Liu K, Shan B, Wei S, Li D, Han H, et al. A genome-wide landscape of mRNAs, lncRNAs, and circRNAs during subcutaneous adipogenesis in pigs. J Anim Sci Biotechnol. 2018; 9:76. pmid:30410752
- 10. Wang H, Zhong J, Zhang C, Chai Z, Cao H, Wang J, et al. The whole-transcriptome landscape of muscle and adipose tissues reveals the ceRNA regulation network related to intramuscular fat deposition in yak. BMC Genomics. 2020; 21(1):347. pmid:32381004
- 11. Zhang Y, Guo X, Pei J, Chu M, Ding X, Wu X, et al. CircRNA Expression Profile during Yak Adipocyte Differentiation and Screen Potential circRNAs for Adipocyte Differentiation. Genes (Basel). 2020;11(4): 414. pmid:32290214
- 12. Huang J, Zhao J, Zheng Q, Wang S, Wei X, Li F, et al. Characterization of Circular RNAs in Chinese Buffalo (Bubalus bubalis) Adipose Tissue: A Focus on Circular RNAs Involved in Fat Deposition. Animals (Basel). 2019; 9(7):403. pmid:31266200
- 13. Zhu Y, Gui W, Lin X, Li H. Knock-down of circular RNA H19 induces human adipose-derived stem cells adipogenic differentiation via a mechanism involving the polypyrimidine tract-binding protein 1. Exp Cell Res. 2020;387(2):111753. pmid:31837293
- 14. Liu Y, Liu H, Li Y, Mao R, Yang H, Zhang Y, et al. Circular RNA SAMD4A controls adipogenesis in obesity through the miR-138-5p/EZH2 axis. Theranostics. 2020; 10(10):4705–19. pmid:32292524
- 15. Qin T, Li J, Zhang KQ. Structure, Regulation, and Function of Linear and Circular Long Non-Coding RNAs. Front Genet. 2020;11:150. pmid:32194627
- 16. Shen M, Wu P, Li T, Wu P, Chen F, Chen L, et al. Transcriptome Analysis of circRNA and mRNA in Theca Cells during Follicular Development in Chickens. Genes (Basel). 2020;11(5): 489. pmid:32365656
- 17. Shen M, Li T, Zhang G, Wu P, Chen F, Lou Q, et al. Dynamic expression and functional analysis of circRNA in granulosa cells during follicular development in chicken. BMC Genomics. 2019; 20(1):96. pmid:30700247
- 18. Zhang L, Xu H, Wang Z, Li T, Guo J, Adu-Asiamah P, et al. Identification and characterization of circular RNAs in chicken hepatocytes. Growth Horm IGF Res. 2019; 46–47:16–23. pmid:31125863
- 19. Ouyang H, Chen X, Wang Z, Yu J, Jia X, Li Z, et al. Circular RNAs are abundant and dynamically expressed during embryonic muscle development in chickens. DNA Res. 2018;25(1):71–86. pmid:29036326
- 20. Shen X, Liu Z, Cao X, He H, Han S, Chen Y, et al. Circular RNA profiling identified an abundant circular RNA circTMTC1 that inhibits chicken skeletal muscle satellite cell differentiation by sponging miR-128-3p. Int J Biol Sci. 2019;15(10):2265–81. pmid:31592238
- 21. Chen B, Yu J, Guo L, Byers MS, Wang Z, Chen X, et al. Circular RNA circHIPK3 Promotes the Proliferation and Differentiation of Chicken Myoblast Cells by Sponging miR-30a-3p. Cells. 2019; 8(2):177. pmid:30791438
- 22. Chen X, Ouyang H, Wang Z, Chen B, Nie Q. A Novel Circular RNA Generated by FGFR2 Gene Promotes Myoblast Proliferation and Differentiation by Sponging miR-133a-5p and miR-29b-1-5p. Cells. 2018;7(11): 199. pmid:30404220
- 23. Ouyang H, Chen X, Li W, Li Z, Nie Q, Zhang X. Circular RNA circSVIL Promotes Myoblast Proliferation and Differentiation by Sponging miR-203 in Chicken. Front Genet. 2018;9:172. pmid:29868120
- 24. Fu S, Zhao Y, Li Y, Li G, Chen Y, Li Z, et al. Characterization of miRNA transcriptome profiles related to breast muscle development and intramuscular fat deposition in chickens. J Cell Biochem. 2018; 119(8):7063–79. pmid:29737555
- 25. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9(4):357–9. pmid:22388286
- 26. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013; 495(7441):333–8. pmid:23446348
- 27. Gao Y, Zhang J, Zhao F. Circular RNA identification based on multiple seed matching. Brief Bioinform. 2018;19(5):803–10. pmid:28334140
- 28. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45. pmid:19541911
- 29. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. pmid:25516281
- 30. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14. pmid:20132535
- 31. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(Web Server issue):W316–22. pmid:21715386
- 32. Zhao J, Wu J, Xu T, Yang Q, He J, Song X. IRESfinder: Identifying RNA internal ribosome entry site in eukaryotic cell using framed k-mer features. J Genet Genomics. 2018; 45(7):403–6. pmid:30054216
- 33. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. pmid:11846609
- 34. Ebbesen KK, Kjems J, Hansen TB. Circular RNAs: Identification, biogenesis and function. Biochim Biophys Acta. 2016; 1859(1):163–8. pmid:26171810
- 35. Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159(1):134–47. pmid:25242744
- 36. Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol cell. 2014;56(1):55–66. pmid:25242144
- 37. Zhang XO, Dong R, Zhang Y, Zhang JL, Luo Z, Zhang J, et al. Diverse alternative back-splicing and alternative splicing landscape of circular RNAs. Genome Res. 2016;26(9):1277–87. pmid:27365365
- 38. Izawa T, Rohatgi N, Fukunaga T, Wang QT, Silva MJ, Gardner MJ, et al. ASXL2 Regulates Glucose, Lipid, and Skeletal Homeostasis. Cell Rep. 2015;11(10):1625–37. pmid:26051940
- 39. Park UH, Seong MR, Kim EJ, Hur W, Kim SW, Yoon SK, et al. Reciprocal regulation of LXRα activity by ASXL1 and ASXL2 in lipogenesis. Biochem Biophys Res Commun. 2014;443(2):489–94. pmid:24321552
- 40. Zhou C, King N, Chen KY, Breslow JL. Activation of PXR induces hypercholesterolemia in wild-type and accelerates atherosclerosis in apoE deficient mice. J Lipid Res. 2009; 50(10):2004–13. pmid:19436068
- 41. Payne VA, Au WS, Gray SL, Nora ED, Rahman SM, Sanders R, et al. Sequential regulation of diacylglycerol acyltransferase 2 expression by CAAT/enhancer-binding protein beta (C/EBPbeta) and C/EBPalpha during adipogenesis. J Biol Chem. 2007;282(29):21005–14. pmid:17504763
- 42. Kulcheski FR, Christoff AP, Margis R. Circular RNAs are miRNA sponges and can be used as a new class of biomarker. J Biotechnol. 2016; 238:42–51. pmid:27671698
- 43. Zhou J, Guo F, Wang G, Wang J, Zheng F, Guan X, et al. miR-20a regulates adipocyte differentiation by targeting lysine-specific demethylase 6b and transforming growth factor-β signaling. Int J Obes(Lond). 2015; 39(8):1282–91. pmid:25817070
- 44. Shen L, Li Q, Wang J, Zhao Y, Niu L, Bai L, et al. miR-144-3p Promotes Adipogenesis Through Releasing C/EBPα From Klf3 and CtBP2. Front Genet. 2018; 9:677. pmid:30619490
- 45. Shi C, Zhang M, Tong M, Yang L, Pang L, Chen L, et al. miR-148a is Associated with Obesity and Modulates Adipocyte Differentiation of Mesenchymal Stem Cells through Wnt Signaling. Sci Rep. 2015; 5:9930. pmid:26001136
- 46. Jin Y, Wang J, Zhang M, Zhang S, Lei C, Chen H, et al. Role of bta-miR-204 in the regulation of adipocyte proliferation, differentiation, and apoptosis. J Cell Physiol. 2019; 234(7):11037–46. pmid:30697738
- 47. Ling HY, Wen GB, Feng SD, Tuo QH, Ou HS, Yao CH, et al. MicroRNA-375 promotes 3T3-L1 adipocyte differentiation through modulation of extracellular signal-regulated kinase signalling. Clin Exp Pharmacol Physiol. 2011; 38(4):239–46. pmid:21291493
- 48. Sun T, Fu M, Bookout AL, Kliewer SA, Mangelsdorf DJ. MicroRNA let-7 regulates 3T3-L1 adipogenesis. Mol Endocrinol. 2009; 23(6):925–31. pmid:19324969
- 49. Li M, Liu Z, Zhang Z, Liu G, Sun S, Sun C. miR-103 promotes 3T3-L1 cell adipogenesis through AKT/mTOR signal pathway with its target being MEF2D. Biol Chem. 2015; 396(3):235–44. pmid:25400071
- 50. Zheng L, Lv GC, Sheng J, Yang YD. Effect of miRNA-10b in regulating cellular steatosis level by targeting PPAR-alpha expression, a novel mechanism for the pathogenesis of NAFLD. J Gastroenterol Hepatol. 2010; 25(1):156–63. pmid:19780876
- 51. Ahn J, Lee H, Jung CH, Jeon TI, Ha TY. MicroRNA-146b promotes adipogenesis by suppressing the SIRT1-FOXO1 cascade. EMBO Mol Med. 2013;5(10):1602–12. pmid:24009212
- 52. Pan S, Yang X, Jia Y, Li Y, Chen R, Wang M, et al. Intravenous injection of microvesicle-delivery miR-130b alleviates high-fat diet-induced obesity in C57BL/6 mice through translational repression of PPAR-γ. J Biomed Sci. 2015; 22:86. pmid:26475357
- 53. Chen C, Peng Y, Peng Y, Peng J, Jiang S. miR-135a-5p inhibits 3T3-L1 adipogenesis through activation of canonical Wnt/β-catenin signaling. J Mol Endocrinol. 2014; 52(3):311–20. pmid:24850830
- 54. Rockstroh D, Löffler D, Kiess W, Landgraf K, Körner A. Regulation of human adipogenesis by miR125b-5p. Adipocyte. 2016; 5(3):283–97. pmid:27617174
- 55. Huang S, Wang S, Bian C, Yang Z, Zhou H, Zeng Y, et al. Upregulation of miR-22 promotes osteogenic differentiation and inhibits adipogenic differentiation of human adipose tissue-derived mesenchymal stem cells by repressing HDAC6 protein expression. Stem Cells Dev. 2012;21(13):2531–40. pmid:22375943
- 56. Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, et al. miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006;3(2):87–98. pmid:16459310
- 57. Yang Z, Bian C, Zhou H, Huang S, Wang S, Liao L, et al. MicroRNA hsa-miR-138 inhibits adipogenic differentiation of human adipose tissue-derived mesenchymal stem cells through adenovirus EID-1. Stem Cells Dev. 2011; 20(2):259–67. pmid:20486779
- 58. Du WW, Liu F, Shan SW, Ma XC, Gupta S, Jin T, et al. Inhibition of Dexamethasone-induced Fatty Liver Development by Reducing miR-17-5p Levels. Mol Ther. 2015; 23(7):1222–33. pmid:25896250
- 59. Zhang M, Li CC, Li F, Li H, Liu XJ, Loor JJ, et al. Estrogen Promotes Hepatic Synthesis of Long-Chain Polyunsaturated Fatty Acids by Regulating ELOVL5 at Post-Transcriptional Level in Laying Hens. Int J Mol Sci. 2017; 18(7):1405. pmid:28665359
- 60. Zhao Q, Li S, Li N, Yang X, Ma S, Yang A, et al. miR-34a Targets HDAC1-Regulated H3K9 Acetylation on Lipid Accumulation Induced by Homocysteine in Foam Cells. J Cell Biochem. 2017; 118(12):4617–27. pmid:28485501
- 61. Davoodian N, Lotfi AS, Soleimani M, Ghaneialvar H. The combination of miR-122 overexpression and Let-7f silencing induces hepatic differentiation of adipose tissue-derived stem cells. Cell Biol Int. 2017; 41(10):1083–92. pmid:28792091
- 62. Zhong D, Huang G, Zhang Y, Zeng Y, Xu Z, Zhao Y, et al. MicroRNA-1 and microRNA-206 suppress LXRα-induced lipogenesis in hepatocytes. Cell Signal. 2013; 25(6):1429–37. pmid:23499676
- 63. Liu S, Yang Y, Wu J. TNFα-induced up-regulation of miR-155 inhibits adipogenesis by down-regulating early adipogenic transcription factors. Biochem Biophys Res Commun. 2011; 414(3):618–24. pmid:21986534
- 64. Jeong BC, Kang IH, Hwang YC, Kim SH, Koh JT. MicroRNA-194 reciprocally stimulates osteogenesis and inhibits adipogenesis via regulating COUP-TFII expression. Cell Death Dis. 2014; 5(11):e1532. pmid:25412310
- 65. Zhu L, Chen Y, Nie K, Xiao Y, Yu H. MiR-101 inhibits cell proliferation and invasion of pancreatic cancer through targeting STMN1. Cancer biomark. 2018; 23(2):301–9. pmid:30198871
- 66. Wang S, Wen X, Han XR, Wang YJ, Shen M, Fan SH, et al. MicroRNA-30d preserves pancreatic islet β-cell function through negative regulation of the JNK signaling pathway via SOCS3 in mice with streptozotocin-induced diabetes mellitus. J Cell Physiol. 2018; 233(9):7343–55. pmid:29663360
- 67. Cui S, Soni CB, Xie J, Li Y, Zhu H, Wu F, et al. MiR-30a-5p accelerates adipogenesis by negatively regulating Sirtuin 1. Int J Clin Exp Pathol. 2018; 11(11):5203–12. pmid:31949600
- 68. Xu B, Shen J, Li D, Ning B, Guo L, Bing H, et al. Overexpression of microRNA-9 inhibits 3T3-L1 cell adipogenesis by targeting PNPLA3 via activation of AMPK. Gene. 2020; 730:144260. pmid:31759991
- 69. Chen Y, Zhao Y, Jin W, Li Y, Zhang Y, Ma X, et al. MicroRNAs and their regulatory networks in Chinese Gushi chicken abdominal adipose tissue during postnatal late development. BMC Genomics. 2019; 20(1):778. pmid:31653195
- 70. Wang Y, Wang Z. Efficient backsplicing produces translatable circular mRNAs. RNA. 2015; 21(2):172–9. pmid:25449546
- 71. Chen CY, Sarnow P. Initiation of protein synthesis by the eukaryotic translational apparatus on circular RNAs. Science. 1995; 268(5209):415–7. pmid:7536344
- 72. Wesselhoeft RA, Kowalski PS, Anderson DG. Engineering circular RNA for potent and stable translation in eukaryotic cells. Nat Commun. 2018; 9(1):2629. pmid:29980667
- 73. Taulli R, Loretelli C, Pandolfi PP. From pseudo-ceRNAs to circ-ceRNAs: a tale of cross-talk and competition. Nat Struct Mol Biol. 2013; 20(5):541–3. pmid:23649362
- 74. Kartha RV, Subramanian S. Competing endogenous RNAs (ceRNAs): new entrants to the intricacies of gene regulation. Front Genet. 2014; 5:8. pmid:24523727