A genome-wide landscape of mRNAs, lncRNAs, and circRNAs during subcutaneous adipogenesis in pigs

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in Small Tail Han sheep with FecB++ genotype

The thyroid gland is an important endocrine gland in animals, which mainly secretes thyroid hormones and acts on various organs of the body. Long-chain non-coding RNA (lncRNA) plays an important role in animal reproduction. However, there is still a lack of understanding of their expression patterns and potential roles in the thyroid of Small Tail Han (STH) sheep. In this study, RNA-seq was used to examine the transcriptome expression patterns of lncRNAs and mRNAs in the follicular phase (ww_FT) and luteal phase (ww_LT) in FecB++ genotype STH Sheep. A total of 17,217 lncRNAs and 39,112 mRNAs were identified including 96 differentially expressed lncRNAs (DELs) and 1054 differentially expressed mRNAs (DEGs). Functional analysis of genes with significant differences in expression level showed that these genes could be enriched in Ras signalling pathway, hedgehog (HH) signalling pathway, ATP-binding cassette (ABC) transporters and other signalling pathways related to animal reproduction. In addition, through correlation analysis for lncRNA-mRNA co-expression and network construction, we found that LNC_009115 and LNC_005796 trans target NIK-related kinase (NRK) and poly(A)-specific ribonuclease (PARN). LNC_007189 and LNC_002045 trans target progesterone-induced blocking factor 1 (PIBF1), LNC_009013 trans targets small mothers against decapentaplegic (SMAD1) are related to animal reproduction. These genes add new resources for elucidating the regulatory mechanisms of reproduction in sheep with different reproductive cycles of the FecB++ genotype STH sheep.

Thyroid transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in small tail han sheep with FecB BB genotype

BACKGROUND

The thyroid gland is an important endocrine gland in animals that secretes thyroid hormones and acts on various organs throughout the body. lncRNAs are long non-coding RNAs that play an important role in animal reproduction; however, there is a lack of understanding of their expression patterns and potential roles in the thyroid gland of the Small Tail Han (STH) sheep. In this study, we used RNA-Seq technology to examine the transcriptome expression pattern of the thyroid from the luteal phase (LP) and follicular phase (FP) of FecB BB (MM) STH sheep.

RESULTS

We identified a total of 122 and 1287 differential expression lncRNAs (DELs) and differential expression mRNAs (DEGs), respectively, which were significantly differentially expressed. These DELs target genes and DEGs can be enriched in several signalling pathways related to the animal reproduction process.

CONCLUSIONS

The expression profiles of DELs and DEGs in thyroid glands provide a more comprehensive resource for elucidating the reproductive regulatory mechanisms of STH sheep.

LncRNA BlncAD1 Modulates Bovine Adipogenesis by Binding to MYH10, PI3K/Akt Signaling Pathway, and miR-27a-5p/CDK6 Axis

Research on adipogenesis will help to improve the meat quality of livestock. Long noncoding RNAs (lncRNAs) are involved in mammalian adipogenesis as epigenetic modulators. In this study, we analyzed lncRNA expression during bovine adipogenesis and detected 195 differentially expressed lncRNAs, including lncRNA BlncAD1, which was significantly upregulated in mature bovine adipocytes. Gain- and loss-of-function experiments confirmed that BlncAD1 promoted the proliferation, apoptosis, and differentiation of bovine preadipocytes. RNA pull-down revealed that the nonmuscle myosin 10 (MYH10) is a potential binding protein of BlncAD1. Then, we elucidated that loss of BlncAD1 caused increased ubiquitination of MYH10, which confirmed that BlncAD1 regulates adipogenesis by enhancing the stability of the MYH10 protein. Western blotting was used to demonstrate that BlncAD1 activated the PI3K/Akt signaling pathway. Bioinformatic analysis and dual-luciferase reporter assays indicated that BlncAD1 competitively absorbed miR-27a-5p. The overexpression and interference of miR-27a-5p in bovine preadipocytes displayed that miR-27a-5p inhibited proliferation, apoptosis, and differentiation. Further results suggested that miR-27a-5p targeted the CDK6 gene and that BlncAD1 controlled the proliferation of bovine preadipocytes by modulating the miR-27a-5p/CDK6 axis. This study revealed the complex mechanisms of BlncAD1 underlying bovine adipogenesis for the first time, which would provide useful information for genetics and breeding improvement of Chinese beef cattle.

Whole-transcriptome analysis of longissimus dorsi muscle in cattle-yaks reveals the regulatory functions of ADAMTS6 gene in myoblasts

Cattle-yak, which is the hybrid F1 generation of cattle and yak, demonstrates better production performance compared to yak. However, there is limited research on the molecular mechanisms responsible for the muscle development of cattle-yak. To address this knowledge gap, a comprehensive transcriptomic survey of the longissimus dorsi muscle in cattle-yak was conducted. Three transcript types, namely lncRNAs, miRNAs, and circRNAs, along with protein-coding genes were characterized at two developmental stages (6 m, 18 m) of cattle-yak. The results revealed significant enrichment of these transcripts into pathways related to myoblast differentiation and muscle development signaling. Additionally, the study identified the TCONS00024465/circHIPK3-bta-miR-499-ADAMTS6 regulatory network, which may play a crucial role in the muscle development of cattle-yak by combining the transcriptome data of yak and constructing the ceRNA co-expression network. HEK 293 T cells were used to validate that TCONS00024465 and circHIPK3 are located upstream of bta-miR-499, and can competitively bind to bta-miR-499 as ceRNA. The study also verified that ADAMTS6 regulates skeletal muscle development by inhibiting myoblast proliferation, promoting myoblast differentiation, and positively regulating the apoptosis of myoblasts. Taken together, this study provides new insights into the advantages of cattle-yak production performance and offers a molecular basis for further research on muscle development.

Expression Analysis of circRNAs in Human Adipogenesis

Purpose

Adipogenesis is one of the major pathways for generating obesity or overweight that can cause a range of metabolic disorders. Circular RNAs (circRNAs), a specific type of RNAs, have a significant influence on metabolic disorders. This study aims to find differentially expressed circRNAs (DECs) during human subcutaneous adipose tissue (SATs) adipogenesis.

Patients and Methods

The human adipose tissue-derived stromal cells (hADSCs) were isolated from human SATs (n = 3), and then induced into adipocytes. Total RNAs were extracted from hADSCs and adipocytes, and he DECs were detected using circRNA microarray. The GO and KEGG pathways of DECs were analyzed by bioinformatic methods, and partial DECs were further validated by quantitative polymerase chain reaction (qPCR).

Results

Our study detected a total of 1987 DECs, among which, 1134 were found upregulated and 853 were downregulated. GO analysis showed that the upregulated DECs have catalytic activity in intracellular organelle and cytoplasms, whereas downregulated DECs are enriched in organelle lumen, and are involved in positive regulation of developmental process. In addition, pathway results demonstrated that upregulated DECs are involved in platinum drug resistance and cellular senescence, and downregulated DECs are enriched in proteoglycans in cancer and focal adhesion pathway. Two circRNAs, namely has_circ_0001600 and has_circ_0001947 were validated to be significantly upregulated in adipocytes compared to hADSCs.

Conclusion

Our study explored DECs between hADSCs derived from SATs and adipocytes, and report that two circRNAs named has_circ_0001600 and has_circ_0001947 might be important factors involved in human adipogenesis, however, the molecular mechanism should be further explored.

Dietary dibutyryl cAMP supplementation regulates the fat deposition in adipose tissues of finishing pigs via cAMP/PKA pathway

This study investigated potential mechanism of dibutyryl-cAMP (db-cAMP) on porcine fat deposition. (1) Exp.1, 72 finishing pigs were allotted to 3 treatments (0, 10 or 20 mg/kg dbcAMP) with 6 replicates. dbcAMP increased the hormone sensitive lipase (HSL) activity and expression of β-adrenergic receptor (β-AR) and growth hormone receptor (GHR), but decreased expression of peroxisome proliferator-activated receptor gamma 2 (PPAR-γ2) and adipocyte fatty acid binding protein (A-FABP) in back fat. dbcAMP upregulated expression of β-AR, GHR, PPAR-γ2 and A-FABP, but decreased insulin receptor (INSR) expression in abdominal fat. Dietary dbcAMP increased HSL activity and expression of G protein-coupled receptor (GPCR), cAMP-response element-binding protein (CREB) and insulin-like growth factor-1 (IGF-1), but decreased fatty acid synthase (FAS) and lipoprotein lipase (LPL) activities, and expression of INSR, cAMP-response element-binding protein (C/EBP-α) and A-FABP in perirenal fat. (2) Exp. 2, dbcAMP suppressed the proliferation and differentiation of porcine preadipocytes in a time- and dose-dependent manner, which might be associated with increased activities of cAMP and protein kinase A (PKA), and expression of GPCR, β-AR, GHR and CREB via inhibiting C/EBP-α and PPAR-γ2 expression. Collectively, dbcAMP treatment may reduce fat deposition by regulating gene expression related to adipocyte differentiation and fat metabolism partially via cAMP-PKA pathway.

