Jideng Ma

Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars

We report the sequencing at 131× coverage, de novo assembly and analyses of the genome of a female Tibetan wild boar. We also resequenced the whole genomes of 30 Tibetan wild boars from six major distributed locations and 18 geographically related pigs in China. We characterized genetic diversity, population structure and patterns of evolution. We searched for genomic regions under selection, which includes genes that are involved in hypoxia, olfaction, energy metabolism and drug response. Comparing the genome of Tibetan wild boar with those of neighboring Chinese domestic pigs further showed the impact of thousands of years of artificial selection and different signatures of selection in wild boar and domestic pig. We also report genetic adaptations in Tibetan wild boar that are associated with high altitudes and characterize the genetic basis of increased salivation in domestic pig.

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.

Identification of differences in microRNA transcriptomes between porcine oxidative and glycolytic skeletal muscles.

BackgroundMicroRNAs (miRNAs) are a type of non-coding small RNA ~22 nucleotides in length that regulate the expression of protein coding genes at the post-transcriptional level. Glycolytic and oxidative myofibers, the two main types of skeletal muscles, play important roles in metabolic health as well as in meat quality and production in the pig industry. Previous expression profile studies of different skeletal muscle types have focused on these aspects of mRNA and proteins; nonetheless, an explanation of the miRNA transcriptome differences between these two distinct muscles types is long overdue.ResultsHerein, we present a comprehensive analysis of miRNA expression profiling between the porcine longissimus doris muscle (LDM) and psoas major muscle (PMM) using a deep sequencing approach. We generated a total of 16.62 M (LDM) and 18.46 M (PMM) counts, which produced 15.22 M and 17.52 M mappable sequences, respectively, and identified 114 conserved miRNAs and 89 novel miRNA*s. Of 668 unique miRNAs, 349 (52.25%) were co-expressed, of which 173 showed significant differences (P < 0.01) between the two muscle types. Muscle-specific miR-1-3p showed high expression levels in both libraries (LDM, 32.01%; PMM, 20.15%), and miRNAs that potentially affect metabolic pathways (such as the miR-133 and -23) showed significant differences between the two libraries, indicating that the two skeletal muscle types shared mainly muscle-specific miRNAs but expressed at distinct levels according to their metabolic needs. In addition, an analysis of the Gene Ontology (GO) terms and KEGG pathway associated with the predicted target genes of the differentially expressed miRNAs revealed that the target protein coding genes of highly expressed miRNAs are mainly involved in skeletal muscle structural development, regeneration, cell cycle progression, and the regulation of cell motility.ConclusionOur study indicates that miRNAs play essential roles in the phenotypic variations observed in different muscle fiber types.

Intrinsic Features in MicroRNA Transcriptomes Link Porcine Visceral Rather than Subcutaneous Adipose Tissues to Metabolic Risk

MicroRNAs (miRNAs) are non-coding small RNA ∼22 nucleotides in length that can regulate the expression of a wide range of coding genes at the post-transcriptional level. Visceral adipose tissues (VATs) and subcutaneous adipose tissues (SATs), the two main fat compartments in mammals, are anatomically, physiologically, metabolically, and clinically distinct. Various studies of adipose tissues have focused mainly on DNA methylation, and mRNA and protein expression, nonetheless little research sheds directly light on the miRNA transcriptome differences between these two distinct adipose tissue types. Here, we present a comprehensive investigation of miRNA transcriptomes across six variant porcine adipose tissues by small RNA-sequencing. We identified 219 known porcine miRNAs, 97 novel miRNA*s, and 124 miRNAs that are conserved to other mammals. A set of universally abundant miRNAs (i.e., miR-148a-3p, miR-143-3p, miR-27b-3p, miR-let-7a-1-5p, and miR-let-7f-5p) across the distinct adipose tissues was found. This set of miRNAs may play important housekeeping roles that are involved in adipogenesis. Clustering analysis indicated significant variations in miRNA expression between the VATs and SATs, and highlighted the role of the greater omentum in responding to potential metabolic risk because of the observed enrichment in this tissue of the immune- and inflammation-related miRNAs, such as the members of miR-17-92 cluster and miR-181 family. Differential expression of the miRNAs between the VATs and SATs, and miRNA target prediction analysis revealed that the VATs-specific enriched miRNAs were associated mainly with immune and inflammation responses. In summary, the differences of miRNA expression between the VATs and SATs revealed some of their intrinsic differences and indicated that the VATs might be closely associated with increased risk of metabolic disorders.

MicroRNA transcriptomes relate intermuscular adipose tissue to metabolic risk.

Intermuscular adipose tissue is located between the muscle fiber bundles in skeletal muscles, and has similar metabolic features to visceral adipose tissue, which has been found to be related to a number of obesity-related diseases. Although various miRNAs are known to play crucial roles in adipose deposition and adipogenesis, the microRNA transcriptome of intermuscular adipose tissue has not, until now, been studied. Here, we sequenced the miRNA transcriptomes of porcine intermuscular adipose tissue by small RNA-sequencing and compared it to a representative subcutaneous adipose tissue. We found that the inflammation- and diabetes-related miRNAs were significantly enriched in the intermuscular rather than in the subcutaneous adipose tissue. A functional enrichment analysis of the genes predicted to be targeted by the enriched miRNAs also indicated that intermuscular adipose tissue was associated mainly with immune and inflammation responses. Our results suggest that the intermuscular adipose tissue should be recognized as a potential metabolic risk factor of obesity.

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.

Mitochondrial DNA Evidence Indicates the Local Origin of Domestic Pigs in the Upstream Region of the Yangtze River

Previous studies have indicated two main domestic pig dispersal routes in East Asia: one is from the Mekong region, through the upstream region of the Yangtze River (URYZ) to the middle and upstream regions of the Yellow River, the other is from the middle and downstream regions of the Yangtze River to the downstream region of the Yellow River, and then to northeast China. The URYZ was regarded as a passageway of the former dispersal route; however, this assumption remains to be further investigated. We therefore analyzed the hypervariable segements of mitochondrial DNA from 513 individual pigs mainly from Sichuan and the Tibet highlands and 1,394 publicly available sequences from domestic pigs and wild boars across Asia. From the phylogenetic tree, most of the samples fell into a mixed group that was difficult to distinguish by breed or geography. The total network analysis showed that the URYZ pigs possessed a dominant position in haplogroup A and domestic pigs shared the same core haplotype with the local wild boars, suggesting that pigs in group A were most likely derived from the URYZ pool. In addition, a region-wise network analysis determined that URYZ contains 42 haplotypes of which 22 are unique indicating the high diversity in this region. In conclusion, our findings confirmed that pigs from the URYZ were domesticated in situ.

Transcriptional Profiling of Swine Lung Tissue after Experimental Infection with Actinobacillus pleuropneumoniae

Porcine pleuropneumonia is a highly contagious respiratory disease that causes great economic losses worldwide. In this study, we aimed to explore the underlying relationship between infection and injury by investigation of the whole porcine genome expression profiles of swine lung tissues post-inoculated with experimentally Actinobacillus pleuropneumoniae. Expression profiling experiments of the control group and the treatment group were conducted using a commercially available Agilent Porcine Genechip including 43,603 probe sets. Microarray analysis was conducted on profiles of lung from challenged versus non-challenged swine. We found 11,929 transcripts, identified as differentially expressed at the p ≤0.01 level. There were 1188 genes annotated as swine genes in the GenBank Data Base. GO term analysis identified a total of 89 biological process categories, 82 cellular components and 182 molecular functions that were significantly affected, and at least 27 biological process categories that were related to the host immune response. Gene set enrichment analysis identified 13 pathways that were significantly associated with host response. Many proinflammatory-inflammatory cytokines were activated and involved in the regulation of the host defense response at the site of inflammation; while the cytokines involved in regulation of the host immune response were suppressed. All changes of genes and pathways of induced or repressed expression not only led to a decrease in antigenic peptides presented to T lymphocytes by APCs via the MHC and alleviated immune response injury induced by infection, but also stimulated stem cells to produce granulocytes (neutrophils, eosinophils, and basophils) and monocyte, and promote neutrophils and macrophages to phagocytose bacterial and foreign antigen at the site of inflammation. The defense function of swine infection with Actinobacillus pleuropneumoniae was improved, while its immune function was decreased.

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.

Detecting mitochondrial signatures of selection in wild Tibetan pigs and domesticated pigs

Abstract Selection in genomic regions is prevalent in mammals; however, the effects of selection on the mitogenome are not clearly understood. We determined the complete mitochondrial DNA (mtDNA) sequences from six wild Tibetan pigs from the Tibetan plateau and four domestic pig breeds from the lowland of neighboring southwest China. Nucleotide diversity analysis using the sliding window method showed that the nucleotide diversity of wild Tibetan pigs in most regions of the mitogenome was higher than that of domestic pigs. The 12 s ribosomal RNA showed relatively lower nucleotide diversity in Tibetan pigs, suggesting purifying selection of these genes during high-altitude adaptation. More non-synonymous nucleotide substitutions in the ATP6 were found in wild Tibetan pigs, indicating adaptive selection in Tibetan pigs. The results suggested distinct impacts of natural selection and artificial selection upon the mitogenome, especially the mitochondrial signatures of adaptive evolution in wild Tibetan pigs under natural selection.