Huashui Ai

Adaptation and possible ancient interspecies introgression in pigs identified by whole-genome sequencing.

Domestic pigs have evolved genetic adaptations to their local environmental conditions, such as cold and hot climates. We sequenced the genomes of 69 pigs from 15 geographically divergent locations in China and detected 41 million variants, of which 21 million were absent from the dbSNP database. In a genome-wide scan, we identified a set of loci that likely have a role in regional adaptations to high- and low-latitude environments within China. Intriguingly, we found an exceptionally large (14-Mb) region with a low recombination rate on the X chromosome that appears to have two distinct haplotypes in the high- and low-latitude populations, possibly underlying their adaptation to cold and hot environments, respectively. Surprisingly, the adaptive sweep in the high-latitude regions has acted on DNA that might have been introgressed from an extinct Sus species. Our findings provide new insights into the evolutionary history of pigs and the role of introgression in adaptation.

Genetic Diversity, Linkage Disequilibrium and Selection Signatures in Chinese and Western Pigs Revealed by Genome-Wide SNP Markers

To investigate population structure, linkage disequilibrium (LD) pattern and selection signature at the genome level in Chinese and Western pigs, we genotyped 304 unrelated animals from 18 diverse populations using porcine 60 K SNP chips. We confirmed the divergent evolution between Chinese and Western pigs and showed distinct topological structures of the tested populations. We acquired the evidence for the introgression of Western pigs into two Chinese pig breeds. Analysis of runs of homozygosity revealed that historical inbreeding reduced genetic variability in several Chinese breeds. We found that intrapopulation LD extents are roughly comparable between Chinese and Western pigs. However, interpopulation LD is much longer in Western pigs compared with Chinese pigs with average r2 0.3 values of 125 kb for Western pigs and only 10.5 kb for Chinese pigs. The finding indicates that higher-density markers are required to capture LD with causal variants in genome-wide association studies and genomic selection on Chinese pigs. Further, we looked across the genome to identify candidate loci under selection using FST outlier tests on two contrast samples: Tibetan pigs versus lowland pigs and belted pigs against non-belted pigs. Interestingly, we highlighted several genes including ADAMTS12, SIM1 and NOS1 that show signatures of natural selection in Tibetan pigs and are likely important for genetic adaptation to high altitude. Comparison of our findings with previous reports indicates that the underlying genetic basis for high-altitude adaptation in Tibetan pigs, Tibetan peoples and yaks is likely distinct from one another. Moreover, we identified the strongest signal of directional selection at the EDNRB loci in Chinese belted pigs, supporting EDNRB as a promising candidate gene for the white belt coat color in Chinese pigs. Altogether, our findings advance the understanding of the genome biology of Chinese and Western pigs.

A comprehensive survey of copy number variation in 18 diverse pig populations and identification of candidate copy number variable genes associated with complex traits.

BackgroundCopy number variation (CNV) is a major source of structural variants and has been commonly identified in mammalian genome. It is associated with gene expression and may present a major genetic component of phenotypic diversity. Unlike many other mammalian genomes where CNVs have been well annotated, studies of porcine CNV in diverse breeds are still limited.ResultHere we used Porcine SNP60 BeadChip and PennCNV algorithm to identify 1,315 putative CNVs belonging to 565 CNV regions (CNVRs) in 1,693 pigs from 18 diverse populations. Total 538 out of 683 CNVs identified in a White Duroc × Erhualian F2 population fit Mendelian transmission and 6 out of 7 randomly selected CNVRs were confirmed by quantitative real time PCR. CNVRs were non-randomly distributed in the pig genome. Several CNV hotspots were found on pig chromosomes 6, 11, 13, 14 and 17. CNV numbers differ greatly among different pig populations. The Duroc pigs were identified to have the most number of CNVs per individual. Among 1,765 transcripts located within the CNVRs, 634 genes have been reported to be copy number variable genes in the human genome. By integrating analysis of QTL mapping, CNVRs and the description of phenotypes in knockout mice, we identified 7 copy number variable genes as candidate genes for phenotypes related to carcass length, backfat thickness, abdominal fat weight, length of scapular, intermuscle fat content of logissimus muscle, body weight at 240 day, glycolytic potential of logissimus muscle, mean corpuscular hemoglobin, mean corpuscular volume and humerus diameter.ConclusionWe revealed the distribution of the unprecedented number of 565 CNVRs in pig genome and investigated copy number variable genes as the possible candidate genes for phenotypic traits. These findings give novel insights into porcine CNVs and provide resources to facilitate the identification of trait-related CNVs.

