Copy number variation in the cattle genome

Impact of STAT5A-CNVs on growth traits in Chinese beef cattle breeds

CNVs, which are a type of structural variation, make a substantial impact on diverse characteristics in multiple species. Q-PCR and data association analysis were used for STAT5A gene copy in this study. This study aimed to investigate the copy number variation (CNV) of the STAT5A gene in seven Chinese cattle breeds, namely Qinchuan cattle, Xianan cattle, Yunling cattle, Ji'an cattle, Jiaxian Red cattle, Qaidam cattle, and Guyuan yellow cattle. Blood samples were collected for CNV typing, and the correlation between CNV type and growth traits was analyzed using SPSS 23.0 software and ANOVA. The findings revealed variations in the distribution of different copy number types among the different cattle breeds. Furthermore, association analysis demonstrated a positive impact of CNV in the STAT5A gene on cattle growth: in the JX, individuals with duplication types exhibited superior performance in terms of rump length (P < 0.05). Conversely, normal GY cattle demonstrated better body height and abdomen circumference (P < 0.05), while QD cattle exhibited a significant correlation between weight and body length with normal individuals (P < 0.05). Moreover, QC bovine duplication individuals outperformed other types, with copy number variation significantly associated with chest depth, chest width, and body length (P < 0.05). The results validate the correlation between copy number variation (CNV) of the STAT5A gene and growth characteristics in five different cattle breeds, providing a reliable benchmark for the purpose of cattle breeding.

Advancements in copy number variation screening in herbivorous livestock genomes and their association with phenotypic traits

Copy number variations (CNVs) have garnered increasing attention within the realm of genetics due to their prevalence in human, animal, and plant genomes. These structural genetic variations have demonstrated associations with a broad spectrum of phenotypic diversity, economic traits, environmental adaptations, epidemics, and other essential aspects of both plants and animals. Furthermore, CNVs exhibit extensive sequence variability and encompass a wide array of genomes. The advancement and maturity of microarray and sequencing technologies have catalyzed a surge in research endeavors pertaining to CNVs. This is particularly prominent in the context of livestock breeding, where molecular markers have gained prominence as a valuable tool in comparison to traditional breeding methods. In light of these developments, a contemporary and comprehensive review of existing studies on CNVs becomes imperative. This review serves the purpose of providing a brief elucidation of the fundamental concepts underlying CNVs, their mutational mechanisms, and the diverse array of detection methods employed to identify these structural variations within genomes. Furthermore, it seeks to systematically analyze the recent advancements and findings within the field of CNV research, specifically within the genomes of herbivorous livestock species, including cattle, sheep, horses, and donkeys. The review also highlighted the role of CNVs in shaping various phenotypic traits including growth traits, reproductive traits, pigmentation and disease resistance etc., in herbivorous livestock. The main goal of this review is to furnish readers with an up-to-date compilation of knowledge regarding CNVs in herbivorous livestock genomes. By integrating the latest research findings and insights, it is anticipated that this review will not only offer pertinent information but also stimulate future investigations into the realm of CNVs in livestock. In doing so, it endeavors to contribute to the enhancement of breeding strategies, genomic selection, and the overall improvement of herbivorous livestock production and resistance to diseases.

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.

Descriptive Statistics and Genome-Wide Copy Number Analysis of Milk Production Traits of Jiangsu Chinese Holstein Cows

Milk production traits are the most important quantitative economic traits in dairy cow production; improving the yield and quality of milk is an important way to ensure the production efficiency of the dairy industry. This study carried out a series of in-depth statistical genetics studies and molecular analyses on the Chinese Holstein cows in the Jiangsu Province, such as descriptive statistics and copy number variation analysis. A genetic correlation, phenotypic correlation, and descriptive statistical analysis of five milk production traits (milk yield, milk fat percentage, milk fat yield, milk protein percentage, and milk protein yield) of the dairy cows were analyzed using the SPSS and DMU software. Through quality control, 4173 cows and their genomes were used for genomic study. Then, SNPs were detected using DNA chips, and a copy number variation (CNV) analysis was carried out to locate the quantitative trait loci (QTL) of the milk production traits by Perl program software Penn CNV and hidden Markov model (HMM). The phenotypic means of the milk yield, milk fat percentage, milk fat mass, milk protein percentage, and milk protein mass at the first trimester were lower than those at the other trimesters by 8.821%, 1.031%, 0.930%, 0.003%, and 0.826%, respectively. The five milk production traits showed a significant phenotypic positive correlation (p < 0.01) and a high genetic positive correlation among the three parities. Based on the GGPBovine 100 K SNP data, QTL-detecting research on the fist-parity milk performance of dairy cows was carried out via the CNV. We identified 1731 CNVs and 236 CNVRs in the 29 autosomes of 984 Holstein dairy cows, and 19 CNVRs were significantly associated with the milk production traits (p < 0.05). These CNVRs were analyzed via a bioinformatics analysis; a total of 13 gene ontology (GO) terms and 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to lipid metabolism, amino acid metabolism, and cellular catabolic processes. This study provided a theoretical basis for the molecular-marker-assisted selection of dairy cows by developing descriptive statistics on the milk production traits of dairy cows and by locating the QTL and functional genes that affect the milk production traits of first-born dairy cows. The results describe the basic status of the milk production traits of the Chinese Holstein cows in Jiangsu and locate the QTL and functional genes that affect the milk production traits of the first-born cows, providing a theoretical basis for the molecular-marker-assisted selection of dairy cows.

