Species composition and environmental adaptation of indigenous Chinese cattle

Temperature adaptation patterns in Chinese cattle revealed by TRPM2 gene mutation analysis

Cattle are sensitive to temperature fluctuations but adapt well to inclement weather conditions. When environmental temperatures exceed specific thresholds, heat stress becomes a critical concern for cattle. The TRPM2 gene, which resides on cattle chromosome 1 encodes a TRP channel protein, holding a unique capacity to sense temperature changes and facilitate rapid response to avoid heat stress. Here, we utilized the Bovine Genome Variation Database (BGVD) (http://animal.omics.pro/code/index.php/BosVar), and identified a missense mutation site, c.805A > G: p. Met269Val (rs527146862), within the TRPM2 gene. To elucidate the functional assessment of this mutation in temperature adaptation attributes of Chinese cattle, we genotyped 407 samples from 20 distinct breeds representing diverse climatic zones across China. The association analysis incorporates three temperature parameters and revealed compelling insights in terms of allele frequency. Interestingly, the prevalence of the wild-type allele A was notably higher among northern cattle breeds and this trend diminished gradually as observed in southern cattle populations. Conversely, the mutant-type allele G demonstrated a contrasting trend. Moreover, southern cattle exhibited markedly higher frequencies of GG and GA genotypes (P < 0.01). The presence of heterozygous and homozygous mutations appears to confer an enhanced capacity for adaptation to elevated temperatures. These results provide unequivocal correlation evidence between TRPM2 genotypes (AA, GA, GG) and environmental temperature parameters and comprehend the genetic mechanisms governing temperature adaptation in cattle. This provides valuable insights for strategic breed selection across diverse climatic regions, thereby aiding livestock production amid evolving climate challenges.

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

Evolutionary adaptations generally reverse phenotypic plasticity to restore ancestral phenotypes during new environment adaptation in cattle

Phenotype plasticity and evolution adaptations are the two main ways in which allow populations to deal with environmental changes, but the potential relationship between them remains controversial. Using a reciprocal transplant approach with cattle adapted to the Tibetan Plateau and adjacent lowlands, we aim to investigate the relative contributions of evolutionary processes and phenotypic plasticity in driving both phenotypic and transcriptomic changes under natural conditions. We observed that while numerous genetic transcriptomic changes were evident during the forward adaptation to highland environments, plastic changes predominantly facilitate the transformation of transcriptomes into a preferred state when Tibetan cattle are reintroduced to lowland habitats. Genes with ancestral plasticity are generally reversed by evolutionary adaptations and show a closer expression level to the ancestral stage in evolved Tibetan cattle. A similar trend was also observed at the phenotypes level, with a majority of biochemical and hemorheology phenotypes showing a tendency to revert to their ancestral patterns, suggesting the restoration of ancestral expression levels is a widespread evolutionary trend during adaptation. The findings of our study contribute to the debate regarding the relative contributions of plasticity and genetic changes in mammal environment adaptation. Furthermore, we highlight that the restoration of ancestral phenotypes represents a general pattern in cattle new environment adaptation.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Genetic factors influencing milk and fat yields in tropically adapted dairy cattle: insights from quantitative trait loci analysis and gene associations

OBJECTIVE

The objective of this study was to identify genes associated with 305-day milk yield (MY) and fat yield (FY) that also influence the adaptability of the Thai multibreed dairy cattle population to tropical conditions.

METHODS

A total of 75,776 imputed and actual single nucleotide polymorphisms (SNPs) from 2,661 animals were used to identify genomic regions associated with MY and FY using the single-step genomic best linear unbiased predictions. Fixed effects included herd-yearseason, breed regression, heterosis regression and calving age regression effects. Random effects were animal additive genetic and residual. Individual SNPs with a p-value smaller than 0.05 were selected for gene mapping, function analysis, and quantitative trait loci (QTL) annotation analysis.

RESULTS

A substantial number of QTLs associated with MY (9,334) and FY (8,977) were identified by integrating SNP genotypes and QTL annotations. Notably, we discovered 17 annotated QTLs within the health and exterior QTL classes, corresponding to nine unique genes. Among these genes, Rho GTPase activating protein 15 (ARHGAP15) and catenin alpha 2 (CTNNA2) have previously been linked to physiological traits associated with tropical adaptation in various cattle breeds. Interestingly, these two genes also showed signs of positive selection, indicating their potential role in conferring tolerance to trypanosomiasis, a prevalent tropical disease.

CONCLUSION

Our findings provide valuable insights into the genetic basis of MY and FY in the Thai multibreed dairy cattle population, shedding light on the underlying mechanisms of tropical adaptation. The identified genes represent promising targets for future breeding strategies aimed at improving milk and fat production while ensuring resilience to tropical challenges. This study significantly contributes to our understanding of the genetic factors influencing milk production and adaptability in dairy cattle, facilitating the development of sustainable genetic selection strategies and breeding programs in tropical environments.

Analysing the genetic diversity and population structure of five native Turkish cattle breeds using SNP data

The conservation and sustainable utilization of cattle genetic resources necessitate a comprehensive understanding of their genetic diversity and population structure. This study provides an analysis of five native Turkish cattle breeds: Anatolian Black (ANB), Turkish Grey (TUR), Anatolian Southern Yellow (ASY), East Anatolian Red (EAR), and South Anatolian Red (SAN) using 50 K SNP data. These breeds were compared with three European breeds, Simmental (SIM), Holstein (HOL), and Jersey (JER), and three Asian Zebu breeds: Arabic Zebu (ZAR), Nelore (NEL), and Red Sindhi (RSI). Genetic diversity indices demonstrated moderate heterogeneity among the breeds, with TUR exhibiting the highest observed heterozygosity (Ho = 0.35). Wright's Fst values indicated significant genetic differentiation, particularly between Turkish breeds and both European (Fst = 0.035-0.071) and Asian breeds (Fst = 0.025-0.150). Principal component analysis distinguished the unique genetic profiles of each breed cluster. Admixture analysis revealed degrees of shared genetic ancestry, suggesting historical gene flow between Turkish, European, and Asian breeds. Analysis of molecular variance (AMOVA) attributed approximately 58% of the variation to population differences. Nei's genetic distances highlighted the closer genetic relatedness within Turkish breeds (distance ranges between 0.032 and 0.046) and suggested a more relative affinity of TUR with European breeds. The study's phylogenetic assessments elucidate the nuanced genetic relationships among these breeds, with runs of homozygosity (ROH) analysis indicating patterns of ancestral relatedness and moderate levels of inbreeding, particularly evident in Turkish breeds. Our findings provide valuable insights into the genetic landscape of Turkish cattle, offering a crucial foundation for informed conservation and breeding strategies aimed at preserving these breeds' genetic integrity and heritage.

Deciphering climate resilience in Indian cattle breeds by selection signature analyses

The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.

Genetic Origin and Introgression Pattern of Pingliang Red Cattle Revealed Using Genome-Wide SNP Analyses

The Pingliang red cattle, an outstanding indigenous resource in China, possesses an exceptional breeding value attributed to its tender meat and superior marbling quality. Currently, research efforts have predominantly concentrated on exploring its maternal origin and conducting conventional phenotypic studies. However, there remains a lack of comprehensive understanding regarding its genetic basis. To address this gap, we conducted a thorough whole-genome analysis to investigate the population structure, phylogenetic relationships, and gene flows of this breed using genomic SNP chip data from 17 bovine breeds. The results demonstrate that Pingliang red cattle have evolved distinct genetic characteristics unique to this breed, clearly distinguishing it from other breeds. Based on the analysis of the population structure and phylogenetic tree, it can be classified as a hybrid lineage between Bos taurus and Bos indicus. Furthermore, Pingliang red cattle display a more prominent B. taurus pedigree in comparison with Jinnan, Qinchuan, Zaosheng, Nanyang, and Luxi cattle. Moreover, this study also revealed closer genetic proximity within the Chinese indigenous cattle breed, particularly Qinchuan cattle, which shares the longest identical by descent (IBD) fragment with Pingliang red cattle. Gene introgression analysis shows that Pingliang red cattle have undergone gene exchange with South Devon and Red Angus cattle from Europe. Admixture analysis revealed that the proportions of East Asian taurine and Chinese indicine in the ancestry of Pingliang red cattle are approximately 52.44% and 21.00%, respectively, while Eurasian taurine, European taurine, and Indian indicine account for approximately 17.55%, 7.27%, and 1.74%. Our findings unveil distinct genetic characteristics in Pingliang red cattle and attribute their origin to B. taurus and B. indicus ancestry, as well as contributions from Qinchuan cattle, South Devon, and Red Angus.

A missense mutation of the WNK1 gene affects cold tolerance in Chinese domestic cattle

Inclement weather conditions, especially cold stress, have threatened the cattle industry. Cattle exposed to cold environments for a longer time suffer developmental delay, immunity decline, and eventually death. WNK1 is a member of With-no-lysine kinases (WNKs), widely expressed in animal organs and tissues. WNK1 and WNK4 are expressed in adipose tissue, and WNK4 promotes adipogenesis. WNK1 does not directly affect adipogenesis but has been shown to promote WNK4 expression in several tissues or organs. One missense mutation NC_037346.1:g.107692244, A > G, rs208265410 in the WNK1 gene was detected from the database of bovine genomic variation (BGVD). Here, we collected 328 individuals of 17 breeds representing four groups of Chinese cattle, northern group cattle, southern group cattle, central group cattle, and special group cattle (Tibetan cattle). We also collected the temperature and humidity data records from their relative locations. The frequencies of the G allele in Chinese breeds increased from northern China to southern China, and the frequencies of the A allele showed an opposite trend. Our results indicate that the WNK1 gene might be a candidate gene marker associated with cold tolerance.

TRPM4 gene variation associated with climatic conditions in Chinese cattle

The transient receptor potential (TRP) superfamily has been reported to play an important role in heat tolerance pathways. Based on the Bovine Genome Variation Database and Selective Signatures, a missense mutation (NC_037345.1: c.2237A > G: p. His746Arg) (rs209689836) was identified in the transient receptor potential cation channel subfamily M member 4 (TRPM4) gene, a member of the TRP family, corresponding to heat tolerance. Here, we explored the prevalence of this variant in 19 native Chinese cattle (comprised of 404 individuals) to determine its possible association with heat tolerance in Chinese cattle by using PCR and DNA sequencing. The distribution of alleles of NC_037345.1: c.2237A > G: p. His746Arg displays significant geographical differences across native Chinese cattle breeds, consistent with the distribution of indicine and taurine cattle in China. Additionally, the association analysis indicated that the G allele was significantly associated with mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) (p < .05), suggesting that cattle carrying allele G were distributed in regions with higher T, RH, and THI. In conclusion, our results suggested that the mutation of the TRPM4 gene in Chinese cattle might be a candidate locus associated with heat tolerance.

