Initial analysis of copy number variations in cattle selected for resistance or susceptibility to intestinal nematodes

Identification of the genetic characteristics of copy number variations in experimental specific pathogen-free ducks using whole-genome resequencing

BACKGROUND

Specific pathogen-free ducks are a valuable laboratory resource for waterfowl disease research and poultry vaccine development. High throughput sequencing allows the systematic identification of structural variants in genomes. Copy number variation (CNV) can explain the variation of important duck genetic traits. Herein, the genome-wide CNVs of the three experimental duck species in China (Jinding ducks (JD), Shaoxing ducks (SX), and Fujian Shanma ducks (SM)) were characterized using resequencing to determine their genetic characteristics and selection signatures.

RESULTS

We obtained 4,810 CNV regions (CNVRs) by merging 73,012 CNVs, covering 4.2% of the duck genome. Functional analysis revealed that the shared CNVR-harbored genes were significantly enriched for 31 gene ontology terms and 16 Kyoto Encyclopedia of Genes and Genomes pathways (e.g., olfactory transduction and immune system). Based on the genome-wide fixation index for each CNVR, growth (SPAG17 and PTH1R), disease resistance (CATHL3 and DMBT1), and thermoregulation (TRPC4 and SLIT3) candidate genes were identified in strongly selected signatures specific to JD, SM, and SX, respectively.

CONCLUSIONS

In conclusion, we investigated the genome-wide distribution of experimental duck CNVs, providing a reference to establish the genetic basis of different phenotypic traits, thus contributing to the management of experimental animal genetic resources.

Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent

Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.

Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection

BACKGROUND

Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding.

RESULTS

CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences.

CONCLUSIONS

Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.

Copy Number Variation Discovery in South African Nguni-Sired and Bonsmara-Sired Crossbred Cattle

Crossbreeding forms part of Climate-Smart beef production and is one of the strategies to mitigate the effects of climate change. Two Nguni-sired and three Bonsmara-sired crossbred animals underwent whole genome sequencing. Following quality control and file preparation, the sequence data were investigated for genome-wide copy number variation (CNV) using the panelcn.MOPS tool. A total of 355 CNVs were identified in the crossbreds, of which 274 were unique in Bonsmara-sired crossbreds and 81 unique in the Nguni-sired crossbreds. Genes that differed in copy number in both crossbreds included genes related to growth (SCRN2, LOC109572916) and fertility-related factors (RPS28, LOC1098562432, LOC109570037). Genes that were present only in the Bonsmara-sired crossbreds included genes relating to lipid metabolism (MAF1), olfaction (LOC109569114), body size (HES7), immunity (LOC10957335, LOC109877039) and disease (DMBT1). Genes that were present only in the Nguni-sired crossbreds included genes relating to ketosis (HMBOX1) and amino acid transport (LOC109572916). Results of this study indicate that Nguni and Bonsmara cattle can be utilized in crossbreeding programs as they may enhance the presence of economically important traits associated with both breeds. This will produce crossbred animals that are good meat producers, grow faster, have high fertility, strong immunity and a better chance of producing in South Africa's harsh climate conditions. Ultimately, this study provides new genetic insights into the adaptability of Nguni and Bonsmara crossbred cattle.

Optical mapping compendium of structural variants across global cattle breeds

Structural variants (SV) have been linked to important bovine disease phenotypes, but due to the difficulty of their accurate detection with standard sequencing approaches, their role in shaping important traits across cattle breeds is largely unexplored. Optical mapping is an alternative approach for mapping SVs that has been shown to have higher sensitivity than DNA sequencing approaches. The aim of this project was to use optical mapping to develop a high-quality database of structural variation across cattle breeds from different geographical regions, to enable further study of SVs in cattle. To do this we generated 100X Bionano optical mapping data for 18 cattle of nine different ancestries, three continents and both cattle sub-species. In total we identified 13,457 SVs, of which 1,200 putatively overlap coding regions. This resource provides a high-quality set of optical mapping-based SV calls that can be used across studies, from validating DNA sequencing-based SV calls to prioritising candidate functional variants in genetic association studies and expanding our understanding of the role of SVs in cattle evolution.

