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

Transcriptome and Metabolome Insights into Key Genes Regulating Fat Deposition and Meat Quality in Pig Breeds

Meat quality is a complex trait that exhibits significant variation across pig breeds, and the regulatory mechanisms governing pork meat quality are not fully elucidated. We compared the transcriptomics and metabolomics of the longissimus dorsi (LD) muscle between the Songliao Black Pig (SBP) and Large White × Landrace Pig (LWLDP) to investigate breed-specific differences in meat quality and underlying regulatory pathways. The results showed that SBP meat had a higher marbling score and backfat thickness, a richer color, a lower shear force, and reduced drip loss. Fatty acid (FA) analysis identified 15 significant FAs in the LWLDP, with docosahexaenoic acid (DHA) in the SBP, while amino acid (AA) analysis revealed no breed-based differences. Transcriptome analysis identified 134 upregulated and 362 downregulated genes in the SBP. Protein-protein interaction (PPI) network analysis found 25 key genes, which are associated with muscle development, fat deposition, and overall meat quality, while genes in the insulin signaling pathway, such as PPP1R3B, PPARGC1A, SOCS1, EIF4E, PRKAR2A, PRKAG2, and FASN, play a crucial role in balancing fat metabolism and catabolism. Metabolomic analysis identified 89 upregulated and 10 downregulated metabolites in the SBP, primarily involved in fructose and mannose metabolism, amino acid biosynthesis, nucleotide sugar metabolism, and glucagon signaling pathways. Gene-metabolite association analysis found that the PPP1R3B gene had a strong association with Thr-Leu, Maltol, D-myo-Inositol-4-phosphate, and Fructose-6-phosphate, while MYOG correlated with Mannose-6-phosphate, Fructose-1-phosphate, Mannose-1-phosphate, and Glucose-6-phosphate. In contrast, NR4A3 and PPARGC1A showed a strong negative correlation with most upregulated metabolites. In conclusion, this study identified functional genes, elucidated the mechanisms associated with meat quality traits, and identified gene-metabolite associations involved in energy metabolism, muscle development, and fat deposition, providing valuable insights into the molecular mechanisms that regulate meat quality between pig breeds.

The long noncoding RNA lncMPD2 inhibits myogenesis by targeting the miR-34a-5p/THBS1 axis

Long noncoding RNAs (lncRNAs) participate in regulating skeletal muscle development. However, little is known about their role in regulating chicken myogenesis. In this study, we identified a novel lncRNA, lncMPD2, through transcriptome sequencing of chicken myoblasts at different developmental stages. Functionally, gain- and loss-of-function experiments showed that lncMPD2 inhibited myoblast proliferation and differentiation. Mechanistically, lncMPD2 directly bound to miR-34a-5p, and miR-34a-5p promoted myoblasts proliferation and differentiation and inhibited the mRNA and protein expression of its target gene THBS1. THBS1 inhibited myoblast proliferation and differentiation in vitro and delayed muscle regeneration in vivo. Furthermore, rescue experiments showed that lncMPD2 counteracted the inhibitory effects of miR-34a-5p on THBS1 and myogenesis-related gene mRNA and protein expression. In conclusion, lncMPD2 regulates the miR-34a-5p/THBS1 axis to inhibit the proliferation and differentiation of myoblasts and skeletal muscle regeneration. This study provides more insight into the molecular regulatory network of skeletal muscle development, identifying novel potential biomarkers for improving chicken quality and increasing chicken yield. In addition, this study provides a potential goal for breeding strategies that minimize muscle damage in chickens.

Research Progress on the Regulating Factors of Muscle Fiber Heterogeneity in Livestock: A Review

The type of muscle fiber plays a crucial role in the growth, development, and dynamic plasticity of animals' skeletal muscle. Additionally, it is a primary determinant of the quality of both fresh and processed meat. Therefore, understanding the regulatory factors that contribute to muscle fibers' heterogeneity is of paramount importance. Recent advances in sequencing and omics technologies have enabled comprehensive cross-verification of research on the factors affecting the types of muscle fiber across multiple levels, including the genome, transcriptome, proteome, and metabolome. These advancements have facilitated deeper exploration into the related biological questions. This review focused on the impact of individual characteristics, feeding patterns, and genetic regulation on the proportion and interconversion of different muscle fibers. The findings indicated that individual characteristics and feeding patterns significantly influence the type of muscle fiber, which can effectively enhance the type and distribution of muscle fibers in livestock. Furthermore, non-coding RNA, genes and signaling pathways between complicated regulatory mechanisms and interactions have a certain degree of impact on muscle fibers' heterogeneity. This, in turn, changes muscle fiber profile in living animals through genetic selection or environmental factors, and has the potential to modulate the quality of fresh meat. Collectively, we briefly reviewed the structure of skeletal muscle tissue and then attempted to review the inevitable connection between the quality of fresh meat and the type of muscle fiber, with particular attention to potential events involved in regulating muscle fibers' heterogeneity.

Long Non-Coding RNAs Differentially Expressed in Swine Fetuses

Long non-coding RNAs (lncRNAs) are non-coding transcripts involved in various biological processes. The Y chromosome is known for determining the male sex in mammals. LncRNAs on the Y chromosome may play important regulatory roles. However, knowledge about their action mechanisms is still limited, especially during early fetal development. Therefore, we conducted this exploratory study aiming to identify, characterize, and investigate the differential expression of lncRNAs between male and female swine fetuses at 35 days of gestation. RNA-Seq libraries from 10 fetuses were prepared and sequenced using the Illumina platform. After sequencing, a data quality control was performed using Trimmomatic, alignment with HISAT2, and transcript assembly with StringTie. The differentially expressed lncRNAs were identified using the limma package of the R software (4.3.1). A total of 871 potentially novel lncRNAs were identified and characterized. Considering differential expression, eight lncRNAs were upregulated in male fetuses. One was mapped onto SSC12 and seven were located on the Y chromosome; among them, one lncRNA is potentially novel. These lncRNAs are involved in diverse functions, including the regulation of gene expression and the modulation of chromosomal structure. These discoveries enable future studies on lncRNAs in the fetal stage in pigs.

Whole-Transcriptome Analysis Sheds Light on the Biological Contexts of Intramuscular Fat Deposition in Ningxiang Pigs

The quality of pork is significantly impacted by intramuscular fat (IMF). However, the regulatory mechanism of IMF depositions remains unclear. We performed whole-transcriptome sequencing of the longissimus dorsi muscle (IMF) from the high (5.1 ± 0.08) and low (2.9 ± 0.51) IMF groups (%) to elucidate potential mechanisms. In summary, 285 differentially expressed genes (DEGs), 14 differentially expressed miRNAs (DEMIs), 83 differentially expressed lncRNAs (DELs), and 79 differentially expressed circRNAs (DECs) were identified. DEGs were widely associated with IMF deposition and liposome differentiation. Furthermore, competing endogenous RNA (ceRNA) regulatory networks were constructed through co-differential expression analyses, which included circRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 47 DECs) and lncRNA-miRNA-mRNA (containing 6 DEMIs, 6 DEGs, 36 DELs) regulatory networks. The circRNAs sus-TRPM7_0005, sus-MTUS1_0004, the lncRNAs SMSTRG.4269.1, and MSTRG.7983.2 regulate the expression of six lipid metabolism-related target genes, including PLCB1, BAD, and GADD45G, through the binding sites of 2-4068, miR-7134-3p, and miR-190a. For instance, MSTRG.4269.1 regulates its targets PLCB1 and BAD via miRNA 2_4068. Meanwhile, sus-TRPM7_0005 controls its target LRP5 through ssc-miR-7134-3P. These findings indicate molecular regulatory networks that could potentially be applied for the marker-assisted selection of IMF to enhance pork quality.

Whole transcriptome profiling reveals a lncMDP1 that regulates myogenesis by adsorbing miR-301a-5p targeting CHAC1

Myoblast proliferation and differentiation are essential for skeletal muscle development. In this study, we generated the expression profiles of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in different developmental stages of chicken primary myoblasts (CPMs) using RNA sequencing (RNA-seq) technology. The dual luciferase reporter system was performed using chicken embryonic fibroblast cells (DF-1), and functional studies quantitative real-time polymerase chain reaction (qPCR), cell counting kit-8 (CCK-8), 5-Ethynyl-2'-deoxyuridine (EdU), flow cytometry cycle, RNA fluorescence in situ hybridization (RNA-FISH), immunofluorescence, and western blotting assay. Our research demonstrated that miR-301a-5p had a targeted binding ability to lncMDP1 and ChaC glutathione-specific gamma-glutamylcyclotransferase 1 (CHAC1). The results revealed that lncMDP1 regulated the proliferation and differentiation of myoblasts via regulating the miR-301a-5p/CHAC1 axis, and CHAC1 promotes muscle regeneration. This study fulfilled the molecular regulatory network of skeletal muscle development and providing an important theoretical reference for the future improvement of chicken meat performance and meat quality.

Comparative transcriptome analysis of longissimus dorsi muscle reveal potential genes affecting meat trait in Chinese indigenous Xiang pig

