Genome-wide copy number variant discovery in dogs using the CanineHD genotyping array

Detection of InDel and CNV of SPAG17 gene and their associations with bovine growth traits

Sperm-associated antigen 17 (SPAG17) gene encodes a central pair protein, which is involved in flagellar motility, male fertility and skeletal growth in ruminants. The insertions/deletions (indels) and copy number variations (CNVs) influence phenotypic traits by altering the sequences and copy numbers of functional genes, respectively. This study identified a novel 8-bp indel of SPAG17 gene in 1520 individuals from eight different cattle breeds, as well as a novel CNV region in 355 animals. The correlation analysis of indel showed that the individuals of ID genotype had superior performance traits such as body height (p = 0.038) and body slanting length (p = 0.041) as compared to other genotypes in Xianan cattle. For the CNV, different copy numbers were closely related to the body height in Qinchuan (p = 0.045) and body weight in Xianan (p = 0.036) breeds. Importantly, significant difference was observed between the 8-bp indel and the copy number loss in Xianan breed (p < 0.01). These findings indicated that the variations within the bovine SPAG17 gene can be considered as an effective DNA molecular marker for beef cattle breeding.

Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds

Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

A novel canine reference genome resolves genomic architecture and uncovers transcript complexity

We present GSD_1.0, a high-quality domestic dog reference genome with chromosome length scaffolds and contiguity increased 55-fold over CanFam3.1. Annotation with generated and existing long and short read RNA-seq, miRNA-seq and ATAC-seq, revealed that 32.1% of lifted over CanFam3.1 gaps harboured previously hidden functional elements, including promoters, genes and miRNAs in GSD_1.0. A catalogue of canine "dark" regions was made to facilitate mapping rescue. Alignment in these regions is difficult, but we demonstrate that they harbour trait-associated variation. Key genomic regions were completed, including the Dog Leucocyte Antigen (DLA), T Cell Receptor (TCR) and 366 COSMIC cancer genes. 10x linked-read sequencing of 27 dogs (19 breeds) uncovered 22.1 million SNPs, indels and larger structural variants. Subsequent intersection with protein coding genes showed that 1.4% of these could directly influence gene products, and so provide a source of normal or aberrant phenotypic modifications.

Polymorphisms of AMY1A gene and their association with growth, carcass traits and feed intake efficiency in chickens