Measurement(s)	Optical Mapping
Technology Type(s)	Optical Mapping
Factor Type(s)	Structural variants
Sample Characteristic - Organism	Bos taurus
Sample Characteristic - Location	United Kingdom • Kenya • Zambia • Uganda • Brazil • Nigeria

Genetic Architecture and Signatures of Selection in the Caqueteño Creole (Colombian Native Cattle)

Evolutionary mechanisms have shaped the genomic architecture of Colombian Creole cattle breeds. The mating and selection processes have impacted several traits, promoting differences within and between populations. Studies of population structure and selection signatures in Colombian Creole breeds are scarce, and need more attention to better understand genetic differentiation, gene flow, and genetic distance. This study aimed to analyze the population structure and identify selection imprints in the Criollo Caqueteño (CAQ) population. It used 127 CAQ animals genotyped with Chip HD 777,000 SNPs. The population structure analyses used discriminant principal component analysis (DAPC), integrated haplotype scoring (iHS), and index-fixing (Fst) methodologies to detect selection signals. We can highlight SNP regions on the genes TMPRSS15, PGAM2, and EGFR, identified by the Fst method. Additionally, the iHS regions for cluster 1 identified candidate genes on BTA 3 (CMPK1 and FOXD2), BTA 11 (RCAN1), and BTA 22 (ARPP21). In group 2, we can highlight the genes on BTA 4 (SLC13A4, BRAF), BTA 9 (ULBP), BTA 14 (CSMD3) and BTA 19 (KRTAP9-2). These candidate genes have been associated with fertility traits, precocity, growth, and environmental and disease resistance, indicating a genetic potential in CAQ animals. All this promotes a better understanding of the diversity and genetic structure in the CAQ population. Based on that, our study can significantly assist the sustainable development and conservation of the breed in the Colombian Amazon.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing

Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

Distribution of Copy Number Variation in SYT11 Gene and Its Association with Growth Conformation Traits in Chinese Cattle

Simple Summary

It is known that many different breeds of cattle are widely distributed in China. However, due to a lengthy selection of draught direction, there are obvious shortcomings in Chinese cattle, such as less meat production, slow weight gain, poor meat quality, and a lack of specialized beef cattle breeds. Animal breeding heavily benefits from molecular technologies, among which molecular genetic markers were widely used to improve the economic traits of beef cattle. Because the copy number variation (CNV) involves a longer DNA sequence or even the entire functional gene, it may have a greater impact on the phenotype. Recent studies have indicated that CNVs are widespread in the Chinese cattle genome. By investigating the effects of CNVs on gene expression and cattle traits, we aim to find those genomic variations which could significantly affect cattle traits, and which could provide a basis for genetic selection and molecular breeding of local Chinese cattle.

Abstract

Currently, studies of the SYT11 gene mainly focus on neurological diseases such as schizophrenia and Parkinson’s disease. However, some studies have shown that the C2B domain of SYT11 can interact with RISC components and affect the gene regulation of miRNA, which is important for cell differentiation, proliferation, and apoptosis, and therefore has an impact on muscle growth and development in animals. The whole-genome resequencing data detected a CNV in the SYT11 gene, and this may affect cattle growth traits. In this study, CNV distribution of 672 individuals from four cattle breeds, Yunling, Pinan, Xianan, and Qinchuan, were detected by qPCR. The relationship between CNV, gene expression and growth traits was further investigated. The results showed that the proportion of multiple copy types was the largest in all cattle breeds, but there were some differences among different breeds. The normal type had higher gene expression than the abnormal copy type. The CNVs of the SYT11 gene were significantly correlated with body length, cannon circumference, chest depth, rump length, and forehead size of Yunling cattle, and was significantly correlated with the bodyweight of Xianan cattle, respectively. These data improve our understanding of the effects of CNV on cattle growth traits. Our results suggest that the CNV of SYT11 gene is a protentional molecular marker, which may be used to improve growth traits in Chinese cattle.