In this study, we compared the transcriptome of longissimus dorsi muscle between Guizhou Xiang pigs (XP) and Western commercial Large White pigs (LW), which show diffirent meat quality between them. In terms of meat quality traits, the pH 45 min, color score, backfat thickness, and intramuscular fat (IMF) content were higher in Xiang pigs than in Large White pigs (P < 0.01), while the drip loss, lean meat percentage, shear force, and longissimus dorsi muscle area of Xiang pigs were lower than that of Large White pigs (P < 0.01). Nutrients such as monounsaturated fatty acid (MUFA), total amino acids (TAA), delicious amino acids (DAA) and essential amino acids (EAA) in Xiang pigs were higher than that in Large White pigs, and the proportion of polyunsaturated fatty acid (PUFA) of Xiang pigs was significantly lower than Large White pigs (P < 0.01). Transcriptome analysis identified 163 up-regulated genes and 88 genes down-regulated in Xiang pigs longissimus dorsi muscle. Combined with the correlation analysis and quantitative trait locis (QTLs) affecting meat quality, a total of 227 DEGs were screened to be significantly associated with meat quality values. Enrichment analysis indicated that numerous members of genes were gathered in muscle development, adipogenesis, amino acid metabolism, fatty acid metabolism and synthesis. Of those, 29 genes were identified to be hub genes that might be related with the meat quality of Xiang pig, such as MYOD1, ACTB, ASNS, FOXO1, ARG2, SLC2A4, PLIN2, and SCD. Thus, we screened and identified the potential functional genes for the formation of meat quality in Xiang pigs, which provides a corresponding theoretical basis for the study of the molecular regulatory mechanism of pork quality and the improvement of pork quality.

Transcriptome Data Revealed the circRNA-miRNA-mRNA Regulatory Network during the Proliferation and Differentiation of Myoblasts in Shitou Goose

CircRNA, a recently characterized non-coding RNA (ncRNA) variant, functions as a molecular sponge, exerting regulatory control by binding to microRNA (miRNA) and modulating the expression of downstream proteins, either promoting or inhibiting their expression. Among poultry species, geese hold significant importance, prized by consumers for their delectable taste and rich nutritional content. Despite the prominence of geese, research on the growth and development of goose muscle, particularly the regulatory role of circRNAs in goose muscle formation, remains insufficiently explored. In this study, we constructed comprehensive expression profiles of circRNAs and messenger RNAs (mRNAs) within the myoblasts and myotubes of Shitou geese. We identified a total of 96 differentially expressed circRNAs (DEcircRNAs) and 880 differentially expressed mRNAs (DEmRNAs). Notably, the parental genes of DEcircRNAs and DEmRNAs exhibited enrichment in the Wnt signaling pathway, highlighting its potential impact on the proliferation and differentiation of goose myoblasts. Employing RNAhybrid and miRDB, we identified circRNA-miRNA pairs and mRNA-miRNA pairs that may play a role in regulating myogenic differentiation or muscle growth. Subsequently, utilizing Cytoscape, we constructed a circRNA-miRNA-mRNA interaction network aimed at unraveling the intricate regulatory mechanisms involved in goose muscle growth and development, which comprises 93 circRNAs, 351 miRNAs, and 305 mRNAs. Moreover, the identification of 10 hub genes (ACTB, ACTN1, BDNF, PDGFRA, MYL1, EFNA5, MYSM1, THBS1, ITGA8, and ELN) potentially linked to myogenesis, along with the exploration of their circRNA-miRNA-hub gene regulatory axis, was also conducted. These competitive endogenous RNA (ceRNA) regulatory networks elucidate the molecular regulatory mechanisms associated with muscle growth in Shitou geese, providing deeper insights into the reciprocal regulation of circRNA, miRNA, and mRNA in the context of goose muscle formation.

Construction of LncRNA-Related ceRNA Networks in Longissimus Dorsi Muscle of Jinfen White Pigs at Different Developmental Stages

The development of skeletal muscle in pigs might determine the quality of pork. In recent years, long non-coding RNAs (lncRNAs) have been found to play an important role in skeletal muscle growth and development. In this study, we investigated the whole transcriptome of the longissimus dorsi muscle (LDM) of Jinfen White pigs at three developmental stages (1, 90, and 180 days) and performed a comprehensive analysis of lncRNAs, mRNAs, and micro-RNAs (miRNAs), aiming to find the key regulators and interaction networks in Jinfen White pigs. A total of 2638 differentially expressed mRNAs (DE mRNAs) and 982 differentially expressed lncRNAs (DE lncRNAs) were identified. Compared with JFW_1d, there were 497 up-regulated and 698 down-regulated DE mRNAs and 212 up-regulated and 286 down-regulated DE lncRNAs in JFW_90d, respectively. In JFW_180d, there were 613 up-regulated and 895 down-regulated DE mRNAs and 184 up-regulated and 131 down-regulated DE lncRNAs compared with JFW_1d. There were 615 up-regulated and 477 down-regulated DE mRNAs and 254 up-regulated and 355 down-regulated DE lncRNAs in JFW_180d compared with JFW_90d. Compared with mRNA, lncRNA has fewer exons, fewer ORFs, and a shorter length. We performed GO and KEGG pathway functional enrichment analysis for DE mRNAs and the potential target genes of DE lncRNAs. As a result, several pathways are involved in muscle growth and development, such as the PI3K-Akt, MAPK, hedgehog, and hippo signaling pathways. These are among the pathways through which mRNA and lncRNAs function. As part of this study, bioinformatic screening was used to identify miRNAs and DE lncRNAs that could act as ceRNAs. Finally, we constructed an lncRNA-miRNA-mRNA regulation network containing 26 mRNAs, 7 miRNAs, and 17 lncRNAs; qRT-PCR was used to verify the key genes in these networks. Among these, XLOC_022984/miR-127/ENAH and XLOC_016847/miR-486/NRF1 may function as key ceRNA networks. In this study, we obtained transcriptomic profiles from the LDM of Jinfen White pigs at three developmental stages and screened out lncRNA-miRNA-mRNA regulatory networks that may provide crucial information for the further exploration of the molecular mechanisms during skeletal muscle development.

Delivery of microRNA-423-5p by exosome from adipose-derived stem/stromal cells inhibits DVL3 to potentiate autologous fat graft survival through adipogenesis and inflammatory response

Autologous fat grafting represents a reconstructive technique but is limited by unstable graft retention. Based on existing reports and bioinformatics prediction, we hypothesized that delivering exosomes from human adipose-derived stem/stromal cells (hADSC-Exo) would increase fat graft survival and further explore the mechanism. hADSC-Exo were extracted and identified. An autologous fat grafting model was established using donor and recipient mice, followed by hADSC-Exo treatment. hADSC-Exo promoted the retention of autologous fat grafts in mice, along with increased adipocyte activity, angiogenesis, and decreased inflammation in grafts. Moreover, hADSC-Exo potentiated the adipose differentiation of 3T3-L1 cells, enhanced the angiogenic and migratory capacity of human umbilical vein endothelial cells, and inhibited the inflammation and viability of RAW 264.7 cells. The therapeutic effect of hADSC-Exo on fat grafting was associated with the delivery of microRNA (miR)-423-5p. Deletion of miR-423-5p in Exo impaired the function of hADSC-Exo on fat retention. miR-423-5p bound to DVL3 to suppress DVL3 expression, and DVL3 deletion promoted adipose differentiation of 3T3-L1 cells. In conclusion, our findings further widen the theoretical basis of the clinical application of hADSC-Exo in autologous fat grafts.

LncRNA and Protein Expression Profiles Reveal Heart Adaptation to High-Altitude Hypoxia in Tibetan Sheep

The Tibetan sheep has an intricate mechanism of adaptation to low oxygen levels, which is influenced by both genetic and environmental factors. The heart plays a crucial role in the adaptation of Tibetan sheep to hypoxia. In the present study, we utilized transcriptomic and proteomic technologies to comprehensively analyze and identify the long non-coding RNAs (lncRNAs), genes, proteins, pathways, and gene ontology (GO) terms associated with hypoxic adaptation in Tibetan sheep at three different altitudes (2500 m, 3500 m, and 4500 m). By integrating the differentially expressed (DE) lncRNA target genes, differentially expressed proteins (DEPs), and differentially expressed genes (DEGs), we were able to identify and characterize the mechanisms underlying hypoxic adaptation in Tibetan sheep. Through this integration, we identified 41 shared genes/proteins, and functional enrichment analyses revealed their close association with lipid metabolism, glycolysis/gluconeogenesis, and angiogenesis. Additionally, significant enrichment was observed in important pathways such as the PPAR signaling pathway, glycolysis/gluconeogenesis, the oxoacid metabolic process, and angiogenesis. Furthermore, the co-expression network of lncRNAs and mRNAs demonstrated that lncRNAs (MSTRG.4748.1, ENSOART00020025894, and ENSOART00020036371) may play a pivotal role in the adaptation of Tibetan sheep to the hypoxic conditions of the plateau. In conclusion, this study expands the existing database of lncRNAs and proteins in Tibetan sheep, and these findings may serve as a reference for the prevention of altitude sickness in humans.