Investigations on the association between chicken traits and genetic variations can provide basic information to improve production performance in chickens. In our previous work, we genotyped 450 male chickens with a 600 K SNP array [1] and found that several SNPs in the genomic regions of the amylase alpha 1A (AMY1A) gene were significantly associated with feed intake efficiency and carcass traits. Given the lower accuracy of the SNP array, we performed direct sequencing with male and female chickens to further test chicken AMY1A polymorphisms and investigate their association with 17 traits in chickens. The results showed that 7 SNPs in the 5' flanking region, exon, intron and 3' UTR (3' untranslated region) of AMY1A, were significantly associated with daily gain (DG), average daily feed intake (ADFI), leg muscle weight (LMW) and abdominal fat (AF) (p < 0.05). Additionally, the haplotypes based on three SNPs, rs15910189, rs314354067 and rs316026696, showed significant associations with DG (p < 0.01), ADFI and AF (p < 0.05). To better understand the transcriptional regulation of AMY1A, we cloned its 5' flanking region and found that the SNPs rs316436216 and rs314213090 which might change the transcriptional regulator binding sites, were in the suppressor and enhancer regions, respectively. In addition, luciferase assays revealed that the SNP rs314613110 in the 3' UTR influenced the binding of the miRNA gga-miR-1764-3p. To validate whether there is any copy number variation in AMY1A in our population, we performed a genome-wide assessment of CNVs through whole-genome resequencing data. However, no CNV was found in AMY1A in our population, which is different from the increased copy number of AMY1A found in humans who consume a high-starch diet. Therefore, the present study provides substantial evidence for the association of AMY1A polymorphisms with growth traits and feed intake efficiency, which might contribute to chicken breeding programs.

Analysis of copy number variation in dogs implicates genomic structural variation in the development of anterior cruciate ligament rupture

Anterior cruciate ligament (ACL) rupture is an important condition of the human knee. Second ruptures are common and societal costs are substantial. Canine cranial cruciate ligament (CCL) rupture closely models the human disease. CCL rupture is common in the Labrador Retriever (5.79% prevalence), ~100-fold more prevalent than in humans. Labrador Retriever CCL rupture is a polygenic complex disease, based on genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers. Dissection of genetic variation in complex traits can be enhanced by studying structural variation, including copy number variants (CNVs). Dogs are an ideal model for CNV research because of reduced genetic variability within breeds and extensive phenotypic diversity across breeds. We studied the genetic etiology of CCL rupture by association analysis of CNV regions (CNVRs) using 110 case and 164 control Labrador Retrievers. CNVs were called from SNPs using three different programs (PennCNV, CNVPartition, and QuantiSNP). After quality control, CNV calls were combined to create CNVRs using ParseCNV and an association analysis was performed. We found no strong effect CNVRs but found 46 small effect (max(T) permutation P<0.05) CCL rupture associated CNVRs in 22 autosomes; 25 were deletions and 21 were duplications. Of the 46 CCL rupture associated CNVRs, we identified 39 unique regions. Thirty four were identified by a single calling algorithm, 3 were identified by two calling algorithms, and 2 were identified by all three algorithms. For 42 of the associated CNVRs, frequency in the population was <10% while 4 occurred at a frequency in the population ranging from 10–25%. Average CNVR length was 198,872bp and CNVRs covered 0.11 to 0.15% of the genome. All CNVRs were associated with case status. CNVRs did not overlap previous canine CCL rupture risk loci identified by GWAS. Associated CNVRs contained 152 annotated genes; 12 CNVRs did not have genes mapped to CanFam3.1. Using pathway analysis, a cluster of 19 homeobox domain transcript regulator genes was associated with CCL rupture (P = 6.6E-13). This gene cluster influences cranial-caudal body pattern formation during embryonic limb development. Clustered genes were found in 3 CNVRs on chromosome 14 (HoxA), 28 (NKX6-2), and 36 (HoxD). When analysis was limited to deletion CNVRs, the association was strengthened (P = 8.7E-16). This study suggests a component of the polygenic risk of CCL rupture in Labrador Retrievers is associated with small effect CNVs and may include aspects of stifle morphology regulated by homeobox domain transcript regulator genes.

The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics

The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.

Whole genome resequencing of the Iranian native dogs and wolves to unravel variome during dog domestication

BACKGROUND

Advances in genome technology have simplified a new comprehension of the genetic and historical processes crucial to rapid phenotypic evolution under domestication. To get new insight into the genetic basis of the dog domestication process, we conducted whole-genome sequence analysis of three wolves and three dogs from Iran which covers the eastern part of the Fertile Crescent located in Southwest Asia where the independent domestication of most of the plants and animals has been documented and also high haplotype sharing between wolves and dog breeds has been reported.

RESULTS

Higher diversity was found within the wolf genome compared with the dog genome. A total number of 12.45 million SNPs were detected in all individuals (10.45 and 7.82 million SNPs were identified for all the studied wolves and dogs, respectively) and a total number of 3.49 million small Indels were detected in all individuals (3.11 and 2.24 million small Indels were identified for all the studied wolves and dogs, respectively). A total of 10,571 copy number variation regions (CNVRs) were detected across the 6 individual genomes, covering 154.65 Mb, or 6.41%, of the reference genome (canFam3.1). Further analysis showed that the distribution of deleterious variants in the dog genome is higher than the wolf genome. Also, genomic annotation results from intron and intergenic regions showed that the proportion of variations in the wolf genome is higher than that in the dog genome, while the proportion of the coding sequences and 3'-UTR in the dog genome is higher than that in the wolf genome. The genes related to the olfactory and immune systems were enriched in the set of the structural variants (SVs) identified in this work.

CONCLUSIONS

Our results showed more deleterious mutations and coding sequence variants in the domestic dog genome than those in wolf genome. By providing the first Iranian dog and wolf variome map, our findings contribute to understanding the genetic architecture of the dog domestication.

Inter- and intra-breed genome-wide copy number diversity in a large cohort of European equine breeds

BACKGROUND

Copy Number Variation (CNV) is a common form of genetic variation underlying animal evolution and phenotypic diversity across a wide range of species. In the mammalian genome, high frequency of CNV differentiation between breeds may be candidates for population-specific selection. However, CNV differentiation, selection and its population genetics have been poorly explored in horses.

RESULTS

We investigated the patterns, population variation and gene annotation of CNV using the Axiom® Equine Genotyping Array (670,796 SNPs) from a large cohort of individuals (N = 1755) belonging to eight European horse breeds, varying from draught horses to several warmblood populations. After quality control, 152,640 SNP CNVs (individual markers), 18,800 segment CNVs (consecutive SNP CNVs of same gain/loss state or both) and 939 CNV regions (CNVRs; overlapping segment CNVs by at least 1 bp) compared to the average signal of the reference (Belgian draught horse) were identified. Our analyses showed that Equus caballus chromosome 12 (ECA12) was the most enriched in segment CNV gains and losses (~ 3% average proportion of the genome covered), but the highest number of segment CNVs were detected on ECA1 and ECA20 (regardless of size). The Friesian horses showed private SNP CNV gains (> 20% of the samples) on ECA1 and Exmoor ponies displayed private SNP CNV losses on ECA25 (> 20% of the samples). The Warmblood cluster showed private SNP CNV gains located in ECA9 and Draught cluster showed private SNP CNV losses located in ECA7. The length of the CNVRs ranged from 1 kb to 21.3 Mb. A total of 10,612 genes were annotated within the CNVRs. The PANTHER annotation of these genes showed significantly under- and overrepresented gene ontology biological terms related to cellular processes and immunity (Bonferroni P-value < 0.05). We identified 80 CNVRs overlapping with known QTL for fertility, coat colour, conformation and temperament. We also report 67 novel CNVRs.

CONCLUSIONS

This work revealed that CNV patterns, in the genome of some European horse breeds, occurred in specific genomic regions. The results provide support to the hypothesis that high frequency private CNVs residing in genes may potentially be responsible for the diverse phenotypes seen between horse breeds.

A Genome-Wide Detection of Copy Number Variations Using SNP Genotyping Arrays in Braque Français Type Pyrénées Dogs

Copy number variants (CNVs) are an important source of genetic variation complementary to single nucleotide polymorphisms (SNPs). Only few studies have been conducted in dogs on CNVs derived from high-density SNP array data, and many canine breeds still remain uncharacterized, e.g., the Braque Français, type Pyrénées breed (BRA). Therefore, in an effort to more comprehensively investigate the canine genome for CNVs, we used a high-density SNP array (170 K) to discover CNVs in BRA. The CNV regions (CNVRs) were identified through the merging of two different CNVRs datasets, obtained separately from SNP data using the PennCNV and SVS software. A total of 45 stringent CNVRs, ranging from 3.5 kb to 458,716 kb in length were detected in 26 dog samples. Results overlapped moderately in comparison with previous studies on CNVs in dogs, leading to the identification of 16 novel CNVRs. A total of 159 genes were annotated in the CNVRs detected with stringent quality criteria in particular high classification stringency and false discovery rate correction. The gene ontology enrichment analysis provided information on biological processes and cellular components related to muscle structure development and muscle cell differentiation. Considering that BRA is a breed used for speed in hunting and retrieval, for the ability to find feathered game, and for pointing, we can hypothesize that selection for such hunting behavior could have driven, at least in part, the presence of these genes into the CNVRs.

Copy number variation and variant discovery in Bullmastiff dogs

Identification of genomic variants within dogs is important for understanding genetic factors contributing to breed diversity and phenotypic traits. This study aimed to identify sources of variation in the Bullmastiff using high-density signal intensity and whole-genome sequence data. Close to 3000 copy number variants (CNVs) were identified in Bullmastiff dogs using Canine HD BeadChip data. When CNVs were collated, 82 CNV regions (CNVRs) were detected, 50% in transcribed regions encompassing 432 genes. Fifty of the CNVRs detected have not been reported in other breeds and represent potential breed-specific variants. A proportion of the CNVR variants with predicted modifying effects on gene pathways may contribute to breed traits. Approximately 5 million putative variants per dog, inclusive of single nucleotide polymorphisms (SNPs), multi-nucleotide polymorphisms (MNPs) and insertion and deletions (INDELs), were identified from DNA sequence data on a small number of animals. Identification of genetic variants in the Bullmastiff highlights sources of variation in the breed and molecular markers that will assist in future trait and disease investigations in dogs.

First genome-wide CNV mapping in FELIS CATUS using next generation sequencing data

Background

Copy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium.

Results

Reads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated.

Conclusion

This study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5297-2) contains supplementary material, which is available to authorized users.