A missense mutation (rs209302038) of KRT9 gene associated with heat stress in Chinese cattle

Type I keratin 9 encoded by the KRT9 gene serves an important special function either in the mature palmar and plantar skin tissue. The changes in skin conditions and thickening of the outer layer of the skin may be affected by environmental variables. A missense mutation rs209302038 (NC_037346.1: g.41782870 G > A) was detected in KRT9, which changing the isoleucine into valine. This study aimed to identify the frequency of allele in this locus in Chinese indigenous cattle, and analyze the connection with heat stress. Our results indicated that the frequency of allele A gradually decreases from south to north, while the frequency of G allele showed the opposite pattern. Further analysis of the association of the different genotypes with three climate factors, which showed that the genotypes (GG, GA, AA) were significantly related to climatic conditions (p < 0.01). Therefore, we speculated that the mutation of the rs209302038 in Chinese indigenous cattle might be a genetic marker to detect heat stress.

A novel SNP of TECPR2 gene associated with heat tolerance in Chinese cattle

Heat stress affects the animal production and causes serious economic losses to the husbandry. Tectonin beta-propeller repeat containing 2 (TECPR2) gene plays an important role in autophagy which may affect the temperature sensation in animals. A missense mutation (XM_024981840.1:c.3989 G > A p.Arg1330His) of the transcripts X4 in the bovine TECPR2 gene was identified. In this study, the c.3989 G > A variant in TECPR2 gene was genotyped in a total of 25 cattle breeds (520 individuals). Our results indicated that the frequency of A allele showed a decreasing pattern from southern cattle to northern cattle, while the frequency of G allele showed the opposite pattern, which was consistent with the climate distribution of China. Compared with the GG genotype, southern cattle carried more the AA and AG genotypes. Furthermore, the association results carried out that the frequencies of genotypes (GG, AG, AA) and the value of climate parameters (mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) were significantly correlated (p < 0.01). Hence, we speculated that the c.3989 G > A variant of TECPR2 gene was associated with the heat tolerance trait in Chinese cattle and the locus may be considered as a molecular marker for Chinese cattle breeding.

Open AccessLack of association between coat color abnormalities in Bali cattle (Bos javanicus) and the coding regions of the MC1R and KIT genes

Background and Aim

Coat color variations in cattle are known to be influenced by the melanocortin 1 receptor (MC1R) and receptor tyrosine kinase (KIT) genes. The presence of coat color abnormalities, such as white spots and albinism, in Bali cattle was the focus of this study. This study aimed to identify single nucleotide polymorphisms (SNPs) in the coding region of MC1R and exons 2 and 3 of KIT associated with coat color abnormalities in Bali cattle.

