Next Article in Journal
Effect of Dextrose Supplementation in the Pre-Ovulatory Sow Diet to Reduce Seasonal Influences on Litter Birth Weight Variation
Previous Article in Journal
The Relationship between Uterine, Fecal, Bedding, and Airborne Dust Microbiota from Dairy Cows and Their Environment: A Pilot Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
*
Authors to whom correspondence should be addressed.
Animals 2019, 9(12), 1008; https://doi.org/10.3390/ani9121008
Submission received: 21 October 2019 / Revised: 16 November 2019 / Accepted: 18 November 2019 / Published: 21 November 2019
(This article belongs to the Section Animal Genetics and Genomics)

Abstract

:

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.

1. Introduction

With large duplications and deletions in gene copy numbers detected in humans and other primates, it is not unreasonable to expect that copy number variation could be a significant source of genetic variation in other animals [1]. In previous studies, many molecular markers have been assessed for improving economic traits in animal breeding, including restriction fragment length polymorphisms (RFLP), single nucleotide polymorphisms (SNP), indels (Insertion/Deletion) and structural variations, such as inversions, duplications and deletions. To date, SNPs are used extensively for complex economic disease traits in livestock [2,3,4,5]. More than 99% of indels are less than 50 bp in length [6,7], and some long insertions and deletions are returned into the concept of structural variation [8]. As a genetic marker, indels have great potential whether in growth and reproduction traits or the treatment of disease [9,10,11]. However, CNVs include duplication and/or deletions of DNA fragments ≥50 kb, and can cause a greater number of differences than SNPs and indels [12], which can be discovered by clone-based array comparative genome hybridization (array-CGH) [13]. Copy number variants can be engendered through four mechanisms: non-allelic homologous recombination (NAHR), non-homologous end joining (NHEJ), fork stalling and template switching (FoSTeS), and L1-mediated retro-transposition [14,15]. This submicroscopic structural variation comprises insertions, inversions and translocations in the genome [13]. Compared with SNPs, CNVs have a higher rate of mutation, with a number of diseases and congenital deficiencies caused by the action or the process of mutation in CNVs. Increasing evidence indicates that CNVs are associated with a number of complex human diseases, for example, the susceptibility to HIV-1 [16], autism spectrum disorders [17] and Tourettes Syndrome [18]. Conversely, many studies in livestock have revealed that CNVs are correlated with some important traits, such as parasite resistance in Angus cattle [19], growth traits in Chinese cattle [20,21], meat production and quality in pig [22] and feed conversion ratios in beef cattle [23]. These results demonstrate that CNVs may have a great influence on variation in phenotype, probably because of the difference in the copy number of specific sequences.
The yak (Bos grunniens) is critical to populations living on the Qinghai-Tibet Plateau and adjacent regions, as they provide materials indispensable for survival such as meat, milk, transportation and fur [24]. The yak’s ability to survive in harsh conditions makes it an ideal model organism to explore the mechanism of adaptation to high elevations in both animals and human beings. The Ashidan yak, a new yak breed bred by Chinese scientists, has typical characteristics such as lack of horns and a mild temperament. It is easier to maintain and feed in the stalls, and thus suited for large-scale intensive breeding in the cold and arid alpine areas of China. In recent years, several functional CNV regions have been discovered and validated in the yak, utilizing an NGS-based method [25,26,27]. In addition to this, in a previous association study, CNVRs throughout the yak genome have been constructed using SNP arrays which indicate that the cell adhesion molecule 2 (CADM2, also known as SynCAM2, Igsf4d and Nectin-like molecule 3) gene overlaps with a portion of CNVR29 and CNVR30 [28]. Variants of the CADM2 gene have been previously identified to play a pivotal role in the human body mass index (BMI) values by means of the central nervous system [29,30,31]. Furthermore, analysis in mice has found that CADM2 is strongly linked to body weight and energy homeostasis via brain activity [32].
In addition, a previous study has demonstrated that the CADM2 gene was located in quantitative trait loci (QTLs) related to body weight both at weaning and in yearlings in cattle [33]. Together, this evidence indicates that CADM2-CNVs in yaks may play a vital role in growth traits that are useful in artificial selection. Nevertheless, research of CADM2-CNV has lagged far behind, in humans and other animals. We hypothesize that copy number variants in the CADM2 gene may affect yak phenotypic traits. The present study investigated the interrelations of CNVs including CNV1 (235,915 bp, exon 1 and partial intron 1) and CNV2 (60,430 bp, intron 9) with the CADM2 gene with growth traits in the yak, by performing quantitative real-time PCR (qPCR).

2. Materials and Methods

2.1. Ethics Statement

All experiments in the current research were approved by the Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences (No. LIHPS-CAAS-2017-115). Blood sample collection and body measurements were performed in strict accordance with the guide for the Care and Use of Laboratory Animals.

2.2. Sample Collection and Growth Traits

Three hundred and fifty unrelated female Ashidan yaks that had been raised with the same levels of nutrition were randomly selected from the Datong yak farm (Qinghai Province, China). None of the animals involved in this study were unhealthy. Phenotypic data (i.e., body weight, withers height, body length, chest girth) were measured in the selected Ashidan yaks at three stages of growth (i.e., at 6 months, 1 year and 18 months of age) for subsequent analysis. Phenotypic measurements were strictly in line with standard methods proposed by Gilbert et al. [34].

2.3. Genomic DNA Extraction and Copy Number Identification

