Congying Chen

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.

Uncovering the composition of microbial community structure and metagenomics among three gut locations in pigs with distinct fatness

Uncovering the phylogenetic composition of microbial community and the potential functional capacity of microbiome in different gut locations is of great importance to pig production. Here we performed a comparative analysis of gut microbiota and metagenomics among jejunum, ileum and cecum in pigs with distinct fatness. 16S rRNA gene sequencing revealed dramatic differences of microbial composition, diversity and species abundance between small intestine and cecum. Clostridium and SMB53 were enriched in the small intestine, while Prevotella, Treponema, Ruminococcus and Faecalibacterium showed a higher abundance in the cecum. Functional capacity analysis of gut microbiome revealed that the microbiome of small intestine plays important roles in the metabolism of small molecule nutrients, while the microbiome of cecum has the stronger ability to degrade xylan, pectin and cellulose. We identified tens of fatness associated-bacterial species including Escherichia spp. that showed a notable increase of relative abundance in all three gut locations of high fatness pigs. We further suggested that the potential pathogens, inflammation process, and microbial metabolism and nutrient sensing are involved in the high fatness of pigs. These results improve our knowledge about microbiota compositions in different gut locations, and give an insight into the effect of gut microbiota on porcine fatness.

Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations

BackgroundFatty acid composition in muscle is an important factor that affects the nutritive value and taste of pork. To investigate the genetic architecture of fatty acid composition of pork, we measured fatty acid contents in longissimus dorsi muscle of 1244 pigs from three divergent populations and conducted genome-wide association studies (GWAS) for fatty acid contents.ResultsWe detected 26 genome-wide significant quantitative trait loci (QTL) on eight chromosomes (SSC for Sus scrofa) for eight fatty acids. These loci not only replicated previously reported QTL for C18:0 on SSC14 and C20:0 on SSC16, but also included several novel QTL such as those for C20:1 on SSC7, C14:0 on SSC9, and C14:0, C16:0 and C16:1 on SSC12. Furthermore, we performed a meta-analysis of GWAS on five populations, including the three populations that were investigated in this study and two additional populations that we had previously examined. This enhanced the strength of the associations detected between fatty acid composition and several marker loci, especially for those for C18:0 on SSC14 and C20:0 on SSC16. The genes ELOVL5, ELOVL6, ELOVL7, FASN, SCD and THRSP, which have functions that are directly relevant to fatty acid metabolism, are proximal to the top associated markers at six significant QTL.ConclusionsThe findings improve our understanding of the genetic architecture of fatty acid composition in pork and contribute to further fine-map and characterize genes that influence fatty acid composition.

Genome-wide association studies for fatty acid metabolic traits in five divergent pig populations

Fatty acid composition profiles are important indicators of meat quality and tasting flavor. Metabolic indices of fatty acids are more authentic to reflect meat nutrition and public acceptance. To investigate the genetic mechanism of fatty acid metabolic indices in pork, we conducted genome-wide association studies (GWAS) for 33 fatty acid metabolic traits in five pig populations. We identified a total of 865 single nucleotide polymorphisms (SNPs), corresponding to 11 genome-wide significant loci on nine chromosomes and 12 suggestive loci on nine chromosomes. Our findings not only confirmed seven previously reported QTL with stronger association strength, but also revealed four novel population-specific loci, showing that investigations on intermediate phenotypes like the metabolic traits of fatty acids can increase the statistical power of GWAS for end-point phenotypes. We proposed a list of candidate genes at the identified loci, including three novel genes (FADS2, SREBF1 and PLA2G7). Further, we constructed the functional networks involving these candidate genes and deduced the potential fatty acid metabolic pathway. These findings advance our understanding of the genetic basis of fatty acid composition in pigs. The results from European hybrid commercial pigs can be immediately transited into breeding practice for beneficial fatty acid composition.

Fine mapping of a QTL for ear size on porcine chromosome 5 and identification of high mobility group AT-hook 2 ( HMGA2 ) as a positional candidate gene

