Genome-wide copy number variation in the bovine genome detected using low coverage sequence of popular beef breeds,

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.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

Identification of copy number variations in the genome of Dairy Gir cattle

Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.

Association of HSF1 gene copy number variation with growth traits in the Ashidan yak

Copy Number Variation (CNV) is the major manner for the variation of genome structure, which is associated with numerous important traits. The heat shock factor 1 (HSF1) gene is a stress response transcriptional regulator. It participates in the heat shock response, simultaneously participated in the development of tissue. The objective of this research was to explore the influence of CNV of the HSF1 gene on the growth traits of the Ashidan yak. In this study, the growth traits (withers height, body weight, chest girth, and body length) of 274 Ashidan yaks were divided into four stages (6, 12, 18, and 30 months old). Moreover, quantitative polymerase chain reaction (qPCR) was exploited for determining the HSF1 gene relative expression level, and SPSS software was utilized for the statistical analysis. The outcomes indicated that HSF1-CNV was significantly associated with body length (p < 0.05) and was extremely significant associated with withers height (p < 0.01) of 18-month-old Ashidan yaks. Besides, the HSF1 relative expression in heart and muscle was higher than that existed in other tissues (p < 0.01). The outcomes suggested that the CNV of HSF1 gene would affect the growth and development of the Ashidan yak, which is conducive to the early breeding of yak.

CNV analysis of VAMP7 gene reveals variation associated with growth traits in Chinese cattle

Copy number variant (CNV), a common genetic polymorphism, is closely related to the phenotypic variation traits of organisms. Vesicle-associated membrane protein 7 gene (VAMP7) codes a protein, which is a member of the SNARE proteins family and plays an important role in the process of intracellular vesicle transport. In this study, a total of four cattle breeds (Yunling cattle, Xianan cattle, Pinan cattle, Jiaxian red cattle) were used to investigate the copy numbers, and we found an association relationship between CNV of VAMP7 gene and growth traits of cattle by SPSS 20.0 software. The results showed that the CNV type of VAMP7 gene in four cattle breeds had the same distribution, Duplication type occupies a dominant position among the four varieties. In Yunling cattle, the Duplication type of VAMP7 is significantly related to the height at the hip cross (p < 0.05), Individuals with Duplication type commonly have less performance on growth and development, which indicates that the Duplication type of the VAMP7 gene may have a negative effect on cattle growth. Individuals with the other two CNV types may become the breeding direction of the VAMP7 gene. This study provided a new perspective and basic material for the molecular genetics of the CNV of the VAMP7 gene, and also promoted the breeding progress of Chinese local cattle.

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.

Copy number variations of the KAT6A gene are associated with body measurements of Chinese sheep breeds

Copy number variation (CNV) is one kind of genomic structure variations and presents as gains and losses of genomic fragments. More recently, we have made an atlas of CNV maps for livestock. In the future, it is a primary focus to determine the phenotypic effects of candidate CNVs. Lysine Acetyltransferase 6 A (KAT6A) is a protein coding gene and plays a critical role in many cellular processes. However, the effects of KAT6A CNVs on sheep body measurements remains unknown. In this study, we performed quantitative polymerase chain reaction (qPCR) to detect the presences and distributions of three CNV regions within KAT6A gene in 672 sheep from four Chinese breeds. Association analysis indicated that the three CNVs of KAT6A gene were significantly associated with body measurement(s) in Small-tailed Han sheep (STH) and Hu sheep (HU) (p < 0.05), while no effects on Large-tailed Han sheep (LTH) were observed (p > 0.05) were observed. Additionally, only one CNV was significantly associated with body measurement (body length) in Chaka sheep (CK) (p < 0.05). Our study provided evidence that the CNV(s) of KAT6A gene could be used as candidate marker(s) for molecular breedings of STH, HU, and CK breeds.

Genome-wide association study between copy number variation regions and carcass- and meat-quality traits in Nellore cattle

