Qianzi Tang

Whole-genome sequencing of Berkshire (European native pig) provides insights into its origin and domestication

Domesticated organisms have experienced strong selective pressures directed at genes or genomic regions controlling traits of biological, agricultural or medical importance. The genome of native and domesticated pigs provide a unique opportunity for tracing the history of domestication and identifying signatures of artificial selection. Here we used whole-genome sequencing to explore the genetic relationships among the European native pig Berkshire and breeds that are distributed worldwide, and to identify genomic footprints left by selection during the domestication of Berkshire. Numerous nonsynonymous SNPs-containing genes fall into olfactory-related categories, which are part of a rapidly evolving superfamily in the mammalian genome. Phylogenetic analyses revealed a deep phylogenetic split between European and Asian pigs rather than between domestic and wild pigs. Admixture analysis exhibited higher portion of Chinese genetic material for the Berkshire pigs, which is consistent with the historical record regarding its origin. Selective sweep analyses revealed strong signatures of selection affecting genomic regions that harbor genes underlying economic traits such as disease resistance, pork yield, fertility, tameness and body length. These discoveries confirmed the history of origin of Berkshire pig by genome-wide analysis and illustrate how domestication has shaped the patterns of genetic variation.

Detection of dietetically absorbed maize-derived microRNAs in pigs

MicroRNAs are a class of small RNAs that are important in post-transcriptional gene regulation in animals and plants. These single-stranded molecules are widely distributed in organisms and influence fundamental biological processes. Interestingly, recent studies have reported that diet-derived plant miRNAs could regulate mammalian gene expression, and these studies have broadened our view of cross-kingdom communication. In the present study, we evaluated miRNA levels in cooked maize-containing chow diets, and found that plant miRNAs were resistant to the harsh cooking conditions to a certain extent. After feeding fresh maize to pigs (7 days), maize-derived miRNAs could be detected in porcine tissues and serum, and the authenticity of these plant miRNAs was confirmed by using oxidization reactions. Furthermore, in vivo and in vitro experiments demonstrated that dietary maize miRNAs could cross the gastrointestinal tract and enter the porcine bloodstream. In the porcine cells, we found that plant miRNAs are very likely to specifically target their endogenous porcine mRNAs and influence gene expression in a fashion similar to that of mammalian miRNAs. Our results indicate that maize-derived miRNAs can cross the gastrointestinal tract and present in pigs, and these exogenous miRNAs have the potential to regulate mammalian gene expression.

Comprehensive variation discovery and recovery of missing sequence in the pig genome using multiple de novo assemblies

Uncovering genetic variation through resequencing is limited by the fact that only sequences with similarity to the reference genome are examined. Reference genomes are often incomplete and cannot represent the full range of genetic diversity as a result of geographical divergence and independent demographic events. To more comprehensively characterize genetic variation of pigs (Sus scrofa), we generated de novo assemblies of nine geographically and phenotypically representative pigs from Eurasia. By comparing them to the reference pig assembly, we uncovered a substantial number of novel SNPs and structural variants, as well as 137.02-Mb sequences harboring 1737 protein-coding genes that were absent in the reference assembly, revealing variants left by selection. Our results illustrate the power of whole-genome de novo sequencing relative to resequencing and provide valuable genetic resources that enable effective use of pigs in both agricultural production and biomedical research.

Transcriptome Analysis of Liangshan Pig Muscle Development at the Growth Curve Inflection Point and Asymptotic Stages Using Digital Gene Expression Profiling

