The spatio-temporal expression patterns of Circular RNA (circRNA) across organs and developmental stages are critical for its function and evolution analysis. However, they remain largely unclear in mammals. Here, we comprehensively analysed circRNAs in nine organs and three skeletal muscles of Guizhou miniature pig (S. scrofa), a widely used biomedical model animal. We identified 5,934 circRNAs and analysed their molecular properties, sequence conservation, spatio-temporal expression pattern, potential function, and interaction with miRNAs. S. scrofa circRNAs show modest sequence conservation with human and mouse circRNAs, are flanked by long introns, exhibit low abundance, and are expressed dynamically in a spatio-temporally specific manner. S. scrofa circRNAs show the greatest abundance and complexity in the testis. Notably, 31% of circRNAs harbour well-conserved canonical miRNA seed matches, suggesting that some circRNAs act as miRNAs sponges. We identified 149 circRNAs potentially associated with muscle growth and found that their host genes were significantly involved in muscle development, contraction, chromatin modification, cation homeostasis, and ATP hydrolysis-coupled proton transport; moreover, this set of genes was markedly enriched in genes involved in tight junctions and the calcium signalling pathway. Finally, we constructed the first public S. scrofa circRNA database, allowing researchers to query comprehensive annotation, expression, and regulatory networks of circRNAs.