Genome-wide detection of copy number variations in Tharparkar cattle

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

A genome-wide scan of copy number variants in three Iranian indigenous river buffaloes

BACKGROUND

In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K.

RESULTS

CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database.

CONCLUSIONS

This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo's natural adaptive history together to a recent selection for milk used as primary food source from this species.

CNV analysis of Meishan pig by next-generation sequencing and effects of AHR gene CNV on pig reproductive traits

Background

Reproductive performance of livestock is an economically important aspect of global food production. The Chinese Meishan pig is a prolific breed, with an average of three to five more piglets per litter than European breeds; however, the genetic basis for this difference is not well understood.

Results

In this study, we investigated copy number variations (CNVs) of 32 Meishan pigs and 29 Duroc pigs by next-generation sequencing. A genome-wide analysis of 61 pigs revealed 12,668 copy number variable regions (CNVRs) that were further divided into three categories based on copy number (CN) of the whole population, i.e., gain (n = 7,638), and loss (n = 5,030) CNVRs. We then compared Meishan and Duroc pigs and identified 17.17 Mb of 6,387 CNVRs that only existing in Meishan pigs CNVRs that overlapped the reproduction-related gene encoding the aryl hydrocarbon receptor (AHR) gene. We found that normal AHR CN was more frequent than CN loss in four different pig breeds. An association analysis showed that AHR CN had a positive effect on litter size (P < 0.05) and that a higher CN was associated with higher total number born (P < 0.05), number born alive (P < 0.05), number of weaned piglets, and birth weight.

Conclusions

The present study provides comprehensive CNVRs for Meishan and Duroc pigs through large-scale population resequencing. Our results provide a supplement for the high-resolution map of copy number variation in the porcine genome and valuable information for the investigation of genomic structural variation underlying traits of interest in pig. In addition, the association results provide evidence for AHR as a candidate gene associated with reproductive traits that can be used as a genetic marker in pig breeding programs.

BGVD: An Integrated Database for Bovine Sequencing Variations and Selective Signatures

Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.

Four single nucleotide polymorphisms (SNPs) are associated with resistance and resilience to Haemonchus contortus in Brazilian Morada Nova sheep

Gastrointestinal nematodes are a major constraint in sheep production. Breeding for resistance has proven to be an effective and feasible approach to address this problem. The use and investigation of genetic markers for resistance traits could accelerate genetic progress and lead to a better understanding of underlying molecular mechanisms. Thus, the aim of this study was to evaluate if five single nucleotide polymorphisms SNPs OAR2_14765360, OAR6_81718546, OAR11_62887032, OAR12_69606944 and OAR15_59871543 are associated with resistance and resilience traits in a flock of the Morada Nova sheep breed. Lambs were submitted to two consecutive parasite challenges by oral infection with 4000 infective larvae L₃) of Haemonchus contortus. Fecal egg counts (FEC), packed cell volume (PVC) and body weight were measured every one or two weeks for 42 days in each trial. DNA samples from 287 lambs, 131 ewes and 4 rams were amplified by ARMS-PCR or PCR-RFLP and genotypes were determined. Analysis of variance (ANOVA) was used for association analyses between genotypes and phenotypes. In case of significant association, the allele substitution effect was calculated based on a linear model. OAR2_14765360 and OAR12_69606944 were associated with FEC, and OAR12_69606944 also had significant effects on PCV and weight gain, showing favourable associations of the CC genotype with all evaluated traits. Both OAR6_81718546 and OAR11_62887032 were associated with weight gain, and OAR6_81718546 had an additional effect on PCV. OAR15_59871543 was not polymorphic in the population. OAR6_81718546 and OAR12_69606944 presented significant allele substitution effects of -1.06 ± 0.52 kg for the T allele on final body weight and 0.74 ± 0.32 for the C allele in PCV of the same sampling date, respectively. This is the first report of SNPs associated with gastrointestinal nematode resistance in this sheep breed. Our findings support the existence of quantitatice trait loci (QTL) for resistance and resilience in linkage disequilibrium with the polymorphic SNPs and suggest their future use for explorations of these traits in Morada Nova sheep.

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

