Identification of genome-wide copy number variations among diverse pig breeds by array CGH
- Yan Li†1,
- Shuqi Mei†2,
- Xuying Zhang1,
- Xianwen Peng2,
- Gang Liu1,
- Hu Tao1,
- Huayu Wu2,
- Siwen Jiang1,
- Yuanzhu Xiong1 and
- Fenge Li1Email author
© Li et al.; licensee BioMed Central Ltd. 2012
Received: 17 August 2012
Accepted: 19 December 2012
Published: 24 December 2012
Recent studies have shown that copy number variation (CNV) in mammalian genomes contributes to phenotypic diversity, including health and disease status. In domestic pigs, CNV has been catalogued by several reports, but the extent of CNV and the phenotypic effects are far from clear. The goal of this study was to identify CNV regions (CNVRs) in pigs based on array comparative genome hybridization (aCGH).
Here a custom-made tiling oligo-nucleotide array was used with a median probe spacing of 2506 bp for screening 12 pigs including 3 Chinese native pigs (one Chinese Erhualian, one Tongcheng and one Yangxin pig), 5 European pigs (one Large White, one Pietrain, one White Duroc and two Landrace pigs), 2 synthetic pigs (Chinese new line DIV pigs) and 2 crossbred pigs (Landrace × DIV pigs) with a Duroc pig as the reference. Two hundred and fifty-nine CNVRs across chromosomes 1–18 and X were identified, with an average size of 65.07 kb and a median size of 98.74 kb, covering 16.85 Mb or 0.74% of the whole genome. Concerning copy number status, 93 (35.91%) CNVRs were called as gains, 140 (54.05%) were called as losses and the remaining 26 (10.04%) were called as both gains and losses. Of all detected CNVRs, 171 (66.02%) and 34 (13.13%) CNVRs directly overlapped with Sus scrofa duplicated sequences and pig QTLs, respectively. The CNVRs encompassed 372 full length Ensembl transcripts. Two CNVRs identified by aCGH were validated using real-time quantitative PCR (qPCR).
Using 720 K array CGH (aCGH) we described a map of porcine CNVs which facilitated the identification of structural variations for important phenotypes and the assessment of the genetic diversity of pigs.
Genetic and archaeological findings suggest that pig domestication began about 9000–10000 years before present (YBP) at multiple sites across Eurasia, followed by their subsequent spread at a worldwide scale . Historically, Europe and China are two major areas of pig breeding . Over the past centuries, pigs have shown marked differences between these two areas, even if many European pig breeds carry far Eastern haplotypes at high frequencies because of an ancient introgression with Chinese swine . The Chinese pigs differ significantly from European pig breeds such as the Large White for many traits including fatness and ear traits [3–5]. Genetic variation within the gene pool which produce the above different phenotypes are selected for or against by evolution. Microsatellites, single nucleotide polymorphisms (SNPs) were the main measures of genetic variations in pigs, producing a USMARC pig SNP map (http://www.marc.usda.gov/genome/swine/marker_list.html) and the PorcineSNP60 Genotyping BeadChip with 62163 SNP probes . Recently, structural variations including insertions, duplications, deletions, inversions and translocations of DNA have been shown to contribute to the major phenotypic variations . Copy number variation (CNV) is described as a segment of DNA >1 kb that is copy number variable when compared with a reference genome . This variation may either be inherited or caused by de novo mutation [9–12]. It has become apparent that CNVs are genome-wide present in the human genome  and the genome of farm animals including cattle [13–16], avian [17–19], sheep , goat . About a range from 5% to 16% of the human genome was covered by CNVs [22, 23]. CNVs can lead to striking phenotypic consequences as a result of altering gene dosage, disrupting coding sequences, or perturbing long-range gene regulation by position effects [24–26]. These striking phenotypic consequences include some common complex diseases such as autism , schizophrenia , auto-immune Addison's disease .
Recently many efforts have been used to detect pig CNVs. By a custom-made tiling oligonucleotide array, 37 CNV regions (CNVRs) across chromosomes 4, 7, 14, and 17 were identified in 12 unrelated Duroc boars . Comparative genome hybridization (CGH) array was also conducted for chromosomes 7 and 8 in 9 different pig populations including Duroc, Large White, Meishan, Pietrain, Hampshire and Wild Boar . By analyzing data from the Porcine SNP60 BeadChip, 49 CNVRs were identified in 55 animals from an Iberian × Landrace cross (IBMAP)  and 382 CNVRs were identified from three purebred populations (Yorkshire, Landrace and Songliao Black) and one Duroc × Erhualian crossbred population . Up until now, few studies have confirmed the genome-wide presence of CNVs in pigs using array CGH (aCGH) with high-density probes. Here we reported the use of high-resolution oligonucleotide aCGH to identify the CNV regions in 12 individual pigs from different pig populations. This analysis provided a high-resolution map of copy number variations in the pig genome with a median probe spacing of 2506 bp relative to the latest porcine genome assembly (Sscrofa9.2).
Results and discussion
The overview of CNVR library
Using the custom tiling oligonucleotide aCGH approach, Fadista et al.  addressed 37 CNVRs on the Sus scrofa chromosomes (SSCs) 4, 7, 14, and 17 of the preliminary assembly of pig genome among 12 Duroc boars. Ramayo-Caldas et al.  detected 49 CNVRs using the Porcine SNP60 BeadChip data of 55 animals from an Iberian × Landrace cross. Wang et al.  detected 382 CNVRs based on the Porcine SNP60 genotyping data of 474 pigs. Two of the 37 CNVRs (5.41%) detected by Fadista et al. , 8 of the 49 CNVRs (16.32%) detected by Ramayo-Caldas et al. , 24 of the 382 CNVRs (6.28%) detected by Wang et al.  were identical or overlapped with the detected CNVRs in this study (Additional file 2). Totally 39 of the presently detected 259 CNVRs (15.06%) were identical or overlapped with those previously reported pig CNVRs (Additional file 2). The main potential reasons for this less well-overlapping result could be the different genetic backgrounds of pig samples, different platforms and various calling algorithms between the present study and other studies.
Compared with PorcineSNP60 Genotyping BeadChip, the detection power of 720 K aCGH was enhanced by dense marker density, uniform distribution of probes along each chromosome [6, 30]. Hence, some small CNVRs can be detected by aCGH technique, as the minimum CNV lengths were 2.30 kb in our present study, and 2.08 kb in the study of Fadista et al. , whereas the minimum CNV length detected by SNP chip were 5.03 kb and 44.65 kb, respectively [30, 31].
CNVRs chromosome distribution and status
Chromosome distribution of CNVRs in pigs
CNVR Size (bp)
Chr. Size (bp)
Concerning copy number status, 93 (35.91%) CNVRs were called as gains, 140 (54.05%) were called as losses and the remaining 26 (10.04%) were called as both gains and losses. Previously, it has been suggested that deletions are under stronger purifying selection than duplications . If so, deletions should be both less frequent and shorter than duplications . However, when we compared the length of gains with losses in the CNVRs, loss regions had slightly larger sizes than gain regions with the average length of 57.39 kb and 45.86 kb respectively (T-test not statistically significant at p value > 0.05). The possible reason was that the aCGH approach might favor the identification of deletions [14, 15, 21, 28]. As the samples were collected from 9 different populations, the considerable number of CNVRs status displaying in ‘both gains and losses’ might be due to the different genetic origins.
Putative population-specific CNVRs and cluster analysis
Some putative population-specific CNVRs were detected. For example, 6 CNVRs including CNVR_132 were purebred Landrace-specific, and CNVR_145 were purebred DIV-specific. CNVR_100 including KIT gene contained amplifications specifically in 8 pigs with dominant white color and a Pietrain pig with black spots, and CNVR_251 contained gains in pigs without dominant white color such as Yangxin, Erhualian, Tongcheng and Pietrain pigs. However, due to the limited samples used in the present study, the putative population-specific CNVRs need future study. And we also found 3 de novo CNVRs, of which CNVR_IDs 36, 149 were present in 2 Landrace × DIV crossbred pigs but not in their parents, while CNVR_259 were absent in 2 Landrace × DIV crossbred pigs but present in their parents.
Duplicated sequences colocalize with CNVRs in the pig genome
Although the exact interpretation of mechanisms responsible for generating CNVs is still unclear, previous studies have noted a four- to twenty-fold enrichment of CNVs near SDs in the other mammalian genomes [22, 32, 36]. Duplicated sequences are typical segments of DNA which range in size from one to hundreds of kb, share a high level of sequence identity (≥ 90%) and occur at more than one site within the genome . Under the same filter criterion, about 66.02% (171/259) of CNV regions directly overlapping with Sus scrofa duplicated sequences were identified through blasting the CNVR sequence against the Ensembl pig genomic sequences. As our present BLAST results did not retain a CNVR overlapping with a duplicated sequence by less than 1000 bp, so the overlaps of CNVs and their targeted duplicated sequences were under reporting. There were 13.5–25.0% CNVRs mapped to duplicated sequences in the previous reports [28, 37]. The difference may be related to differences in samples. CNVRs overlapping duplicated sequences were significantly different in average size (87.12 kb versus 22.23 kb, t-test p < 0.01) with the CNVRs that did not overlap duplicated sequences, consistent with previous CNV studies reporting a stronger association between duplicated sequences and long CNVRs [9, 11].
Gene contents of pig CNV regions
Enriched GO terms and KEGG pathway associated with the CNV regions (Modified Fisher Exact P-value ≤ 0.05)
Involved genes/ total genes (%)
GO:0007608 ~ sensory perception of smell
GO:0007606 ~ sensory perception of chemical stimulus
GO:0007600 ~ sensory perception
GO:0050890 ~ cognition
GO:0007186 ~ G-protein coupled receptor protein signaling pathway
GO:0050877 ~ neurological system process
GO:0007166 ~ cell surface receptor linked signal transduction
GO:0004984 ~ olfactory receptor activity
GO:0047961 ~ glycine N-acyltransferase activity
GO:0016021 ~ integral to membrane
GO:0031224 ~ intrinsic to membrane
Evolutionary rates of pig monomorphic and CNVR genes compared with human and mouse
Pig CNVRs overlapped with QTL regions
We queried the animal QTL database that held publicly available QTL data on livestock species. Retrieving all the porcine QTLs (http://www.animalgenome.org/cgi-bin/QTLdb/SS/download?file=gbpSS_10.2) within 2 Mb of our CNVRs resulted that 34 CNVRs overlapped with QTLs for several important traits including average daily gain (ADG) (Additional file 4). However, as the pig QTLs are not fully defined, the contribution of these QTL-overlapping CNVRs to complex traits needs further study.
Validation of CNVRs by real-time quantitative (qPCR)
qPCR was performed to validate 2 CNVRs (CNVR_IDs 100 and 215) detected by the aCGH experiment. Thirteen DNA samples including the reference used in aCGH were used for qPCR analysis. CNVR_100 and CNVR_215 were validated (Additional file 5) with the p threshold values 0.05 as the previous reports .
In summary, we described a map of porcine CNVs between breeds by a high-resolution array CGH, which was confirmed to be a very valid method to detect porcine genome-wide CNVs. With a stringent CNV calling criterion, 259 highly reliable CNV regions were reported here among diverse pig breeds. Future studies are required to assess the function of CNVs on pig important phenotypes. Our results facilitated the identification of structural variations for important phenotypes and the assessment of the genetic diversity in pigs.
All animal procedures were performed according to protocols approved by the Biological Studies Animal Care and Use Committee of Hubei Province, PR China. Twelve pigs including one White Duroc pig (♀), one Chinese Yangxin pig (♂), one Chinese Erhualian pig (♀), one Chinese Tongcheng pig (♀), one Large White pig (♀), one Pietrain pig (♂), two Landrace pigs (♂), two DIV pigs (♀) and two Landrace × DIV pigs (♀, ♂) were selected to function as test animals. Chinese Erhualian pigs were a strain of Chinese Taihu pig breed. Synthetic Line DIV was a result of cross of Landrace, Large White, Tongcheng or Taihu pigs. An unrelated female Duroc pig was selected as the common reference. The genomic DNA of 13 pig samples was extracted and purified from semen, whole blood or ear notch.
A 3 × 720 K whole genome tiling aCGH (NCBI GEO accession no. GPL16165) was designed (NimbleGen Systems, http://www.nimblegen.com) from the Sscrofa9.2 release (http://www.sanger.ac.uk/Projects/S_scrofa/), which was the new release at the time of the experiment. The probe design fundamentals were described in the NimbleGen technical note (http://www.nimblegen.com/products/lit/probe_design_2008_06_04.pdf). The probes with length of 50–60 bp were integrated into an array design using ArrayScribeTM, which resulted in a design with a median probe spacing of 2506 bp. Test DNA and reference DNA samples were independently labeled with either Cy3 or Cy5 dyes. Labeled DNA was co-hybridized to the custom-made NimbleGen CGH array (3 × 720 K). The array format included 3 arrays on single slides containing 719,336 probes. The arrays were scanned using a 5 μm scanner, and NimbleScan software (Roche NimbleGen) was used to retrieve fluorescent intensity raw data from the scanned images of the oligonucleotide tiling arrays. For each spot on the array, log2 ratios of the Cy3-labeled test sample versus Cy5-labeled reference sample were computed. Before normalization and segmentation analysis, spatial correction was applied. Specifically, locally weighted polynomial regression (LOESS) was used to adjust signal intensities based on X, Y feature position . Normalization was then performed using the q-spline method followed by segmentation using the CNV calling algorithm segMNT included in NimbleScan software . CNVRs were called as the segments with at least 5 consecutive probes, a mean |log2 ratio| of >0.50 and detected in two or more animals . Since the CNV calling pipeline requires at least 5 consecutive probes, our theoretical resolution for CNV detection is 10299 bp (median spacing × 4 + median oligo length × 5). As females had two copies of X-linked genes and males only had one copy, male–female aCGH resulted in an excess of female signals for X-linked genes that can be used to calibrate the threshold values and detection methods . aCGH data have been submitted to the GenBank gene expression omnibus database under the accession number GSE41488. The dendrogram were generated by average linkage clustering algorithm of Cluster 3.0 software .
In order to check if the CNVRs overlapped any duplicated sequence, BLAST was used to query the CNVRs sequences against the Sus scrofa genome sequence (Sscrofa9.2). Sequences were retained as duplicated sequences if they had ≥ 1 kb and ≥ 90% identity and occurred at more than one site within the genome.
Gene contents in the identified CNVRs were retrieved from the Sscrofa9.2 assembly using the BioMart (http://www.biomart.org/) . Gene content of pig CNV regions was assessed using Ensembl transcripts. The DAVID functional annotation tool (http://david.abcc.ncifcrf.gov/) was used to perform GO classification and KEGG pathway annotation of CNV mRNAs. Functional annotation terms from the ontologies of "biological processes", "molecular function" and "cellular component" were recorded. Since only a limited number of genes in the pig genome have been annotated, we converted the pig Ensembl transcripts IDs to orthologous mouse and human Ensembl gene IDs by BioMart, then carried out the GO and pathway analyses, as described previously .
All the porcine QTLs data were downloaded from pig QTL database (http://www.animalgenome.org/cgi-bin/QTLdb/SS/download?file=gbpSS_10.2) . The CNVRs were considered to be overlapping pig QTLs if they were within 2 Mb of pig QTLs .
Validation of CNVRs by qPCR
Determination of CNVRs by qPCR was performed using the Roche LightCycler® 480 Detection System and obtained the crossing thresholds (Ct) value following the guidelines of the manufacturer. The primers were designed using the Primer Premier 5 software and were available in the Additional file 6. As previously reported , the copy number of each CNVR was normalized against the Col10 region, a control region in the genome that did not vary in copy number between the pigs. Triplicate wells of reactions (15 μL) contained 7.5 μL SYBR Green Real-time PCR Master Mix, 1 μL of 10–20 ng/μL gDNA, 0.3 μL 5 μM of each primer and 0.1 μL ROX. The cycling conditions consisted of 1 cycle at 95°C for 10 min, followed by 40 cycles at 94°C for 20 sec, 60°C for 20 sec, and 72°C for 20 sec, with fluorescence acquisition at 74°C in single mode. The specific PCR products were confirmed by the results of melting curve analysis and agarose gel electrophoresis. Analysis of resultant crossing thresholds (Ct) was performed using the -ΔΔCt method .
Copy number variation
Polymerase chain reaction
Comparative genome hybridization
Real-time quantitative PCR
Relative quantification value
Quantitative trait locus
Tyrosine-protein kinase Kit
Cytochrome P450 gene family
Single nuclotide polymorphism
Homo sapiens chromosome
Sus scrofa chromosome
The database for annotation, visualization and integrated discovery
kyoto encyclopedia of genes and genomes
locally weighted polynomial regression
We thank the anonymous reviewers for critical reading and discussions of the manuscript. We are grateful to Prof. Alan Archibald (The Roslin Institute) for the suggestions for this study, and to CapitalBio Corporation for the technical assistance with NimbleGen CGH analysis. The authors also acknowledge the farmers for providing pig samples.
- Amills M, Clop A, Ramı′rez O, Pe′rez-Enciso M: Origin and genetic diversity of pig breeds. Encyclopedia of Life Sciences (ELS). 2010, John Wiley & Sons, Ltd, ChichesterGoogle Scholar
- Megens HJ, Crooijmans RP, San Cristobal M, Hui X, Li N, Groenen MA: Biodiversity of pig breeds from China and Europe estimated from pooled DNA samples: differences in microsatellite variation between two areas of domestication. Genet Sel Evol. 2008, 40: 103-128.PubMed CentralPubMedGoogle Scholar
- Haley CS, Agaro E, Ellis M: Genetic components of growth and ultrasonic fat depth traits in Meishan and Large White pigs and their reciprocal crosses. Anim Prod. 1992, 54: 105-115. 10.1017/S0003356100020626.View ArticleGoogle Scholar
- Haley CS, Lee GJ, Ritchie M: Comparative reproductive-performance in Meishan and Large White pigs and threir crosses. Anim Sci. 1995, 60: 259-267. 10.1017/S1357729800008420.View ArticleGoogle Scholar
- Wei WH, de Koning DJ, Penman JC, Finlayson HA, Archibald AL, Haley CS: QTL modulating ear size and erectness in pigs. Anim Genet. 2007, 38: 222-226. 10.1111/j.1365-2052.2007.01591.x.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Sebat J: Major changes in our DNA lead to major changes in our thinking. Nat Genet. 2007, 39: S3-S5. 10.1038/ng2095.View ArticlePubMedGoogle Scholar
- Feuk L, Carson AR, Scherer SW: Structural variation in the human genome. Nat Rev Genet. 2006, 7: 85-97.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
- Greenway SC, Pereira AC, Lin JC, DePalma SR, Israel SJ, Mesquita SM, Ergul E, Conta JH, Korn JM, McCarroll SA, Gorham JM, Gabriel S, Altshuler DM, Quintanilla-Dieck Mde L, Artunduaga MA, Eavey RD, Plenge RM, Shadick NA, Weinblatt ME, De Jager PL, Hafler DA, Breitbart RE, Seidman JG, Seidman CE: De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot. Nat Genet. 2009, 41: 931-935. 10.1038/ng.415.PubMed CentralView ArticlePubMedGoogle Scholar
- Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, Yamrom B, Yoon S, Krasnitz A, Kendall J, Leotta A, Pai D, Zhang R, Lee YH, Hicks J, Spence SJ, Lee AT, Puura K, Lehtimäki T, Ledbetter D, Gregersen PK, Bregman J, Sutcliffe JS, Jobanputra V, Chung W, Warburton D, King MC, Skuse D, Geschwind DH, Gilliam TC, Ye K, Wigler M: Strong association of De Novo copy number mutations with Autism. Science. 2007, 316: 445-449. 10.1126/science.1138659.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu B, Roos JL, Levy S, van Rensburg EJ, Gogos JA, Karayiorgou M: Strong association of de novo copy number mutations with sporadic schizophrenia. Nat Genet. 2008, 40: 880-885. 10.1038/ng.162.View ArticlePubMedGoogle Scholar
- Bae JS, Cheong HS, Kim LH, NamGung S, Park TJ, Chun JY, Kim JY, Pasaje CF, Lee JS, Shin HD: Identification of copy number variations and common deletion polymorphisms in cattle. BMC Genomics. 2010, 11: 232-10.1186/1471-2164-11-232.PubMed CentralView ArticlePubMedGoogle Scholar
- Fadista J, Thomsen B, Holm LE, Bendixen C: Copy number variation in the bovine genome. BMC Genomics. 2010, 11: 284-10.1186/1471-2164-11-284.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu GE, Hou Y, Zhu B, Cardone MF, Jiang L, Cellamare A, Mitra A, Alexander LJ, Coutinho LL, Dell'Aquila ME, Gasbarre LC, Lacalandra G, Li RW, Matukumalli LK, Nonneman D, Regitano LC, Smith TP, Song J, Sonstegard TS, Van Tassell CP, Ventura M, Eichler EE, McDaneld TG, Keele JW: Analysis of copy number variations among diverse cattle breeds. Genome Res. 2010, 20: 693-703. 10.1101/gr.105403.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu GE, Van Tassel CP, Sonstegard TS, Li RW, Alexander LJ, Keele JW, Matukumalli LK, Smith TP, Gasbarre LC: Detection of germline and somatic copy number variations in cattle. Dev Biol (Basel). 2008, 132: 231-237.Google Scholar
- Griffin DK, Robertson LB, Tempest HG, Vignal A, Fillon V, Crooijmans RP, Groenen MA, Deryusheva S, Gaginskaya E, Carré W, Waddington D, Talbot R, Völker M, Masabanda JS, Burt DW: Whole genome comparative studies between chicken and turkey and their implications for avian genome evolution. BMC Genomics. 2008, 9: 168-10.1186/1471-2164-9-168.PubMed CentralView ArticlePubMedGoogle Scholar
- Skinner BM, Robertson LB, Tempest HG, Langley EJ, Ioannou D, Fowler KE, Crooijmans RP, Hall AD, Griffin DK, Völker M: Comparative genomics in chicken and Pekin duck using FISH mapping and microarray analysis. BMC Genomics. 2009, 10: 357-10.1186/1471-2164-10-357.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Nahashon S, Feaster TK, Bohannon-Stewart A, Adefope N: An initial map of chromosomal map of chromosomal segmental copy number variations in the chicken. BMC Genomics. 2010, 11: 351-10.1186/1471-2164-11-351.PubMed CentralView ArticlePubMedGoogle Scholar
- Fontanesi L, Beretti F, Martelli PL, Colombo M, Dall'olio S, Occidente M, Portolano B, Casadio R, Matassino D, Russo V: A first comparative map of copy number variations in the sheep genome. Genomics. 2011, 97: 158-165. 10.1016/j.ygeno.2010.11.005.View ArticlePubMedGoogle Scholar
- Fontanesi L, Martelli PL, Beretti F, Riggio V, Dall'Olio S, Colombo M, Casadio R, Russo V, Portolano B: An initial comparative map of copy number variations in the goat (Capra hircus) genome. BMC Genomics. 2010, 11: 639-10.1186/1471-2164-11-639.PubMed CentralView ArticlePubMedGoogle Scholar
- Itsara A, Cooper GM, Baker C, Girirajan S, Li J, Absher D, Krauss RM, Myers RM, Ridker PM, Chasman DI, Mefford H, Ying P, Nickerson DA, Eichler EE: Population analysis of large copy number variants and hotspots of human genetic disease. Am J Hum Genet. 2009, 84: 148-161. 10.1016/j.ajhg.2008.12.014.PubMed CentralView ArticlePubMedGoogle Scholar
- McCarroll SA, Kuruvilla FG, Korn JM, Cawley S, Nemesh J, Wysoker A, Shapero MH, de Bakker PI, Maller JB, Kirby A, Elliott AL, Parkin M, Hubbell E, Webster T, Mei R, Veitch J, Collins PJ, Handsaker R, Lincoln S, Nizzari M, Blume J, Jones KW, Rava R, Daly MJ, Gabriel SB, Altshuler D: Integrated detection and population-genetic analysis of SNPs and copy number variation. Nat Genet. 2008, 40: 1166-1174. 10.1038/ng.238.View 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
- Orozco LD, Cokus SJ, Ghazalpour A, Ingram-Drake L, Wang S, van Nas A, Che N, Araujo JA, Pellegrini M, Lusis AJ: Copy number variation influences gene expression and metabolic traits in mice. Hum Mol Genet. 2009, 18: 4118-4129. 10.1093/hmg/ddp360.PubMed CentralView ArticlePubMedGoogle Scholar
- Butler MW, Hackett NR, Salit J, Strulovici-Barel Y, Omberg L, Mezey J, Crystal RG: Glutathione S-transferase copy number variation alters lung gene expression. Eur Respir J. 2011, 38: 15-28. 10.1183/09031936.00029210.PubMed CentralView ArticlePubMedGoogle Scholar
- Brønstad I, Wolff AS, Løvås K, Knappskog PM, Husebye ES: Genome-wide copy number variation (CNV) in patients with autoimmune Addison's disease. BMC Med Genet. 2011, 12: 111-10.1186/1471-2350-12-111.PubMed CentralView ArticlePubMedGoogle 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 A: Structural And Copy Number Variation In The Pig Genome. Book Structural And Copy Number Variation In The Pig Genome. 2010, Plant & Animal Genomes XVIII Conference. Town, City, January 9-13, 2010Google Scholar
- Ramayo-Caldas Y, Castelló A, Pena RN, Alves E, Mercadé A, Souza CA, Fernández 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
- Wang J, Jiang J, Fu W, Jiang L, Ding X, Liu J, Zhang Q: A genome-wide detection of copy number variations using SNP genotyping arrays in swine. BMC Genomics. 2012, 13: 273-10.1186/1471-2164-13-273.PubMed CentralView 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
- Armengol L, Pujana MA, Cheung J, Scherer SW, Estivill X: Enrichment of segmental duplications in regions of breaks of synteny between the human and mouse genomes suggest their involvement in evolutionary rearrangements. Hum Mol Genet. 2003, 12: 2201-2208. 10.1093/hmg/ddg223.View ArticlePubMedGoogle Scholar
- Kemkemer C, Kohn M, Cooper DN, Froenicke L, Högel J, Hameister H, Kehrer-Sawatzki H: Gene synteny comparisons between different vertebrates provide new insights into breakage and fusion events during mammalian karyotype evolution. BMC Evol Biol. 2009, 9: 84-10.1186/1471-2148-9-84.PubMed CentralView ArticlePubMedGoogle Scholar
- Locke DP, Sharp AJ, McCarroll SA, McGrath SD, Newman TL, Cheng Z, Schwartz S, Albertson DG, Pinkel D, Altshuler DM, Eichler EE: Linkage disequilibrium and heritability of CNPs within duplicated regions of the human genome. Am J Hum Genet. 2006, 79: 275-290. 10.1086/505653.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharp AJ, Locke DP, McGrath SD, Cheng Z, Bailey JA, Vallente RU, Pertz LM, Clark RA, Schwartz S, Segraves R, Oseroff VV, Albertson DG, Pinkel D, Eichler EE: Segmental duplications and copy-number variation in the human genome. Am J Hum Genet. 2005, 77: 78-88. 10.1086/431652.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
- 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
- Kondrashov FA, Kondrashov AS: Role of selection in fixation of gene duplications. J Theor Biol. 2006, 239: 141-151. 10.1016/j.jtbi.2005.08.033.View ArticlePubMedGoogle Scholar
- Nguyen DQ, Webber C, Ponting CP: Bias of selection on human copy-number variants. PLoS Genet. 2006, 2: e20-10.1371/journal.pgen.0020020.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohno S: Evolution by gene duplication. 1970, Springer-Verlag, New York Heidelberg Berlin, 1View ArticleGoogle Scholar
- Nguyen DQ, Webber C, Hehir-Kwa J, Pfundt R, Veltman J, Ponting CP: Reduced purifying selection prevails over positive selection in human copy number variant evolution. Genome Res. 2008, 18: 1711-1723. 10.1101/gr.077289.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Belle G, Fisher LD, Heagerty PJ, Lumley T: Biostatistics: A Methodology For the Health Sciences. 2004, Wiley, New Jersey, 291-356. Association and prediction: linear models with one predictor variable, 2, 9.View ArticleGoogle Scholar
- Giuffra E, Evans G, Törnsten A, Wales R, Day A, Looft H, Plastow G, Andersson L: The Belt mutation in pigs is an allele at the Dominant white (I/KIT) locus. Mamm Genome. 1999, 10: 1132-1136. 10.1007/s003359901178.View ArticlePubMedGoogle Scholar
- Giuffra E, Törnsten A, Marklund S, Bongcam-Rudloff E, Chardon P, Kijas JM, Anderson SI, Archibald AL, Andersson L: A large duplication associated with dominant white color in pigs originated by homologous recombination between LINE elements flanking KIT. Mamm Genome. 2002, 13: 569-577. 10.1007/s00335-002-2184-5.View ArticlePubMedGoogle Scholar
- Wehrle-Haller B: The role of Kit-ligand in melanocyte development and epidermal homeostasis. Pigment Cell Res. 2003, 16: 287-296. 10.1034/j.1600-0749.2003.00055.x.View ArticlePubMedGoogle Scholar
- Hutt KJ, McLaughlin EA, Holland MK: Kit ligand and c-Kit have diverse roles during mammalian oogenesis and folliculogenesis. Mol Hum Reprod. 2006, 12: 61-69. 10.1093/molehr/gal010.View ArticlePubMedGoogle Scholar
- Smyth GK, Speed T: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.View ArticlePubMedGoogle Scholar
- Zhou J, Lemos B, Dopman EB, Hartl DL: Copy-number variation: the balance between gene dosage and expression in Drosophila melanogaster. Genome Biol Evol. 2011, 3: 1014-1024. 10.1093/gbe/evr023.PubMed CentralView ArticlePubMedGoogle Scholar
- de Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20: 1453-1454. 10.1093/bioinformatics/bth078.View ArticlePubMedGoogle Scholar
- Guberman JM, Ai J, Arnaiz O, Baran J, Blake A, Baldock R, Chelala C, Croft D, Cros A, Cutts RJ, Di Génova A, Forbes S, Fujisawa T, Gadaleta E, Goodstein DM, Gundem G, Haggarty B, Haider S, Hall M, Harris T, Haw R, Hu S, Hubbard S, Hsu J, Iyer V, Jones P, Katayama T, Kinsella R, Kong L, Lawson D, Liang Y, Lopez-Bigas N, Luo J, Lush M, Mason J, Moreews F, Ndegwa N, Oakley D, Perez-Llamas C, Primig M, Rivkin E, Rosanoff S, Shepherd R, Simon R, Skarnes B, Smedley D, Sperling L, Spooner W, Stevenson P, Stone K, Teague J, Wang J, Wang J, Whitty B, Wong DT, Wong-Erasmus M, Yao L, Youens-Clark K, Yung C, Zhang J, Kasprzyk A: BioMart Central Portal: an open database network for the biological community. Database (Oxford). 2011, 18: bar041-Google Scholar
- Hu ZL, Fritz ER, Reecy JM: AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond. Nucleic Acids Res. 2007, 35: D604-D609. 10.1093/nar/gkl946.PubMed CentralView ArticlePubMedGoogle Scholar
- Graubert TA, Cahan P, Edwin D, Selzer RR, Richmond TA, Eis PS, Shannon WD, Li X, McLeod HL, Cheverud JM, Ley TJ: A high-resolution map of segmental DNA copy number variation in the mouse genome. PLoS Genet. 2007, 3: e3-10.1371/journal.pgen.0030003.PubMed CentralView ArticlePubMedGoogle Scholar
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.