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

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

BACKGROUND

Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.

RESULTS

Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.

CONCLUSIONS

The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

The Identification of Goat KCNJ15 Gene Copy Number Variation and Its Association with Growth Traits

(1) Background: Copy number variation (CNV) is a critical component of genome structural variation and has garnered significant attention. High-throughput screening of the KCNJ15 gene has revealed a correlation between the CNV region and the growth traits of goats. We aimed to identify the CNV of the KCNJ15 gene in five goat breeds and analyze its association with growth characteristics. (2) Methods: We utilized 706 goats from five breeds: Guizhou black goat (GZB), Guizhou white goat (GZW), Bohuai goat (BH), Huai goat (HH), and Taihang goat (TH). To evaluate the number of copies of the KCNJ15 gene using qPCR, we analyzed the correlation between the CNV and growth characteristics and then used a universal linear model. The findings revealed variations in the distribution of different copy number types among the different goat breeds. (3) Results: Association analysis revealed a positive influence of the CNV in the KCNJ15 gene on goat growth. In GZB, individuals with duplication types exhibited superior performance in terms of cannon bone circumference (p < 0.05). In HH, individuals with duplication types exhibited superior performance in terms of body slanting length (p < 0.05). Conversely, normal TH demonstrated better body height and body weight (p < 0.05), while in GZW, when CN = 3, it performed better than other types in terms of body weight and chest circumference (p < 0.05). However, in BH, it had no significant effect on growth traits. (4) Conclusions: We confirmed that the CNV in the KCNJ15 gene significantly influences the growth characteristics of four distinct goat breeds. The correlation between KCNJ15 gene CNVs and goat growth traits offers valuable insights to breeders, enabling them to employ precise and efficient breeding methods that enhance livestock welfare, productivity, and overall economic benefits in the industry.

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

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

Copy Number Variations in the MICALL2 and MOGAT2 Genes Are Associated with Ashidan Yak Growth Traits

Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.

Genome-Wide Association Study of Body Weight Trait in Yaks

The yak is the largest meat-producing mammal around the Tibetan Plateau, and it plays an important role in the economic development and maintenance of the ecological environment throughout much of the Asian highlands. Understanding the genetic components of body weight is key for future improvement in yak breeding; therefore, genome-wide association studies (GWAS) were performed, and the results were used to mine plant and animal genetic resources. We conducted whole genome sequencing on 406 Maiwa yaks at 10 × coverage. Using a multiple loci mixed linear model (MLMM), fixed and random model circulating probability unification (FarmCPU), and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK), we found that a total of 25,000 single-nucleotide polymorphisms (SNPs) were distributed across chromosomes, and seven markers were identified as significantly (p-values < 3.91 × 10−7) associated with the body weight trait,. Several candidate genes, including MFSD4, LRRC37B, and NCAM2, were identified. This research will help us achieve a better understanding of the genotype−phenotype relationship for body weight.

Copy number variation (CNV) of the AHR gene in the Ashidan yak and its association with growth traits

Copy number variation (CNV) is a principal genomic structure variation affecting the gene expression through the dose-effect and change of gene regulatory region. It plays an important role in regulating the various complex traits of vertebrates. The aromatic hydrocarbon receptor (AHR) is a member of ligand-dependent transcription factors which belong to the alkaline helix-loop-helix PASS family. It is used as a conservative environmental sensor during biological evolution. This study, tracked the growth data (body weight, withers height, body length, chest girth) of 332 yaks in four stages (6, 12, 18, and 30 months) were tracked. The CNV of the yaks was analyzed using real-time quantitative PCR, and the correlation between CNV of AHR and yak growth traits was analyzed using the SPSS and R software. The AHR gene expression profiles were assessed in different tissues of the 18-month-old yak. The statistical analysis indicated the AHR-CNV of the Ashidan yak to significantly correlate with the body length (P < 0.05), and was found to be correlated with the withers height at 18 months old (P < 0.01) with extreme significance. To sum up, this study for the first time discussed the relationship between AHR-CNV and the growth traits of the Ashidan yak. The results indicated that the AHR gene might become a new molecular marker in the breeding yak.

Genome-Wide Placental Gene Methylations in Gestational Diabetes Mellitus, Fetal Growth and Metabolic Health Biomarkers in Cord Blood

