Whole-genome resequencing reveals world-wide ancestry and adaptive introgression events of domesticated cattle in East Asia

Recent selection and introgression facilitated high-altitude adaptation in cattle

During the past 3000 years, cattle on the Qinghai-Xizang Plateau have developed adaptive phenotypes under the selective pressure of hypoxia, ultraviolet (UV) radiation, and extreme cold. The genetic mechanism underlying this rapid adaptation is not yet well understood. Here, we present whole-genome resequencing data for 258 cattle from 32 cattle breeds/populations, including 89 Tibetan cattle representing eight populations distributed at altitudes ranging from 3400 m to 4300 m. Our genomic analysis revealed that Tibetan cattle exhibited a continuous phylogeographic cline from the East Asian taurine to the South Asian indicine ancestries. We found that recently selected genes in Tibetan cattle were related to body size (HMGA2 and NCAPG) and energy expenditure (DUOXA2). We identified signals of sympatric introgression from yak into Tibetan cattle at different altitudes, covering 0.64%-3.26% of their genomes, which included introgressed genes responsible for hypoxia response (EGLN1), cold adaptation (LRP11), DNA damage repair (LATS1), and UV radiation resistance (GNPAT). We observed that introgressed yak alleles were associated with noncoding variants, including those in present EGLN1. In Tibetan cattle, three yak introgressed SNPs in the EGLN1 promoter region reduced the expression of EGLN1, suggesting that these genomic variants enhance hypoxia tolerance. Taken together, our results indicated complex adaptation processes in Tibetan cattle, where recently selected genes and introgressed yak alleles jointly facilitated rapid adaptation to high-altitude environments.

Evolution and legacy of East Asian aurochs

Aurochs (Bos primigenius), once widely distributed in Afro-Eurasia, became extinct in the early 1600 s. However, their phylogeography and relative contributions to domestic cattle remain unknown. In this study, we analyzed 16 genomes of ancient aurochs and three mitogenomes of ancient bison (Bison priscus) excavated in East Asia, dating from 43,000 to 3,590 years ago. These newly generated data with previously published genomic information on aurochs as well as ancient/extant domestic cattle worldwide through genome analysis. Our findings revealed significant genetic divergence between East Asian aurochs and their European, Near Eastern, and African counterparts on the basis of both mitochondrial and nuclear genomic data. Furthermore, we identified evidence of gene flow from East Asian aurochs into ancient and present-day taurine cattle, suggesting their potential role in facilitating the environmental adaptation of domestic cattle.

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Genome-wide local ancestry and the functional consequences of admixture in African and European cattle populations

Bos taurus (taurine) and Bos indicus (indicine) cattle diverged at least 150,000 years ago and, since that time, substantial genomic differences have evolved between the two lineages. During the last two millennia, genetic exchange in Africa has resulted in a complex tapestry of taurine-indicine ancestry, with most cattle populations exhibiting varying levels of admixture. Similarly, there are several Southern European cattle populations that also show evidence for historical gene flow from indicine cattle, the highest levels of which are found in the Central Italian White breeds. Here we use two different software tools (MOSAIC and ELAI) for local ancestry inference (LAI) with genome-wide high- and low-density SNP array data sets in hybrid African and residually admixed Southern European cattle populations and obtained broadly similar results despite critical differences in the two LAI methodologies used. Our analyses identified genomic regions with elevated levels of retained or introgressed ancestry from the African taurine, European taurine, and Asian indicine lineages. Functional enrichment of genes underlying these ancestry peaks highlighted biological processes relating to immunobiology and olfaction, some of which may relate to differing susceptibilities to infectious diseases, including bovine tuberculosis, East Coast fever, and tropical theileriosis. Notably, for retained African taurine ancestry in admixed trypanotolerant cattle we observed enrichment of genes associated with haemoglobin and oxygen transport. This may reflect positive selection of genomic variants that enhance control of severe anaemia, a debilitating feature of trypanosomiasis disease, which severely constrains cattle agriculture across much of sub-Saharan Africa.

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.

The genomic natural history of the aurochs

Now extinct, the aurochs (Bos primigenius) was a keystone species in prehistoric Eurasian and North African ecosystems, and the progenitor of cattle (Bos taurus), domesticates that have provided people with food and labour for millennia1. Here we analysed 38 ancient genomes and found 4 distinct population ancestries in the aurochs-European, Southwest Asian, North Asian and South Asian-each of which has dynamic trajectories that have responded to changes in climate and human influence. Similarly to Homo heidelbergensis, aurochsen first entered Europe around 650 thousand years ago2, but early populations left only trace ancestry, with both North Asian and European B. primigenius genomes coalescing during the most recent glaciation. North Asian and European populations then appear separated until mixing after the climate amelioration of the early Holocene. European aurochsen endured the more severe bottleneck during the Last Glacial Maximum, retreating to southern refugia before recolonizing from Iberia. Domestication involved the capture of a small number of individuals from the Southwest Asian aurochs population, followed by early and pervasive male-mediated admixture involving each ancestral strain of aurochs after domestic stocks dispersed beyond their cradle of origin.

Deep Ancestral Introgressions between Ovine Species Shape Sheep Genomes via Argali-Mediated Gene Flow