Whole-Transcriptome RNA Sequencing Uncovers the Global Expression Changes and RNA Regulatory Networks in Duck Embryonic Myogenesis

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

Functional analysis of differentially expressed circular RNAs in sheep subcutaneous fat

BACKGROUND

Circular RNAs (circRNAs), as important non-coding RNAs (ncRNAs), are involved in many biological activities. However, the exact chemical mechanism behind fat accumulation is unknown. In this paper, we obtained the expression profiles of circRNAs using high-throughput sequencing and investigated their differential expression in subcutaneous fat tissue of Duolang and Small Tail Han sheep.

RESULTS

From the transcriptomic analysis, 141 differentially expressed circRNAs were identified, comprising 61 up-regulated circRNAs and 80 down-regulated circRNAs. These host genes were primarily enriched in the MAPK and AMPK signaling pathways which is closely associated with fat deposition regulation. We identified circRNA812, circRNA91, and circRNA388 as vital genes in fat deposition by miRNA-circRNA target gene prediction. The functional annotation results of target genes of key circRNAs showed that the signaling pathways mainly included PI3K-Akt and AMPK. We constructed the competing endogenous RNA (ceRNA) regulatory network to study the role of circRNAs in sheep lipid deposition, and circRNA812, circRNA91, and circRNA388 can adsorb more miRNAs. NC_040253.1_5757, as the source of miRNA response element (MRE) among the three, may play an important role during the process of sheep fat deposition.

CONCLUSIONS

Our study gives a systematic examination of the circRNA profiles expressed in sheep subcutaneous fat. These results from this study provide some new basis for understanding circRNA function and sheep fat metabolism.

Comprehensive analysis on the regulation of differentially expressed of mRNA and ncRNA in different ovarian stages of ark shell Scapharca broughtonii

BACKGROUND

Ovarian development is an important prerequisite and basis for animal reproduction. In many vertebrates, it is regulated by multiple genes and influenced by sex steroid hormones and environmental factors. However, relative information is limited in shellfish. To explore the biological functions and molecular mechanisms of mRNA and non-coding RNA that regulate ovarian development in Scapharca broughtonii, we performed whole transcriptome sequencing analysis on ovaries at three developmental stages. Furthermore, the biological processes involved in the differential expression of mRNA and ncRNA were analyzed.

RESULTS

A total of 11,342 mRNAs, 6897 lncRNAs, 135 circRNAs, and 275 miRNAs were differentially expressed. By mapping the differentially expressed RNAs from the three developmental stages of Venn diagram, multiple groups of shared mRNAs and lncRNAs were found to be associated with ovarian development, with some mRNA and ncRNA functions associated with steroid hormone. In addition, we constructed and visualized the lncRNA/circRNA-miRNA-mRNA network based on ceRNA targeting relationships.

CONCLUSIONS

These findings may facilitate our further understanding the mRNA and ncRNAs roles in the regulation of shellfish reproduction.

Integrated analysis of circRNA, lncRNA, miRNA and mRNA to reveal the ceRNA regulatory network of postnatal skeletal muscle development in Ningxiang pig

Introduction: The development of skeletal muscle is regulated by regulatory factors of genes and non-coding RNAs (ncRNAs). Methods: The objective of this study was to understand the transformation of muscle fiber type in the longissimus dorsi muscle of male Ningxiang pigs at four different growth stages (30, 90, 150, and 210 days after birth, n = 3) by histological analysis and whole transcriptome sequencing. Additionally, the study investigated the expression patterns of various RNAs involved in muscle fiber transformation and constructed a regulatory network for competing endogenous RNA (ceRNA) that includes circular RNA (circRNA)/long non-coding RNA (lncRNA)-microRNA (miRNA)-messenger RNA (mRNA). Results: Histomorphology analysis showed that the diameter of muscle fiber reached its maximum at 150 days after birth. The slow muscle fiber transformation showed a pattern of initial decrease followed by an increase. 29,963 circRNAs, 2,683 lncRNAs, 986 miRNAs and 22,411 mRNAs with expression level ≥0 were identified by whole transcriptome sequencing. Furthermore, 642 differentially expressed circRNAs (DEc), 505 differentially expressed lncRNAs (DEl), 316 differentially expressed miRNAs (DEmi) and 6,090 differentially expressed mRNAs (DEm) were identified by differential expression analysis. Functions of differentially expressed mRNA were identified by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). GO enrichment analysis indicates that 40 known genes and 6 new genes are associated with skeletal muscle development. Additionally, KEGG analysis shows that these genes regulate skeletal muscle development via MAPK, FoxO, Hedgehog, PI3K-Akt, Notch, VEGF and other signaling pathways. Through protein-protein interaction (PPI) and transcription factor prediction (TFP), the action mode of skeletal muscle-related genes was explored. PPI analysis showed that there were stable interactions among 19 proteins, meanwhile, TFP analysis predicted 22 transcription factors such as HMG20B, MYF6, MYOD1 and MYOG, and 12 of the 19 interacting proteins were transcription factors. The regulatory network of ceRNA related to skeletal muscle development was constructed based on the correlation of various RNA expression levels and the targeted binding characteristics with miRNA. The regulatory network included 31 DEms, 59 miRNAs, 667 circRNAs and 224 lncRNAs. conclusion: Overall, the study revealed the role of ceRNA regulatory network in the transformation of skeletal muscle fiber types in Ningxiang pigs, which contributes to the understanding of ceRNA regulatory network in Ningxiang pigs during the skeletal muscle development period.

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.

Population Genetic Analysis of Six Chinese Indigenous Pig Meta-Populations Based on Geographically Isolated Regions

The diversification of indigenous pig breeds in China has resulted from multiple climate, topographic, and human cultural influences. The numerous indigenous pig breeds can be geographically divided into six meta-populations; however, their genetic relationships, contributions to genetic diversity, and genetic signatures remain unclear. Whole-genome SNP data for 613 indigenous pigs from the six Chinese meta-populations were obtained and analyzed. Population genetic analyses confirmed significant genetic differentiation and a moderate mixture among the Chinese indigenous pig meta-populations. The North China (NC) meta-population had the largest contribution to genetic and allelic diversity. Evidence from selective sweep signatures revealed that genes related to fat deposition and heat stress response (EPAS1, NFE2L2, VPS13A, SPRY1, PLA2G4A, and UBE3D) were potentially involved in adaptations to cold and heat. These findings from population genetic analyses provide a better understanding of indigenous pig characteristics in different environments and a theoretical basis for future work on the conservation and breeding of Chinese indigenous pigs.

CircHOMER1 inhibits porcine adipogenesis via the miR-23b/SIRT1 axis

Fat deposition is critical to pork quality. However, the mechanism of fat deposition remains to be elucidated. Circular RNAs (circRNAs) are ideal biomarkers and are involved in adipogenesis. Here, we investigated the effect and mechanism of circHOMER1 on porcine adipogenesis in vitro and in vivo. Western blotting, Oil red O staining, and HE staining were used to assess the function of circHOMER1 in adipogenesis. The results showed that circHOMER1 inhibited adipogenic differentiation of porcine preadipocytes and suppressed adipogenesis in mice. Dual-luciferase reporter gene, RIP, and pull-down assays demonstrated that miR-23b directly bound to circHOMER1 and the 3'-UTR of SIRT1. Rescue experiments further illustrated the regulatory relationship among circHOMER1, miR-23b, and SIRT1. Conclusively, we demonstrate that circHOMER1 plays an inhibitory role in porcine adipogenesis through miR-23b and SIRT1. The present study revealed the mechanism of porcine adipogenesis, which may be helpful to improve pork quality.

Genome-Wide Landscape of mRNAs, lncRNAs, and circRNAs during Testicular Development of Yak

Testicular development is a tightly regulated process in mammals. Understanding the molecular mechanisms of yak testicular development will benefit the yak breeding industry. However, the roles of different RNAs, such as mRNA, lncRNA, and circRNA in the testicular development of yak, are still largely unclear. In this study, transcriptome analyses were performed on the expression profiles of mRNAs, lncRNAs, and circRNAs in testis tissues of Ashidan yak at different developmental stages, including 6-months-old (M6), 18-months-old (M18), and 30-months-old (M30). A total of 30, 23, and 277 common differentially expressed (DE) mRNAs, lncRNAs, and circRNAs were identified in M6, M18, and M30, respectively. Furthermore, functional enrichment analysis showed that the common DE mRNAs during the entire developmental process were mainly involved in gonadal mesoderm development, cell differentiation, and spermatogenesis processes. Additionally, co-expression network analysis identified the potential lncRNAs related to spermatogenesis, e.g., TCONS_00087394 and TCONS_00012202. Our study provides new information about changes in RNA expression during yak testicular development, which greatly improves our understanding of the molecular mechanisms regulating testicular development in yaks.

Integrated transcriptome analysis reveals roles of long non-coding RNAs (lncRNAs) in caprine skeletal muscle mass and meat quality

