Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. We previously reported an initial analysis of copy number variations (CNVs) in Angus cattle selected for resistance or susceptibility to gastrointestinal nematodes. In this study, we performed a large-scale analysis of CNVs using SNP genotyping data from 472 animals of the same population. We detected 811 candidate CNV regions, which represent 141.8 Mb (~4.7%) of the genome. To investigate the functional impacts of CNVs, we created 2 groups of 100 individual animals with extremely low or high estimated breeding values of eggs per gram of feces and referred to these groups as parasite resistant (PR) or parasite susceptible (PS), respectively. We identified 297 (~51 Mb) and 282 (~48 Mb) CNV regions from PR and PS groups, respectively. Approximately 60% of the CNV regions were specific to the PS group or PR group of animals. Selected PR- or PS-specific CNVs were further experimentally validated by quantitative PCR. A total of 297 PR CNV regions overlapped with 437 Ensembl genes enriched in immunity and defense, like WC1 gene which uniquely expresses on gamma/delta T cells in cattle. Network analyses indicated that the PR-specific genes were predominantly involved in gastrointestinal disease, immunological disease, inflammatory response, cell-to-cell signaling and interaction, lymphoid tissue development, and cell death. By contrast, the 282 PS CNV regions contained 473 Ensembl genes which are overrepresented in environmental interactions. Network analyses indicated that the PS-specific genes were particularly enriched for inflammatory response, immune cell trafficking, metabolic disease, cell cycle, and cellular organization and movement.     

Signatures of selection in sheep bred for resistance or susceptibility to gastrointestinal nematodes

Background

Gastrointestinal nematodes are one of the most serious causes of disease in domestic ruminants worldwide. There is considerable variation in resistance to gastrointestinal nematodes within and between sheep breeds, which appears to be due to underlying genetic diversity. Through selection of resistant animals, rapid genetic progress has been demonstrated in both research and commercial flocks. Recent advances in genome sequencing and genomic technologies provide new opportunities to understand the ovine host response to gastrointestinal nematodes at the molecular level, and to identify polymorphisms conferring nematode resistance.

Results

Divergent lines of Romney and Perendale sheep, selectively bred for high and low faecal nematode egg count, were genotyped using the Illumina® Ovine SNP50 BeadChip. The resulting genome-wide SNP data were analysed for selective sweeps on loci associated with resistance or susceptibility to gastrointestinal nematode infection. Population differentiation using F_ST and Peddrift revealed sixteen regions, which included candidate genes involved in chitinase activity and the cytokine response. Two of the sixteen regions identified were contained within previously identified QTLs associated with nematode resistance.

Conclusions

In this study we identified fourteen novel regions associated with resistance or susceptibility to gastrointestinal nematodes. Results from this study support the hypothesis that host resistance to internal nematode parasites is likely to be controlled by a number of loci of moderate to small effects.

Electronic supplementary material

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

拷贝数变异的全基因组关联分析

基因组拷贝数变异(copy number variations,CNVs)是指与基因组参考序列相比,基因组中≥1 kb的DNA片段插入、缺失和/或扩增,及其互相组合衍生出的复杂变异．由于其具有分布范围广、可遗传、相对稳定和高度异质性等特点,目前认为,CNVs是一种新的可以作为疾病易感标志的基因组DNA多态性,其变异引起的基因剂量改变可以导致表型改变．最近,一种基于CNVs的新的疾病易感基因鉴定策略——CNV全基因组关联分析开始出现,这一策略和传统的基于单核苷酸多态性的关联分析具有互补性,通过认识基因组结构变异可以认识复杂疾病的分子机制和遗传基础．     

New copy number variations in schizophrenia

Genome-wide screenings for copy number variations (CNVs) in patients with schizophrenia have demonstrated the presence of several CNVs that increase the risk of developing the disease and a growing number of large rare CNVs; the contribution of these rare CNVs to schizophrenia remains unknown. Using Affymetrix 6.0 arrays, we undertook a systematic search for CNVs in 172 patients with schizophrenia and 160 healthy controls, all of Italian origin, with the aim of confirming previously identified loci and identifying novel schizophrenia susceptibility genes. We found five patients with a CNV occurring in one of the regions most convincingly implicated as risk factors for schizophrenia: NRXN1 and the 16p13.1 regions were found to be deleted in single patients and 15q11.2 in 2 patients, whereas the 15q13.3 region was duplicated in one patient. Furthermore, we found three distinct patients with CNVs in 2q12.2, 3q29 and 17p12 loci, respectively. These loci were previously reported to be deleted or duplicated in patients with schizophrenia but were never formally associated with the disease. We found 5 large CNVs (>900 kb) in 4q32, 5q14.3, 8q23.3, 11q25 and 17q12 in five different patients that could include some new candidate schizophrenia susceptibility genes. In conclusion, the identification of previously reported CNVs and of new, rare, large CNVs further supports a model of schizophrenia that includes the effect of multiple, rare, highly penetrant variants.     