Genetic Diversity of Seven Cattle Breeds Inferred Using Copy Number Variations

Copy number variations (CNVs) comprise deletions, duplications, and insertions found within the genome larger than 50 bp in size. CNVs are thought to be primary role-players in breed formation and adaptation. South Africa boasts a diverse ecology with harsh environmental conditions and a broad spectrum of parasites and diseases that pose challenges to livestock production. This has led to the development of composite cattle breeds which combine the hardiness of Sanga breeds and the production potential of the Taurine breeds. The prevalence of CNVs within these respective breeds of cattle and the prevalence of CNV regions (CNVRs) in their diversity, adaptation and production is however not understood. This study therefore aimed to ascertain the prevalence, diversity, and correlations of CNVRs within cattle breeds used in South Africa. Illumina Bovine SNP50 data and PennCNV were utilized to identify CNVRs within the genome of 287 animals from seven cattle breeds representing Sanga, Taurine, Composite, and cross breeds. Three hundred and fifty six CNVRs of between 36 kb to 4.1 Mb in size were identified. The null hypothesis that one CNVR loci is independent of another was tested using the GENEPOP software. One hunded and two and seven of the CNVRs in the Taurine and Sanga/Composite cattle breeds demonstrated a significant (p ≤ 0.05) association. PANTHER overrepresentation analyses of correlated CNVRs demonstrated significant enrichment of a number of biological processes, molecular functions, cellular components, and protein classes. CNVR genetic variation between and within breed group was measured using phiPT which allows intra-individual variation to be suppressed and hence proved suitable for measuring binary CNVR presence/absence data. Estimate PhiPT within and between breed variance was 2.722 and 0.518 respectively. Pairwise population PhiPT values corresponded with breed type, with Taurine Holstein and Angus breeds demonstrating no between breed CNVR variation. Phylogenetic trees were drawn. CNVRs primarily clustered animals of the same breed type together. This study successfully identified, characterized, and analyzed 356 CNVRs within seven cattle breeds. CNVR correlations were evident, with many more correlations being present among the exotic Taurine breeds. CNVR genetic diversity of Sanga, Taurine and Composite breeds was ascertained with breed types exposed to similar selection pressures demonstrating analogous incidences of CNVRs.