Copy number variation of bovine S100A7 as a positional candidate affected body measurements

Beef production is closely related to the national economy and the attention has been paid to the improvement of beef cattle by molecular markers associated. Copy number variations (CNVs) recently have been gained many researches and recognized as an important source of genetic variation. Extensive studies have indicated that CNVs have effects on a large range of economic traits by a wide range of gene copy number alteration. S100A7 is a member of S100 family which is a famous family of Ca2+-binding proteins. S100A7 plays a crucial role in many important phenotypes (progress) including inflammatory diseases, psoriasis, obesity, etc. The aim of our study was to explore the phenotypic effects of CNV located in the S100A7 gene of bovine chromosome 3. We detected S100A7 CNV by qPCR in different cattle breeds, including Qinchuan cattle, Yunling cattle, Xianan cattle and a crossbred group Pinan. The copy number was identified as gain, normal and loss type, our results showed that the gain type was the main type in three types of S100A7 CNV of the whole tested breeds. After CNV detection, association analysis between S100A7 CNV and growth traits was carried out in four cattle breeds. We found significant effects of the CNV on cattle growth traits with differently preferred CNV types such as gain type with better chest depth (p = 0.043) in QC, loss type with better body length (p = 0.008) and rump width (p = 0.014) in YL, normal with better chest girth (p = 0.001), gain with better waist width (p = 0.001) and rump width (p = 0.044) in PN. These results suggested that the S100A7 CNV could affect the phenotypic traits and be used as a promising genetic marker for cattle molecular breeding.

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.

Genome-Wide Detection of Copy Number Variations Associated with Miniature Features in Horses

Copy number variations (CNVs) are crucial structural genomic variants affecting complex traits in humans and livestock animals. The current study was designed to conduct a comprehensive comparative copy number variation analysis among three breeds, Debao (DB), Baise (BS), and Warmblood (WB), with a specific focus on identifying genomic regions associated with miniature features in horses. Using whole-genome next-generation resequencing data, we identified 18,974 CNVs across 31 autosomes. Among the breeds, we found 4279 breed-specific CNV regions (CNVRs). Baise, Debao, and Warmblood displayed 2978, 986, and 895 distinct CNVRs, respectively, with 202 CNVRs shared across all three breeds. After removing duplicates, we obtained 1545 CNVRs from 26 horse genomes. Functional annotation reveals enrichment in biological functions, including antigen processing, cell metabolism, olfactory conduction, and nervous system development. Debao horses have 970 genes overlapping with CNVRs, possibly causing their small size and mountainous adaptations. We also found that the genes GHR, SOX9, and SOX11 may be responsible for the miniature features of the Debao horse by analyzing their overlapping CNVRs. Overall, this study offers valuable insights into the widespread presence of CNVs in the horse genome. The findings contribute to mapping horse CNVs and advance research on unique miniature traits observed in the Debao horse.

Identification and characterisation of de novo germline structural variants in two commercial pig lines using trio-based whole genome sequencing

BACKGROUND

De novo mutations arising in the germline are a source of genetic variation and their discovery broadens our understanding of genetic disorders and evolutionary patterns. Although the number of de novo single nucleotide variants (dnSNVs) has been studied in a number of species, relatively little is known about the occurrence of de novo structural variants (dnSVs). In this study, we investigated 37 deeply sequenced pig trios from two commercial lines to identify dnSVs present in the offspring. The identified dnSVs were characterised by identifying their parent of origin, their functional annotations and characterizing sequence homology at the breakpoints.

RESULTS

We identified four swine germline dnSVs, all located in intronic regions of protein-coding genes. Our conservative, first estimate of the swine germline dnSV rate is 0.108 (95% CI 0.038-0.255) per generation (one dnSV per nine offspring), detected using short-read sequencing. Two detected dnSVs are clusters of mutations. Mutation cluster 1 contains a de novo duplication, a dnSNV and a de novo deletion. Mutation cluster 2 contains a de novo deletion and three de novo duplications, of which one is inverted. Mutation cluster 2 is 25 kb in size, whereas mutation cluster 1 (197 bp) and the other two individual dnSVs (64 and 573 bp) are smaller. Only mutation cluster 2 could be phased and is located on the paternal haplotype. Mutation cluster 2 originates from both micro-homology as well as non-homology mutation mechanisms, where mutation cluster 1 and the other two dnSVs are caused by mutation mechanisms lacking sequence homology. The 64 bp deletion and mutation cluster 1 were validated through PCR. Lastly, the 64 bp deletion and the 573 bp duplication were validated in sequenced offspring of probands with three generations of sequence data.