Previous studies revealed extensive genetic introgression between Ovis species, which affects genetic adaptation and morphological traits. However, the exact evolutionary scenarios underlying the hybridization between sheep and allopatric wild relatives remain unknown. To address this problem, we here integrate the reference genomes of several ovine and caprine species: domestic sheep, argali, bighorn sheep, snow sheep, and domestic goats. Additionally, we use 856 whole genomes representing 169 domestic sheep populations and their six wild relatives: Asiatic mouflon, urial, argali, snow sheep, thinhorn sheep, and bighorn sheep. We implement a comprehensive set of analyses to test introgression among these species. We infer that the argali lineage originated ∼3.08 to 3.35 Mya and hybridized with the ancestor of Pachyceriforms (e.g. bighorn sheep and snow sheep) at ∼1.56 Mya. Previous studies showed apparent introgression from North American Pachyceriforms into the Bashibai sheep, a Chinese native sheep breed, despite of their wide geographic separation. We show here that, in fact, the apparent introgression from the Pachyceriforms into Bashibai can be explained by the old introgression from Pachyceriforms into argali and subsequent recent introgression from argali into Bashibai. Our results illustrate the challenges of estimating complex introgression histories and provide an example of how indirect and direct introgression can be distinguished.

[Runs of Homozygosity Decipher Genetic Diversity in Cattle Breed Dwelling in the Colder Regions of the World]

BACKGROUND

Our study focuses on Yakutian cattle, a Siberian native breed, examining its inbreeding and diversity through genome-wide analysis of runs of homozygosity (ROHs). Yakutian cattle are adapted to Siberia's harsh sub-arctic conditions, enduring temperatures below -70°C. However, the population genetics studies on this breed are scanty, to document the genetic uniqueness in these cattle.

RESULTS

We analyzed 40 Yakutian cattle with strict quality control for ROH detection yielding 683 homozygous segments, averaging 17 per individual with an average length of 9 Mb. ROH regions were found to be involved in important pathways pertaining to cold adaptation. Autozygosity ranged from 1% to 12% of the genome, with a relatively low average inbreeding coefficient (FROH) of 0.057, as compared to other breeds. Also, the different diversity indicators, namely, principal component analysis, heterozygosity, and effective population size analysis, revealed the prevalence of genetic diversity within the breed.

CONCLUSION

Our findings on ROH are the first of its kind in Yakutian cattle that support their adaptability to colder environments, as evidenced by low inbreeding and high genetic diversity.

Development and validation of a 5K low-density SNP chip for Hainan cattle

BACKGROUND

This study aimed to design and develop a 5K low-density liquid chip for Hainan cattle utilizing targeted capture sequencing technology. The chip incorporates a substantial number of functional single nucleotide polymorphism (SNP) loci derived from public literature, including SNP loci significantly associated with immunity, heat stress, meat quality, reproduction, and other traits. Additionally, SNPs located in the coding regions of immune-related genes from the Bovine Genome Variation Database (BGVD) and Hainan cattle-specific SNP loci were included.

RESULTS

A total of 5,293 SNPs were selected, resulting in 9,837 DNA probes with a coverage rate of 85.69%, thereby creating a Hainan cattle-specific 5K Genotyping by Target Sequencing (GBTS) liquid chip. Evaluation with 152 cattle samples demonstrated excellent clustering performance and a detection rate ranging from 96.60 to 99.07%, with 94.5% of SNP sites exhibiting polymorphism. The chip achieved 100% gender coverage and displayed a heterozygosity rate between 14.20% and 29.65%, with a repeatability rate of 99.65-99.85%. Analyses using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed the potential regulatory roles of exonic SNPs in immune response pathways.

CONCLUSION

The development and validation of the 5K GBTS liquid chip for Hainan cattle represent a valuable tool for genome analysis and genetic diversity assessment. Furthermore, it facilitates breed identification, gender determination, and kinship analysis, providing a foundation for the efficient utilization and development of local cattle genetic resources.

Introgression drives adaptation to the plateau environment in a subterranean rodent

BACKGROUND

Introgression has repeatedly been shown to play an important role in the adaptation of species to extreme environments, yet how introgression enables rodents with specialized subterranean lifestyle to acclimatize to high altitudes is still unclear. Myospalacinae is a group of subterranean rodents, among which the high-altitude plateau zokors (Eospalax baileyi) and the low-altitude Gansu zokors (E. cansus) are sympatrically distributed in the grassland ecosystems of the Qinghai-Tibet Plateau (QTP). Together, they provide a model for the study of the role of introgression in the adaptation of low-altitude subterranean rodents to high altitudes.

RESULTS

Applying low-coverage whole-genome resequencing and population genetics analyses, we identified evidence of adaptive introgression from plateau zokors into Gansu zokors, which likely facilitated the adaptation of the latter to the high-altitude environment of the QTP. We identified positively selected genes with functions related to energy metabolism, cardiovascular system development, calcium ion transport, and response to hypoxia which likely made critical contributions to adaptation to the plateau environment in both plateau zokors and high-altitude populations of Gansu zokors.

CONCLUSIONS

Introgression of genes associated with hypoxia adaptation from plateau zokors may have played a role in the adaptation of Gansu zokors to the plateau environment. Our study provides new insights into the understanding of adaptive evolution of species on the QTP and the importance of introgression in the adaptation of species to high-altitude environments.

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing

Structural variants play an important role in evolutionary processes. Besides, they constitute a large source of inter individual genetic variation that might represent a major factor in the aetiology of complex, multifactorial traits. Their importance in adaptation is becoming increasingly evident in literature. Yet, the characterization of the genomic landscape of structural variants in local breeds remains scarce to date. Herein, we investigate patterns and gene annotation of structural variants in the Creole cattle from Guadeloupe breed using whole genome sequences from 23 bulls representative of the population. In total, we detected 32821 ascertained SV defining 15258 regions, representing ~ 17% of the Creole cattle genome. Among these, 6639 regions have not been previously reported in the Database of Genomic Variants archive. Average number of structural variants detected per individual in the studied population is in the same order of magnitude of that observed in indicine populations and higher than that reported in taurine breeds. We observe an important within-individual variability where approximately half of the detected structural variants have low frequency (MAF < 0.25). Most of the detected structural variants (55%) occurred in intergenic regions. Genic structural variants overlapped with 7793 genes and the predicted effect of most of them is ranked as "modifier". Among the structural variants that were predicted to have a high functional impact on the protein, a 5.5 Kb in length, highly frequent deletion on chromosome 2, affects ALPI, a gene associated with the interaction between gut microbiota and host immune system. The 6639 newly identified structural variants regions include three deletions and three duplications shared by more than 80% of individuals that are significantly enriched for genes related to tRNA threonylcarbamoyladenosine metabolic process, important for temperature adaptation in thermophilic organisms, therefore suggesting a potential role in the thermotolerance of Creole cattle from Guadeloupe cattle to tropical climate. Overall, highly frequent structural variants that are specific to the Creole cattle population encompass olfactory receptor and immunity genes as well as genes involved in muscle tone, muscle development and contraction. Beyond mapping and characterizing structural variants in the Creole cattle from Guadeloupe breed, this study provides valuable information for a better understanding of the potential role of chromosomal rearrangements in adaptive traits in cattle.

A Whole-Genome Scan Revealed Genomic Features and Selection Footprints of Mengshan Cattle

(1) Background: Mengshan cattle from the Yimeng mountainous region in China stand out as a unique genetic resource, known for their adaptive traits and environmental resilience. However, these cattle are currently endangered and comprehensive genomic characterization remains largely unexplored. This study aims to address this gap by investigating the genomic features and selection signals in Mengshan cattle. (2) Methods: Utilizing whole-genome resequencing data from 122 cattle, including 37 newly sequenced Mengshan cattle, we investigated population structure, genetic diversity, and selection signals. (3) Results: Our analyses revealed that current Mengshan cattle primarily exhibit European taurine cattle ancestry, with distinct genetic characteristics indicative of adaptive traits. We identified candidate genes associated with immune response, growth traits, meat quality, and neurodevelopment, shedding light on the genomic features underlying the unique attributes of Mengshan cattle. Enrichment analysis highlighted pathways related to insulin secretion, calcium signaling, and dopamine synapse, further elucidating the genetic basis of their phenotypic traits. (4) Conclusions: Our results provide valuable insights for further research and conservation efforts aimed at preserving this endangered genetic resource. This study enhances the understanding of population genetics and underscores the importance of genomic research in informing genetic resources and conservation initiatives for indigenous cattle breeds.

Genome-Wide Analysis of Genetic Diversity and Selection Signatures in Zaobei Beef Cattle

This investigation provides a comprehensive analysis of genomic diversity and selection signatures in Zaobei beef cattle, an indigenous breed known for its adaptation to hot and humid climates and superior meat quality. Whole-genome resequencing was conducted on 23 Zaobei cattle, compared with 46 Simmental cattle to highlight genetic distinctions. Population structure analysis confirmed the genetic uniqueness of Zaobei cattle. Using methods such as DASDC v1.01, XPEHH, and θπ ratio, we identified 230, 232, and 221 genes through DASDC, including hard sweeps, soft sweeps, and linkage sweeps, respectively. Coincidentally, 109 genes were identified when using XPEHH and θπ ratio methods. Together, these analyses revealed eight positive selection genes (ARHGAP15, ZNF618, USH2A, PDZRN4, SPATA6, ROR2, KCNIP3, and VWA3B), which are linked to critical traits such as heat stress adaptation, fertility, and meat quality. Moreover, functional enrichment analyses showed pathways related to autophagy, immune response, energy metabolism, and muscle development. The comprehensive genomic insights gained from this study provide valuable knowledge for breeding programs aimed at enhancing the beneficial traits in Zaobei cattle.

Whole-genome resequencing deciphers patterns of genetic diversity, phylogeny, and evolutionary dynamics in Kashmir cattle

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Illuminating Genetic Diversity and Selection Signatures in Matou Goats through Whole-Genome Sequencing Analysis

(1) Background: Matou goats, native to Hunan and Hubei provinces in China, are renowned for their exceptional meat and skin quality. However, a comprehensive whole-genome-based exploration of the genetic architecture of this breed is scant in the literature. (2) Methods: To address this substantial gap, we used whole-genome sequences of 20 Matou goats and compared them with published genomic data of 133 goats of different breeds across China. This comprehensive investigation sought to assess genetic diversity, population structure, and the presence of genomic selection signals. (3) Results: The whole genome of Matou goat populations yielded a substantial catalog of over 19 million single nucleotide polymorphisms (SNPs), primarily distributed within intergenic and intron regions. The phylogenetic tree analysis revealed distinct clades corresponding to each goat population within the dataset. Notably, this analysis positioned Matou goats in a closer genetic affinity with Guizhou White goats, compared to other recognized goat breeds. This observation was corroborated by principal component analysis (PCA) and admixture analysis. Remarkably, Matou goats exhibited diminished genetic diversity and a notable degree of inbreeding, signifying a reduced effective population size. Moreover, the study employed five selective sweep detection methods (including PI, CLR, PI-Ratio, Fst, and XP-EHH) to screen top signal genes associated with critical biological functions, encompassing cardiomyocytes, immunity, coat color, and meat quality. (4) Conclusions: In conclusion, this study significantly advances our understanding of the current genetic landscape and evolutionary dynamics of Matou goats. These findings underscore the importance of concerted efforts in resource conservation and genetic enhancement for this invaluable breed.