Genomic DNA was extracted using an EasyPure Blood Genomic DNA kit (TransGen Biotech Co., Ltd., Beijing, China) and its concentration and quality measured using a Nanodrop 2000 spectrophotometer (ThermoFisher Scientific Inc., Waltham, MA, USA). Quantified DNA samples were stored at −20 °C.
The copy number variation regions CADM2-CNV1 and CADM2-CNV2 are located at 33,143,512–33,379,427 bp and 33,569,316–33,629,746 bp from the reference genome sequence NC_037328.1, respectively (Figure 1).

2.4. Primer Design

The National Center for Biotechnology Information (NCBI) primer-BLAST web tool was utilized to design the primers for the analysis of CADM2-CNVs (Table 1). A polymerase chain reaction (PCR) was utilized to verify the quality and optimum temperature of primers. Products were analyzed on 1% agarose gel electrophoresis.

2.5. Quantitative PCR Analysis

A LightCyler® 96 quantitative polymerase chain reaction (qPCR) system (Roche, Basel, Switzerland) with SYBR Green dye was used in this study because of its convenience and accuracy. The bovine basic transcription factor 3 (BTF3) gene was chosen as the diploid internal reference gene [35]. The total volume of the reaction mixture was 20 μL, comprising 40 ng genomic DNA with 10 μL of SYBR Premix Ex Taq II (TaKaRa, Dalian, China), 2 μL of amplification primers and 7 μL of ddH2O. The q-PCR procedure involved a 30 s preincubation at 95 °C then a two-step amplification with 45 cycles of 95 °C for 5 s followed by 62 °C for 30 s. After completion of the cycling, the reaction was terminated by incubation at 95 °C for 5 s then 65 °C for 60 s then 95 °C continuously. All experiments were conducted in triplicates, and the results represented as mean values ± standard deviation (SD).

2.6. Statistical Analysis

Copy number values of the yak CADM2 gene were determined using the formula: 2 × 2 Δ Δ c t [36]. Cycle threshold (Ct) values for the target gene were normalized using the formula: Δ C t = C t t a r g e t C t r e f e r e n c e , with Δ C t of the test samples normalized using the formula: Δ Δ C t = Δ C t t e s t Δ C t c o n t r o l . As the precondition of the analysis of variance, the homogeneity of variance test and the normality test of each trait in different copy number variation (CNV) regions were done by Bartlett test and Shapiro-Wilk test before association analysis. One-way analysis of variance (ANOVA) and a nonparametric test were used to analyze the association of CADM2 gene copy number variant types with growth traits in the Ashidan yak. The following model was used to analyze the CNV effects on growth traits: Y i j = μ + C N V i j + e i j , where i represents the ith region, j the jth CNV type, Y i j the observed growth trait, μ the overall mean of each trait, with C N V i j representing the effect of ith CNV region and jth CNV type of the CADM2 gene. e i j is a random error. Assessment of differences between means was performed by one-way ANOVA with a post hoc least significant digit (LSD) multiple comparison test. Differences in the distribution of CNV types for different CNV regions were analyzed using a chi-square test in terms of the frequency of CNV types [37]. Based on the results of association analysis, a Kendall Tau correlation test was carried out. Other than these, the combination and interaction effect of CNV1 and CNV2 also detected by two-way analysis of variance (ANOVA).
All the statistical analyses were carried out using IBM SPSS Statistics software (Version 23, IBM, New York, NY, USA). Significance level was set up at p < 0.05 in this experiment.

3. Results

3.1. Quantitative PCR and the Distribution of CADM2-CNVs

Specific primers for the two CNVs are displayed in Table 1. Cycle threshold (Ct) values were used subsequently for quantitative analysis. The amplification curves and melt peaks (Figure 2) were used to determine primer specificity. In addition, the distribution frequency of different CNV types in CNV1 and CNV2 were measured (Figure 3 and Table 2).
The majority of individual animals had the loss-type variant in CNV1, the proportion exceeding 50%. However, in CNV2, the largest number of individuals had the normal variant, followed by the loss-type. The gain-type variant was similarly distributed within CNV1 and CNV2. The chi-square test indicated a significant difference in the distribution of CNV types between CNV1 and CNV2 (p < 0.01), suggesting that the CNV type is specific for different loci (Table 3).

3.2. Association Analysis

Association analysis indicated a significant correlation between copy number type and phenotypic measurements in the Ashidan yak population. Bartlett test showed the homogeneity of variance of each trait within each locus (Supplementary Table S1). And the Shapiro-Wilk test showed the normality of each trait in different CNV regions (Supplementary Table S2). One-way ANOVA indicated that the CNV2 mutation had a significant effect on body weight (p < 0.05) in 6th month yaks (Table 4 and Table 5). Multiple comparisons were conducted based on the ANOVA described above. Interestingly, the gain-type copy number variant of CADM2-CNV2 exhibited clear phenotypic superiority in comparison with the loss and normal types, with greater body weight values in 6th month Ashidan yaks (p < 0.05). Furthermore, our correlation test verified the significant correlation between CNV2 and body weight in 6th month (p < 0.01) (Table 6).
No significant difference was observed in the traits at 12 and 18 months. Association analysis of the CADM2 gene CNVs with growth traits is presented in Table 4 and Table 5.
Moreover, a two-way analysis of variance (ANOVA) was also carried out to test the joint effect of CADM2-CNV1 and CADM2-CNV2. The results revealed a significant association in 18th month withers height both in CNV1 (p < 0.01) and CNV2(p < 0.05) (Table 7). But no significant association was found in their combination.

