An imputation-based genome-wide association study for growth and fatness traits in Sujiang pigs

Sujiang pigs are a synthetic breed derived from Jiangquhai, Fengjing, and Duroc pigs. In this study, we sequenced the genome of 62 pigs with a coverage depth of 10× to 20×, including 27 Sujiang and 35 founder breed pigs, and we collected 360 global pigs’ genome sequence data from public databases including 39 Duroc pigs. We obtained a high-quality variant dataset of 365 Sujiang pigs by imputing the porcine 80 K single nucleotide polymorphism (SNP) Beadchip to the whole-genome scale with a total of 422 pigs as a reference panel. A dataset of 365 imputated Sujiang pigs was used to perform single-trait genome-wide association study (GWAS) and meta-analyses for growth and fatness traits. Single-trait GWAS identified 1 907, 18, and 14 SNPs surpassing the suggestively significant threshold for backfat thickness, chest circumference, and chest width, respectively. Meta-analyses identified 2 400 genome-wide significant SNPs and 520 suggestively significant SNPs for backfat thickness and chest circumference, and 719 genome-wide significant SNPs and 1 225 suggestively significant SNPs for all seven traits. According to the meta-analysis of backfat thickness and chest circumference, a remarkable region of 2.69 Mb on Sus scrofa chromosome 4 containing FAM110B, IMPAD1, LYN, MOS, PENK, PLAG1, SDR16C5 and XKR4 was identified as a candidate region. The haplotype heat map of the 2.69 Mb region verified that Sujiang pigs were derived from Duroc and Chinese indigenous pigs, especially Jiangquhai pigs. The Kruskal-Wallis test showed that haplotypes of the 2.69 Mb region significantly affected backfat thickness and chest circumference traits. We then focused on PLAG1, an important growth-related gene, and identified two synonymous SNPs with obvious differences among different breeds in the PLAG1 gene. We then performed genotyping of 365 Sujiang, 150 Duroc, 95 Jiangquhai, and 100 Fengjing pigs to confirm the above result and verified that the two variants significantly affected phenotypes of growth and fatness traits. Our findings not only provide insights into the genetic architecture of porcine growth and fatness traits but also provide potential markers for selective breeding of these traits in Sujiang pigs.     

Genome‐wide detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.     

A genome-wide association study identifies genomic loci associated with backfat thickness,carcass weight,and body weight in two commercial pig populations

Growth and fatness traits are economically important in the pig industry. To dissect the genetic architecture of these traits in commercial pigs, we conducted a genome-wide association study (GWAS) for carcass weight, backfat thickness, and body weight in two commercial populations: Duroc × (Landrace × Yorkshire) (DLY) and Duroc populations. To enhance the detection power, three GWAS approaches including single-trait GWAS, multi-trait GWAS and meta-analysis were used in this study. A total of 13 suggestive loci were identified on nine chromosomes. The most significant locus was detected at 272.05 Mb on SSC1, and it was associated with backfat thickness at the first rib in the DLY population. Three genes at the identified loci (TBC1D1, BAAT and PHLPP1) were highlighted as functionally plausible candidate genes for pig growth and fatness traits. Genome-wide significant locus was not evidenced in this study, indicating that large populations are required to identify QTL with minor effects on growth and fatness traits in commercial pig populations, in which intensively artificial selections have been imposed on these traits and small genetical variances usually retain in these traits.     

Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations

Background

Recently, genome-wide association studies (GWAS) have been reported on various pig traits. We performed a GWAS to analyze 22 traits related to growth and fatness on two pig populations: a White Duroc × Erhualian F₂ intercross population and a Chinese Sutai half-sib population.

Results

We identified 14 and 39 loci that displayed significant associations with growth and fatness traits at the genome-wide level and chromosome-wide level, respectively. The strongest association was between a 750 kb region on SSC7 (SSC for Sus scrofa) and backfat thickness at the first rib. This region had pleiotropic effects on both fatness and growth traits in F₂ animals and contained a promising candidate gene HMGA1 (high mobility group AT-hook 1). Unexpectedly, population genetic analysis revealed that the allele at this locus that reduces fatness and increases growth is derived from Chinese indigenous pigs and segregates in multiple Chinese breeds. The second strongest association was between the region around 82.85 Mb on SSC4 and average backfat thickness. PLAG1 (pleiomorphic adenoma gene 1), a gene under strong selection in European domestic pigs, is proximal to the top SNP and stands out as a strong candidate gene. On SSC2, a locus that significantly affects fatness traits mapped to the region around the IGF2 (insulin-like growth factor 2) gene but its non-imprinting inheritance excluded IGF2 as a candidate gene. A significant locus was also detected within a recombination cold spot that spans more than 30 Mb on SSCX, which hampered the identification of plausible candidate genes. Notably, no genome-wide significant locus was shared by the two experimental populations; different loci were observed that had both constant and time-specific effects on growth traits at different stages, which illustrates the complex genetic architecture of these traits.

Conclusions

We confirm several previously reported QTL and provide a list of novel loci for porcine growth and fatness traits in two experimental populations with Chinese Taihu and Western pigs as common founders. We showed that distinct loci exist for these traits in the two populations and identified HMGA1 and PLAG1 as strong candidate genes on SSC7 and SSC4, respectively.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0089-5) contains supplementary material, which is available to authorized users.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

丙酸通路基因多态性与猪肉质及胴体性状的关联分析

为了挖掘新的猪肉品质及胴体性状的候选基因,揭示猪肉质及胴体性状的遗传机制,文章将丙酸代谢通路作为候选通路,将通路内基因与猪肉质及胴体性状进行关联分析。实验采用37头三元杂交商品猪作为研究对象,首次针对丙酸通路中7个基因的36个SNP位点利用SNaPshot方法进行基因分型,分别用最小二乘模型及MB-MDR模型与肉质及胴体性状进行关联分析。结果发现,基因PCCB、MUT、MCEE及ACSS2上的4个SNP位点分别与肌内脂肪含量、背膘厚等性状显著相关(P<0.05),ACSS2与猪脂肪含量显著相关;MCEE及MUT与猪的背膘厚显著相关;PCCB基因与脂重显著相关。通过MB-MDR方法检测到多个SNP位点具有互作效应,并与背膘厚、水分含量、脂肪含量显著相关(P<0.05)。另外,丙酸代谢通路中的基因间的互作效应对猪肉品质有显著影响。     

Detection of novel SNPs and mapping of the fatness QTL on pig chromosome 7q1.1-1.4 region

Many QTLs for fatness traits have been mapped on pig chromosome 7q1.1-1.4 in various pig resource populations. Eight novel markers, including seven SNPs and one insertion or deletion within BTNL1, COL21A1, PPARD, GLP1R, MDFI, GNMT, ABCC10, and PLA2G7 genes, as well as two previously reported SNPs in SLC39A7 and HMGA1 genes, were genotyped in Large White and Meishan pig breeds. Except for two SNPs in HMGA1 and ABCC10 genes, allele frequencies of the other eight markers are highly significant different between Chinese indigenous Meishan breeds and Large White pig breeds. Eight polymorphic sites were then used for linkage and QTL mapping to refine the fatness QTL in a Large White × Meishan F(2) resource population. Five chromosome-wise significant QTLs were detected, of which the QTLs for leaf fat weight, backfat thickness at 6-7th rib and rump, and mean backfat thickness were narrowed to the interval between PPARD and GLP1R genes and the QTL for backfat thickness at thorax-waist between GNMT and PLA2G7 genes on SSC7p1.1-q1.4.     

Imprinting analysis of porcine <Emphasis Type="Italic">DIO3</Emphasis> gene in two fetal stages and association analysis with carcass and meat quality traits

Imprinted genes play important roles in mammalian growth, development and behavior. In this study, we obtained 1568 bp mRNA sequence of porcine DIO3 (deiodinase, iodothyronine, type III), and also identified its imprinting status during porcine fetal development. The complete open reading frame (ORF) encoding 278 amino acids. The porcine DIO3 mRNA was expressed predominantly in backfat, mildly in liver, uterus, kidney, heart, small intestine, muscle and stomach, and almost absent in spleen and lung. A single nucleotide polymorphism in exon (A/C ⁶⁸⁷) was used to investigate the allele frequencies in different pig breeds and the imprinting status in porcine embryonic tissues. The results indicate that DIO3 was imprinted in all the tested tissues. Statistical analysis showed the DIO3 gene polymorphism was significantly associated with almost all the fat deposition and carcass traits, including lean meat percentage (LMP), fat meat percentage (FMP), ratio of lean to fat (RLF), shoulder fat thickness (SFT), sixth–seventh rib fat thickness (RFT), buttock fat thickness (BFT), loin eye area (LEA), and intramuscular fat (IMF).     

Genomic regions affecting backfat thickness and cannon bone circumference identified by genome‐wide association study in a Duroc pig population

We performed a genome‐wide association study using the porcine 60K SNP array to detect QTL regions for nine traits in a three‐generational Duroc samples (n = 651), viz. generations 1, 2 and 3 from a population selected over five generations using a closed nucleus breeding scheme. We applied a linear mixed model for association mapping to detect SNP effects, adjusting for fixed effects (sex and season) and random polygenic effects (reflecting genetic relatedness), and derived a likelihood ratio statistic for each SNP using the efficient mixed‐model association method. We detected a region on SSC6 for backfat thickness (BFT) and on SSC7 for cannon bone circumference (CANNON), with a genome‐wide significance of P < 0.01 after Bonferroni correction. These regions had been detected previously in other pig populations. Six genes are located in the BFT‐associated region, while the CANNON‐associated region includes 66 genes. In the future, significantly associated SNPs, derived by sequencing the coding regions of the six genes in the BFT region, can be used in marker‐assisted selection of BFT, whereas haplotypes constructed from the SSC7 region with strong LD can be used to select for the CANNON trait in our resource family.     

Effects of genetic variability of the dairy goat growth hormone releasing hormone receptor (<Emphasis Type="Italic">GHRHR</Emphasis>) gene on growth traits

Growth hormone-releasing hormone receptor (GHRHR) plays a critical role in growth hormone (GH) synthesis, release and regulation of pituitary somatotroph expansion in vertebrates. The objective of this study was to investigate variations in goat GHRHR gene and their associations with growth traits in 668 dairy goats. The results showed four novel single nucleotide polymorphisms (SNPs): NC_007302:g.5203C>T, 7307C>G, 9583G>A and 9668A>C. In detail, the novel SNP C>T in the 5203rd nucleotide identified a missense mutation: CCC (Pro)>TCC (Phe) at position 116aa of the goat GHRHR (423aa). Besides, 9583G>A and 9668A>C polymorphism were in complete linkage disequilibrium. The genetic diversity analysis revealed that the Guanzhong dairy goat possessed intermediate genetic diversity in P3 and P7 loci, and the Xinong Sannen dairy goat belonged to poor genetic diversity in P4 locus. Significant associations between the genotypes of P3 locus and body length, body height and chest circumference was observed in Guanzhong goat (P < 0.05). However, in Xinong saanen population, significant statistical difference was only found in body height and body length (P < 0.05). In P4 and P7 loci, no significant associations were detected between any variant sites and body length, body height and chest circumference, as well as for the milk traits (P > 0.05). These results strongly suggested that the goat GHRHR gene is a candidate gene that influences growth traits in dairy goat.     

Identification of three novel SNPs and association with carcass traits in porcine <Emphasis Type="Italic">TNNI1</Emphasis> and <Emphasis Type="Italic">TNNI2</Emphasis>

In this study, two novel SNPs (EU743939:g.5174T>C in intron 4 and EU743939:g.8350C>A in intron 7) in TNNI1 and one SNP (EU696779:g.1167C>T in intron 3) in TNNI2 were identified by PCR–RFLP (PCR restriction fragment length polymorphism) using XbaI, MspI and SmaI restriction enzyme, respectively. The allele frequencies of three novel SNPs were determined in the genetically diverse pig breeds including ten Chinese indigenous pigs and three Western commercial pig breeds. Association analysis of the SNPs with the carcass traits were conducted in a Large White × Meishan F₂ pig population. The linkage of two SNPs (g.5174T>C and g.8350C>A) in TNNI1 gene had significant effect on fat percentage. Besides these, the g.5174T>C polymorphism was also significantly associated with skin percentage (P < 0.05), shoulder fat thickness (P < 0.05) and backfat thickness between sixth and seventh ribs (P < 0.05). The significant effects of g.1167C>T polymorphism in TNNI2 gene on fat percentage (P < 0.01), lean meat percentage (P < 0.05), lion eye area (P < 0.05), thorax–waist backfat thickness (P < 0.01) and average backfat thickness (P < 0.05) were also found.     

Molecular characterization,chromosomal location,alternative splicing and polymorphism of porcine <Emphasis Type="Italic">GFAT1</Emphasis> gene

Glutamine: fructose-6-phosphate amidotransferase (GFAT) is the rate-limiting enzyme of the hexosamine synthesis pathway, which plays important roles in insulin resistance and glucose toxicity. GFAT1 is one of the two isoenzymes of GFAT. In the present study, we cloned cDNA sequence of the porcine GFAT1 gene and identified a GFAT1 splice variant (designed GFAT1-L) that contains a 54 bp insertion within the coding region. Nested RT–PCR revealed that GFAT1 was ubiquitously expressed in all tested tissues, but GFAT1-L was only expressed in skeletal muscle and heart, not in liver, spleen, lung, kidney, small intestine, stomach and fat tissue, suggested that GFAT1-L was selectively expressed in striate muscle in pig. Using both the somatic cell hybrid panel and radiation hybrid panel, the GFAT1 gene was mapped to porcine chromosome 3q21-q27, in which several significant QTLs for carcass traits were found. Among the SNPs we found in porcine GFAT1 gene, only the g. 101A>G polymorphism which located in intron 8 was polymorphic in two pig populations we investigated in the study. Association analyses revealed that the g. 101A>G polymorphism has a significant effect on lean meat percentage (P < 0.05), corrected backfat thickness (P < 0.05) and backfat at the rump (P < 0.05).     

Analysis of the genetic architecture of maize kernel size traits by combined linkage and association mapping

Kernel size‐related traits are the most direct traits correlating with grain yield. The genetic basis of three kernel traits of maize, kernel length (KL), kernel width (KW) and kernel thickness (KT), was investigated in an association panel and a biparental population. A total of 21 single nucleotide polymorphisms (SNPs) were detected to be most significantly (P < 2.25 × 10^?6) associated with these three traits in the association panel under four environments. Furthermore, 50 quantitative trait loci (QTL) controlling these traits were detected in seven environments in the intermated B73 × Mo17 (IBM) Syn10 doubled haploid (DH) population, of which eight were repetitively identified in at least three environments. Combining the two mapping populations revealed that 56 SNPs (P < 1 × 10^?3) fell within 18 of the QTL confidence intervals. According to the top significant SNPs, stable‐effect SNPs and the co‐localized SNPs by association analysis and linkage mapping, a total of 73 candidate genes were identified, regulating seed development. Additionally, seven miRNAs were found to situate within the linkage disequilibrium (LD) regions of the co‐localized SNPs, of which zma‐miR164e was demonstrated to cleave the mRNAs of Arabidopsis CUC1, CUC2 and NAC6 in vitro. Overexpression of zma‐miR164e resulted in the down‐regulation of these genes above and the failure of seed formation in Arabidopsis pods, with the increased branch number. These findings provide insights into the mechanism of seed development and the improvement of molecular marker‐assisted selection (MAS) for high‐yield breeding in maize.     

Nucleotide variability and haplotype heterogeneity at the porcine fat mass and obesity‐associated (FTO) gene

A large number of studies have confirmed that variants within the fat mass and obesity‐associated (FTO) gene are associated with higher obesity risk in humans. We and others have shown that FTO polymorphisms are associated with fat deposition and related traits in several pig populations, thus confirming the role of this gene in fatness across species. However, some differences observed in different pig populations may be derived, at least in part, from genetic heterogeneity at this locus. Here, we characterise the nucleotide variability and haplotype diversity of the porcine FTO gene in breeds having different predispositions to fat deposition traits. We resequenced 4749 bp of coding and non‐coding regions of the porcine FTO gene in 44 pigs of eight different breeds and identified 27 single nucleotide polymorphisms (SNPs) and four insertions/deletions. A positive Tajima's D‐value (P < 0.10) obtained in Italian Duroc pigs may be compatible with putative balancing selection. From the sequenced pig panel, 20 haplotypes were inferred, some of which clustered according to the breed of origin (Meishan and Italian Duroc). Genetic heterogeneity at this locus could complicate the dissection of the effects of this gene on fat deposition and production traits in pigs. This situation resembles, to some extent, what has been reported in humans, thus making the study of the porcine FTO gene variability especially interesting, as it could be used as a model to understand the complex and elusive role of this gene in human obesity.     

Genetics of fat tissue accumulation in pigs: a comparative approach

Fatness traits are important in pig production since they influence meat quality and fattening efficiency. On the other hand, excessive fat accumulation in humans has become a serious health problem due to worldwide spread of obesity. Since the pig is also considered as an animal model for numerous human diseases, including obesity and metabolic syndrome, comparative genomic studies may bring new insights into genetics of fatness/obesity. Input of genetic factors into phenotypic variability of these traits is rather high and the heritability coefficient (h ²) of these traits oscillates around 0.5. Genome scanning revealed the presence of more than 500 QTLs for fatness in the pig genome. In addition to QTL studies, many candidate gene polymorphisms have been analyzed in terms of their associations with pig fatness, including genes encoding leptin (LEP) and its receptor (LEPR), insulin-like growth factor 2 (IGF-2), fatty acid-binding proteins (FABP3 andFABP4), melanocortin receptor type 4 (MC4R), and theFTO (fat mass and obesity-associated) gene. Among them, a confirmed effect on pig fatness was found for a well-known polymorphism of theIGF-2 gene. In humans the strongest association with predisposition to obesity was shown for polymorphism of theFTO gene, while in pigs such an association seems to be doubtful. The development of functional genomics has revealed a large number of genes whose expression is associated with fat accumulation and lipid metabolism, so far not studied extensively in terms of the association of their polymorphism with pig fatness. Recently, epigenomic mechanisms, mainly RNA interference, have been considered as a potential source of information on genetic input into the fat accumulation process. The rather limited progress in studies focused on the identification of gene polymorphism related with fatness traits shows that their genetic background is highly complex.     

Molecular characterization and association analysis of porcine adipose triglyceride lipase (<Emphasis Type="Italic">PNPLA2</Emphasis>) gene

The adipose triglyceride lipase (PNPLA2, also known as ATGL) is a novel triacylglycerol (TG) lipase which specifically removes the first fatty acid from the triglyceride molecule generating free fatty acid and diglyceride (DG) in mammalian cells. Here we describe the molecular characterization of the porcine ATGL gene. The full-length cDNA sequence contains a 1,461 bp open reading frame encoding a protein of 486 amino acids with a calculated molecular mass of 53.2 kDa and an isoelectric point of 7.90. The porcine ATGL protein shares high identity with other mammalian ATGL. The ATGL gene contains 9 coding exons, spans approximately 6 kb. The porcine ATGL mRNA was expressed predominantly in backfat, mildly in muscle, small intestine and heart, and almost absent in liver, spleen, lung, stomach, kidney and ovary. Statistical analysis showed the ATGL gene polymorphism (G/A³⁹²) was different between Chinese indigenous and introduced commercial western pig breeds, and was highly associated with almost all the fat deposition and carcass traits, including subcutaneous fat thickness, viscera adipose tissue, lean percentage, loin eye traits and even rib numbers.