Long non-coding RNAs (lncRNAs) play important roles in the growth and development of skeletal muscle. However, there is limited information on goats. In this study, expression profiles of lncRNAs in Longissimus dorsi muscle from Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with divergent meat yield and meat quality were compared using RNA-sequencing. Based on our previous microRNA (miRNA) and mRNA profiles obtained from the same tissues, the target genes and binding miRNAs of differentially expressed lncRNAs were obtained. Subsequently, lncRNA-mRNA interaction networks and a ceRNA network of lncRNA-miRNA-mRNA were constructed. A total of 136 differentially expressed lncRNAs were identified between the two breeds. Fifteen cis target genes and 143 trans target genes were found for differentially expressed lncRNAs, and they were enriched in muscle contraction, muscle system process, muscle cell differentiation, and p53 signaling pathway. A total of 69 lncRNA-trans target gene pairs were constructed, with close relationship with muscle development, intramuscular fat deposition, and meat tenderness. A total of 16 lncRNA-miRNA-mRNA ceRNA pairs were identified, of which some reportedly associated with skeletal muscle development and fat deposition were found. The study will provide an improved understanding of the roles of lncRNAs in caprine meat yield and meat quality.

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.

Comparison of meat quality and glycolysis potential of two hybrid pigs in three-way hybrid model

With the improvement of consumers' requirements for pork quality, the method of crossbreeding with excellent local pig breeds to improve meat quality is popular. Saba pig has high reproduction rate, good meat quality and high utilization rate of roughage, but its excellent characteristics have not been fully developed and utilized. To promote the development and utilization of Saba pigs and production of high-quality pork, the meat quality traits and glycolysis potential of Duroc × (Landrace × Yorkshire) (DLY), Berkshire × (Duroc × Saba) (BDS), and Duroc × (Berkshire × Saba) (DBS) three-way crossbred pigs were compared. The results showed that DLY had the highest live weight, carcass weight, lean meat percentage, drip loss, glycolysis potential, muscle diameter, and relative mRNA expression levels of type IIb muscle fibers as well as the lowest ultimate pH (p < 0.05). The lightness value of DBS was the highest (p < 0.05). Among the three crossbred pigs, myristic, arachidic, palmitoleic, and eicosenoic acids were the highest in BDS. These results indicated that the carcass traits of local crossbred pigs were worse than those of DLY pigs, but meat quality was markedly higher, with BDS showing the best meat quality.

Genome-Wide Analysis of Circular RNAs Reveals circCHRNG Regulates Sheep Myoblast Proliferation via miR-133/SRF and MEF2A Axis

As relatively new members of the non-coding RNA family, circRNAs play important roles in a variety of biological processes. However, the temporal expression pattern and the function of circRNAs during sheep skeletal muscle development remains unclear. This study aimed to identify circRNAs related to sheep skeletal muscle development and explore their roles in myoblast proliferation. The circRNA expression profiles of longissimus dorsi of sheep from F90, L30, and A3Y were obtained by the RNA-seq method. The function and mechanisms of the novel circCHRNG in muscle satellite cell proliferation were explored using CCK-8 assay, Western blot, qPCR, and dual-luciferase reporter assay. We identified 12,375 circRNAs, including 476, 133, and 233 DEcircRNAs found among three comparative groups. KEGG results showed that DEcircRNAs were enriched in muscle contraction, the regulation of cell proliferation, and the AMPK, insulin, and PI3K-Akt signaling pathways. Notably, a novel circRNA, termed circRNA CHRNG, acts as a miR-133 sponge to promote skeletal muscle satellite cell proliferation. Our study provides a systematic description of circRNAs of ovine skeletal muscle across fetal, lamb, and adult stages. GO and KEGG analyses showed that DEcircRNAs were enriched in multiple pathways associated with muscle development, such as the PI3K-Akt and AMPK signaling pathways. In addition, we propose that circCHRNG acts as a miR-133 sponge to upregulate the expression levels of SRF and MEF2A, thereby promoting myoblast proliferation.

Integrated analysis of microRNAs, circular RNAs, long non-coding RNAs, and mRNAs revealed competing endogenous RNA networks involved in brown adipose tissue whitening in rabbits

BACKGROUND

The brown adipose tissue (BAT) is a target for treating obesity. BAT losses thermogenic capacity and gains a "white adipose tissue-like" phenotype ("BAT whitening") under thermoneutral environments, which is a potential factor causing a low curative effect in BAT-related obesity treatments. Circular RNAs (circRNAs) and long non-coding RNAs (lncRNAs) can act as competing endogenous RNAs (ceRNA) to mRNAs and function in various processes by sponging shared microRNAs (miRNAs). However, the roles of circRNA- and lncRNA-related ceRNA networks in regulating BAT whitening remain litter known.

RESULTS

In this study, BATs were collected from rabbits at day0 (D0), D15, D85, and 2 years (Y2). MiRNA-seq was performed to investigate miRNA changes during BAT whitening. Then, a combined analysis of circRNA-seq and whole-transcriptome sequencing was used for circRNA assembly and quantification during BAT whitening. Our data showed that 1187 miRNAs and 6204 circRNAs were expressed in the samples, and many of which were identified as significantly changed during BAT whitening. Target prediction showed that D0-selective miRNAs were significantly enriched in the Ras, MAPK, and PI3K-Akt signaling pathways, and Y2-selective miRNAs were predicted to be involved in cell proliferation. The cyclization of several circRNAs could form novel response elements of key thermogenesis miRNAs at the back-splicing junction (BSJ) sites, and in combination with a dual-luciferase reporter assay confirmed the binding between the BSJ site of novel_circ_0013792 and ocu-miR-378-5p. CircRNAs and lncRNAs have high cooperativity in sponging miRNAs during BAT whitening. Both circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA triple networks were significantly involved in immune response-associated biological processes. The D15-selective circRNA might promote BAT whitening by increasing the expression of IDH2. The Y2-selective circRNA-related ceRNA network and lncRNA-related ceRNA network might regulate the formation of the WAT-like phenotype of BAT via MAPK and Ras signaling pathways, respectively.

CONCLUSIONS

Our work systematically revealed ceRNA networks during BAT whitening in rabbits and might provide new insight into BAT-based obesity treatments.

Identification and characterization of circRNAs related to meat quality during embryonic development of the longissimus dorsi muscle in two pig breeds

Meat quality, an important economic trait, is regulated by many factors, especially by genetic factors, including coding genes, miRNAs, and lncRNAs. Recent studies have elucidated that circRNAs also play a key role in muscle development and lipid deposition. However, the functions and regulatory mechanisms of circRNAs in meat quality remain mostly unknown. The circRNA expression profiles between Huainan pigs (Chinese indigenous pigs, fat-type, Huainan HN) and Large White pigs (Western commercial pigs, lean-type, LW) in the longissimus dorsi (LD) muscle at 38, 58, and 78 days post conception (dpc) were compared by sequencing. In total, 39,887 circRNAs were identified in 18 samples, and 60, 78, and 86 differentially expressed circRNAs (DECs) were found at the three stages mentioned above between these two breeds. The parent genes of DECs were enriched in myogenesis, proliferation, adipogenesis and muscle fiber-type transition. The circRNA-miRNA interaction networks included 38 DECs and 47 miRNAs, and these miRNAs were involved in muscle development and lipid metabolism. Two shared DECs (circ_0030593 and circ_0032760) of these three stages were selected, their head-to-tail junction sites were validated by Sanger sequencing, and RT‒qPCR results suggested that these two DECs might be involved in intramuscular fat deposition. These findings provide a basis for understanding the role of circRNAs in meat quality.

Integrated analysis of lncRNA and gene expression in longissimus dorsi muscle at two developmental stages of Hainan black goats

It is deemed that meat quality of kids’ is better than that of adults’ for Hainan black goat. Generally, meat quality is affected by many indicators, such as intramuscular fat (IMF) content, muscle fiber diameter and shear force. It is indicated that long non-coding RNAs (lncRNAs) play essential roles in meat quality of goats. However, it is unclear whether and how lncRNAs and genes play their roles in meat quality of Hainan Black goats. Here, we firstly compared the meat quality between two-month-old kids (kids) and adult goats (adults). Then, the lncRNA-seq and RNA-seq data were integrated and analyzed to explore the potential functions of lncRNAs and genes. The results showed that adults’ IMF content and muscle fiber diameter were extremely significantly higher than that of kids (P<0.01). For the sequenced data, average 84,970,398, and 83,691,250 clean reads were obtained respectively for Kids and adults, among which ~96% were mapped to the reference genome of goats. Through analyzing, 18,242 goat annotated genes, 1,429 goat annotated lncRNAs and 2,967 novel lncRNAs were obtained. Analysis of differential expression genes (DEGs) and lncRNAs (DELs) showed that 328 DEGs and 98 DELs existed between kids and adults. Furthermore, functional enrichment analysis revealed that a number of DEGs and DELs were mainly associated with IMF. Primarily, DGAT2 expressed higher in adults than that in kids and CPT1A expressed higher in kids than that in adults. Both of them were overlapped by DEGs and targets of DELs, suggesting the two DEGs and the DELs targeted by the two DEGs might be the potential regulators of goat IMF deposition. Taken together, our results provide basic support for further understanding the function and mechanism of lncRNAs and genes in meat quality of Hainan black goats.