Genome-wide associations and functional gene analyses for endoparasite resistance in an endangered population of native German Black Pied cattle

Background

Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections.

Results

A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (p_Bonf = 4.47 × 10^− 7), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10^− 6 on BTA 1 to 2.18 × 10^− 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections.

Conclusions

We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5659-4) contains supplementary material, which is available to authorized users.

Copy number variation study in Japanese quail associated with stress related traits using whole genome re-sequencing data

Copy number variation (CNV) is a major driving factor for genetic variation and phenotypic diversity in animals. To detect CNVs and understand genetic components underlying stress related traits, we performed whole genome re-sequencing of pooled DNA samples of 20 birds each from High Stress (HS) and Low Stress (LS) Japanese quail lines using Illumina HiSeq 2×150 bp paired end method. Sequencing data were aligned to the quail genome and CNVnator was used to detect CNVs in the aligned data sets. The depth of coverage for the data reached to 41.4x and 42.6x for HS and LS birds, respectively. We identified 262 and 168 CNV regions affecting 1.6 and 1.9% of the reference genome that completely overlapped 454 and 493 unique genes in HS and LS birds, respectively. Ingenuity pathway analysis showed that the CNV genes were significantly enriched to phospholipase C signaling, neuregulin signaling, reelin signaling in neurons, endocrine and nervous development, humoral immune response, and carbohydrate and amino acid metabolisms in HS birds, whereas CNV genes in LS birds were enriched in cell-mediated immune response, and protein and lipid metabolisms. These findings suggest CNV genes identified in HS and LS birds could be candidate markers responsible for stress responses in birds.

Genomic predictions combining SNP markers and copy number variations in Nellore cattle

Background

Due to the advancement in high throughput technology, single nucleotide polymorphism (SNP) is routinely being incorporated along with phenotypic information into genetic evaluation. However, this approach often cannot achieve high accuracy for some complex traits. It is possible that SNP markers are not sufficient to predict these traits due to the missing heritability caused by other genetic variations such as microsatellite and copy number variation (CNV), which have been shown to affect disease and complex traits in humans and other species.

Results

In this study, CNVs were included in a SNP based genomic selection framework. A Nellore cattle dataset consisting of 2230 animals genotyped on BovineHD SNP array was used, and 9 weight and carcass traits were analyzed. A total of six models were implemented and compared based on their prediction accuracy. For comparison, three models including only SNPs were implemented: 1) BayesA model, 2) Bayesian mixture model (BayesB), and 3) a GBLUP model without polygenic effects. The other three models incorporating both SNP and CNV included 4) a Bayesian model similar to BayesA (BayesA+CNV), 5) a Bayesian mixture model (BayesB+CNV), and 6) GBLUP with CNVs modeled as a covariable (GBLUP+CNV). Prediction accuracies were assessed based on Pearson’s correlation between de-regressed EBVs (dEBVs) and direct genomic values (DGVs) in the validation dataset. For BayesA, BayesB and GBLUP, accuracy ranged from 0.12 to 0.62 across the nine traits. A minimal increase in prediction accuracy for some traits was noticed when including CNVs in the model (BayesA+CNV, BayesB+CNV, GBLUP+CNV).

Conclusions

This study presents the first genomic prediction study integrating CNVs and SNPs in livestock. Combining CNV and SNP marker information proved to be beneficial for genomic prediction of some traits in Nellore cattle.

Electronic supplementary material

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

CNVs are associated with genomic architecture in a songbird

Background

Understanding variation in genome structure is essential to understand phenotypic differences within populations and the evolutionary history of species. A promising form of this structural variation is copy number variation (CNV). CNVs can be generated by different recombination mechanisms, such as non-allelic homologous recombination, that rely on specific characteristics of the genome architecture. These structural variants can therefore be more abundant at particular genes ultimately leading to variation in phenotypes under selection. Detailed characterization of CNVs therefore can reveal evolutionary footprints of selection and provide insight in their contribution to phenotypic variation in wild populations.

Results