Context Indicine breeds are the main source of beef products in tropical and subtropical regions. However, genetic improvement for carcass- and meat-quality traits in zebu cattle have been limited and genomics studies concerning structural variations that influence these traits are essential. Aim The aim of this study was to perform a genome-wide association study between copy number variation regions (CNVRs) and carcass- and meat quality-traits in Nellore cattle. Methods In total, 3794 animals, males and females included, were genotyped using a 777962 single-nucleotide polymorphism platform of BovineHD BeadChip (777k; Illumina Inc.). Of these, 1751 Nellore bulls were slaughtered at 24 months of age for further carcass beef analysis. The following traits were studied: beef tenderness, marbling, rib-eye area, backfat thickness and meat colour (lightness, redness and yellowness). The CNV detection was conducted through PennCNV software. The association analyses were performed using CNVRuler software. Key results Several identified genomic regions were linked to quantitative trait loci associated with fat deposition (FABP7) and lipid metabolism (PPARA; PLA2 family; BCHE), extracellular matrix (INS; COL10A1), contraction (SLC34A3; TRDN) and muscle development (CAPZP). The gene-enrichment analyses highlighted biological mechanisms directly related to the metabolism and synthesis of lipids and fatty acids. Conclusions The large number of potential candidate genes identified within the CNVRs, as well as the functions and pathways identified, should help better elucidate the genetic mechanisms involved in the expression of beef and carcass traits in Nellore cattle. Several CNVRs harboured genes that might have a functional impact to improve the beef and carcass traits. Implications The results obtained contribute to upgrade the sensorial and organoleptic attributes of Nellore cattle and make feasible the genetic improvement of carcass- and meat-quality traits.

Genomic sequencing analysis reveals copy number variations and their associations with economically important traits in beef cattle

Copy number variation (CNV) represents a major source of genetic variation, which may have potentially large effects, including alternating gene regulation and dosage, as well as contributing to gene expression and risk for normal phenotypic variability. We carried out a comprehensive analysis of CNV based on whole genome sequencing in Chinese Simmental beef cattle. Totally, we found 9313 deletion and 234 duplication events, covering 147.5 Mb autosomal regions. Within them, 257 deletion events of high frequency overlapped with 193 known RefGenes. Among these genes, we observed several genes were related to economically important traits, like residual feed intake, immune responding, pregnancy rate and muscle differentiation. Using a locus-based analysis, we identified 11 deletions and 1 duplication, which were significantly associated with three traits including carcass weight, tenderloin and longissimus muscle area. Our sequencing-based study provided important insights into investigating the association of CNVs with important traits in beef cattle.

High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data

Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.

A large interactive visual database of copy number variants discovered in taurine cattle

BACKGROUND

Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations.

RESULTS

CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle.

CONCLUSIONS

We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.

A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

Copy number variation (CNV) in the IGF1R gene across four cattle breeds and its association with economic traits

The insulin-like growth factor 1 receptor (IGF1R) plays a vital role in immunomodulation and muscle and bone growth. The copy number variation (CNV) is believed to the reason for many complex phenotypic variations. In this paper, we statistically analyzed the copy number and the expression profiling in different tissue types of the IGF1R gene using the 422 samples from four Chinese beef cattle breeds, and the mRNA of IGF1R was widely expressed in nine tissue types of adult cattle (heart, liver, kidney, muscle, fat, stomach, spleen, lung and testis). Results of CNV and growth traits indicated that the IGF1R CNV was significantly associated with body weight and body height of Jinnan (JN) cattle and was significantly associated with body height and hucklebone width of Qinchuan (QC) cattle, making IGF1R CNV a promising molecular marker to improve meat production in beef cattle breeding. Bioinformatics predictions show that the CNV region is highly similar to the human genome, and there are a large number of transcription factors, DNase I hypersensitive sites, and high levels of histone acetylation, suggesting that this region may play a role in transcriptional regulation, providing directions for further study of the role of bovine CNV and economic traits.

The genetic background of clinical mastitis in Holstein-Friesian cattle

Mastitis is an inflammatory disease of the mammary gland, which has a significant economic impact and is an animal welfare concern. This work examined the association between single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) with the incidence of clinical mastitis (CM). Using information from 16 half-sib pairs of Holstein-Friesian cows (32 animals in total) we searched for genomic regions that differed between a healthy (no incidence of CM) and a mastitis-prone (multiple incidences of CM) half-sib. Three cows with average sequence depth of coverage below 10 were excluded, which left 13 half-sib pairs available for comparisons. In total, 191 CNV regions were identified, which were deleted in a mastitis-prone cow, but present in its healthy half-sib and overlapped in at least nine half-sib pairs. These regions overlapped with exons of 46 genes, among which APP (BTA1), FOXL2 (BTA1), SSFA2 (BTA2), OTUD3 (BTA2), ADORA2A (BTA17), TXNRD2 (BTA17) and NDUFS6 (BTA20) have been reported to influence CM. Moreover, two duplicated CNV regions present in nine healthy individuals and absent in their mastitis-affected half-sibs overlapped with exons of a cholinergic receptor nicotinic α 10 subunit on BTA15 and a novel gene (ENSBTAG00000008519) on BTA27. One CNV region deleted in nine mastitis-affected sibs overlapped with two neighbouring long non-coding RNA sequences located on BTA12. Single nucleotide polymorphisms with differential genotypes between a healthy and a mastitis-affected sib included 17 polymorphisms with alternate alleles in eight affected and healthy half-sib families. Three of these SNPs were located introns of genes: MET (BTA04), RNF122 (BTA27) and WRN (BTA27). In summary, structural polymorphisms in form of CNVs, putatively play a role in susceptibility to CM. Specifically, sequence deletions have a greater effect on reducing resistance against mastitis, than sequence duplications have on increasing resistance against the disease.