Genome-Wide Identification and Evolutionary and Mutational Analysis of the Bos taurus Pax Gene Family

Bos taurus is known for its tolerance of coarse grains, adaptability, high temperature, humidity, and disease resistance. Primarily, cattle are raised for their meat and milk, and pinpointing genes associated with traits relevant to meat production can enhance their overall productivity. The aim of this study was to identify the genome, analyze the evolution, and explore the function of the Pax gene family in B. taurus to provide a new molecular target for breeding in meat-quality-trait cattle. In this study, 44 Pax genes were identified from the genome database of five species using bioinformatics technology, indicating that the genetic relationships of bovids were similar. The Pax3 and Pax7 protein sequences of the five animals were highly consistent. In general, the Pax gene of the buffalo corresponds to the domestic cattle. In summary, there are differences in affinity between the Pax family genes of buffalo and domestic cattle in the Pax1/9, Pax2/5/8, Pax3/7, and Pax4/6 subfamilies. We believe that Pax1/9 has an effect on the growth traits of buffalo and domestic cattle. The Pax3/7 gene is conserved in the evolution of buffalo and domestic animals and may be a key gene regulating the growth of B. taurus. The Pax2/5/8 subfamily affects coat color, reproductive performance, and milk production performance in cattle. The Pax4/6 subfamily had an effect on the milk fat percentage of B. taurus. The results provide a theoretical basis for understanding the evolutionary, structural, and functional characteristics of the Pax family members of B. taurus and for molecular genetics and the breeding of meat-production B. taurus species.

Multi-omic Analyses Shed Light on The Genetic Control of High-altitude Adaptation in Sheep

Sheep were domesticated in the Fertile Crescent and then spread globally, where they have been encountering various environmental conditions. The Tibetan sheep has adapted to high altitudes on the Qinghai-Tibet Plateau over the past 3000 years. To explore genomic variants associated with high-altitude adaptation in Tibetan sheep, we analyzed Illumina short-reads of 994 whole genomes representing ∼ 60 sheep breeds/populations at varied altitudes, PacBio High fidelity (HiFi) reads of 13 breeds, and 96 transcriptomes from 12 sheep organs. Association testing between the inhabited altitudes and 34,298,967 variants was conducted to investigate the genetic architecture of altitude adaptation. Highly accurate HiFi reads were used to complement the current ovine reference assembly at the most significantly associated β-globin locus and to validate the presence of two haplotypes A and B among 13 sheep breeds. The haplotype A carried two homologous gene clusters: (1) HBE1, HBE2, HBB-like, and HBBC, and (2) HBE1-like, HBE2-like, HBB-like, and HBB; while the haplotype B lacked the first cluster. The high-altitude sheep showed highly frequent or nearly fixed haplotype A, while the low-altitude sheep dominated by haplotype B. We further demonstrated that sheep with haplotype A had an increased hemoglobin-O2 affinity compared with those carrying haplotype B. Another highly associated genomic region contained the EGLN1 gene which showed varied expression between high-altitude and low-altitude sheep. Our results provide evidence that the rapid adaptive evolution of advantageous alleles play an important role in facilitating the environmental adaptation of Tibetan sheep.

Genome-wide association study provided insights into the polled phenotype and polled intersex syndrome (PIS) in goats

BACKGROUND

Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals.

METHODS

The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats.

RESULTS

The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats.

CONCLUSION

We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

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.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle

BACKGROUND

The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data.

RESULT

The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1).

CONCLUSION

We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Unraveling genomic diversity and positive selection signatures of Qaidam cattle through whole-genome re-sequencing

Qaidam cattle are a typical Chinese native breed inhabiting northwest China. They bear the characteristics of high cold and roughage tolerance, low-oxygen adaptability and good meat quality. To analyze the genetic diversity of Qaidam cattle, 60 samples were sequenced using whole-genome resequencing technology, along with 192 published sets of whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, North Chinese cattle breeds, East Asian taurine cattle, Eurasian taurine cattle and European taurine cattle as controls. It was found that Qaidam cattle have rich genetic diversity in Bos taurus, but the degree of inbreeding is also high, which needs further protection. The phylogenetic analysis, principal component analysis and ancestral component analysis showed that Qaidam cattle mainly originated from East Asian taurine cattle. Qaidam cattle had a closer genetic relationship with the North Chinese cattle breeds and the least differentiation from Mongolian cattle. Annotating the selection signals obtained by composite likelihood ratio, nucleotide diversity analysis, integrated haplotype score, genetic differentiation index, genetic diversity ratio and cross-population extended haplotype homozygosity methods, several genes associated with immunity, reproduction, meat, milk, growth and adaptation showed strong selection signals. In general, this study provides genetic evidence for understanding the germplasm characteristics of Qaidam cattle. At the same time, it lays a foundation for the scientific and reasonable protection and utilization of genetic resources of Chinese local cattle breeds, which has great theoretical and practical significance.

Whole-genome resequencing uncovers diversity and selective sweep in Kazakh cattle

The Kazakh cattle in the Xinjiang Uygur Autonomous Region of China are highly adaptable and have multiple uses, including milk and meat production, and use as draft animals. They are an excellent original breed that could be enhanced by breeding and hybrid improvement. However, the genomic diversity and signature of selection underlying the germplasm characteristics require further elucidation. Herein, we evaluated 26 Kazakh cattle genomes in comparison with 103 genomes of seven other cattle breeds from regions around the world to assess the Kazakh cattle genetic variability. We revealed that the relatively low linkage disequilibrium at large SNP distances was strongly correlated with the largest effective population size among Kazakh cattle. Using population structural analysis, we next demonstrated a taurine lineage with restricted Bos indicus introgression among Kazakh cattle. Notably, we identified putative selected genes associated with resistance to disease and body size within Kazakh cattle. Together, our findings shed light on the evolutionary history and breeding profile of Kazakh cattle, as well as offering indispensable resources for germplasm resource conservation and crossbreeding program implementation.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