Materials and Methods

The study included the analysis of 48 Bali cattle, including 20 individuals with standard coat color, 15 with white spots, and 13 with albinism. Total DNA was extracted using a DNA Extraction Kit, and MC1R (coding region) and KIT (exons 2 and 3) gene amplifications were analyzed using forward and reverse primers with polymerase chain reaction product lengths of 1071, 234, and 448 bp, respectively. The determination of MC1R and KIT gene diversity was analyzed through direct sequencing. Melanocortin 1 receptor and KIT gene sequence data were analyzed using BioEdit and MEGA6 to identify SNPs associated with standard and abnormal coat color phenotypes (white-spotted and albino) in Bali cattle.

Results

No SNPs associated with coat color abnormalities were found in the coding region of MC1R and exons 2 and 3 of KIT genes in Bali cattle. However, the intron two regions of KIT contained the SNP g.70208534A>G, which showed a high degree of diversity. The AA genotype frequency was highest in albino Bali cattle, whereas the G allele frequency was highest and the A allele frequency was lowest in white-spotted Bali cattle.

Conclusion

The results indicated that standard, white-spotted, and albino coat colors in Bali cattle could not be distinguished by analyzing the MC1R and KIT genes.

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Deciphering local adaptation of native Indian cattle (Bos indicus) breeds using landscape genomics and in-silico prediction of deleterious SNP effects on protein structure and function

India has 50 registered breeds of native cattle (Bos indicus) which are locally adapted to diverse environmental conditions. This study aimed to investigate the genomic basis of adaptation of native Indian cattle and to predict the impact of key SNPs on the amino acid changes that affect protein function. The Illumina 777 K BovineHD BeadChip was used to genotype 178 native cattle belonging to contrasting landscapes and agro-climatic conditions. The genotype-environment association was investigated with R. SamBada, using 5,74,382 QC passed SNPs and 11 predictor variables (10 multi-collinearity controlled environmental variables and 1 variable as "score of PCA" on ancestry coefficients of individuals). In total, 1,12,780 models were selected as significant (q < 0.05) based on G score. The pathway ontology of the annotated genes revealed many important pathways and genes having a direct and indirect role in cold and hot adaptation. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Adaptive alleles identified in the present investigation might be responsible for the local adaptation of these cattle breeds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03493-3.

Metformin attenuates depressive-like behaviour of methamphetamine withdrawal in mice: A mechanistic approach

OBJECTIVES

Methamphetamine (METH) as a potent psychostimulant drug with a high potency of dependence rate that results in neurotoxicity has become a major drug of abuse in many parts of the world. Unfortunately, there is limited evidence regarding treatment of METH withdrawal syndrome. Therefore, we aimed to investigate whether metformin mitigate the methamphetamine (METH) withdrawal syndrome in male mice. Based on the literature, depression and anxiety are the major METH withdrawal symptoms.

METHODS

Here, METH (2 mg/kg) was administered to mice twice a day for 14 constitutive days to induce animal model of METH-induced withdrawal syndrome. To do this, mice in control group and those with METH withdrawal syndrome were divided into treatment (receiving metformin in 3 doses of 50, 100 and 200 mg/kg for 10 days) and non-treatment sub-groups. Following the behavioural test, the animals were sacrificed; their hippocampus was dissected to measure oxidative stress parameters and expression of cellular energy homeostasis and immune-inflammatory genes.

RESULTS

Our data revealed that metformin provoked antidepressant effects in behavioural tests through AMPK overexpression as an important mitochondrial energetic sensor and inhibition of Tlr4 overexpression in the immune system gene expression. In addition, metformin was able to improve oxidative stress biomarkers and neuronal damage in the hippocampus and restore cellular energy homeostasis and immune system gene expression.

CONCLUSIONS

The data suggested that metformin can influence the hippocampus through targeting mitochondria and their performance, and consequently, neuroinflammation responses and brain metabolic changes. It is supposed to be a new therapeutic option in clinical trials of depression and anxiety following METH withdrawal treatment.

Favored single nucleotide variants identified using whole genome Re-sequencing of Austrian and Chinese cattle breeds

Cattle have been essential for the development of human civilization since their first domestication few thousand years ago. Since then, they have spread across vast geographic areas following human activities. Throughout generations, the cattle genome has been shaped with detectable signals induced by various evolutionary processes, such as natural and human selection processes and demographic events. Identifying such signals, called selection signatures, is one of the primary goals of population genetics. Previous studies used various selection signature methods and normalized the outputs score using specific windows, in kbp or based on the number of SNPs, to identify the candidate regions. The recent method of iSAFE claimed for high accuracy in pinpointing the candidate SNPs. In this study, we analyzed whole-genome resequencing (WGS) data of ten individuals from Austrian Fleckvieh (Bos taurus) and fifty individuals from 14 Chinese indigenous breeds (Bos taurus, Bos taurus indicus, and admixed). Individual WGS reads were aligned to the cattle reference genome of ARS. UCD1.2 and subsequently undergone single nucleotide variants (SNVs) calling pipeline using GATK. Using these SNVs, we examined the population structure using principal component and admixture analysis. Then we refined selection signature candidates using the iSAFE program and compared it with the classical iHS approach. Additionally, we run Fst population differentiation from these two cattle groups. We found gradual changes of taurine in north China to admixed and indicine to the south. Based on the population structure and the number of individuals, we grouped samples to Fleckvieh, three Chinese taurines (Kazakh, Mongolian, Yanbian), admixed individuals (CHBI_Med), indicine individuals (CHBI_Low), and a combination of admixed and indicine (CHBI) for performing iSAFE and iHS tests. There were more significant SNVs identified using iSAFE than the iHS for the candidate of positive selection and more detectable signals in taurine than in indicine individuals. However, combining admixed and indicine individuals decreased the iSAFE signals. From both within-population tests, significant SNVs are linked to the olfactory receptors, production, reproduction, and temperament traits in taurine cattle, while heat and parasites tolerance in the admixed individuals. Fst test suggests similar patterns of population differentiation between Fleckvieh and three Chinese taurine breeds against CHBI. Nevertheless, there are genes shared only among the Chinese taurine, such as PAX5, affecting coat color, which might drive the differences between these yellowish coated breeds, and those in the greater Far East region.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Cloning and Expression Analysis of HAT1 and HDAC1 in the Testes of Mature Yaks and Their Sterile Hybrids

Simple Summary