Commingled black and white cows (Bos taurus; Angus and Charolais) in high-elevation rangeland are differentially parasitised by Haematobia irritans

Globally horn flies (Haematobia irritans) are one of the most economically damaging parasites of beef cattle. These obligate blood-feeding external parasites take blood meals from cattle leading to blood loss, annoyance avoidance behaviours, and reductions in animal performance. Development of chemical resistance by horn flies suggests that novel management strategies are needed. More in-depth understanding of parasitism relative to hide colour and temperature, especially in a changing climate, may enhance animal production. In peak parasitism periods of 2016 and 2017, we measured horn fly loads on commingled black Angus and white Charolais cows in a cold high-elevation rangeland in Wyoming, USA. We tested how breed, year, and interactions explained horn fly parasitism and economic thresholds. In 2016 we also measured ambient and external cow temperatures to further elucidate thermal ecology mechanisms explaining horn fly hide colour preferences. Mean annual horn fly infestations were always four times greater or more on black cows than white cattle both years, but not all cattle reached economic thresholds all years and the breed by year interaction was not significant. Difference in horn fly preference for black cattle over white cattle in our cold high-elevation environment may be explained by greater absolute and relative external surface temperatures of black hided cows. Host colour and thermal preferences of horn flies could be incorporated into integrated pest management strategies that only treat darker hided cattle and producers in cold high-elevation environments conduct real-time monitoring to determine if treatments are even needed on a year-by-year basis.

Evaluation and Recommendations for Routine Genotyping Using Skim Whole Genome Re-sequencing in Canola

Whole genome sequencing offers genome wide, unbiased markers, and inexpensive library preparation. With the cost of sequencing decreasing rapidly, many plant genomes of modest size are amenable to skim whole genome resequencing (skim WGR). The use of skim WGR in diverse sample sets without the use of imputation was evaluated in silico in 149 canola samples representative of global diversity. Fastq files with an average of 10x coverage of the reference genome were used to generate skim samples representing 0.25x, 0.5x, 1x, 2x, 3x, 4x, and 5x sequencing coverage. Applying a pre-defined list of SNPs versus de novo SNP discovery was evaluated. As skim WGR is expected to result in some degree of insufficient allele sampling, all skim coverage levels were filtered at a range of minimum read depths from a relaxed minimum read depth of 2 to a stringent read depth of 5, resulting in 28 list-based SNP sets. As a broad recommendation, genotyping pre-defined SNPs between 1x and 2x coverage with relatively stringent depth filtering is appropriate for a diverse sample set of canola due to a balance between marker number, sufficient accuracy, and sequencing cost, but depends on the intended application. This was experimentally examined in two sample sets with different genetic backgrounds: 1x coverage of 1,590 individuals from 84 Australian spring type four-parent crosses aimed at maximizing diversity as well as one commercial F1 hybrid, and 2x coverage of 379 doubled haploids (DHs) derived from a subset of the four-parent crosses. To determine optimal coverage in a simpler genetic background, the DH sample sequence coverage was further down sampled in silico. The flexible and cost-effective nature of the protocol makes it highly applicable across a range of species and purposes.

In silico identification of conserved miRNAs and their selective target gene prediction in indicine (Bos indicus) cattle

The modern cattle was domesticated from aurochs, sharing its physiological traits into two subspecies Bos taurus and Bos indicus. MicroRNAs (miRNAs) are a class of non-coding short RNAs of ~22nt which have a key role in the regulation of many cellular and physiological processes in the animal. The current study was aimed to predict and annotate the potential mutations in indicine miRNAs throughout the genome using de novo and homology-based in silico approaches. Genome-wide mapping was performed in available indicine assembly by the homology-based approach and 768 miRNAs were recovered out of 808 reported taurine miRNAs belonging to 521 unique mature miRNA families. While 42 precursors were dropped due to lack of secondary miRNA structure, increasing stringency or decreasing similarity between the two genomes' miRNA. Increasing tendency of miRNAs incidence was observed on chr5, chr7, chr8, chr12 and chr21 with 19 polycistronic miRNA within 1-kilobase distance throughout the indicine genome. Notably, 12 miRNAs showed copy number variation. Eighteen miRNAs showed a mutation in their mature sequences in which eight were found in their seed region. Whilst in de novo based approach, 12 novel potential miRNAs on Y chromosome in indicine cattle along with a new miRNA (bind-miR-1264) on chrX were found. The final data set is annotated and explains the impending target genes that are responsible for enhanced immunity, heat tolerance and disease tolerance regulation in indicine. The study conforms to better understanding and perceptive approach towards indicine genome.