Animal growth curves can provide essential information for animal breeders to optimize feeding and management strategies. However, the genetic mechanism underlying the phenotypic differentiation between the inflection point and asymptotic stages of the growth curve is not well characterized. Here, we employed Liangshan pigs in stages of growth at the inflection point (under inflection point: UIP) and the two asymptotic stages (before the inflection point: BIP, after the inflection point: AIP) as models to survey global gene expression in the longissimus dorsi muscle using digital gene expression (DGE) tag profiling. We found Liangshan pigs reached maximum growth rate (UIP) at 163.6 days of age and a weight of 134.6 kg. The DGE libraries generated 117 million reads of 5.89 gigabases in length. 21,331, 20,996 and 20,139 expressed transcripts were identified BIP, UIP and AIP, respectively. Among them, we identified 757 differentially expressed genes (DEGs) between BIP and UIP, and 271 DEGs between AIP and UIP. An enrichment analysis of DEGs proved the immune system was strengthened in the AIP stage. Energy metabolism rate, global transcriptional activity and bone development intensity were highest UIP. Meat from Liangshan pigs had the highest intramuscular fat content and most favorable fatty acid composition in the AIP. Three hundred eighty (27.70%) specific expression genes were highly enriched in QTL regions for growth and meat quality traits. This study completed a comprehensive analysis of diverse genetic mechanisms underlying the inflection point and asymptotic stages of growth. Our findings will serve as an important resource in the understanding of animal growth and development in indigenous pig breeds.

Transcriptomic analysis provides insight into high-altitude acclimation in domestic goats

Domestic goats are distributed in a wide range of habitats and have acclimated to their local environmental conditions. To investigate the gene expression changes of goats that are induced by high altitude stress, we performed RNA-seq on 27 samples from the three hypoxia-sensitive tissues (heart, lung, and skeletal muscle) in three indigenous populations from distinct altitudes (600 m, 2000 m, and 3000 m). We generated 129Gb of high-quality sequencing data (~4Gb per sample) and catalogued the expression profiles of 12,421 annotated hircine genes in each sample. The analysis showed global similarities and differences of high-altitude transcriptomes among populations and tissues as well as revealed that the heart underwent the most high-altitude induced expression changes. We identified numerous differentially expressed genes that exhibited distinct expression patterns, and nonsynonymous single nucleotide variant-containing genes that were highly differentiated between the high- and low-altitude populations. These genes have known or potential roles in hypoxia response and were enriched in functional gene categories potentially responsible for high-altitude stress. Therefore, they are appealing candidates for further investigation of the gene expression and associated regulatory mechanisms related to high-altitude acclimation.

Exosomal microRNAs in giant panda (Ailuropoda melanoleuca) breast milk: potential maternal regulators for the development of newborn cubs

The physiological role of miRNAs is widely understood to include fine-tuning the post-transcriptional regulation of a wide array of biological processes. Extensive studies have indicated that exosomal miRNAs in the bodily fluids of various organisms can be transferred between living cells for the delivery of gene silencing signals. Here, we illustrated the expression characteristics of exosomal miRNAs in giant panda breast milk during distinct lactation periods and highlighted the enrichment of immune- and development-related endogenous miRNAs in colostral and mature giant panda milk. These miRNAs are stable, even under certain harsh conditions, via the protection of extracellular vesicles. These findings indicate that breast milk may facilitate the dietary intake of maternal miRNAs by infants for the regulation of postnatal development. We also detected exogenous plant miRNAs from the primary food source of the giant panda (bamboo) in the exosomes of giant panda breast milk that were associated with regulatory roles in basic metabolism and neuron development. This result suggested that dietary plant miRNAs are absorbed by host cells and subsequently secreted into bodily fluids as potential cross-kingdom regulators. In conclusion, exosomal miRNAs in giant panda breast milk may be crucial maternal regulators for the development of intrinsic ‘slink’ newborn cubs.

Genomic analysis reveals selection in Chinese native black pig

Identification of genomic signatures that help reveal mechanisms underlying desirable traits in domesticated pigs is of significant biological, agricultural and medical importance. To identify the genomic footprints left by selection during domestication of the Enshi black pig, a typical native and meat-lard breed in China, we generated about 72-fold coverage of the pig genome using pools of genomic DNA representing three different populations of Enshi black pigs from three different locations. Combining this data with the available whole genomes of 13 Chinese wild boars, we identified 417 protein-coding genes embedded in the selected regions of Enshi black pigs. These genes are mainly involved in developmental and metabolic processes, response to stimulus, and other biological processes. Signatures of selection were detected in genes involved in body size and immunity (RPS10 and VASN), lipid metabolism (GSK3), male fertility (INSL6) and developmental processes (TBX19). These findings provide a window into the potential genetic mechanism underlying development of desirable phenotypes in Enshi black pigs during domestication and subsequent artificial selection. Thus, our results illustrate how domestication has shaped patterns of genetic variation in Enshi black pigs and provide valuable genetic resources that enable effective use of pigs in agricultural production.