Identification of Candidate Genes and Regulatory Competitive Endogenous RNA (ceRNA) Networks Underlying Intramuscular Fat Content in Yorkshire Pigs with Extreme Fat Deposition Phenotypes

Intramuscular fat (IMF) content is vital for pork quality, serving an important role in economic performance in pig industry. Non-coding RNAs, with mRNAs, are involved in IMF deposition; however, their functions and regulatory mechanisms in porcine IMF remain elusive. This study assessed the whole transcriptome expression profiles of the Longissimus dorsi muscle of pigs with high (H) and low (L) IMF content to identify genes implicated in porcine IMF adipogenesis and their regulatory functions. Hundreds of differentially expressed RNAs were found to be involved in fatty acid metabolic processes, lipid metabolism, and fat cell differentiation. Furthermore, combing co-differential expression analyses, we constructed competing endogenous RNAs (ceRNA) regulatory networks, showing crosstalk among 30 lncRNAs and 61 mRNAs through 20 miRNAs, five circRNAs and 11 mRNAs through four miRNAs, and potential IMF deposition-related ceRNA subnetworks. Functional lncRNAs and circRNAs (such as MSTRG.12440.1, ENSSSCT00000066779, novel_circ_011355, novel_circ_011355) were found to act as ceRNAs of important lipid metabolism-related mRNAs (LEP, IP6K1, FFAR4, CEBPA, etc.) by sponging functional miRNAs (such as ssc-miR-196a, ssc-miR-200b, ssc-miR10391, miR486-y). These findings provide potential regulators and molecular regulatory networks that can be utilized for research on IMF traits in pigs, which would aid in marker-assisted selection to improve pork quality.

Analyzing the interactions of mRNAs, miRNAs and lncRNAs to predict ceRNA networks in bovine cystic follicular granulosa cells

Cross-talk between competitive endogenous RNAs (ceRNAs) may play a critical role in revealing potential mechanism of bovine follicular cysts. Ovarian cyst has always been an intractable scientific problem and has led to considerable economic losses to bovine breeding industry. However, its pathogenesis and molecular mechanisms are still not well understood. Here, this study aimed to investigate the role of non-coding RNAs (ncRNAs) and the ceRNA networks in bovine follicular cyst. Whole transcriptome sequencing of bovine follicular granulosa cells (GCs) was conducted to obtain the expression profiles of mRNAs, lncRNAs and miRNAs. The results for the identified expressions of 8,003 mRNAs, 579 lncRNAs and 205 miRNAs were often altered between cystic and normal follicular GCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed on these differentially expressed mRNAs. Furthermore, the ceRNA network combining mRNAs, miRNAs, and lncRNAs using several bioinformatics methods based on co-expression analysis between the differentially expressed RNAs was conducted. Finally, the lncRNA NONBTAT027373.1-miR-664b-HSD17B7 pathway was verified by dual-luciferase reporting assay and RNA binding protein immunoprecipitation (RIP) assay. LncRNA NONBTAT027373.1 sponged miR-664b in GCs and prevented miR-664b from binding to the HSD17B7 3′-UTR. These results indicated that genes and lncRNAs related to steroid hormone synthesis and energy metabolism could play important roles in the formation of bovine cystic follicles through the ceRNA mechanism and represent candidate targets for further research. This can be used as a practical guideline for promoting healthy and highly efficient development in the bovine industry.

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.

The comprehensive detection of miRNA and circRNA in the regulation of intramuscular and subcutaneous adipose tissue of Laiwu pig

circRNAs, as miRNA sponges, participate in many important biological processes. However, it remains unclear whether circRNAs can regulate lipid metabolism. This study aimed to explore the competing endogenouse RNA (ceRNA) regulatory network that affects the difference between intramuscular fat (IMF) and subcutaneous fat (SCF) deposition, and to screen key circRNAs and their regulatory genes. In this experiment, we identified 265 differentially expressed circRNAs, of which 187 up-regulated circRNA and 78 down-regulated circRNA in IMF. Subsequently, we annotated the function of DEcircRNA's host genes, and found that DEcircRNA's host genes were mainly involved in GO terms (including cellular response to fatty acids, lysophosphatidic acid acyltransferase activity, R-SMAD binding, etc.) and signaling pathways (fatty acid biosynthesis, Citrate cycle, TGF- β Signal pathway) related to adipogenesis, differentiation and lipid metabolism. By constructing a circRNA-miRNA network, we screened out DEcircRNA that can competitively bind to more miRNAs as key circRNAs (circRNA_06424 and circRNA_08840). Through the functional annotation of indirect target genes and protein network analysis, we found that circRNA_06424 affects the expression of PPARD, MMP9, UBA7 and other indirect target genes by competitively binding to miRNAs such as ssc-miR-339-5p, ssc-miR-744 and ssc-miR-328, and participates in PPAR signaling pathway, Wnt signaling pathway, unsaturated fatty acid and other signaling pathways, resulting in the difference of fat deposition between IMF and SCF. This study provide a theoretical basis for further research investigating the differences of lipid metabolism in different adipose tissues, providing potential therapeutic targets for ectopic fat deposition and lipid metabolism diseases.

Identification of critical lncRNAs for milk fat metabolism in dairy cows using WGCNA and the construction of a ceRNAs network

As key regulators, long non-coding RNAs (lncRNAs) play a crucial role in the ruminant mammary gland. However, the function of lncRNAs in milk fat synthesis from dairy cows is largely unknown. In this study, we used the weighted gene co-expression network analysis (WGCNA) to comprehensive analyze the expression profile data of lncRNAs from the group's previous Illumina PE150 sequencing results based on bovine mammary epithelial cells from high- and low-milk-fat-percentage (MFP) cows, and identify core_lncRNAs significantly associated with MFP by module membership (MM) and gene significance (GS). Functional enrichment analysis (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes) of core_lncRNA target genes (co-localization and co-expression) was performed to screen potential lncRNAs regulating milk fat metabolism and further construct an interactive regulatory network of lipid metabolism-related competing endogenous RNAs (ceRNAs). A total of 4876 lncRNAs were used to construct the WGCNA. The MEdarkturquoise module among the 19 modules obtained was significantly associated with MFP (r = 0.78, p-value <0.05) and contained 64 core_lncRNAs (MM > 0.8, GS > 0.4). Twenty-four lipid metabolism-related lncRNAs were identified by core_lncRNA target gene enrichment analysis. TCONS_00054233, TCONS_00152292, TCONS_00048619, TCONS_00033839, TCONS_00153791 and TCONS_00074642 were key candidate lncRNAs for regulating milk fat synthesis. The 22 ceRNAs most likely to be involved in milk fat metabolism were constructed by interaction network analysis, and TCONS_00133813 and bta-miR-2454-5p were located at the network's core. TCONS_00133813_bta-miR-2454-5p_TNFAIP3, TCONS_00133813_bta-miR-2454-5p_ARRB1 and TCONS_00133813_bta-miR-2454-5p_PIK3R1 are key candidate ceRNAs associated with milk fat metabolism. This study provides a framework for the co-expression module of MFP-related lncRNAs in ruminants, identifies several major lncRNAs and ceRNAs that influence milk fat synthesis, and provides a new understanding of the complex biology of milk fat synthesis in dairy cows.

Comparative Transcriptomic Analysis of mRNAs, miRNAs and lncRNAs in the Longissimus dorsi Muscles between Fat-Type and Lean-Type Pigs

In pigs, meat quality and production are two important traits affecting the pig industry and human health. Compared to lean-type pigs, fat-type pigs contain higher intramuscular fat (IMF) contents, better taste and nutritional value. To uncover genetic factors controlling differences related to IMF in pig muscle, we performed RNA-seq analysis on the transcriptomes of the Longissimus dorsi (LD) muscle of Laiwu pigs (LW, fat-type pigs) and commercial Duroc × Landrace × Yorkshire pigs (DLY, lean-type pigs) at 150 d to compare the expression profiles of mRNA, miRNA and lncRNA. A total of 225 mRNAs, 12 miRNAs and 57 lncRNAs were found to be differentially expressed at the criteria of |log2(foldchange)| > 1 and q < 0.05. The mRNA expression of LDHB was significantly higher in the LD muscle of LW compared to DLY pigs with log2(foldchange) being 9.66. Using protein interaction prediction method, we identified more interactions of estrogen-related receptor alpha (ESRRA) associated with upregulated mRNAs, whereas versican (VCAN) and proenkephalin (PENK) were associated with downregulated mRNAs in LW pigs. Integrated analysis on differentially expressed (DE) mRNAs and miRNAs in the LD muscle between LW and DLY pigs revealed two network modules: between five upregulated mRNA genes (GALNT15, FKBP5, PPARGC1A, LOC110258214 and LOC110258215) and six downregulated miRNA genes (ssc-let-7a, ssc-miR190-3p, ssc-miR356-5p, ssc-miR573-5p, ssc-miR204-5p and ssc-miR-10383), and between three downregulated DE mRNA genes (IFRD1, LOC110258600 and LOC102158401) and six upregulated DE miRNA genes (ssc-miR1379-3p, ssc-miR1379-5p, ssc-miR397-5p, ssc-miR1358-5p, ssc-miR299-5p and ssc-miR1156-5p) in LW pigs. Based on the mRNA and ncRNA binding site targeting database, we constructed a regulatory network with miRNA as the center and mRNA and lncRNA as the target genes, including GALNT15/ssc-let-7a/LOC100523888, IFRD1/ssc-miR1379-5p/CD99, etc., forming a ceRNA network in the LD muscles that are differentially expressed between LW and DLY pigs. Collectively, these data may provide resources for further investigation of molecular mechanisms underlying differences in meat traits between lean- and fat-type pigs.