The emerging roles of circRNAs in traits associated with livestock breeding

Many indicators can be used to evaluate the productivity and quality of livestock, such as meat and milk production as well as fat deposition. Meat and milk production are measures of livestock performance, while fat deposition affects the taste and flavor of the meat. The circRNAs, are non-coding RNAs, that are involved in the regulation of all these three traits. We review the functions and mechanisms of circRNAs in muscle and fat development as well as lactation to provide a theoretical basis for circRNA research in animal husbandry. Various phenotypic changes presented in livestock may be produced by different circRNAs. Our current concern is how to use the roles played by circRNAs to our advantage to produce the best possible livestock. Hence, we describe the advantages and disadvantages of knockout techniques for circRNAs. In addition, we also put forward our thoughts regarding the mechanism and network of circRNAs to provide researchers with novel ideas of how molecular biology can help us advance our goals in animal farming. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Transcriptional Regulation Associated with Subcutaneous Adipogenesis in Porcine ACSL1 Gene

Long-chain acyl-CoA synthetase 1 (ACSL1) plays an important role in fatty acid metabolism and fat deposition. The transcription of the ACSL1 gene is regulated specifically among cells and physiological processes, and transcriptional regulation of ACSL1 in adipogenesis remains elusive. Here, we characterize transcription factors (TFs) associated with adipogenesis in the porcine ACSL1 gene. CCAAT-enhancer binding protein (C/EBP)α, a well-known adipogenic marker, was found to enhance the expression of the ACSL1 gene via binding two tandem motifs in the promoter. Further, we demonstrate that ACSL1 mediates C/EBPα effects on adipogenesis in preadipocytes cultured from subcutaneous fat tissue of pigs via gain- and loss-of-function analyses. The cAMP-response element binding protein, another TF involved in adipogenesis, was also identified in the regulation of ACSL1 gene expression. Additionally, single nucleotide polymorphisms (SNPs) were screened in the promoter of ACSL1 among four breeds including the Chinese indigenous Min, and Duroc, Berkshire, and Yorkshire pigs through sequencing of PCR products. Two tightly linked SNPs, -517G>T and -311T>G, were found exclusively in Min pigs. The haplotype mutation decreases promoter activity in PK-15 and ST cells, and in vivo the expression of ACSL1, illustrating a possible role in adipogenesis regulated by C/EBPα/ACSL1 axis. Additionally, a total of 24 alternative splicing transcripts were identified, indicating the complexity of alternative splicing in the ACSL1 gene. The results will contribute to further revealing the regulatory mechanisms of ACSL1 during adipogenesis and to the characterization of molecular markers for selection of fat deposition in pigs.

Comprehensive analysis of miRNAs, lncRNAs and mRNAs profiles in backfat tissue between Daweizi and Yorkshire pigs

OBJECTIVE

Daweizi (DWZ) is a famous indigenous pig breed in China and characterized by tender meat and high fat percentage. However, the expression profiles and functions of transcripts in DWZ pigs is still in infancy. The object of this study was to depict the transcript profiles in DWZ pigs and screen the potential pathway influence adipogenesis and fat deposition.

METHODS

Histological analysis of backfat tissue was firstly performed between DWZ and lean-type Yorkshire pigs, and then RNA sequencing technology was utilized to explore miRNAs, lncRNAs and mRNAs profiles in backfat tissue. 18 differentially expressed (DE) transcripts were randomly selected for quantitative real-time polymerase chain reaction (QPCR) to validate the reliability of the sequencing results. Finally, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis were conducted to investigate the potential pathways influence adipocyte differentiation, adipogenesis and lipid metabolism, and a schematic model was further proposed.

RESULTS

A total of 1,625 differentially expressed transcripts were identified in DWZ pigs, including 27 upregulated and 45 downregulated miRNAs, 64 upregulated and 119 downregulated lncRNA, 814 upregulated and 556 downregulated mRNAs. QPCR analysis exhibited strong consistency with the sequencing data. GO and KEGG analysis elucidated that the differentially expressed transcripts were mainly associated with cell growth and death, signal transduction, peroxisome proliferator-activated receptors (PPAR), AMP-activated protein kinase (AMPK), PI3K-Akt, adipocytokine and foxo signaling pathways, all of which are strongly involved in cell development, lipid metabolism and adipogenesis. Further analysis indicated that the BGIR9823_87926/miR-194a-5p/AQP7 network may be effective in the process of adipocyte differentiation or adipogenesis.

CONCLUSION

Our study provides comprehensive insights into the regulatory network of backfat deposition and lipid metabolism in pigs from the point of view of miRNAs, lncRNAs and mRNAs.

Circular RNA- and microRNA-Mediated Post-Transcriptional Regulation of Preadipocyte Differentiation in Adipogenesis: From Expression Profiling to Signaling Pathway

Adipogenesis is an indispensable cellular process that involves preadipocyte differentiation into mature adipocyte. Dysregulated adipogenesis contributes to obesity, diabetes, vascular conditions and cancer-associated cachexia. This review aims to elucidate the mechanistic details on how circular RNA (circRNA) and microRNA (miRNA) modulate post-transcriptional expression of targeted mRNA and the impacted downstream signaling and biochemical pathways in adipogenesis. Twelve adipocyte circRNA profiling and comparative datasets from seven species are analyzed using bioinformatics tools and interrogations of public circRNA databases. Twenty-three circRNAs are identified in the literature that are common to two or more of the adipose tissue datasets in different species; these are novel circRNAs that have not been reported in the literature in relation to adipogenesis. Four complete circRNA-miRNA-mediated modulatory pathways are constructed via integration of experimentally validated circRNA-miRNA-mRNA interactions and the downstream signaling and biochemical pathways involved in preadipocyte differentiation via the PPARγ/C/EBPα gateway. Despite the diverse mode of modulation, bioinformatics analysis shows that the circRNA-miRNA-mRNA interacting seed sequences are conserved across species, supporting mandatory regulatory functions in adipogenesis. Understanding the diverse modes of post-transcriptional regulation of adipogenesis may contribute to the development of novel diagnostic and therapeutic strategies for adipogenesis-associated diseases and in improving meat quality in the livestock industries.

Recent Progress on Circular RNAs in the Development of Skeletal Muscle and Adipose Tissues of Farm Animals

Circular RNAs (circRNAs) are a highly conserved and specifically expressed novel class of covalently closed non-coding RNAs. CircRNAs can function as miRNA sponges, protein scaffolds, and regulatory factors, and play various roles in development and other biological processes in mammals. With the rapid development of high-throughput sequencing technology, thousands of circRNAs have been discovered in farm animals; some reportedly play vital roles in skeletal muscle and adipose development. These are critical factors affecting meat yield and quality. In this review, we have highlighted the recent advances in circRNA-related studies of skeletal muscle and adipose in farm animals. We have also described the biogenesis, properties, and biological functions of circRNAs. Furthermore, we have comprehensively summarized the functions and regulatory mechanisms of circRNAs in skeletal muscle and adipose development in farm animals and their effects on economic traits such as meat yield and quality. Finally, we propose that circRNAs are putative novel targets to improve meat yield and quality traits during animal breeding.

Identification and characterization of extrachromosomal circular DNA in Wei and Large White pigs by high-throughput sequencing