Cattle-yak is the hybrid between male cattle (Bos taurus) and female yak (Bos grunniens). Male cattle-yak can not produce normal sperm. The mechanisms that underlie cattle-yak male sterility have not been elucidated. Histone acetylation is a common regulation mode that plays an important role in the development of gametes. The objective of this study was to explore the molecular mechanism of male sterility in yak hybrids based on histone acetyltransferase 1 (HAT1) and histone deacetylase 1 (HDAC1), two enzymes that regulate histone acetylation. The mRNA and protein expression levels of HAT1 in the testes of adult cattle-yaks were significantly lower than in adult yaks, and the protein expression levels of HDAC1 were significantly higher than in yaks. In addition, H3K9 acetylation levels in cattle-yak testes were significantly lower than in yaks. These results suggest that male sterility in cattle-yaks might be associated with decreased histone acetylation levels in the testes.

Abstract

The objective of this study was to explore the molecular mechanism of male sterility in yak hybrids based on HAT1 and HDAC1. Total RNA was extracted from the testes of adult yaks (n = 11) and sterile cattle-yaks (n = 11) followed by reverse transcription. The coding sequence (CDS) of yak HAT1 and HDAC1 were obtained by conventional polymerase chain reaction (PCR) and gene cloning. The testicular mRNA and protein levels of HAT1 and HDAC1 in yaks and cattle-yaks were detected by quantitative PCR (qPCR) and Western blotting, respectively, and the histone H3 lysine 9 (H3K9) histone acetylation level in the testes of yaks and cattle-yaks was assayed using enzyme linked immunosorbent assay (ELISA). The results showed that the CDS of HAT1 and HDAC1 were 1242 bp and 1449 bp in length, encoding 413 and 482 amino acids, respectively; yaks had a similar mRNA sequence as cattle in both genes. The testicular mRNA and protein levels of HAT1 of cattle-yaks were significantly lower than those of yaks, and the protein level of HDAC1 was significantly higher than that of yaks. ELISA showed that the acetylation level of testicular H3K9 was significantly lower in yak hybrids than that of yaks. The present results suggest that the decreased level of HAT1 and increased level of HDAC1 may result in the decreased H3K9 acetylation in cattle-yaks and might be associated with their sterility.

Whole-Genome Analysis Deciphers Population Structure and Genetic Introgression Among Bovine Species

It is known that throughout history and presently, taurine (Bos taurus) and indicine/zebu (Bos indicus) cattle were crossed with other bovine species (e.g., gayal, gaur, banteng, yak, wisent, and bison). Information on the role of interspecific hybridization to facilitate faster adaptation of the newly arrived domestic species to new environments is poorly known. Herein, we collected 266 samples of bovine species of the taurine, zebu, yak, and gaur from West Europe, Southwest China, Indian subcontinent, and Southeast Asia to conduct the principal component analysis (PCA), admixture, gene flow, and selection signature analyses by using SNPs distributed across the bovine autosomes. The results showed that the genetic relationships between the zebu, yak, and gaur mirrored their geographical origins. Three ancestral components of the European taurine, East Asian taurine, and Indian zebu were found in domestic cattle, and the bidirectional genetic introgression between the Diqing cattle and Zhongdian yak was also detected. Simultaneously, the introgressed genes from the Zhongdian yak to the Diqing cattle were mainly enriched with immune-related pathways, and the ENPEP, FLT1, and PIK3CA genes related to the adaptation of high-altitude hypoxia were detected. Additionally, we found the genetic components of the Zhongdian yak had introgressed into Tibetan cattle. The 30 selected genes were detected in Tibetan cattle, which were significantly enriched in the chemokine signaling pathway. Interestingly, some genes (CDC42, SLC39A2, and EPAS1) associated with hypoxia response were discovered, in which CDC42 and SLC39A2 played important roles in angiogenesis and erythropoiesis, and heart function, respectively. This result showed that genetic introgression was one of the important ways for the environmental adaptation of domestic cattle.

Distribution of a missense mutation (rs525805167) within the SLC45A2 gene associated with climatic conditions in Chinese cattle

SLC45A2 is involved in the synthesis of melanin transporters. We investigated the association between single nucleotide polymorphisms (SNPs) of the SLC45A2 gene and humidity and hot conditions in indigenous cattle habitat. According to the Bovine Genome Variation Database and Selective Signatures (BGVD), we explored the frequency distribution of a missense mutation (NC_037347.1: c.1543A > G, p.ser515gly) in the SLC45A2 gene in Chinese indigenous cattle. This variation from serine to glycine caused a significant change in the protein modeling structure. PCR and partial DNA sequencing were used to genotype 541 individuals, including 28 Chinese indigenous cattle breeds as well as Angus and zebu. From our results, the mutant allele frequency of this SNP in Chinese native cattle increases gradually from north to south, which is consistent with the distribution of climatic conditions in China. In addition, according to association analysis of a missense mutation (NC_037347.1: c.1543A > G) (rs525805167) in Chinese cattle, it is closely related to the annual average temperature (T), relative humidity (RH), temperature and humidity index (THI) and solar radiation time (P < 0.01). Based on the statistical analysis of the data, we assumed that rs525805167 was associated with heat tolerance traits. Simple Summary: The characteristics of Chinese indigenous cattle are closely related to their climatic environment. In China, Bos taurus is mainly distributed in the northern regions; Bos indicus is mainly distributed in southern China. In addition, the average temperature is higher in the south than in the north, and there are many mixed ancestry breeds of B. taurus and B. indicus in the middle area. The SLC45A2 gene is related to melanin synthesis, which may be closely related to heat tolerance in cattle. The purpose of our study was to investigate whether the SLC45A2 gene is related to heat tolerance in Chinese indigenous cattle.

A novel missense mutation (rs464874590) within BoLA-DOB gene associated with the heat-resistance in Chinese cattle

Bovine leukocyte antigen, class II, DO beta (BoLA-DOB) is related to antigen presentation, which can triggered by multicul factors. And the condition of immune function determines how much cattle load to heat stress. To evaluate the relationship between heat-resistance and single nucleotide polymorphisms (SNPs) in BoLA-DOB gene, our study has taken further analysis in Chinese indigenous cattle for the first time. A missense single nucleotide polymorphism (rs464874590) was detected in BoLA-DOB gene. We directly sequenced rs464874590 (NM_001013600.1 g.7122762 A > G) in BoLA-DOB gene of 522 individuals of 26 cattle breeds. The frequency of allele G gradually decreases from south to north with distinct climatic distribution characteristics. Further association analysis was carried out between different genotypes and environmental parameters, including annual mean temperature (T), relative humidity (RH), and temperature-humidity index (THI). The result showed that three genotypes were significantly correlated with T, H, and THI (P < 0.01), indicating that GG genotype was distributed in areas with hot and moist conditions. Therefore, our results suggested that the rs464874590 could be applied as a genetic marker to detect the heat-resistance of Chinese indigenous cattle.

TAS2R16 introgression from banteng into indigenous Chinese cattle