CONCLUSIONS

Our estimate of 0.108 dnSVs per generation in the swine germline is conservative, due to our small sample size and restricted possibilities of dnSV detection from short-read sequencing. The current study highlights the complexity of dnSVs and shows the potential of breeding programs for pigs and livestock species in general, to provide a suitable population structure for identification and characterisation of dnSVs.

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.

Genome-wide detection of copy number variation in American mink using whole-genome sequencing

BACKGROUND

Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta.

RESULTS

A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3).

CONCLUSIONS

This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.

Comparative analysis of allele frequencies for DNA polymorphisms associated with disease and economically important traits in the genomes of Russian and foreign cattle breeds

The genetic makeup of a breed including its genetic differences from other breeds determines its appearance and characteristics, including economically important traits and resistance to pathologies. To date, many loci controlling significant phenotypes have been identified, which is successfully used in the world practice of marker-assisted selection to improve breed properties. The aim of this study was a comparative analysis of frequencies for known causative nucleotide substitutions, insertions and deletions associated with disease and economically important traits in Russian and foreign cattle breeds. As a result, we identified frequencies of these DNA polymorphisms in the populations of Russian cattle breeds, compared them with those of foreign populations of the same breed, as well as other foreign breeds. Our results indicate similarities in frequencies for most of such alleles within breeds (populations of Russian and foreign breeding), as well as the relationship between the causative allele prevalence and the presence of phenotypic traits under the effect. We also found an excess of some undesirable alleles in the Russian cattle populations, which should be paid attention to when designing breeding programs. We found that the alleles increasing fertility in the Hereford breed have a higher frequency in the Russian Hereford population compared to the foreign counterpart. Interestingly, unlike for the European breeds, for Asian Turano-Mongolian Wagyu and Yakut cattle, there was a less clear link between phenotypic traits and frequencies of known causative alleles. Our work points to specific genetic variants that could be used to improve and/or maintain the performance of certain cattle breeds bred in the Russian Federation.

The copy number variation of DMBT1 gene effects body traits in two Chinese cattle breeds

Copy number variations (CNVs) belong to mutations in the genome level with loci in the region of genic or intergenic. It is through different effects (such as position effect and dose effect) that influence complex traits and diseases. Deleted in Malignant Brain Tumors 1 (DMBT1) gene is a member of the scavenger receptor cysteine-rich super family. In cattle, this gene has been associated with the susceptibility to bovine tuberculosis. In this study, a new CNV was found in DMBT1 gene of Chinese cattle breeds and tested in two different Chinese cattle breeds (Jiaxian red and Pinan) for frequency distribution analysis. Besides, the body size data such as body length, body height, chest girth, chest width, rump length, and rump girth for Jiaxian (JX) and Pinan (PN) cattle were collected and associated with the newly identified CNV. The CNV was significantly associated with the body length and chest girth of JX cattle, and the rump length of PN cattle (P < 0.05). Furthermore, the expression profile of the DMBT1 gene was tested in calves' tissues and the myoblasts differentiation. It was found that the DMBT1 gene expression was high in tuberculosis susceptible tissues (liver and lungs) at the calf stage and high in myoblast early differentiation. These tests were done using the qPCR method. As the result, the CNV of DMBT1 gene could be used as a candidate marker for bovine growth and health in marker-assisted selection (MAS) breeding.

The relationship between MFN1 copy number variation and growth traits of beef cattle

Copy number variation, as a kind of genetic submicroscopic structural variation, refers to the deletion or repetition of a large segment of genomic DNA, involving a segment size ranging from 50 bp to several MB. Mitochondrial fusion protein (MFN1) gene regulates the fusion of mitochondrial outer membrane in cells and maintains the dynamic needs of reticular mitochondria in cells. In this study, we conducted to tested the dstribution characteristics of MFN1-CNV in 522 cattles across Xianan cattle (XN), Pinan cattle (PN), Qinchuan cattle (QC), Jiaxian cattle (JX), Yunling cattle (YL), and correlated it with phenotypic traits. Then we observed the expression of MFN1 in various tissues of QC cattle (n = 3), and the expression levels were higher in lung and muscle. The results showed that there was significant correlation between MFN1 gene CNV and hucklebone width of QC cattle, hip width and height at sacrum of JX red cattle, chest width and rump length of YL cattle (P < 0.05). Individuals with duplication type were better than the type of normal or deletion in phenotypic traits. In conclusion, our data showed the correlation between MFN1 gene and growth traits of Chinese cattle. MFN1 gene can be used as a molecular marker for cattle selection and breeding, and accelerate the improvement of Chinese 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.

[Signatures of selection and candidate genes for adaptation to extreme environmental factors in the genomes of Turano-Mongolian cattle breeds]

