Analysis of copy number variations in Holstein cows identify potential mechanisms contributing to differences in residual feed intake

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. In this study, we performed an initial analysis of copy number variations (CNVs) using BovineHD SNP genotyping data from 147 Holstein cows identified as having high or low feed efficiency as estimated by residual feed intake (RFI). We detected 443 candidate CNV regions (CNVRs) that represent 18.4?Mb (0.6?%) of the genome. To investigate the functional impacts of CNVs, we created two groups of 30 individual animals with extremely low or high estimated breeding values (EBVs) for RFI, and referred to these groups as low intake (LI; more efficient) or high intake (HI; less efficient), respectively. We identified 240 (~9.0?Mb) and 274 (~10.2?Mb) CNVRs from LI and HI groups, respectively. Approximately 30–40?% of the CNVRs were specific to the LI group or HI group of animals. The 240 LI CNVRs overlapped with 137 Ensembl genes. Network analyses indicated that the LI-specific genes were predominantly enriched for those functioning in the inflammatory response and immunity. By contrast, the 274 HI CNVRs contained 177 Ensembl genes. Network analyses indicated that the HI-specific genes were particularly involved in the cell cycle, and organ and bone development. These results relate CNVs to two key variables, namely immune response and organ and bone development. The data indicate that greater feed efficiency relates more closely to immune response, whereas cattle with reduced feed efficiency may have a greater capacity for organ and bone development.     

Fine mapping of copy number variations on two cattle genome assemblies using high density SNP array

ABSTRACT: BACKGROUND: Btau_4.0 and UMD3.1 are two distinct cattle reference genome assemblies. In our previous study using the low density BovineSNP50 array, we reported a copy number variation (CNV) analysis on Btau_4.0 with 521 animals of 21 cattle breeds, yielding 682 CNV regions with a total length of 139.8 megabases. RESULTS: In this study using the high density BovineHD SNP array, we performed high resolution CNV analyses on both Btau_4.0 and UMD3.1 with 674 animals of 27 cattle breeds. We first compared CNV results derived from these two different SNP array platforms on Btau_4.0. With two thirds of the animals shared between studies, on Btau_4.0 we identified 3,346 candidate CNV regions representing 142.7 megabases (~4.70%) of the genome. With a similar total length but 5 times more event counts, the average CNVR length of current Btau_4.0 dataset is significantly shorter than the previous one (42.7kb vs. 205 kb). Although subsets of these two results overlapped, 64% (91.6 megabases) of current dataset was not present in the previous study. We also performed similar analyses on UMD3.1 using these BovineHD SNP array results. Approximately 50% more and 20% longer CNVs were called on UMD3.1 as compared to those on Btau_4.0. However, a comparable result of CNVRs (3,438 regions with a total length 146.9 megabases) was obtained. We suspect that these results are due to that UMD3.1's efforts of placing unplaced contigs and removing unmerged alleles. Selected CNVs were further experimentally validated, achieving a 73% PCR validation rate, which is considerably higher than the previous validation rate. About 20-45% of CNV regions overlapped with cattle RefSeq genes and Ensembl genes. Panther and IPA analyses indicated that these genes provide a wide spectrum of biological processes involving immune system, lipid metabolism, cell, organism and system development. CONCLUSION: We present a comprehensive result of cattle CNVs at a higher resolution and sensitivity. We identified over 3,000 candidate CNV regions on both Btau_4.0 and UMD3.1, further compared current datasets with previous results, and examined the impacts of genome assemblies on CNV calling.     

A genome‐wide scan for copy number variations using high‐density single nucleotide polymorphism array in Simmental cattle

Copy number variations (CNVs) have recently been identified as promising sources of genetic variation, complementary to single nucleotide polymorphisms (SNPs). As a result, detection of CNVs has attracted a great deal of attention. In this study, we performed genome‐wide CNV detection using Illumina Bovine HD BeadChip (770k) data on 792 Simmental cattle. A total of 263 CNV regions (CNVRs) were identified, which included 137 losses, 102 gains and 24 regions classified as both loss and gain, covering 35.48 Mb (1.41%) of the bovine genome. The length of these CNVRs ranged from 10.18 kb to 1.76 Mb, with an average length of 134.78 kb and a median length of 61.95 kb. In 136 of these regions, a total of 313 genes were identified related to biological functions such as transmembrane activity and olfactory transduction activity. To validate the results, we performed quantitative PCR to detect nine randomly selected CNVRs and successfully confirmed seven (77.6%) of them. Our results present a map of cattle CNVs derived from high‐density SNP data, which expands the current CNV map of the cattle genome and provides useful information for investigation of genomic structural variation in cattle.     

