The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

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.

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.

MicroRNAs are evolutionarily conserved small RNAs that post-transcriptionally regulate gene expression and have emerged as critical regulators of skeletal muscle development. Here, we identified miR-148a as a novel myogenic microRNA that mediated myogenic differentiation. The expression levels of miR-148a increased during C2C12 myoblast differentiation. Overexpression of miR-148a significantly promoted myogenic differentiation of both C2C12 myoblast and primary muscle cells. Blocking the function of miR-148a with a 2'-O-methylated antisense oligonucleotide inhibitor repressed C2C12 myoblast differentiation. Using a bioinformatics approach, we identified Rho-associated coiled-coil containing protein kinase 1 (ROCK1), a known inhibitor of myogenesis, as a target of miR-148a. A dual-luciferase reporter assay was used to demonstrate that miR-148a directly targeted the 3'-UTR of ROCK1. In addition, the overexpression of miR-148a decreased the protein expression of ROCK1 in C2C12 myoblast and primary muscle cells. Furthermore, ROCK1 inhibition with specific siRNA leaded to accelerated myogenic differentiation progression, underscoring a negative regulatory function of ROCK1 in myogenesis. Therefore, our results revealed a novel mechanism in which miR-148a positively regulates myogenic differentiation via ROCK1 down-regulation.

Background: Obese and lean pig breeds show obvious differences in muscle growth; however, the molecular mechanism underlying phenotype variation remains unknown. Prenatal muscle development programs postnatal performance. Here, we describe a genome-wide analysis of differences in prenatal skeletal muscle between Tongcheng (a typical indigenous Chinese breed) and Landrace (a leaner Western breed) pigs. Results: We generated transcriptome profiles of skeletal muscle from Tongcheng and Landrace pigs at 33, 65 and 90 days post coitus (dpc), using long serial analysis of gene expression (LongSAGE). We sequenced 317,115 LongSAGE tags and identified 1,400 and 1,201 differentially expressed transcripts during myogenesis in Tongcheng and Landrace pigs, respectively. From these, the Gene Ontology processes and expression patterns of these differentially expressed genes were constructed. Most of the genes showed different expression patterns in the two breeds. We also identified 532, 653 and 459 transcripts at 33, 65 and 90 dpc, respectively, that were differentially expressed between the two breeds. Growth factors, anti-apoptotic factors and genes involved in the regulation of protein synthesis were up-regulated in Landrace pigs. Finally, 12 differentially expressed genes were validated by quantitative PCR. Conclusion: Our data show that gene expression phenotypes differ significantly between the two breeds. In particular, a slower muscle growth rate and more complicated molecular changes were found in Tongcheng pigs, while genes responsible for increased cellular growth and myoblast survival were up-regulated in Landrace pigs. Our analyses will assist in the identification of candidate genes for meat production traits and elucidation of the development of prenatal skeletal muscle in mammals.

MicroRNAs (miRNAs) are short, single-stranded non-coding RNAs that repress their target genes by binding their 3′ UTRs. These RNAs play critical roles in myogenesis. To gain knowledge about miRNAs involved in the regulation of myogenesis, porcine longissimus muscles were collected from 18 developmental stages (33-, 40-, 45-, 50-, 55-, 60-, 65-, 70-, 75-, 80-, 85-, 90-, 95-, 100- and 105-day post-gestation fetuses, 0 and 10-day postnatal piglets and adult pigs) to identify miRNAs using Solexa sequencing technology. We detected 197 known miRNAs and 78 novel miRNAs according to comparison with known miRNAs in the miRBase (release 17.0) database. Moreover, variations in sequence length and single nucleotide polymorphisms were also observed in 110 known miRNAs. Expression analysis of the 11 most abundant miRNAs were conducted using quantitative PCR (qPCR) in eleven tissues (longissimus muscles, leg muscles, heart, liver, spleen, lung, kidney, stomach, small intestine and colon), and the results revealed that ssc-miR-378, ssc-miR-1 and ssc-miR-206 were abundantly expressed in skeletal muscles. During skeletal muscle development, the expression level of ssc-miR-378 was low at 33 days post-coitus (dpc), increased at 65 and 90 dpc, peaked at postnatal day 0, and finally declined and maintained a comparatively stable level. This expression profile suggested that ssc-miR-378 was a new candidate miRNA for myogenesis and participated in skeletal muscle development in pigs. Target prediction and KEGG pathway analysis suggested that bone morphogenetic protein 2 (BMP2) and mitogen-activated protein kinase 1 (MAPK1), both of which were relevant to proliferation and differentiation, might be the potential targets of miR-378. Luciferase activities of report vectors containing the 3′UTR of porcine BMP2 or MAPK1 were downregulated by miR-378, which suggested that miR-378 probably regulated myogenesis though the regulation of these two genes.