Association study and expression analysis of CYP4A11 gene copy number variation in Chinese cattle

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.

Distribution of DGAT1 copy number variation in Chinese goats and its associations with milk production traits

Goat is an important sector for meat and dairy products. Diacylglycerol O-acyltransferase 1 (DGAT1), which is a key gene in milk production, has been recently detected to overlap with a novel copy number variation (CNV) in goats. CNVs could be genetic markers providing new insights into the genetic basis of phenotypic variation. Up to now, there are no reports on the DGAT1-related CNV (DGAT1 CNV) in Chinese goats. This study first detected the distribution of the DGAT1 CNV in Chinese seven goat breeds, finding substantial differences among dairy, meat, and fiber goats (P < 0.01). The association analysis between the DGAT1 CNV and milk production traits revealed significant associations: Xinong Sannen (XS) dairy goat with copy number loss type had higher freezing point depression (FPD) (P < 0.01) and milk solids-not-fat (SNF) content (P < 0.05). Overall, our study unraveled the distribution of DGAT1 CNV in Chinese goats for the first time and found the potential role of this CNV in the marker-assisted selection of dairy goat breeding.

Copy Number Variants in Two Northernmost Cattle Breeds Are Related to Their Adaptive Phenotypes

Copy number variations (CNVs) are genomic structural variants with potential functional and evolutionary effects on phenotypes. In this study, we report the identification and characterization of CNVs from the whole-genome resequencing data of two northernmost cattle breeds from Russia: the Yakut and Kholmogory cattle and their phylogenetically most related breeds, Hanwoo and Holstein, respectively. Comparisons of the CNV regions (CNVRs) among the breeds led to the identification of breed-specific CNVRs shared by cold-adapted Kholmogory and Yakut cattle. An investigation of their overlap with genes, regulatory domains, conserved non-coding elements (CNEs), enhancers, and quantitative trait loci (QTLs) was performed to further explore breed-specific biology and adaptations. We found CNVRs enriched for gene ontology terms related to adaptation to environments in both the Kholmogory and Yakut breeds and related to thermoregulation specifically in Yakut cattle. Interestingly, the latter has also been supported when exploring the enrichment of breed-specific CNVRs in the regulatory domains and enhancers, CNEs, and QTLs implying the potential contribution of CNVR to the Yakut and Kholmogory cattle breeds’ adaptation to a harsh environment.

Detection of InDel and CNV of SPAG17 gene and their associations with bovine growth traits

Sperm-associated antigen 17 (SPAG17) gene encodes a central pair protein, which is involved in flagellar motility, male fertility and skeletal growth in ruminants. The insertions/deletions (indels) and copy number variations (CNVs) influence phenotypic traits by altering the sequences and copy numbers of functional genes, respectively. This study identified a novel 8-bp indel of SPAG17 gene in 1520 individuals from eight different cattle breeds, as well as a novel CNV region in 355 animals. The correlation analysis of indel showed that the individuals of ID genotype had superior performance traits such as body height (p = 0.038) and body slanting length (p = 0.041) as compared to other genotypes in Xianan cattle. For the CNV, different copy numbers were closely related to the body height in Qinchuan (p = 0.045) and body weight in Xianan (p = 0.036) breeds. Importantly, significant difference was observed between the 8-bp indel and the copy number loss in Xianan breed (p < 0.01). These findings indicated that the variations within the bovine SPAG17 gene can be considered as an effective DNA molecular marker for beef cattle breeding.

Two Different Copy Number Variations of the CLCN2 Gene in Chinese Cattle and Their Association with Growth Traits