Gestational diabetes mellitus (GDM) "program" an elevated risk of metabolic syndrome in the offspring. Epigenetic alterations are a suspected mechanism. GDM has been associated with placental DNA methylation changes in some epigenome-wide association studies. It remains unclear which genes or pathways are affected, and whether any placental differential gene methylations are correlated to fetal growth or circulating metabolic health biomarkers. In an epigenome-wide association study using the Infinium MethylationEPIC Beadchip, we sought to identify genome-wide placental differentially methylated genes and enriched pathways in GDM, and to assess the correlations with fetal growth and metabolic health biomarkers in cord blood. The study samples were 30 pairs of term placentas in GDM vs. euglycemic pregnancies (controls) matched by infant sex and gestational age at delivery in the Shanghai Birth Cohort. Cord blood metabolic health biomarkers included insulin, C-peptide, proinsulin, IGF-I, IGF-II, leptin and adiponectin. Adjusting for maternal age, pre-pregnancy BMI, parity, mode of delivery and placental cell type heterogeneity, 256 differentially methylated positions (DMPs,130 hypermethylated and 126 hypomethylated) were detected between GDM and control groups accounting for multiple tests with false discovery rate <0.05 and beta-value difference >0.05. WSCD2 was identified as a differentially methylated gene in both site- and region-level analyses. We validated 7 hypermethylated (CYP1A2, GFRA1, HDAC4, LIMS2, NAV3, PAX6, UPK1B) and 10 hypomethylated (DPP10, CPLX1, CSMD2, GPR133, NRXN1, PCSK9, PENK, PRDM16, PTPRN2, TNXB) genes reported in previous epigenome-wide association studies. We did not find any enriched pathway accounting for multiple tests. DMPs in 11 genes (CYP2D7P1, PCDHB15, ERG, SIRPB1, DKK2, RAPGEF5, CACNA2D4, PCSK9, TSNARE1, CADM2, KCNAB2) were correlated with birth weight (z score) accounting for multiple tests. There were no significant correlations between placental gene methylations and cord blood biomarkers. In conclusions, GDM was associated with DNA methylation changes in a number of placental genes, but these placental gene methylations were uncorrelated to the observed metabolic health biomarkers (fetal growth factors, leptin and adiponectin) in cord blood. We validated 17 differentially methylated placental genes in GDM, and identified 11 differentially methylated genes relevant to fetal growth.

Genome-Wide Association Study on Reproduction-Related Body-Shape Traits of Chinese Holstein Cows

Reproduction is an important production activity for dairy cows, and their reproductive performance can directly affect the level of farmers' income. To better understand the genomic regions and biological pathways of reproduction-related traits of dairy cows, in the present study, three body shape traits-Loin Strength (LS), Rump Angle (RA), and Pin Width (PW)-were selected as indicators of the reproductive ability of cows, and we conducted genome-wide association analyses on them. The heritability of these three traits was medium, ranging from 0.20 to 0.38. A total of 11 significant single-nucleotide polymorphisms (SNPs) were detected associated with these three traits. Bioinformatics analysis was performed on genes close to the significant SNPs (within 200 Kb) of LS, RA, and PW, and we found that these genes were totally enriched in 20 gene ontology terms and six KEGG signaling pathways. Finally, the five genes CDH12, TARP, PCDH9, DTHD1, and ARAP2 were selected as candidate genes that might affect LS. The six genes LOC781835, FSTL4, ATG4C, SH3BP4, DMP1, and DSPP were selected as candidate genes that might affect RA. The five genes USP6NL, CNTN3, LOC101907665, UPF2, and ECHDC3 were selected as candidate genes that might affect the PW of Chinese Holstein cows. Our results could provide useful biological information for the improvement of body shape traits and contribute to the genomic selection of Chinese Holstein cows.

Genome-wide selective sweep analysis of the high-altitude adaptability of yaks by using the copy number variant

The domestic yak (Bos grunniens) from the Qinghai-Tibet Plateau is an important animal model in high-altitude adaptation studies. Here, we performed the genome-wide selective sweep analysis to identify the candidate copy number variation (CNV) for the high-altitude adaptation of yaks. A total of 531 autosomal CNVs were determined from 29 yak genome-wide resequencing data (15 high- and 14 low-altitude distributions) by using a CNV caller with a CNV identification interval > 5 kb, CNV silhouette score > 0.7, and minimum allele frequency > 0.05. Most high-frequency CNVs were located at the exonic (44.63%) and intergenic (46.52%) regions. In accordance with the results of the selective sweep analysis, 7 candidate CNVs were identified from the interaction of the top 20 CNVs with highest divergence from the F _ST and V _ST between the low (LA) and high (HA) altitudes. Five genes (i.e., GRIK4, IFNLR1, LOC102275985, GRHL3, and LOC102275713) were also annotated from the seven candidate CNVs and their upstream and downstream ranges at 300 kb. GRIK4, IFNLR1, and LOC102275985 were enriched in five known signal pathways, namely, glutamatergic synapse, JAK-STAT signaling pathway, cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and olfactory transduction. These pathways are involved in the environmental adaptability and various physiological functions of animals, especially the physiological regulation under a hypoxic environment. The results of this study advanced the understanding of CNV as an important genomic structure variant type that contributes to HA adaptation and helped further explain the molecular mechanisms underlying the altitude adaptability of yaks.