- Research article
- Open Access
An initial comparative map of copy number variations in the goat (Capra hircus) genome
BMC Genomicsvolume 11, Article number: 639 (2010)
The goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals. Despite its importance, studies on the goat genome are still in their infancy compared to those in other farm animal species. Comparative mapping between cattle and goat showed only a few rearrangements in agreement with the similarity of chromosome banding. We carried out a cross species cattle-goat array comparative genome hybridization (aCGH) experiment in order to identify copy number variations (CNVs) in the goat genome analysing animals of different breeds (Saanen, Camosciata delle Alpi, Girgentana, and Murciano-Granadina) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome.
We identified a total of 161 CNVs (an average of 17.9 CNVs per goat), with the largest number in the Saanen breed and the lowest in the Camosciata delle Alpi goat. By aggregating overlapping CNVs identified in different animals we determined CNV regions (CNVRs): on the whole, we identified 127 CNVRs covering about 11.47 Mb of the virtual goat genome referred to the bovine genome (0.435% of the latter genome). These 127 CNVRs included 86 loss and 41 gain and ranged from about 24 kb to about 1.07 Mb with a mean and median equal to 90,292 bp and 49,530 bp, respectively. To evaluate whether the identified goat CNVRs overlap with those reported in the cattle genome, we compared our results with those obtained in four independent cattle experiments. Overlapping between goat and cattle CNVRs was highly significant (P < 0.0001) suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Genes with environmental functions were over-represented in goat CNVRs as reported in other mammals.
We describe a first map of goat CNVRs. This provides information on a comparative basis with the cattle genome by identifying putative recurrent interspecies CNVs between these two ruminant species. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.
The goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals and about 560 breeds, which constitute approximately 12% of the total number of recorded domesticated mammalian livestock breeds of the world . The diffusion of this species is mainly due to its capacity to supply milk, meat, and fibers for human consumption and use in a large number of different environments, including those poor of vegetation. In general, goat breeding represents an essential support for marginal economies in most developed and developing countries.
Despite the importance of this species, studies on the goat genome are still in their infancy compared to those in other farm animal species. A first and a second generation genetic maps of the goat genome have been obtained by Vaiman et al.  and Schibler et al.  by mapping a few hundred microsatellite markers in half-sib paternal goat families, with about 90% genome coverage. The relatively short evolutionary time separating the goat from the cattle and sheep [4–6] made it possible to use microsatellites developed in these two species to successfully genotype goats, even though interspecific priming often resulted in a marked loss of heterozygosity [2, 3]. A comparative cytogenetic map of the goat genome has been developed using cattle and sheep BAC clones . This map has been improved adding many other physically mapped genes, as recently reviewed in a compiled list including 268 genes and 144 microsatellites, roughly including 65% of the goat chromosome bands . Comparative mapping between cattle and goat (both species have 2n = 60) has shown only few rearrangements in agreement with the similarity of chromosome banding [7, 8].
Analysis of the goat genome provided a few important findings including the positional cloning of the polled intersex syndrome (PIS) locus located on goat chromosome (CHI) 1q43 and determined by a deletion of 11.7 kb containing mainly repetitive sequences [9, 10]. Additional studies have been focused on milk protein gene polymorphisms and their effects on milk production traits (e.g. [11–13]). Polymorphisms in the goat PRNP gene have been associated with susceptibility to scrapie in different breeds . Few other studies reported QTL for milk and fleece production traits and disease resistance [15–20]. Investigations of genes affecting coat colour identified polymorphisms associated with this phenotypic trait [21, 22]. In particular, the Agouti locus in goat was shown to be highly variable including missense mutations and copy number variation (CNV) .
Recent studies have shown that copy number variants, defined as intraspecific gains or losses of ≥ 1 kb of genomic DNA [23, 24], represent an important source of variability of mammalian genomes (~0.4-25% of the genome) as reported in human (e.g. [25–34]), chimpanzee [35, 36], rhesus macaque , mouse [38–42], rat [43, 44], dog [45, 46], pig , and cattle [48–52]. CNVs can change gene structure and dosage, can regulate gene expression and function and for these reasons they have potentially more effects than the most frequent single nucleotide polymorphisms (SNPs) in determining phenotypic differences [43, 53–56]. CNVs can represent benign polymorphic variants even if in many other cases they are associated with human Mendelian and complex genetic disorders (reviewed in [57, 58]). In farm animals, several traits are caused by CNV affecting genes or gene regions. For example, the Dominant white locus in pigs includes alleles determined by duplications of the KIT gene [59, 60]. CNV also affects the Agouti locus in sheep and goats and contributes to the variability of coat colour in these two species [22, 61]. CNV in intron 1 of the SOX5 gene causes the pea-comb phenotype in chicken  and the late feathering locus in this avian species includes a partial duplication of the PRLR and SPEF2 genes .
Genome-wide discovery and frequency evaluation of CNVs have been possible with the development of high-resolution array comparative genome hybridisation (aCGH) and, subsequently, with data analysis of high-density SNP platforms and paired end and deep sequencing approaches [64–69]. An advantage of aCGH is that hybridisation can be performed using heterologous DNA, i.e. genomic DNA of a different species but close to that used to develop the array, taking advantages from completely sequenced, assembled, and richly annotated genomes. Cross species aCGH experiments have been successfully applied using human arrays to analyse CNVs in chimpanzee and other primates [35, 36, 70], and using chicken based arrays to identify CNVs in turkey , duck , and zebra-finch  genomes.
Here we designed a cross species cattle-goat aCGH experiment in order to identify CNVs in goats of different breeds (both cosmopolitan and local) using information of the cattle genome and we obtained a first comparative map of CNVs of the Capra hircus genome.
Results and discussion
Identification of goat CNVs and comparative analysis between goat and cattle CNVRs
The goat genome has not been sequenced yet and the cattle is the closest species to the goat for which an assembled genome is available [74, 75]. In order to give a preliminary evaluation of the extent of sequence identity between these two ruminant species, we compared goat genomic sequences longer than 1 kb available in EMBL database with homologous cattle genomic regions. The average sequence identity between these two species was 82.79% (over 166.8 kb of aligned sequences), that increased at 93.77% if only exonic sequences (32.5 kb) were considered (data not shown). Even if this rough evaluation cannot give a complete picture of sequence divergence between cattle and goat, it indicates that cattle-goat cross-species DNA hybridization is possible as also demonstrated in physical mapping of cattle BAC and YAC clones to the goat genome and vice versa as well as by interspecific use of microsatellites [2, 3, 7, 76]. Therefore, to identify CNVs in goats, we carried out a cross-species aCGH experiment using a cattle (Bos taurus) custom tiling array including ~385,000 oligonucleotide probes and goat genomic DNA obtained from 9 goats of four different breeds (3 Saanen, 1 Camosciata delle Alpi, 3 Girgentana, and 2 Murciano-Granadina goats), chosen according to their differences in production and morphological traits, and origin (Figure 1). Saanen is the most popular dairy breed that was originated in Switzerland, highly selected for heavy milk production and with completely white/cream coat colour. Camosciata delle Alpi is an Alp mountain breed of the Chamois group. Girgentana is a Sicilian breed in an endangered status that is probably of Afghan and Himalayan origin, well adapted to the dry Sicilian environment and with very mild behavior. Girgentana goats have long corkscrew horns and cream/light-grey coat colour with, usually, a few small red spots around eyes and ears. Murciano-Granadina is a native Spanish breed with two colour types, solid black or solid brown (caoba).
Specific criteria were used to call CNVs in this cross-species aCGH experiment. CNVs were reported using 10 different algorithms developed for data segmentation and averaging their results as implemented in the CGHweb server . The averaged log2 ratio threshold used to call CNVs was empirically determined considering the number of false positives called in the reference DNA self-self hybridization and the validation obtained by semiquantitative fluorescent multiplex PCR (SQF-PCR) (see Methods and data reported below for details). Applying these criteria, on the whole we identified 161 CNVs (Table 1 and Additional file 1). The largest number of CNVs was identified in the Saanen breed and the lowest number was reported for the Camosciata delle Alpi analysed goat. About 74% of the identified CNVs were observed in only one breed. On average each sample contained 17.9 CNVs. This result is comparable to what has been obtained in similar aCGH experiments that analysed the cattle , dog , and chicken  genomes, in which 11.6, 17.2 and 9.6 CNVs were evidenced for each animal in the three species, respectively (Table 2). Two CNVs (not overlapped to any other CNV) were also called in the self-self hybridization providing a rough estimation of the false discovery rate (FDR) among the set of identified CNVs (FDR = 11%). However, technical issues, like sequence divergence between the reference genome and the hybridized DNA and heterogeneity of DNA quality among different samples, make it difficult to precisely estimate the experimental FDR. The estimated FDR using the self-self hybridization data in this cross-species experiment is a little bit larger than that obtained with homologous DNA hybridization in aCGH. Indeed, in an aCGH study in dog, FDR was about 3% due to only one CNV detected in the self-self hybridization . However, in another experiment carried out in humans, FDR was estimated to be 8-24% .
CNV regions (CNVRs) were determined by aggregating overlapping CNVs identified in different animals as previously reported [28, 51] and considering a conservative approach due to the specificity of our experiment (see Methods for details). On the whole, we detected 126 CNVRs covering about 11.39 Mb of the virtual goat genome referred to the bovine genome, version Btau_4.0 (Figure 2 and Additional file 2). This fraction corresponds to 0.432% of the latter genome, considered adding bases in the 29 autosomes and the X chromosome assembled in the Btau_4.0 version (11.39 Mb/2634 Mb). The chrUnAll (unassembled scaffolds) of the Btau_4.0 version was not included in the tiling array due to difficulties in interpreting the results that might have been derived by the short assembled fragments and mapping uncertainty . In addition, the tiling arrays included 4,673 oligonucleotides designed on a portion of BTA13 (from nucleotide position 48 Mb to nucleotide position 78 Mb) derived from the UMD Bos taurus v. 2.0 assembly . This additional BTA13 portion was included in the tiling arrays because the agouti signaling protein (ASIP) gene was not correctly assembled in the Btau_4.0 version of BTA13 and was reported in unassembled scaffolds only. This UMD 2.0 BTA13 region was added as an internal control because we previously demonstrated that the ASIP gene and the close S-adenosylhomocysteine hydrolase (AHCY) gene are included in a goat CNVR . This CNVR may represent a recurrent interspecies CNVR since the same two genes are involved in a large duplicated region in sheep . Considering this additional CNVR, on the whole we identified 127 CNVRs for a total of 11.47 Mb (Figure 2 and Additional file 2). These 127 CNVRs included 86 losses and 41 gains, whereas none reported both events. Of these CNVRs, 14 were found in multiple animals of different breeds (no. = 7) or in multiple animals of the same and different breeds (no. = 6) or in multiple animals of the same breed only (no. = 1). All other events (no. = 113) were found only in one animal. The regions that were affected by these gain or loss events in goats ranged from 24,605 bp (BTA22, CNVR no. 101) to 1,075,778 bp (BTA17, CNVR no. 90) with mean and median equal to 90,292 bp and 49,530 bp, respectively (Additional file 2). Using a similar aCGH experiment Liu et al.  reported that cattle CNVRs are on average a little bit larger than those we reported in goats (mean and median equal to 159,031 bp and 89,053 bp, respectively). This is in contrast to Fadista et al.  who, using more dense arrays in another aCGH experiment reported a median size of CNVRs equal to 16.7 kb. Mean and median differences between goat and cattle experiments might be due to i) the cross species experiment performed in goat that might not be able to correctly identify the borders of the CNVRs, ii) the oligonucleotide density in the aCGH experiments, iii) true differences between species/breeds. It is interesting to note that the number of loss events was about two fold the number of gain events in both cattle and goat. This finding might be derived by both biological and technical reasons. Non-allelic homologous recombination, which seems to be one of the most important mechanisms responsible for CNV formation, has been shown to generate more deletions than duplications . On the other hand, aCGH detection methods seem to favor the identification of deletions as reported by several other studies [28, 47, 51, 52].
Additional file 3 reports the extension of CNVRs distributed for the different bovine chromosomes used in the comparative analysis with the goat genome. In only three chromosomes (BTA11, BTA20, and BTA21) we did not identify any CNVRs. BTA5 included the largest number of CNVRs (no. = 11), whereas BTA17, BTA10, and BTA18 included the largest extension of regions affected by CNVs (1.6%, 1.3% and 1.0% of their length, respectively) (Figure 2). In cattle a similar aCGH experiment (that however included a larger number of animals ) showed the greatest enrichment for CNVRs on BTA5, BTA15, BTA18, BTA27, BTA29, and BTAX. Comparative mapping and chromosome banding similarities between cattle and goat indicate highly conserved synteny between these two ruminant species even if a few rearrangements have been evidenced mainly on CHI14 containing a small BTA9q11-q13 segment, and some other gene order rearrangements for chromosomes CHI2 compared to BTA2, CHI19 compared to BTA19, and CHIX compared to BTAX [7, 8]. CNVRs were evidenced in these interested chromosome regions. However, it is not possible to evidence if CNVs are precisely positioned in rearranged regions because of the low resolution of the rearrangements so far described between goat and cattle chromosomes due to the few mapped genes.
To evaluate if CNVRs we identified in goats overlap with CNVRs reported in cattle, we compared our results with those obtained in four independent cattle experiments [49–52] carried out i) using aCGH including ~385,000 tiling oligonucleotides (177 CNVRs ) or ii) including 6.3 million of probes (304 CNVRs ) and iii) using the Illumina BovineSNP50 BeadChip containing about 50K SNPs as reported by Matukumalli et al. (; 79 CNVRs) and by Bae et al. (; 368 CNVRs) (Additional file 4 and Additional file 5). Overlapping between aCGH results obtained in goat and cattle was highly significant (P < 0.0001) for both cattle experiments (17 and 11 goat CNVRs overlapped with cattle CNVRs identified by Liu et al.  and Fadista et al. , respectively). Only two goat CNVRs matched cattle CNVRs identified with the SNP panel [49, 50], therefore overlapping was not significant. A similar bias on the common CNVRs between aCGH and SNP genotyping experiments was also evident comparing the cattle data obtained by the two methods (Additional file 5 and ). This could be due to resolution and genome coverage differences between the two platforms . However, merging all cattle CNVRs reported in the four different experiments (on the whole 764 unique cattle CNVRs were obtained [49–52], Additional file 4), overlap with goat CNVRs was highly significant (P < 0.0001), confirming the results obtained considering the different cattle datasets separately (Additional file 5).
As the goat genome is not sequenced yet, we could not evaluate if the goat CNVRs have similar sequence characteristics in goat and cattle. Segmental duplications have been shown to significantly overlap with CNVRs in cattle [51, 52] as well as in several other species [31, 32, 34, 35, 37, 39, 46]. As segmental duplications might facilitate non-allelic homologous recombination, it is likely that they are as well involved in the mechanisms that produce CNVs in goats. The overlapping CNVRs between goat and cattle might represent homoplastic recurrent interspecies CNVRs probably driven by genomic regions prone to instability present in the cattle-goat common ancestor that might have been retained in the genomes of both extant species. Indeed, cattle and goat share a common ancestor in the early Miocene about 17-23 Million years ago [4–6]. Similar reasoning could be considered for the CNVR that includes the ASIP gene  for which a recurrent CNVR has been reported in sheep , but not in cattle. Sheep and goat lineages separated about 6-14 millions of years before present . Significant overlap of CNVRs among different species has been also observed comparing the human with both chimpanzee and rhesus macaque genomes [35–37]. These two non-human primate species diverged from the human lineage about 6 and 25 million years ago [81, 82]. Together these results suggest that certain genomic regions are prone to recurrent CNV formation and instability in both the primate and the Artiodactyla evolutionary lineages. However, a change in the formation process of CNVs and segmental duplications in the human genome might be occurred in recent evolutionary history . It will be interesting to evaluate if this has occurred in other lineages. A comparative analysis of CNVRs identified in cattle, goat, and sheep can open perspectives to evaluate the evolutionary mechanisms determining CNV formation during the mammalian evolution.
Validation and gene content of CNVRs
Several results suggest most of our CNVs are correctly identified. First of all, the number of CNVs was lower in the analysed Camosciata delle Alpi goat (no. = 8; Table 1). This was expected because the aCGH reference was a sample of genomic DNA of another goat of the same breed. Similar results were also reported in mouse and dog aCGH studies that used a reference genomic DNA of an animal of the same breed/line of others that have been analysed for CNV discovery [39, 45]. In addition, an internal control included in the design of the tiling arrays was represented by a portion of BTA13 derived from the UMD 2.0 assembly that contains the ASIP and AHCY genes that we previously showed to be affected by CNV in goats . All goats that were shown to have multiple copies of the ASIP and AHCY genes with different methods  reported evidence of gain in the aCGH experiment, compared to the reference Camosciata delle Alpi genomic DNA. The results for this region were also used to set an empirical threshold to call CNVs by using the CGHweb platform and multiple algorithms  (see Methods).
Other three CNVRs identified on the goat chromosomes corresponding to BTA4, BTA10, and BTA17 (Table 3 and Additional file 6) were directly validated with SQF-PCR as reported for the ASIP and AHCY genes . The CNVR of BTA4 (CNVR no. 25) included a gene coding for the GTPase, IMAP family member 1 (GIMAP1; ENSBTAG00000001198) that was used to design the primers for PCR validation. Semiquantitative fluorescent multiplex PCR confirmed the gain of copies already reported in the goats used in the aCGH experiment and in additional goats of the Saanen, Girgentana and Murciano-Granadina breeds (Figure 3A) (Additional file 7). The evaluated CNVR of BTA10 (CNVR no. 61) was originally identified in a Murciano-Granadina goat by aCGH (Additional file 1 and Additional file 2). The bovine sequence coding for an unknown transcript (ENSBTAG00000027170) positioned in this region was used to design PCR primers for validation. However, analyzing the bovine genome by BLAT using the ENSBTAG00000027170 sequence, a large number of significant hits (E-value < 4E-233, identity > 89%) were obtained, and the hits included only sequences assembled in the BTA10 encompassing the goat CNVRs no. 61 and 62 or included in unassigned scaffolds (data not shown). These regions of BTA10 have been already reported to contain CNVs in different cattle breeds [51, 52] (Additional file 4). Analyzing by SQF-PCR goats other than those used for aCGH, we identified that this complex region (CNVRs no. 61 and 62) is affected by both gain and loss of DNA copies that did not occur specifically in the evaluated goat breeds (Additional file 7). Therefore, even if these CNVRs were classified as gain (Figure 3B) (Additional file 2), they should be considered together as gain/loss CNVRs. CVNR no. 90 on BTA17 (the largest we identified) includes the nuclear receptor subfamily 3, group C, member 2 (NR3C2) gene (ENSBTAG00000027182). PCR primers were tested on additional animals as reported for the previous two CNVRs. Only and all Girgentana goats showed gain of DNA copies in this region (Figure 3C), including the three Girgentana goats used for aCGH analyses (Additional file 7), even if the signal for one of them did not trespass the averaged log2 ratio threshold of 0.175. The NR3C2 gene encodes the mineralocorticoid receptor, which mediates aldosterone actions on salt and water balance within distal nephron cells, with crucial effects on blood pressure and potassium homeostasis. The protein acts as a ligand-dependent transcription factor that binds to mineralocorticoid response elements in order to transactivate target genes. Mutations in this gene in humans cause autosomal dominant pseudohypoaldosteronism type I, a disorder characterized by renal resistance to aldosterone as well as salt wasting, dehydration, hyperkalemia, metabolic acidosis and failure to thrive in the newborn . Defects in this gene are also associated with early onset hypertension , whereas overexpression of this gene in forebrain decreases anxiety-like behavior . It is tempting to speculate that additional copies of the NR3C2 gene (if functional) could contribute to specific adaptation traits to harsh and dry environments and to the very mild behavior of the Girgentana goats compared to other breeds.
Considering all 127 goat CNVRs, 78 of them (61.4%) partially or completely spanned cattle Ensembl annotated genes (Btau_4.0 version), including 249 Ensembl peptides, corresponding to 199 unique Ensembl cattle genes, representing 261, 126, and 864 gene ontology (GO) categories for molecular function, cellular components and biological processes, and 870 different PANTHER terms (Additional file 8 and Additional file 9). For example, considering GO annotation for biological processes, several GO categories are significantly (P < 0.001) over-represented in goat CNVRs (Table 4; Additional file 10). A few of these GO terms (response to stimulus; defense response) have been also already reported to be over-represented in CNVRs of the mouse and human genomes [39, 87]. In general, goat CNVRs resulted particularly enriched in "environmental" functions (Table 4; Additional file 10 and Additional file 11) as already reported in other species [39, 46, 50–52, 87]. This indication might be important in understanding the evolutionary and selection processes that determined the occurrence and diffusion of this variability in the goat genome. It could be possible that positive selection on duplications (based on functional effects) has occurred for these particular gene categories or that these enrichments might instead have been arisen from nonuniform negative selection on gene copy changes. This could be due to the fact that duplication or deletion of nonenvironmental genes might be more frequently deleterious than copy number changes of environmental genes . The complex domestication process and artificial selection for breeding purposes (including bottleneck and expansion) that largely contributed to establish goat breed differences  might further complicate the evaluation of these two different hypotheses.
Considering the 199 unique Ensembl cattle genes included in the CNVRs identified in goats, for 119 we retrieved a human orthologous gene. Mutations in only 8 of these genes cause Mendelian disorders or are associated with genetic diseases in human (Additional file 12). None of these 199 cattle genes is involved in any reported genetic disease in goat, sheep or cattle. Among the goat genes already mapped  only one (TTN) is included in a CNVR (Additional file 12).
Few QTL studies have been reported in goats so far [15–20] and all used breeds/populations not included in our CNV study, except one in which a few Saanen crosses were analysed . In addition, the identified QTL have very large confidence intervals including most of the CNVRs we observed in several chromosomes (data not shown). Therefore a comparison among CNVRs we identified in goats, their gene content and QTL regions is not very informative. However, it is interesting to note that Bolormaa et al.  reported putative QTL for faecal worm egg and eosinophil counts on CHI23, in the region including the major histocompatibility complex (MHC), in which we reported two CNVRs, one affecting a MHC class I antigen gene (CNVR no. 108, gain of copies) and the other one including MHC class II alpha and beta chain genes (CNVR no. 107, loss of copies) (Additional file 8 and Additional file 9). In sheep, several reports have identified polymorphisms in the class I and class II regions of the MHC being associated with resistance to nematodes . It will be interesting to evaluate if these CNVs we identified in goats are associated with resistance to nematode infection and other diseases.
We provide a first comparative map of CNVRs in the goat genome using a cross-species aCGH experiment based on the cattle genome. Considering the limited number of analysed animals and breeds and the fact that the cross-species hybridization could have limited the detection power of this study, the reported goat CNVRs largely underestimate the true number of this kind of variation in the goat genome. Additional studies including other breeds should be carried out in order to better evaluate the extension and distribution of CNVs in the genome of this farm animal species. However, it appeared that possible evolutionary conserved mechanisms might be the causative factors of putative recurrent interspecies CNVs between cattle and goat. Using this cross-species design it seems possible to analyse several other genomes of the Bovidae family in order to obtain comparative CNV maps in other species close to the cattle and provide additional evidence on the co-occurrence of CNVs in particular chromosome regions. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.
We analysed CNVs in the goat genome by means of a cross-species aCGH experiment using the Roche NimbleGen platform (Roche NimbleGen Inc., Madison, WI; http://www.nimblegen.com) based on custom tiling arrays designed on the bovine (Bos taurus) genome, Btau_4.0 version , including a fraction of BTA13 of the University of Maryland (UMD) Bos taurus v. 2.0 assembly (ftp://ftp.cbcb.umd.edu/pub/data/Bos_taurus/Bos_taurus_UMD_2.0/). Arrays contained ~385,000 probes on a single slide to provide an evenly distributed coverage with an average interval of ~6 kb for the Btau_4.0 genome. The BTA13 of the UMD v. 2.0 assembly was included from nucleotide position 48 M bp to nucleotide position 78 M bp (4,673 oligonucleotides and average spacing of ~6 kb). This chromosome region was analysed as internal control because it contains the ASIP gene, not assembled in the BTA13 of the Btau_4.0 version. We previously showed that this goat gene is included in a CNVRs in different goat breeds .
Goat genomic DNA was extracted from blood of 2 Camosciata delle Alpi, 3 Girgentana, 3 Saanen, 1 black and 1 brown Murciano-Granadina goats using the Wizard® Genomic DNA Purification kit (Promega Corporation, Madison, WI). All analysed animals were females. Reference DNA sample of one Camosciata delle Alpi goat (C1) was labeled with Cy5 and co-hybridised with the other test DNA samples labeled with Cy3 on 9 different arrays. A self hybridisation (reference labeled by both Cy5 and Cy3) was carried out in another array. Hybridization and array scanning were performed by Roche NimbleGen as previously described . Data normalization was conducted using the normalize.qsline method from the Bioconductor package in R . Then data were analysed for each hybridization using normalized log2 ratios using the CGHweb server (http://compbio.med.harvard.edu/CGHweb/) that includes multiple algorithms. We used the self-self hybridisation and the BTA13 control region to define a suitable threshold to apply to the CGHweb calls in order to minimize false positives. Specifically we retained predicted CNVs if it had at least five consecutives datapoints supporting it (considering an average of probe values inside a smoothing window of five), thus limiting the minimum CNV size to about 30 kb, even if this resolution can vary in different regions depending on the relative distance of the probes that can be different from the averaged ~6 kb. Pointwise averaging of all computed profiles and maps of gains/losses for smoothed/segmented obtained from several algorithms (Lowess, Wavelet, Quantreg, ruavg, CBS, CGHseg, BioHMM, cghFLasso, GLAD, and FASeg) and summary data were generated. Pointwise averaging was shown to have good performances in calling alteration of copy number  and was chosen to compensate possible signal differences among probes in the cattle-goat heterologous experiment. Summary data were considered to call gain/loss in a chromosome region and to compile a high confidence set of CNVs. Then CNVs were called considering a conservative approach joining regions of at least 4-5 contiguous probes with CNV signal separated by up to three probes without CNV signal in the same individual (Additional file 1). This adjustment was applied in order to overcome possible signal losses or hybridisation problems in the cross-species aCGH experiment. CNVRs were reported aggregating overlapping or partially overlapping CNVs in different animals as previously reported [28, 51] and applying the same criteria for CNVs within individuals (Additional file 2). The false discovery rate (FDR) was estimated based on the observation of 2 false positives in the self-self hybridisation. A rough estimate of the FDR is the expected number of false positives per array (n. 2) times the number of total arrays divided by the total number of CNVs (n. 161), resulting in an estimated FDR of 11%. This calculation should be considered only an approximation because it does not consider the potential for varying false positive rates across arrays. Based on these criteria the averaged log2 ratio threshold to call gains and losses  was empirically established at 0.175 considering the results obtained for the ASIP gene region. Four goats out of five with independent validated CNV  in this gene reported an averaged log2 ratio > 0.175, therefore this value was used as threshold even if in another goat the averaged log2 ratio for the ASIP region was 0.156. However, even if this latter value did not change the self-self FDR results, we used as threshold the value of 0.175 because the self-self hybridization could not fully reflect the variance of our 9 test experiments and we preferred a low false-positive rate even at the expense of having more false negatives in our dataset.
Validation of CNVs
Validation of CNVs was performed by semiquantitative fluorescent multiplex PCR (SQF-PCR) as previously reported [22, 92] using genomic DNA of the same goats analysed in the aCGH experiment and genomic DNA of additional goats (additional 8 Saanaen, 12 Girgentana, and 14 Murciano-Granadina) extracted as reported above. Briefly, two internal control regions known to have no CNV (DGAT1 and MC1R gene fragments) and CNVRs of interest were co-amplified in multiplex PCR under quantitative PCR conditions (with forward primers labelled in 5' with 6FAM) and the products were separated by capillary electrophoresis using an ABI3100 Avant sequencer (Applied Biosystems, Foster City, CA, USA) . Peak heights of regions of interest were normalized against those of the internal controls. Primer pairs for control gene fragments and analysed CNVRs are reported in Additional file 6. SQF-PCR was performed in a total volume of 10 μL using 1-6 pmol of each primer and the conditions reported in Additional file 6. PCR profile was as follows: 5 min at 95°C; 20-22 amplification cycles of 30 sec at 95°C, 30 sec at 59°C, 30 sec at 72°C; 5 min at 72°C. Capillary electrophoresis was performed using 1 μL of reaction product, diluted in 10 μL of Hi-Di formamide (Applied Biosystems), and added with 0.1 μL of Rox labelled DNA ladder (500HD Rox, Applied Biosystems). Peak heights were obtained using GeneScan software v. 3.7 (Applied Biosystems). DNA dosages were calculated by comparing the normalized peak height ratios of goats of interest with the average normalized ratios of the reference Camosciata delle Alpi goat as follows: the peak height of a fragment of interest was divided by the peak height of the internal control; the averaged value obtained from at least two analyses for each sample/target region was divided by the same averaged value obtained for the control goat DNA. We adopted the theoretical values of 1.5, 2.0, 2.5, and so on for a gain of multiple of one, two, three or other copies, respectively, compared to the copy content (unknown) of the reference DNA sample. Similarly, a loss of one set of copies (or one copy in case of a simple duplication) would theoretically result in a value of 0.5. These values should be considered only approximation of the copy number content as the objective was to validate the results obtained with aCGH and not to precisely estimate the number of copies of the analysed DNA fragments.
Bioinformatic and computational analyses
Capra hircus genomic sequences longer than 1 kb and including complete coding sequences were retrieved from EMBL database (Sept. 2010). Sequences were clustered with BLASTclust http://blast.ncbi.nlm.nih.gov/ on the basis of their identity (> 98%) resulting in 30 sequences covering on the whole 196,665 bp. These sequences were aligned with homologous cattle transcript regions identified using BLASTN on the basis of the best hits. The global sequence alignment without end-gap penalty was performed with LALIGN program http://www.ch.embnet.org/software/LALIGN_form.html. Exonic regions in goat sequences were defined according to the cattle annotation of the Btau_4.0 genome version http://www.ensembl.org/Bos_taurus/Info/Index. The goat CNVRs were mapped on the Btau_4.0 version of the bovine genome. To determine whether goat and cattle CNVRs occur in orthologous regions more often than expected by chance we considered the data reported for cattle in four different experiments [49–52]. The data reported in these four studies were considered separately due to differences in the methods and populations used for CNV detection. A merged list of the CNVRs reported in these investigations was also compiled (Additional file 4). In one of these cattle studies , CNVs were reported with reference to the Btau_3.0 version, therefore the LiftOver tool at the UCSC Genome Bioinformatics Site http://genome.ucsc.edu/cgi-bin/hgLiftOver was used to map CNVs coordinates on the Btau_4.0 version. In this case, only 45 out of the reported 79 CNVs were successfully re-mapped. Within each experiment, overlapping CNVs were fused to define CNVRs. These procedures ended up with 37, 368, 177, and 266 CNVRs for Matukumalli et al. , Bae et al. , Liu et al. , and Fadista et al.  experiments, respectively, for a total of 764 combined CNVRs (Additional file 4). The number of overlapping segments between each pair of CNVR sets was computed and the overlap significance was evaluated with a permutation test . For each experiment, 10,000 artificial random rearrangements of the CNVRs were generated and mapped on the Btau_4.0 bovine genome. The CNVR length distribution was preserved in each random rearrangement. In order to evaluate the significance of the overlap between two CNVR sets, we computed the distribution of the number of overlapping segments between one of the CNVR sets and the 10,000 random rearrangements of the other one. The reported P-value is the fraction of random CNVR rearrangements that obtain at least the same number of overlapping segments as the real one.
Goat CNVRs superimposing with cattle transcripts annotated in the Btau_4.0 version were determined on the basis of the genome coordinates, without imposing a minimum overlap threshold. Gene ontology terms associated with bovine transcripts were downloaded with the Ensembl BioMart retrieval system http://www.ensembl.org/biomart/index.html and the complete annotation was obtained by reconstructing the complete list of ancestors of each term in the directed acyclic graph described by the OBO file downloaded from the Gene Ontology web site on May 2010 http://www.geneontology.org/. The GOTermFinder tool was adopted for this task http://search.cpan.org/dist/GO-TermFinder/. We computed the occurrence of each term in the set of transcripts overlapping with goat CNVRs and we compared it with the occurrence of the same term in the whole bovine genome (Btau_4.0 version). The Fisher exact test was adopted to assess the significance of the overrepresentation of the terms in the set of transcripts overlapping with the goat CNVRs. The multiple-hypothesis correction  was adopted for discriminating the significant terms at different False Discovery Rates (FDR): 0.001, 0.01, 0.05, and 0.1.
To supplement the functional annotation, PANTHER Molecular Function terms were assigned to all bovine transcripts using the Hidden Markov Model scoring tools of the PANTHER (Protein ANalysis THrough Evolutionary Relationships) Classification System version 6.1 http://panther6.ai.sri.com/tools/hmmScoreForm.jsp. Similarly to the GO annotation, the distribution of the PANTHER terms in the set of transcripts overlapping with goat CNVRs was compared with the occurrence in the whole genome and the significance of the overrepresentation was evaluated with the Fisher exact test adopting the multiple-hypothesis correction.
aCGH data have been submitted to the gene expression omnibus http://www.ncbi.nlm.nih.gov/geo/ under the accession number GSE24436.
FAO: The State of the World's Animal Genetic Resources for Food and Agriculture, Rome. 2007
Vaiman D, Schibler L, Bourgeois F, Oustry A, Amigues Y, Cribiu EP: A genetic linkage map of the male goat genome. Genetics. 1996, 144: 279-305.
Schibler L, Vaiman D, Oustry A, Giraud-Delville C, Cribiu EP: Comparative gene mapping: a fine-scale survey of chromosome rearrangements between ruminants and humans. Genome Res. 1998, 8: 901-915.
Allard MW, Miyamoto MM, Jarecki L, Kraus F, Tennant MR: DNA systematics and evolution of the artiodactyl family Bovidae. Proc Natl Acad Sci USA. 1992, 89: 3972-3976. 10.1073/pnas.89.9.3972.
Gatesy J, Yelon D, DeSalle R, Vrba ES: Phylogeny of the Bovidae (Artiodactyla, Mammalia), based on mitochondrial ribosomal DNA sequences. Mol Biol Evol. 1992, 9: 433-446.
Matthee CA, Davis SK: Molecular insights into the evolution of the family Bovidae: a nuclear DNA perspective. Mol Biol Evol. 2001, 18: 1220-1230.
Schibler L, Di Meo GP, Cribiu EP, Iannuzzi L: Molecular cytogenetics and comparative mapping in goats (Capra hircus, 2n = 60). Cytogenet Genome Res. 2009, 126: 77-85. 10.1159/000245908.
Cribiu EP, Di Berardino D, Di Meo GP, Eggen A, Gallagher DS, Gustavsson I, Hayes H, Iannuzzi L, Popescu CP, Rubes J, Schmutz S, Stranzinger G, Vaiman A, Womack J: International System for Chromosome Nomenclature of Domestic Bovids (ISCNDB 2000). Cytogenet Cell Genet. 2001, 92: 283-299. 10.1159/000056917.
Schibler L, Cribiu EP, Oustry-Vaiman A, Furet JP, Vaiman D: Fine mapping suggests that the goat Polled Intersex Syndrome and the human Blepharophimosis Ptosis Epicanthus Syndrome map to a 100-kb homologous region. Genome Res. 2000, 10: 311-318. 10.1101/gr.10.3.311.
Pailhoux E, Vigier B, Chaffaux S, Servel N, Taourit S, Furet JP, Fellous M, Grosclaude F, Cribiu EP, Cotinot C, Vaiman D: A 11.7-kb deletion triggers intersexuality and polledness in goats. Nat Genet. 2001, 29: 453-458. 10.1038/ng769.
Barbieri ME, Manfredi E, Elsen JM, Ricordeau G, Bouillon J, Grosclaude F, Mahé MF, Bibé B: Influence du locus de la caséine αs1 sur les performances laitières et les paramètres génétiques des chèvres de race Alpine. Genet Sel Evol. 1995, 27: 437-450. 10.1186/1297-9686-27-5-437.
Hayes B, Hagesaether N, Adnøy T, Pellerud G, Berg PR, Lien S: Effects on production traits of haplotypes among casein genes in Norwegian goats and evidence for a site of preferential recombination. Genetics. 2006, 174: 455-464. 10.1534/genetics.106.058966.
Gigli I, Maizon DO, Riggio V, Sardina MT, Portolano B: Short communication: casein haplotype variability in Sicilian dairy goat breeds. J Dairy Sci. 2008, 91: 3687-3692. 10.3168/jds.2008-1067.
Vaccari G, Panagiotidis CH, Acin C, Peletto S, Barillet F, Acutis P, Bossers A, Langeveld J, van Keulen L, Sklaviadis T, Badiola JJ, Andreéoletti O, Groschup MH, Agrimi U, Foster J, Goldmann W: State-of-the-art review of goat TSE in the European Union, with special emphasis on PRNP genetics and epidemiology. Vet Res. 2009, 40: 48-10.1051/vetres/2009031.
Roldán DL, Rabasa AE, Saldaño S, Holgado F, Poli MA, Cantet RJ: QTL detection for milk production traits in goats using a longitudinal model. J Anim Breed Genet. 2008, 125: 187-193. 10.1111/j.1439-0388.2007.00714.x.
Marrube G, Cano EM, Roldan DL, Bidinost F, Abad M, Allain D, Vaiman D, Taddeo H, Poli MA: QTL affecting conformation traits in Angora goats. Small Rum Res. 2007, 71: 255-263. 10.1016/j.smallrumres.2006.07.008.
Cano E, Marrube G, Roldan D, Bidinost F, Abad M, Allain D, Vaiman D, Taddeo H, Poli M: QTL affecting fleece traits in Angora goats. Small Rum Res. 2007, 71: 158-164. 10.1016/j.smallrumres.2006.06.002.
Cano EM, Debenedetti S, Abad M, Allain D, Taddeo HR, Poli MA: Chromosomal segments underlying quantitative trait loci for mohair production in Angora goats. Anim Genet Res Information. 2009, 45: 107-112. 10.1017/S1014233909990423.
Mohammad Abadi MR, Askari N, Baghizadeh A, Esmailizadeh AK: A directed search around caprine candidate loci provided evidence for microsatellites linkage to growth and cashmere yield in Rayini goats. Small Rum Res. 2009, 81: 146-151. 10.1016/j.smallrumres.2008.12.012.
Bolormaa S, van der Werf JH, Walkden-Brown SW, Marshall K, Ruvinsky A: A quantitative trait locus for faecal worm egg and blood eosinophil counts on chromosome 23 in Australian goats. J Anim Breed Genet. 2010, 127: 207-214. 10.1111/j.1439-0388.2009.00824.x.
Fontanesi L, Beretti F, Riggio V, Dall'Olio S, González EG, Finocchiaro R, Davoli R, Russo V, Portolano B: Missense and nonsense mutations in melanocortin 1 receptor (MC1R) gene of different goat breeds: association with red and black coat colour phenotypes but with unexpected evidences. BMC Genet. 2009, 10: 47-10.1186/1471-2156-10-47.
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.
Feuk L, Carson AR, Scherer SW: Structural variation in the human genome. Nat Rev Genet. 2006, 7: 85-97. 10.1038/nrg1767.
Freeman JL, Perry GH, Feuk L, Redon R, McCarroll SA, Altshuler DM, Aburatani H, Jones KW, Tyler-Smith C, Hurles ME, Carter NP, Scherer SW, Lee C: Copy number variation: new insights in genome diversity. Genome Res. 2006, 16: 949-961. 10.1101/gr.3677206.
Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C: Detection of large-scale variation in the human genome. Nat Genet. 2004, 36: 949-951. 10.1038/ng1416.
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.
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.
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.
Qiao Y, Liu X, Harvard C, Nolin SL, Brown WT, Koochek M, Holden JJ, Lewis ME, Rajcan-Separovic E: Large-scale copy number variants (CNVs): distribution in normal subjects and FISH/real-time qPCR analysis. BMC Genomics. 2007, 8: 167-10.1186/1471-2164-8-167.
Wong KK, deLeeuw RJ, Dosanjh NS, Kimm LR, Cheng Z, Horsman DE, MacAulay C, Ng RT, Brown CJ, Eichler EE, Lam WL: A comprehensive analysis of common copy-number variations in the human genome. Am J Hum Genet. 2007, 80: 91-104. 10.1086/510560.
Perry GH, Ben-Dor A, Tsalenko A, Sampas N, Rodriguez-Revenga L, Tran CW, Scheffer A, Steinfeld I, Tsang P, Yamada NA, Park HS, Kim JI, Seo JS, Yakhini Z, Laderman S, Bruhn L, Lee C: The fine-scale and complex architecture of human copy-number variation. Am J Hum Genet. 2008, 82: 685-695. 10.1016/j.ajhg.2007.12.010.
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.
Yim SH, Kim TM, Hu HJ, Kim JH, Kim BJ, Lee JY, Han BG, Shin SH, Jung SH, Chung YJ: Copy number variations in East-Asian population and their evolutionary and functional implications. Hum Mol Genet. 2010, 19: 1001-1008. 10.1093/hmg/ddp564.
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AW, Robson S, Stirrups K, Valsesia A, Walter K, Wei J, Wellcome Trust Case Control Consortium, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464: 704-712. 10.1038/nature08516.
Perry GH, Tchinda J, McGrath SD, Zhang J, Picker SR, Cáceres AM, Iafrate AJ, Tyler-Smith C, Scherer SW, Eichler EE, Stone AC, Lee C: Hotspots for copy number variation in chimpanzees and humans. Proc Natl Acad Sci USA. 2006, 103: 8006-8011. 10.1073/pnas.0602318103.
Perry GH, Yang F, Marques-Bonet T, Murphy C, Fitzgerald T, Lee AS, Hyland C, Stone AC, Hurles ME, Tyler-Smith C, Eichler EE, Carter NP, Lee C, Redon R: Copy number variation and evolution in humans and chimpanzees. Genome Res. 2008, 18: 1698-1710. 10.1101/gr.082016.108.
Lee AS, Gutiérrez-Arcelus M, Perry GH, Vallender EJ, Johnson WE, Miller GM, Korbel JO, Lee C: Analysis of copy number variation in the rhesus macaque genome identifies candidate loci for evolutionary and human disease studies. Hum Mol Genet. 2008, 17: 1127-1136. 10.1093/hmg/ddn002.
Adams DJ, Dermitzakis ET, Cox T, Smith J, Davies R, Banerjee R, Bonfield J, Mullikin JC, Chung YJ, Rogers J, Bradley A: Complex haplotypes, copy number polymorphisms and coding variation in two recently divergent mouse strains. Nat Genet. 2005, 37: 532-536. 10.1038/ng1551.
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.
Egan CM, Sridhar S, Wigler M, Hall IM: Recurrent DNA copy number variation in the laboratory mouse. Nat Genet. 2007, 39: 1384-1389. 10.1038/ng.2007.19.
Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD: Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007, 17: 1743-1754. 10.1101/gr.6754607.
She X, Cheng Z, Zöllner S, Church DM, Eichler EE: Mouse segmental duplication and copy number variation. Nat Genet. 2008, 40: 909-914. 10.1038/ng.172.
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.
Charchar FJ, Kaiser M, Bingham AJ, Fotinatos N, Ahmady F, Tomaszewski M, Samani NJ: Whole genome survey of copy number variation in the spontaneously hypertensive rat: relationship to quantitative trait loci, gene expression, and blood pressure. Hypertension. 2010, 55: 1231-1238. 10.1161/HYPERTENSIONAHA.109.141663.
Chen WK, Swartz JD, Rush LJ, Alvarez CE: Mapping DNA structural variation in dogs. Genome Res. 2009, 19: 500-509. 10.1101/gr.083741.108.
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. 10.1101/gr.084715.108.
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.
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.
Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O'Connell J, Moore SS, Smith TP, Sonstegard TS, Van Tassell CP: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4: e5350-10.1371/journal.pone.0005350.
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.
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.
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.
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.
Cahan P, Li Y, Izumi M, Graubert TA: The impact of copy number variation on local gene expression in mouse hematopoietic stem and progenitor cells. Nat Genet. 2009, 41: 430-437. 10.1038/ng.350.
Henrichsen CN, Vinckenbosch N, Zöllner S, Chaignat E, Pradervand S, Schütz F, Ruedi M, Kaessmann H, Reymond A: Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009, 41: 424-429. 10.1038/ng.345.
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.
Zhang F, Gu W, Hurles ME, Lupski JR: Copy number variation in human health, disease, and evolution. Annu Rev Genomics Hum Genet. 2009, 10: 451-481. 10.1146/annurev.genom.9.081307.164217.
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.
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.
Fontanesi L, D'Alessandro E, Scotti E, Liotta L, Crovetti A, Chiofalo V, Russo V: Genetic heterogeneity and selection signature at the KIT gene in pigs showing different coat colours and patterns. Anim Genet. 2010, 41: 478-492. 10.1111/j.1365-2052.2010.02054.x.
Norris BJ, Whan VA: A gene duplication affecting expression of the ovine ASIP gene is responsible for white and black sheep. Genome Res. 2008, 18: 1282-1293. 10.1101/gr.072090.107.
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.
Elferink MG, Vallée AA, Jungerius AP, Crooijmans RP, Groenen MA: Partial duplication of the PRLR and SPEF2 genes at the late feathering locus in chicken. BMC Genomics. 2008, 9: 391-10.1186/1471-2164-9-391.
Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Kowbel D, Collins C, Kuo WL, Chen C, Zhai Y, Dairkee SH, Ljung BM, Gray JW, Albertson DG: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays. Nat Genet. 1998, 20: 207-211. 10.1038/2524.
Ishkanian AS, Malloff CA, Watson SK, DeLeeuw RJ, Chi B, Coe BP, Snijders A, Albertson DG, Pinkel D, Marra MA, Ling V, MacAulay C, Lam WL: A tiling resolution DNA microarray with complete coverage of the human genome. Nat Genet. 2004, 36: 299-303. 10.1038/ng1307.
Korbel JO, Urban AE, Affourtit JP, Godwin B, Grubert F, Simons JF, Kim PM, Palejev D, Carriero NJ, Du L, Taillon BE, Chen Z, Tanzer A, Saunders AC, Chi J, Yang F, Carter NP, Hurles ME, Weissman SM, Harkins TT, Gerstein MB, Egholm M, Snyder M: Paired-end mapping reveals extensive structural variation in the human genome. Science. 2007, 318: 420-426. 10.1126/science.1149504.
Carter NP: Methods and strategies for analyzing copy number variation using DNA microarrays. Nat Genet. 2007, 39: S16-S21. 10.1038/ng2028.
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.
Park H, Kim JI, Ju YS, Gokcumen O, Mills RE, Kim S, Lee S, Suh D, Hong D, Kang HP, Yoo YJ, Shin JY, Kim HJ, Yavartanoo M, Chang YW, Ha JS, Chong W, Hwang GR, Darvishi K, Kim H, Yang SJ, Yang KS, Kim H, Hurles ME, Scherer SW, Carter NP, Tyler-Smith C, Lee C, Seo JS: Discovery of common Asian copy number variants using integrated high-resolution array CGH and massively parallel DNA sequencing. Nat Genet. 2010, 42: 400-405. 10.1038/ng.555.
Dumas L, Kim YH, Karimpour-Fard A, Cox M, Hopkins J, Pollack JR, Sikela JM: Gene copy number variation spanning 60 million years of human and primate evolution. Genome Res. 2007, 17: 1266-1277. 10.1101/gr.6557307.
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.
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.
Völker M, Backström N, Skinner BM, Langley EJ, Bunzey SK, Ellegren H, Griffin DK: Copy number variation, chromosome rearrangement, and their association with recombination during avian evolution. Genome Res. 2010, 20: 503-511. 10.1101/gr.103663.109.
Bovine Genome Sequencing and Analysis Consortium: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324: 522-528. 10.1126/science.1169588.
Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, Hanrahan F, Pertea G, Van Tassell CP, Sonstegard TS, Marçais G, Roberts M, Subramanian P, Yorke JA, Salzberg SL: A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009, 10: R42-10.1186/gb-2009-10-4-r42.
Hayes H, Elduque C, Gautier M, Schibler L, Cribiu E, Eggen A: Mapping of 195 genes in cattle and updated comparative map with man, mouse, rat and pig. Cytogenet Genome Res. 2003, 102: 16-24. 10.1159/000075719.
Lai W, Choudhary V, Park PJ: CGHweb: a tool for comparing DNA copy number segmentations from multiple algorithms. Bioinformatics. 2008, 24: 1014-1015. 10.1093/bioinformatics/btn067.
Wang X, Nahashon S, Feaster TK, Bohannon-Stewart A, Adefope N: An initial map of chromosomal segmental copy number variations in the chicken. BMC Genomics. 2010, 11: 351-10.1186/1471-2164-11-351.
Turner DJ, Miretti M, Rajan D, Fiegler H, Carter NP, Blayney ML, Beck S, Hurles ME: Germline rates of de novo meiotic deletions and duplications causing several genomic disorders. Nat Genet. 2008, 40: 90-95. 10.1038/ng.2007.40.
Lalueza-Fox C, Castresana J, Sampietro L, Marquès-Bonet T, Alcover JA, Bertranpetit J: Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal. BMC Evol Biol. 2005, 5: 70-10.1186/1471-2148-5-70.
Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392: 917-920. 10.1038/31927.
Steiper ME, Young NM, Sukarna TY: Genomic data support the hominoid slowdown and an Early Oligocene estimate for the hominoid-cercopithecoid divergence. Proc Natl Acad Sci USA. 2004, 101: 17021-17026. 10.1073/pnas.0407270101.
Kim PM, Lam HY, Urban AE, Korbel JO, Affourtit J, Grubert F, Chen X, Weissman S, Snyder M, Gerstein MB: Analysis of copy number variants and segmental duplications in the human genome: Evidence for a change in the process of formation in recent evolutionary history. Genome Res. 2008, 18: 1865-1874. 10.1101/gr.081422.108.
Geller DS, Rodriguez-Soriano J, Vallo Boado A, Schifter S, Bayer M, Chang SS, Lifton RP: Mutations in the mineralocorticoid receptor gene cause autosomal dominant pseudohypoaldosteronism type I. Nat Genet. 1998, 19: 279-281. 10.1038/966.
Martinez F, Mansego ML, Escudero JC, Redon J, Chaves FJ: Association of a mineralocorticoid receptor gene polymorphism with hypertension in a Spanish population. Am J Hypertens. 2009, 22: 649-655. 10.1038/ajh.2009.39.
Rozeboom AM, Akil H, Seasholtz AF: Mineralocorticoid receptor overexpression in forebrain decreases anxiety-like behavior and alters the stress response in mice. Proc Natl Acad Sci USA. 2007, 104: 4688-4693. 10.1073/pnas.0606067104.
Nguyen DQ, Webber C, Ponting CP: Bias of selection on human copy-number variants. PLoS Genet. 2006, 2: e20-10.1371/journal.pgen.0020020.
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.
Naderi S, Rezaei HR, Pompanon F, Blum MG, Negrini R, Naghash HR, Balkiz O, Mashkour M, Gaggiotti OE, Ajmone-Marsan P, Kence A, Vigne JD, Taberlet P: The goat domestication process inferred from large-scale mitochondrial DNA analysis of wild and domestic individuals. Proc Natl Acad Sci USA. 2008, 105: 17659-17664. 10.1073/pnas.0804782105.
Dukkipati VS, Blair HT, Garrick DJ, Murray A: Ovar-Mhc - ovine major histocompatibility complex: role in genetic resistance to diseases. N Z Vet J. 2006, 54: 153-160.
Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, Vivanco I, Lee JC, Huang JH, Alexander S, Du J, Kau T, Thomas RK, Shah K, Soto H, Perner S, Prensner J, Debiasi RM, Demichelis F, Hatton C, Rubin MA, Garraway LA, Nelson SF, Liau L, Mischel PS, Cloughesy TF, Meyerson M, Golub TA, Lander ES, Mellinghoff IK, Sellers WR: Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci USA. 2007, 104: 20007-20012. 10.1073/pnas.0710052104.
Le Maréchal C, Masson E, Chen JM, Morel F, Ruszniewski P, Levy P, Férec C: Hereditary pancreatitis caused by triplication of the trypsinogen locus. Nat Genet. 2006, 38: 1372-1374. 10.1038/ng1904.
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.
We thank Drs. Charles Lee, Omer Gökçümen (Harvard Medical School, USA), and Giuliano Galimberti (University of Bologna, Italy) for advices on statistical evaluations of goat-cattle CNVR overlaps. We thank farmers and people who helped in obtaining goat samples. We also thank the three anonymous reviewers for their useful comments to the early version of the manuscript. This work was funded by the Assessorato Agricoltura e Foreste of the Regione Siciliana - U.O.B. 108, SOAT n. 69 Aragona (AG), by the Italian MiPAAF SELMOL project and by the University of Bologna RFO funds and was associated with the MIUR-PNR 2003 project (FIRB art.8) termed LIBI-Laboratorio Internazionale di BioInformatica.
LF conceived and designed the study, analysed data, contributed to the sampling, coordinated and organized the laboratory work, drafted the manuscript and partially funded the study. PLM conducted bioinformatic and computational analyses. FB carried out the laboratory work and prepared data for analyses. VRi and SD contributed to the sampling and data. MC collaborated in the laboratory work. RC and VRu supervised the work and were involved in the design of the study. BP provided samples, supervised, took part in designing the work and partially funded the study. All authors reviewed the manuscript and accepted the final version.