Mapping of genome-wide copy number variations in the Iranian indigenous cattle using a dense SNP data set

The objective of this study was to present the first map of the copy number variations (CNVs) in Iranian indigenous cattle based on a high-density single nucleotide polymorphism (SNP) dataset. A total of 90 individuals were genotyped using the Illumina BovineHD BeadChip containing 777 962 SNPs. The QuantiSNP algorithm was used to perform a genome-wide CNV detection across autosomal genome. After merging the overlapping CNV, a total of 221 CNV regions were identified encompassing 36.4 Mb or 1.44% of the bovine autosomal genome. The length of the CNV regions ranged from 3.5 to 2252.8 Kb with an average of 163.8 Kb. These regions included 147 loss (66.52%) and 74 gain (33.48%) events containing a total of 637 annotated Ensembl genes. Gene ontology analysis revealed that most of genes in the CNV regions were involved in environmental responses, disease susceptibility and immune system functions. Furthermore, 543 of these genes corresponded to the human orthologous genes, which involved in a wide range of biological functions. Altogether, 73% of the 221 CNV regions overlapped either completely or partially with those previously reported in other cattle studies. Moreover, novel CNV regions involved several quantitative trait loci (QTL)-related to adaptative traits of Iranian indigenous cattle. These results provided a basis to conduct future studies on association between CNV regions and phenotypic variations in the Iranian indigenous cattle.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing

BACKGROUND

Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids.

RESULTS

Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts.

CONCLUSIONS

This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

Gene expression patterns of chicken neuregulin 3 in association with copy number variation and frameshift deletion

Background

Neuregulin 3 (NRG3) plays a key role in central nervous system development and is a strong candidate for human mental disorders. Thus, genetic variation in NRG3 may have some impact on a variety of phenotypes in non-mammalian vertebrates. Recently, genome-wide screening for short insertions and deletions in chicken (Gallus gallus) genomes has provided useful information about structural variation in functionally important genes. NRG3 is one such gene that has a putative frameshift deletion in exon 2, resulting in premature termination of translation. Our aims were to characterize the structure of chicken NRG3 and to compare expression patterns between NRG3 isoforms.

Results