A Missense Mutation in PPARD Causes a Major QTL Effect on Ear Size in Pigs

Chinese Erhualian is the most prolific pig breed in the world. The breed exhibits exceptionally large and floppy ears. To identify genes underlying this typical feature, we previously performed a genome scan in a large scale White Duroc × Erhualian cross and mapped a major QTL for ear size to a 2-cM region on chromosome 7. We herein performed an identical-by-descent analysis that defined the QTL within a 750-kb region. Historically, the large-ear feature has been selected for the ancient sacrificial culture in Erhualian pigs. By using a selective sweep analysis, we then refined the critical region to a 630-kb interval containing 9 annotated genes. Four of the 9 genes are expressed in ear tissues of piglets. Of the 4 genes, PPARD stood out as the strongest candidate gene for its established role in skin homeostasis, cartilage development, and fat metabolism. No differential expression of PPARD was found in ear tissues at different growth stages between large-eared Erhualian and small-eared Duroc pigs. We further screened coding sequence variants in the PPARD gene and identified only one missense mutation (G32E) in a conserved functionally important domain. The protein-altering mutation showed perfect concordance (100%) with the QTL genotypes of all 19 founder animals segregating in the White Duroc × Erhualian cross and occurred at high frequencies exclusively in Chinese large-eared breeds. Moreover, the mutation is of functional significance; it mediates down-regulation of β-catenin and its target gene expression that is crucial for fat deposition in skin. Furthermore, the mutation was significantly associated with ear size across the experimental cross and diverse outbred populations. A worldwide survey of haplotype diversity revealed that the mutation event is of Chinese origin, likely after domestication. Taken together, we provide evidence that PPARD G32E is the variation underlying this major QTL.

Genome-Wide Association Study of Meat Quality Traits in a White Duroc×Erhualian F2 Intercross and Chinese Sutai Pigs

Thousands of QTLs for meat quality traits have been identified by linkage mapping studies, but most of them lack precise position or replication between populations, which hinder their application in pig breeding programs. To localize QTLs for meat quality traits to precise genomic regions, we performed a genome-wide association (GWA) study using the Illumina PorcineSNP60K Beadchip in two swine populations: 434 Sutai pigs and 933 F2 pigs from a White Duroc×Erhualian intercross. Meat quality traits, including pH, color, drip loss, moisture content, protein content and intramuscular fat content (IMF), marbling and firmness scores in the M. longissimus (LM) and M. semimembranosus (SM) muscles, were recorded on the two populations. In total, 127 chromosome-wide significant SNPs for these traits were identified. Among them, 11 SNPs reached genome-wise significance level, including 1 on SSC3 for pH, 1 on SSC3 and 3 on SSC15 for drip loss, 3 (unmapped) for color a*, and 2 for IMF each on SSC9 and SSCX. Except for 11 unmapped SNPs, 116 significant SNPs fell into 28 genomic regions of approximately 10 Mb or less. Most of these regions corresponded to previously reported QTL regions and spanned smaller intervals than before. The loci on SSC3 and SSC7 appeared to have pleiotropic effects on several related traits. Besides them, a few QTL signals were replicated between the two populations. Further, we identified thirteen new candidate genes for IMF, marbling and firmness, on the basis of their positions, functional annotations and reported expression patterns. The findings will contribute to further identification of the causal mutation underlying these QTLs and future marker-assisted selection in pigs.

Genome-wide QTL mapping for three traits related to teat number in a White Duroc × Erhualian pig resource population

