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.

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F2 intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F2 intercross and an F19 Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F2 animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F2 and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related 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.

Carcass and meat quality traits are economically important in pigs. In this study, 17 carcass composition traits and 23 meat quality traits were recorded in 1028 F(2) animals from a White Duroc x Erhualian resource population. All pigs in this experimental population were genotyped for 194 informative markers covering the entire porcine genome. Seventy-seven genome-wide significant quantitative trait loci (QTL) for carcass traits and 68 for meat quality were mapped to 34 genomic regions. These results not only confirmed many previously reported QTL but also revealed novel regions associated with the measured traits. For carcass traits, the most prominent QTL was identified for carcass length and head weight at 57 cM on SSC7, which explained up to 50% of the phenotypic variance and had a 95% confidence interval of only 3 cM. Moreover, QTL for kidney and spleen weight and lengths of cervical vertebrae were reported for the first time in pigs. For meat quality traits, two significant QTL on SSC5 and X were identified for both intramuscular fat content and marbling score in the longissimus muscle, while three significant QTL on SSC1 and SSC9 were found exclusively for IMF. Both LM and the semimembranous muscle showed common QTL for colour score on SSC4, 5, 7, 8, 13 and X and discordant QTL on other chromosomes. White Duroc alleles at a majority of QTL detected were favourable for carcass composition, while favourable QTL alleles for meat quality originated from both White Duroc and Erhualian.

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.