POLB Regulates Proliferation and Apoptosis of Bovine Primary Myocytes

DNA polymerase β (DNA polymerase beta (POLB)) belongs to a member of the DNA polymerase X family, mainly involved in various biological metabolic processes, such as eukaryotic DNA replication, DNA damage repair, gene recombination, and cell cycle regulation. In this study, the muscle development-related gene POLB was screened by selection signature and RNA-seq analysis and then validated for the proliferation and apoptosis of bovine primary myocytes. It was also found that overexpression of the POLB gene had a pro-apoptosis effect, but interfering with the expression of the gene had no significant effect on cells. Then, the analysis of related apoptotic genes revealed that POLB overexpression affected CASP9 gene expression.

Research Progress and Applications of Bovine Genome in the Tribe Bovini

Various bovine species have been domesticated and bred for thousands of years, and they provide adequate animal-derived products, including meat, milk, and leather, to meet human requirements. Despite the review studies on economic traits in cattle, the genetic basis of traits has only been partially explained by phenotype and pedigree breeding methods, due to the complexity of genomic regulation during animal development and growth. With the advent of next-generation sequencing technology, genomics projects, such as the 1000 Bull Genomes Project, Functional Annotation of Animal Genomes project, and Bovine Pangenome Consortium, have advanced bovine genomic research. These large-scale genomics projects gave us a comprehensive concept, technology, and public resources. In this review, we summarize the genomics research progress of the main bovine species during the past decade, including cattle (Bos taurus), yak (Bos grunniens), water buffalo (Bubalus bubalis), zebu (Bos indicus), and gayal (Bos frontalis). We mainly discuss the development of genome sequencing and functional annotation, focusing on how genomic analysis reveals genetic variation and its impact on phenotypes in several bovine species.

Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data

OBJECTIVE

The objective of this study was to investigate the genetic diversity, population structure and whole-genome selection signatures of Luxi cattle to reveal its genomic characteristics in terms of meat and carcass traits, skeletal muscle development, body size, and other traits.

METHODS

To further analyze the genomic characteristics of Luxi cattle, this study sequenced the whole-genome of 16 individuals from the core conservation farm in Shandong region, and collected 174 published genomes of cattle for conjoint analysis. Furthermore, three different statistics (pi, Fst, and XP-EHH) were used to detect potential positive selection signatures related to selection in Luxi cattle. Moreover, gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were performed to reveal the potential biological function of candidate genes harbored in selected regions.

RESULTS

The results showed that Luxi cattle had high genomic diversity and low inbreeding levels. Using three complementary methods (pi, Fst, and XP-EHH) to detect the signatures of selection in the Luxi cattle genome, there were 2,941, 2,221 and 1,304 potentially selected genes identified, respectively. Furthermore, there were 45 genes annotated in common overlapping genomic regions covered 0.723 Mb, including PLAG1 zinc finger (PLAG1), dedicator of cytokinesis 3 (DOCK3), ephrin A2 (EFNA2), DAZ associated protein 1 (DAZAP1), Ral GTPase activating protein catalytic subunit alpha 1 (RALGAPA1), mediator complex subunit 13 (MED13), and decaprenyl diphosphate synthase subunit 2 (PDSS2), most of which were enriched in pathways related to muscle growth and differentiation and immunity.

CONCLUSION

In this study, we provided a series of genes associated with important economic traits were found in positive selection regions, and a scientific basis for the scientific conservation and genetic improvement of Luxi cattle.

Genomic insights into the population history and adaptive traits of Latin American Criollo cattle

Criollo cattle, the descendants of animals brought by Iberian colonists to the Americas, have been the subject of natural and human-mediated selection in novel tropical agroecological zones for centuries. Consequently, these breeds have evolved distinct characteristics such as resistance to diseases and exceptional heat tolerance. In addition to European taurine (Bos taurus) ancestry, it has been proposed that gene flow from African taurine and Asian indicine (Bos indicus) cattle has shaped the ancestry of Criollo cattle. In this study, we analysed Criollo breeds from Colombia and Venezuela using whole-genome sequencing (WGS) and single-nucleotide polymorphism (SNP) array data to examine population structure and admixture at high resolution. Analysis of genetic structure and ancestry components provided evidence for African taurine and Asian indicine admixture in Criollo cattle. In addition, using WGS data, we detected selection signatures associated with a myriad of adaptive traits, revealing genes linked to thermotolerance, reproduction, fertility, immunity and distinct coat and skin coloration traits. This study underscores the remarkable adaptability of Criollo cattle and highlights the genetic richness and potential of these breeds in the face of climate change, habitat flux and disease challenges. Further research is warranted to leverage these findings for more effective and sustainable cattle breeding programmes.

Whole-genome resequencing provides insights into the diversity and adaptation to desert environment in Xinjiang Mongolian cattle

BACKGROUND

Xinjiang Mongolian cattle is an indigenous breed that inhabits the Taklimakan Desert and is characterized by its small body size. However, the genomic diversity, origin, and genetic basis underlying the adaptation to the desert environment have been poorly studied.

RESULTS