Wei pig (WP) and Large White pig (LP) are fatty and lean breeds, respectively. Extrachromosomal circular DNA (eccDNA) plays an important role in regulating signaling pathway processes of cell. However, there are few reports regarding the eccDNA and ecDNA profiles in WP and LP. The present work aimed to investigate the eccDNA and ecDNA profiles between WP and LP. Three WPs and three LPs (100 ± 1.3 kg) were selected for analysis of eccDNA and ecDNA in the ear samples. Results showed that there were 39,686,953,656-58,411,217,258 and 53,824,168,657-58,311,810,737 clean data for WP and LP, respectively. Sequencing yielded 15,587-25,479 and 71,123-79,605 eccDNAs from the ear samples of WP and LP, respectively. There were 15,111 and 22,594 eccDNA-derived genes in the WP and LP, respectively, and 13,807 eccDNA-derived genes were common in the ear samples of both pigs. Sequencing yielded 13-19 and 27-43 ecDNAs in the ears of WP and LP, respectively. There were 1,005 and 1,777 ecDNA-derived genes in WP and LP, respectively, and 351 ecDNA-derived genes were common in the ear samples of both pigs. The most significant KEGG pathways of eccDNA-derived genes were axon guidance, focal adhesion, metabolic pathways, MAPK signaling pathway, Hedgehog signaling pathway, microRNAs in cancer, tight junction, phospholipase D signaling pathway, endocytosis, and sphingolipid signaling pathway. Furthermore, the most significant KEGG pathways of ecDNA-derived genes were olfactory transduction, B cell receptor signaling pathway, and chemical carcinogenesis. The eccDNA00044301 was lower abundance, while the ecDNA00000060 was higher abundance in WP compared with that in LP. Summary, we found that eccDNAs and ecDNAs are common in WP and LP and occur in sizes large enough to carry one or several partial or complete genes. These findings have expanded the knowledge repertoire of circular DNA in pig and will provide a reference for the use of pigs as a medical model and help discovery of new genetic markers to select high-quality breeds.

Genomic and transcriptomic analyses enable the identification of important genes associated with subcutaneous fat deposition in Holstein cows

Subcutaneous fat deposition has many important roles in dairy cattle, including immunological defense and mechanical protection. The main objectives of this study are to identify key candidate genes regulating subcutaneous fat deposition in high-producing dairy cows by integrating genomic and transcriptomic datasets. A total of 1,654 genotyped Holstein cows are used to perform a genome-wide association study (GWAS) aiming to identify genes associated with subcutaneous fat deposition. Subsequently, weighted gene co-expression network analyses (WGCNA) are conducted based on RNA-sequencing data of 34 cows and de-regressed estimated breeding values of subcutaneous fat deposition. Lastly, differentially expressed (DE) mRNA, lncRNA, and differentially alternative splicing genes are obtained for 12 Holstein cows with extreme and divergent phenotypes for subcutaneous fat deposition. Forty-six protein-coding genes are identified as candidate genes regulating subcutaneous fat deposition in Holstein cattle based on the GWAS. Eleven overlapping genes are identified based on the analyses of DE genes and WGCNA. Furthermore, the candidate genes identified based on the GWAS, WGCNA, and analyses of DE genes are significantly enriched for pathways involved in metabolism, oxidative phosphorylation, thermogenesis, fatty acid degradation, and glycolysis/gluconeogenesis pathways. Integrating all findings, the NID2, STARD3, UFC1, DEDD, PPP1R1B, and USP21 genes are considered to be the most important candidate genes influencing subcutaneous fat deposition traits in Holstein cows. This study provides novel insights into the regulation mechanism underlying fat deposition in high-producing dairy cows, which will be useful when designing management and breeding strategies.

Differences in histomorphology and expression of key lipid regulated genes of four adipose tissues from Tibetan pigs

Tibetan pigs, an indigenous pig breed in China, have high overall fat deposition and flavorful and tasty meat. They are thus good models for studying adipogenesis. Few studies have been conducted focusing on expression of lipid regulated genes in different adipose tissues of Tibetan pigs. Therefore, we compared the difference of histomorphology and expression level of lipid regulated genes through qPCR and western blot in subcutaneous fat, perirenal fat, omental adipose tissue, and inguinal fat of Tibetan pigs. Our results showed that the area of subcutaneous adipocytes in Tibetan pigs was smaller, while the other three adipose tissues (perirenal fat, greater omentum fat, inguinal fat) had cell areas of similar size. The gene expression of FABP4, FASN, FABP3, and ME1 in subcutaneous fat was significantly higher than that in perirenal fat. Furthermore, the protein expression of FABP4 was significantly lower in subcutaneous fat than in perirenal fat (p < 0.05), and the expression of FASN was higher in greater omentum fat than in subcutaneous fat (p = 0.084). The difference in adipocyte cell size and expression of lipid-regulated genes in adipose tissues from the various parts of the pig body is likely due to the different cellular lipid metabolic processes. Specially, FABP4 and FASN may be involved in the regulation of fat deposition in different adipose tissues of Tibetan pigs.

Comprehensive Analysis of the lncRNA-miRNA-mRNA Regulatory Network for Intramuscular Fat in Pigs

Intramuscular fat (IMF) is an essential trait closely related to meat quality. The IMF trait is a complex quantitative trait that is regulated by multiple genes. In order to better understand the process of IMF and explore the key factors affecting IMF deposition, we identified differentially expressed mRNA, miRNA, and lncRNA in the longissimus dorsi muscle (LD) between Songliao Black (SL) pigs and Landrace pigs. We obtained 606 differentially expressed genes (DEGs), 55 differentially expressed miRNAs (DEMs), and 30 differentially expressed lncRNAs (DELs) between the SL pig and Landrace pig. Enrichment results from GO and KEGG indicate that DEGs are involved in fatty acid metabolism and some pathways related to glycogen synthesis. We constructed an lncRNA-miRNA-mRNA interaction network with 18 DELs, 11 DEMs, and 42 DEGs. Finally, the research suggests that ARID5B, CPT1B, ACSL1, LPIN1, HSP90AA1, IRS1, IRS2, PIK3CA, PIK3CB, and PLIN2 may be the key genes affecting IMF deposition. The LncRNAs MSTRG.19948.1, MSTRG.13120.1, MSTRG.20210.1, and MSTRG.10023.1, and the miRNAs ssc-miRNA-429 and ssc-miRNA-7-1, may play a regulatory role in IMF deposition through their respective target genes. Our research provides a reference for further understanding the regulatory mechanism of IMF.

The landscape of the long non-coding RNAs and circular RNAs of the abdominal fat tissues in the chicken lines divergently selected for fatness

BACKGROUND

Excessive deposition of abdominal fat poses serious problems in broilers owing to rapid growth. Recently, the evolution of the existing knowledge on long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) have established their indispensable roles in multiple physiological metabolic processes, including adipogenesis and fat deposition. However, not much has been explored on their profiles in the abdominal fat tissues of broilers to date. In the study, we aimed to characterize the vital candidates of lncRNAs and circRNAs and their underlying regulations for abdominal fat deposition in broilers.

RESULTS

The present study sequenced the lncRNAs and circRNAs expression profiles in the abdominal fat tissues isolated from 7-week-old broilers, who were divergently selected for their fatness. It identified a total of 3359 lncRNAs and 176 circRNAs, demonstrating differential expressed (DE) 30 lncRNAs and 17 circRNAs between the fat- and lean-line broilers (|log2FC| ≥ 1, P < 0.05). Subsequently, the 20 cis-targets and 48 trans-targets of the candidate DE lncRNAs were identified for depositing abdominal fat by adjacent gene analysis and co-expression analysis, respectively. In addition, the functional enrichment analysis showed the DE lncRNAs targets and DE circRNAs host genes to be mainly involved in the cellular processes, amino/fatty acid metabolism, and immune inflammation-related pathways and GO terms. Finally, the vital 16 DE lncRNAs located in cytoplasm and specifically expressed in fat/lean line and their targets were used to construct the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) regulatory network, comprising 7 DE lncRNAs, 28 miRNAs, 11 DE mRNAs. Notably, three lncRNAs including XR_001468036.2, XR_003077610.1 and XR_001466431.2 with the most connected degrees might play hub regulatory roles in abdominal fat deposition of broilers.

CONCLUSIONS

This study characterized the whole expression difference of lncRNAs and circRNAs between the two lines broilers with divergently ability of abdominal fat. The vital candidate DE lncRNAs/circRNAs and ceRNA regulations were identified related to the deposition of abdominal fat in chicken. These results might further improve our understanding of regulating the non-coding RNAs in obesity.

Transcriptome-Wide Study of mRNAs and lncRNAs Modified by m6A RNA Methylation in the Longissimus Dorsi Muscle Development of Cattle-Yak