Changes in the environment force populations of organisms to adapt to new conditions, either through phenotypic plasticity or through genetic or epigenetic changes. Signatures of selection, such as specific changes in the frequency of alleles and haplotypes, as well as the reduction or increase in genetic diversity, help to identify changes in the cattle genome in response to natural and artificial selection, as well as loci and genetic variants directly affecting adaptive and economically important traits. Advances in genetics and biotechnology enable a rapid transfer of unique genetic variants that have originated in local cattle breeds in the process of adaptation to local environments into the genomes of cosmopolitan high-performance breeds, in order to preserve their outstanding performance in new environments. It is also possible to use genomic selection approach to increase the frequency of already present adaptive alleles in cosmopolitan breeds. The review examines recent work on the origin and evolution of Turano-Mongolian cattle breeds, adaptation of Turano-Mongolian cattle to extreme environments, and summarizes available information on potential candidate genes for climate adaptation of Turano-Mongolian breeds, including cold resistance genes, immune response genes, and high-altitude adaptation genes. The authors conclude that the current literature data do not provide preference to one of the two possible scenarios of Turano-Mongolian breed origins: as a result of the domestication of a wild aurochs at East Asia or as a result of the migration of taurine proto-population from the Middle East. Turano-Mongolian breeds show a high degree of adaptation to extreme climatic conditions (cold, heat, lack of oxygen in the highlands) and parasites (mosquitoes, ticks, bacterial and viral infections). As a result of high-density genotyping and sequencing of genomes and transcriptomes, prospective candidate genes and genetic variants involved in adaptation to environmental factors have recently been identified.

Population differentiated copy number variation of Bos taurus, Bos indicus and their African hybrids

Background

CNV comprises a large proportion in cattle genome and is associated with various traits. However, there were few population-scale comparison studies on cattle CNV.

Results

Here, autosome-wide CNVs were called by read depth of NGS alignment result and copy number variation regions (CNVRs) defined from 102 Eurasian taurine (EAT) of 14 breeds, 28 Asian indicine (ASI) of 6 breeds, 22 African taurine (AFT) of 2 breeds, and 184 African humped cattle (AFH) of 17 breeds. The copy number of every CNVRs were compared between populations and CNVRs with population differentiated copy numbers were sorted out using the pairwise statistics V_ST and Kruskal-Wallis test. Three hundred sixty-two of CNVRs were significantly differentiated in both statistics and 313 genes were located on the population differentiated CNVRs.

Conclusion

For some of these genes, the averages of copy numbers were also different between populations and these may be candidate genes under selection. These include olfactory receptors, pathogen-resistance, parasite-resistance, heat tolerance and productivity related genes. Furthermore, breed- and individual-level comparison was performed using the presence or copy number of the autosomal CNVRs. Our findings were based on identification of CNVs from short Illumina reads of 336 individuals and 39 breeds, which to our knowledge is the largest dataset for this type of analysis and revealed important CNVs that may play a role in cattle adaption to various environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07808-7.

Detection of copy number variants in African goats using whole genome sequence data

Background

Copy number variations (CNV) are a significant source of variation in the genome and are therefore essential to the understanding of genetic characterization. The aim of this study was to develop a fine-scaled copy number variation map for African goats. We used sequence data from multiple breeds and from multiple African countries.

Results

A total of 253,553 CNV (244,876 deletions and 8677 duplications) were identified, corresponding to an overall average of 1393 CNV per animal. The mean CNV length was 3.3 kb, with a median of 1.3 kb. There was substantial differentiation between the populations for some CNV, suggestive of the effect of population-specific selective pressures. A total of 6231 global CNV regions (CNVR) were found across all animals, representing 59.2 Mb (2.4%) of the goat genome. About 1.6% of the CNVR were present in all 34 breeds and 28.7% were present in all 5 geographical areas across Africa, where animals had been sampled. The CNVR had genes that were highly enriched in important biological functions, molecular functions, and cellular components including retrograde endocannabinoid signaling, glutamatergic synapse and circadian entrainment.

Conclusions

This study presents the first fine CNV map of African goat based on WGS data and adds to the growing body of knowledge on the genetic characterization of goats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07703-1.

Insights into adaption and growth evolution: a comparative genomics study on two distinct cattle breeds from Northern and Southern China

Mongolian cattle (MG, Bos taurus) and Minnan cattle (MN, Bos indicus) are two different breeds of Chinese indigenous cattle, representing North type and South type, respectively. However, their value and potential have not yet been discovered at the genomic level. In this study, 26 individuals of MN and MG were sequenced for the first time at an average of 13.9- and 12.8-fold, respectively. Large numbers of different variations were identified. In addition, the analyses of phylogenetic and population structure showed that these two cattle breeds are distinct from each other, and results of linkage disequilibrium analysis revealed that these two cattle breeds have undergone various degrees of intense natural or artificial selection. Subsequently, 496 and 306 potential selected genes (PSRs) were obtained in MN and MG, containing 1,096 and 529 potential selected genes (PSGs), respectively. These PSGs, together with the analyzed copy number variation (CNV)-related genes, showed potential relations with their phenotypic characteristics, including environmental adaptability (e.g., DVL2, HSPA4, CDHR4), feed efficiency (e.g., R3HDM1, PLAG1, XKR4), and meat/milk production (e.g., PDHB, LEMD3, APOF). The results of this study help to gain new insights into the genetic characteristics of two distinct cattle breeds and will contribute to future cattle breeding.