Here we use genotypic data from a long-term population of great tits (Parus major), a widely studied passerine bird in ecology and evolution, to detect CNVs and identify genomic features prevailing within these regions. We used allele intensities and frequencies from high-density SNP array data from 2,175 birds. We detected 41,029 CNVs concatenated into 8,008 distinct CNV regions (CNVRs). We successfully validated 93.75% of the CNVs tested by qPCR, which were sampled at different frequencies and sizes. A mother-daughter family structure allowed for the evaluation of the inheritance of a number of these CNVs. Thereby, only CNVs with 40 probes or more display segregation in accordance with Mendelian inheritance, suggesting a high rate of false negative calls for smaller CNVs. As CNVRs are a coarse-grained map of CNV loci, we also inferred the frequency of coincident CNV start and end breakpoints. We observed frequency-dependent enrichment of these breakpoints at homologous regions, CpG sites and AT-rich intervals. A gene ontology enrichment analyses showed that CNVs are enriched in genes underpinning neural, cardiac and ion transport pathways.

Conclusion

Great tit CNVs are present in almost half of the genes and prominent at repetitive-homologous and regulatory regions. Although overlapping genes under selection, the high number of false negatives make neutrality or association tests on CNVs detected here difficult. Therefore, CNVs should be further addressed in the light of their false negative rate and architecture to improve the comprehension of their association with phenotypes and evolutionary history.

Electronic supplementary material

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

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.

Distribution and Functionality of Copy Number Variation across European Cattle Populations

Copy number variation (CNV), which is characterized by large-scale losses or gains of DNA fragments, contributes significantly to genetic and phenotypic variation. Assessing CNV across different European cattle populations might reveal genetic changes responsible for phenotypic differences, which have accumulated throughout the domestication history of cattle as consequences of evolutionary forces that act upon them. To explore pattern of CNVs across European cattle, we genotyped 149 individuals, that represent different European regions, using the Illumina Bovine HD Genotyping array. A total of 9,944 autosomal CNVs were identified in 149 samples using a Hidden Markov Model (HMM) as employed in PennCNV. Animals originating from several breeds of British Isles, and Balkan and Italian regions, on average, displayed higher abundance of CNV counts than Dutch or Alpine animals. A total of 923 CNV regions (CNVRs) were identified by aggregating CNVs overlapping in at least two animals. The hierarchical clustering of CNVRs indicated low differentiation and sharing of high-frequency CNVRs between European cattle populations. Various CNVRs identified in the present study overlapped with olfactory receptor genes and genes related to immune system. In addition, we also detected a CNV overlapping the Kit gene in English longhorn cattle which has previously been associated with color-sidedness. To conclude, we provide a comprehensive overview of CNV distribution in genome of European cattle. Our results indicate an important role of purifying selection and genomic drift in shaping CNV diversity that exists between different European cattle populations.

β2-microglobulin gene duplication in cetartiodactyla remains intact only in pigs and possibly confers selective advantage to the species

Several β2-microglobulin (B2M) -bound protein complexes undertake key roles in various immune system pathways, including the neonatal Fc receptor (FcRn), cluster of differentiation 1 (CD1) protein, non-classical major histocompatibility complex (MHC), and well-known MHC class I molecules. Therefore, the duplication of B2M may lead to an increase in the biological competence of organisms to the environment. Based on the pig genome assembly SSC10.2, a segmental duplication of ~45.5 kb, encoding the entire B2M protein, was identified in pig chromosome 1. Through experimental validation, we confirmed the functional duplication of the B2M gene with a completely identical coding sequence between two copies in pigs. Considering the importance of B2M in the immune system, we performed the phylogenetic analysis of B2M duplication in ten mammalian species, confirming the presence of B2M duplication in cetartioldactyls, like cattle, sheep, goats, pigs and whales, but non-cetartiodactyl species, like mice, cats, dogs, horses, and humans. The density of long interspersed nuclear element (LINE) at the edges of duplicated blocks (39 to 66%) was found to be 2 to 3-fold higher than the average (20.12%) of the pig genome, suggesting its role in the duplication event. The B2M mRNA expression level in pigs was 12.71 and 7.57 times (2-ΔΔCt values) higher than humans and mice, respectively. However, we were unable to experimentally demonstrate the difference in the level of B2M protein because species specific anti-B2M antibodies are not available. We reported, for the first time, the functional duplication of the B2M gene in animals. The identification of partially remaining duplicated B2M sequences in the genomes of only cetartiodactyls indicates that the event was lineage specific. B2M duplication could be beneficial to the immune system of pigs by increasing the availability of MHC class I light chain protein, B2M, to complex with the proteins encoded by the relatively large number of MHC class I heavy chain genes in pigs. Further studies are necessary to address the biological meaning of increased expression of B2M.

Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems

The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently.

Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases

Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.

Identification of genomic loci associated with Rhodococcus equi susceptibility in foals

Pneumonia caused by Rhodococcus equi is a common cause of disease and death in foals. Although agent and environmental factors contribute to the incidence of this disease, the genetic factors influencing the clinical outcomes of R. equi pneumonia are ill-defined. Here, we performed independent single nucleotide polymorphism (SNP)- and copy number variant (CNV)-based genome-wide association studies to identify genomic loci associated with R. equi pneumonia in foals. Foals at a large Quarter Horse breeding farm were categorized into 3 groups: 1) foals with R. equi pneumonia (clinical group [N = 43]); 2) foals with ultrasonographic evidence of pulmonary lesions that never developed clinical signs of pneumonia (subclinical group [N = 156]); and, 3) foals without clinical signs or ultrasonographic evidence of pneumonia (unaffected group [N = 49]). From each group, 24 foals were randomly selected and used for independent SNP- and CNV-based genome-wide association studies (GWAS). The SNP-based GWAS identified a region on chromosome 26 that had moderate evidence of association with R. equi pneumonia when comparing clinical and subclinical foals. A joint analysis including all study foals revealed a 3- to 4-fold increase in odds of disease for a homozygous SNP within the associated region when comparing the clinical group with either of the other 2 groups of foals or their combination. The region contains the transient receptor potential cation channel, subfamily M, member 2 (TRPM2) gene, which is involved in neutrophil function. No associations were identified in the CNV-based GWAS. Collectively, these data identify a region on chromosome 26 associated with R. equi pneumonia in foals, providing evidence that genetic factors may indeed contribute to this important disease of foals.

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.

The challenges and importance of structural variation detection in livestock

Recent studies in humans and other model organisms have demonstrated that structural variants (SVs) comprise a substantial proportion of variation among individuals of each species. Many of these variants have been linked to debilitating diseases in humans, thereby cementing the importance of refining methods for their detection. Despite progress in the field, reliable detection of SVs still remains a problem even for human subjects. Many of the underlying problems that make SVs difficult to detect in humans are amplified in livestock species, whose lower quality genome assemblies and incomplete gene annotation can often give rise to false positive SV discoveries. Regardless of the challenges, SV detection is just as important for livestock researchers as it is for human researchers, given that several productive traits and diseases have been linked to copy number variations (CNVs) in cattle, sheep, and pig. Already, there is evidence that many beneficial SVs have been artificially selected in livestock such as a duplication of the agouti signaling protein gene that causes white coat color in sheep. In this review, we will list current SV and CNV discoveries in livestock and discuss the problems that hinder routine discovery and tracking of these polymorphisms. We will also discuss the impacts of selective breeding on CNV and SV frequencies and mention how SV genotyping could be used in the future to improve genetic selection.

Copy Number Variation in Chickens: A Review and Future Prospects

DNA sequence variations include nucleotide substitution, deletion, insertion, translocation and inversion. Deletion or insertion of a large DNA segment in the genome, referred to as copy number variation (CNV), has caught the attention of many researchers recently. It is believed that CNVs contribute significantly to genome variability, and thus contribute to phenotypic variability. In chickens, genome-wide surveys with array comparative genome hybridization (aCGH), SNP chip detection or whole genome sequencing have revealed a large number of CNVs. A large portion of chicken CNVs involves protein coding or regulatory sequences. A few CNVs have been demonstrated to be the determinant factors for single gene traits, such as late-feathering, pea-comb and dermal hyperpigmentation. The phenotypic effects of the majority of chicken CNVs are to be delineated.