Copy number variation (CNV) can affect gene function and even individual phenotypic traits by changing the transcription and translation level of related genes, and it also plays an important role in species evolution. Chloride voltage-gated channel 2 (CLCN2) encodes a voltage-gated chloride channel (CLC-2), which has a wide organ distribution and is ubiquitously expressed. Based on previous studies, we hypothesize that CLCN2 could be a candidate gene involved in cell volume regulation, transepithelial transport and cell proliferation. This study aimed to explore CNVs in the CLCN2 gene and investigate its association with growth traits in four Chinese cattle breeds (Yunling cattle, Xianan cattle, Qinchuan cattle and Pinan cattle). We identified there are two copy number variation regions (CNV1: 3600 bp, including exon 2-11; CNV2: 4800 bp, including exon 21-22) of the CLCN2 gene. The statistical analysis showed that the CNV1 mutation in the YL cattle population was significantly associated with cannon circumference (p < 0.01). The CNV2 mutation in the XN cattle population had a significant effect on body slanting length, chest girth and body weight (p < 0.05). In the YL cattle, the association analysis of CLCN2 gene CNV1 and CNV2 combination with cannon circumference was significant (p < 0.01). Our results provide evidence that CNV1 and CNV2 in CLCN2 are associated with growth traits in two different cattle populations and could be used as candidate markers for cattle molecular breeding.

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.

Copy number variation of the HPGDS gene in the Ashidan yak and its associations with growth traits

Copy number variation (CNV) is a structural variation at the submicroscopic level of the genome, which can affect gene-related phenotypes by changing genes dosage and transcript structure. Hematopoietic prostaglandin D synthase (HPGDS) is a member whose functions are closely related to weight gain and inflammatory diseases of the glutathione S-transferase (GSTs) family. In this study, the growth characteristics (body weight, withers height, body length, and chest girth) of 336 Ashidan yaks were monitored at four stages (6 months, 12 months, 18 months, and 30 months). In addition, CNV of the HPGDS gene was detected, discovered relationships of CNV with growth traits, and explored the level of gene expression. Based on the statistical analysis by IBM SPSS software, significant correlations were observed between HPGDS-CNV and body weight in 12-month-old yak (P < 0.01), 18-month-old yak (P < 0.001) and 30-month-old yak (P < 0.001) and body length in 18-month-old yak (P < 0.05) and 30-month-old yak (P < 0.05), respectively. Additionally, the individuals with gain copy number type performed better in body weight and body length than those with normal or loss copy number type. To our best of knowledge, this is the first time to make efforts to probe into the role of HPGDS-CNV and its interaction with livestock growth traits. Our results suggested that the CNV of the HPGDS gene may be an active candidate gene for the marker-assisted selection (MAS) of yaks.

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

Background

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

Results

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

Conclusions

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

Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle

Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.

Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak

Simple Summary

Cell adhesion molecule 2 (CADM2), also known as synaptic cell adhesion molecule 2 (SYNCAM2), is the mediator of synaptic signals enriched in the brain. Overlaps between copy number variation (CNV) regions in CADM2 and quantitative trait loci (QTL) related to body weight have been clarified in a previous study. In this study, two loci were amplified in the CADM2 gene (CNV1: 235,915 bp, exon 1 and partial intron 1; CNV2: 60,430 bp, intron 9) to explore the relationship between CNV types in the CADM2 gene and growth traits in 350 Ashidan yaks. Association analysis illustrated that no significant effect was found on growth traits in CNV1. However, the CNV2 mutation had a significant effect on body weight at the sixth month (p < 0.05). Individuals with the gain-type copy number variation CNV2 were significantly superior to those with loss- or normal-type in terms of body weight (p < 0.05). In summary, this study confirmed that CADM2-CNVs affect growth traits in yaks, and may be candidate genes for successful yak breeding and genetics projects.

Abstract

Copy number variation (CNV) is currently accepted as a common source of genetic variation. It is reported that CNVs may influence the resistance to disease and complex economic traits, such as residual feed intake, muscle formation, and fat deposition in livestock. Cell adhesion molecule 2 (CADM2) is expressed widely in the brain and adipose tissue and can regulate body weight through the central nervous system. Growth traits are important economic traits for animal selection. In this study, we aimed to explore the effect of CADM2 gene copy number variants on yak growth traits. Here, two CNVs in the CADM2 gene were investigated using the quantitative polymerase chain reaction (qPCR), and the association of the CNVs with growth traits in yak was analyzed using statistical methods by SPSS software. Differences were considered significant if the p value was < 0.05. Statistical analysis indicated significant association of CADM2-CNV2 with the body weight of the Chinese Ashidan yak. A significant effect of CNV2 (p < 0.05) was found on body weight at 6 months. In CNV2, the gain-type copy number variation exhibited greater performance than the other variants, with greater body weight observed at 6 months (p < 0.05). To the best of our knowledge, this is the first attempt to investigate the function of CADM2-CNVs and their association with growth traits in animals. This may be a useful candidate marker in marker-assisted selection of yaks.