Copy number variation of CADM2 gene revealed its association with growth traits across Chinese Capra hircus (goat) populations

Copy number variations (CNVs) are the wide structural variations ranging from 50 bp to several Mb at genome which can affect gene expression and further impacting growth and development traits of livestock. Comparing with single nucleotide polymorphisms (SNPs), CNVs can better explain the genetic and phenotypic diversity, are increasingly important in biological research. As a member of immunoglobulin super-family, cell adhesion molecule 2 (CADM2) plays a vital role in cancer development and metabolic regulation. Here, we tested the CNV of CADM2 gene in 443 goats across five breeds (Guizhou white goat, GZW; Guizhou black goat, GZB; Africa Nubian goat, AN; Boer goat × Huai goat, BH; Boer goat, BG) and detected its association with phenotypic traits. Subsequently, we analyzed the CADM2 gene expression level in different tissues of NB goats (n = 3, Nubian × Black) and the transcriptional expression in lung is much higher than others. The results showed that the CNV of CADM2 has a significant association with withers height and body length in GZB goat (P < 0.01), in which individuals with type of deletion were superior to those with duplication or normal type in term of body hight and body length (P < 0.01). In summary, this study confirmed the association between CNV of CADM2 gene and growth traits, and our research data indicated the CADM2-CNV may considered as a prospective candidate for the molecular marker-assisted selection breeding of goat growth traits, which conducived to accelerating the genetic amelioration in Chinese goats.

The Association of the Copy Number Variation of the MLLT10 Gene with Growth Traits of Chinese Cattle

Simple Summary

Copy number variation is a common genetic polymorphism, mainly represented by submicroscopic levels of deletion and duplication, caused by rearrangement of the genome. It is well-known that the copy number variation of a gene is associated with growth traits of livestock. In this study, we detected the correlation between the copy number variation of the the MLLT10 gene and the growth traits of Chinese yellow cattle. We found that the copy number variation of the MLLT10 gene has a significant influence on hip width, rump length, hucklebone width, and cannon bone circumference of some Chinese yellow cattle breeds. The results provide preliminary suggestions for Chinese yellow cattle breeding and new insights about the future of copy number variation (CNV) as a new promising molecular marker in animal breeding.

Abstract

Copy number variation is a part of genomic structural variation and has caused widespread concern. According to the results of high-throughput screening of the MLLT10 gene, we found that the copy number variation region of the MLLT10 gene was correlated with bovine growth traits. We aimed to detect the MLLT10 gene copy number variation and provide materials for the Chinese yellow cattle breed. In this study, the SPSS software was used to analyze the correlation among the copy number type of six different cattle breeds (i.e., Qinchuan, Xianan, Jiaxian, Yanbian, Sinan, Yunling) and the corresponding growth traits. The results showed the following: In Qinchuan cattle, the copy number duplication type was greater than the deletion and normal types; in Xianan cattle, the copy number duplication and normal types were less as compared with the deletion type; and in Yunling cattle, the frequency of the duplication type was dominant among the three types of copy number variants. The correlation analysis result showed that there is a significant correlation between the copy number variation (CNV) of the MLLT10 gene and the growth traits of three cattle breeds. Furthermore, correlation analysis showed that MLLT10 CNV had positive effects on growth traits such as hip width, rump length, hucklebone width, and cannon bone circumference (p < 0.05). This study provides a basis for the molecular-assisted marker breeding of cattle and contributes to the breeding of cattle.

Copy number variation detection in Chinese indigenous cattle by whole genome sequencing

Copy number variation (CNV) refers to a kind of structural variation, having functional and evolutionary effects on phenotypes. Thus far, further elucidation of the CNVs in different Chinese indigenous cattle breeds by whole genome sequencing have yet not been done. In this study, a comprehensive genomic analysis was performed on 75 cattle individuals including six Chinese indigenous cattle breeds and two non-native specialized beef cattle breeds. Based on the 11,486 CNVRs discovered, population analysis was performed, showed that all the cattle breeds clustered in to three clades, consistent with their lineages Bos taurus, Bos taurus × Bos indicus and Bos indicus. Importantly, a set of CNVRs related genes were found to be associated with the traits of interest, which include meat production or quality (CAST, ACTC1, etc.), adaption (BLA-DQB, EGLN2, etc.) and coat color (KIT, MITF, etc.). These results provide valuable full genome variation resources for Chinese bovine genome research and would be helpful for cattle breeding and selection programs in the future.

A genome-wide single nucleotide polymorphism and copy number variation analysis for number of piglets born alive

Background