Breeding cattle and sheep for resistance to gastrointestinal nematodes

Gastrointestinal nematodes are an important cause of reduced production of meat, milk and wool in domestic livestock. It is generally believed that problems caused by these parasites have increased owing to the intensification of animal husbandry(1-3) of resistance to anthelmintics, current research is focussed on alternative control strategies that do not rely on anthelmintics. Here, Bram Kloosterman, Henk Parmentier and Harm Ploeger review work on the genetic resistance of domestic ruminants to these nematodes and discuss the practicality of breeding programmes.     

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.     

Modeling genetic inheritance of copy number variations

Copy number variations (CNVs) are being used as genetic markers or functional candidates in gene-mapping studies. However, unlike single nucleotide polymorphism or microsatellite genotyping techniques, most CNV detection methods are limited to detecting total copy numbers, rather than copy number in each of the two homologous chromosomes. To address this issue, we developed a statistical framework for intensity-based CNV detection platforms using family data. Our algorithm identifies CNVs for a family simultaneously, thus avoiding the generation of calls with Mendelian inconsistency while maintaining the ability to detect de novo CNVs. Applications to simulated data and real data indicate that our method significantly improves both call rates and accuracy of boundary inference, compared to existing approaches. We further illustrate the use of Mendelian inheritance to infer SNP allele compositions in each of the two homologous chromosomes in CNV regions using real data. Finally, we applied our method to a set of families genotyped using both the Illumina HumanHap550 and Affymetrix genome-wide 5.0 arrays to demonstrate its performance on both inherited and de novo CNVs. In conclusion, our method produces accurate CNV calls, gives probabilistic estimates of CNV transmission and builds a solid foundation for the development of linkage and association tests utilizing CNVs.     

Analysis of copy number variants in the cattle genome

Copy number variation (CNV) is likely to be an important component of heritable variation in livestock. To characterise CNVs in cattle, we performed a genome wide survey to determine the number, location and gene content of these genomic features. A tiling oligonucleotide array with ~385,000 probes was used for comparative genomic hybridisation of both taurine and zebu cattle. Using a conservative set of calling criteria, a total of 51 CNV were detected that collectively spanned approximately half of one percent of the bovine genome. The size of the average CNV within each animal ranged from 213 kb up to 335 kb. Half of the CNV were detected in a single animal only, whilst the remainder was independently identified in multiple individuals. Analysis was performed to determine the gene content for each CNV region. This revealed that the majority of CNV (82%) spanned at least one gene, with a number of CNV containing genes which are known to control aspects of phenotypic variation in cattle. Whilst additional studies are required to determine the impact of individual CNV, this study confirmed them as an important class of genomic variation in cattle.     

Genome-wide patterns of copy number variations in Spodoptera litura

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.     

A map of copy number variations in Chinese populations

It has been shown that the human genome contains extensive copy number variations (CNVs). Investigating the medical and evolutionary impacts of CNVs requires the knowledge of locations, sizes and frequency distribution of them within and between populations. However, CNV study of Chinese minorities, which harbor the majority of genetic diversity of Chinese populations, has been underrepresented considering the same efforts in other populations. Here we constructed, to our knowledge, a first CNV map in seven Chinese populations representing the major linguistic groups in China with 1,440 CNV regions identified using Affymetrix SNP 6.0 Array. Considerable differences in distributions of CNV regions between populations and substantial population structures were observed. We showed that ～35% of CNV regions identified in minority ethnic groups are not shared by Han Chinese population, indicating that the contribution of the minorities to genetic architecture of Chinese population could not be ignored. We further identified highly differentiated CNV regions between populations. For example, a common deletion in Dong and Zhuang (44.4% and 50%), which overlaps two keratin-associated protein genes contributing to the structure of hair fibers, was not observed in Han Chinese. Interestingly, the most differentiated CNV deletion between HapMap CEU and YRI containing CCL3L1 gene reported in previous studies was also the highest differentiated regions between Tibetan and other populations. Besides, by jointly analyzing CNVs and SNPs, we found a CNV region containing gene CTDSPL were in almost perfect linkage disequilibrium between flanking SNPs in Tibetan while not in other populations except HapMap CHD. Furthermore, we found the SNP taggability of CNVs in Chinese populations was much lower than that in European populations. Our results suggest the necessity of a full characterization of CNVs in Chinese populations, and the CNV map we constructed serves as a useful resource in further evolutionary and medical studies.     

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.     