4. Discussion

Researchers have previously given considerable attention to the analysis of molecular markers associated with economic traits in a number of species. For example, genetic variation in the ovine PROP1 gene was identified in 670 New Zealand Romney sheep, which were related to traits such as tailing weight and growth rate [38]. One 43-bp indel polymorphism in HS6ST3 was associated with shank length, chest depth and other growth traits in 860 chickens [39]. SNPs within the coding region of the KDM4D gene were explored and linked with testis morphology traits of male pigs [40]. However, copy number variance gives rise to more substantial mutations and more extensive genetic effects than other molecular markers [41]. It has been established that the total number of nucleotides affected by copy number variants is larger than the total number of SNPs [42]. Recently, CNV maps in cattle [35], sheep [43] and pigs [44] have been constructed, indicating that structural polymorphisms within them constitute an important proportion of their genomes. Furthermore, CNVs have been associated with various growth traits, such as body length in cattle [45] and pigs [46] and pigmentation traits in goats [47] and horses [48], in addition to reduced immunological resistance to diseases [49]. Consistent with these results, it is widely believed that this form of genetic variation in livestock is of great importance to phenotypic traits.
Quantitative polymerase chain reaction (q-PCR) is widely used, especially for the measurement of copy number variation because of its accuracy and simplicity. To investigate the probable effects of CADM2-CNVs in the yak growth process, we conducted an association study between 12 phenotypic traits and different CNV types using a general linear model for the first time. Significant association was observed between 6th body weight and CNV2. The statistics demonstrate that weaning weights are significantly correlated with CNV2 type. Interestingly, this is in agreement with the conclusion of previous studies [28]. This phenomenon may be due to the function of the CADM2 gene, which plays an important role in maintaining the synaptic circuitry of the central nervous system [50], and can mediate leptin sensitivity, thermogenesis and energy metabolism [51]. Notably, the gain-types in CNV2 caught our attention due to the outstanding body weight in 6th month animals. Here, no significant association between CNV1 and CNV2 and the body weight in 12 month-olds was observed, while in cattle a QTL harboring CADM2 gene was thought to be associated with the yearling weight. This was most likely cause by the following points. First of all, a quantitative trait like body weight is regulated by multiple genes with small effect. Though a 350 yak sample size is an appropriate sample to examine the variance, the marker’s effect is too insufficient to be detected. A larger population is needed to be studied in the future. Secondly, it is probably due to the yak’s special physiological regulation mechanism due to its particularity harsh living conditions. Unlike the cattle industry, yak has to struggle to survive in a cold season lasting for six months after weaning, which is extremely cold along with severe forage shortage.
A 6 months hay period they will endure after weaning, which is different from cattle. This phenomenon may lead to a unique physiological regulation of growth in yaks. Here, we did not find any associations between CNV types and phenotype at 12 and 18 months.
To explain why this CNV has such an uncommon influence on phenotype, firstly the CADM2 gene had been reported to overlap a QTL region that is correlated with body weight (weaning and yearling) [28]. Secondly, body weight and energy metabolism are maintained by the central nervous system and previous studies have identified that CADM2 is expressed more abundantly in different regions of the brain than in any other tissue [52]. Furthermore, plenty of studies have reported that CADM2 can significantly improve body weight and decrease leptin sensitivity via the brain and adipose tissue [32,53]. Moreover, genome-wide association studies in humans have demonstrated clear evidence that the CADM2 gene is strongly associated with the body mass index through the encoding of membrane proteins that mediate synaptic assembly [31,52,54]. Because of the importance of both the CADM2 gene and copy number variance in animals, we propose that CADM2-CNVs can be used as indispensable molecular markers of genetic improvement of the yak. To the best of our knowledge, copy number variance of CADM2 has not been previously investigated in livestocks such as cattle, sheep, pigs, etc. Therefore, we speculate that the CADM2 gene might play a crucial role in growth traits in the yak.
Growth traits are an economically important factor in animal husbandry that can improve production performance by selective breeding [55]. Consistent with this research, numerous studies have illustrated the association between gene copy number type and growth traits. Shi et al. demonstrated that the CNV of the LEPR gene exhibited significant correlation with growth traits in Chinese cattle [56]. Jiang et al. clarified that CNV of the sheep SHE gene was greatly linked with economic traits, such as body length and chest width [57]. Lin et al. suggested that the CNV of SOX6 is responsible for muscle development in chickens [58]. We propose that these molecular markers in various livestock offer the theoretical promise of animal selection.
The yak is a specialized breed famous for its adaptation to life on the plateau and hypoxia. Distinctive regulatory mechanisms require exploration, particularly their growth mechanism. In the present study, the two CNVs of the CADM2 gene in the Ashidan yak have been characterized and validated. The results demonstrate the effect of the CADM2 CNVs on phenotypic traits in different CNV types. From this evidence, we have established that CADM2 CNVs can act as molecular markers in Ashidan yaks for early selection.

5. Conclusions

In conclusion, this study described the distribution of CADM2 gene copy number variant types in 350 Ashidan yaks. CADM2-CNV2 is significantly associated with body weight in 6 months in the Ashidan yak. Further studies, such as the mechanism of physiological regulation with the CADM2 gene in yaks need to be carried out. Also a larger population is required to validate the association of this CADM2 gene with growth traits in the yak.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/9/12/1008/s1, Table S1: The homogeneity of variance test of each trait, Table S2: The test of normality of each trait in different CNV regions.