In this study we integrated the CNV (copy number variation) and WssGWAS (weighted single-step approach for genome-wide association) analyses to increase the knowledge about number of piglets born alive, an economically important reproductive trait with significant impact on production efficiency of pigs.

Results

A total of 3892 samples were genotyped with the Porcine SNP80 BeadChip. After quality control, a total of 57,962 high-quality SNPs from 3520 Duroc pigs were retained. The PennCNV algorithm identified 46,118 CNVs, which were aggregated by overlapping in 425 CNV regions (CNVRs) ranging from 2.5 Kb to 9718.4 Kb and covering 197 Mb (~ 7.01%) of the pig autosomal genome. The WssGWAS identified 16 genomic regions explaining more than 1% of the additive genetic variance for number of piglets born alive. The overlap between CNVR and WssGWAS analyses identified common regions on SSC2 (4.2–5.2 Mb), SSC3 (3.9–4.9 Mb), SSC12 (56.6–57.6 Mb), and SSC17 (17.3–18.3 Mb). Those regions are known for harboring important causative variants for pig reproductive traits based on their crucial functions in fertilization, development of gametes and embryos. Functional analysis by the Panther software identified 13 gene ontology biological processes significantly represented in this study such as reproduction, developmental process, cellular component organization or biogenesis, and immune system process, which plays relevant roles in swine reproductive traits.

Conclusion

Our research helps to improve the understanding of the genetic architecture of number of piglets born alive, given that the combination of GWAS and CNV analyses allows for a more efficient identification of the genomic regions and biological processes associated with this trait in Duroc pigs. Pig breeding programs could potentially benefit from a more accurate discovery of important genomic regions.

Electronic supplementary material

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

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.

Genetic Diversity of Seven Cattle Breeds Inferred Using Copy Number Variations

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

Integrating CNVs into meta-QTL identified GBP4 as positional candidate for adult cattle stature

Copy number variation (CNV) of DNA sequences, functionally significant but yet fully ascertained, is believed to confer considerable increments in unexplained heritability of quantitative traits. Identification of phenotype-associated CNVs (paCNVs) therefore is a pressing need in CNV studies to speed up their exploitation in cattle breeding programs. Here, we provided a new avenue to achieve this goal that is to project the published CNV data onto meta-quantitative trait loci (meta-QTL) map which connects causal genes with phenotypes. Any CNVs overlapping meta-QTL therefore will be potential paCNVs. This study reported potential paCNVs in Bos taurus autosome 3 (BTA3). Notably, overview indexes and CNVs both highlighted a narrower region (BTA3 54,500,000-55,000,000 bp, named BTA3_INQTL_6) within one constructed meta-QTL. Then, we ascertained guanylate-binding protein 4 (GBP4) among the nine positional candidate genes was significantly associated with adult cattle stature, including body weight (BW, P < 0.05) and withers height (WHT, P < 0.05), fitting GBP4 CNV either with three levels or with six levels in the model. Although higher copy number downregulated the mRNA levels of GBP2 (P < 0.05) and GBP4 (P < 0.05) in 1-Mb window (54.0-55.0 Mb) in muscle and adipose, additional analyses will be needed to clarify the causality behind the ascertained association.

Detection of breed-specific copy number variations in domestic chicken genome

Copy number variations (CNVs) are important large-scale variants. They are widespread in the genome and may contribute to phenotypic variation. Detection and characterization of CNVs can provide new insights into the genetic basis of important traits. Here, we perform whole-genome short read sequence analysis to identify CNVs in two indigenous and commercial chicken breeds to evaluate the impact of the identified CNVs on breed-specific traits. After filtration, a total of 12 955 CNVs spanning (on average) about 9.42% of the chicken genome were found that made up 5467 CNV regions (CNVRs). Chicken quantitative trait loci (QTL) datasets and Ensembl gene annotations were used as resources for the estimation of potential phenotypic effects of our CNVRs on breed-specific traits. In total, 34% of our detected CNVRs were also detected in earlier CNV studies. These CNVRs partly overlap several previously reported QTL and gene ontology terms associated with some important traits, including shank length QTL in Creeper-specific CNVRs and body weight and egg production characteristics, as well as muscle and body organ growth, in the Arian commercial breed. Our findings provide new insights into the genomic structure of the chicken genome for an improved understanding of the potential roles of CNVRs in differentiating between breeds or lines.

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves-an asthmalike inflammatory disease in horses. Analysis of 16 RAO-susceptible (cases) and six RAO-resistant (control) horses on a custom-made whole-genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO-affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.

Genomic population structure and prevalence of copy number variations in South African Nguni cattle

Background

Copy number variations (CNVs) are modifications in DNA structure comprising of deletions, duplications, insertions and complex multi-site variants. Although CNVs are proven to be involved in a variety of phenotypic discrepancies, the full extent and consequence of CNVs is yet to be understood. To date, no such genomic characterization has been performed in indigenous South African Nguni cattle. Nguni cattle are recognized for their ability to sustain harsh environmental conditions while exhibiting enhanced resistance to disease and parasites and are thought to comprise of up to nine different ecotypes.

