A comprehensive survey of copy number variation in 18 diverse pig populations and identification of candidate copy number variable genes associated with complex traits
- Congying Chen†1,
- Ruimin Qiao†1,
- Rongxing Wei†1, 2,
- Yuanmei Guo1,
- Huashui Ai1,
- Junwu Ma1,
- Jun Ren1 and
- Lusheng Huang1Email author
© Chen et al.; licensee BioMed Central Ltd. 2012
Received: 28 February 2012
Accepted: 15 December 2012
Published: 27 December 2012
Copy 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.
Here 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.
We 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.
KeywordsCopy number variation Copy number variable gene Complex trait QTL Pig
Copy number variation (CNV), a major type of structural variants, is defined as DNA segments that vary from one kilobase to several megabases in length and present at variable copy numbers in comparison with a reference genome[1, 2]. Following the completion of pig whole-genome sequencing, genome-wide polymorphisms including CNV, SNP and deletion/insertion will be well annotated in the near future. Currently, Porcine 60K SNP BeadChips are commercially available for genome-wide analyses of 62,163 SNPs. However, a comprehensive study of genome-scale CNVs in pigs remains unexplored.
Tilling oligonucleotide array and comparative genomic hybridization (CGH) array have been commonly used to detect whole-genome CNVs. In recent years, whole-genome SNP genotyping arrays offer alternative methods for CNV detection. Computer programs like PennCNV, QuantiSNP, cnvPartition (Illumina) and CNV Workshop have been developed to identify CNVs from SNP array data. The unique ability to integrate family relationships from parent-offspring trios, total signal intensity and allelic intensity ratio at each SNP marker, and the allele frequency of SNP makes PennCNV algorithm have a moderate power and the lowest false positive rate. Winchester et al. (2009) reported that PennCNV is the most accurate program in the prediction of CNVs for the Illumina’s platform by comparing different algorithms for CNV detection.
CNVs have been commonly identified in humans, rats, dogs, cattle and horse, and occupy about 3.7%, 1.4%, 4.2%, 4.6% and 3.6% of their assembled genome, respectively[8–12]. It has been estimated that CNVs account for at least 17.7% of heritable variation of gene expression in a variety of ways including gene dosage effects, disruption of gene coding region and deletion or duplication of regulatory elements. In humans, SNP-tagged CNVs are enriched for expression quantitative trait loci (eQTL). CNVs have been confirmed to associate with Mendelian diseases and complex genetic disorders in humans, such as schizophrenia, body mass index, Crohn’s disease and intellectual disability and various congenital defects. Similarly, in livestock, more and more studies evidenced that CNVs play causative effects on phenotypic variations, such as CNV in intron 1 of SOX5 causing the pea-comb phenotype in chickens, a 4.6-kb intronic duplication in STX17 for hair greying and melanoma in horses, duplication of FGF3, FGF4, FGF19 and ORAOV1 resulting in hair ridge and predisposition to dermoid sinus in Ridgeback dogs, and CNV and missense mutations of the agouti signaling protein (ASIP) gene leading to different coat colors in goats. But in pigs, few such examples have been reported at present. KIT is the first pig gene that has been confirmed that gene duplication and a splice mutation leading the skipping of exon 17 are responsible for the dominant white phenotype and peripheral blood cell[23, 24].
Until now, to our knowledge, there are only three studies on pig CNV discovery. Fadista et al. (2008) found 37 CNVRs across chromosome (SSC) 4, 7, 14 and 17 using a custom tilling oligonucleotide array. Tang et al. (2010) investigated the CNV distribution on SSC7 and SSC8 by CGH array. More recently, Ramayo-Caldas et al. (2010) detected 49 CNVRs in porcine autosomal chromosomes in 55 animals from an Iberian × Landrace cross with Porcine SNP60 BeadChip. However, the distribution of CNVs in large scale and diverse pig populations remains largely unknown.
We herein used Porcine SNP60 BeadChip and PennCNV algorithm to identify porcine autosomal CNVs in 1,693 animals from 18 populations, and analyzed the CNV distribution in pig genome and different populations. We compared the identified CNVs with the reported porcine CNV call sets and investigated the copy number variable genes. Especially, we used a large scale White Duroc × Erhualian F2 intercross, in which the QTL were mapped for 422 traits. The intercross allowed us to systematically investigate the effects of pig CNVs on phenotypic variations.
Results and discussion
CNV discovery, distribution and validation
Identification of CNVs in 18 diverse pig populations
Number of animals identified CNV
Number of CNVs
Status of CNVs
White Duroc × Erhualian F2 intercross
Bamaxiang × Erhualian F2 intercross
Rongchang × Erhualian F2 intercross
Shazilin × Erhualian F2 intercross
Tongcheng × Erhualian F2 intercross
Just as the cases in human and cattle, we found that the CNVRs were non-randomly distributed across the pig genome. Chromosome 4, for instance, has only 2.02% of sequences showing copy number variable, while chromosome 18 has > 18.01% of sequences with copy number variation. Several “hotspots” of copy number variation were obsevered in this study, such as SSC6: 138.40-145.74 Mb, SSC11: 43.58-46.62 Mb, SSC13:213.45-215.94 Mb, SSC14: 2.74-7.42 Mb and SSC17: 8.25-10.28 Mb. These regions contain clusters of four to six CNVRs, indicating CNV hotspots (Additional file1: Table S1; Figure1).
Porcine CNV frequencies among breeds
Like CNV frequencies in humans, CNV numbers differ greatly among different pig populations. The average number of CNVs per population was 117.89, ranging from 23 (White Duroc) to 683 (White Duroc × Erhualian F2 intercross). The most number of CNVs per sample was detected in Duroc pigs (12.40 CNVs per animal on average), in comparison with the least number of 0.91 CNVs per animal in the White Duroc × Erhualian F2 intercross (Table1).
Of the 565 CNVRs, only 20 CNVRs were detected in more than 50% of populations (9 populations). Similar to the finding in humans, most CNVRs (72.87%) were restricted to one population. This should be due to sampling variances or the fact that they were recent evolution events. Among 18 diverse populations, the largest number of 310 CNVRs was identified in the White Duroc × Erhualian F2 intercross (Additional file4: Table S2). It is most likely due to the fact that this population was a cross between European and Asian divergent breeds. Of these 310 CNVRs, 215 were unique to this resource population. As two Western pig breeds (Duroc and White Duroc) used in this study, we found that all 5 CNVRs identified in White Duroc were included within 42 CNVRs detected in Duroc. It is consistent with the fact that White Duroc is one of specialized lines of Duroc. Total 26 CNVRs were unique to the two Western breeds. Among Chinese indigenous pig breeds, the most number of 30 CNVRs was identified in both Jinhua and Shazilin. Moreover, Jinhua had the highest average number of CNVs per individual (5.36 per individual; Additional file4: Table S2).
Gene content of CNV regions
The BioMart gene database based on the Sscrofa 10.2 reference genome assembly was used to retrieve genes within the detected CNVRs. A total of 1,764 transcripts were annotated within 320 out of 565 putative CNVRs, including 1,546 transcripts completely located within the CNVRs and 218 transcripts overlapping with the CNVRs. These 1,764 transcripts were composed of 1,587 protein coding genes, 37 pseudogenes, 1 retrotransposed genes, 42 miRNAs, 7 rRNAs, 36 snoRNAs, 47 snRNAs and 7 miscRNAs. No annotated transcripts were identified within the other 235 CNVRs. Of the 1,587 protein coding genes, 1,055 genes were well annotated in pigs including 634 genes that have CNVs in the human genome (Additional file5: Table S4). The average gene number per Mb in CNVRs was 13.26, 14.74 and 6.01, respectively, for gain, loss and gain/loss. However, compared to the average gene number per Mb in whole genome, the CNVRs have higher gene density (11.21 vs. 8.37). This result was consistent with the findings in other species in which a higher gene content was discovered in CNVRs[9, 29, 30].
We found that many CNV-associated genes appear to be certain gene clusters or families, such as olfactory receptor family, solute carrier family, apolipoprotein gene family, myosin gene family, CD gene family, cytochrome C oxidase gene family, interleukin gene family, protocadherin gamma gene cluster, beta defensin protein family, neuroblastoma breakpoint family, zinc finger protein family and ring finger protein family (Additional file5: Table S4). Some of these gene families have been well characterized and play important roles in biological processes[27, 35]. For example, olfactory receptor family is the most well characterized CNV-related genes in humans. Over 400 human olfactory receptor genes are reported to be variable in copy number. In this study, total 84 olfactory receptor genes were located within CNVRs. More than 15 members of solute carrier family were detected in CNVRs in this study. Solute carrier family encodes membrane transport proteins including over 300 members organized into 51 families. The family proteins play important roles in transportation and exchange of ion, amino acid and other substance which take part in important biological process.
Identification of copy number variable genes as potential candidate genes for complex traits
Identification of potential candidate CNV genes for complex traits in the White Duroc × Erhualian F 2 population
Trait of overlapped QTL*
Phenotype in knockout mice
decreased body size
decreased subcutaneous adipose tissue
abdominal fat weight
decreased retroperitoneal fat pad weight
length of scapular
decreased length of long bones
intermuscle fat content of LD
decreased white adipose tissue amount
body weight at 240 days
abnormal developmental patterning
Glycolytic potential of LD
decreased skeletal muscle glycogen level
mean corpuscular volume
increased mean corpuscular volume
mean corpuscular hemoglobin
decreased hemoglobin content
diluted coat color
increased bone size and stiffnes
We further chose KIT as a proof-of-principle example as the confirmed association between KIT duplication and MCH and MCV. Previous studies have confirmed that KIT regulatory mutations including the gene duplication and splice mutation are responsible for the dominant white phenotype in pigs and have pleiotropic effects on peripheral blood cell measures in Western commercial pigs[23, 24]. A significant QTL for MCH and MCV at day 240 was detected at SSC8: 43,550,231 in the White Duroc × Erhualian F2 resource population, which fell within the genomic region of CNVR268 (SSC8: 43,425,758-43,955,459; Table2). More than 378 F2 animals from the intercross inherited this CNV variant. Two genes of KIT and KDR were located within this CNVR. The KIT knockout mice exhibited phenotypes of increased mean corpuscular volume, decreased hemoglobin content and diluted coat color. Furthermore, association analysis showed that CNV268 was significantly associated with MCH and MCV in the White Duroc × Erhualian F2 resource population (P = 1.17 × 10-5). And if we included the CNV268 as fixed effect in QTL mapping, the QTL on SSC8 for MCH and MCV was never detected again.
Moreover, the result obtained in this study was also consistent with the causative relation between KIT duplication and dominant white coat color identified before. The CNVR268 harboring KIT was detected only in the solid white breed White Duroc. It was absent in all other pigs from diverse populations having colored phenotypes. It is noteworthy that the CNVR was either not found in Chinese belt-like breeds including Bamaxiang, Dongshan, Shanggao, Jinhua, Shazilin and Tongcheng, or in Rongchang pigs with the white coat color. This was in agreement with our previous conclusions that the belt-like and white coat colors in Chinese pigs were not caused by the dominant white allele of KIT[57, 58].
Although some identified copy number variable genes were not overlapped with our reported QTL, they have been reported to associate with complex traits in pigs, humans or mice. For instance, we detected an 836.67-kb CNV in SSC6: 49,802,217-50,638,891 in the White Duroc × Erhualian resource population. This region contains the alpha-1-fucosyltransferase (FUT1) gene. FUT1 has been identified as a strong candidate gene encoding the intestinal Escherichia coli F18 receptor that determines susceptibility to oedema disease, post-weaning diarrhoea in Western piglets and total number of born piglets[59, 60]. APOE/C4/C2 gene cluster is located within a 427.46-kb CNV region on SSC6: 46,893,592-47,321,053. The APOE/C1/C4/C2 gene cluster variation in humans is associated with plasma lipids, particularly low density lipoprotein (LDL) level and coronary heart disease. CNVR128 (SSC3: 120,293,272-120,344,603) harbors HADHA and HADHB, encoding the mitochondria trifunctional protein (MTP) alpha- and beta-subunits, respectively. The HADHA HADHB deficient mice had a decreased weight gain and cardiac arrhythmias.
In this study, we revealed the distribution of the unprecedented number of 565 CNVRs in 1,327 pigs from 18 diverse populations, and found that CNVRs were non-randomly distributed in the porcine genome. CNV numbers differ greatly among diverse populations. More than 72.87% of CNVRs were restricted to one population. The main functional categories of CNV-related genes were similar to those of in other mammals. With the QTL mapping data and the identified CNVRs in the White Duroc × Erhualian F2 intercross, and the description of phenotypes in knockout mice, we identified 7 copy number variable genes as potential candidate genes for porcine complex traits. These findings give novel insights into porcine CNVs and provide resources to facilitate the further identification of trait-related CNVs.
A total of 1,693 animals from 18 pig populations including 10 Chinese indigenous breeds, 2 Western commercial breeds, 1 wild boar and 5 F2 resource populations comprising White Duroc × Erhualian, Bamaxiang × Erhualian, Rongchang × Erhualian, Shazilin × Erhualian and Tongcheng × Erhualian were used in this study (Table1). Thereinto, 1,021 animals were obtained from the large scale White Duroc × Erhualian F2 resource population in which two White Duroc boars and 17 Erhualian sows were crossed as founder animals to produce F1 animals, and 59 F1 sows were randomly mated with 9 F1 boars to generate 1,912 F2 individuals. Total 422 traits related to growth, meat quality, body composition, blood physiological and biochemical parameters, reproduction and immune capacity were well phenotyped in this intercross and genome-wide QTL mapping was carried out for these traits[38–44]. All animal procedures were conducted according to the guidelines for the care and use of experimental animals established by the Ministry of Agriculture of China.
Genome-wide SNP genotyping
Genomic DNA was extracted from ear or spleen tissues with the routine phenol/chloroform extraction method. All 1,693 animals were genotyped with Porcine SNP60 BeadChip using the Infinium HD Assay Ultra protocol (Illumina Inc., San Diego, USA). The position of each SNP in the pig genome assembly (Sscrofa10.2) was determined by SOAP2 software. The quality control of genotypes was performed with GenABEL procedure in R. The SNPs in sex chromosomes and those not mapped or mapped to multi-positions in the Sscrofa10.2 assembly were discarded. A final set of 5,2596 SNPs on 18 autosomes with a unique position in Sscrofa10.2 was used for further analysis.
PennCNV was used to CNV calling. The software integrates a Hidden Markov Model (HMM) for high resolution copy number variation detection with whole-genome SNP genotyping data. The signal intensity data of log R Ratio (LRR) and B allele frequency (BAF) were exported from Illumina BeadStudio software. Individual-based CNV calling was performed using the default parameters of the HMM model by integrating Log R Ratio, BAF, population allele frequency and the SNP distance.
To reduce the false discover rate in CNV calling, we used a calling criteria requiring that the standard deviation (SD) of LRR must be under or less 0.35, the CNV contained three or more consecutive SNPs and the length of CNV region must be more than 50 kb at the calling-level. We set the “-qcnumcnv 50” argument in the command line to treat any samples with > 50 CNV calls as low quality samples and eliminated them from analysis. GC model file was used to adjust the signal intensity values for CNV calling. For F2 individuals, the “-trio” argument was employed in CNV calling to make use of the family information. The CNVs whose 50.0% of sequence overlapped with the telomere region and those detected in only one individual and not overlapped with any other discovered CNVs were also removed from analysis. All putative CNVs identified in this study were pooled across breeds. We aggregated the overlapping CNVs identified across all samples to determine CNVRs following the previously published protocols.
Quantitative real time PCR
Seven CNVRs were randomly selected for validation by quantitative real time PCR (qPCR). The 2-△△Ct method was used to estimate relative quantification (RQ) of CNVRs. This comparative method used a target assay for the CNV region and a reference assay of β-ACTIN as an internal control. The test and control primers were verified for their amplification efficiency. Six DNA samples were randomly selected including those with or without copy number variant for each CNVR. Primers and TaqMan probes labelled with FAM for each CNVR were designed with Allele 6.0 software (Applied Biosystems, Foster City, USA) and are listed in Additional file2: Table S3. qPCR was carried out in a total volume of 20 μl mixture containing 1 × Premix Ex Taq™ (TaKaRa, Dalian, China), 0.2 μmol/L each primer, 1 × ROX Reference Dye II and 100 ng genomic DNA in an ABI 7500 FAST instrument (Applied Biosystems, Foster City, USA). The thermal cycle parameter was: 30 sec at 95°C, and 40 cycles of 3 sec at 95°C and 30 sec at 60°C. Each sample was analyzed in triplicate. Results were analyzed with the ABI7500 software v2.0.5 (Applied Biosystems, Foster City, USA).
Gene identification and functional classification
Genes within porcine CNVRs were annotated by BioMart. These genes were tested for enrichment of molecule function, cell component and biological process in gene ontology (GO) terms in DAVID Bioinformatics Resources 6.7. Considering the limited number of pig genes assigned to GO terms, the human annotated genes that were homologous to pig genes were used as the background. Multiple tests were corrected by FDR corrections and enrichment threshold was set as EASE score of adjusted FDR P ≤ 0.05.
Association of the CNV268 with MCH and MCV was analyzed with a mixed linear model. Gender and batch were considered as fixed effects. The QTL for MCH and MCV was re-mapped with the CNV268 as the fixed effect. All the analyses were performed with R package.
Sex determining region Y-box 5
V-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog
Hydroxyacyl-CoA dehydrogenase, alpha subunit
Hydroxyacyl-CoA dehydrogenase, belta subunit
Kinase insert domain protein receptor
Acidic (leucine-rich) nuclear phosphoprotein 32 family, member B
Berardinelli-Seip congenital lipodystrophy 2 (seipin)
Latent transforming growth factor beta binding protein 3
Growth differentiation factor 3
Glycogen synthase 1
We thank the colleagues in Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China for their assistance in the management of resource population and sample collection. This work was supported by National Natural Science Foundation of China (31160225), National Key Basic Research Program of China (2012CB124702, 2012CB722502) and Program for Changjiang Scholars and Innovative Research Team in University (IRT1136).
- Stankiewicz P, Lupski JR: Structural variation in the human genome and its role in disease. Annu Rev Med. 2010, 61: 437-455. 10.1146/annurev-med-100708-204735.View ArticlePubMedGoogle Scholar
- Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, Cho EK, Dallaire S, Freeman JL, González JR, Gratacòs M, Huang J, Kalaitzopoulos D, Komura D, MacDonald JR, Marshall CR, Mei R, Montgomery L, Nishimura K, Okamura K, Shen F, Somerville MJ, Tchinda J, Valsesia A, Woodwark C, Yang F, Zhang J, Zerjal T, Zhang J, Armengol L, Conrad DF, Estivill X, Tyler-Smith C, Carter NP, Aburatani H, Lee C, Jones KW, Scherer SW, Hurles ME: Global variation in copy number in the human genome. Nature. 2006, 444: 444-454. 10.1038/nature05329.PubMed CentralView ArticlePubMedGoogle Scholar
- Marenne G, Rodríguez-Santiago B, Closas MG, Pérez-Jurado L, Rothman N, Rico D, Pita G, Pisano DG, Kogevinas M, Silverman DT, Valencia A, Real FX, Chanock SJ, Génin E, Malats N: Assessment of copy number variation using the Illumina Infinium 1 M SNP-array: a comparison of methodological approaches in the Spanish Bladder Cancer/EPICURO study. Hum Mutat. 2011, 32: 240-248. 10.1002/humu.21398.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Li M, Hadley D, Liu R, Glessner J, Grant SF, Hakonarson H, Bucan M: PennCNV: an integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007, 17: 1665-1674. 10.1101/gr.6861907.PubMed CentralView ArticlePubMedGoogle Scholar
- Colella S, Yau C, Taylor JM, Mirza G, Butler H, Clouston P, Bassett AS, Seller A, Holmes CC, Ragoussis J: QuantiSNP: an Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data. Nucleic Acids Res. 2007, 35: 2013-2025. 10.1093/nar/gkm076.PubMed CentralView ArticlePubMedGoogle Scholar
- Gai X, Perin JC, Murphy K, O'Hara R, D'arcy M, Wenocur A, Xie HM, Rappaport EF, Shaikh TH, White PS: CNV Workshop: an integrated platform for high-throughput copy number variation discovery and clinical diagnostics. BMC Bioinformatics. 2010, 11: 74-10.1186/1471-2105-11-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Winchester L, Yau C, Ragoussis J: Comparing CNV detection methods for SNP arrays. Brief Funct Genomic Proteomic. 2009, 8: 353-366. 10.1093/bfgp/elp017.View ArticlePubMedGoogle Scholar
- Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, et al: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464: 704-712. 10.1038/nature08516.PubMed CentralView ArticlePubMedGoogle Scholar
- Guryev V, Saar K, Adamovic T, Verheul M, van Heesch SA, Cook S, Pravenec M, Aitman T, Jacob H, Shull JD, Hubner N, Cuppen E: Distribution and functional impact of DNA copy number variation in the rat. Nat Genet. 2008, 40: 538-545. 10.1038/ng.141.View ArticlePubMedGoogle Scholar
- Nicholas TJ, Cheng Z, Ventura M, Mealey K, Eichler EE, Akey JM: The genomic architecture of segmental duplications and associated copy number variants in dogs. Genome Res. 2009, 19: 491-499.PubMed CentralView ArticlePubMedGoogle Scholar
- Hou Y, Liu GE, Bickhart DM, Cardone MF, Wang K, Kim ES, Matukumalli LK, Ventura M, Song J, VanRaden PM, Sonstegard TS, Van Tassell CP: Genomic characteristics of cattle copy number variations. BMC Genomics. 2011, 12: 127-10.1186/1471-2164-12-127.PubMed CentralView ArticlePubMedGoogle Scholar
- Doan R, Cohen N, Harrington J, Veazy K, Juras R, Cothran G, McCue ME, Skow L, Dindot SV: Identification of copy number variants in horses. Genome Res. 2012, 22: 899-907. 10.1101/gr.128991.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, Lee C, Tyler-Smith C, Carter N, Scherer SW, Tavaré S, Deloukas P, Hurles ME, Dermitzakis ET: Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007, 315: 848-853. 10.1126/science.1136678.PubMed CentralView ArticlePubMedGoogle Scholar
- Gamazon ER, Nicolae DL, Cox NJ: A study of CNVs as trait-associated polymorphisms and as expression quantitative trait loci. PLoS Genet. 2011, 7: e1001292-10.1371/journal.pgen.1001292.PubMed CentralView ArticlePubMedGoogle Scholar
- Tam GW, Redon R, Carter NP, Grant SG: The role of DNA copy number variation in schizophrenia. Biol Psychiatry. 2009, 66: 1005-1012. 10.1016/j.biopsych.2009.07.027.View ArticlePubMedGoogle Scholar
- Sha BY, Yang TL, Zhao LJ, Chen XD, Guo Y, Chen Y, Pan F, Zhang ZX, Dong SS, Xu XH, Deng HW: genome-wide association study suggested copy number variation may be associated with body mass index in the Chinese population. J Hum Genet. 2009, 54: 199-202. 10.1038/jhg.2009.10.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaschl H, Aitman TJ, Vyse TJ: Copy number variation in the human genome and its implication in autoimmunity. Clin Exp Immunol. 2009, 156: 12-16. 10.1111/j.1365-2249.2008.03865.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper GM, Nickerson DA, Eichler EE: Mutational and selective effects on copy-number variants in the human genome. Nat Genet. 2007, 39: S22-S29. 10.1038/ng2054.View ArticlePubMedGoogle Scholar
- Wright D, Boije H, Meadows JR, Bed'hom B, Gourichon D, Vieaud A, Tixier-Boichard M, Rubin CJ, Imsland F, Hallböök F, Andersson L: Copy number variation in intron 1 of SOX5 causes the Pea-comb phenotype in chickens. PLoS Genet. 2009, 5: e1000512-10.1371/journal.pgen.1000512.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosengren Pielberg G, Golovko A, Sundstrom E, Curik I, Lennartsson J, Seltenhammer MH, Druml T, Binns M, Fitzsimmons C, Lindgren G, Sandberg K, Baumung R, Vetterlein M, Strömberg S, Grabherr M, Wade C, Lindblad-Toh K, Pontén F, Heldin CH, Sölkner J, Andersson L: A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse. Nat Genet. 2008, 40: 1004-1009. 10.1038/ng.185.View ArticlePubMedGoogle Scholar
- Salmon Hillbertz NH, Isaksson M, Karlsson EK, Hellmén E, Pielberg GR, Savolainen P, Wade CM, von Euler H, Gustafson U, Hedhammar A, Nilsson M, Lindblad-Toh K, Andersson L, Andersson G: Duplication of FGF3, FGF4, FGF19 and ORAOV1 causes hair ridge and predisposition to dermoid sinus in Ridgeback dogs. Nat Genet. 2007, 39: 1318-1320. 10.1038/ng.2007.4.View ArticlePubMedGoogle Scholar
- Fontanesi L, Beretti F, Riggio V, Gómez González E, Dall'Olio S, Davoli R, Russo V, Portolano B: Copy number variation and missense mutations of the agouti signaling protein (ASIP) gene in goat breeds with different coat colors. Cytogenet Genome Res. 2009, 126: 333-347. 10.1159/000268089.View ArticlePubMedGoogle Scholar
- Johansson A, Pielberg G, Andersson L, Edfors-Lilja I: Polymorphism at the porcine Dominant white/KIT locus influence coat colour and peripheral blood cell measures. Anim Genet. 2005, 36: 288-296. 10.1111/j.1365-2052.2005.01320.x.View ArticlePubMedGoogle Scholar
- Pielberg G, Olsson C, Syvänen AC, Andersson L: Unexpectedly high allelic diversity at the KIT locus causing dominant white color in the domestic pig. Genetics. 2002, 160: 305-311.PubMed CentralPubMedGoogle Scholar
- Fadista J, Nygaard M, Holm LE, Thomsen B, Bendixen C: A snapshot of CNVs in the pig genome. PLoS One. 2008, 3: e3916-10.1371/journal.pone.0003916.PubMed CentralView ArticlePubMedGoogle Scholar
- Tang H, Li F, Finlayson HA, Smith S, Lu Z, Langford C, Archibald AL: Structural and copy number variation in the pig genome. Proceedings of the Plant & Animal Genomes XVIII Conference: 9–13 January 2010; San Diego [abstract]. Edited by: Gerard L, Vicoria Carollo B, David G. 2010, 609-Google Scholar
- Ramayo-Caldas Y, Castello A, Pena RN, Alves E, Mercade A, Souza CA, Fernandez AI, Perez-Enciso M, Folch JM: Copy number variation in the porcine genome inferred from a 60 k SNP BeadChip. BMC Genomics. 2010, 11: 593-10.1186/1471-2164-11-593.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo Y, Mao H, Ren J, Yan X, Duan Y, Yang G, Ren D, Zhang Z, Yang B, Ouyang J, Brenig B, Haley C, Huang L: A linkage map of the porcine genome from a large-scale White Duroc x Erhualian resource population and evaluation of factors affecting recombination rates. Anim Genet. 2009, 40: 47-52. 10.1111/j.1365-2052.2008.01802.x.View ArticlePubMedGoogle Scholar
- Ramos AM, Crooijmans RPMA, 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 TPL, Schnabel RD, Van Tassell CP, Taylor JF, Wiedmann RT, Schook LB, Groenen MAM: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Månér S, Massa H, Walker M, Chi M, Navin N, Lucito R, Healy J, Hicks J, Ye K, Reiner A, Gilliam TC, Trask B, Patterson N, Zetterberg A, Wigler M: Large-scale copy number polymorphism in the human genome. Science. 2004, 305: 525-528. 10.1126/science.1098918.View ArticlePubMedGoogle Scholar
- Chen W, Hayward C, Wright AF, Hicks AA, Vitart V, Knott S, Wild SH, Pramstaller PP, Wilson JF, Rudan I, Porteous DJ: Copy number variation across European populations. PLoS One. 2011, 6: e23087-10.1371/journal.pone.0023087.PubMed CentralView ArticlePubMedGoogle Scholar
- Database of Genomic Variants.http://projects.tcag.ca/variation/,
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Young JM, Endicott RM, Parghi SS, Walker M, Kidd JM, Trask BJ: Extensive copy-number variation of the human olfactory receptor gene family. Am J Hum Genet. 2008, 83: 228-242. 10.1016/j.ajhg.2008.07.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Hussain A, Saraiva LR, Korsching SI: Positive Darwinian selection and the birth of an olfactory receptor clade in teleosts. PNAS. 2009, 106: 4313-4318. 10.1073/pnas.0803229106.PubMed CentralView ArticlePubMedGoogle Scholar
- Hediger MA, Romero MF, Peng JB, Rolfs A, Takanaga H, Bruford EA: The ABCs of solute carriers: physiological, pathological and therapeutic implications of human membrane transport proteins: Introduction. Pflugers Arch. 2004, 447: 465-468. 10.1007/s00424-003-1192-y.View ArticlePubMedGoogle Scholar
- Chen R, Ren J, Li W, Huang X, Yan X, Yang B, Zhao Y, Guo Y, Mao H, Huang L: A genome-wide scan for quantitative trait loci affecting serum glucose and lipids in a White Duroc x Erhualian intercross F(2) population. Mamm Genome. 2009, 20: 386-392. 10.1007/s00335-009-9190-9.View ArticlePubMedGoogle Scholar
- Ma J, Ren J, Guo Y, Duan Y, Ding N, Zhou L, Li L, Yan X, Yang K, Huang L, Song Y, Xie J, Milan D, Huang L: Genome-wide identification of quantitative trait loci for carcass composition and meat quality in a large-scale White Duroc x Chinese Erhualian resource population. Anim Genet. 2009, 40: 637-647. 10.1111/j.1365-2052.2009.01892.x.View ArticlePubMedGoogle Scholar
- Ai H, Ren J, Zhang Z, Ma J, Guo Y, Yang B, Huang L: Detection of quantitative trait loci for growth- and fatness-related traits in a large-scale White Duroc × Erhualian intercross pig population. Anim Genet. 2012, 43: 383-391. 10.1111/j.1365-2052.2011.02282.x.View ArticlePubMedGoogle Scholar
- Ren DR, Ren J, Xing YY, Guo YM, Wu YB, Yang GC, Mao HR, Huang LS: A genome scan for quantitative trait loci affecting male reproductive traits in a White Duroc x Chinese Erhualian resource population. J Anim Sci. 2009, 87: 17-23.View ArticlePubMedGoogle Scholar
- Li K, Ren J, Xing Y, Zhang Z, Ma J, Guo Y, Huang L: Quantitative trait loci for litter size and prenatal loss in a White Duroc x Chinese Erhualian resource population. Anim Genet. 2009, 40: 963-966. 10.1111/j.1365-2052.2009.01931.x.View ArticlePubMedGoogle Scholar
- Zou Z, Ren J, Yan X, Huang X, Yang S, Zhang Z, Yang B, Li W, Huang L: Quantitative trait loci for porcine baseline erythroid traits at three growth ages in a White Duroc x Erhualian F(2) resource population. Mamm Genome. 2008, 19: 640-646. 10.1007/s00335-008-9142-9.View ArticlePubMedGoogle Scholar
- Mao H, Guo Y, Yang G, Yang B, Ren J, Liu S, Ai H, Ma J, Brenig B, Huang L: A genome-wide scan for quantitative trait loci affecting limb bone lengths and areal bone mineral density of the distal femur in a White Duroc x Erhualian F2 population. BMC Genet. 2008, 9: 63-PubMed CentralView ArticlePubMedGoogle Scholar
- Pig QTL database.http://www.animalgenome.org/cgi-bin/QTLdb/SS/index,
- Reilly PT, Afzal S, Gorrini C, Lui K, Bukhman YV, Wakeham A, Haight J, Ling TW, Cheung CC, Elia AJ, Turner PV, Mak TW: Acidic nuclear phosphoprotein 32 kDa (ANP32) B-deficient mouse reveals a hierarchy of ANP32 importance in mammalian development. Proc Natl Acad Sci USA. 2011, 108: 10243-10248. 10.1073/pnas.1106211108.PubMed CentralView ArticlePubMedGoogle Scholar
- Cui X, Wang Y, Tang Y, Liu Y, Zhao L, Deng J, Xu G, Peng X, Ju S, Liu G, Yang H: Seipin ablation in mice results in severe generalized lipodystrophy. Hum Mol Genet. 2011, 20: 3022-3030. 10.1093/hmg/ddr205.View ArticlePubMedGoogle Scholar
- Dabovic B, Chen Y, Colarossi C, Obata H, Zambuto L, Perle MA, Rifkin DB: Bone abnormalities in latent TGF-[beta] binding protein (Ltbp)-3-null mice indicate a role for Ltbp-3 in modulating TGF-[beta] bioavailability. J Cell Biol. 2002, 156: 227-232. 10.1083/jcb.200111080.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen JJ, Huang L, Li L, Jorgez C, Matzuk MM, Brown CW: Deficiency of growth differentiation factor 3 protects against diet-induced obesity by selectively acting on white adipose. Mol Endocrinol. 2009, 23: 113-123.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen C, Ware SM, Sato A, Houston-Hawkins DE, Habas R, Matzuk MM, Shen MM, Brown CW: The Vg1-related protein Gdf3 acts in a Nodal signaling pathway in the pre-gastrulation mouse embryo. Development. 2006, 133: 319-329.View ArticlePubMedGoogle Scholar
- Bouskila M, Hunter RW, Ibrahim AF, Delattre L, Peggie M, van Diepen JA, Voshol PJ, Jensen J, Sakamoto K: Allosteric regulation of glycogen synthase controls glycogen synthesis in muscle. Cell Metab. 2010, 12: 456-466. 10.1016/j.cmet.2010.10.006.View ArticlePubMedGoogle Scholar
- Magnol L, Chevallier MC, Nalesso V, Retif S, Fuchs H, Klempt M, Pereira P, Riottot M, Andrzejewski S, Doan BT, Panthier JJ, Puech A, Beloeil JC, de Angelis MH, Hérault Y: KIT is required for hepatic function during mouse post-natal development. BMC Dev Biol. 2007, 7: 81-10.1186/1471-213X-7-81.PubMed CentralView ArticlePubMedGoogle Scholar
- Waskow C, Paul S, Haller C, Gassmann M, Rodewald HR: Viable c-Kit(W/W) mutants reveal pivotal role for c-kit in the maintenance of lymphopoiesis. Immunity. 2002, 17: 277-288. 10.1016/S1074-7613(02)00386-2.View ArticlePubMedGoogle Scholar
- Tsujimura T, Koshimizu U, Katoh H, Isozaki K, Kanakura Y, Tono T, Adachi S, Kasugai T, Tei H, Nishimune Y: Mast cell number in the skin of heterozygotes reflects the molecular nature of c-kit mutation. Blood. 1993, 81: 2530-2538.PubMedGoogle Scholar
- Rubin J, Schwartz Z, Boyan BD, Fan X, Case N, Sen B, Drab M, Smith D, Aleman M, Wong KL, Yao H, Jo H, Gross TS: Caveolin-1 knockout mice have increased bone size and stiffness. J Bone Miner Res. 2007, 22: 1408-1418. 10.1359/jbmr.070601.View ArticlePubMedGoogle Scholar
- Johansson Moller M, Chaudhary R, Hellmén E, Höyheim B, Chowdhary B, Andersson L: Pigs with the dominant white coat color phenotype carry a duplication of the KIT gene encoding the mast/stem cell growth factor receptor. Mamm Genome. 1996, 7: 822-830. 10.1007/s003359900244.View ArticlePubMedGoogle Scholar
- Xu GL, Ren J, Ding NS, Ai HS, Guo YM, Chen CY, Huang LS: Genetic analysis of the KIT and MC1R genes in Chinese indigenous pigs with belt-like coat color phenotypes. Anim Genet. 2006, 37: 518-519. 10.1111/j.1365-2052.2006.01504.x.View ArticlePubMedGoogle Scholar
- Lai F, Ren J, Ai H, Ding N, Ma J, Zeng D, Chen C, Guo Y, Huang L: Chinese white Rongchang pig does not have the dominant white allele of KIT but has the dominant black allele of MC1R. J Hered. 2007, 98: 84-87.View ArticlePubMedGoogle Scholar
- Meijerink E, Fries R, Vogeli P, Masabanda J, Wigger G, Stricker C, Neuenschwander S, Bertschinger HU, Stranzinger G: Two a(1,2)fucosyltransferase genes on porcine chromosome 6q11 are closely linked to the blood group inhibitor (S) and Escherichia coli F18 receptor (ECF18R) loci. Mamm Genome. 1997, 8: 736-741. 10.1007/s003359900556.View ArticlePubMedGoogle Scholar
- Horák P, Urban T, Dvorák J: The FUT1 and ESR genes–their variability and associations with reproduction in Prestice Black-Pied sows. J Anim Breed Genet. 2005, 122: 210-213. 10.1111/j.1439-0388.2005.00502.x.View ArticlePubMedGoogle Scholar
- Ken-Dror G, Talmud PJ, Humphries SE, Drenos F: APOE/C1/C4/C2 gene cluster genotypes, haplotypes and lipid levels in prospective coronary heart disease risk among UK healthy men. Mol Med. 2010, 16: 389-399.PubMed CentralView ArticlePubMedGoogle Scholar
- Kao HJ, Cheng CF, Chen YH, Hung SI, Huang CC, Millington D, Kikuchi T, Wu JY, Chen YT: ENU mutagenesis identifies mice with cardiac fibrosis and hepatic steatosis caused by a mutation in the mitochondrial trifunctional protein beta-subunit. Hum Mol Genet. 2006, 15: 3569-3577. 10.1093/hmg/ddl433.View ArticlePubMedGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25: 1966-1967. 10.1093/bioinformatics/btp336.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- DAVID Bioinformatics Resources 6.7.http://david.abcc.ncifcrf.gov/,
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.