Author Contributions

Conceptualization, F.G., C.J. and P.Y.; methodology, F.G.; software, F.G. and C.J.; validation, F.G., C.J. and P.Y.; formal analysis, F.G.; investigation, F.G.; resources, C.L., C.J. and M.C.; data curation, F.G. and C.J.; writing—original draft preparation, F.G.; writing—review and editing, C.J., C.L. and M.C.; visualization, P.Y.; supervision, P.Y.; project administration, P.Y.; funding acquisition, P.Y.

Funding

This research was funded by the Innovation Project of Chinese Academy of Agriculture Sciences (grant number CAAS-ASTIP-2014-LIHPS-01). The APC was funded by the Lanzhou Institute of Animal Husbandry and Pharmaceutical Science.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sebat, J.; Lakshmi, B.; Troge, J.; Alexander, J.; Young, J.; Lundin, P.; Månér, S.; Massa, H.; Walker, M.; Chi, M.; et al. Large-scale copy number polymorphism in the human genome. Science 2004, 305, 525–528. [Google Scholar] [CrossRef]
  2. Chung, E.R.; Kim, W.T. Association of SNP marker in IGF-I and MYF5 candidate genes with growth traits in Korean cattle. Asian Australas. J. Anim. Sci. 2005, 18, 1061–1065. [Google Scholar] [CrossRef]
  3. Wang, A.; Zhang, Y.; Li, M.; Lan, X.; Wang, J.; Chen, H. SNP identification in FBXO32 gene and their associations with growth traits in cattle. Gene 2013, 515, 181–186. [Google Scholar] [CrossRef]
  4. Wu, J.; Qiao, L.; Liu, J.; Yuan, Y.; Liu, W. SNP variation in ADRB3 gene reflects the breed difference of sheep populations. Mol. Biol. Rep. 2012, 39, 8395–8403. [Google Scholar] [CrossRef]
  5. Wang, X.; Yang, Q.; Zhang, S.; Zhang, X.; Pan, C.; Chen, H.; Zhu, H.; Lan, X. Genetic effects of single nucleotide polymorphisms in the goat GDF9 gene on prolificacy: True or false positive? Animals 2019, 9, 886. [Google Scholar] [CrossRef]
  6. Britten, R.J.; Rowen, L.; Williams, J.; Cameron, R.A. Majority of divergence between closely related DNA samples is due to indels. Proc. Natl. Acad. Sci. USA 2003, 100, 4661–4665. [Google Scholar] [CrossRef]
  7. Alkan, C.; Coe, B.P.; Eichler, E.E. Genome structural variation discovery and genotyping. Nat. Rev. Genet. 2011, 12, 363–376. [Google Scholar] [CrossRef]
  8. Tuzun, E.; Sharp, A.J.; Bailey, J.A.; Kaul, R.; Morrison, V.A.; Pertz, L.M.; Haugen, E.; Hayden, H.; Albertson, D.; Pinkel, D. Fine-scale structural variation of the human genome. Nat. Genatics 2005, 37, 727–732. [Google Scholar] [CrossRef]
  9. Han, S.H.; Shin, K.Y.; Lee, S.S.; Ko, M.S.; Cho, I.C. SINE indel polymorphism of AGL gene association with growth and carcass traits in Landrace X Jeju black pig F2 population. Mol. Biol. Rep. 2010, 37, 467–471. [Google Scholar] [CrossRef]
  10. Zhang, W.; Fang, M.; Li, Y.; Nie, Q.; Zhang, X. Identification of TDRP1 gene and its association with pig reproduction traits. DNA Cell Biol. 2012, 31, 371–377. [Google Scholar] [CrossRef]
  11. Dong-Mei, X.I.; Liu, Q.; Hong-Man, Y.U.; Yang, Y.A.; Mao, H.M.; Deng, W.D. Research progress in the relationship of polymorphism of prion protein gene (PRNP) with disease resistance for sheep and cattle. J. Yunnan Agric. Univ. 2011, 26, 418–425. [Google Scholar]
  12. Iafrate, A.J.; Feuk, L.; Rivera, M.N.; Listewnik, M.L.; Donahoe, P.K.; Qi, Y.; Scherer, S.W.; Lee, C. Detection of large-scale variation in the human genome. Nat Genet 2004, 36, 949–951. [Google Scholar] [CrossRef]
  13. Lars, F.; Carson, A.R.; Scherer, S.W. Structural variation in the human genome. Nat. Rev. Genet. 2006, 7, 85–97. [Google Scholar]
  14. Kidd, J.M.; Cooper, G.M.; Donahue, W.F.; Hayden, H.S.; Nick, S.; Tina, G.; Nancy, H.; Brian, T.; Can, A.; Francesca, A. Mapping and sequencing of structural variation from eight human genomes. Nature 2008, 453, 56–64. [Google Scholar] [CrossRef]
  15. Gu, W.; Zhang, F.; Lupski, J.R. Mechanisms for human genomic rearrangements. Pathogenetics 2008, 1, 4. [Google Scholar] [CrossRef]
  16. Gonzalez, E.; Kulkarni, H.; Bolivar, H.; Mangano, A.; Sanchez, R.; Catano, G.; Nibbs, R.J.; Freedman, B.I.; Quinones, M.P.; Bamshad, M.J. The influence of CCL3L1 gene-containing segmental duplications on HIV-1/AIDS susceptibility. Science 2005, 307, 1434–1440. [Google Scholar] [CrossRef]
  17. Dalila, P.; Pagnamenta, A.T.; Lambertus, K.; Richard, A.; Daniele, M.; Regina, R.; Judith, C.; Magalhaes, T.R.; Catarina, C.; Abrahams, B.S. Functional impact of global rare copy number variation in autism spectrum disorders. Nature 2010, 466, 368–372. [Google Scholar]
  18. Bassett, A.S.; Scherer, S.W. Copy number variation in tourette syndrome. Neuron 2017, 94, 1041. [Google Scholar] [CrossRef]
  19. Hou, Y.; Liu, G.E.; Bickhart, D.M.; Matukumalli, L.K.; Li, C.; Song, J.; Gasbarre, L.C.; Tassell, C.P.V.; Sonstegard, T.S. Genomic regions showing copy number variations associate with resistance or susceptibility to gastrointestinal nematodes in Angus cattle. Funct. Integr. Genom. 2012, 12, 81–92. [Google Scholar] [CrossRef]
  20. Zhang, G.M.; Zheng, L.; He, H.; Song, C.C.; Zhang, Z.J.; Cao, X.K.; Lei, C.Z.; Lan, X.Y.; Qi, X.L.; Chen, H. Associations of GBP2 gene copy number variations with growth traits and transcriptional expression in Chinese cattle. Gene 2018, 647, 101–106. [Google Scholar] [CrossRef]
  21. Yao, X.; Tao, S.; Hanfang, C.; Yang, Z.; Xianyong, L.; Chunlei, Z.; Chuzhao, L.; Xinglei, Q.; Hong, C. Associations of MYH3 gene copy number variations with transcriptional expression and growth traits in Chinese cattle. Gene 2014, 535, 106–111. [Google Scholar]
  22. Chen, C.; Qiao, R.; Wei, R.; Guo, Y.; Ai, H.; Ma, J.; Ren, J.; Huang, L. A comprehensive survey of copy number variation in 18 diverse pig populations and identification of candidate copy number variable genes associated with complex traits. BMC Genom. 2012, 13, 733. [Google Scholar] [CrossRef]
  23. Santana, M.H.D.A.; Junior, G.A.O.; Cesar, A.S.M.; Freua, M.C.; Gomes, R.D.C.; Silva, S.D.L.E.; Leme, P.R.; Fukumasu, H.; Carvalho, M.E.; Ventura, R.V. Copy number variations and genome-wide associations reveal putative genes and metabolic pathways involved with the feed conversion ratio in beef cattle. J. Appl. Genet. 2016, 57, 1–10. [Google Scholar]
  24. Wiener, G.; Han, J.L.; Long, R.J.; Wiener, G.; Han, J.L.; Long, R.J. The yak. RAP Publ. 2011, 44, 57–58. [Google Scholar]
  25. Zhang, X.; Wang, K.; Wang, L.; Yang, Y.; Ni, Z.; Xie, X.; Shao, X.; Han, J.; Wan, D.; Qiu, Q. Genome-wide patterns of copy number variation in the Chinese yak genome. BMC Genom. 2016, 17, 379. [Google Scholar] [CrossRef]
  26. Goshu, H.; Wu, X.; Chu, M.; Bao, P.; Ding, X.; Yan, P. Copy number variations of KLF6 modulate gene transcription and growth traits in Chinese Datong Yak (Bos Grunniens). Animals 2018, 8, 145. [Google Scholar] [CrossRef]
  27. Goshu, H.A.; Chu, M.; Xiaoyun, W.; Pengjia, B.; Zhi, D.X.; Yan, P. Genomic copy number variation of the CHKB gene alters gene expression and affects growth traits of Chinese domestic yak (Bos grunniens) breeds. Mol. Genet. Genom. 2019, 294, 549–561. [Google Scholar] [CrossRef]
  28. Jia, C.; Wang, H.; Li, C.; Wu, X.; Zan, L.; Ding, X.; Guo, X.; Bao, P.; Pei, J.; Chu, M. Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip. BMC Genom. 2019, 20, 376. [Google Scholar] [CrossRef]
  29. Speliotes, E.K.; Willer, C.J.; Berndt, S.I.; Monda, K.L.; Thorleifsson, G.; Jackson, A.U.; Allen, H.L.; Lindgren, C.M.; Luan, J.; Magi, R.; et al. Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat. Genet. 2010, 42, 937–953. [Google Scholar] [CrossRef]
  30. Speakman, J.R.; Loos, R.J.F.; O’Rahilly, S.; Hirschhorn, J.N.; Allison, D.B. GWAS for BMI: A treasure trove of fundamental insights into the genetic basis of obesity. Int. J. Obes. 2018, 42, 1524–1531. [Google Scholar] [CrossRef]
  31. Locke, A.E.; Bratati, K.; Berndt, S.I.; Justice, A.E.; Pers, T.H.; Day, F.R.; Corey, P.; Sailaja, V.; Buchkovich, M.L.; Jian, Y. Genetic studies of body mass index yield new insights for obesity biology. Nature 2015, 518, 197–206. [Google Scholar] [CrossRef] [PubMed]
  32. Yan, X.; Wang, Z.; Schmidt, V.; Gauert, A.; Willnow, T.E.; Heinig, M.; Poy, M.N. Cadm2 regulates body weight and energy homeostasis in mice. Mol. Metab. 2017, 8, 180–188. [Google Scholar] [CrossRef] [PubMed]
  33. Snelling, W.M.; Allan, M.F.; Keele, J.W.; Kuehn, L.A.; Mcdaneld, T.; Smith, T.P.L.; Sonstegard, T.S.; Thallman, R.M.; Bennett, G.L. Genome-wide association study of growth in crossbred beef cattle. J. Anim. Sci. 2010, 88, 837–848. [Google Scholar] [CrossRef]
  34. Gilbert, R.P.; Bailey, D.R.; Shannon, N.H. Linear body measurements of cattle before and after 20 years of selection for postweaning gain when fed two different diets. J. Anim. Sci. 1993, 71, 1712. [Google Scholar] [CrossRef]
  35. Bickhart, D.M.; Hou, Y.L.; Schroeder, S.G.; Alkan, C.; Cardone, M.F.; Matukumalli, L.K.; Song, J.Z.; Schnabe, R.D.; Ventura, M.; Taylor, J.F.; et al. Copy number variation of individual cattle genomes using next-generation sequencing. Genome Res. 2012, 22, 778–790. [Google Scholar] [CrossRef] [PubMed]
  36. Bae, J.S.; Cheong, H.S.; Kim, L.H.; Namgung, S.; Park, T.J.; Chun, J.Y.; Kim, J.Y.; Pasaje, C.F.A.; Jin, S.L.; Shin, H.D. Identification of copy number variations and common deletion polymorphisms in cattle. BMC Genom. 2010, 11, 232. [Google Scholar] [CrossRef]
  37. Yoav, B. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar]
  38. Ekegbu, U.J.; Burrows, L.; Amirpour-Najafabadi, H.; Zhou, H.; Hickford, J.G. Gene polymorphisms in PROP1 associated with growth traits in sheep. Gene 2019, 683, 41–46. [Google Scholar] [CrossRef]
  39. Wang, X.; Li, Z.; Guo, Y.; Wang, Y.; Sun, G.; Jiang, R.; Kang, X.; Han, R. Identification of a novel 43-bp insertion in the heparan sulfate 6-O-sulfotransferase 3 (HS6ST3) gene and its associations with growth and carcass traits in chickens. Anim. Biotechnol. 2019, 30, 252–259. [Google Scholar] [CrossRef]
  40. Zhou, T.; Wei, H.; Li, D.; Yang, W.; Cui, Y.; Gao, J.; Yu, T.; Lv, X.; Pan, C. A novel missense mutation within the domain of lysine demethylase 4D (KDM4D) gene is strongly associated with testis morphology traits in pigs. Anim. Biotechnol. 2019, 1–7. [Google Scholar] [CrossRef]
  41. Henrichsen, C.N.; Vinckenbosch, N.; Zöllner, S.; Chaignat, E.; Pradervand, S.; Schütz, F.; Ruedi, M.; Kaessmann, H.; Reymond, A. Segmental copy number variation shapes tissue transcriptomes. Nat. Genet. 2009, 41, 424. [Google Scholar] [CrossRef] [PubMed]
  42. Conrad, D.F.; Christine, B.; Ben, B.; Sarah, L.; Lira, M.; Charles, L.; Turner, D.J.; Hurles, M.E. Mutation spectrum revealed by breakpoint sequencing of human germline CNVs. Nat. Genet. 2010, 42, 385–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, J.; Zhang, L.; Xu, L.; Ren, H.; Lu, J.; Zhang, X.; Zhang, S.; Zhou, X.; Wei, C.; Zhao, F. Analysis of copy number variations in the sheep genome using 50K SNP BeadChip array. BMC Genom. 2013, 14, 229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Jiang, J.; Wang, J.; Wang, H.; Zhang, Y.; Kang, H.; Feng, X.; Wang, J.; Yin, Z.; Bao, W.; Zhang, Q. Global copy number analyses by next generation sequencing provide insight into pig genome variation. BMC Genom. 2014, 15, 593. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, M.; Lv, J.; Zhang, L.; Li, M.; Zhou, Y.; Lan, X.; Lei, C.; Chen, H. Association study and expression analysis of CYP4A11 gene copy number variation in Chinese cattle. Sci. Rep. 2017, 7, 46599. [Google Scholar] [CrossRef] [Green Version]
  46. Revilla, M.; Puig-Oliveras, A.; Castelló, A.; Crespo-Piazuelo, D.; Paludo, E.; Fernández, A.I.; Ballester, M.; Folch, J.M. A global analysis of CNVs in swine using whole genome sequence data and association analysis with fatty acid composition and growth traits. PLoS ONE 2017, 12, e0177014. [Google Scholar] [CrossRef]
  47. Zhang, R.-Q.; Wang, J.-J.; Zhang, T.; Zhai, H.-L.; Shen, W. Copy-number variation in goat genome sequence: A comparative analysis of the different litter size trait groups. Gene 2019, 696, 40–46. [Google Scholar] [CrossRef]
  48. Sundström, E.; Imsland, F.; Mikko, S.; Wade, C.; Sigurdsson, S.; Pielberg, G.R.; Golovko, A.; Curik, I.; Seltenhammer, M.H.; Sölkner, J. Copy number expansion of the STX17 duplication in melanoma tissue from Grey horses. BMC Genom. 2012, 13, 365. [Google Scholar] [CrossRef] [Green Version]
  49. Cantsilieris, S.; White, S.J. Correlating multiallelic copy number polymorphisms with disease susceptibility. Hum. Mutat. 2013, 34, 1–13. [Google Scholar] [CrossRef]
  50. Thomas, L.A.; Akins, M.R.; Thomas, B. Expression and adhesion profiles of SynCAM molecules indicate distinct neuronal functions. J. Comp. Neurol. 2010, 510, 47–67. [Google Scholar] [CrossRef] [Green Version]
  51. Rathjen, T.; Yan, X.; Kononenko, N.L.; Ku, M.C.; Song, K.; Ferrarese, L.; Tarallo, V.; Puchkov, D.; Kochlamazashvili, G.; Brachs, S. Regulation of body weight and energy homeostasis by neuronal cell adhesion molecule 1. Nat. Neurosci. 2017, 20, 1096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ibrahim-Verbaas, C.A.; Bressler, J.; Debette, S.; Schuur, M.; Smith, A.V.; Bis, J.C.; Davies, G.; Trompet, S.; Smith, J.A.; Wolf, C. GWAS for executive function and processing speed suggests involvement of the CADM2 gene. Mol. Psychiatry 2016, 21, 189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Morris, J.; Bailey, M.E.; Baldassarre, D.; Cullen, B.; de Faire, U.; Ferguson, A.; Gigante, B.; Giral, P.; Goel, A.; Graham, N. Genetic variation in CADM2 as a link between psychological traits and obesity. Sci. Rep. 2019, 9, 7339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lamiquiz-Moneo, I.; Mateo-Gallego, R.; Bea, A.M.; Dehesa-García, B.; Pérez-Calahorra, S.; Marco-Benedí, V.; Baila-Rueda, L.; Laclaustra, M.; Civeira, F.; Cenarro, A. Genetic predictors of weight loss in overweight and obese subjects. Sci. Rep. 2019, 9, 10770. [Google Scholar] [CrossRef] [Green Version]
  55. Rout, P.K.; Matika, O.; Kaushik, R.; Dass, G.; Singh, M.K.; Dige, M.S.; Bhusan, S. Genetic analysis of growth parameters and survival potential of Jamunapari goats in semiarid tropics. Small Rumin. Res. 2018, 165, 124–130. [Google Scholar] [CrossRef]
  56. Shi, T.; Xu, Y.; Yang, M.; Huang, Y.; Lan, X.; Lei, C.; Qi, X.; Yang, X.; Chen, H. Copy number variations at LEPR gene locus associated with gene expression and phenotypic traits in Chinese cattle. Anim. Sci. J. 2016, 87, 336–343. [Google Scholar] [CrossRef]
  57. Jiang, R.; Cheng, J.; Cao, X.-K.; Ma, Y.-L.; Chaogetu, B.; Huang, Y.-Z.; Lan, X.-Y.; Lei, C.-Z.; Hu, L.-Y.; Chen, H. Copy number variation of the SHE gene in sheep and its association with economic traits. Animals 2019, 9, 531. [Google Scholar] [CrossRef] [Green Version]
  58. Wright, D.; Boije, H.; Meadows, J.R.; Bed’Hom, B.; Gourichon, D.; Vieaud, A.; Tixier-Boichard, M.; Rubin, C.-J.; Imsland, F.; Hallböök, F. Copy number variation in intron 1 of SOX5 causes the Pea-comb phenotype in chickens. PLoS Genet. 2009, 5, e1000512. [Google Scholar] [CrossRef]