Genome-Wide Identification of Copy Number Variations in Chinese Holstein

Recent studies of mammalian genomes have uncovered the vast extent of copy number variations (CNVs) that contribute to phenotypic diversity. Compared to SNP, a CNV can cover a wider chromosome region, which may potentially incur substantial sequence changes and induce more significant effects on phenotypes. CNV has been becoming an alternative promising genetic marker in the field of genetic analyses. Here we firstly report an account of CNV regions in the cattle genome in Chinese Holstein population. The Illumina Bovine SNP50K Beadchips were used for screening 2047 Holstein individuals. Three different programes (PennCNV, cnvPartition and GADA) were implemented to detect potential CNVs. After a strict CNV calling pipeline, a total of 99 CNV regions were identified in cattle genome. These CNV regions cover 23.24 Mb in total with an average size of 151.69 Kb. 52 out of these CNV regions have frequencies of above 1%. 51 out of these CNV regions completely or partially overlap with 138 cattle genes, which are significantly enriched for specific biological functions, such as signaling pathway, sensory perception response and cellular processes. The results provide valuable information for constructing a more comprehensive CNV map in the cattle genome and offer an important resource for investigation of genome structure and genomic variation underlying traits of interest in cattle.     

A first comparative map of copy number variations in the sheep genome

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.     

Inherited deletion of 9p22.3-p24.3 and duplication of 18p11.31-p11.32 associated with neurodevelopmental delay: Phenotypic matching of involved genes

We describe a 3.5-year-old Iranian female child and her affected 10-month-old brother with a maternally inherited derivative chromosome 9 [der(9)]. The postnatally detected rearrangement was finely characterized by aCGH analysis, which revealed a 15.056 Mb deletion of 9p22.3-p24.3p22.3 encompassing 14 OMIM morbid genes such as DOCK8, KANK1, DMRT1 and SMARCA2, and a gain of 3.309 Mb on 18p11.31-p11.32 encompassing USP14, THOC1, COLEC12, SMCHD1 and LPIN2. We aligned the genes affected by detected CNVs to clinical and functional phenotypic features using PhenogramViz. In this regard, the patient's phenotype and CNVs data were entered into PhenogramViz. For the 9p deletion CNV, 53 affected genes were identified and 17 of them were matched to 24 HPO terms describing the patient's phenotypes. Also, for CNV of 18p duplication, 22 affected genes were identified and six of them were matched to 13 phenotypes. Moreover, we used DECIPHER for in-depth characterization of involved genes in detected CNVs and also comparison of patient phenotypes with 9p and 18p genomic imbalances. Based on our filtration strategy, in the 9p22.3-p24.3 region, approximately 80 pathogenic/likely pathogenic/uncertain overlapping CNVs were in DECIPHER. The size of these CNVs ranged from 12.01 kb to 18.45 Mb and 52 CNVs were smaller than 1 Mb in size affecting 10 OMIM morbid genes. The 18p11.31-p11.32 region overlapped 19 CNVs in the DECIPHER database with the size ranging from 23.42 kb to 1.82 Mb. These CNVs affect eight haploinsufficient genes.     

Detection of copy number variants in the horse genome and examination of their association with recurrent laryngeal neuropathy

We used the data from a recently performed genome‐wide association study using the Illumina Equine SNP50 beadchip for the detection of copy number variants (CNVs) and examined their association with recurrent laryngeal neuropathy (RLN), an important equine upper airway disease compromising performance. A total of 2797 CNVs were detected for 477 horses, covering 229 kb and seven SNPs on average. Overlapping CNVs were merged to define 478 CNV regions (CNVRs). CNVRs, particularly deletions, were shown to be significantly depleted in genes. Fifty‐two of the 67 common CNVRs (frequency ≥ 1%) were validated by association mapping, Mendelian inheritance, and/or Mendelian inconsistencies. None of the 67 common CNVRs were significantly associated with RLN when accounting for multiple testing. However, a duplication on chromosome 10 was detected in 10 cases (representing three breeds) and two unphenotyped parents but in none of the controls. The duplication was embedded in an 8‐Mb haplotype shared across breeds.     

Genome-Wide Detection of Copy Number Variations among Diverse Horse Breeds by Array CGH

Recent studies have found that copy number variations (CNVs) are widespread in human and animal genomes. CNVs are a significant source of genetic variation, and have been shown to be associated with phenotypic diversity. However, the effect of CNVs on genetic variation in horses is not well understood. In the present study, CNVs in 6 different breeds of mare horses, Mongolia horse, Abaga horse, Hequ horse and Kazakh horse (all plateau breeds) and Debao pony and Thoroughbred, were determined using aCGH. In total, seven hundred CNVs were identified ranging in size from 6.1 Kb to 0.57 Mb across all autosomes, with an average size of 43.08 Kb and a median size of 15.11 Kb. By merging overlapping CNVs, we found a total of three hundred and fifty-three CNV regions (CNVRs). The length of the CNVRs ranged from 6.1 Kb to 1.45 Mb with average and median sizes of 38.49 Kb and 13.1 Kb. Collectively, 13.59 Mb of copy number variation was identified among the horses investigated and accounted for approximately 0.61% of the horse genome sequence. Five hundred and eighteen annotated genes were affected by CNVs, which corresponded to about 2.26% of all horse genes. Through the gene ontology (GO), genetic pathway analysis and comparison of CNV genes among different breeds, we found evidence that CNVs involving 7 genes may be related to the adaptation to severe environment of these plateau horses. This study is the first report of copy number variations in Chinese horses, which indicates that CNVs are ubiquitous in the horse genome and influence many biological processes of the horse. These results will be helpful not only in mapping the horse whole-genome CNVs, but also to further research for the adaption to the high altitude severe environment for plateau horses.     

Gene-based copy number variation study reveals a microdeletion at 12q24 that influences height in the Korean population

Height is a classic polygenic trait with high heritability (h² = 0.8). Recent genome-wide association studies have revealed many independent loci associated with human height. In addition, although many studies have reported an association between copy number variation (CNV) and complex diseases, few have explored the relationship between CNV and height. Recent studies reported that single nucleotide polymorphisms (SNPs) are highly correlated with common CNVs, suggesting that it is warranted to survey CNVs to identify additional genetic factors affecting heritable traits such as height.This study tested the hypothesis that there would be CNV regions (CNVRs) associated with height nearby genes from the GWASs known to affect height. We identified regions containing > 1% copy number deletion frequency from 3667 population-based cohort samples using the Illumina HumanOmni1-Quad BeadChip. Among the identified CNVRs, we selected 15 candidate regions that were located within 1 Mb of 283 previously reported genes. To assess the effect of these CNVRs on height, statistical analyses were conducted with samples from a case group of 370 taller (upper 10%) individuals and a control group of 1828 individuals (lower 50%).We found that a newly identified 17.7 kb deletion at chromosomal position 12q24.33, approximately 171.6 kb downstream of GPR133, significantly correlated with height; this finding was validated using quantitative PCR. These results suggest that CNVs are potentially important in determining height and may contribute to height variation in human populations.     

Copy number variations identified in the chicken using a 60K SNP BeadChip

Copy number variation (CNV) is considered an important genetic variation, contributing to many economically important traits in the chicken. Although CNVs can be detected using a comparative genomic hybridization array, the high‐density SNP array has provided an alternative way to identify CNVs in the chicken. In the current study, a chicken 60K SNP BeadChip was used to identify CNVs in two distinct chicken genetic lines (White Leghorn and dwarf) using the penncnv program. A total of 209 CNV regions were identified, distributing on chromosomes 1–22 and 24–28 and encompassing 13.55 Mb (1.42%) of chicken autosomal genome area. Three of seven selected CNVs (73.2% individuals) were completely validated by quantitative PCR. To our knowledge, this is the first report in the chicken identifying CNVs using a SNP array. Identification of 190 new identified CNVs illustrates the feasibility of the chicken 60K SNP BeadChip to detect CNVs in the chicken, which lays a solid foundation for future analyses of associations of CNVs with economically important phenotypes in chickens.     

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Background

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F_ST. The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

Results

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F_ST, potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F_ST signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F_ST values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

Conclusions

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-876) contains supplementary material, which is available to authorized users.     

The application of gene co-expression network reconstruction based on CNVs and gene expression microarray data in breast cancer