Methods

Illumina BovineSNP50 Beadchip data was utilized to investigate genomic population structure and the prevalence of CNVs in 492 South African Nguni cattle. PLINK, ADMIXTURE, R, gPLINK and Haploview software was utilized for quality control, population structure and haplotype block determination. PennCNV hidden Markov model identified CNVs and genes contained within and 10 Mb downstream from reported CNVs. PANTHER and Ensembl databases were subsequently utilized for gene annotation analyses.

Results

Population structure analyses on Nguni cattle revealed 5 sub-populations with a possible sub-structure evident at K equal to 8. Four hundred and thirty three CNVs that formed 334 CNVRs ranging from 30 kb to 1 Mb in size are reported. Only 231 of the 492 animals demonstrated CNVRs. Two hundred and eighty nine genes were observed within CNVRs identified. Of these 149, 28, 44, 2 and 14 genes were unique to sub-populations A, B, C, D and E respectively. Gene ontology analyses demonstrated a number of pathways to be represented by respective genes, including immune response, response to abiotic stress and biological regulation processess.

Conclusions

CNVs may explain part of the phenotypic diversity and the enhanced adaptation evident in Nguni cattle. Genes involved in a number of cellular components, biological processes and molecular functions are reported within CNVRs identified. The significance of such CNVRs and the possible effect thereof needs to be ascertained and may hold interesting insight into the functional and adaptive consequence of CNVs in cattle.

Electronic supplementary material

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

Genetic Susceptibility to Rhodococcus equi

Rhodococcus equi pneumonia is a major cause of morbidity and mortality in neonatal foals. Much effort has been made to identify preventative measures and new treatments for R. equi with limited success. With a growing focus in the medical community on understanding the genetic basis of disease susceptibility, investigators have begun to evaluate the interaction of the genetics of the foal with R. equi. This review describes past efforts to understand the genetic basis underlying R. equi susceptibility and tolerance. It also highlights the genetic technology available to study horses and describes the use of this technology in investigating R. equi. This review provides readers with a foundational understanding of candidate gene approaches, single nucleotide polymorphism‐based, and copy number variant‐based genome‐wide association studies, and next generation sequencing (both DNA and RNA).

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

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

Genome-wide patterns of copy number variation in the diversified chicken genomes using next-generation sequencing

BACKGROUND

Copy number variation (CNV) is important and widespread in the genome, and is a major cause of disease and phenotypic diversity. Herein, we performed a genome-wide CNV analysis in 12 diversified chicken genomes based on whole genome sequencing.

RESULTS

A total of 8,840 CNV regions (CNVRs) covering 98.2 Mb and representing 9.4% of the chicken genome were identified, ranging in size from 1.1 to 268.8 kb with an average of 11.1 kb. Sequencing-based predictions were confirmed at a high validation rate by two independent approaches, including array comparative genomic hybridization (aCGH) and quantitative PCR (qPCR). The Pearson's correlation coefficients between sequencing and aCGH results ranged from 0.435 to 0.755, and qPCR experiments revealed a positive validation rate of 91.71% and a false negative rate of 22.43%. In total, 2,214 (25.0%) predicted CNVRs span 2,216 (36.4%) RefSeq genes associated with specific biological functions. Besides two previously reported copy number variable genes EDN3 and PRLR, we also found some promising genes with potential in phenotypic variation. Two genes, FZD6 and LIMS1, related to disease susceptibility/resistance are covered by CNVRs. The highly duplicated SOCS2 may lead to higher bone mineral density. Entire or partial duplication of some genes like POPDC3 may have great economic importance in poultry breeding.

CONCLUSIONS

Our results based on extensive genetic diversity provide a more refined chicken CNV map and genome-wide gene copy number estimates, and warrant future CNV association studies for important traits in chickens.

Copy number variation in the horse genome

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

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

Detection of genome-wide copy number variations in two chicken lines divergently selected for abdominal fat content

Background

The chicken (Gallus gallus) is an important model organism that bridges the evolutionary gap between mammals and other vertebrates. Copy number variations (CNVs) are a form of genomic structural variation widely distributed in the genome. CNV analysis has recently gained greater attention and momentum, as the identification of CNVs can contribute to a better understanding of traits important to both humans and other animals. To detect chicken CNVs, we genotyped 475 animals derived from two broiler chicken lines divergently selected for abdominal fat content using chicken 60 K SNP array, which is a high-throughput method widely used in chicken genomics studies.

Results