Identification and Functional Prediction of Circular RNAs Related to Growth Traits and Skeletal Muscle Development in Duroc pigs

Porcine skeletal muscle is a highly heterogeneous tissue type, and the Longissimus Dorsi muscle (LDM), as the most economical and physiologically metabolized skeletal muscle in pigs, has always been the focus of research and improvement in pig molecular breeding. Circular RNA, as an important new member of regulatory non-coding RNA after microRNA, has become a frontier hot spot in life science research. This study aims to explore candidate circRNAs related to growth, meat quality, and skeletal muscle development among Duroc pigs with different average daily gain (ADG). Eight pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG), followed by RNA-Seq analysis to identify circRNAs. The results showed that backfat at 6-7 rib (BF) and Intramuscular fat (IMF) content in the H group was significantly lower than L group, but ribeye area (REA) in the H group was higher than in the L group. In RNA-seq, 296 Differentially expressed (DE) circRNAs (157 upregulated and 139 downregulated) were identified and exons flanking long introns are easier to circularize to produce circRNAs. Most of the DE circRNAs were enriched in Quantitative trait locus (QTL) regions related to meat quality and growth traits. In addition, a gene can produce one or more circRNA transcripts. It was also found that the source genes of DE circRNAs were enriched in MAPK, FoXO, mTOR, PI3K-Akt, and Wnt signaling pathways. The results showed that different ADG, carcass, and meat quality traits among half-sibling Duroc pigs with the same diet may be due to the DE circRNAs related to skeletal muscle growth and development.

Comprehensive analysis of differently expression mRNA and non-coding RNAs, and their regulatory mechanisms on relationship in thiram-induced tibial dyschondroplasia in chicken

Thiram pollution is one of the main causes of tibial dyschondroplasia (TD) induced by feed sources. Several studies have speculated that miRNA, circRNA and lncRNA may have significant impact on the development of TD, however, the specific mRNAs and noncoding RNAs and their respective regulatory mechanisms and functions in the development of TD have not been explored. Therefore, in this present study, we screened the differentially expressed mRNA, miRNA, circRNA and lncRNA by whole-transcriptome sequencing (RNA-seq) and differentially expressed genes (DEGs) enrichment, as well as constructed the interaction network among the mRNA-miRNA, mRNA-lncRNA and mRNA-miRNA-circRNA. The sequencing results were verified by fluorescence real-time quantitative PCR (RT-qPCR). The results obtained in this study, revealed that the cells were atrophied and disordered in the TD group, and the expression of BMP6, TGF-β and VEGF were significantly reduced. A total of 141 mRNAs, 10 miRNAs, 23 lncRNAs and 35 circRNAs of DEGs were obtained (p<0.05) Theses DEGs were enriched in the adhere junction and insulin signaling pathways. In addition, the mRNA-miRNA-circRNA network suggested that several pivotal ceRNA showed a regulatory relationship between the transcripts with miRNA, circRNA or lncRNA. Taken together, the results in the present study, represent an insight for further functional research on the ceRNA regulatory mechanism of TD in broilers.

CircRNA, lncRNA, and mRNA profiles of umbilical cord blood exosomes from preterm newborns showing bronchopulmonary dysplasia

Bronchopulmonary dysplasia (BPD) represents a multifactorial chronic pulmonary pathology and a major factor causing premature illness and death. The therapeutic role of exosomes in BPD has been feverishly investigated. Meanwhile, the potential roles of exosomal circRNAs, lncRNAs, and mRNAs in umbilical cord blood (UCB) serum have not been studied. This study aimed to detect the expression profiles of circRNAs, lncRNAs, and mRNAs in UCB-derived exosomes of infants with BPD. Microarray analysis was performed to compare the RNA profiles of UCB-derived exosomes of a preterm newborn with (BPD group) and without (non-BPD, NBPD group) BPD. Then, circRNA/lncRNA-miRNA-mRNA co-expression networks were built to determine their association with BPD. In addition, cell counting kit-8 (CCK-8) assay was used to evaluate the proliferation of lipopolysaccharide (LPS)-induced human bronchial epithelial cells (BEAS-2B cells) and human umbilical vein endothelial cells (HUVECs). The levels of tumor necrosis factor (TNF)-α and interleukin (IL)-1β in LPS-induced BEAS-2B cells and HUVECs were assessed through Western blot analysis. Then, quantitative reverse transcription-polymerase chain reaction assay was used to evaluate the expression levels of four differentially expressed circRNAs (hsa_circ_0086913, hsa_circ_0049170, hsa_circ_0087059, and hsa_circ_0065188) and two lncRNAs (small nucleolar RNA host gene 20 (SNHG20) and LINC00582) detected in LPS-induced BEAS-2B cells or HUVECs. A total of 317 circRNAs, 104 lncRNAs, and 135 mRNAs showed significant differential expression in UCB-derived exosomes of preterm infants with BPD compared with those with NBPD. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted to examine differentially expressed exosomal circRNAs, lncRNAs, and mRNAs. The results showed that the GO terms and KEGG pathways mostly involving differentially expressed exosomal RNAs were closely associated with endothelial or epithelial cell development. In vitro, CCK-8 and Western blot assays revealed that LPS remarkably inhibited the viability and promoted inflammatory responses (TNF-α and IL-1β) of BEAS-2B cells or HUVECs. The expression levels of circRNAs hsa_circ_0049170 and hsa_circ_0087059 were upregulated in LPS-induced BEAS-2B cells; the expression level of hsa_circ_0086913 was upregulated and that of hsa_circ_0065188 was downregulated in LPS-induced HUVECs. Moreover, the expression level of lncRNA SNHG20 was upregulated and that of LINC00582 was downregulated in LPS-induced BEAS-2B cells. Further, 455 circRNA/lncRNA-miRNA-mRNA interaction networks were predicted, including hsa_circ_0086913/hsa-miR-103a-3p/transmembrane 4 L six family member 1 (TM4SF1) and lncRNA-SNHG20/hsa-miR-6720-5p/spermine synthase (SMS) networks, which may take part in BPD.

CONCLUSION

This study provided a systematic perspective on UCB-derived exosomal circRNAs and lncRNAs and laid an important foundation for further investigating the potential biological functions of exosomal circRNAs and lncRNAs in BPD.

WHAT IS KNOWN

• BPD represents a multifactorial chronic pulmonary pathology and a major factor causing premature illness and death. • The therapeutic role of exosomes in BPD has been feverishly investigated, and exosomal RNAs were ignored.

WHAT IS NEW

• The profiles of UCB-derived exosomal circRNAs, lncRNAs, and mRNAs were performed. • Several differentially expressed circRNAs and lncRNAs were identified in LPS-induced BEAS-2B cells and HUVECs.

Genome-Wide Identification and Characterization of Circular RNAs during Skeletal Muscle Development in Meat Rabbits

Simple Summary

Our knowledge of circRNAs regulating skeletal muscle development remains largely unknown in meat rabbits. Therefore, we collected the leg muscle tissues of ZIKA rabbits at three key growth stages. A combination of circRNA assembly from a circRNA-seq library and the whole-transcriptome sequencing data identified credible circRNAs in our samples. We found these circRNAs were more conserved between rabbits and humans than between rabbits and mice. A prediction of circRNA–microRNA–mRNAs networks revealed that circRNAs might be the regulators that mainly functioned in rabbits’ muscle neuron development and metabolic processes. Our work provides a catalog of circRNAs regulating skeletal muscle development at key growth stages in rabbits and might give a new insight into rabbit breeding.

Abstract