Deciphering the microRNA transcriptome of skeletal muscle during porcine development

MicroRNAs (miRNAs) play critical roles in many important biological processes, such as growth and development in mammals. Various studies of porcine muscle development have mainly focused on identifying miRNAs that are important for fetal and adult muscle development; however, little is known about the role of miRNAs in middle-aged muscle development. Here, we present a comprehensive investigation of miRNA transcriptomes across five porcine muscle development stages, including one prenatal and four postnatal stages. We identified 404 known porcine miRNAs, 118 novel miRNAs, and 101 miRNAs that are conserved in other mammals. A set of universally abundant miRNAs was found across the distinct muscle development stages. This set of miRNAs may play important housekeeping roles that are involved in myogenesis. A short time-series expression miner analysis indicated significant variations in miRNA expression across distinct muscle development stages. We also found enhanced differentiation- and morphogenesis-related miRNA levels in the embryonic stage; conversely, apoptosis-related miRNA levels increased relatively later in muscle development. These results provide integral insight into miRNA function throughout pig muscle development stages. Our findings will promote further development of the pig as a model organism for human age-related muscle disease research.

Detection of genetic diversity and selection at the coding region of the melanocortin receptor 1 (MC1R) gene in Tibetan pigs and Landrace pigs.

Domestication and subsequent selective pressures have produced a large variety of pig coat colors in different regions and breeds. The melanocortin 1 receptor (MC1R) gene plays a crucial role in determining coat color of mammals. Here, we investigated genetic diversity and selection at the coding region of the porcine melanocortin receptor 1 (MC1R) in Tibetan pigs and Landrace pigs. By contrast, genetic variability was much lower in Landrace pigs than in Tibetan pigs. Meanwhile, haplotype analysis showed that Tibetan pigs possessed shared haplotypes, suggesting a possibility of recent introgression event by way of crossbreeding with neighboring domestic pigs or shared ancestral polymorphism. Additionally, we detected positive selection at the MC1R in both Tibetan pigs and Landrace pigs through the dN/dS analysis. These findings suggested that novel phenotypic change (dark coat color) caused by novel mutations may help Tibetan pigs against intensive solar ultraviolet (UV) radiation and camouflage in wild environment, whereas white coat color in Landrace were intentionally selected by human after domestication. Furthermore, both the phylogenetic analysis and the network analysis provided clues that MC1R in Asian and European wild boars may have initially experienced different selective pressures, and MC1R alleles diversified in modern domesticated pigs.

The miRNA Transcriptome Directly Reflects the Physiological and Biochemical Differences between Red, White, and Intermediate Muscle Fiber Types

MicroRNAs (miRNAs) are small non-coding RNAs that can regulate their target genes at the post-transcriptional level. Skeletal muscle comprises different fiber types that can be broadly classified as red, intermediate, and white. Recently, a set of miRNAs was found expressed in a fiber type-specific manner in red and white fiber types. However, an in-depth analysis of the miRNA transcriptome differences between all three fiber types has not been undertaken. Herein, we collected 15 porcine skeletal muscles from different anatomical locations, which were then clearly divided into red, white, and intermediate fiber type based on the ratios of myosin heavy chain isoforms. We further illustrated that three muscles, which typically represented each muscle fiber type (i.e., red: peroneal longus (PL), intermediate: psoas major muscle (PMM), white: longissimus dorsi muscle (LDM)), have distinct metabolic patterns of mitochondrial and glycolytic enzyme levels. Furthermore, we constructed small RNA libraries for PL, PMM, and LDM using a deep sequencing approach. Results showed that the differentially expressed miRNAs were mainly enriched in PL and played a vital role in myogenesis and energy metabolism. Overall, this comprehensive analysis will contribute to a better understanding of the miRNA regulatory mechanism that achieves the phenotypic diversity of skeletal muscles.