Comparative Analysis of CNV Calling Algorithms: Literature Survey and a Case Study Using Bovine High-Density SNP Data

Copy number variations (CNVs) are gains and losses of genomic sequence between two individuals of a species when compared to a reference genome. The data from single nucleotide polymorphism (SNP) microarrays are now routinely used for genotyping, but they also can be utilized for copy number detection. Substantial progress has been made in array design and CNV calling algorithms and at least 10 comparison studies in humans have been published to assess them. In this review, we first survey the literature on existing microarray platforms and CNV calling algorithms. We then examine a number of CNV calling tools to evaluate their impacts using bovine high-density SNP data. Large incongruities in the results from different CNV calling tools highlight the need for standardizing array data collection, quality assessment and experimental validation. Only after careful experimental design and rigorous data filtering can the impacts of CNVs on both normal phenotypic variability and disease susceptibility be fully revealed.

Copy number variants in locally raised Chinese chicken genomes determined using array comparative genomic hybridization

Background

Copy number variants contribute to genetic variation in birds. Analyses of copy number variants in chicken breeds had focused primarily on those from commercial varieties with nothing known about the occurrence and diversity of copy number variants in locally raised Chinese chicken breeds. To address this deficiency, we characterized copy number variants in 11 chicken breeds and compared the variation among these breeds.

Results

We presented a detailed analysis of the copy number variants in locally raised Chinese chicken breeds identified using a customized comparative genomic hybridization array. We identified 833 copy number variants contained within 308 copy number variant regions. The median and mean sizes of the copy number variant regions were 14.6 kb and 35.1 kb, respectively. Of the copy number variant regions, 138 (45%) involved gain of DNA, 159 (52%) involved loss of DNA, and 11 (3%) involved both gain and loss of DNA. Principal component analysis and agglomerative hierarchical clustering revealed the close relatedness of the four locally raised chicken breeds, Shek-Ki, Langshan, Qingyuan partridge, and Wenchang. Biological process enrichment analysis of the copy number variant regions confirmed the greater variation among the four aforementioned varieties than among the seven other breeds studied.

Conclusion

Our description of the distribution of the copy number variants and comparison of the differences among the copy number variant regions of the 11 chicken breeds supplemented the information available concerning the copy number variants of other Chinese chicken breeds. In addition to its relevance for functional analysis, our results provided the first insight into how chicken breeds can be clustered on the basis of their genomic copy number variation.

Genome-wide copy number variant analysis in inbred chickens lines with different susceptibility to Marek's disease

Breeding of genetically resistant chickens to Marek’s disease (MD) is a vital strategy to poultry health. To find the markers underlying the genetic resistance to MD, copy number variation (CNV) was examined in inbred MD-resistant and -susceptible chicken lines. A total of 45 CNVs were found in four lines of chickens, and 28 were potentially involved in immune response and cell proliferation, etc. Importantly, two CNVs related with MD resistance were transmitted to descendent recombinant congenic lines that differ in susceptibility to MD. Our findings may lead to better strategies for genetic improvement of disease resistance in poultry.

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

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.7 kb 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 the UMD3.1 assembly'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 detection of copy number variations using SNP genotyping arrays in swine

Background

Copy Number Variations (CNVs) have been shown important in both normal phenotypic variability and disease susceptibility, and are increasingly accepted as another important source of genetic variation complementary to single nucleotide polymorphism (SNP). Comprehensive identification and cataloging of pig CNVs would be of benefit to the functional analyses of genome variation.

Results