Using PennCNV algorithm, we detected 438 and 291 CNVs in the lean and fat lines, respectively, corresponding to 271 and 188 CNV regions (CNVRs), which were obtained by merging overlapping CNVs. Out of these CNVRs, 99% were confirmed also by the CNVPartition program. These CNVRs covered 40.26 and 30.60 Mb of the chicken genome in the lean and fat lines, respectively. Moreover, CNVRs included 176 loss, 68 gain and 27 both (i.e. loss and gain within the same region) events in the lean line, and 143 loss, 25 gain and 20 both events in the fat line. Ten CNVRs were chosen for the validation experiment using qPCR method, and all of them were confirmed in at least one qPCR assay. We found a total of 886 genes located within these CNVRs, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed they could play various roles in a number of biological processes. Integrating the results of CNVRs, known quantitative trait loci (QTL) and selective sweeps for abdominal fat content suggested that some genes (including SLC9A3, GNAL, SPOCK3, ANXA10, HELIOS, MYLK, CCDC14, SPAG9, SOX5, VSNL1, SMC6, GEN1, MSGN1 and ZPAX) may be important for abdominal fat deposition in the chicken.

Conclusions

Our study provided a genome-wide CNVR map of the chicken genome, thereby contributing to our understanding of genomic structural variations and their potential roles in abdominal fat content in the chicken.

Electronic supplementary material

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

Enhancing genome-wide copy number variation identification by high density array CGH using diverse resources of pig breeds

Copy number variations (CNVs) are important forms of genomic variation, and have attracted extensive attentions in humans as well as domestic animals. In the study, using a custom-designed 2.1 M array comparative genomic hybridization (aCGH), genome-wide CNVs were identified among 12 individuals from diverse pig breeds, including one Asian wild population, six Chinese indigenous breeds and two modern commercial breeds (Yorkshire and Landrace), with one individual of the other modern commercial breed, Duroc, as the reference. A total of 1,344 CNV regions (CNVRs) were identified, covering 47.79 Mb (∼1.70%) of the pig genome. The length of these CNVRs ranged from 3.37 Kb to 1,319.0 Kb with a mean of 35.56 Kb and a median of 11.11 Kb. Compared with similar studies reported, most of the CNVRs (74.18%) were firstly identified in present study. In order to confirm these CNVRs, 21 CNVRs were randomly chosen to be validated by quantitative real time PCR (qPCR) and a high rate (85.71%) of confirmation was obtained. Functional annotation of CNVRs suggested that the identified CNVRs have important function, and may play an important role in phenotypic and production traits difference among various breeds. Our results are essential complementary to the CNV map in the pig genome, which will provide abundant genetic markers to investigate association studies between various phenotypes and CNVs in pigs.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

Sequencing and annotated analysis of the Holstein cow genome

The aim of our study was to create a high-quality Holstein cow genome reference sequence and describe the different types of variations in this genome compared to the reference Hereford breed. We generated one fragment and three mate-paired libraries from genomic DNA. Raw files were mapped and paired to the reference cow (Bos taurus) genome assemblies bosTau6/UMD_3.1. BioScope (v1.3) software was used for mapping and variant analysis. Initial sequencing resulted in 2,842,744,008 of 50-bp reads. Average mapping efficiency was 78.4 % and altogether 2,168,425,497 reads and 98,022,357,422 bp were successfully mapped, resulting in 36.7X coverage. Tertiary analysis found 5,923,230 SNPs in the bovine genome, of which 3,833,249 were heterozygous and 2,089,981 were homozygous variants. Annotation revealed that 4,241,000 of all discovered SNPs were annotated in the dbSNP database and 1,682,230 SNPs were considered as novel. Large indel variations accounted for 48,537,190 bp of the entire genome and there were 138,504 of them. The largest deletion was 18,594 bp and the largest insertion was 13,498 bp. Another group of variants, small indels (n = 458,061), accounted for the total variation of 1,839,872 nucleotides in the genome. Only 92,115 small indels were listed in the dbSNP and therefore 365,946 small indels were novel. Finally, we identified 1,876 inversions in the bovine genome. In conclusion, this is another description of the Holstein cow genome and, similar to previous studies, we found a large amount of novel variations. Better knowledge of these variations could explain significant phenotypic differences (e.g., health, production, reproduction) between different breeds.

Identification of genome-wide copy number variations among diverse pig breeds using SNP genotyping arrays

Copy number variations (CNVs) are important forms of genetic variation complementary to SNPs, and can be considered as promising markers for some phenotypic and economically important traits or diseases susceptibility in domestic animals. In the present study, we performed a genome-wide CNV identification in 14 individuals selected from diverse populations, including six types of Chinese indigenous breeds, one Asian wild boar population, as well as three modern commercial foreign breeds. We identified 63 CNVRs in total, which covered 9.98 Mb of polymorphic sequence and corresponded to 0.36% of the genome sequence. The length of these CNVRs ranged from 3.20 to 827.21 kb, with an average of 158.37 kb and a median of 97.85 kb. Functional annotation revealed these identified CNVR have important molecular function, and may play an important role in exploring the genetic basis of phenotypic variability and disease susceptibility among pigs. Additionally, to confirm these potential CNVRs, we performed qPCR for 12 randomly selected CNVRs and 8 of them (66.67%) were confirmed successfully. CNVs detected in diverse populations herein are essential complementary to the CNV map in the pig genome, which provide an important resource for studies of genomic variation and the association between various economically important traits and CNVs.