Copy number variations (CNVs) are one type of the human genetic variations and are pervasive in the human genome. It has been confirmed that they can play a causal role in complex diseases. Previous studies of CNVs focused more on identifying the disease-specific CNV regions or candidate genes on these CNV regions, but less on the synergistic actions between genes on CNV regions and other genes. Our research combined the CNVs with related gene co-expression to reconstruct gene co-expression network by using single nucleotide polymorphism microarray datasets and gene microarray datasets of breast cancer, and then extracted the modules which connected densely inside and analyzed the functions of modules. Interestingly, all of these modules’ functions were related to breast cancer according to our enrichment analysis, and most of the genes in these modules have been reported to be involved in breast cancer. Our findings suggested that integrating CNVs and gene co-expressed relations was an available way to analyze the roles of CNV genes and their synergistic genes in breast cancer, and provided a novel insight into the pathological mechanism of breast cancer.     

Tropically adapted cattle of Africa: perspectives on potential role of copy number variations

Africa is host to diverse and locally adapted cattle breeds that are expected to survive the harsh and extreme tropical environments associated with diseases and parasite infections, heat stress and episodes of feed and water scarcity. Genomic copy number variations (CNVs) are considered to be primary role players in cattle breed formation and adaptation where isolation and genetic drift together with subsequent mutations have created an enormous diversity of local populations. CNVs are modifications in DNA structure comprising deletions, duplications and insertions that are >1 kb in size. Despite attracting much attention, the frequency and pattern of bovine CNV events, especially in African cattle breeds, are for the most part largely unknown. Characterization of genetic variation in the indigenous cattle of Africa will be a vital step toward dissecting the molecular mechanisms underlying phenotypic variation and local adaptation. This review therefore aims to describe the current knowledge regarding bovine CNVs and the implications and potentials they encompass for dissecting genetic adaptation and the genotypic skeleton of tropical African cattle populations.     

Disruption of Contactin 4 in two subjects with autism in Chinese population

Autism is a heterogeneous childhood neurodevelopmental disorder that is characterised by deficits in verbal communication, impaired social interactions, restricted interests and repetitive behaviours. Using an Illumina HumanCNV370-Quad BeadChip, we identified two Han Chinese individuals with autism and large duplications (~1.6 Mb and ~2.4 Mb) disrupting the same CNTN4 gene. CNTN4 encodes a protein that functions as a cell-adhesion molecule and may play an essential role in the formation of axon connections in the developing nervous system. The disruption of this gene has been reported to be the cause of the 3p deletion syndrome and also a possible susceptibility factor for autism spectrum disorders (ASDs). Our results suggest that rare copy number variations (CNVs) in CNTN4 may also influence autism susceptibility in Asian populations. Interestingly, a comparison of the clinical phenotypes between the two subjects revealed that the subject with the 2.4 Mb CNV (involving several other genes) presented with a more severe phenotype than the subject with the 1.6 Mb CNV (disrupting only CNTN4 and CNTN6). This suggests that other genes in the nearby region may contribute to the pathogenesis.     

猪13号染色体上拷贝数变异区内基因信息发掘及遗传规律分析

拷贝数变异(Copy number variation, CNV)是染色体上发生的一种微结构变异, 已引起越来越多研究者的关注。本课题组前期已获得猪13号染色体上的32个CNV区域(CNV region, CNVR), 为了发掘CNVR内的基因信息, 文章在线检索了上述CNVR内的基因并进行基因本体(Gene Ontology)分析。结果共发现236个基因, 其中有注释基因169个, 主要参与蛋白质水解、细胞粘附、大分子降解等生物过程。为了探索这些基因拷贝数变异的遗传规律, 文章选择RCAN1(Regulators of calcineurin 1)基因为候选基因, 利用QPCR方法在莱芜猪群中检测了该基因的拷贝数, 并分析了CNV在莱芜猪3个家系中的遗传规律。结果表明, RCAN1基因在莱芜猪群体中存在拷贝数的缺失、重复现象, 其拷贝数变异的遗传规律符合孟德尔遗传方式。     

Copy number variation in the cattle genome