Skeletal muscle development plays a vital role in muscle quality and yield in meat rabbits. Circular RNAs (circRNAs) are a new type of single-stranded endogenous non-coding RNAs involved in different biological processes. However, our knowledge of circRNAs regulating skeletal muscle development remains largely unknown in meat rabbits. In this study, we collected the leg muscle tissues of ZIKA rabbits at three key growth stages. By performing whole-transcriptome sequencing, we found the sequential expression of day 0- (D0-), D35-, and D70-selective mRNAs mainly functioned in muscle development, nervous development, and immune response during skeletal muscle development, respectively. Then, a combination of circRNA assembly from a circRNA-seq library and the whole-transcriptome sequencing data identified 6845 credible circRNAs in our samples. Most circRNAs were transcribed from exons of known genes, contained few exons, and showed short length, and these circRNAs were more conserved between rabbits and humans than between rabbits and mice. The upregulated circRNAs, which were synchronously changed with host genes, primarily played roles in MAPK signaling pathways and fatty acid biosynthesis. The prediction of circRNA–microRNA–mRNAs networks revealed that circRNAs might be the regulators that mainly functioned in rabbits’ muscle neuron development and metabolic processes. Our work provides a catalog of circRNAs regulating skeletal muscle development at key growth stages in rabbits and might give a new insight into rabbit breeding.

Identification of potential miRNA-mRNA regulatory network and the key miRNAs in intramuscular and subcutaneous adipose

Intramuscular fat (IMF) is an important indicator for evaluating meat quality. Breeds with high IMF content are often accompanied by high subcutaneous fat (SCF), severely affecting the meat rate of pigs. Studying the mechanisms of miRNAs related to lipogenesis and lipid metabolism has important implications for pig breeding. We constructed two small RNA libraries from intramuscular and subcutaneous fat to evaluate the patterns of lipogenesis in Laiwu pig, a Chinese breed. A total of 286 differentially expressed miRNAs (DEmiRNAs), including 193 known miRNA and 93 novel miRNAs, were identified from two types of adipose. GO and KEGG enrichment analysis for DEmiRNAs showed that their target genes involved in many adipogenesis and lipid metabolism biological processes and signaling pathways, such as Wnt signaling pathway,MAPK signaling pathway, Hippo signaling pathway, PI3K-Akt signaling pathway, Melanogenesis, Signaling pathways regulating pluripotency of stem cells and so on. Then, we constructed a miRNA-mRNA interaction network to find out which miRNAs were the key miRNAs of regulation in Wnt signaling pathway. In this pathway, miR-331-3p, miR-339-5p, miR-874 and novel346_mature target PPARD, WNT10B, RSPO3, WNT2B. This study provides a theoretical basis for further understanding the post-transcriptional regulation mechanism of meat quality formation and predicting and treating diseases caused by ectopic fat.

lncRNA-miRNA-mRNA ceRNA Network Involved in Sheep Prolificacy: An Integrated Approach

Understanding the molecular pattern of fertility is considered as an important step in breeding of different species, and despite the high importance of the fertility, little success has been achieved in dissecting the interactome basis of sheep fertility. However, the complex mechanisms associated with prolificacy in sheep have not been fully understood. Therefore, this study aimed to use competitive endogenous RNA (ceRNA) networks to evaluate this trait to better understand the molecular mechanisms responsible for fertility. A competitive endogenous RNA (ceRNA) network of the corpus luteum was constructed between Romanov and Baluchi sheep breeds with either good or poor genetic merit for prolificacy using whole-transcriptome analysis. First, the main list of lncRNAs, miRNAs, and mRNA related to the corpus luteum that alter with the breed were extracted, then miRNA−mRNA and lncRNA−mRNA interactions were predicted, and the ceRNA network was constructed by integrating these interactions with the other gene regulatory networks and the protein−protein interaction (PPI). A total of 264 mRNAs, 14 lncRNAs, and 34 miRNAs were identified by combining the GO and KEGG enrichment analyses. In total, 44, 7, 7, and 6 mRNAs, lncRNAs, miRNAs, and crucial modules, respectively, were disclosed through clustering for the corpus luteum ceRNA network. All these RNAs involved in biological processes, namely proteolysis, actin cytoskeleton organization, immune system process, cell adhesion, cell differentiation, and lipid metabolic process, have an overexpression pattern (Padj < 0.01). This study increases our understanding of the contribution of different breed transcriptomes to phenotypic fertility differences and constructed a ceRNA network in sheep (Ovis aries) to provide insights into further research on the molecular mechanism and identify new biomarkers for genetic improvement.

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.

circMEF2D Negatively Regulated by HNRNPA1 Inhibits Proliferation and Differentiation of Myoblasts via miR-486-PI3K/AKT Axis

Circular RNA (circRNA) is a form of endogenous RNA that can regulate gene expression and participate in the regulation of myogenesis. However, the molecular mechanisms and potential roles of circRNAs in bovine muscle development remain largely unknown. Nevertheless, the RNA splicing factors regulating the biogenesis of bovine circRNA have not yet been characterized. In this study, we identified a novel circRNA, circMEF2D, formed by back-splicing of constitutive exons (exons 5-7) of the bovine MEF2D gene. Functional assays showed that circMEF2D inhibited the proliferation and differentiation of bovine myoblasts. Importantly, we showed that circMEF2D regulated the PI3K-AKT signaling pathway through direct and competitive binding to miR-486. Furthermore, to explore the formation mechanism of circMEF2D, we explored the MEF2D gene alternative splicing progress. Four alternative linear variants of MEF2D were found. Due to its role in alternative splicing, the RNA-binding protein HNRNPA1 was selected for further study and the modulation of HNRNPA1 levels showed that it negatively regulated both back-splicing and linear splicing of MEF2D gene. Overall, in addition to the characterization of bovine circRNAs, these findings revealed the crucial role of HNRNPA1 in MEF2D gene alternative splicing and demonstrated a regulatory circMEF2D-miR-486-PI3K-AKT axis.

Effects of swine microRNA mimics on lipopolysaccharide (LPS) induced inflammatory changes in 3D4/21 cells

There have been limited studies focused on validation of swine microRNAs (miRNA) with mRNA targets. The objective of this study was to validate a defined set of targets using artificial miRNA mimics transfected into cell lines to confirm specific targets of endogenous miRNAs after administration of Escherichia coli lipopolysaccharide (LPS). Sixteen hours after mimic transfection of 3D4/21 cell lines, the cells were stimulated with 1 μg/ml LPS or phosphate-buffered saline (PBS). The cells were harvested and collected at 0, 1, 3, and 8 h post administration. The selected genes DAD1, IL8, and ESR, which are involved in known pathways of inflammation. and are predicted or validated human targets of either miR-146a, let-7a, or miR-22-3p. These were then evaluated by quantitative real-time-PCR (qRT-PCR) to verify microRNA-mRNA interaction in swine. Using the ROX reference dye, mRNA changes in expression were assessed using the comparative CT Method (ΔΔCT method) for normalization against the PBS control group. DAD1 and ESR1 were negatively regulated by miR-22-3p and miR-146a-5p, respectively in 3D4/21 cells after LPS stimulation. However, miR-146a-5p may play an indirect positive regulatory role of both DAD1 and IL8 mRNA expression. Furthermore, we found an inverse relationship between LPS stimulation compared with the let-7a-5p overexpression with DAD1. Our inflammation study provides new evidence on the roles and predicted targets of miR-146a, let-7a, and miR-22-3p in swine.

Insight into mechanisms of pig lncRNA FUT3-AS1 regulating E. coli F18-bacterial diarrhea

Escherichia coli F18 is a common conditional pathogen that is associated with a variety of infections in humans and animals. LncRNAs have emerged as critical players in pathogen infection, but their role in the resistance of the host to bacterial diarrhea remains unknown. Here, we used piglets as animal model and identified an antisense lncRNA termed FUT3-AS1 as a host regulator related to E. coli F18 infection by RNA sequencing. Downregulation of FUT3-AS1 expression contributed to the enhancement of E. coli F18 resistance in IPEC-J2 cells. FUT3-AS1 knockdown reduced FUT3 expression via decreasing the H4K16ac level of FUT3 promoter. Besides, the FUT3-AS1/miR-212 axis could act as a competing endogenous RNA to regulate FUT3 expression. Functional analysis demonstrated that target FUT3 plays a vital role in the resistance of IPEC-J2 cells to E. coli F18 invasion. A Fut3-knockout mice model was established and Fut3-knockout mice obviously improved the ability of resistance to bacterial diarrhea. Interestingly, FUT3 could enhance E. coli F18 susceptibility by activating glycosphingolipid biosynthesis and toll-like receptor signaling which are related to receptor formation and immune response, respectively. In summary, we have identified a novel biomarker FUT3-AS1 that modulates E. coli F18 susceptibility via histone H4 modifications or miR-212/FUT3 axis, which will provide theoretical guidance to develop novel strategies for combating bacterial diarrhea in piglets.

Effect of a Plateau Environment on the Oxidation State of the Heart and Liver through AMPK/p38 MAPK/Nrf2-ARE Signaling Pathways in Tibetan and DLY Pigs