Historical hybridization between southern indigenous Chinese cattle and banteng has been well-documented and has resulted in gene introgression. Bitter taste receptors were reported in indigenous cattle as a result of introgression from banteng. To determine the level of introgression of the taste 2 receptor member 16 (TAS2R16) gene from banteng into Chinese cattle, two missense mutations in the bovine TAS2R16 gene were examined. Here, we explored the prevalence of the two variants in 28 indigenous Chinese cattle and banteng breeds (comprising 750 individuals) to determine the influence of banteng introgressions on Chinese cattle based on PCR and DNA sequencing. In our study, the two mutant alleles had a higher frequency distribution in southern China with strong geographic distribution, especially in the south-central and southeast areas. In conclusion, this study examines the impact of introgression on the frequency distributions of mutations in variable regions and the subsequent adaptation of Chinese indigenous cattle to different environmental conditions.

Unraveling Admixture, Inbreeding, and Recent Selection Signatures in West African Indigenous Cattle Populations in Benin

The Dwarf Lagune and the Savannah Somba cattle in Benin are typical representatives of the endangered West African indigenous Shorthorn taurine. The Lagune was previously exported to African and European countries and bred as Dahomey cattle, whereas the Somba contributed to the formation of two indigenous hybrids known as Borgou and Pabli cattle. These breeds are affected by demographic, economic, and environmental pressures in local production systems. Considering current and historical genomic data, we applied a formal test of admixture, estimated admixture proportions, and computed genomic inbreeding coefficients to characterize the five breeds. Subsequently, we unraveled the most recent selection signatures using the cross-population extended haplotype homozygosity approach, based on the current and historical genotypes. Results from principal component analyses and high proportion of Lagune ancestry confirm the Lagune origin of the European Dahomey cattle. Moreover, the Dahomey cattle displayed neither indicine nor European taurine (EUT) background, but they shared on average 40% of autozygosity from common ancestors, dated approximately eight generations ago. The Lagune cattle presented inbreeding coefficients larger than 0.13; however, the Somba and the hybrids (Borgou and Pabli) were less inbred (≤0.08). We detected evidence of admixture in the Somba and Lagune cattle, but they exhibited a similar African taurine (AFT) ancestral proportion (≥96%) to historical populations, respectively. A moderate and stable AFT ancestral proportion (62%) was also inferred for less admixed hybrid cattle including the Pabli. In contrast, the current Borgou samples displayed a lower AFT ancestral proportion (47%) than historical samples (63%). Irrespective of the admixture proportions, the hybrid populations displayed more selection signatures related to economic traits (reproduction, growth, and milk) than the taurine. In contrast, the taurine, especially the Somba, presented several regions known to be associated with adaptive traits (immunity and feed efficiency). The identified subregion of bovine leukocyte antigen (BoLA) class IIb (including DSB and BOLA-DYA) in Somba cattle is interestingly uncommon in other African breeds, suggesting further investigations to understand its association with specific adaptation to endemic diseases in Benin. Overall, our study provides deeper insights into recent evolutionary processes in the Beninese indigenous cattle and their aptitude for conservation and genetic improvement.

Bos taurus haplotypes segregating in Nellore (Bos indicus) cattle

Brazil is the largest exporter of beef in the world, and most of that beef derives from Nellore cattle. Although considered a zebu breed (Bos indicus), the history of Nellore cattle in Brazil is marked by the importation of bulls from India, the use of a Creole taurine (Bos taurus) maternal lineage to quickly expand the herds and backcrossing to Nellore bulls to recover zebu ancestry. As a consequence, the current Brazilian Nellore population carries an average taurine ancestry of approximately 1%. Although that percentage seems small, some taurine variants deviate substantially from that average, with the better-known cases being the PLAG1-Q haplotype involved with body size variation and the Guarani (P_G ) polled variant producing hornless animals. Here, we report taurine haplotypes in 9074 Nellore animals genotyped for 539 657 imputed SNP markers. Apart from PLAG1-Q and P_G , our analysis further revealed common taurine haplotypes (>3%) spanning genes related to immunity, growth, reproduction and hair and skin phenotypes. Using data from 22 economically important traits, we showed that many of the major QTL previously reported in the breed are at least partially driven by taurine haplotypes. As B. taurus and B. indicus haplotypes are highly divergent, presenting widely different sets of functional variants, our results provide promising targets for future scrutiny in Nellore 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.

Genome analyses revealed genetic admixture and selection signatures in Bos indicus