Copy number variations (CNVs) are gains and losses of genomic sequence greater than 50?bp between two individuals of a species. While single nucleotide polymorphisms (SNPs) are more frequent, CNVs impact a higher percentage of genomic sequence and have potentially greater effects, including the changing of gene structure and dosage, altering gene regulation and exposing recessive alleles. In particular, segmental duplications (SDs) were shown to be one of the catalysts and hotspots for CNV formation. Substantial progress has been made in understanding CNVs in mammals, especially in humans and rodents. CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility. Recently, interest in CNV study has extended into domesticated animals, including cattle. Multiple genome-wide cattle CNV studies have been carried out using both microarray and next generation sequencing technologies. Integration of SD and CNV results with SNP and other datasets are beginning to reveal impacts of CNVs on cattle domestication, health, and production traits.     

A genome-wide survey reveals a deletion polymorphism associated with resistance to gastrointestinal nematodes in Angus cattle

Gastrointestinal (GI) nematode infections are a worldwide threat to human health and animal production. In this study, we performed a genome-wide association study between copy number variations (CNVs) and resistance to GI nematodes in an Angus cattle population. Using a linear regression analysis, we identified one deletion CNV which reaches genome-wide significance after Bonferroni correction. With multiple mapped human olfactory receptor genes but no annotated bovine genes in the region, this significantly associated CNV displays high population frequencies (58.26 %) with a length of 104.8 kb on chr7. We further investigated the linkage disequilibrium (LD) relationships between this CNV and its nearby single nucleotide polymorphisms (SNPs) and genes. The underlining haplotype blocks contain immune-related genes such as ZNF496 and NLRP3. As this CNV co-segregates with linked SNPs and associated genes, we suspect that it could contribute to the detected variations in gene expression and thus differences in host parasite resistance.     

Copy number polymorphism in plant genomes

Copy number variants (CNVs) are genomic rearrangements resulting from gains or losses of DNA segments. Typically, the term refers to rearrangements of sequences larger than 1 kb. This type of polymorphism has recently been shown to be a key contributor to intra-species genetic variation, along with single-nucleotide polymorphisms and short insertion-deletion polymorphisms. Over the last decade, a growing number of studies have highlighted the importance of copy number variation (CNV) as a factor affecting human phenotype and individual CNVs have been linked to risks for severe diseases. In plants, the exploration of the extent and role of CNV is still just beginning. Initial genomic analyses indicate that CNVs are prevalent in plants and have greatly affected plant genome evolution. Many CNV events have been observed in outcrossing and autogamous species. CNVs are usually found on all chromosomes, with CNV hotspots interspersed with regions of very low genetic variation. Although CNV is mainly associated with intergenic regions, many CNVs encompass protein-coding genes. The collected data suggest that CNV mainly affects the members of large families of functionally redundant genes. Thus, the effects of individual CNV events on phenotype are usually modest. Nevertheless, there are many cases in which CNVs for specific genes have been linked to important traits such as flowering time, plant height and resistance to biotic and abiotic stress. Recent reports suggest that CNVs may form rapidly in response to stress.     

Genome‐wide copy number profiling using high‐density SNP array in chickens

Phenotypic diversity is a direct consequence resulting mainly from the impact of underlying genetic variation, and recent studies have shown that copy number variation (CNV) is emerging as an important contributor to both phenotypic variability and disease susceptibility. Herein, we performed a genome‐wide CNV scan in 96 chickens from 12 diversified breeds, benefiting from the high‐density Affymetrix 600 K SNP arrays. We identified a total of 231 autosomal CNV regions (CNVRs) encompassing 5.41 Mb of the chicken genome and corresponding to 0.59% of the autosomal sequence. The length of these CNVRs ranged from 2.6 to 586.2 kb with an average of 23.4 kb, including 130 gain, 93 loss and eight both gain and loss events. These CNVRs, especially deletions, had lower GC content and were located particularly in gene deserts. In particular, 102 CNVRs harbored 128 chicken genes, most of which were enriched in immune responses. We obtained 221 autosomal CNVRs after converting probe coordinates to Galgal3, and comparative analysis with previous studies illustrated that 153 of these CNVRs were regarded as novel events. Furthermore, qPCR assays were designed for 11 novel CNVRs, and eight (72.73%) were validated successfully. In this study, we demonstrated that the high‐density 600 K SNP array can capture CNVs with higher efficiency and accuracy and highlighted the necessity of integrating multiple technologies and algorithms. Our findings provide a pioneering exploration of chicken CNVs based on a high‐density SNP array, which contributes to a more comprehensive understanding of genetic variation in the chicken genome and is beneficial to unearthing potential CNVs underlying important traits of chickens.