Common copy number variations in fifty radiosensitive cell lines

Hypersensitivity to radiation exposure is a major challenge to radiotherapy in the treatment of cancer patients. Copy number variations (CNVs) are believed to identify genomic regions of functional significance for radiosensitivity (RS) but have yet to be systematically investigated. We used Affymetrix 6.0 SNP arrays to survey common CNVs in a cohort of 50 radiosensitive lymphoblastoid cell lines (RS-LCLs) derived from patients with undiagnosed diseases. A total of 317 CNVs that were present in at least 10% of the studied cell lines were identified. Three hundred and eight CNVs overlapped with polymorphic CNVs, 13 of which were significantly enriched in the RS-LCLs compared to the reference. The remaining 9 CNVs were novel. The majority of these enriched and novel CNVs were chromosomal gains. The dominance of the chromosomal gains over losses is inconsistent with the traditional concept of molecular basis of RS and suggests more complex genetic mechanisms for RS.     

Mx mRNA expression and RFLP analysis of rainbow trout Oncorhynchus mykiss genetic crosses selected for susceptibility or resistance to IHNV

Three interferon-inducible Mx genes have been identified in rainbow trout Oncorhynchus mykiss and their roles in virus resistance have yet to be determined. In mice, expression of the Mx1 protein is associated with resistance to influenza virus. We report a study to determine whether there was a correlation between the expression of Mx in rainbow trout and resistance to a fish rhabdovirus, infectious hematopoietic necrosis virus (IHNV). A comparison of Mx mRNA expression was made between different families of cultured rainbow trout selected for resistance or for susceptibility to IHNV. A trout-specific Mx cDNA gene probe was used to determine whether there was a correlation between Mx mRNA expression and resistance to the lethal effects of IHNV infection. Approximately 99% of trout injected with a highly virulent strain of the fish rhabdovirus, IHNV, were able to express full length Mx mRNA at 48 h post infection. This is markedly different from the expression of truncated, non-functional Mx mRNA found in most laboratory strains of mice, and the ability of only 25% of wild mice to express functional Mx protein. A restriction fragment length polymorphism (RFLP) assay was developed to compare the Mx locus between individual fish and between rainbow trout genetic crosses bred for IHNV resistance or susceptibility. The assay was able to discriminate 7 distinct RFLP patterns in the rainbow trout crosses. One cross was identified that showed a correlation between homozygosity at the Mx locus and greater susceptibility to IHN-caused mortality.     

BoLA-DRB3 gene polymorphism and FMD resistance or susceptibility in Wanbei cattle

For the further characterization of foot-and-mouth disease virus (FMDV)-induced foot-and-mouth disease, we investigated the association between polymorphism of BoLA-DRB3 gene and FMD resistance/susceptibility of Wanbei cattle challenged with FMDV. One hundred cattle were challenged with FMDV and exon 2 of BoLA-DRB3 genes was amplified by hemi-nested polymerase chain reaction from asymptomatic animals and from animals with FMD. PCR products were characterized by the RFLP technique using restriction enzymes Hae III. The results revealed extensive polymorphisms, 6 RFLP patterns were identified. By analyzing alleles and genotypic frequencies between healthy and infection with FMD cattle, we found that allele Hae III A was associated with susceptibility to FMD in Wanbei cattle (P < 0.05), whereas Hae III C was associated with resistance to FMD (P < 0.01) and may have a strong protective effect against FMD. Hae IIICC and Hae III BC genotype were associated with resistance to FMD (P < 0.01). By contrast, Hae III AA genotype was associated with susceptibility to FMD (P < 0.01). Sequence analysis show that 89 amino acids were translated in exon 2 of BoLA-DRB3 and 13.70 % of nucleotide mutated, which resulted in 14.61 % of amino acid change. One PKC, one Tyr and one CAMP phosphorylation were increased; the hydrophobicity and secondary structure of proteins produced change after amino acid substitution. These results revealed that Wanbei cattle had the ability of resistance to disease by mutation which result changes of the protein structure to perform the regulation of the cell using different signaling pathways in the long process of choice evolution.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

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