This study evaluated the effect of a plateau environment on the heart and liver oxidation state of Tibetan pigs (TPs) and DLY pigs through analyzing AMPK, p38 MAPK, and Nrf2 signaling pathways. Twelve 120-day-old TPs and twelve 120-day-old DLY pigs were randomly divided into two groups in a plateau environment for three weeks. Exposed to a plateau environment, TPs exhibited a lower (p < 0.05) malondialdehyde level but higher (p < 0.05) glutathione, superoxide dismutase (SOD), glutathione peroxidase (GSH-Px), and total antioxidant capacity (T-AOC) activities in the liver and heart than those observed in DLY pigs. TPs also showed higher (p < 0.05) mRNA levels of SOD and GSH-Px in the liver and heart compared with those of DLY pigs. The TPs showed higher (p < 0.05) mRNA and protein levels of AMPK and Nrf2 in the liver and heart compared with those of DLY pigs. Furthermore, TPs showed higher (p < 0.05) mRNA and protein levels of p38 MAPK in the heart and higher mRNA levels of p38 MAPK in the liver compared with those of DLY pigs under a plateau environment. In summary, TPs possess a stronger antioxidant capacity in the heart and liver than that of DLY pigs in a plateau environment through AMPK/p38 MAPK/Nrf2-ARE signaling pathways.

Comprehensive Transcriptome Analysis of Follicles from Two Stages of the Estrus Cycle of Two Breeds Reveals the Roles of Long Intergenic Non-Coding RNAs in Gilts

Simple Summary

This study provides new perspectives about the roles of lincRNAs in the estrus expression of gilts, which is correlated with ovarian steroid hormone and follicular development. Follicular tissues from two stages of the estrus cycle of Large White and Mi gilts were used for RNA-seq. Some genes and lincRNAs related to estrus expression in pigs were discovered. PPI and ceRNA networks related to the estrus expression were constructed. These results suggest that the estrus expression may be affected by lincRNAs and their target genes.

Abstract

Visible and long-lasting estrus expression of gilts and sows effectively sends a mating signal. To reveal the roles of Long Intergenic Non-coding RNAs (lincRNAs) in estrus expression, RNA-seq was used to investigate the lincRNAs expression of follicular tissues from Large White gilts at diestrus (LD) and estrus (LE), and Chinese Mi gilts at diestrus (MD) and estrus (ME). Seventy-three differentially expressed lincRNAs (DELs) were found in all comparisons (LE vs. ME, LD vs. LE, and MD vs. ME comparisons). Eleven lincRNAs were differentially expressed in both LD vs. LE and MD vs. ME comparisons. Fifteen DELs were mapped onto the pig corpus luteum number Quantitative Trait Loci (QTL) fragments. A protein–protein interaction (PPI) network that involved estrus expression using 20 DEGs was then constructed. Interestingly, three predicted target DEGs (PTGs) (CYP19A1 of MSTRG.10910, CDK1 of MSTRG.10910 and MSTRG.23984, SCARB1 of MSTRG.1559) were observed in the PPI network. A competitive endogenous RNA (ceRNA) network including three lincRNAs, five miRNAs, and five genes was constructed. Our study provides new insight into the lincRNAs associated with estrus expression and follicular development in gilts.

Spatiotemporal Regulation of Circular RNA Expression during Liver Development of Chinese Indigenous Ningxiang Pigs

BACKGROUND

There have been many studies on the relationship between circRNAs and fat deposition. Although the liver is a central organ for fat metabolism, there are few reports on the relationship between circRNAs in the liver and fat deposition.

METHODS

In this study, we systematically analyzed circular RNAs in the liver of Ningxiang pigs, at four time points after birth (30 days, 90 days, 150 days and 210 days).

RESULTS

A total of 3705 circRNAs were coexpressed in four time periods were found, and KEGG analysis showed that the significantly upregulated pathways were mainly enriched in lipid metabolism and amino acid metabolism, while significantly downregulated pathways were mainly related to signal transduction, such as ECM-receptor interaction, MAPK signaling pathway, etc. Short time-series expression miner (STEM) analysis showed multiple model spectra that were significantly enriched over time in the liver. By constructing a competing endogenous RNA (ceRNA) regulatory network, 9187 pairs of networks related to the change in development time were screened.

CONCLUSIONS

The expression profiles of circRNAs in Ningxiang pig liver were revealed at different development periods, and it was determined that there is differential coexpression. Through enrichment analysis of these circRNAs, it was revealed that host genes were involved in metabolism-related signaling pathways and fatty acid anabolism. Through STEM analysis, many circRNAs involved in fat metabolism, transport, and deposition pathways were screened, and the first circRNA-miRNA-mRNA regulation network map in Ningxiang pig liver was constructed. The highly expressed circRNAs related to fat deposition were verified and were consistent with RNA-Seq results.

Knockdown of CDR1as Decreases Differentiation of Goat Skeletal Muscle Satellite Cells via Upregulating miR-27a-3p to Inhibit ANGPT1

Myogenesis is a complex process controlled by several coding and non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) that are known to function as endogenous microRNAs (miRNAs) sponges. Cerebellar Degeneration-Related protein 1 antisense (CDR1as) is the most spotlighted circRNA that is known as an miR-7 sponge, which has bloomed circRNAs’ research in animal disease and physiology. Here, we screened for miRNAs and mRNA associated with CDR1as and further characterized their regulatory function during muscle differentiation. We found that a total of 43 miRNAs (including miR-107-3p, miR-125b-5p, miR-140-5p, miR-29a-3p, and miR-27a-3p upregulated) and 789 mRNAs (including ANGPT1, E2F2, CCN1, FGFR1, and MEF2C downregulated) were differentially expressed in goat skeletal muscle satellite cells (SMSCs). Further, knockdown of CDR1as and ANGPT1 inhibited SMSCs differentiation. miR-27a-3p was differentially upregulated after the knockdown of CDR1as in SMSCs. Overexpressed miR-27a-3p decreased SMSCs differentiation. Via RNAhybrid and luciferase, miR-27a-3p was identified to regulate ANGPT1. We discovered that miR-27a-3p has an inverse relationship with CDR1as and decreases the expression level of ANGPT1 during SMSCs differentiation. In summary, our study demonstrates that siCDR1as inhibits myoblast differentiation by downregulating ANGPT1 mRNA via miR-27a-3p in SMSCs.

Transcriptomic Profiling of Circular RNAs in the Goat Rumen During Fetal and Prepubertal Period

Circular RNAs (circRNAs) are key regulatory factors with vital functions in various biological activities. However, little has been reported concerning the genetic regulation of circRNAs during rumen development in goats. The aim of this study was to identify the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. Histological morphology showed that from the fetal period (days 60 and 135 of gestation) to the prepuberal period (days 60 and 150 of age) the rumen papilla developed gradually, and the thickness of the rumen muscular layer increased. A total of 11,149 circRNAs were identified in the four development stages by RNA-sequencing. From this, 1,518 were differentially expressed circRNAs (DECs). Fifty-eight DECs were up-regulated from 60 to 135 days of gestation, and 93 from day 135 of pregnancy to 30 days after birth. A large proportion (598) of DECs were down-regulated from day 135 of gestation to 30 days after birth. The expression levels of six randomly selected circRNAs were validated by qPCR, and their back-splicing junction (BSJ) sites were also confirmed. Ontology and pathway analyses revealed that the parental genes of DECs were mainly involved in the signaling pathways related to cell proliferation and apoptosis. The interaction network of circRNAs with their target miRNAs showed its involvement in cell proliferation and apoptosis signaling pathways. In conclusion, we identified the genome-wide expression profiles of circRNAs in the rumen of goats during fetal development and before and after weaning. These results provide a basis for further study on the regulatory effect of circRNAs on the development of rumen tissues.

The Competitive Endogenous RNA (ceRNA) Regulation in Porcine Alveolar Macrophages (3D4/21) Infected by Swine Influenza Virus (H1N1 and H3N2)

H1N1 and H3N2 are the two most common subtypes of swine influenza virus (SIV). They not only endanger the pig industry, but are also a huge risk of zoonotic diseases. However, the molecular mechanism and regulatory network of pigs (hosts) against influenza virus infection are still unclear. In this study, porcine alveolar macrophage cell (3D4/21) models infected by swine influenza virus (H1N1 and H3N2) were constructed. The expression profiles of miRNAs, mRNAs, lncRNAs and circRNAs after H1N1 and H3N2 infected 3D4/21 cells were revealed in this study. Then, two ceRNAs (TCONS_00166432-miR10391-MAN2A1 and novel_circ_0004733-miR10391-MAN2A1) that regulated H1N1 and H3N2 infection in 3D4/21 cells were verified by the methods of bioinformatics analysis, gene overexpression, gene interference, real-time quantitative PCR (qPCR), dual luciferase activity assay and RNA immunoprecipitation (RIP). In addition, the important candidate molecules (miR-10391, TCONS_00166432, and novel_circ_0004733) were identified by qPCR and enzyme linked immunosorbent assay (ELISA). Finally, the regulatory effect and possible molecular mechanism of the target gene MAN2A1 were identified by the methods of gene interference, qPCR, Western blot and ELISA. The results of this study suggested that TCONS_00166432 and novel_circ_0004733 could competitively bind miR-10391 to target the MAN2A1 gene to regulate swine influenza virus infecting 3D4/21 cells. This study reported for the first time the ceRNA networks involved in the regulation of the swine influenza virus infecting 3D4/21 cells, which provided a new insight into the molecular mechanism of 3D4/21 cells against swine influenza virus infection.

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.