Depending on the presence or absence of the 2-bp deletion in chicken NRG3, 3 breeds (red junglefowl [RJF], Boris Brown [BB], and Hinai-jidori [HJ]) were genotyped using flanking primers. In the commercial breeds (BB and HJ), approximately 45% of individuals had at least one exon 2 allele with the 2-bp deletion, whereas there was no deletion allele in RJF. The lack of a homozygous mutant indicated the existence of duplicated NRG3 segments in the chicken genome. Indeed, highly conserved elements consisting of exon 1, intron 1, exon 2, and part of intron 2 were found in the reference RJF genome, and quantitative PCR detected copy number variation (CNV) between breeds as well as between individuals. The copy number of conserved elements was significantly higher in chicks harboring the 2-bp deletion in exon 2. We identified 7 novel transcript variants using total mRNA isolated from the amygdala. Novel isoforms were found to lack the exon 2 cassette, which probably harbored the premature termination codon. The relative transcription levels of the newly identified isoforms were almost the same between chick groups with and without the 2-bp deletion, while chicks with the deletion showed significant suppression of the expression of previously reported isoforms.

Conclusions

A putative frameshift deletion and CNV in chicken NRG3 are structural mutations that occurred before the establishment of commercial chicken lines. Our results further suggest that the putative frameshift deletion in exon 2 may potentially affect the expression level of particular isoforms of chicken NRG3.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0537-z) contains supplementary material, which is available to authorized users.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits

Background

Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance.

Results

Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively.

Conclusions

Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2752-9) contains supplementary material, which is available to authorized users.

Absolute quantification reveals the stable transmission of a high copy number variant linked to autoinflammatory disease

BACKGROUND

Dissecting the role copy number variants (CNVs) play in disease pathogenesis is directly reliant on accurate methods for quantification. The Shar-Pei dog breed is predisposed to a complex autoinflammatory disease with numerous clinical manifestations. One such sign, recurrent fever, was previously shown to be significantly associated with a novel, but unstable CNV (CNV_16.1). Droplet digital PCR (ddPCR) offers a new mechanism for CNV detection via absolute quantification with the promise of added precision and reliability. The aim of this study was to evaluate ddPCR in relation to quantitative PCR (qPCR) and to assess the suitability of the favoured method as a genetic test for Shar-Pei Autoinflammatory Disease (SPAID).

RESULTS

One hundred and ninety-six individuals were assayed using both PCR methods at two CNV positions (CNV_14.3 and CNV_16.1). The digital method revealed a striking result. The CNVs did not follow a continuum of alleles as previously reported, rather the alleles were stable and pedigree analysis showed they adhered to Mendelian segregation. Subsequent analysis of ddPCR case/control data confirmed that both CNVs remained significantly associated with the subphenotype of fever, but also to the encompassing SPAID complex (p < 0.001). In addition, harbouring CNV_16.1 allele five (CNV_16.1|5) resulted in a four-fold increase in the odds for SPAID (p < 0.001). The inclusion of a genetic marker for CNV_16.1 in a genome-wide association test revealed that this variant explained 9.7 % of genetic variance and 25.8 % of the additive genetic heritability of this autoinflammatory disease.

CONCLUSIONS

This data shows the utility of the ddPCR method to resolve cryptic copy number inheritance patterns and so open avenues of genetic testing. In its current form, the ddPCR test presented here could be used in canine breeding to reduce the number of homozygote CNV_16.1|5 individuals and thereby to reduce the prevalence of disease in this breed.

Genome-wide identification of copy number variation using high-density single-nucleotide polymorphism array in Japanese Black cattle

Background

Copy number variation (CNV) is an important source of genetic variability associated with phenotypic variation and disease susceptibility. Comprehensive genome-wide CNV maps provide valuable information for genetic and functional studies. To identify CNV in Japanese Black cattle, we performed a genome-wide autosomal screen using genomic data from 1,481 animals analyzed with the Illumina Bovine High-Density (HD) BeadChip Array (735,293 single-nucleotide polymorphisms (SNPs) with an average marker interval of 3.4 kb on the autosomes).

Results

We identified a total of 861 CNV regions (CNVRs) across all autosomes, which covered 43.65 Mb of the UMD3.1 genome assembly and corresponded to 1.74 % of the 29 bovine autosomes. Overall, 35 % of the CNVRs were present at a frequency of > 1 % in 1,481 animals. The estimated lengths of CNVRs ranged from 1.1 kb to 1.4 Mb, with an average of 50.7 kb. The average number of CNVR events per animal was 35. Comparisons with previously reported cattle CNV showed that 72 % of the CNVR calls detected in this study were within or overlapped with known CNVRs. Experimentally, three CNVRs were validated using quantitative PCR, and one CNVR was validated using PCR with flanking primers for the deleted region. Out of the 861 CNVRs, 390 contained 717 Ensembl-annotated genes significantly enriched for stimulus response, cellular defense response, and immune response in the Gene Ontology (GO) database. To associate genes contained in CNVRs with phenotypes, we converted 560 bovine Ensembl gene IDs to their 438 orthologous associated mouse gene IDs, and 195 of these mouse orthologous genes were categorized into 1,627 phenotypes in the Mouse Genome Informatics (MGI) database.