BackgroundEar size and shape are distinct conformation characteristics of pig breeds. Previously, we identified a significant quantitative trait locus (QTL) influencing ear surface on pig chromosome 5 in a White Duroc × Erhualian F2 resource population. This QTL explained more than 17% of the phenotypic variance.MethodsFour new markers on pig chromosome 5 were genotyped across this F2 population. RT-PCR was performed to obtain expression profiles of different candidate genes in ear tissue. Standard association test, marker-assisted association test and F-drop test were applied to determine the effects of single nucleotide polymorphisms (SNP) on ear size. Three synthetic commercial lines were also used for the association test.ResultsWe refined the QTL to an 8.7-cM interval and identified three positional candidate genes i.e. HMGA2, SOX5 and PTHLH that are expressed in ear tissue. Seven SNP within these three candidate genes were selected and genotyped in the F2 population. Of the seven SNP, HMGA2 SNP (JF748727: g.2836 A > G) showed the strongest association with ear size in the standard association test and marker-assisted association test. With the F-drop test, F value decreased by more than 97% only when the genotypes of HMGA2 g.2836 A > G were included as a fixed effect. Furthermore, the significant association between g.2836 A > G and ear size was also demonstrated in the synthetic commercial Sutai pig line. The haplotype-based association test showed that the phenotypic variance explained by HMGA2 was similar to that explained by the QTL and at a much higher level than by SOX5. More interestingly, HMGA2 is also located within the dog orthologous chromosome region, which has been shown to be associated with ear type and size.ConclusionsHMGA2 was the closest gene with a potential functional effect to the QTL or marker for ear size on chromosome 5. This study will contribute to identify the causative gene and mutation underlying this QTL.

Unraveling the Fecal Microbiota and Metagenomic Functional Capacity Associated with Feed Efficiency in Pigs

Gut microbiota plays fundamental roles in energy harvest, nutrient digestion, and intestinal health, especially in processing indigestible components of polysaccharides in diet. Unraveling the microbial taxa and functional capacity of gut microbiome associated with feed efficiency can provide important knowledge to improve pig feed efficiency in swine industry. In the current research, we studied the association of fecal microbiota with feed efficiency in 280 commercial Duroc pigs. All experimental pigs could be clustered into two enterotype-like groups. Different enterotypes showed the tendency of association with the feed efficiency (P = 0.07). We further identified 31 operational taxonomic units (OTUs) showing the potential associations with porcine feed efficiency. These OTUs were mainly annotated to the bacteria related to the metabolisms of dietary polysaccharides. Although we did not identify the RFI-associated bacterial species at FDR < 0.05 level, metagenomic sequencing analysis did find the distinct function capacities of gut microbiome between the high and low RFI pigs (FDR < 0.05). The KEGG orthologies related to nitrogen metabolism, amino acid metabolism, and transport system, and eight KEGG pathways including glycine, serine, and threonine metabolism were positively associated with porcine feed efficiency. We inferred that gut microbiota might improve porcine feed efficiency through promoting intestinal health by the SCFAs produced by fermenting dietary polysaccharides and improving the utilization of dietary protein. The present results provided important basic knowledge for improving porcine feed efficiency through modulating gut microbiome.

A meta analysis of genome-wide association studies for limb bone lengths in four pig populations

BackgroundLimb bone length is an economically important trait in pigs, because it is negatively correlated with backfat thickness, and is also a determinant to the yield of hip and loin. Moreover, abnormal growth of the limb bone leads to leg structural weakness. Until now, the genetic architecture of the pig lime bone length remains poorly understood. The object of this study was to map genomic loci for limb bone length by genome-wide association study (GWAS) on 4 pig populations.ResultsWe measured the lengths of five limb bones including scapula, humerus, ulna, femur and tibia that were dissected from the right-side carcass of 925, 331, 314 and 434 animals from White Duroc × Erhualian F2 intercross, Erhualian, Laiwu and Sutai populations, respectively. We genotyped the 2004 pigs for 62,163 single nucleotide polymorphisms (SNPs) on the Porcine SNP60 BeadChip, and performed GWAS and a GWAS meta analysis in the 4 populations. In total, we identified 12 and 4 loci associated with the limb bone lengths at suggestive and genome-wide significant levels respectively, of which 4 loci were reported for the first time. The most prominent locus was identified in a 924-kb (kilo base pairs) linkage disequilibrium block on Sus Scrofa chromosome (SSC) 7, and High Mobility Group AT-hook 1 (HMGA1) appears to be a strong candidate gene in this region. Another promising locus is located in the middle of SSC4, and Pleiomorphic Adenoma Gene 1 (PLAG1) is a functionally plausible candidate gene underlying the locus. Because the lengths of the 5 limb bones are highly correlated to each other, most of significant loci were associated with all of the 5 traits; however, several loci showed specific effect on the length of one limb bone, such as the locus at the proximal end of SSC2 associated with only the scapula length.ConclusionTo our knowledge, this study was the first GWAS meta analysis for limb bone lengths in pigs. As expected, the meta analysis is more powerful to identify genomic loci. A total of 16 loci were identified in this study, including four novel loci. HMGA1 and PLAG1 are two appearing candidate genes for pig limb bone lengths, which warrant further investigations.