Misidentification of runs of homozygosity islands in cattle caused by interference with copy number variation or large intermarker distances

Background

Runs of homozygosity (ROH) islands are stretches of homozygous sequence in the genome of a large proportion of individuals in a population. Algorithms for the detection of ROH depend on the similarity of haplotypes. Coverage gaps and copy number variants (CNV) may result in incorrect identification of such similarity, leading to the detection of ROH islands where none exists. Misidentified hemizygous regions will also appear as homozygous based on sequence variation alone. Our aim was to identify ROH islands influenced by marker coverage gaps or CNV, using Illumina BovineHD BeadChip (777 K) single nucleotide polymorphism (SNP) data for Austrian Brown Swiss, Tyrol Grey and Pinzgauer cattle.

Methods

ROH were detected using clustering, and ROH islands were determined from population inbreeding levels for each marker. CNV were detected using a multivariate copy number analysis method and a hidden Markov model. SNP coverage gaps were defined as genomic regions with intermarker distances on average longer than 9.24 kb. ROH islands that overlapped CNV regions (CNVR) or SNP coverage gaps were considered as potential artefacts. Permutation tests were used to determine if overlaps between CNVR with copy losses and ROH islands were due to chance. Diversity of the haplotypes in the ROH islands was assessed by haplotype analyses.

Results

In Brown Swiss, Tyrol Grey and Pinzgauer, we identified 13, 22, and 24 ROH islands covering 26.6, 389.0 and 35.8 Mb, respectively, and we detected 30, 50 and 71 CNVR derived from CNV by using both algorithms, respectively. Overlaps between ROH islands, CNVR or coverage gaps occurred for 7, 14 and 16 ROH islands, respectively. About 37, 44 and 52% of the ROH islands coverage in Brown Swiss, Tyrol Grey and Pinzgauer, respectively, were affected by copy loss. Intersections between ROH islands and CNVR were small, but significantly larger compared to ROH islands at random locations across the genome, implying an association between ROH islands and CNVR. Haplotype diversity for reliable ROH islands was lower than for ROH islands that intersected with copy loss CNVR.

Conclusions

Our findings show that a significant proportion of the ROH islands in the bovine genome are artefacts due to CNV or SNP coverage gaps.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0414-x) contains supplementary material, which is available to authorized users.

Detection of genome-wide structural variations in the Shanghai Holstein cattle population using next-generation sequencing

Objective

The Shanghai Holstein cattle breed is susceptible to severe mastitis and other diseases due to the hot weather and long-term humidity in Shanghai, which is the main distribution centre for providing Holstein semen to various farms throughout China. Our objective was to determine the genetic mechanisms influencing economically important traits, especially diseases that have huge impact on the yield and quality of milk as well as reproduction.

Methods

In our study, we detected the structural variations of 1,092 Shanghai Holstein cows by using next-generation sequencing. We used the DELLY software to identify deletions and insertions, cn.MOPS to identify copy-number variants (CNVs). Furthermore, we annotated these structural variations using different bioinformatics tools, such as gene ontology, cattle quantitative trait locus (QTL) database and ingenuity pathway analysis (IPA).

Results

The average number of high-quality reads was 3,046,279. After filtering, a total of 16,831 deletions, 12,735 insertions and 490 CNVs were identified. The annotation results showed that these mapped genes were significantly enriched for specific biological functions, such as disease and reproduction. In addition, the enrichment results based on the cattle QTL database showed that the number of variants related to milk and reproduction was higher than the number of variants related to other traits. IPA core analysis found that the structural variations were related to reproduction, lipid metabolism, and inflammation. According to the functional analysis, structural variations were important factors affecting the variation of different traits in Shanghai Holstein cattle. Our results provide meaningful information about structural variations, which may be useful in future assessments of the associations between variations and important phenotypes in Shanghai Holstein cattle.

Conclusion

Structural variations identified in this study were extremely different from those of previous studies. Many structural variations were found to be associated with mastitis and reproductive system diseases; these results are in accordance with the characteristics of the environment that Shanghai Holstein cattle experience.