BackgroundTeat number is an important fertility trait for pig production, reflecting the mothering ability of sows. It is also a discrete and often canalized trait presenting bilateral symmetry with minor differences between the two sides, providing a potential power to evaluate fluctuating asymmetry and developmental instability. The knowledge of its genetic control is still limited. In this study, a genome-wide scan was performed with 183 microsatellites covering the pig genome to identify quantitative trait loci (QTL) for three traits related to teat number including the total teat number (TTN), the teat number at the left (LTN) and right (RTN) sides in a large scale White Duroc × Erhualian resource population.ResultsA sex-average linkage map with a total length of 2350.3 cM and an average marker interval of 12.84 cM was constructed. Eleven genome-wide significant QTL for TTN were detected on 8 autosomes including pig chromosomes (SSC) 1, 3, 4, 5, 6, 7, 8 and 12. Six suggestive QTL for this trait were detected on SSC6, 9, 13, 14 and 16. Eight chromosomal regions each on SSC1, 3, 4, 5, 6, 7, 8 and 12 showed significant associations with LTN. These regions were also evidenced as significant QTL for RTN except for those on SSC6 and SSC8. The most significant QTL for the 3 traits were all located on SSC7. Erhualian alleles at most of the identified QTL had positive additive effects except for three QTL on SSC1 and SSC7, at which White Duroc alleles increased teat numbers. On SSC1, 6, 9, 13 and 16, significant dominance effects were observed on TTN, and predominant imprinting effect on TTN was only detected on SSC12.ConclusionThe results not only confirmed the QTL regions from previous experiments, but also identified five new QTL for the total teat number in swine. Minor differences between the QTL regions responsible for LTN and RTN were validated. Further fine mapping should be focused on consistently identified regions with small confidence intervals, such as those on SSC1, SSC7 and SSC12.

Population history and genomic signatures for high-altitude adaptation in Tibetan pigs

BackgroundThe Tibetan pig is one of domestic animals indigenous to the Qinghai-Tibet Plateau. Several geographically isolated pig populations are distributed throughout the Plateau. It remained an open question if these populations have experienced different demographic histories and have evolved independent adaptive loci for the harsh environment of the Plateau. To address these questions, we herein investigated ~ 40,000 genetic variants across the pig genome in a broad panel of 678 individuals from 5 Tibetan geographic populations and 34 lowland breeds.ResultsUsing a series of population genetic analyses, we show that Tibetan pig populations have marked genetic differentiations. Tibetan pigs appear to be 3 independent populations corresponding to the Tibetan, Gansu and Sichuan & Yunnan locations. Each population is more genetically similar to its geographic neighbors than to any of the other Tibetan populations. By applying a locus-specific branch length test, we identified both population-specific and -shared candidate genes under selection in Tibetan pigs. These genes, such as PLA2G12A, RGCC, C9ORF3, GRIN2B, GRID1 and EPAS1, are involved in high-altitude physiology including angiogenesis, pulmonary hypertension, oxygen intake, defense response and erythropoiesis. A majority of these genes have not been implicated in previous studies of highlanders and high-altitude animals.ConclusionTibetan pig populations have experienced substantial genetic differentiation. Historically, Tibetan pigs likely had admixture with neighboring lowland breeds. During the long history of colonization in the Plateau, Tibetan pigs have developed a complex biological adaptation mechanism that could be different from that of Tibetans and other animals. Different Tibetan pig populations appear to have both distinct and convergent adaptive loci for the harsh environment of the Plateau.

A Further Look at Porcine Chromosome 7 Reveals VRTN Variants Associated with Vertebral Number in Chinese and Western Pigs

The number of vertebrae is an economically important trait that affects carcass length and meat production in pigs. A major quantitative trait locus (QTL) for thoracic vertebral number has been repeatedly identified on pig chromosome (SSC) 7. To dissect the genetic basis of the major locus, we herein genotyped a large sample of animals from 3 experimental populations of Chinese and Western origins using 60K DNA chips. Genome-wide association studies consistently identified the locus across the 3 populations and mapped the locus to a 947-Kb region on SSC7. An identical-by-descent sharing assay refined the locus to a 100-Kb segment that harbors only two genes including VRTN and SYNDIG1L. Of them, VRNT has been proposed as a strong candidate of the major locus in Western modern breeds. Further, we resequenced the VRTN gene using DNA samples of 35 parental animals with known QTL genotypes by progeny testing. Concordance tests revealed 4 candidate causal variants as their genotypes showed the perfect segregation with QTL genotypes of the tested animals. An integrative analysis of evolutional constraints and functional elements supported two VRTN variants in a complete linkage disequilibrium phase as the most likely causal mutations. The promising variants significantly affect the number of thoracic vertebrae (one vertebra) in large scale outbred animals, and are segregating at rather high frequencies in Western pigs and at relatively low frequencies in a number of Chinese breeds. Altogether, we show that VRTN variants are significantly associated with the number of thoracic vertebrae in both Chinese and Western pigs. The finding advances our understanding of the genetic architecture of the vertebral number in pigs. Furthermore, our finding is of economical importance as it provides a robust breeding tool for the improvement of vertebral number and meat production in both Chinese indigenous pigs and Western present-day commercial pigs.

Susceptibility towards Enterotoxigenic Escherichia coli F4ac Diarrhea Is Governed by the MUC13 Gene in Pigs

Enterotoxigenic Escherichia coli (ETEC) F4ac is a major determinant of diarrhea and mortality in neonatal and young pigs. Susceptibility to ETEC F4ac is governed by the intestinal receptor specific for the bacterium and is inherited as a monogenic dominant trait. To identify the receptor gene (F4acR), we first mapped the locus to a 7.8-cM region on pig chromosome 13 using a genome scan with 194 microsatellite markers. A further scan with high density markers on chromosome 13 refined the locus to a 5.7-cM interval. Recombination breakpoint analysis defined the locus within a 2.3-Mb region. Further genome-wide mapping using 39,720 informative SNPs revealed that the most significant markers were proximal to the MUC13 gene in the 2.3-Mb region. Association studies in a collection of diverse outbred populations strongly supported that MUC13 is the most likely responsible gene. We characterized the porcine MUC13 gene that encodes two transcripts: MUC13A and MUC13B. Both transcripts have the characteristic PTS regions of mucins that are enriched in distinct tandem repeats. MUC13B is predicated to be heavily O-glycosylated, forming the binding site of the bacterium; while MUC13A does not have the O-glycosylation binding site. Concordantly, 127 independent pigs homozygous for MUC13A across diverse breeds are all resistant to ETEC F4ac, and all 718 susceptible animals from the broad breed panel carry at least one MUC13B allele. Altogether, we conclude that susceptibility towards ETEC F4ac is governed by the MUC13 gene in pigs. The finding has an immediate translation into breeding practice, as it allows us to establish an efficient and accurate diagnostic test for selecting against susceptible animals. Moreover, the finding improves our understanding of mucins that play crucial roles in defense against enteric pathogens. It revealed, for the first time, the direct interaction between MUC13 and enteric bacteria, which is poorly understood in mammals.

Detection of quantitative trait loci for growth- and fatness-related traits in a large-scale White Duroc × Erhualian intercross pig population

Growth and fatness are economically important traits in pigs. In this study, a genome scan was performed to detect quantitative trait loci (QTL) for 14 growth and fatness traits related to body weight, backfat thickness and fat weight in a large-scale White Duroc × Erhualian F(2) intercross. A total of 76 genome-wide significant QTL were mapped to 16 chromosomes. The most significant QTL was found on pig chromosome (SSC) 7 for fatness with unexpectedly small confidence intervals of ∼2 cM, providing an excellent starting point to identify causal variants. Common QTL for both fatness and growth traits were found on SSC4, 5, 7 and 8, and shared QTL for fat deposition were detected on SSC1, 2 and X. Time-series analysis of QTL for body weight at six growth stages revealed the continuously significant effects of the QTL on SSC4 at the fattening period and the temporal-specific expression of the QTL on SSC7 at the foetus and fattening stages. For fatness traits, Chinese Erhualian alleles were associated with increased fat deposition except that at the major QTL on SSC7. For growth traits, most of White Duroc alleles enhanced growth rates except for those at three significant QTL on SSC6, 7 and 9. The results confirmed many previously reported QTL and also detected novel QTL, revealing the complexity of the genetic basis of growth and fatness in pigs.