We analyzed patterns of Xinjiang Mongolian cattle genetic variation by sequencing 20 genomes together with seven previously sequenced genomes and comparing them to the 134 genomes of nine representative breeds worldwide. Among the breeds of Bos taurus, we found the highest nucleotide diversity (0.0024) associated with the lower inbreeding coefficient (2.0110-6), the lowest linkage disequilibrium (r2 = 0.3889 at distance of 10 kb), and the highest effective population size (181 at 20 generations ago) in Xinjiang Mongolian cattle. The genomic diversity pattern could be explained by a limited introgression of Bos indicus genes. More importantly, similarly to desert-adapted camel and same-habitat sheep, we also identified signatures of selection including genes, GO terms, and/or KEGG pathways controlling water reabsorption and osmoregulation, metabolic regulation and energy balance, as well as small body size in Xinjiang Mongolian cattle.

CONCLUSIONS

Our results imply that Xinjiang Mongolian cattle might have acquired distinct genomic diversity by virtue of the introgression of Bos indicus, which helps understand the demographic history. The identification of selection signatures can provide novel insights into the genomic basis underlying the adaptation of Xinjiang Mongolian cattle to the desert environment.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Genome-wide survey reveals the genetic background of Xinjiang Brown cattle in China

Introduction: Xinjiang Brown cattle are a famous dual-purpose (dairy-beef) cultivated breed in China that occupy a pivotal position within the cattle breeding industry in Xinjiang, China. However, little information is available on the genetic background of this breed. To fill this research gap, we conducted a whole-genome screen using specific-locus amplified fragment sequencing to examine the genetic structure and diversity of 130 Xinjiang Brown cattle-grazing type (XBG, traditional type) cattle. Methods: A subsequent joint analysis incorporating two ancestral breeds, specifically 19 Brown Swiss (BS) foreign and nine Kazakh (KZ) Chinese cattle, as well as 20 Xinjiang Brown cattle-housing type (XBH) cattle, was used to explore the genetic background of the Xinjiang Brown cattle. Results: The results showed that, after nearly a century of crossbreeding, XBG cattle formed a single population with a stable genetic performance. The genetic structure, genetic diversity, and selection signature analysis of the two ancestral types showed highly different results compared to that of XBH cattle. Local ancestry inference showed that the average proportions of XGB cattle within the BS and KZ cattle lineages were 37.22% and 62.78%, respectively, whereas the average proportions of XBH cattle within the BS and KZ cattle lineages were 95.14% and 4.86%, respectively. Thus, XGB cattle are more representative of all Xinjiang Brown cattle, in line with their breeding history, which involves crossbreeding. Two complementary approaches, fixation index and mean nucleotide diversity, were used to detect selection signals in the four aforementioned cattle breeds. Finally, the analysis of 26 candidate genes in Xinjiang Brown cattle revealed significant enrichment in 19 Gene Ontology terms, and seven candidate genes were enriched in three pathways related to disease resistance (CDH4, SIRPB1, and SIRPα) and the endocrine system (ADCY5, ABCC8, KCNJ11, and KCNMA1). Finally, development of the core SNPs in XBG cattle yielded 8,379 loci. Conclusion: The results of this study detail the evolutionary process of crossbreeding in Xinjiang Brown cattle and provide guidance for selecting and breeding new strains of this species.

Hypoxia Inducible Factor pathway proteins in high-altitude mammals

Humans and other mammals inhabit hypoxic high-altitude locales. In many of these species, genes under positive selection include ones in the Hypoxia Inducible Factor (HIF) pathway. One is PHD2 (EGLN1), which encodes for a key oxygen sensor. Another is HIF2A (EPAS1), which encodes for a PHD2-regulated transcription factor. Recent studies have provided insights into mechanisms for these high-altitude alleles. These studies have (i) shown that selection can occur on nonconserved, unstructured regions of proteins, (ii) revealed that high altitude-associated amino acid substitutions can have differential effects on protein-protein interactions, (iii) provided evidence for convergent evolution by different molecular mechanisms, and (iv) suggested that mutations in different genes can complement one another to produce a set of adaptive phenotypes.

Evidence for early domestic yak, taurine cattle, and their hybrids on the Tibetan Plateau

Domestic yak, cattle, and their hybrids are fundamental to herder survival at high altitudes on the Tibetan Plateau. However, little is known about their history. Bos remains are uncommon in this region, and ancient domestic yak have not been securely identified. To identify Bos taxa and investigate their initial management, we conducted zooarchaeological analyses of 193 Bos specimens and sequenced five nuclear genomes from recently excavated assemblages at Bangga. Morphological data indicated that more cattle than yak were present. Ancient mitochondrial DNA and nuclear genome sequences identified taurine cattle and provided evidence for domestic yak and yak-cattle hybridization ~2500 years ago. Reliance on diverse Bos species and their hybrid has increased cattle adaptation and herder resilience to plateau conditions. Ancient cattle and yak at Bangga were closely related to contemporary livestock, indicating early herder legacies and the continuity of cattle and yak husbandry on the Tibetan Plateau.

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.

Inference of Admixture Origins in Indigenous African Cattle

Present-day African cattle retain a unique genetic profile composed of a mixture of the Bos taurus and Bos indicus populations introduced into the continent at different time periods. However, details of the admixture history and the exact origins of the source populations remain obscure. Here, we infer the source of admixture in the earliest domestic cattle in Africa, African taurine. We detect a significant contribution (up to ∼20%) from a basal taurine lineage, which might represent the now-extinct African aurochs. In addition, we show that the indicine ancestry of African cattle, although most closely related to so-far sampled North Indian indicine breeds, has a small amount of additional genetic affinity to Southeast Asian indicine breeds. Our findings support the hypothesis of aurochs introgression into African taurine and generate a novel hypothesis that the origin of indicine ancestry in Africa might be different indicine populations than the ones found in North India today.

The genetic secrets of adaptation: decoding the significance of the 30-bp insertion in the KRT77 gene for Chinese cattle

The KRT77 gene is a type II epithelial cell α-keratin gene family member that plays a crucial role in animal epidermal and coat formation. This study aimed to investigate the relationship between the KRT77 gene and the adaptability of Chinese cattle in varying environments by exploring the distribution of an exon insertion of the KRT77 gene in different cattle populations. Our analysis involved amplifying and sequencing DNA samples from 362 individuals from 24 cattle breeds in China. Our findings reveal a gradual increase in the frequency of insertion from the northwest to the southeast population. We conducted an association analysis between the genotypes and climate data, revealing a correlation between the insertion and local annual mean temperature, relative humidity, and temperature humidity index. The study highlights the significance of the newly identified KRT77 gene insertion as a variation associated with environmental adaptation in Chinese cattle.This insertion variation increased insights into the genetic mechanisms that drive adaptation in Chinese cattle, emphasizing the importance of the 30-bp insertion in the KRT77 gene. Our findings facilitate further research to improve cattle breeding strategies for adaptability to changing environments from the northwest to the southeast population. In conclusion, this study provides value.

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.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

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.

Assessing genomic diversity and selective pressures in Bashan cattle by whole-genome sequencing data

Specific ecological environments and domestication have continuously influenced the physiological characteristics of Chinese indigenous cattle. Among them, Bashan cattle belongs to one of the indigenous breeds. However, the genomic diversity of Bashan cattle is still unknown. Published whole-genome sequencing (WGS) data of 13 Bashan cattle and 48 worldwide cattle were used to investigate the genetic composition and selection characteristics of Bashan cattle. The population structure analysis revealed that Bashan cattle harbored ancestries with East Asian taurine and Chinese indicine. Genetic diversity analysis implied the relatively high genomic diversity in Bashan cattle. Through the identification of containing >5 nsSNPs or frameshift mutations genes in Bashan cattle, a large number of pathways related to sensory perception were discovered. CLR, θπ ratio, FST, and XP-EHH methods were used to detect the candidate signatures of positive selection in Bashan cattle. Among the identified genes, most of the enriched signal pathways were related to environmental information processing, biological systems, and metabolism. We mainly reported genes related to the nervous system (HCN1, KATNA1, FSTL1, GRIK2, and CPLX2), immune (CD244, SLAMF1, LY9, and CD48), and reproduction (AKR1C1, AKR1C3, AKR1C4, and TUSC3). Our findings will be significant in understanding the molecular basis underlying phenotypic variation of breed-related traits and improving productivity in Bashan cattle.

Selection signatures of Qinchuan cattle based on whole-genome sequences

Qinchuan cattle has gradually improved in body shape and growth rate in the long-term breeding process from the draft cattle to beef cattle. As the head of the five local yellow cattle in China, the Qinchuan cattle has been designated as a specialized beef cattle breed. We investigated the selection signatures using whole genome sequencing data in Qinchuan cattle. Based on Fst, we detected hundreds of candidate genes under selection across Qinchuan, Red Angus, and Japanese Black cattle. Through protein-protein interaction analysis and functional annotation of candidate genes, the results revealed that KMT2E, LTBP1 and NIPBL were related to brain size, body characteristics, and limb development, respectively, suggesting that these potential genes may affect the growth and development traits in Qinchuan cattle. ARIH2, DACT1 and DNM2, et al. are related to meat quality. Meanwhile, TBXA2R can be used as a gene associated with reproductive function, and USH2A affect coat color. This provided a glimpse into the formation of breeds and molecular genetic breeding. Our findings will promote genome-assisted breeding to improve animal production and health.

Genome-wide identification and analysis of TMT-based proteomes in longissimus dorsi tissue from Kazakh cattle and Xinjiang brown cattle

With the gradual completion of the human genome project, proteomes have gained extremely important value in the fields of human disease and biological process research. In our previous research, we performed transcriptomic analyses of longissimus dorsi tissue from Kazakh cattle and Xinjiang brown cattle and conducted in-depth studies on the muscles of both species through epigenetics. However, it is unclear whether differentially expressed proteins in Kazakh cattle and Xinjiang brown cattle regulate muscle production and development. In this study, a proteomic analysis was performed on Xinjiang brown cattle and Kazakh cattle by using TMT markers, HPLC classification, LC/MS and bioinformatics analysis. A total of 13,078 peptides were identified, including 11,258 unique peptides. We identified a total of 1874 proteins, among which 1565 were quantifiable. Compared to Kazakh cattle, Xinjiang brown cattle exhibited 75 upregulated proteins and 44 downregulated proteins. These differentially expressed proteins were enriched for the functions of adrenergic signaling in cardiomyocytes, fatty acid degradation and glutathione metabolism. In our research, we found differentially expressed proteins in longissimus dorsi tissue between Kazakh cattle and Xinjiang brown cattle. We predict that these proteins regulate muscle production and development through select enriched signaling pathways. This study provides novel insights into the roles of proteomes in cattle genetics and breeding.

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.

Structural variation and introgression from wild populations in East Asian cattle genomes confer adaptation to local environment

BACKGROUND