A Genome-Wide Investigation of Expression Characteristics of Natural Antisense Transcripts in Liver and Muscle Samples of Pigs

Natural antisense transcripts are endogenous transcripts that are complementary to the sense-strand of DNA. These transcripts have been identified in various eukaryotic species and are involved in a broad range of regulatory events and biological processes. However, their general biological functions, expression characteristics and regulatory mechanisms are still unclear. In this study, 497 liver and 586 muscle samples were harvested from a White Duroc×Erhualian F2 resource population. The expression profiles of sense and antisense transcripts were determined by tag-based RNA sequencing. We identified 33.7% and 20.4% of transcripts having both sense and antisense expression, and 12.5% and 6.1% of transcripts only expressing antisense transcripts in liver and muscle, respectively. More than 32.2% of imprinting or predicted imprinting genes in the geneimprint database were detected with both sense and antisense expression. The correlations between sense and antisense expression in sense-antisense pairs were diverse in both liver and muscle, showing positive, negative or absent correlation. Antisense expression increases gene expression variability. More interestingly, compared to eQTL mapping of sense transcripts in which more than one eQTL was mapped for a transcript, only one eQTL was identified for each antisense transcript, and the percentage of cis-eQTL in antisense eQTL was higher than that in sense eQTL. This suggests that the expressions of antisense transcripts tend to be cis-regulated by a single genomic locus. To our knowledge, this study is the first systematical investigation of antisense transcription in pigs. The findings improve our understanding of the complexity of porcine transcriptome.

Evaluating the Contribution of Gut Microbiota to the Variation of Porcine Fatness with the Cecum and Fecal Samples

Microbial community in gastrointestinal tract participates in the development of the obesity as well as quite a few metabolic diseases in human. However, there are few studies about the relationship between gut microbiota and porcine fatness. Here, we used high-throughput sequencing to perform 16S rRNA gene analysis in 256 cecum luminal samples from Erhualian pigs and 244 stools from Bamaxiang pigs, and adopted a two-part model statistical method to evaluate the association of gut microbes with porcine fatness. As the results, we identified a total of 6 and 108 operational taxonomic units (OTUs), and 9 and 10 bacterial taxa which showed significant associations with fatness traits in the stool and cecum samples, respectively. Cross-validation analysis indicated that gut microbiome showed the largest effect on abdominal adipose by explaining 2.73% phenotypic variance of abdominal fat weight. Significantly more fatness-associated OTUs were identified in the cecum samples than that in the stools, suggesting that cecum luminal samples were better used for identification of fatness-associated microbes than stools. The fatness-associated OTUs were mainly annotated to Lachnospiraceae, Ruminococcaceae, Prevotella, Treponema, and Bacteroides. These microbes have been reported to produce short-chain fatty acids by fermenting dietary indigested polysaccharide and pectin. The short-chain fatty acids can regulate host body energy homeostasis, protect host from inflammation and inhibit fat mass development. Our findings suggested that the gut microbiome may be an important factor modulating fatness in pigs.

Genome-Wide Association Analysis for Blood Lipid Traits Measured in Three Pig Populations Reveals a Substantial Level of Genetic Heterogeneity

Serum lipids are associated with myocardial infarction and cardiovascular disease in humans. Here we dissected the genetic architecture of blood lipid traits by applying genome-wide association studies (GWAS) in 1,256 pigs from Laiwu, Erhualian and Duroc × (Landrace × Yorkshire) populations, and a meta-analysis of GWAS in more than 2,400 pigs from five diverse populations. A total of 22 genomic loci surpassing the suggestive significance level were detected on 11 pig chromosomes (SSC) for six blood lipid traits. Meta-analysis of GWAS identified 5 novel loci associated with blood lipid traits. Comparison of GWAS loci across the tested populations revealed a substantial level of genetic heterogeneity for porcine blood lipid levels. We further evaluated the causality of nine polymorphisms nearby or within the APOB gene on SSC3 for serum LDL-C and TC levels. Of the 9 polymorphisms, an indel showed the most significant association with LDL-C and TC in Laiwu pigs. But the significant association was not identified in the White Duroc × Erhualian F2 resource population, in which the QTL for LDL-C and TC was also detected on SSC3. This indicates that population-specific signals may exist for the SSC3 QTL. Further investigations are warranted to validate this assumption.