In this study, we performed a genome-wide CNV detection based on the Porcine SNP60 genotyping data of 474 pigs from three pure breed populations (Yorkshire, Landrace and Songliao Black) and one Duroc × Erhualian crossbred population. A total of 382 CNV regions (CNVRs) across genome were identified, which cover 95.76Mb of the pig genome and correspond to 4.23% of the autosomal genome sequence. The length of these CNVRs ranged from 5.03 to 2,702.7kb with an average of 250.7kb, and the frequencies of them varied from 0.42 to 20.87%. These CNVRs contains 1468 annotated genes, which possess a great variety of molecular functions, making them a promising resource for exploring the genetic basis of phenotypic variation within and among breeds. To confirmation of these findings, 18 CNVRs representing different predicted status and frequencies were chosen for validation via quantitative real time PCR (qPCR). Accordingly, 12 (66.67%) of them was successfully confirmed.

Conclusions

Our results demonstrated that currently available Porcine SNP60 BeadChip can be used to capture CNVs efficiently. Our study firstly provides a comprehensive map of copy number variation in the pig genome, which would be of help for understanding the pig genome and provide preliminary foundation for investigating the association between various phenotypes and CNVs.

Molecular cloning and characterisation of ovine dual oxidase 2

The dual oxidases (DUOX1 and DUOX2) are NADPH-dependent hydrogen peroxide-producing enzymes that are reported to function in a physiological capacity and as a component of the mucosal immune response. We have previously reported increased expression of the DUOX2 gene in the gut mucosa of sheep in response to gastrointestinal nematode (GIN) challenge. In this paper, we report the cloning of the full-length ovine DUOX2 transcript, using a PCR based strategy. The ovine DUOX2 transcript includes an ORF of 4644 bases, and encodes a protein with 97% identity to the bovine sequence. We also cloned a fragment of DUOX1 (encompassing nucleotides 2692-2829), and the proximal promoter sequence of DUOX2. Through analysis of sequence data we have confirmed that DUOX1 and DUOX2 are co-located in a head to tail arrangement conserved across many species. Alignment of the sequences to the ovine genome predicts a location of this gene cluster on ovine chromosome 7. We quantified the expression of ovine DUOX1 and DUOX2 transcripts in 24 different sheep tissues, and discovered tissue specific expression signatures. DUOX2 was found to be most highly expressed in tissues of the gastrointestinal tract, while expression of DUOX1 predominated in the bladder. Rapid amplification of cDNA ends (RACE) analysis identified the existence of multiple 5' UTR variants in DUOX2, ranging in size from 32 to 242 nucleotides, with 3 distinct transcribed regions. Real time PCR quantification of the DUOX2 UTR variants revealed that these were differentially expressed between tissues, and at various stages of the response to GIN parasite infection. The collective evidence suggested a complex regulation of DUOX2, prompting a bioinformatic analysis of the proximal promoter regions of ovine DUOX2 to identify potential transcription factor binding sites (TFBS) that may explain the differences in the observed expression of the transcript variants of DUOX2. Possible transcription factor families that may regulate this process were identified as Kruppel-like factors (KLF), ETS-factors, erythroid growth receptor factors (EGRF) and myogenic differentiation factors (MYOD).

Copy number variation in the genomes of domestic animals

Copy number variation (CNV) might be one of the main contributors to phenotypic diversity and evolutionary adaptation in animals and plants, employing a wide variety of mechanisms, such as gene dosage and transcript structure alterations, to modulate organismal plasticity. In the past 4 years, considerable advances have been made in the characterization of the genomic architecture of CNV in domestic species. First, low-resolution CNV maps were produced for cattle, goat, sheep, pig, dog, chicken, duck and turkey, showing that these structural polymorphisms comprise a significant part of these genomes. Furthermore, CNVs have been associated with several pigmentation (white coat in horse, pig and sheep) and morphological (late feathering and pea comb in chicken) traits, as well as with susceptibility to a wide array of diseases and developmental disorders, for example osteopetrosis, anhidrotic ectodermal dysplasia, copper toxicosis, intersexuality, cone degeneration, periodic fever and dermoid sinus, among others. In the future, development of high-resolution tools for CNV detection and typing combined with the implementation of databases integrating CNV, QTL and gene expression data will be essential to identify and measure the impact of this source of structural variation on the many phenotypes that are relevant to animal breeders and veterinary practitioners.

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.