Associations of ORMDL1 gene copy number variations with growth traits in four Chinese sheep breeds

Copy number variations (CNVs) are gains and losses of genomic sequence of more than 50 bp between two individuals of a species. Also, CNV is considered to be one of the main elements affecting the phenotypic diversity and evolutionary adaptation of animals. ORMDL sphingolipid biosynthesis regulator 1 (ORMDL1) is a protein-coding gene associated with diseases and development. In our study, the polymorphism of ORMDL1 gene copy numbers in four Chinese sheep breeds (abbreviated CK, HU, STH, and LTH) was detected. In addition, we analyzed the transcriptional expression level of ORMDL1 gene in different tissues of sheep and examined the association of ORMDL1 CNV with growth traits. The statistical analysis revealed that ORMDL1 CNV was remarkably correlated with body height, heart girth, and circumference of cannon bone in HU sheep (P<0.05), and there are significant effects on body weight, body height, body length, chest depth, and height of hip cross in STH sheep (P<0.05). In conclusion, our results provide a basis for the relationship between CNV of ORMDL1 gene and sheep growth traits, suggesting that ORMDL1 CNV may be considered a promising marker for the molecular breeding of Chinese sheep.

Copy number variation of bovine SHH gene is associated with body conformation traits in Chinese beef cattle

Sonic Hedgehog (Shh) regulates many key developmental processes during vertebrate limb development, fat formation, and skeletal tissue regeneration. Current whole genome sequencing data have identified a copy number variation mapping to bovine Sonic Hedgehog gene (SHH-CNV). The object of this study was to characterize the SHH-CNV distributions in 648 individuals from 11 Chinese cattle populations and further to investigate the associations of the copy number changes with gene expression and cattle growth traits. The SHH-CNV showed a high variance within Chinese indigenous yellow cattle. Compared to yak and dairy cattle, the beef cattle like Luxi and Xianan breed had significantly higher median copy numbers, suggesting the diversity of SHH-CNV in beef cattle selections. The negative correlation of SHH-CNV with SHH transcriptional level in adult adipose tissue (P < 0.01) indicated the dosage effects of SHH-CNV related to bovine fat formation. Association analysis of SHH-CNV and body size traits was conducted in five breeds. The results revealed that the copy number gain type of SHH-CNV exhibited significantly better chest depth in 24 months old Qinchuan cattle, and better body weight, body length, and chest girth in 18 months old Nanyang cattle, whereas the normal copy number had superior chest girth and body weight in adult Jinnan cattle (P < 0.05 or P < 0.01). In summary, this research uncovered meaningful effects of SHH-CNV on gene expression and cattle phenotypic traits, indicating its potential applications for genetic improvement of beef cattle.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins

BACKGROUND

Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.

RESULTS

We identified 1043 CNV regions (CNVRs) from array comparative genomic hybridization (CGH) data of 47 Holstein bulls. Using a probe-based CNV association approach, we detected 87 CNVRs significantly (Bonferroni-corrected P value < 0.05) associated with at least one out of 41 complex traits. Within them, 39 CNVRs were simultaneously associated with at least 2 complex traits. Notably, 24 CNVRs were markedly related to daughter pregnancy rate (DPR). For example, CNVR661 containing CYP4A11 and CNVR213 containing CTR9, respectively, were associated with DPR and other traits related to reproduction, production, and body conformation. CNVR758 was also significantly related to DPR, with a nearby gene CAPZA3, encoding one of F-actin-capping proteins which play a role in determining sperm architecture and male fertility. We corroborated these CNVRs by examining their overlapped quantitative trait loci and comparing with previously published CNV results.

CONCLUSION

To our knowledge, this is one of the first genome-wide association studies based on CNVs called by array CGH in Holstein cattle. Our results contribute substantial information about the potential CNV impacts on reproduction, health, production, and body conformation traits, which lay the foundation for incorporating CNV into the future dairy cattle breeding program.

Distribution and association study in copy number variation of KCNJ12 gene across four Chinese cattle populations

Copy number variation is a large genome variation which usually happens in the noncoding-region, and it may occur at the locus associated with the functional gene to further influence the phenotype. Potassium inwardly-rectifying channel, subfamily J 12 (KCNJ12) gene expressed widely in cardiomyocytes and neurons, plays an important role in tumor therapy and muscle movement regulation. In this study, we detected the distribution of CNVs for KCNJ12 gene in 404 individuals belonging to four Chinese cattle breeds (NY, JX, JA and GF). We also investigated the KCNJ12 gene expression in different tissues of JX cattle. Additionally, we examined the association of two CNV regions (CNV1: 1,600 bp, intron 1; CNV2: 4,800 bp, intergenic) with growth traits. The statistical analyses indicated that the CNV1 is associated with the body length, rump length and weight in JX cattle population (P < 0.05); and there has a significant association with the body length, chest circumference, and body weight in GF cattle (P < 0.05).The CNV2 had a significant effect on the body length and body weight in JX cattle (P < 0.05); the body length, chest circumference, rump length and body weight in GF cattle (P < 0.01 or P < 0.05). The copy numbers of KCNJ12 gene presented the negative correlations with the transcript level of gene in skeletal muscles (P < 0.05). Our results provide evidence that CNV1 and CNV 2 in KCNJ12 are associated with growth traits in two cattle populations and may be used as candidates for marker-assisted selection and breeding management in cattle.

Copy Number Variations of KLF6 Modulate Gene Transcription and Growth Traits in Chinese Datong Yak (Bos Grunniens)

Copy number variation (CNV) is a significant marker of the genetic and phenotypic diversity among individuals that accounts for complex quantitative traits of phenotype and diseases via modulating gene dosage and disrupting coding regions in the genome. Biochemically, Kruppel-like factor 6 (KLF6) genes plays a significant role in the regulation of cell differentiation and proliferation and muscle development. The aim of this study was to detect the distributions of KLF6 copy number variations (CNVs) in five breeds of domestic yak and to explore their effect on growth traits and gene expression. The data were analyzed by real-time quantitative PCR (qPCR). Our results elucidated that a decreased CNV in the KLF6 gene is more highly associated (p < 0.05) with various growth traits than increased or normal CNVs in six-month-old and five-year-old Datong yak. Nevertheless, negative correlations between the DNA copy number and KLF6 gene expression were observed in the skeletal muscle of adult Datong yak. These results suggest that CNVs of the KLF6 gene could be crucial genomic markers for growth phenotypes of Chinese Datong yak breeds and this finding constitutes the first evidence of the biological role of KLF6 CNVs in Chinese Datong yak breeds.

Copy number variations of MTHFSD gene across pig breeds and its association with litter size traits in Chinese indigenous Xiang pig

Copy number variation (CNV) is a major proportion of genetic variation, which changes the gene structure and dosage and affects gene expression and function. To validate the presence and the function of CNV in pig, we used real-time quantitative polymerase chain reaction (qPCR) method to validate a 496 kb CNV region comprising MTHFSD gene on chromosome 6 of Xiang pig detected by single nucleotide polymorphism (SNP) array. Then we investigated the distribution of the MTHFSD CNV in a total of 545 pigs in four breeds. About 46.2% and 32.7% individuals in the four pig breeds were detected to be types of loss and gain of MTHFSD locus. The relative copy numbers of MTHFSD gene showed the largest variation range (0-55 copies) in the Xiang pig population. The copy numbers of MTHFSD gene presented the positive correlations with the transcript level of MTHFSD gene in adult ovaries. Statistical analysis indicated that CNVs of MTHFSD gene was significantly changed the litter size traits of Xiang pigs, and the individuals with CNV gain showed more litter size than the CNV loss pigs. We have reasons to believe that the MTHFSD as RNA-binding protein play an important role in pig reproduction as a result of regulating MTHFS mRNA metabolism.