Cattle-yak is a hybrid F₁ generation of cattle and yak, which has a history of more than 3000 years and has shown better production performance and higher economic benefits than those of yaks. However, up to now, there has been no study on the transcriptome-wide m⁶A methylation profile of bovine skeletal muscle and its potential biological function during muscle development. Here, we observed significant changes in the expression levels of muscle-related marker genes and methylation-related enzymes during the development of cattle-yak, and the overall m⁶A content in the Longissimus dorsi muscle of 18-month-old cattle-yak decreased significantly. A total of 36,602 peaks, 11,223 genes and 8388 lncRNAs were identified in the two groups, including 2989 differential peaks (427 up-regulated peaks and 2562 down-regulated peaks), 1457 differentially expressed genes (833 up-regulated genes and 624 down-regulated genes) and 857 differentially expressed lncRNAs (293 up-regulated lncRNAs and 564 down-regulated lncRNAs). GO and KEGG analysis revealed that they were significantly enriched in some muscle-related pathways (Wnt signaling pathway and MAPK signaling pathway) and high-altitude adaptation-related pathway (HIF-1 signaling pathway). Moreover, m⁶A abundance was positively correlated with gene expression levels, while it was negatively correlated with lncRNA expression levels. This indicates that m⁶A modification played an important role in the Longissimus dorsi muscle development of cattle-yak; however, the regulation mechanism of m⁶A-modified mRNA and lncRNA may be different. This study was the first report of transcriptome-wide m⁶A-modified mRNAs and lncRNAs atlas in the Longissimus dorsi muscle development of cattle-yak, one which will provide new perspectives for genetic improvement in bovines.

A genome-wide landscape of mRNAs, lncRNAs, circRNAs and miRNAs during intramuscular adipogenesis in cattle

Background

Intramuscular preadipocyte differentiation plays a critical role in bovine intramuscular fat (IMF) deposition. However, the roles of different RNAs, including mRNAs, circRNAs, lncRNAs and miRNAs, in regulating the adipogenic differentiation of intramuscular preadipocytes remain largely unclear.

Results

In the present study, a whole transcriptome sequencing and analysis, including the analysis of mRNAs, circRNAs, lncRNAs and miRNAs, during different differentiation stages (0, 3, 6, and 9 d) of intramuscular preadipocytes from Qinchuan cattle was performed. All samples were prepared with 3 biological replicates. Here, a total of 27,153 mRNAs, 14,070 circRNAs, 7035 lncRNAs, and 427 miRNAs were annotated. Among them, we identified 4848 differentially expressed mRNAs (DEMs), 181 DE circRNAs (DECs), 501 DE lncRNAs (DELs) and 77 DE miRNAs (DEmiRs) between 0 d and other differentiation days (3, 6, and 9 d). GO and KEGG functional enrichment analyses showed that these differentially expressed genes were mainly enriched in cell differentiation, fat metabolism and adipogenesis-related pathways. Furthermore, weighted gene coexpression network analysis (WGCNA) and co-expression network analysis screened out multiple important mRNAs, circRNAs and lncRNAs related to intramuscular adipogenesis. Based on the competing endogenous RNA (ceRNA) regulatory mechanism, we finally identified 24 potential ceRNA networks and 31 potential key genes, including FOXO1/miR-330/circRNA2018/MSTRG.20301, GPAM/miR-27b/ciRNA489 and SESN3/miR-433/circRNA2627MSTRG.20342.

Conclusions

This study provides new insights into the differential expression patterns of different transcript types (i.e., mRNAs, circRNAs, lncRNAs and miRNAs) in intramuscular preadipocyte differentiation. Our findings provide data support for studying the molecular mechanism of key mRNAs and noncoding RNAs in IMF deposition, and provide new candidate markers for the molecular breeding of beef cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08911-z.

Comparative proteomics reveals genetic mechanisms of body weight in Hu sheep and Dorper sheep

Body weight (BW) is a critical economic trait for meat production in sheep, and it is a complex trait affected by numerous elements. The aim of this study was to investigate the genetic mechanisms of sheep BW by a label-free proteomics approach. The result showed, a total of 27, 14, 61, and 65 differentially abundant proteins (DAPs) were identified in the Hu_HBW vs. Hu_LBW, DP_HBW vs. DP_LBW, Hu_HBW vs. DP_HBW, and Hu_LBW vs. DP_LBW comparisons, respectively. Five proteins (including ILK, AHCYL2, MLIP, CYB5A, and SMTNL1) related to fat synthesis and muscle development were detected in the Hu sheep group. In the Dorper sheep group, the screened DAPs strictly related to muscle development and fat synthesis were significantly enriched in MAP kinase activity (MAPK12), Arachidonic acid metabolism, and Steroid hormone biosynthesis (PGFS, LOC101107119) pathways. Several DAPs related to immune responses (SERPINA1, FGG, SERPINC1, and LOC101108131), fat deposition (APOH, GC, AHSG, SKP1, ACSL1, ACAT1, and ACADS), and muscle development (LMOD3 and LRRC39) were detected in the Hu vs. Dorper sheep comparison. These analyses indicated that the BW of sheep is regulated via a variety of pathways, and these DAPs can be further investigated as candidate markers for predicting the BW of sheep. SIGNIFICANCE: Body weight is one of the key traits in sheep and involves multiple coordinated regulatory mechanisms, but the genetic mechanism of BW is still unclear in sheep. In the current study, the label-free method was used to identify the proteins and pathways related to BW using LT muscle of Hu sheep and Dorper sheep with different BW. These findings will provide new candidate proteins and vital pathways into the molecular mechanisms involved growth traits in sheep.

C/EBPα promotes porcine pre-adipocyte proliferation and differentiation via mediating MSTRG.12568.2/FOXO3 trans-activation for STYX

As a key adipogenic marker, C/EBPα (CCAAT/enhancer binding protein α) is also an important factor in regulating targets containing CCAAT element for transcription, whose products include coding and non-coding RNAs (ncRNAs). However, knowledge of the mechanism of C/EBPα affecting pre-adipocyte proliferation and adipogenesis through regulating ncRNA is still limited. In this study, we firstly conducted an investigation concerning the impact of C/EBPα knockdown on porcine pre-adipocytes by using RNA sequencing (RNA-Seq) to identify the role of key ncRNAs, especially lncRNAs and their correlated mRNAs in regulating proliferation and differentiation of porcine pre-adipocytes. 97 differentially expressed (DE) mRNAs and 4 DE lncRNAs were identified in si-C/EBPα groups compared with the si-NC groups. Meanwhile, we found C/EBPα directly target the promoter of a novel lncRNA, namely MSTRG.12568.2, which was trans-correlated with STYX (serine/threonine/tyrosine interacting protein), an important candidate gene for regulating cell proliferation. Moreover, FOXO3 (forkhead box O3) was identified as a co-regulator with MSTR.12568.2 for STYX. Overexpression and knockdown of any of the MSTRG.12568.2, STYX, and FOXO3 increased and decreased the levels of pre-adipocyte proliferation and differentiation, respectively, which demonstrated that they played a positive role in adipogenesis of pre-adipocytes. Moreover, our results revealed that FOXO3 was necessary for MSTRG.12568.2 to trans-activate STYX. We revealed that C/EBPα regulated pre-adipocyte proliferation and differentiation through mediating trans-activation of MSTRG.12568.2-FOXO3 to STYX. These results provide a novel regulation signal for C/EBPα to influence porcine pre-adipocyte proliferation and differentiation and greatly benefit to our understanding of molecular mechanism regulating subcutaneous adipogenesis.

Weighted Gene Co-expression Network Analysis Revealed That CircMARK3 Is a Potential CircRNA Affects Fat Deposition in Buffalo

Background

Buffalo meat is increasingly widely accepted for consumption as it shares several quality attributes with cattle meat (beef). Hence, there is a huge opportunity for growth in the buffalo meat industry. However, buffalo meat has relatively low intramuscular fat (IMF) content, affecting its flavor, tenderness and juiciness. As there is a dearth of information on factors that control fat deposition, this study was undertaken to provide new candidate factor associated with buffalo fat deposition. Circular RNA (circRNA) is a novel class of non-coding RNA with a closed-loop structure, and play an important role in fat deposition.

Methods

In this study, weighted gene co-expression network analysis (WGCNA) was used to construct a circRNA co-expression network and revealed a candidate circRNA that may affect the IMF deposition of buffalo as determined by RT-qPCR, semiquantitative PCR and gain-of-function experiments.

Results