Tissue-specific miRNAs (TS miRNA) specifically expressed in particular tissues play an important role in tissue identity, differentiation and function. However, transcription factor (TF) and TS miRNA regulatory networks across multiple tissues have not been systematically studied. Here, we manually extracted 116 TS miRNAs and systematically investigated the regulatory network of TF-TS miRNA in 12 human tissues. We identified 2,347 TF-TS miRNA regulatory relations and revealed that most TF binding sites tend to enrich close to the transcription start site of TS miRNAs. Furthermore, we found TS miRNAs were regulated widely by non-tissue specific TFs and the tissue-specific expression level of TF have a close relationship with TF-genes regulation. Finally, we describe TSmiR (http://bioeng.swjtu.edu.cn/TSmiR), a novel and web-searchable database that houses interaction maps of TF-TS miRNA in 12 tissues. Taken together, these observations provide a new suggestion to better understand the regulatory network and mechanisms of TF-TS miRNAs underlying different tissues.

Despite modest sequence conservation and rapid evolution, long non-coding RNAs (lncRNAs) appear to be conserved in expression pattern and function. However, analysis of lncRNAs across tissues and developmental stages remains largely uncharacterized in mammals. Here, we systematically investigated the lncRNAs of the Guizhou miniature pig (Sus scrofa), which was widely used as biomedical model. We performed RNA sequencing across 9 organs and 3 developmental skeletal muscle, and developed a filtering pipeline to identify 10,813 lncRNAs (9,075 novel). Conservation patterns analysis revealed that 57% of pig lncRNAs showed homology to humans and mice based on genome alignment. 5,455 lncRNAs exhibited typical hallmarks of regulatory molecules, such as high spatio-temporal specificity. Notably, conserved lncRNAs exhibited higher tissue specificity than pigspecific lncRNAs and were significantly enriched in testis and ovary. Weighted co-expression network analysis revealed a set of conserved lncRNAs that are likely involved in postnatal muscle development. Based on the high degree of similarity in the structure, organization, and dynamic expression of pig lncRNAs compared with human and mouse lncRNAs, we propose that these lncRNAs play an important role in organ physiology and development in mammals. Our results provide a resource for studying animal evolution, morphological complexity, breeding, and biomedical research.Intensive transcriptome sequencing, also known as deep sequencing, has led to the discovery that mammalian genomes encode a vast range of non-protein-coding RNAs (ncRNAs) that differ in size and level of conservation 1,2 . The proportion of ncRNAs in an organism's genome has a direct correlation with its developmental complexity 3 . ncRNAs are generally classified into many different RNA types, including microRNA (miRNAs), Piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), small interfering RNAs (siRNAs), and long noncoding RNA (lncRNAs). LncRNAs are defined as transcribed RNA fragments > 200 bp, and they do not have open reading frames of > 100 amino acids. Several recent studies have shown mammalian lncRNAs to be heterogeneous and diverse as well as critically important in cellular function, development, and disease via their transcriptional and posttranscriptional regulation of gene expression 4 . In this study, we aimed to further elucidate the origin, evolution, and function of mammalian ncRNAs, the genome of S. scrofa 5 , a species widely used in medical research, by analyzing the content and function of its lncRNA.LncRNAs have long been ascribed an important role in the evolution of complex traits, and recent studies have revealed that thousands of lncRNAs are evolutionarily conserved in mammals, though not to the same extent as

DNA methylation is important for the epigenetic regulation of gene expression and plays a critical role in mammalian development. However, the dynamic regulation of genome-wide DNA methylation in skeletal muscle development remains largely unknown. Here, we generated the first single-base resolution DNA methylome and transcriptome maps of porcine skeletal muscle across 27 developmental stages. The overall methylation level decreased from the embryo to the adult, which was highly correlated with the downregulated expression of DNMT1 and an increase in partially methylated domains. Notably, we identified over 40 000 developmentally differentially methylated CpGs (dDMCs) that reconstitute the developmental trajectory of skeletal muscle and associate with muscle developmental genes and transcription factors (TFs). The dDMCs were significantly under-represented in promoter regulatory regions but strongly enriched as enhancer histone markers and in chromatin-accessible regions. Integrative analysis revealed the negative regulation of both promoter and gene body methylation in genes associated with muscle contraction and insulin signaling during skeletal muscle development. Mechanistically, DNA methylation affected the expression of muscle-related genes by modulating the accessibly of upstream myogenesis TF binding, indicating the involvement of the DNA methylation/SP1/IGF2BP3 axis in skeletal myogenesis. Our results highlight the function and regulation of dynamic DNA methylation in skeletal muscle development.

MicroRNAs (miRNAs), which are short (22–24 base pairs), non-coding RNAs, play critical roles in myogenesis. Using Solexa deep sequencing, we detected the expression levels of 229 and 209 miRNAs in swine skeletal muscle at 90 days post-coitus (E90) and 100 days postnatal (D100), respectively. A total of 138 miRNAs were up-regulated on E90, and 31 were up-regulated on D100. Of these, 9 miRNAs were selected for the validation of the small RNA libraries by quantitative RT-PCR (RT-qPCR). We found that miRNA-21 was down-regulated by 17-fold on D100 (P<0.001). Bioinformatics analysis suggested that the transforming growth factor beta-induced (TGFβI) gene was a potential target of miRNA-21. Both dual luciferase reporter assays and western blotting demonstrated that the TGFβI gene was regulated by miRNA-21. Co-expression analysis revealed that the mRNA expression levels of miRNA-21 and TGFβI were negatively correlated (r = -0.421, P = 0.026) in skeletal muscle during the 28 developmental stages. Our results revealed that more miRNAs are expressed in prenatal than in postnatal skeletal muscle. The miRNA-21 is a novel myogenic miRNA that is involved in skeletal muscle development and regulates PI3K/Akt/mTOR signaling by targeting the TGFβI gene.

No abstract

MicroRNAs (miRNAs) play a vital role in muscle development by binding to messenger RNAs (mRNAs). Based on prenatal skeletal muscle at 33, 65 and 90 days post-coitus (dpc) from Landrace, Tongcheng and Wuzhishan pigs, we carried out integrated analysis of miRNA and mRNA expression profiling. We identified 33, 18 and 67 differentially expressed miRNAs and 290, 91 and 502 mRNA targets in Landrace, Tongcheng and Wuzhishan pigs, respectively. Subsequently, 12 mRNAs and 3 miRNAs differentially expressed were validated using quantitative real-time PCR (qPCR), and 5 predicted miRNA targets were confirmed via dual luciferase reporter or western blot assays. We identified a set of miRNAs and mRNA genes differentially expressed in muscle development. Gene ontology (GO) enrichment analysis suggests that the miRNA targets are primarily involved in muscle contraction, muscle development and negative regulation of cell proliferation. Our data indicated that more mRNAs are regulated by miRNAs at earlier stages than at later stages of muscle development. Landrace and Tongcheng pigs also had longer phases of myoblast proliferation than Wuzhishan pigs. This study will be helpful to further explore miRNA-mRNA interactions in myogenesis and aid to uncover the molecular mechanisms of muscle development and phenotype variance in pigs.

The pig is an important source of animal protein, and is also an ideal model for human disease. There are significant differences in growth rate, muscle mass, and meat quality between different breeds. To understand the molecular mechanisms underlying porcine skeletal muscle phenotypes, we performed mRNA and miRNA profiling of muscle from three different breeds of pig, Landrace (lean-type), Tongcheng (obese-type), and Wuzhishan (mini-type) by Solexa sequencing. Forty-three genes and 106 miRNAs were differentially expressed between Landrace and Tongcheng pigs, 92 genes and 151 miRNAs were differentially expressed between Tongcheng and Wuzhishan pigs, and 145 genes and 156 miRNAs were differential expressed between Landrace and Wuzhishan pigs. Gene ontology analysis suggested that genes differentially expressed between Landrace and Tongcheng pigs were mainly involved in the biological processes of oxidative stress and muscle organ development. Meanwhile, for Tongcheng vs Wuzhishan and Landrace vs Wuzhishan pigs, the differentially expressed genes were involved in fatty acid metabolism, oxidative stress, muscle contraction, and muscle organ development, processes that are closely related to meat quality. To investigate the molecular mechanisms underlying meat quality diversity based on differentially expressed genes and miRNAs, interaction networks were constructed, according to target prediction results and integration analysis of up-regulated genes with down-regulated miRNAs or down-regulated genes with up-regulated miRNAs. Our findings identify candidate genes and miRNAs associated with muscle development and indicate their potential roles in muscle phenotype variance between different pig breeds. These results serve as a foundation for further studies on muscle development and molecular breeding.