The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (F_ST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.

Assessing Genomic Diversity and Productivity Signatures in Dianzhong Cattle by Whole-Genome Scanning

Dianzhong cattle is a classic Chinese indigenous cattle breed with historical records dating back to 200 BC. But with its genomic differences having not been clearly elucidated, the quest for genomic characterization will be an essential step towards understanding the genomic basis of productivity and adaptation to survival under Chinese farming systems. Here we compared 10 Dianzhong cattle (four newly sequenced and six downloaded) with 29 published genomes of three underlying ancestral populations (Chinese zebu, Indian zebu, and Yanbian cattle) to characterize the genomic variations of Dianzhong cattle. Dianzhong cattle has a high nucleotide diversity (0.0034), second only to Chinese zebu. Together with analyses of linkage disequilibrium decay and runs of homozygosity, Dianzhong cattle displayed higher genomic diversity and weaker artificial selection compared with Yanbian cattle. From a selective sweep analysis by four methods (Fst, π-ratio, XP-CLR, and XP-EHH), the positive selective signals were mainly manifested in candidate genes and pathways related to heat resistance, growth and development, fat deposition, and male reproduction. Missense mutations were detected in candidate genes, SDS (c.944C > A and p.Ala315Glu), PDGFD (c.473A > G and p.Lys158Arg), and DDX4 (rs460251486, rs722912933, and rs517668236), which related to heat resistance, fat deposition, and spermatogenesis, respectively. Our findings unravel, at the genome-wide level, the unique diversity of Dianzhong cattle while emphasizing the opportunities for improvement of livestock productivity in further breeding programs.

Genomic diversity and breed composition of Vietnamese smallholder dairy cows

Vietnamese smallholder dairy cows (VDC) are the result of crossbreeding between different zebu (ZEB) and taurine dairy breeds through many undefined generations. Thus, the predominant breed composition of VDC is currently unknown. This study aimed to evaluate the level of genetic diversity and breed composition of VDC. The SNP data of 344 animals from 32 farms located across four dairy regions of Vietnam were collected and merged with genomic reference data, which included three ZEB breeds: Red Sindhi, Sahiwal and Brahman, three taurine breeds: Holstein (HOL), Jersey (JER) and Brown Swiss (BSW), and a composite breed: Chinese Yellow cattle. Diversity and admixture analyses were applied to the merged data set. The VDC were not excessively inbred, as indicated by very low inbreeding coefficients (Wright's F_IS ranged from -0.017 to 0.003). The genetic fractions in the test herds suggested that the VDC are primarily composed of HOL (85.0%); however, JER (6.0%), BSW 5.3%) and ZEB (4.5%) had also contributed. Furthermore, major genotype groupings in the test herds were pure HOL (48%), B3:15/16HOL_1/16ZEB (22%) and B2:7/8HOL_1/8ZEB (12%). The genetic makeup of the VDC is mainly components of various dairy breeds but also has a small percentage of ZEB; thus, the VDC could be a good genetic base for selecting high milk-producing cows with some degree of adaptation to tropical conditions.

The distribution of a missense mutation in SDK1 gene across native Chinese breeds

Sidekick cell adhesion molecule 1 (SDK1) gene belonging to the immunoglobulin superfamily (IgSF) is reported to be associated with disease resistance. A novel missense mutation (XM_015469413.1: c.2678 A > G: p. Asp893Gly) in bovine SDK1 gene was observed in the Bovine Genome Variation Database (BGVD). The purpose of the current study was to determine the allelic frequency distribution of XM_015469413.1: c.2678 A > G: p. Asp893Gly and analyze its association with disease resistance in native Chinese cattle. Here, we explored the allele frequency of the missense mutation (M_015469413.1: c.2678 A > G: p. Asp893Gly) in 542 individuals from 27 Chinese indigenous cattle breeds using PCR and DNA sequencing methods. The frequency of the mutant allele (G) gradually decreased from the southern cattle groups to the northern cattle groups, whereas the frequency of the wild-type allele A showed an opposite pattern, consistent with the distribution of indicine and taurine cattle in China. Thus, our results proved that the bovine SDK1 gene might be candidate genes associated with disease resistance in Chinese cattle.

Genome-wide SNP analysis clearly distinguished the Belarusian Red cattle from other European cattle breeds

Local breeds can serve as an important source of genetic variability in domestic animal species. This study aimed to assess the genetic diversity and population structure of Belarusian Red cattle and their differentiation from other European cattle populations based on genome-wide SNP genotypes. Twenty pedigree-recorded non-closely related cows of Belarusian Red cattle were genotyped using the Illumina BovineHD BeadChip. Genotypes of 22 other European cattle breeds were included in the study for comparison. A total of 28 562 SNPs passed through the quality control checks and were selected for analysis. The Belarusian Red cattle displayed a moderate level of genetic variability (_U H_E  = 0.341, H_O  = 0.368), and the highest heterozygote excess (_U F_IS  = -0.066), among the studied breeds; this reflects the contribution of multiple breeds to their formation. The principal component analysis, F_ST -based Neighbor-Net tree and Admixture clustering, clearly distinguished the Belarusian Red cattle from the other European cattle breeds. Moreover, the presence of ancestral genomic components of Danish Red and Brown Swiss breeds were clearly visible, which agrees with the breed's history and its recent development. Our study highlights the importance of maintaining the specific genomic components, which makes a significant contribution to the global genetic diversity in the modern population of Belarusian Red cattle, allowing us to consider them a valuable national genetic resource. Our research results will be useful for the development of conservation programs for this local cattle breed.

Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains

Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.

Mitochondrial genomes from modern and ancient Turano-Mongolian cattle reveal an ancient diversity of taurine maternal lineages in East Asia

Turano-Mongolian cattle are a group of taurine cattle from Northern and Eastern Asia with distinct morphological traits, which are known for their ability to tolerate harsh environments, such as the Asian steppe and the Tibetan plateau. Through the analysis of 170 mitogenomes from ten modern breeds, two sub-lineages within T3 (T3₁₁₉ and T3₀₅₅) were identified as specific of Turano-Mongolian cattle. These two T3 sub-lineages, together with the previously identified T4, were also present in six Neolithic samples, dated to ~3900 years BP, which might represent the earliest domestic taurine stocks from Southwest Asia. The rare haplogroup Q, found in three Tibetan cattle, testifies for the legacy of ancient migrations from Southwest Asia and suggests that the isolated Tibetan Plateau preserved unique prehistoric genetic resources. These findings confirm the geographic substructure of Turano-Mongolian cattle breeds, which have been shaped by ancient migrations and geographic barriers.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.

Introgression, admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.

Assessing the genetic background and genomic relatedness of red cattle populations originating from Northern Europe

Background

Local cattle breeds need special attention, as they are valuable reservoirs of genetic diversity. Appropriate breeding decisions and adequate genomic management of numerically smaller populations are required for their conservation. At this point, the analysis of dense genome-wide marker arrays provides encompassing insights into the genomic constitution of livestock populations. We have analyzed the genetic characterization of ten cattle breeds originating from Germany, The Netherlands and Denmark belonging to the group of red dairy breeds in Northern Europe. The results are intended to provide initial evidence on whether joint genomic breeding strategies of these populations will be successful.

Results

Traditional Danish Red and Groningen White-Headed were the most genetically differentiated breeds and their populations showed the highest levels of inbreeding. In contrast, close genetic relationships and shared ancestry were observed for the populations of German Red and White Dual-Purpose, Dutch Meuse-Rhine-Yssel, and Dutch Deep Red breeds, reflecting their common histories. A considerable amount of gene flow from Red Holstein to German Angler and to German Red and White Dual-Purpose was revealed, which is consistent with frequent crossbreeding to improve productivity of these local breeds. In Red Holstein, marked genomic signatures of selection were reported on chromosome 18, suggesting directed selection for important breeding goal traits. Furthermore, tests for signatures of selection between Red Holstein, Red and White Dual-Purpose, and Meuse-Rhine-Yssel uncovered signals for all investigated pairs of populations. The corresponding genomic regions, which were putatively under different selection pressures, harboured various genes which are associated with traits such as milk and beef production, mastitis and female fertility.

Conclusions

This study provides comprehensive knowledge on the genetic constitution and genomic connectedness of divergent red cattle populations in Northern Europe. The results will help to design and optimize breeding strategies. A joint genomic evaluation including some of the breeds studied here seems feasible.

Genetic Diversity and Signatures of Selection for Thermal Stress in Cattle and Other Two Bos Species Adapted to Divergent Climatic Conditions

Understanding the biological mechanisms of climatic adaptation is of paramount importance for the optimization of breeding programs and conservation of genetic resources. The aim of this study was to investigate genetic diversity and unravel genomic regions potentially under selection for heat and/or cold tolerance in thirty-two worldwide cattle breeds, with a focus on Chinese local cattle breeds adapted to divergent climatic conditions, Datong yak (Bos grunniens; YAK), and Bali (Bos javanicus) based on dense SNP data. In general, moderate genetic diversity levels were observed in most cattle populations. The proportion of polymorphic SNP ranged from 0.197 (YAK) to 0.992 (Mongolian cattle). Observed and expected heterozygosity ranged from 0.023 (YAK) to 0.366 (Sanhe cattle; SH), and from 0.021 (YAK) to 0.358 (SH), respectively. The overall average inbreeding (±SD) was: 0.118 ± 0.028, 0.228 ± 0.059, 0.194 ± 0.041, and 0.021 ± 0.004 based on the observed versus expected number of homozygous genotypes, excess of homozygosity, correlation between uniting gametes, and runs of homozygosity (ROH), respectively. Signatures of selection based on multiple scenarios and methods (F ST, HapFLK, and ROH) revealed important genomic regions and candidate genes. The candidate genes identified are related to various biological processes and pathways such as heat-shock proteins, oxygen transport, anatomical traits, mitochondrial DNA maintenance, metabolic activity, feed intake, carcass conformation, fertility, and reproduction. This highlights the large number of biological processes involved in thermal tolerance and thus, the polygenic nature of climatic resilience. A comprehensive description of genetic diversity measures in Chinese cattle and YAK was carried out and compared to 24 worldwide cattle breeds to avoid potential biases. Numerous genomic regions under positive selection were detected using three signature of selection methods and candidate genes potentially under positive selection were identified. Enriched function analyses pinpointed important biological pathways, molecular function and cellular components, which contribute to a better understanding of the biological mechanisms underlying thermal tolerance in cattle. Based on the large number of genomic regions identified, thermal tolerance has a complex polygenic inheritance nature, which was expected considering the various mechanisms involved in thermal stress response.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.

A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species

BACKGROUND

Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available.

RESULTS

The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep).

CONCLUSIONS

We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.

Discovery of Genomic Characteristics and Selection Signatures in Southern Chinese Local Cattle

Chinese local cattle with a high level of genetic diversity mainly originate from two subspecies; the cattle in northern China are primarily Bos Taurus, and the cattle in southern China are primarily Bos indicus. Cattle from southern China are characterized by a specific phenotype and adapted to the local environment. This study explored the genetic diversity, degree of admixture, and selection signature in eight local cattle breeds in southern China. The lowest level of heterozygosity was found in Hainan and Nandan cattle from Hainan and Guangxi province, respectively, whereas the highest level of heterozygosity was detected in Zhaotong cattle from Yunnan province. A neighbor-joining phylogenetic tree analysis clearly separated Lufeng cattle from other breeds, whereas Leiqiong and Hainan cattle have some crossover. Based on linkage disequilibrium-filtered single nucleotide polymorphisms (SNPs), the admixture analysis revealed two clusters corresponding to the taurine and indicine cattle lineages, and the local cattle breeds from southern China showed a certain degree of admixture. When K = 4 and 9, we found a slight separation among Leiqiong, Lufeng, and Hainan cattle. Meanwhile, we performed a selection signature analysis in Hainan, Leiqiong, and Lufeng cattle distributed in the extreme south of China, using the integrated haplotype score (iHS), Rsb statistic, and BayeScan software. Using the iHS approach, we identified 251, 270, and 256 candidate regions in Lufeng, Leiqiong, and Hainan cattle, respectively. Moreover, we identified 184, 174, and 146 candidate regions in pairwise comparisons of Leiqiong vs. Lufeng, Leiqiong vs. Hainan, and Hainan vs. Lufeng cattle using the Rsb approach. In addition, we identified 76 loci with a total of 48 genes under selection, based on the F_ST approach. Several candidate genes under selection were found to be related to meat quality, immunity, and adaptation to the local environment in southern China. Our results provide significant information about the genetic differences among the cattle breeds from southern China and the possible cause of difference in breed-specific characteristics. Selection signature analysis identified a few candidate SNPs and genes related to certain important traits of these cattle. In general, our results provide valuable insights into the genetic basis of specific traits under selection in certain local cattle breeds.

A novel missense mutation within KRT75 gene strongly affects heat stress in Chinese cattle

The KRT75 gene (Keratin 75) is a member of the type II epithelial α-keratin gene family which plays a key role in hair and nail formation. And the coat conformation affects heat tolarence in mammals. Therefore, the aim of this study was to identify novel single nucleotide polymorphisms (SNPs) of the KRT75 gene and further evaluate its relation to heat stress in Chinese cattle. A missense mutation (NC_037332.1: g.1052 T > C) of the bovine KRT75 was identified using the Bovine Genome Variation Database (BGVD). The g.1052 T > C variant was then genotyped in 519 individuals of 22 cattle breeds. Further analyses showed that the frequency of T allele in Chinese indigenous cattle breeds gradually diminished from northern groups to southern groups, whereas the frequency of C allele displayed a contrary patternl. Simultaneously, the frequency of the CC and CT genotype for southern groups was much greater than that of the TT genotype. Additionally, association analysis showed the genotypes were remarkably associated with mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) (P < 0.01). Our results demonstrated that the KRT75 gene might be a candidate gene associated with the heat stress in Chinese cattle.

A novel SNP of MYO1A gene associated with heat-tolerance in Chinese cattle

With the advent of global climate change, heat-tolerance is becoming more and more important to the sustainability of animal husbandry production systems. Previous studies have shown that MYO1A gene associated with pigmentation may be closely related to heat-tolerance in cattle. In this study, a novel missense mutation (NC_037332.1 g.56390345 A > G) was first detected in MYO1A in 891 individuals of 35 cattle breeds, which transformed the amino acid isoleucine into valine. The purpose of this study was to determine the allele frequencies distribution of this locus in Chinese indigenous cattle and to analyze the relationship between this locus and heat-tolerance. Further analysis showed that frequency of wild allele A decreased gradually from northern cattle to southern cattle, whereas frequency of mutant type allele G showed the opposite pattern, which was consistent with the distribution of various climatic conditions of China. Additionally, association analysis was carried out between genotypes and four climatic conditions (annual mean temperature (T), relative humidity (H), temperature-humidity index (THI) and average annual sunshine hours (100-cloudiness) (SR)). Analysis results showed that genotypes were significantly correlated with climatic conditions. Therefore, our results suggest that the novel SNP (rs209559414) is related to heat-tolerance trait of Chinese indigenous cattle.

The distribution of four missense mutations in SPTBN5 gene across native Chinese breeds

Modern cattle belong to two subspecies, Bos taurus and Bos indicus. Since divergence, cattle types have accumulated different genetic variations, which have contributed to highly differentiated phenotypes. The mammalian inner ear possesses functional and morphological innovations that contribute to its unique hearing capacities. The spectrin beta, non-erythrocytic 5 (SPTBN5) gene has been shown to play an important function in the inner ear. Four missense mutations: rs522333459 (c.7232G > C:p.Cys2411Ser), rs718838405 (c.6568A > C:p.Met2190Leu), rs516536785 (c.6283C > T:p.Leu2095Phe) and rs480278206 (c.4201T > C:p.Cys1401Arg) were identified in the bovine SPTBN5 gene by whole genome resequencing (http://animal.nwsuaf.edu.cn/code/index.php/BosVar), which might be candidate mutations related with hearing of both taurine and indicine cattle. In our study, PCR and DNA sequencing were used to explore the allele frequencies of four mutations of 971 individuals belonging to 38 native Chinese cattle breeds. We find that four mutant alleles showing strong geographic distribution, consisting with the ancestry distribution of taurine and indicine in China. In addition, we identified four mutations of SPTBN5 were diverged in taurine and indicine cattle showing signatures of adaptive evolution in two subspecies, which might participate in bovine inner ear development.

Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species - taurine cattle, zebu and yak - and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.

Distribution of the variant at AKIRIN2: c.*188G > A in Chinese cattle

Marbling score (MS), is an economically important trait in cattle. Previous results showed that a SNP (c.*188G > A) of akirin 2 (AKIRIN2) gene was associated with MS in Japanese Black cattle and Korean cattle. However, the distribution of the genotypic frequency of the single nucleotide polymorphism (SNP) has not been explored in Chinese cattle. In this study, we used polymerase chain reaction (PCR) and DNA sequencing to detect the variation in 1296 individuals from 39 Chinese cattle breeds, one semi-wild bovine species (Dulong) and three introduced breeds (Angus, Holstein and Brahman). Our study found the frequency of the A allele at this locus roughly diminished from north to south in Chinese cattle, and we detected statistically significant differences between Angus and Brahman (p < 0.05), Dulong and another two breeds (Angus and Holstein; p < 0.01) using Chi-Square Independence Test. Our results reflected the variation of AKIRIN2: c.*188G > A in Chinese cattle, which would help us better understand Chinese cattle genetic resources and provide reference for further research.

Exploring of InDel in bovine PSAP gene and their association with growth traits in different development stages

PSAP (prosaposin) is widely expressed in different organs, and plays an important role in fat deposit. Insertion/Deletion (InDel) is a relatively simple and effective DNA marker. However, the association of molecular marker at different stages of animal development has not received enough attention, especially fat deposition related traits. Therefore, eight cattle breeds were used to explore novel InDels variants within bovine PSAP gene, and to evaluate their effects on growth traits in different development stages. Herein, two novel InDels (P5:NC037355.1g.27974439-27974440 ins AGTGTGGTTAATGTCAAC and P8:NC037355.1g.27980734-27980752 del GTCAAAAAATCAGGGGAAAC) within the bovine PSAP gene were found, and their association with growth traits in different development stages were analyzed. Interestingly, the dominant genotype was different in different development stages both in NY cattle and JX cattle for daily gain and body weight. PSAP Gene expression patterns were analyzed in this study, high expression in the middle stage of adipocytes differentiation suggests that it plays a certain role in fat development. It reveals that InDels could affect phenotype in different development stages, which depend on the expression pattern of the host gene and their function in different tissues. These findings could provide a new way for molecular marker studies in bovine breeding and genetics.

A genetically unique Chinese cattle population shows evidence of common ancestry with wild species when analysed with a reduced ascertainment bias SNP panel

In Hong Kong, there is a cattle population of ~1,200 individuals of uncertain origin and genetic diversity. This population shows heterogeneous morphology, both in body type and pigmentation. Once used as draught animals by the local farmers, they were abandoned around the 1970s due to changes in the economy, and since then have lived as feral populations. To explore the origins of these cattle, we analysed ~50k genotype data of 21 Hong Kong feral cattle, along with data from 703 individuals of 36 cattle populations of European, African taurine, and Asian origin, the wild x domestic hybrid gayal, plus two wild bovine species, gaur and banteng. To reduce the effect of ascertainment bias ~4k loci that are polymorphic in the two wild species were selected for further analysis. The stringent SNP selection we applied resulted in increased heterozygosity across all populations studies, compared with the full panel of SNP, thus reducing the impact of ascertainment bias and facilitating the comparison of divergent breeds of cattle. Our results showed that Hong Kong feral cattle have relatively high levels of genetic distinctiveness, possibly due to the low level of artificial selection, and a likely common ancestry with wild species. We found signs of a putative taurine introgression, probably dating to the import of north European breeds during the British colonialism of Hong Kong. We showed that Hong Kong feral cattle, are distinct from Bos taurus and Bos indicus breeds. Our results highlight the distinctiveness of Hong Kong feral cattle and stress the conservation value of this indigenous breed that is likely to harbour adaptive genetic variation, which is a fundamental livestock resource in the face of climate change and diversifying market demands.

Effect direction meta-analysis of GWAS identifies extreme, prevalent and shared pleiotropy in a large mammal

In genome-wide association studies (GWAS), variants showing consistent effect directions across populations are considered as true discoveries. We model this information in an Effect Direction MEta-analysis (EDME) to quantify pleiotropy using GWAS of 34 Cholesky-decorrelated traits in 44,000+ cattle with sequence variants. The effect-direction agreement between independent bull and cow datasets was used to quantify the false discovery rate by effect direction (FDRed) and the number of affected traits for prioritised variants. Variants with multi-trait p < 1e–6 affected 1∼22 traits with an average of 10 traits. EDME assigns pleiotropic variants to each trait which informs the biology behind complex traits. New pleiotropic loci are identified, including signals from the cattle FTO locus mirroring its bystander effects on human obesity. When validated in the 1000-Bull Genome database, the prioritized pleiotropic variants consistently predicted expected phenotypic differences between dairy and beef cattle. EDME provides robust approaches to control GWAS FDR and quantify pleiotropy.

Xiang et al. developed an Effect Direction Meta-analysis (EDME) approach to identify true pleiotropy. They used Cholesky-transformation to decorrelate the traits and identified many pleiotropic variants that consistently predicted phenotypic differences in cattle.

DGAT1 K232A polymorphism is associated with milk production traits in Chinese cattle

The production traits of cattle, especially milk trait, are of great significance to human life. A quantitative trait loci (QTL) associated with milk fat content was detected in the centromeric region of cattle chromosome 14. This QTL harbors a strong candidate gene called DGAT1 responsible for the milk quality. A non-conservative substitution of lysine by alanine (K232A) was found in DGAT1 gene producing a strong effect on milk composition and yield. The lysine (K allele) is associated with increased milk fat content, while the decreased milk fat content is linked to the alanine (A allele) amino acid. To estimate the frequencies of the DGAT1 K232A polymorphism in Chinese cattle breeds, PCR and DNA sequencing methods were used to investigate the polymorphism of DGAT1 K232A in a total of 682 individuals, including 655 Chinese cattle and 27 Holstein cattle. The results demonstrated that the frequency of K allele gradually elevated from the northern group to the southern group of native Chinese cattle, whereas the frequency of A allele showed a contrary pattern, displaying a significant geographical difference across native Chinese cattle breeds. Our results confirm that the southern cattle group has higher milk fat content than that of the northern group.