Figure 1. The two copy number variation regions of the CADM2 gene. The yellow boxes represent the coding regions and red boxes represent the CNVRs.
Figure 1. The two copy number variation regions of the CADM2 gene. The yellow boxes represent the coding regions and red boxes represent the CNVRs.
Animals 09 01008 g001
Figure 2. Amplification curves and melt peaks of the CADM2 and BTF3 genes.
Figure 2. Amplification curves and melt peaks of the CADM2 and BTF3 genes.
Animals 09 01008 g002
Figure 3. Distribution of different copy number variation (CNV) types in CNV1 and CNV2 in the Ashidan Yak.
Figure 3. Distribution of different copy number variation (CNV) types in CNV1 and CNV2 in the Ashidan Yak.
Animals 09 01008 g003
Table 1. Details of primers.
Table 1. Details of primers.
GenePrimer Pairs Sequence (5′-3′)Amplicon Length (bp)Tm (°C)
CADM2-P1F: GACTTCCCAGGATTGCCTGT
R: CCCTGGGAGCACAGTTGTTT
18662 °C
CADM2-P2F: GGCTGTCACGTTCTTCTCTCA
R: AGGGTTCATCCTGGAGGCTT
19662 °C
BTF3F: AACCAGGAGAAACTCGCCAA
R: TTCGGTGAAATGCCCTCTCG
16662 °C
F: forward primer; R: reverse primer.
Table 2. Distribution and frequency of CNV types in CNV1 and CNV2.
Table 2. Distribution and frequency of CNV types in CNV1 and CNV2.
TypeLossNormalGain
CNV1176 (0.50)161 (0.46)13 (0.04)
CNV2101 (0.29)215 (0.61)34 (0.1)
Figures represent numbers of animals (and proportion of the total population studied).
Table 3. Chi-square test of CNV type distribution in CNV1 and CNV2.
Table 3. Chi-square test of CNV type distribution in CNV1 and CNV2.
TypeCNV1CNV2
CNV1 37.445 (p = 0.000 **)
CNV2
Chi-square values ( χ 2 ) for differences between CNVs. ** Represent statistically significant differences at a level of p < 0.01.
Table 4. Statistical analysis of CADM2-CNV1 growth traits in the Ashidan yak.
Table 4. Statistical analysis of CADM2-CNV1 growth traits in the Ashidan yak.
AgeGrowth TraitCNV Type (Mean ± SE)p-Value
Loss (0.50)Normal (0.46)Gain (0.04)
6 monthsBody weight (kg)83.14 ± 0.773 a85.60 ± 0.833 b81.15 ± 2.207 a,b0.052
Withers height (cm)94.50 ± 0.37694.27 ± 0.43595.31 ± 1.4910.765
Body length (cm)91.46 ± 0.55792.53 ± 0.59490.54 ± 1.2280.325
Chest girth(cm)124.21 ± 0.567123.66 ± 0.662124.38 ± 1.6070.802
12 monthsBody weight (kg)82.17 ± 0.78483.14 ± 0.88180.64 ± 3.0580.584
Withers height (cm)90.62 ± 0.32490.28 ± 0.33590.31 ± 0.8940.755
Body length (cm)95.98 ± 0.38095.74 ± 0.38796.92 ± 1.7190.682
Chest girth(cm)117.35 ± 0.379116.93 ± 0.428116.92 ± 0.8430.747
18 monthsBody weight (kg)121.82 ± 1.119122.59 ± 1.258126.25 ± 3.3710.516
Withers height (cm)101.71 ± 0.452101.73 ± 0.476105.54 ± 3.0980.180
Body length (cm)101.22 ± 0.464102.07 ± 0.460101.69 ± 1.7840.439
Chest girth(cm)137.97 ± 0.840138.57 ± 0.795141.54 ± 3.5800.464
a,b Represent statistically significant differences at a level of p < 0.05.
Table 5. Statistical analysis of CADM2-CNV2 growth traits in the Ashidan yak.
Table 5. Statistical analysis of CADM2-CNV2 growth traits in the Ashidan yak.
AgeGrowth Trait CNV Type (Mean ± SE)p-Value
Loss (0.29)Normal (0.61)Gain (0.1)
6 monthsBody weight (kg)82.22 ± 1.004 a84.75 ± 0.723 b86.61 ± 1.559 c0.048 *
Withers height (cm)94.27 ± 0.54494.28 ± 0.35795.82 ± 0.7850.263
Body length (cm)90.93 ± 0.74392.30 ± 0.49992.74 ± 1.2700.278
Chest girth(cm)123.56 ± 0.877 a123.66 ± 0.518 a127.08 ± 1.013 b0.051
12 monthsBody weight (kg)82.53 ± 1.10782.43 ± 0.72684.53 ± 2.0130.579
Withers height (cm)90.09 ± 0.45690.54 ± 0.27890.94 ± 0.7290.534
Body length (cm)95.64 ± 0.49196.09 ± 0.33395.47 ± 1.0380.656
Chest girth(cm)116.74 ± 0.525117.24 ± 0.355117.65 ± 0.7780.606
18 monthsBody weight (kg)122.71 ± 1.510121.71 ± 1.076125.19 ± 1.9960.428
Withers height (cm)101.49 ± 0.640102.31 ± 0.414100.34 ± 1.2990.180
Body length (cm)101.11 ± 0.667101.79 ± 0.396102.13 ± 0.9130.567
Chest girth(cm)137.99 ± 1.165138.45 ± 0.698139.19 ± 1.9230.844
Association of CNVR30 of the CADM2 gene with yak growth traits. a,b,c Represent statistically significant differences at a level of p < 0.05; * p < 0.05.
Table 6. Correlation test of body weight in 6 months and CNV2.
Table 6. Correlation test of body weight in 6 months and CNV2.
Kendall Tau correlation coefficient0.261
Significance0.007 **
** Represent statistically significant differences at a level of p < 0.01.
Table 7. Two-way analysis of variance (ANOVA) of the effects of CADM2-CNVs in the Ashidan yak.
Table 7. Two-way analysis of variance (ANOVA) of the effects of CADM2-CNVs in the Ashidan yak.
AgeGrowth TraitCNVs (P Value)
CNV1CNV2CNV1 × CNV2
6 monthsBody weight (kg)0.1290.4850.884
Withers height (cm)0.4390.2950.418
Body length (cm)0.4700.9600.370
Chest girth(cm)0.9050.6520.387
12 monthsBody weight (kg)0.4470.9580.859
Withers height (cm)0.8700.9440.735
Body length (cm)0.6310.8910.996
Chest girth(cm)0.9970.9950.845
18 monthsBody weight (kg)0.6220.9410.354
Withers height (cm)0.003 **0.039 *0.190
Body length (cm)0.5080.8210.750
Chest girth(cm)0.1220.2530.168
* Represent statistically significant differences at a level of p < 0.05; ** Represent statistically significant differences at a level of p < 0.01.

Share and Cite

MDPI and ACS Style

Ge, F.; Jia, C.; Chu, M.; Liang, C.; Yan, P. Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak. Animals 2019, 9, 1008. https://doi.org/10.3390/ani9121008

AMA Style

Ge F, Jia C, Chu M, Liang C, Yan P. Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak. Animals. 2019; 9(12):1008. https://doi.org/10.3390/ani9121008

Chicago/Turabian Style

Ge, Fei, Congjun Jia, Min Chu, Chunnian Liang, and Ping Yan. 2019. "Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak" Animals 9, no. 12: 1008. https://doi.org/10.3390/ani9121008

APA Style

Ge, F., Jia, C., Chu, M., Liang, C., & Yan, P. (2019). Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak. Animals, 9(12), 1008. https://doi.org/10.3390/ani9121008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop