A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes
BMC Genomics volume 13, Article number: 583 (2012)
Carcass fatness is an important trait in most pig breeding programs. Following market requests, breeding plans for fresh pork consumption are usually designed to reduce carcass fat content and increase lean meat deposition. However, the Italian pig industry is mainly devoted to the production of Protected Designation of Origin dry cured hams: pigs are slaughtered at around 160 kg of live weight and the breeding goal aims at maintaining fat coverage, measured as backfat thickness to avoid excessive desiccation of the hams. This objective has shaped the genetic pool of Italian heavy pig breeds for a few decades. In this study we applied a selective genotyping approach within a population of ~ 12,000 performance tested Italian Large White pigs. Within this population, we selectively genotyped 304 pigs with extreme and divergent backfat thickness estimated breeding value by the Illumina PorcineSNP60 BeadChip and performed a genome wide association study to identify loci associated to this trait.
We identified 4 single nucleotide polymorphisms with P≤5.0E-07 and additional 119 ones with 5.0E-07<P≤5.0E-05. These markers were located throughout all chromosomes. The largest numbers were found on porcine chromosomes 6 and 9 (n=15), 4 (n=13), and 7 (n=12) while the most significant marker was located on chromosome 18. Twenty-two single nucleotide polymorphisms were in intronic regions of genes already recognized by the Pre-Ensembl Sscrofa10.2 assembly. Gene Ontology analysis indicated an enrichment of Gene Ontology terms associated with nervous system development and regulation in concordance with results of large genome wide association studies for human obesity.
Further investigations are needed to evaluate the effects of the identified single nucleotide polymorphisms associated with backfat thickness on other traits as a pre-requisite for practical applications in breeding programs. Reported results could improve our understanding of the biology of fat metabolism and deposition that could also be relevant for other mammalian species including humans, confirming the role of neuronal genes on obesity.
Fat deposition is a key biological process that has important similarities between humans and pigs, potentially useful to elucidate mechanisms determining human obesity. This trait has practical and economical implications in pig breeding as it indirectly affects feeding efficiency and determines carcass value and consumers’ acceptance of pork.
Following consumer demands, breeding goals for fresh pork generally aim at reducing carcass fatness and increasing lean meat content which has adversely affected pork quality (e.g.).
The Italian pig breeding industry is mainly devoted to the production of high quality Protected Designation of Origin (PDO) dry cured hams for which pigs are raised until they reach about 160 kg live weight and appropriate fat coverage of the hams is required[2, 3]. Therefore, breeding objectives aim at maintaining fat coverage measured as backfat thickness (BFT). This objective has shaped the genetic pool of Italian heavy pig breeds for a few decades.
To investigate molecular genetic aspects of fat deposition in these pigs, we have recently applied a systematic candidate gene approach and have identified tens of single nucleotide polymorphisms (SNPs) associated with BFT and/or intermuscular fat content in Italian Large White and Italian Duroc pigs[4–10]. For example, a list of more than 30 SNPs has been associated with BFT in Italian Large White, including SNPs already found by other authors in IGF2, MC4R, TBC1D1, PPARG genes or newly identified in the PCSK1, ACP2, CALR, JAK3, and NT5E, among several other genes. Moreover, many other SNPs in additional candidate genes have been shown to explain a proportion of genetic variability of fat deposition traits in pigs[4, 15, 16]. In addition, a large number of QTLs for a variety of fat deposition and related traits have been already reported and listed in the Pig QTL database[17, 18].
Recently, with the development of a commercial high throughput SNP genotyping tool in pig (PorcineSNP60 BeadChip), a number of genome wide association (GWA) studies have been carried out in this species focusing on reproduction[20, 21], boar taint[22, 23], disease resistance, structural and body composition, including BFT.
We have herein applied a selective genotyping approach in the Italian Large White pig breed and genotyped extreme and divergent pigs for BFT estimated breeding value (EBV) by the Illumina PorcineSNP60 BeadChip (http://www.illumina.com) tool to identify chromosome regions and markers associated with BFT.
Results and discussion
A total of 304 performance tested Italian Large White pigs were genotyped with the Illumina PorcineSNP60 BeadChip, interrogating 62,163 loci. One pig was excluded from further analysis due to a call rate below 90%. A call rate ≥0.90 was obtained for 58,680 SNPs (for 2,293 SNP, call rate was 0.0; 1,190 SNP had 0.0<call rate<0.90). About 15.8% (9,287 SNPs) of these potentially useful SNPs had a minor allele frequency <0.05 and were discarded. The remaining 49,393 SNPs were re-mapped on the Sscrofa10.2 genome assembly.
Genome wide association (GWA) results
Only individuals with extreme phenotypes were genotyped for association study. Several authors have shown that this approach allows to attain the same power with less genotyped individuals (e.g.[26–28]). A recent GWA study for human obesity showed that this design can obtain very similar results to previous studies on general body mass index performed on unselected cohorts of tens of thousands of subjects.
In our study, genotyped pigs had extreme and divergent EBV for BFT: 151 had the lowest (thinnest BFT) and 152 the highest (thickest BFT) EBV. These animals were two generation unrelated gilts taken from the performance test of the National selection program of the Italian Large White breed carried out by the National Pig Breeders Association (ANAS).
Figure1 reports a Manhattan plot showing significant (P≤5.0E-07) and suggestively significant (5.0E-07<P≤5.0E-05) SNPs (Pnominal value thresholds for significant results were those indicated by the Wellcome Trust Case Control Consortium, WTCCC). Using these values, 4 SNPs were significantly associated (Table1) whereas 119 SNPs were suggestively associated with BFT (Additional file1: Table S1).
The WTCCC criteria to reduce the number of false positive are rather conservative: the suggestive threshold for significance (P=5.0E-05) corresponds to a FDR of 0.02. Had we assumed a false discovery rate (FDR) of 0.05 (Pnominal value = 0.000412), a total of 410 SNPs would had been below this threshold and therefore considered at least “suggestively associated” (data not shown).
Single nucleotide polymorphisms with P≤5.0E-05 were located in all porcine autosomal chromosomes (SCC), and on SSCX, and 5 SNPs were in unassembled scaffolds of the Sscrofa10.2 genome version. Among the mapped SNPs, the largest number was on SSC6 and SSC9 (n = 15), SSC4 (n = 13), SSC7 (n = 12) and SSC1 (n = 11) (Additional file1: Table S1). Twenty-two SNPs were in intronic regions of recognized genes in the Pre-Ensembl Sscrofa10.2 assembly. The closest gene for the remaining mapped SNPs (n = 96) was located in a range from 481 bp to 4.69 Mb (mean = 287.6 kb ± 580 kb, median = 88 kb).
The most significant SNP (ALGA0098168; P=3.07E-11) was on SSC18 (Table1). This SNP was localized in intron 3 of the phosphodiesterase 1C, calmodulin-dependent 70kDa (PDE1C) gene. PDEC1, highly expressed in brain and heart, is involved in the regulation of the cellular level of adenosine 3′,5′-cyclic monophosphate (cAMP) and guanosine 3′,5′-cyclic monophosphate (cGMP) that play critical roles in signal transduction. The second most significant SNP (M1GA0010276; P=1.45E-08) was localized on SSC7 at about 3.6 kb from the cysteine-rich secretory protein 1 (CRISP1) gene whose known function in reproduction processes is not directly linked to any fat or energy related biological function. The third top SNP (ALGA0109557; P=3.81E-07) was mapped on SSC15 at about 460 kb from the signal transducer and activator of transcription 4 (STAT4) gene. STAT4 is a member of the STAT family of transcription factors that transduces interleukin and type 1 interferon cytokine signals in T cells and monocytes, leading to important immunological functions. Reduction of STAT4 activation has been proposed to control obesity-induced inflammation. The fourth most significant marker (ALGA0069549; P=3.87E-07) was located on SSC13 at about 20 kb from the stabilin 1 (STAB1) gene. Another close marker (ALGA0109216; at position 38330168 of SSC13; Additional file1: Table S1) was suggestively significant (P=1.01E-06). The SSC13 region bracketed by these two SNPs includes the STAB1-nischarin (NISCH) gene interval that in human has been shown to be associated with waist-hip ratio (a measure of body fat distribution).
Several other genes close or within the additional suggestively significant SNPs (Additional file1: Table S1) have been already involved in obesity related biological mechanisms. Among this list it is worth mentioning: ATP-binding cassette, sub-family B (MDR/TAP), member 1 (ABCB1) gene (SSC9; ALGA0109564, P=9.01E-07) whose altered function contributes to steatosis and obesity in mice and a polymorphism in this gene has been associated with obesity risk in Japanese subjects; galanin receptor 3 (GALR3) gene (SSC5; M1GA0007458, P=1.25E-06) that is upregulated in adipose tissues of mice fed a high fat diet, and whose function is to bind galanin, a neuropeptide that regulates food intake, neurogenesis, memory, and gut secretion; olfactory receptor genes (two genes on SSC9, OR52N2 and OR56A3) have been associated with eating behaviour and adiposity in humans; Parkinson protein 2 (PARK2) gene on SSC1 (ALGA0108518, P=5.48E-06) that is regulated in a lipid-dependent manner and modulates systemic fat uptake via ubiquitin ligase-dependent effects; phosphodiesterase 4B, cAMP-specific (PDE4B) on SSC6 (ALGA0109354, P=5.95E-06) that has been already shown to be associated with BFT in pigs as well as with obesity in humans; vacuolar protein sorting 13 homolog B (yeast) (VPS13B) on SSC4 (ALGA0024658, P=3.00E-05) that causes Cohen syndrome, characterized by truncal obesity; iroquois homeobox 3 (IRX3) gene on SSC6 (M1GA0008432, P=4.66E-05), that is involved in the stress response after fat loss and could be linked to obesity and type 2 diabetes through its pancreatic function. Interestingly the second closest gene to this latter SSC6 SNP was FTO, that is well known to affect human obesity (i.e.).
Even though the annotation of the pig genome available at present in Pre-Ensembl should be considered preliminary, we further evaluated the potential functional role of regions around associated or suggestively associated SNPs with BFT in our pig population (Additional file1: Table S1). For this evaluation we used Gene Ontology (GO) information of their corresponding closest genes. Table2 reports GO terms enriched in this dataset. Interestingly, most of the statistically significant GO terms were related to nervous system development and regulation. This indication might support and extend the role of the nervous system in the genetic predisposition of fat accumulation in mammals, as in part reported in large GWA studies in humans and, subsequently, in pigs. Among the genes listed in these neuronal GO categories (Table2), few have been already reported to be indirectly associated or involved in obesity related traits. Apart from those already described above (IRX2 and PARK2), it is interesting to mention the delta-like 1 (Drosophila) gene (DLL1) as this gene is located in a quite large region (~1 Mb) associated to type 1 diabetes on human chromosome 6. For several other genes involved in neuronal processes, at present, there is no direct reported link with obesity or fat metabolism. It would be important to further explore their role in affecting the investigated phenotype as a possible strategy to identify new pathways and mechanisms affecting fat deposition. For example, it could be possible to speculatively suggest a relationship between dysbindin (DTNBP1), involved in the modulation of glutamatergic neurotransmission in the brain, schizophrenia and obesity.
Comparison with other studies in pigs
We compared our GWA results with results obtained in our previous candidate gene studies for BFT in pigs[4–10] and those obtained by other GWA and QTL mapping studies. In our previous studies[5, 10], the IGF2 intron3-g.3072G>A mutation was the most significant marker (P < 1.00E-50 by selective genotyping). As the IGF2 gene is not assembled in the Sscrofa10.2 genome version, it was not possible to obtain a direct comparison with results obtained for SNPs mapped on SSC2 included in the Illumina PorcineSNP60 BeadChip. However, no SNP in the region where IGF2 is likely to be found (0–10 Mb) reached the significance level of P<5.0E-05 (Additional file1: Table S1). Only one SNP (ASGA0008884, position 9139348; P=2.12E-04) was included in the list of markers with PFDR<0.05. Several other SSC2 SNPs were suggestively significant (Additional file1: Table S1) indicating that they might pick up other regions affecting fatness as already reported by QTL studies (e.g.[47, 48]) or candidate gene studies[5, 10, 49, 50].
The second most significant marker of our previous candidate gene investigation was the MC4R p.Asp298Asn substitution. In the current GWA study, no significant or suggestively significant SNPs were located in the SSC1 region around the MC4R gene, even if a few markers had P<1.0E-3 (data not shown). The GWA study by Fan et al., conducted on gilts of a commercial breeding stock, showed that markers around MC4R were significantly associated with 10th rib and last rib backfat. These slight differences in terms of level of significance of the markers between the two studies might be due to different MC4R haplotype structures in the two pig populations (Fontanesi et al. submitted) or to different positions in the pig body where BFT measurements were taken. However, in general, few results we obtained confirmed those previously obtained by Fan et al. in their GWA study on BFT. This could be due to different experimental designs, incomplete power in the two studies, and/or to differences between the investigated populations. Other results we previously obtained in candidate gene studies (i.e.) could be confirmed if we relaxed the significance threshold up to FDR <0.05 (data not shown).
QTLs for fat deposition traits can be found over all pig chromosomes. Many different studies have repeatedly reported the presence of complex QTL patterns for fat related traits in SSC1, SSC2, SSC4, SSC6 and SSC7. In the present GWA study, SSC4, SSC6, SSC7, and also SSC9 resulted to be rich in significant or suggestively significant markers (SSC4: expected proportion = 0.068, observed = 0.110; SSC6: expected = 0.059, observed = 0.127; SSC7: expected = 0.063, observed = 0.102; SSC9: expected = 0.061, observed = 0.127). These results seem to indicate these chromosomes to support an important proportion of genetic variability for BFT in the Italian Large White breed. In particular, two markers below the suggestive significance threshold were located both on IGSF3 or close to PKN2 on SSC4 and a few close blocks of SNPs with P<5.0E-05 (from about 65.1 - 65.4, 70.6 - 72.5, and 100.7 - 101.8 Mb) were located on SSC6 (Additional file1: Table S1). As mentioned above, FTO is close to the marker at position 28215213 on SSC6. Single marker analysis using a few FTO SNPs in our previous large association study with BFT in Italian Large White pigs did not produce significant results. However, subsequent haplotype analysis at this locus tended to confirm FTO as an important locus affecting fat deposition also in this pig breed.
This study is the first genome wide association analysis for BFT in Italian heavy pigs. The targeted trait is of paramount importance for the Italian pig breeding industry that is devoted to the production of high quality dry-cured hams for which fat coverage is a key factor during the processing and curing steps[2, 3]. The genetic dissection of BFT could open new perspectives to improve selection efficiency. In this study we applied a selective genotyping approach within the Italian Large White pig population to reduce the cost of genotyping without losing much power[26–31]. We took advantage of the large number of pigs that have been performance tested and genetically evaluated under the National selection program for this breed. The association analysis that compared SNP genotype frequencies between low BFT-EBV vs. high BFT-EBV groups identified 123 SNPs with P<5.0E-5 that were more densely represented in a few chromosomes known to harbor important QTLs for fat deposition traits. The quite large number of markers below this threshold (spread in different chromosome regions) might indirectly support the fact that many genes, each with a small-medium contribution, are involved in determining BFT, according to the classical definition of a quantitative trait.
Several significant or suggestively significant SNPs were close to genes whose function might be directly or indirectly related to energy metabolism and fat deposition. Many other cannot be easily linked to the targeted trait and might provide, if confirmed in following up studies, new evidence on this matter. Even if the annotation available in Pre-Ensembl for Sscrofa10.2 is preliminary, GO enrichment analysis indicated that neuronal genes might affect fat deposition in pig confirming and enlarging previous indications reported in humans.
Summarizing, as more information is becoming available in pigs on biological aspects of fat metabolism and deposition, it is more and more clear that this species could represent an attractive biomedical model for human obesity and associated diseases. Data here reported could give an insight over genetic mechanisms of fat metabolism and deposition that could be helpful in understanding also biology aspects of human obesity.
Animals and phenotypic traits
All animals used in this study were kept according to Italian and European legislation for pig production and all procedures described were in compliance with national and European Union regulations for animal care and slaughtering.
The national selection program of the Italian Large White breed is based on triplets of siblings from the same litter (two females and one castrated male) that are individually performance tested at the Central Test Station of the National Pig Breeder Association (ANAS) for the genetic evaluation of a boar from the same litter (sib-testing). Performance evaluation starts when the pigs are 30 to 45 days of age and it ends when the animals reach 155 ± 5 kg live weight. The nutritive level is quasi ad libitum, meaning that about 60% of the pigs are able to ingest the entire supplied ration. At the end of test, animals are transported to a commercial abattoir where they are slaughtered following standard procedures. Then, backfat thickness is measured on the carcasses at the level of Musculus gluteus medius.
The association study was conducted following a selective genotyping approach (e.g.[26–31]). In this study we genotyped two extreme and divergent groups of Italian Large White gilts of these triplets (one female per triplet), performance tested in the period 1996–2007. Two-generation unrelated females (i.e. gilts with different and unique parents) were chosen according to their EBV for BFT (152 with most negative and 152 with most positive EBV) within a performance tested population of ~12000 pigs (details of EBV calculation are reported below). The two extreme groups were chosen ranking the animals according to their BFT EBV and then taking only the first unrelated gilts in the list (with the most positive or the most negative BFT EBV). BFT EBV used to choose the animals were recalculated for the whole performance tested population in 2007. Average BFT EBV in the negative and positive selected groups of pigs were −9.8 ± 1.6 mm and +6.6 ± 2.3 mm, respectively. Genotyped pigs were a subset of the 560 two-generation unrelated pigs used in our previous candidate gene association study.
Genomic DNA was extracted from dried-blood by standard protocols. Based on quality control, all animals were used for genotyping using the PorcineSNP60 BeadChip developed by Illumina according to manufacturer’s protocol.
Estimated breeding values for BFT were calculated in the population using a BLUP-Multiple Trait-Animal Model that included the fixed effect of sex (considering the triplets of pigs from the same litter), batch on trial, inbreeding coefficient of the animal, interaction of sex by age at slaughtering, date of slaughtering and random effect of litter and animal. Three criteria were used to filter animals and SNP before association analysis: call rate >0.9 both at the 1) animal and 2) SNP level, and 3) MAF>0.05. Animals and SNPs that passed these filters were taken for association analysis treating the two groups as cases and controls. Full pedigree information available was used to obtain a kinship matrix. In order to correct for possible family-based stratification (see Additional file2: Figure S1), the EIGENSTRAT method was applied including the kinship matrix, and association tests were performed. All analyses were performed in R, using an option of the package GenABEL for computing the test-statistics according to the EIGENSTRAT method, and the package kinship for building the pedigree kinship matrix.
For n animals, the first K<n principal components, c1, …, c K , of the kinship matrix among the animals were used as axes of genetic variation. Let g ij and p j be the genotype at SNP i (g ij = 0, 1 or 2) and the phenotype of animal j, respectively, a PC-based adjustment was performed on genotypes and phenotypes according to the following formulas:
where c kj is the score of the k-th component on animal j, β ki and γ k are the partial regression coefficients for predicting the i-th genotype and the phenotype, respectively, on the basis of the k-th component (with k = 1, ., K).
The association test-statistic is computed as (n − K − 1)r i 2, where
is the squared correlation coefficient between the i-th PC-adjusted genotype and PC-adjusted phenotype. As noted by Price et al., this statistic is a generalization of the Armitage trend statistic for discrete genotypes and phenotypes.
Wellcome Trust Case Control Consortium significance thresholds, whose definition depends on the prior odds and power, were used in this study. In addition, correction for multiple testing was achieved by using a False Discovery Rate approach. For each chromosome, the expected proportion of SNPs with P<5.0E-5 was computed under the assumption of uniform distribution from the informative SNPs over the chromosome. This proportion was compared to the proportion of significant or suggestively significant markers actually observed on the same chromosome.
Mapping of the PorcineSNP60 BeadChip SNPs was obtained by using BWA on the Sscrofa9.2 and Sscrofa10.2 genome assemblies as previously described and confirmed using the BLAT analysis available at the Ensembl (http://www.ensembl.org/Sus_scrofa/Info/Index) and Pre-Ensembl (http://pre.ensembl.org/Sus_scrofa/Info/Index) databases (February 2012). Coordinates for the Sscrofa10 genome preliminary version (September 2010) were downloaded from the Animal Genome repository web sitehttp://www.animalgenome.org/repository/. Identification of the closest genes to SNPs with P<5.0E-05 was obtained using Pre-Ensembl annotation of Sscrofa10.2 genome version and verified using Ensembl Sscrofa9.2 genome version (February 2012). Starting from the corresponding protein sequences retrieved from these databases, the corresponding gene symbols were extracted from NCBI Gene section (http://www.ncbi.nlm.nih.gov/) and/or Uniprot (http://www.uniprot.org/) databases (February 2012). Gene annotation was verified by BLAST analysis (http://blast.ncbi.nlm.nih.gov/). Gene Ontology analysis was carried out using DAVID Bioinformatics Resources 6.7 (http://david.abcc.ncifcrf.gov/).
Lonergan SM, Huff-Lonergan E, Rowe LJ, Kuhlers DL, Jungst SB: Selection for lean growth efficiency in Duroc pigs influences pork quality. J Anim Sci. 2001, 79: 2075-2085.
Russo V, Nanni Costa L: Suitability of pig meat for salting and the production of quality processed products. Pig News Inform. 1995, 16: 7N-26N.
Bosi P, Russo V: The production of the heavy pig for high quality processed products. Ital J Anim Sci. 2004, 3: 309-321.
Fontanesi L, Scotti E, Buttazzoni L, Davoli R, Russo V: The porcine fat mass and obesity associated (FTO) gene is associated with fat deposition in Italian Duroc pigs. Anim Genet. 2009, 40: 90-93. 10.1111/j.1365-2052.2008.01777.x.
Fontanesi L, Speroni C, Buttazzoni L, Scotti E, Dall'Olio S, Nanni Costa L, Davoli R, Russo V: The insulin-like growth factor 2 (IGF2) gene intron3-g.3072G>A polymorphism is not the only Sus scrofa chromosome 2p mutation affecting meat production and carcass traits in pigs: evidence from the effects of a cathepsin D (CTSD) gene polymorphism. J Anim Sci. 2010, 88: 2235-2245. 10.2527/jas.2009-2560.
Fontanesi L, Colombo M, Scotti E, Buttazzoni L, Bertolini F, Dall’Olio S, Davoli R, Russo V: The porcine tribbles homolog 3 (TRIB3) gene: identification of a missense mutation and association analysis with meat quality and production traits in Italian heavy pigs. Meat Sci. 2010, 86: 806-813.
Fontanesi L, Scotti E, Speroni C, Buttazzoni L, Russo V: A selective genotyping approach identifies single nucleotide polymorphisms in porcine chromosome 2 genes associated with production and carcass traits in Italian heavy pigs. Ital J Anim Sci. 2011, 10: e15-
Fontanesi L, Colombo M, Tognazzi L, Scotti E, Buttazzoni L, Dall’Olio S, Davoli R, Russo V: The porcine TBC1D1 gene: mapping, SNP identification, and association study with meat, carcass and production traits in Italian heavy pigs. Mol Biol Rep. 2011, 38: 1425-1431. 10.1007/s11033-010-0247-3.
Fontanesi L, Speroni C, Buttazzoni L, Scotti E, Dall’Olio S, Davoli R, Russo V: Association between polymorphisms in cathepsin and cystatin genes with meat production and carcass traits in Italian Duroc pigs: confirmation of the effects of a cathepsin L (CTSL) gene marker. Mol Biol Rep. 2012, 39: 109-115. 10.1007/s11033-011-0715-4.
Fontanesi L, Galimberti G, Calò DG, Fronza R, Martelli PL, Scotti E, Colombo M, Schiavo G, Casadio R, Buttazzoni L, Russo V: Identification and association analysis of several hundred single nucleotide polymorphisms within candidate genes for back fat thickness in Italian Large White pigs using a selective genotyping approach. J Anim Sci. 2012, 90: 2450-2464. 10.2527/jas.2011-4797.
Van Laere A-S, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, Archibald AL, Haley CS, Buys N, Tally M, Andersson G, Georges M, Andersson L: A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature. 2003, 425: 832-836. 10.1038/nature02064.
Kim KS, Larsen N, Short T, Plastow G, Rothschild MF: A missense variant of the porcine melanocortin-4 receptor (MC4R) gene is associated with fatness, growth, and feed intake traits. Mamm Genome. 2000, 11: 131-135. 10.1007/s003350010025.
Fan B, Onteru SK, Mote BE, Serenius T, Stalder KJ, Rothschild MF: Large-scale association study for structural soundness and leg locomotion traits in the pig. Genet Sel Evol. 2009, 41: 14-10.1186/1297-9686-41-14.
Fontanesi L, Bertolini F, Scotti E, Trevisi P, Buttazzoni L, Dall’Olio S, Davoli R, Bosi P, Russo V: Polymorphisms in an obesity related gene (PCSK1) are associated with fat deposition and production traits in Italian heavy pigs. Animal. 2012, 6: 1913-1924. 10.1017/S1751731112001280.
Fan B, Du ZQ, Rothschild MF: The fat mass and obesity-associated (FTO) gene is associated with intramuscular fat content and growth rate in the pig. Anim Biotechnol. 2009, 20: 58-70. 10.1080/10495390902800792.
Switonski M, Stachowiak M, Cieslak J, Bartz M, Grzes M: Genetics of fat tissue accumulation in pigs: a comparative approach. J Appl Genet. 2010, 51: 153-168. 10.1007/BF03195724.
Rothschild MF, Hu Z-L, Jiang Z: Advances in QTL mapping in pigs. Int J Biol Sci. 2007, 3: 192-197.
Hu Z-L, Reecy JM: Animal QTLdb: Beyond a repository – A public platform for QTL comparisons and integration with diverse types of structural genomic information. Mamm Genome. 2007, 18: 1-4. 10.1007/s00335-006-0105-8.
Ramos AM, Crooijmans RP, Affara NA, Amaral AJ, Archibald AL, Beever JE, Bendixen C, Churcher C, Clark R, Dehais P, Hansen MS, Hedegaard J, Hu ZL, Kerstens HH, Law AS, Megens HJ, Milan D, Nonneman DJ, Rohrer GA, Rothschild MF, Smith TP, Schnabel RD, Van Tassell CP, Taylor JF, Wiedmann RT, Schook LB, Groenen MA: Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology. PLoS One. 2009, 4: e6524-10.1371/journal.pone.0006524.
Uimari P, Sironen A, Sevon-Aimonen ML: Whole-genome SNP association analysis of reproduction traits in the Finnish Landrace pig breed. Genet Sel Evol. 2011, 43: 42-10.1186/1297-9686-43-42.
Onteru SK, Fan B, Du ZQ, Garrick DJ, Stalder KJ, Rothschild MF: A whole-genome association study for pig reproductive traits. Anim Genet. 2012, 43: 18-26. 10.1111/j.1365-2052.2011.02213.x.
Grindflek E, Lien S, Hamland H, Hansen MH, Kent M, van Son M, Meuwissen TH: Large scale genome-wide association and LDLA mapping study identifies QTLs for boar taint and related sex steroids. BMC Genomics. 2011, 12: 362-10.1186/1471-2164-12-362.
Gregersen VR, Conley LN, Sorensen KK, Guldbrandtsen B, Velander IH, Bendixen C: Genome-wide association scan and phased haplotype construction for quantitative trait loci affecting boar taint in three pig breeds. BMC Genomics. 2012, 13: 22-10.1186/1471-2164-13-22.
Boddicker N, Waide EH, Rowland RR, Lunney JK, Garrick DJ, Reecy JM, Dekkers JC: Evidence for a major QTL associated with host response to Porcine Reproductive and Respiratory Syndrome virus challenge. J Anim Sci. 2012, 10.2527/jas.2011-4464.
Fan B, Onteru SK, Du ZQ, Garrick DJ, Stalder KJ, Rothschild MF: Genome-wide association study identifies Loci for body composition and structural soundness traits in pigs. PLoS One. 2011, 6: e14726-10.1371/journal.pone.0014726.
Darvasi A, Soller M: Selective genotyping for determination of linkage between a marker locus and a quantitative trait locus. Theor Appl Genet. 1992, 85: 353-359.
Van Gestel S, Houwing-Duistermaat JJ, Adolfsson R, van Duijn CM, Van Broeckhoven C: Power of selective genotyping in genetic association analyses of quantitative traits. Behav Genet. 2000, 30: 141-146. 10.1023/A:1001907321955.
Zhang G, Nebert DW, Chakraborty R, Jin L: Statistical power of association using the extreme discordant phenotype design. Pharmacogenet Genomics. 2006, 16: 401-413. 10.1097/01.fpc.0000204995.99429.0f.
Paternoster L, Evans DM, Nohr EA, Holst C, Gaborieau V, Brennan P, Gjesing AP, Grarup N, Witte DR, Jørgensen T, Linneberg A, Lauritzen T, Sandbaek A, Hansen T, Pedersen O, Elliott KS, Kemp JP, St Pourcain B, McMahon G, Zelenika D, Hager J, Lathrop M, Timpson NJ, Smith GD, Sørensen TI: Genome-wide population-based association study of extremely overweight young adults - the GOYA study. PLoS One. 2011, 6: e24303-10.1371/journal.pone.0024303.
Wellcome Trust Case Control Consortium: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447: 661-683. 10.1038/nature05911.
Lakics V, Karran EH, Boess FG: Quantitative comparison of phosphodiesterase mRNA distribution in human brain and peripheral tissues. Neuropharmacology. 2010, 59: 367-374. 10.1016/j.neuropharm.2010.05.004.
Chakrabarti SK, Wen Y, Dobrian AD, Cole BK, Ma Q, Pei H, Williams MD, Bevard MH, Vandenhoff GE, Keller SR, Gu J, Nadler JL: Evidence for activation of inflammatory lipoxygenase pathways in visceral adipose tissue of obese Zucker rats. Am J Physiol Endocrinol Metab. 2011, 300: E175-E187. 10.1152/ajpendo.00203.2010.
Heid IM, Jackson AU, Randall JC, Winkler TW, Qi L, Steinthorsdottir V, Thorleifsson G, Zillikens MC, Speliotes EK, Mägi R, Workalemahu T, White CC, Bouatia-Naji N, Harris TB, Berndt SI, Ingelsson E, Willer CJ, Weedon MN, Luan J, Vedantam S, Esko T, Kilpeläinen TO, Kutalik Z, Li S, Monda KL, Dixon AL, Holmes CC, Kaplan LM, Liang L, Min JL, et al, et al: Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. 2010, 42: 949-960. 10.1038/ng.685.
Foucaud-Vignault M, Soayfane Z, Ménez C, Bertrand-Michel J, Martin PG, Guillou H, Collet X, Lespine A: P-glycoprotein dysfunction contributes to hepatic steatosis and obesity in mice. PLoS One. 2011, 6: e23614-10.1371/journal.pone.0023614.
Ichihara S, Yamada Y, Kato K, Hibino T, Yokoi K, Matsuo H, Kojima T, Watanabe S, Metoki N, Yoshida H, Satoh K, Aoyagi Y, Yasunaga A, Park H, Tanaka M, Nozawa Y: Association of a polymorphism of ABCB1 with obesity in Japanese individuals. Genomics. 2008, 91: 512-516. 10.1016/j.ygeno.2008.03.004.
Kim A, Park T: Diet-induced obesity regulates the galanin-mediated signalling cascade in the adipose tissue of mice. Mol Nutr Food Res. 2010, 54: 1361-1370. 10.1002/mnfr.200900317.
Choquette AC, Bouchard L, Drapeau V, Lemieux S, Tremblay A, Bouchard C, Vohl MC, Pérusse L: Association between olfactory receptor genes, eating behaviour traits and adiposity: results from the Quebec Family Study. Physiol Behav. 2012, 105: 772-776. 10.1016/j.physbeh.2011.10.015.
Kim KY, Stevens MV, Akter MH, Rusk SE, Huang RJ, Cohen A, Noguchi A, Springer D, Bocharov AV, Eggerman TL, Suen DF, Youle RJ, Amar M, Remaley AT, Sack MN: Parkin is a lipid-responsive regulator of fat uptake in mice and mutant human cells. J Clin Invest. 2011, 121: 3701-3712. 10.1172/JCI44736.
Lee KT, Byun MJ, Kang KS, Park EW, Lee SH, Cho S, Kim H, Kim KW, Lee T, Park JE, Park W, Shin D, Park HS, Jeon JT, Choi BH, Jang GW, Choi SH, Kim DW, Lim D, Park HS, Park MR, Ott J, Schook LB, Kim TH, Kim H: Neuronal genes for subcutaneous fat thickness in human and pig are identified by local genomic sequencing and combined SNP association study. PLoS One. 2011, 6: e16356-10.1371/journal.pone.0016356.
Balikova I, Lehesjoki AE, de Ravel TJ, Thienpont B, Chandler KE, Clayton-Smith J, Träskelin AL, Fryns JP, Vermeesch JR: Deletions in the VPS13B (COH1) gene as a cause of Cohen syndrome. Hum Mutat. 2009, 30: E845-E854. 10.1002/humu.21065.
Dankel SN, Fadnes DJ, Stavrum AK, Stansberg C, Holdhus R, Hoang T, Veum VL, Christensen BJ, Våge V, Sagen JV, Steen VM, Mellgren G: Switch from stress response to homeobox transcription factors in adipose tissue after profound fat loss. PLoS One. 2010, 5: e11033-10.1371/journal.pone.0011033.
Ragvin A, Moro E, Fredman D, Navratilova P, Drivenes O, Engström PG, Alonso ME, de la Calle Mustienes E, Gómez Skarmeta JL, Tavares MJ, Casares F, Manzanares M, van Heyningen V, Molven A, Njølstad PR, Argenton F, Lenhard B, Becker TS: Long-range gene regulation links genomic type 2 diabetes and obesity risk regions to HHEX, SOX4, and IRX3. Proc Natl Acad Sci USA. 2010, 107: 775-780. 10.1073/pnas.0911591107.
Dina C, Meyre D, Gallina S, Durand E, Körner A, Jacobson P, Carlsson LM, Kiess W, Vatin V, Lecoeur C, Delplanque J, Vaillant E, Pattou F, Ruiz J, Weill J, Levy-Marchal C, Horber F, Potoczna N, Hercberg S, Le Stunff C, Bougnères P, Kovacs P, Marre M, Balkau B, Cauchi S, Chèvre JC, Froguel P: Variation in FTO contributes to childhood obesity and severe adult obesity. Nat Genet. 2007, 39: 724-726. 10.1038/ng2048.
Willer CJ, Speliotes EK, Loos RJ, Li S, Lindgren CM, Heid IM, Berndt SI, Elliott AL, Jackson AU, Lamina C, Lettre G, Lim N, Lyon HN, McCarroll SA, Papadakis K, Qi L, Randall JC, Roccasecca RM, Sanna S, Scheet P, Weedon MN, Wheeler E, Zhao JH, Jacobs LC, Prokopenko I, Soranzo N, Tanaka T, Timpson NJ, Almgren P, Bennett A, Genetic Investigation of ANthropometric Traits Consortium: Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet. 2009, 41: 25-34. 10.1038/ng.287.
Bradfield JP, Qu HQ, Wang K, Zhang H, Sleiman PM, Kim CE, Mentch FD, Qiu H, Glessner JT, Thomas KA, Frackelton EC, Chiavacci RM, Imielinski M, Monos DS, Pandey R, Bakay M, Grant SF, Polychronakos C, Hakonarson H: A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci. PLoS Genet. 2011, 7: e1002293-10.1371/journal.pgen.1002293.
Numakawa T, Yagasaki Y, Ishimoto T, Okada T, Suzuki T, Iwata N: Evidence of novel neuronal functions of dysbindin, a susceptibility gene for schizophrenia. Hum Mol Genet. 2004, 13: 2699-26708. 10.1093/hmg/ddh280.
Lee SS, Chen Y, Moran C, Cepica S, Reiner G, Bartenschlager H, Moser G, Geldermann H: Linkage and QTL mapping for Sus scrofa chromosome 2. J Anim Breed Genet. 2003, 120 (Suppl. 1): 11-19.
Tortereau F, Gilbert H, Heuven HC, Bidanel JP, Groenen MA, Riquet J: Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees. Genet Sel Evol. 2011, 43: 11-10.1186/1297-9686-43-11.
Cepica S, Ovilo C, Masopust M, Knoll A, Fernandez A, Lopez A, Rohrer GA, Nonneman D: Four genes located on a SSC2 meat quality QTL region are associated with different meat quality traits in Landrace × Chinese-European crossbred population. Anim Genet. 2012, 43: 333-336. 10.1111/j.1365-2052.2011.02252.x.
Fontanesi L, Buttazzoni L, Scotti E, Russo V: Confirmation of the association between a single nucleotide polymorphism in the porcine LDHA gene and average daily gain and correlated traits in Italian Large White pigs. Anim Genet. 2012, 43: 649-650. 10.1111/j.1365-2052.2012.02355.x.
Fontanesi L, Russo V: Nucleotide variability and haplotype heterogeneity at the porcine fat mass and obesity associated (FTO) gene. Anim Genet. 10.1111/j.1365-2052.2012.02352.x.
Fontanesi L, Davoli R, Nanni Costa L, Beretti F, Scotti E, Tazzoli M, Tassone F, Colombo M, Buttazzoni L, Russo V: Investigation of candidate genes for glycolytic potential of porcine skeletal muscle: association with meat quality and production traits in Italian Large White pigs. Meat Sci. 2008, 80: 780-787. 10.1016/j.meatsci.2008.03.022.
Steemers FJ, Gunderson KL: Whole genome genotyping technologies on the BeadArray platform. Biotechnol J. 2007, 2: 41-49. 10.1002/biot.200600213.
Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.
R Development Core Team: R: A language and environment for statistical computing.http://www.R-project.org/,
Aulchenko YS, Ripke S, Isaacs A, van Duijn CM: GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007, 23: 1294-1296. 10.1093/bioinformatics/btm108.
Atkinson B, Therneau T: kinship: mixed-effects Cox models, sparse matrices, and modeling data from large pedigrees. R package version 1.1.0-23, 2009:http://CRAN.R-project.org/package=kinship,
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc Series B. 1995, 57: 289-300.
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.
Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.
We thank ANAS for providing samples and data. This work was supported by University of Bologna RFO and FAGenomicH projects and by Italian MiPAAF INNOVAGEN funds.
The authors declare that they have no competing interests.
LF conceived and coordinated the study, analysed data and drafted the manuscript. GS, GG, DGC, and PLM performed statistical and bioinformatics analyses. ES carried out laboratory activities. LB, RC and VR coordinated and conceived the study. All authors reviewed and contributed to draft the manuscript. All authors read and approved the final manuscript.
Giuseppina Schiavo, Giuliano Galimberti, Daniela Giovanna Calò contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1. Suggestively significant SNPs (5.0E-07< P≤5.0E-05), their chromosome positions and their closest genes in Sscrofa10.2 (Pre-Ensembl). Notes are the same as those reported for Table1. (JPEG 86 KB)
Additional file 2: Figure S1. Two-dimensional graphical representation of relatedness among animals based on a multidimensional scaling representation of pedigree-based kinship matrix. Different symbols are used to denote pigs with positive or negative backfat thickness estimated breeding values. (XLS 46 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Fontanesi, L., Schiavo, G., Galimberti, G. et al. A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes. BMC Genomics 13, 583 (2012). https://doi.org/10.1186/1471-2164-13-583