Conclusions

We identified 861 CNVRs in 1,481 Japanese Black cattle using the Illumina BovineHD BeadChip Array. The genes contained in CNVRs were characterized using GO analysis and the mouse orthologous genes were characterized using the MGI database. The comprehensive genome-wide CNVRs map will facilitate identification of genetic variation and disease-susceptibility alleles in Japanese Black cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0335-z) contains supplementary material, which is available to authorized users.

Copy number variation in the region harboring SOX9 gene in dogs with testicular/ovotesticular disorder of sex development (78,XX; SRY-negative)

Although the disorder of sex development in dogs with female karyotype (XX DSD) is quite common, its molecular basis is still unclear. Among mutations underlying XX DSD in mammals are duplication of a long sequence upstream of the SOX9 gene (RevSex) and duplication of the SOX9 gene (also observed in dogs). We performed a comparative analysis of 16 XX DSD and 30 control female dogs, using FISH and MLPA approaches. Our study was focused on a region harboring SOX9 and a region orthologous to the human RevSex (CanRevSex), which was located by in silico analysis downstream of SOX9. Two highly polymorphic copy number variable regions (CNVRs): CNVR1 upstream of SOX9 and CNVR2 encompassing CanRevSex were identified. Although none of the detected copy number variants were specific to either affected or control animals, we observed that the average number of copies in CNVR1 was higher in XX DSD. No copy variation of SOX9 was observed. Our extensive studies have excluded duplication of SOX9 as the common cause of XX DSD in analyzed samples. However, it remains possible that the causative mutation is hidden in highly polymorphic CNVR1.

Copy number variation in the horse genome

We constructed a 400K WG tiling oligoarray for the horse and applied it for the discovery of copy number variations (CNVs) in 38 normal horses of 16 diverse breeds, and the Przewalski horse. Probes on the array represented 18,763 autosomal and X-linked genes, and intergenic, sub-telomeric and chrY sequences. We identified 258 CNV regions (CNVRs) across all autosomes, chrX and chrUn, but not in chrY. CNVs comprised 1.3% of the horse genome with chr12 being most enriched. American Miniature horses had the highest and American Quarter Horses the lowest number of CNVs in relation to Thoroughbred reference. The Przewalski horse was similar to native ponies and draft breeds. The majority of CNVRs involved genes, while 20% were located in intergenic regions. Similar to previous studies in horses and other mammals, molecular functions of CNV-associated genes were predominantly in sensory perception, immunity and reproduction. The findings were integrated with previous studies to generate a composite genome-wide dataset of 1476 CNVRs. Of these, 301 CNVRs were shared between studies, while 1174 were novel and require further validation. Integrated data revealed that to date, 41 out of over 400 breeds of the domestic horse have been analyzed for CNVs, of which 11 new breeds were added in this study. Finally, the composite CNV dataset was applied in a pilot study for the discovery of CNVs in 6 horses with XY disorders of sexual development. A homozygous deletion involving AKR1C gene cluster in chr29 in two affected horses was considered possibly causative because of the known role of AKR1C genes in testicular androgen synthesis and sexual development. While the findings improve and integrate the knowledge of CNVs in horses, they also show that for effective discovery of variants of biomedical importance, more breeds and individuals need to be analyzed using comparable methodological approaches.

Genomes of individuals in a species vary in many ways, one of which is DNA copy number variation (CNV). This includes deletions, duplications, and complex rearrangements typically larger than 50 base-pairs. CNVs are part of normal genetic variation contributing to phenotypic diversity but can also be pathogenic and associated with diseases and disorders. In order to distinguish between the two, detailed knowledge about CNVs in the species of interest is needed. Here we studied the genomes of 38 normal horses of 16 diverse breeds, and identified 258 CNV regions. We integrated our findings with previously published horse CNVs and generated a composite dataset of ∼1400 CNVRs. Despite this large number, our analysis shows that CNV research in horses needs further improvement because the current data are based on 10% of horse breeds and that most CNVRs are study-specific and require validation. Finally, we analyzed CNVs in horses with disorders of sexual development and found in two male pseudo-hermaphrodites a large deletion disrupting a group of genes involved in sex hormone metabolism and sexual differentiation. The findings underline the possible role of CNVs in complex disorders such as development and reproduction.