Landscape of copy number variations in Bos taurus: individual - and inter-breed variability

Background

The number of studies of Copy Number Variation in cattle has increased in recent years. This has been prompted by the increased availability of data on polymorphisms and their relationship with phenotypes. In addition, livestock species are good models for some human phenotypes. In the present study, we described the landscape of CNV driven genetic variation in a large population of 146 individuals representing 13 cattle breeds, using whole genome DNA sequence.

Results

A highly significant variation among all individuals and within each breed was observed in the number of duplications (P < 10⁻¹⁵) and in the number of deletions (P < 10⁻¹⁵). We also observed significant differences between breeds for duplication (P = 0.01932) and deletion (P = 0.01006) counts. The same variation CNV length - inter-individual and inter-breed differences were significant for duplications (P < 10⁻¹⁵) and deletions (P < 10⁻¹⁵). Moreover, breed-specific variants were identified, with the largest proportion of breed-specific duplications (9.57%) found for Fleckvieh and breed-specific deletions found for Brown Swiss (5.00%). Such breed-specific CNVs were predominantly located in intragenic regions, however in Simmental, one deletion present in five individuals was found in the coding sequence of a novel gene ENSBTAG00000000688 on chromosome 18. In Brown Swiss, Norwegian Red and Simmental breed-specific deletions were located within KIT and MC1R genes, which are responsible for a coat colour. The functional annotation of coding regions underlying the breed-specific CNVs showed that in Norwegian Red, Guernsey, and Simmental significantly under- and overrepresented GO terms were related to chemical stimulus involved in sensory perception of smell and the KEGG pathways for olfactory transduction. In addition, specifically for the Norwegian Red breed, the dopaminergic synapse KEGG pathway was significantly enriched within deleted parts of the genome.

Conclusions

The CNV landscape in Bos taurus genome revealed by this study was highly complex, with inter-breed differences, but also a significant variation within breeds. The former, may explain some of the phenotypic differences among analysed breeds, and the latter contributes to within-breed variation available for selection.

Electronic supplementary material

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

Identification of copy number variation in French dairy and beef breeds using next-generation sequencing

Background

Copy number variations (CNV) are known to play a major role in genetic variability and disease pathogenesis in several species including cattle. In this study, we report the identification and characterization of CNV in eight French beef and dairy breeds using whole-genome sequence data from 200 animals. Bioinformatics analyses to search for CNV were carried out using four different but complementary tools and we validated a subset of the CNV by both in silico and experimental approaches.

Results

We report the identification and localization of 4178 putative deletion-only, duplication-only and CNV regions, which cover 6% of the bovine autosomal genome; they were validated by two in silico approaches and/or experimentally validated using array-based comparative genomic hybridization and single nucleotide polymorphism genotyping arrays. The size of these variants ranged from 334 bp to 7.7 Mb, with an average size of ~ 54 kb. Of these 4178 variants, 3940 were deletions, 67 were duplications and 171 corresponded to both deletions and duplications, which were defined as potential CNV regions. Gene content analysis revealed that, among these variants, 1100 deletions and duplications encompassed 1803 known genes, which affect a wide spectrum of molecular functions, and 1095 overlapped with known QTL regions.

Conclusions

Our study is a large-scale survey of CNV in eight French dairy and beef breeds. These CNV will be useful to study the link between genetic variability and economically important traits, and to improve our knowledge on the genomic architecture of cattle.

Electronic supplementary material

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

Evolutionary and Functional Features of Copy Number Variation in the Cattle Genome

Genomic structural variations are an important source of genetic diversity. Copy number variations (CNVs), gains and losses of large regions of genomic sequence between individuals of a species, have been associated with a wide variety of phenotypic traits. However, in cattle, as well as many other species, relatively little is understood about CNV, including frequency of CNVs in the genome, sizes, and locations, chromosomal properties, and evolutionary processes acting to shape CNV. In this work, we focused on copy number variation in the bovine genome, with the aim to detect CNVs in Bos taurus coding sequence and explore potential evolutionary mechanisms shaping these CNV. We identified and characterized CNV regions by utilizing exome sequence from 175 influential sires used in the Germplasm Evaluation project, representing 10 breeds. We examined various evolutionary and functional aspects of these CNVs, including selective constraint on CNV-overlapped genes, centrality of CNV genes in protein-protein interaction networks, and tissue-specific expression of CNV genes. Patterns of CNV in the Bos taurus genome reveal that reduced functional constraint and mutational bias may play a prominent role in shaping this type of structural variation.