Structural variations (SVs) in individual genomes are major determinants of complex traits, including adaptability to environmental variables. The Mongolian and Hainan cattle breeds in East Asia are of taurine and indicine origins that have evolved to adapt to cold and hot environments, respectively. However, few studies have investigated SVs in East Asian cattle genomes and their roles in environmental adaptation, and little is known about adaptively introgressed SVs in East Asian cattle.

RESULTS

In this study, we examine the roles of SVs in the climate adaptation of these two cattle lineages by generating highly contiguous chromosome-scale genome assemblies. Comparison of the two assemblies along with 18 Mongolian and Hainan cattle genomes obtained by long-read sequencing data provides a catalog of 123,898 nonredundant SVs. Several SVs detected from long reads are in exons of genes associated with epidermal differentiation, skin barrier, and bovine tuberculosis resistance. Functional investigations show that a 108-bp exonic insertion in SPN may affect the uptake of Mycobacterium tuberculosis by macrophages, which might contribute to the low susceptibility of Hainan cattle to bovine tuberculosis. Genotyping of 373 whole genomes from 39 breeds identifies 2610 SVs that are differentiated along a "north-south" gradient in China and overlap with 862 related genes that are enriched in pathways related to environmental adaptation. We identify 1457 Chinese indicine-stratified SVs that possibly originate from banteng and are frequent in Chinese indicine cattle.

CONCLUSIONS

Our findings highlight the unique contribution of SVs in East Asian cattle to environmental adaptation and disease resistance.

Ancient Mongolian aurochs genomes reveal sustained introgression and management in East Asia

Societies in East Asia have utilized domesticated cattle for over 5000 years, but the genetic history of cattle in East Asia remains understudied. Genome-wide analyses of 23 ancient Mongolian cattle reveal that East Asian aurochs and ancient East Asian taurine cattle are closely related, but neither are closely related to any modern East Asian breeds. We observe binary variation in aurochs diet throughout the early Neolithic, and genomic evidence shows millennia of sustained male-dominated introgression. We identify a unique connection between ancient Mongolian aurochs and the European Hereford breed. These results point to the likelihood of human management of aurochs in Northeast Asia prior to and during the initial adoption of taurine cattle pastoralism.

Local ancestry and selection in admixed Sanjiang cattle

The majority of native cattle are taurine × indicine cattle of diverse phenotypes in the central region of China. Sanjiang cattle, a typical breed in the central region, play a central role in human livelihood and have good adaptability, including resistance to dampness, heat, roughage, and disease, and are thus regarded as an important genetic resource. However, the genetic history of the successful breed remains unknown. Here, we sequenced 10 Sanjiang cattle genomes and compared them to the 70 genomes of 5 representative populations worldwide. We characterized the genomic diversity and breed formation process of Sanjiang cattle and found that Sanjiang cattle have a mixed ancestry of indicine (55.6%) and taurine (33.2%) dating to approximately 30 generations ago, which has shaped the genome of Sanjiang cattle. Through ancestral fragment inference, selective sweep and transcriptomic analysis, we identified several genes linked to lipid metabolism, immune regulation, and stress reactions across the mosaic genome of Sanjiang cattle showing an excess of taurine or indicine ancestry. Taurine ancestry might contribute to meat quality, and indicine ancestry is more conducive to adaptation to hot climate conditions, making Sanjiang cattle a valuable genetic resource for the central region of China. Our results will help us understand the evolutionary history and ancestry components of Sanjiang cattle, which will provide a reference for resource conservation and selective breeding of Chinese native cattle.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

Ancient DNA confirms diverse origins of early post-Columbian cattle in the Americas

Before the arrival of Europeans, domestic cattle (Bos taurus) did not exist in the Americas, and most of our knowledge about how domestic bovines first arrived in the Western Hemisphere is based on historical documents. Sixteenth-century colonial accounts suggest that the first cattle were brought in small numbers from the southern Iberian Peninsula via the Canary archipelago to the Caribbean islands where they were bred locally and imported to other circum-Caribbean regions. Modern American heritage cattle genetics and limited ancient mtDNA data from archaeological colonial cattle suggest a more complex story of mixed ancestries from Europe and Africa. So far little information exists to understand the nature and timing of the arrival of these mixed-ancestry populations. In this study we combine ancient mitochondrial and nuclear DNA from a robust sample of some of the earliest archaeological specimens from Caribbean and Mesoamerican sites to clarify the origins and the dynamics of bovine introduction into the Americas. Our analyses support first arrival of cattle from diverse locales and potentially confirm the early arrival of African-sourced cattle in the Americas, followed by waves of later introductions from various sources over several centuries.

Whole genome sequencing reveals signals of adaptive admixture in Creole cattle

The Creole cattle from Guadeloupe (GUA) are well adapted to the tropical environment. Its admixed genome likely played an important role in such adaptation. Here, we sought to detect genomic signatures of selection in the GUA genome. For this purpose, we sequenced 23 GUA individuals and combined our data with sequenced genomes of 99 animals representative of European, African and indicine groups. We detect 17,228,983 single nucleotide polymorphisms (SNPs) in the GUA genome, providing the most detailed exploration, to date, of patterns of genetic variation in this breed. We confirm the higher level of African and indicine ancestries, compared to the European ancestry and we highlight the African origin of indicine ancestry in the GUA genome. We identify five strong candidate regions showing an excess of indicine ancestry and consistently supported across the different detection methods. These regions encompass genes with adaptive roles in relation to immunity, thermotolerance and physical activity. We confirmed a previously identified horn-related gene, RXFP2, as a gene under strong selective pressure in the GUA population likely owing to human-driven (socio-cultural) pressure. Findings from this study provide insight into the genetic mechanisms associated with resilience traits in livestock.