Herein, WGCNA determined that one module (turquoise module) is significantly associated with the growth and development stages of buffalo. Further analysis revealed a total of 191 overlapping circRNAs among differentially expressed (DE) circRNAs and the co-expression module. A candidate circRNA was found, 21:6969877|69753491 (circRNA_ID), with a reported involvement in lipid metabolism. This circRNA is stably expressed and originates from the MARK3 gene, hence the name circMARK3. circMARK3 is highly expressed in adipose tissue and mature adipocytes and is located in the cytoplasm. Gain-of-function experiments demonstrated that circMARK3 promoted adipogenic differentiation of buffalo adipocytes and 3T3-L1 cells by up-regulating the expression levels of adipogenic marker genes PPARG, C/EBPα and FABP4.

Conclusion

These results indicate that circMARK3 is a potential factor that promotes fat deposition by regulating adipocyte differentiation and adipogenesis in buffalo.

Sus_circPAPPA2 Regulates Fat Deposition in Castrated Pigs through the miR-2366/GK Pathway

CircRNAs play an important role in fat deposition, and testosterone-deficient boars exhibit significantly increased fat deposition; however, the mechanism by which testosterone regulates fat deposition through circRNAs remains unclear. In this study, circRNA-seq of backfat and abdominal fat from castrated and intact full-sib Yorkshire pigs was performed. The GO and KEGG enrichment analyses revealed that the host genes of the dorsal DE circRNAs were mainly involved in fatty acid transport, while in abdominal tissues, these genes were mainly involved in adipogenesis and inflammation. The interaction among sus_circPAPPA2, ssc-miR-2366 and GK was verified by dual fluorescence experiments and in porcine preadipocytes. The overexpression of sus_circPAPPA2 significantly inhibited the differentiation of preadipocytes. The expression of sus_circPAPPA2 was increased after adding 100 nM of testosterone, and preadipocyte differentiation was significantly inhibited. Testosterone can affect preadipocyte differentiation by upregulating the expression of sus_circPAPPA2, sponging miR-2366 and regulating the expression of genes, such as GK. These results indicate that testosterone can regulate the expression of adipocyte differentiation- and lipid metabolism-related genes by regulating the expression of circRNA, and ceRNA networks are different in the testosterone regulation of adipose deposition in different parts. This study provides basic data enhancing the understanding of the interaction between the hormone environment and mir-2366/GK to regulate trait performance in pigs.

Thyroid Transcriptomic Profiling Reveals the Follicular Phase Differential Regulation of lncRNA and mRNA Related to Prolificacy in Small Tail Han Sheep with Two FecB Genotypes

Long non-coding RNA (lncRNA) accounts for a large proportion of RNA in animals. The thyroid gland has been established as an important gland involved in animal reproduction, however, little is known of its gene expression patterns and potential roles in the sheep. Herein, RNA-Seq was used to detect reproduction-related differentially expressed lncRNAs (DELs) and mRNAs (DEGs) in the follicular phase (FT) FecBBB (MM) and FecB++ (ww) genotypes of Small Tail Han (STH) sheep thyroids. Overall, 29 DELs and 448 DEGs in thyroid between MM and ww sheep were screened. Moreover, GO and KEGG enrichment analysis showed that targets of DELs and DEGs were annotated in biological transitions, such as cell cycle, oocyte meiosis and methylation, which in turn affect reproductive performance in sheep. In addition, we constructed co-expression and networks of lncRNAs-mRNAs. Specifically, XLOC_075176 targeted MYB, XLOC_014695 targeted VCAN, 106991527 targeted CASR, XLOC_075176 targeted KIFC1, XLOC_360232 targeted BRCA2. All these differential lncRNAs and mRNAs expression profiles in the thyroid provide a new resource for elucidating the regulatory mechanism underlying STH sheep prolificacy.

Circular RNA Profiling Identifies Novel circPPARA that Promotes Intramuscular Fat Deposition in Pigs

Intramuscular fat (IMF) content plays an important role in pork quality. Circular RNAs (circRNAs) implicate various biological processes; however, the regulatory mechanisms and functions of circRNAs in porcine IMF remains elusive. Hence, the study assessed the circRNA expression profiling in the longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to unravel their regulatory functions in improving meat quality. The RNA sequencing analysis identified 29,732 circRNAs from six sampled pigs, most of which were exon-derived. In the muscle, 336 were differentially expressed (DE) between the H and L IMF groups; 196 circRNAs were upregulated, and 140 were downregulated. Subsequent qRT-PCR validation of 10 DE circRNAs revealed expression patterns consistent with the RNA-seq data. Gene ontology and KEGG enrichment analysis revealed that most significantly enriched DE circRNAs' host genes were linked to lipid metabolism and adipogenesis processes. The circRNA-miRNA regulatory network analysis found several circRNAs targeting miRNAs associated with adipogenesis. Finally, a novel circRNA, circPPARA, was identified with the expression positively correlated with the IMF content. Detailed analysis revealed that circPPARA was formed via head-to-tail splicing and was more stable than the linear PPARA, predominantly located in the cytoplasm. Functional studies using overexpression and siRNA constructs demonstrated that circPPARA promotes differentiation and hinders the proliferation of porcine intramuscular preadipocytes. Moreover, the dual-luciferase assay revealed that circPPARA adsorbed miR-429 and miR-200b, thereby promoting intramuscular adipogenesis in pigs. Our results identified a candidate circRNA, circPPARA, that affects porcine IMF content. The study provides knowledge of the regulatory functions of circRNAs in intramuscular adipogenesis and abundant resource for future research on circRNAs in pigs.

MiR-218-5p Affects Subcutaneous Adipogenesis by Targeting ACSL1, A Novel Candidate for Pig Fat Deposition

As a centre enzyme in fatty acid activation, acyl-CoA synthetase long-chain family member 1 (ACSL1) plays an important role in body lipid homeostasis. However, the functions of ACSL1 in the subcutaneous adipogenesis of pigs are largely unknown. In the present study, we found that the expression of ACSL1 significantly increased during the process of porcine preadipocyte differentiation. Moreover, silencing of ACSL1 in preadipocytes decreased levels of triglyceride and adipogenic-related markers, including FABP4, APOE, and FASN (p < 0.01), and simultaneously increased levels of lipolytic-related markers, such as ATGL and HSL (p < 0.05). Conversely, overexpression of ACSL1 in preadipocytes increased levels of triglyceride and FABP4, APOE, and FASN (p < 0.01), and reduced levels of ATGL and HSL (p < 0.05). Luciferase reporter assays revealed that ACSL1 is a target of miR-218-5p, which can reduce the mRNA and protein levels of ACSL1 by directly binding the 3′ untranslated region of ACSL1. Furthermore, miR-218-5p has an inhibition role in porcine preadipocyte differentiation by suppressing ACSL1 expression. Taken together, these data provide insights into the mechanism of the miR-218-5p/ACSL1 axis in regulating subcutaneous fat deposition of pigs.

Genome-wide analysis of mRNAs, lncRNAs, and circRNAs during intramuscular adipogenesis in Chinese Guizhou Congjiang pigs

Intramuscular fat content is an important determinant of meat quality, and preadipocyte differentiation plays a critical role in intramuscular fat deposition in pigs. However, many types of RNA differentiation, including messenger RNA (mRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA) remain unreported despite their crucial roles in regulating adipogenesis. Chinese Guizhou Congjiang pigs are raised in the Guizhou province of China for their high-quality meat. Therefore, it is important for breeders to explore the mechanisms of proliferation and differentiation of intramuscular adipocytes from the longissimus dorsi muscle of these pigs. In the present study, a transcriptome analysis of intramuscular preadipocytes from Chinese Guizhou Congjiang pigs, including analyses of mRNAs, lncRNAs, and circRNAs at days 0 (D0), 4 (D4), and 8 (D8) was performed. A total of 1,538, 639, and 445 differentially expressed (DE) mRNAs, 479, 192, and 126 DE lncRNAs, and 360, 439, and 304 DE circRNAs were detected between D4 and D0, D8 and D0, and D8 and D4, respectively. Functional analyses identified many significantly enriched RNAs related to lipid deposition, cell differentiation, metabolism processes, and obesity-related diseases, biological processes, and pathways. We identified two lncRNAs (TCONS_00012086 and TCONS_00007245) closely related to fat deposition according to their target genes and tissue expression profiles. Subcellular distribution analysis using quantitative real-time PCR (qRT-PCR) revealed that both TCONS_00012086 and TCONS_00007245 are cytoplasmic lncRNAs. These data provide a genome-wide resource for mRNAs, lncRNAs, and circRNAs potentially involved in Chinese Guizhou Congjiang pig fat metabolism, thus improving our understanding of their function in adipogenesis.

Comprehensive Analysis of mRNA, lncRNA, circRNA, and miRNA Expression Profiles and Their ceRNA Networks in the Longissimus Dorsi Muscle of Cattle-Yak and Yak

Cattle-yak, as the hybrid offspring of cattle (Bos taurus) and yak (Bos grunniens), demonstrates obvious heterosis in production performance. Male hybrid sterility has been focused on for a long time; however, the mRNAs and non-coding RNAs related to muscle development as well as their regulatory networks remain unclear. The phenotypic data showed that the production performance (i.e., body weight, withers height, body length, and chest girth) of cattle-yak was significantly better than that of the yak, and the economic benefits of the cattle-yak were higher under the same feeding conditions. Then, we detected the expression profiles of the longissimus dorsi muscle of cattle-yak and yak to systematically reveal the molecular basis using the high-throughput sequencing technology. Here, 7,126 mRNAs, 791 lncRNAs, and 1,057 circRNAs were identified to be differentially expressed between cattle-yaks and yaks in the longissimus dorsi muscle. These mRNAs, lncRNA targeted genes, and circRNA host genes were significantly enriched in myoblast differentiation and some signaling pathways related to muscle development (such as HIF-1 signaling pathway and PI3K-Akt signaling pathway). We constructed a competing endogenous RNA (ceRNA) network and found that some non-coding RNAs differentially expressed may be involved in the regulation of muscle traits. Taken together, this study may be used as a reference tool to provide the molecular basis for studying muscle development.

Castration induced circRNA expressional changes in subcutaneous adipose tissue of male pigs

Circular RNAs (circRNAs) participated in regulation of lipid metabolism; however, its functional role on castration-induced lipid deposition has not been deeply researched. So in this research, we firstly compared circRNAs expressional differences in subcutaneous adipose tissue between intact and castrated male Huainan pigs. A total of 6116 differentially expressed circRNAs (DECs) were detected between these two groups (|log₂ foldchange| ≥ 1 and p_adj  ≤ 0.05); GO and KEGG analysis showed that their parent genes were mainly enriched in metabolism-related pathway. And TGF-beta, insulin, AMPK, and MAPK pathways might play vital role in castration-induced lipid deposition. The miRNAs enriched in the constructed circRNA-miRNA network were mainly participated in adipogenesis and lipid metabolism, such as miR-143a-3p, miR-378, and miR-195. And it was verified that testosterone upregulated miR-181a but downregulated circ_0005912 expression in a dose-dependent manner in porcine intramuscular adipocytes, and overexpression of miR-181a inhibited circ_0005912. Taken together, these DECs may participate in the regulation of lipid metabolism after castration by reaction with miRNAs, which indicated the novel role of circRNAs in castration-induced lipid deposition.

Identification of Circular RNA Expression Profiles in White Adipocytes and Their Roles in Adipogenesis

Obesity and its related metabolic diseases have become great public health threats worldwide. Although accumulated evidence suggests that circRNA is a new type of non-coding RNAs regulating various physiological and pathological processes, little attention has been paid to the expression profiles and functions of circRNAs in white adipose tissue. In this study, 3,771 circRNAs were detected in three stages of white adipogenesis (preadipocyte, differentiating preadipocyte, and mature adipocyte) by RNA-seq. Experimental validation suggested that the RNA-seq results are highly reliable. We found that nearly 10% of genes which expressed linear RNAs in adipocytes could also generate circRNAs. In addition, 40% of them produced multiple circRNA isoforms. We performed correlation analysis and found that a great deal of circRNAs (nearly 50%) and their parental genes were highly correlated in expression levels. A total of 41 differential expression circRNAs (DECs) were detected during adipogenesis and an extremely high ratio of them (80%) were correlated with their parental genes, indicating these circRNAs may potentially play roles in regulating the expression of their parental genes. KEGG enrichment and GO annotation of the parental genes suggesting that the DECs may participate in several adipogenesis-related pathways. Following rigorous selection, we found that many up-regulated circRNAs contain multiple miRNAs binding sites, such as miR17, miR-30c, and miR-130, indicating they may potentially facilitate their regulatory functions by acting as miRNA sponges. These results suggest that plenty of circRNAs are expressed in white adipogenesis and the DECs may serve as new candidates for future adipogenesis regulation.

Comprehensive analysis of lncRNAs involved in skeletal muscle development in ZBED6-knockout Bama pigs

BACKGROUND

The mutation of insulin-like growth factor 2 (IGF2 mutation) that a single-nucleotide substitution (G→A) in the third intron of IGF2 abrogates the interaction with zinc finger BED-type containing 6 (ZBED6) and leads to increased muscle mass in pigs. IGF2 mutation knock-in (IGF2 KI) and ZBED6 knockout (ZBED6 KO) lead to changes in IGF2 expression and increase muscle mass in mice and pigs. Long noncoding RNAs (lncRNAs) may participate in numerous biological processes, including skeletal muscle development. However, the role of the ZBED6-lncRNA axis in skeletal muscle development is poorly characterized.

RESULTS

In this study, we assembled transcriptomes using RNA-seq data published in previous studies by our group and identified 11,408 known lncRNAs and 2269 potential lncRNAs in seven tissues, heart, longissimus dorsi, gastrocnemius muscle, liver, spleen, lung and kidney, of ZBED6 KO (lean mass model) and WT Bama pigs. ZBED6 affected the expression of 1570 lncRNAs (differentially expressed lncRNAs [DE-lncRNAs]; log2-fold change ≥ 1, nominal p-value ≤ 0.05) in the seven examined tissues. The expressed lncRNAs (FPKM > 0.1) exhibited tissue-specific patterns in WT pigs. Specifically, 3410 lncRNAs were expressed exclusively in only one tissue. Potential functions of lncRNAs were indirectly predicted by searching their target cis- and trans-regulated protein-coding genes. LncRNAs with tissue-specific expression influence numerous genes related to tissue functions. Weighted gene coexpression network analysis (WGCNA) of 1570 DE-lncRNAs between WT and ZBED6 KO pigs was used to define the following six lncRNA modules specific to different tissues: skeletal muscle, heart, lung, spleen, kidney and liver modules. Furthermore, by conjoint analysis of longissimus dorsi data (tissue-specific expression, muscle module and DE-lncRNAs) and ChIP-PCR revealed NONSUSG002145.1 (adjusted p-values = 0.044), which is coexpressed with the IGF2 gene and binding with ZBED6, may play important roles in ZBED6 KO pig skeletal muscle development.

CONCLUSIONS

These findings indicate that the identified lncRNAs may play essential roles in tissue function and regulate the mechanism of ZBED6 action in skeletal muscle development in pigs. To our knowledge, this is the first study describing lncRNAs in ZBED6 KO pigs. These results may open new research directions leading to a better understanding of the global functions of ZBED6 and of lncRNA functions in skeletal muscle development in pigs.

Comparison of long non-coding RNAs in adipose and muscle tissues between seven indigenous Chinese and the Yorkshire pig breeds

lncRNAs play crucial roles in fat metabolism in animals. Previously, we have compared the mRNA transcriptome profiles between seven fat-type Chinese pig breeds and one lean-type Western breed (Yorkshire, YY). The associations between differentially expressed (DE) genes and phenotypical traits were investigated. In the present study, to further explore the underlying regulatory mechanisms, lncRNAs were sequenced and compared between YY and Chinese indigenous breeds. The results showed 9114 and 7538 DE lncRNAs between at least one Chinese breed and the YY breed in the adipose and muscle tissue respectively. KEGG enrichment analysis revealed that the target genes of these DE lncRNAs mainly influenced the glucolipid metabolism, which is an important process affecting meat quality. Correlation analyses between the DE lncRNA and DE mRNA genes related to meat quality and growth traits were performed. The results showed that LTCONS_00073280 was associated with intramuscular fat content. Four lncRNAs (LTCONS_00101781, LTCONS_00037879, LTCONS_00088260 and LTCONS-00128343) might mediate backfat thickness. Overall, this study provides candidate lncRNAs that potentially affect meat quality, which might be useful for molecular breeding of pig breeds in future.

Dynamic Expression Profiles of Circular RNAs during Brown to White Adipose Tissue Transformation in Goats (Capra hircus)

Simple Summary

In our study, we launched RNA-seq in order to investigate the potential functions of circRNA during brown adipose tissue (BAT) to white adipose tissue (WAT) transformation. As a result, 6610 circRNAs and 61 differentially expressed circRNAs (DEcircRNAs) were identified. Moreover, 65 miRNAs were detected that could potentially interact with DEcircRNAs. The present study provides a detailed circRNA expression landscape and evidence for circRNA functions in the transformation from BAT to WAT.

Abstract

Adipose tissues are mainly divided into brown adipose tissue (BAT) and white adipose tissue (WAT). WAT mainly functions to buffer excess calories, whereas BAT plays a role in the non-shivering thermogenesis to maintain body temperature and energy balance. Moreover, circRNAs play important roles in various biological processes. However, knowledge of the expression profile and function of circRNAs from BAT to WAT remains largely unknown. In this study, a total of 6610 unique circRNAs were identified in the perirenal adipose tissues of 1-day, 30-days, and 1-year goats. Functional annotation revealed that host genes of circRNAs were involved in some BAT-related pathways, such as the thyroid hormone signaling pathway, MAPK signaling pathway, and VEGF signaling pathway. Furthermore, a total of 61 DEcircRNAs were detected across three stages. Additionally, five selected circRNAs were validated by RNase R assay, qPCR, and Sanger sequencing. Finally, the circRNA–miRNA network was constructed between the DEcircRNAs and their miRNA binding sites.

Characteristics and expression profiles of circRNAs during abdominal adipose tissue development in Chinese Gushi chickens

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.

Pineal gland transcriptomic profiling reveals the differential regulation of lncRNA and mRNA related to prolificacy in STH sheep with two FecB genotypes

BACKGROUND

Long noncoding RNA (lncRNA) has been identified as important regulator in hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, little is known of their expression pattern and potential roles in the pineal gland of sheep. Herein, RNA-Seq was used to detect transcriptome expression pattern in pineal gland between follicular phase (FP) and luteal phase (LP) in FecB^BB (MM) and FecB⁺⁺ (ww) STH sheep, respectively, and differentially expressed (DE) lncRNAs and mRNAs associated with reproduction were identified.

RESULTS

Overall, 135 DE lncRNAs and 1360 DE mRNAs in pineal gland between MM and ww sheep were screened. Wherein, 39 DE lncRNAs and 764 DE mRNAs were identified (FP vs LP) in MM sheep, 96 DE lncRNAs and 596 DE mRNAs were identified (FP vs LP) in ww sheep. Moreover, GO and KEGG enrichment analysis indicated that the targets of DE lncRNAs and DE mRNAs were annotated to multiple biological processes such as phototransduction, circadian rhythm, melanogenesis, GSH metabolism and steroid biosynthesis, which directly or indirectly participate in hormone activities to affect sheep reproductive performance. Additionally, co-expression of lncRNAs-mRNAs and the network construction were performed based on correlation analysis, DE lncRNAs can modulate target genes involved in related pathways to affect sheep fecundity. Specifically, XLOC_466330, XLOC_532771, XLOC_028449 targeting RRM2B and GSTK1, XLOC_391199 targeting STMN1, XLOC_503926 targeting RAG2, XLOC_187711 targeting DLG4 were included.

CONCLUSION

All of these differential lncRNAs and mRNAs expression profiles in pineal gland provide a novel resource for elucidating regulatory mechanism underlying STH sheep prolificacy.

Comparative analysis of gene expression profiles in differentiated subcutaneous adipocytes between Jiaxing Black and Large White pigs

Background

Chinese domestic pig breeds are reputed for pork quality, but their low ratio of lean-to-fat carcass weight decreases production efficiency. A better understanding of the genetic regulation network of subcutaneous fat tissue is necessary for the rational selection of Chinese domestic pig breeds. In the present study, subcutaneous adipocytes were isolated from Jiaxing Black pigs a Chinese indigenous pig breed with redundant subcutaneous fat deposition and Large White pigs a lean-type pig breed with relatively low subcutaneous fat deposition. The expression profiles of mRNAs and lncRNAs were compared by RNA-seq analysis to identify biomarkers correlated with the differences of subcutaneous fat deposition between the two breeds.

Results

A total of 1058 differentially expressed genes and 221 differentially expressed lncRNAs were identified in subcutaneous adipocytes between Jiaxing Black and Large White pigs, which included 275 up-regulated mRNAs, 783 down-regulated mRNAs, 118 up-regulated lncRNAs and 103 down-regulated lncRNAs. Gene Ontology and KEGG pathway enrichment analyses revealed that the differentially expressed genes and differentially expressed lncRNAs were mainly involved in the immune response, cell fate determination, PI3K-Akt signaling pathway and MAPK signaling pathway, which are known to be related to adipogenesis and lipid metabolism. The expression levels of differentially expressed genes and differentially expressed lncRNAs according to the RNA-seq data were verified by quantitative PCR, which showed 81.8% consistency. The differences in MAPK pathway activity between Jiaxing Black and Large White pigs was confirmed by western blot analysis, which revealed elevated p38 phosphorylation in Jiaxing Black pigs.

Conclusions

This study offers a detailed characterization of mRNAs and lncRNAs in fat- and lean-type pig breeds. The activity of the MAPK signaling pathway was found to be associated with subcutaneous adipogenesis. These results provide new targets for further investigation of the molecular mechanisms regulating subcutaneous fat deposition in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07361-9.

The expression profile and role of non-coding RNAs in obesity

Latest years have experienced a dramatic upsurge in the knowledge about the function of non-coding transcripts in the determination of diverse human phenotypes including obesity. Several miRNAs and lncRNAs participate in the regulation of metabolic pathways leading to obesity. Several lncRNAs such as Mist, lincIRS2, lncRNA-p5549, H19, GAS5 and SNHG9 have been shown to be down-regulated in adipose tissues or other biological samples in the obese human or animal subjects. On the other hand, Meg3, Plnc1, Blnc1, AC092834.1, TINCR and PVT1 are among up-regulated lncRNAs in the obese subjects. Tens of miRNAs have differential expression between obese and non-obese subjects or between mature adipocytes and pre-adipocytes. Understanding the molecular mechanism of involvement of non-coding RNAs in the pathobiology of obesity would simplify design of therapeutic choices for protecting against obesity and its related comorbidities. We explain the available literature on the function of these transcripts in the pathobiology of obesity.

Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.

Myostatin regulates fatty acid desaturation and fat deposition through MEF2C/miR222/SCD5 cascade in pigs

Myostatin (MSTN), associated with the “double muscling” phenotype, affects muscle growth and fat deposition in animals, whereas how MSTN affects adipogenesis remains to be discovered. Here we show that MSTN can act through the MEF2C/miR222/SCD5 cascade to regulate fatty acid metabolism. We generated MSTN-knockout (KO) cloned Meishan pigs, which exhibits typical double muscling trait. We then sequenced transcriptome of subcutaneous fat tissues of wild-type (WT) and MSTN-KO pigs, and intersected the differentially expressed mRNAs and miRNAs to predict that stearoyl-CoA desaturase 5 (SCD5) is targeted by miR222. Transcription factor binding prediction showed that myogenic transcription factor 2C (MEF2C) potentially binds to the miR222 promoter. We hypothesized that MSTN-KO upregulates MEF2C and consequently increases the miR222 expression, which in turn targets SCD5 to suppress its translation. Biochemical, molecular and cellular experiments verified the existence of the cascade. This novel molecular pathway sheds light on new targets for genetic improvements in pigs.

Ren, Xiao et al. identify a mechanism by which myostatin regulates adipogenesis, using myostatin-knockout pigs. Myostatin deficiency upregulates MEF2C that binds to the promoter of miR222. miR222 in turn downregulates stearoyl-CoA desaturase 5. This study provides potential targets that can be engineered to generate a new pig variety that has high leanness while maintaining its high intramuscular fat content.

Regulatory roles of circRNAs in adipogenesis and lipid metabolism: emerging insights into lipid-related diseases

Disorder of lipid metabolism has become an urgent health problem that brings about a variety of metabolic syndromes, including hepatic steatosis, adipose tissue dysfunction, diabetes and obesity. Circular RNAs (circRNAs), a class of emerging RNA molecules with unique structure and extensive effects, have been verified to participate in various biological programs through distinct mechanisms, especially in lipid-related processes. In this review, the biogenesis, characteristics, and functional mechanisms of circRNAs are discussed. Furthermore, the methods for circRNA identification and expression profiles of circRNAs associated with adipogenesis and lipid metabolism are described. Additionally, we emphasize the regulatory roles of circRNAs in adipogenesis, lipid metabolism, and lipid-related diseases. Finally, the diagnostic and therapeutic potential of circRNAs is highlighted, showing potential for the clinical application of circRNAs in the treatment of lipid-related diseases in the near future.

Exploration of Long Non-coding RNAs and Circular RNAs in Porcine Milk Exosomes

RNA in milk exosomes can be absorbed in the mammalian intestinal tract and function in gene expression regulations. Our previous work demonstrated that porcine milk exosomes (PME) contain large amounts of miRNAs and mRNAs. Increasing evidence suggests that long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are of particular interest, given their key role in diverse biological processes of animals. However, the expression profiles and the potential functions of lncRNAs and circRNAs in PME are still unknown. In the present study, we isolated PME by ultracentrifugation and performed a comprehensive analysis of lncRNA and circRNA in PME by using RNA sequencing. As a result, 2,466 novel lncRNAs, 809 annotated lncRNAs, and 61 circRNAs were identified in PME. The lncRNAs shared similar characteristics with other mammals in terms of length, exon number, and open reading frames. However, lncRNAs showed a higher level compared with mRNAs. Eight lncRNAs and five circRNAs in PME were selected for PCR identification. A functional enrichment analysis on the target genes of lncRNAs indicated that these genes were involved in cellular macromolecule metabolic, RNA metabolic, and immune processes. The circRNAs host genes were enriched in nucleic acid binding and adherence junction. We also evaluated the potential interaction targets between miRNAs and PME lncRNAs or circRNAs, and the results showed that the PME lncRNAs and the circRNAs have a high density of miRNA target sites. The top 20 highly expressed lncRNAs were found to interact with the proliferation-related miRNAs, and the circRNAs potentially targeted many miRNAs that are associated with the intestinal barrier. This study is the first to provide a resource for lncRNA and circRNA research of porcine milk. Moreover, the potential interaction between lncRNA/circRNA and miRNA is revealed. The present study expands our knowledge of non-coding RNAs in milk, and additional research is necessary to confirm their exactly physiological functions.

The Emerging Role and Promise of Circular RNAs in Obesity and Related Metabolic Disorders

Circular RNAs (circRNAs) are genome transcripts that are produced from back-splicing of specific regions of pre-mRNA. These single-stranded RNA molecules are widely expressed across diverse phyla and many of them are stable and evolutionary conserved between species. Growing evidence suggests that many circRNAs function as master regulators of gene expression by influencing both transcription and translation processes. Mechanistically, circRNAs are predicted to act as endogenous microRNA (miRNA) sponges, interact with functional RNA-binding proteins (RBPs), and associate with elements of the transcriptional machinery in the nucleus. Evidence is mounting that dysregulation of circRNAs is closely related to the occurrence of a range of diseases including cancer and metabolic diseases. Indeed, there are several reports implicating circRNAs in cardiovascular diseases (CVD), diabetes, hypertension, and atherosclerosis. However, there is very little research addressing the potential role of these RNA transcripts in the occurrence and development of obesity. Emerging data from in vitro and in vivo studies suggest that circRNAs are novel players in adipogenesis, white adipose browning, obesity, obesity-induced inflammation, and insulin resistance. This study explores the current state of knowledge on circRNAs regulating molecular processes associated with adipogenesis and obesity, highlights some of the challenges encountered while studying circRNAs and suggests some perspectives for future research directions in this exciting field of study.

Comprehensive Analysis of Differentially Expressed Profiles of mRNA, lncRNA, and circRNA in the Uterus of Seasonal Reproduction Sheep

Photoperiod is one of the important factors leading to seasonal reproduction of sheep. However, the molecular mechanisms underlying the photoperiod regulation of seasonal reproduction remain poorly understood. In this study, we compared the expression profiles of mRNAs, lncRNAs, and circRNAs in uterine tissues from Sunite sheep during three different photoperiods, namely, the short photoperiod (SP), short transfer to long photoperiod (SLP), and long photoperiod (LP). The results showed that 298, 403, and 378 differentially expressed (DE) mRNAs, 171, 491, and 499 DE lncRNAs, and 124, 270, and 400 DE circRNAs were identified between SP and LP, between SP and SLP, and between LP and SLP, respectively. Furthermore, functional enrichment analysis showed that the differentially expressed RNAs were mainly involved in the GnRH signaling pathway, thyroid hormone synthesis, and thyroid hormone signaling pathway. In addition, co-expression networks of lncRNA–mRNA were constructed based on the correlation analysis between the differentially expressed RNAs. Our study provides new insights into the expression changes of RNAs in different photoperiods, which might contribute to understanding the molecular mechanisms of seasonal reproduction in sheep.

Transcriptome Analysis of Landrace Pig Subcutaneous Preadipocytes during Adipogenic Differentiation

Fat deposition in pigs, which significantly contributes to meat quality, fattening efficiency, reproductive performance, and immunity, is critically affected by preadipocyte adipogenic differentiation. We elucidated adipogenesis in pigs using transcriptome analysis. Preadipocytes from subcutaneous adipose tissue (SAT) of Landrace piglets were differentiated into adipocytes in vitro. RNA sequencing (RNA-seq) used to screen differentially expressed genes (DEGs) during preadipocyte differentiation up to day 8 revealed 15,918 known and 586 novel genes. We detected 21, 144, and 394 DEGs, respectively, including 16 genes differentially expressed at days 2, 4 and 8 compared to day 0. Th number of DEGs increased time-dependently. Lipid metabolism, cell differentiation and proliferation, peroxisome proliferator-activated receptor (PPAR), wingless-type MMTV integration site (Wnt), tumor necrosis factor (TNF) signaling, and steroid biosynthesis were significant at days 2, 4, and 8 compared to day 0 (adjusted p < 0.05). Short time-series expression miner (STEM) analysis obtained 26 clusters of differential gene expression patterns, and nine were significant (p < 0.05). Functional analysis showed many significantly enriched lipid deposition- and cellular process-related biological processes and pathways in profiles 9, 21, 22, and 24. Glycerolipid and fatty-acid metabolism, PPAR signaling, fatty-acid degradation, phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt), and TNF signaling were observed during preadipocyte differentiation in vitro. These findings will facilitate the comprehension of preadipocyte differentiation and fat deposition in pigs.

GROWTH AND DEVELOPMENT SYMPOSIUM: STEM AND PROGENITOR CELLS IN ANIMAL GROWTH: Long noncoding RNAs in adipogenesis and adipose development of meat animals12

Sequencing technology, especially next-generation RNA sequencing, has greatly facilitated the identification and annotation of long noncoding RNAs (lncRNAs). In mammals, a large number of lncRNAs have been identified, which regulate various biological processes. An increasing number of lncRNAs have been identified which could function as key regulators of adipogenesis (adipocyte formation), a key step of the development of adipose tissue. Because proper adipose tissue development is a key factor affecting animal growth efficiency, lean/fat ratio, and meat quality, summarizing the roles and recent advances of lncRNAs in adipogenesis is needed in order to develop strategies to effectively manage fat deposition. In this review, we updated lncRNAs contributed to the regulation of adipogenesis, focusing on their roles in fat development of farm animals.

Comprehensive Analysis of Differentially Expressed mRNA, lncRNA and circRNA and Their ceRNA Networks in the Longissimus Dorsi Muscle of Two Different Pig Breeds

Circular RNA (circRNA) and long non-coding RNA (lncRNA) are known to participate in adipogenesis and myogenic differentiation, but their impact on porcine muscle traits is not well understood. We compared their expressional profiles in the longissimus dorsi muscle of Chinese Huainan pigs (HN, the fat type) and Western commercial Duroc×(Landrace×Yorkshire) (DLY, the thin type) pigs, and 854 mRNAs, 233 lncRNAs, and 66 circRNAs (p < 0.05 and ｜log₂FoldChange｜>1) were found to be differentially expressed. The differentially expressed mRNA and circRNA parental genes were enriched in the Wnt signaling pathway (adipogenesis), the transition between fast and slow fibers (myogenic differentiation), and alanine, aspartate and glutamate metabolism (pork flavor). The potential lncRNAs/circRNAs-miRNAs-mRNAs regulatory networks shared MYOD1, PPARD, miR-423-5p and miR-874, which were associated with skeletal muscle muscular proliferation, differentiation/regeneration and adipogenesis. Taken together, these differentially expressed non-coding RNAs may be involved in the molecular basis of muscle traits, acting as the competing endogenous RNA (ceRNA) for miRNAs.