The sequence and de novo assembly of the wild yak genome

Joint exploration of network pharmacology and metabolomics on the effects of traditional Chinese medicine compounds in weaned yaks

Introduction

Chinese herbal medicines are relatively inexpensive and have fewer side effects, making them an effective option for improving health and treating diseases. As a result, they have gained more attention in recent years. The weaning period is a critical stage in the life of yaks, often inducing stress in calves. Weaning stress, along with dietary changes, can lead to a decline in physical fitness and immune function, making yaks more susceptible to diarrhea and resulting in high mortality rates during this period. Therefore, our study aimed to address this issue by incorporating traditional Chinese medicine (TCM) formulas into the diet of yaks during the weaning period.

Methods

Following a dialectical analysis, three TCM formulas, mainly composed of Paeonia lactiflora, Coptis chinensis, and Dandelion, were identified for their anti-inflammatory, antioxidant, and immune enhancing potentials. We explored the possible molecular mechanisms of these TCM formulas using network pharmacology analysis and investigated their effects on the physiology of yaks through metabolomics.

Results

Network pharmacology analysis revealed several key target proteins in the protein-protein interaction (PPI) network between three formulas and immune-related genes, including PIK3R1, PIK3CA, JAK2, PTK2, and PYPN11. The key target proteins in the PPI network associated with metabolism-related genes included ENPP1, CYP1A1, PTGS1, members of the CYP1 family, and EPHX2. GO analysis of co-targets revealed highly enriched pathways such as protein phosphorylation, plasma membrane, and one-carbon metabolic processes. Metabolomics revealed significant changes in the abundance of metabolites including dimethyl sulfoxide, tyrphostin A25, and thromboxane A2 in the intestines of weaned yaks supplemented with these Chinese herbal compounds. Significant changes were also observed in pathways such as vitamin A metabolism, chloroalkane, and chloroalkene degradation.

Discussion

Based on these findings, it can be inferred that TCM formulas improve the physical fitness of weaned yaks by enhancing antioxidant capacity, boosting immunity, and reducing intestinal inflammation. This study preliminarily elucidates the pharmacological mechanisms by which TCM formulas prevent diarrhea and improve physical fitness in weaned yaks through metabolomics and network pharmacology, paving the way for further evaluation of the effectiveness of these three formulas.

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Metagenomic analysis reveals high diversity of gut viromes in yaks (Bos grunniens) from the Qinghai-Tibet Plateau

The Qinghai-Tibet Plateau (QTP), renowned for its exceptional biological diversity, is home to numerous endemic species. However, research on the virology of vulnerable vertebrates like yaks remains limited. In this study, our objective was to use metagenomics to provide a comprehensive understanding of the diversity and evolution of the gut virome in yak populations across different regions of the QTP. Our findings revealed a remarkably diverse array of viruses in the gut of yaks, including those associated with vertebrates and bacteriophages. Notably, some vertebrate-associated viruses, such as astrovirus and picornavirus, showed significant sequence identity across diverse yak populations. Additionally, we observed differences in the functional profiles of genes carried by the yak gut virome across different regions. Moreover, the virus-bacterium symbiotic network that we discovered holds potential significance in maintaining the health of yaks. Overall, this research expands our understanding of the viral communities in the gut of yaks and highlights the importance of further investigating the interactions between viruses and their hosts. These data will be beneficial for revealing the crucial role that viruses play in the yak gut ecology in future studies.

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.

Evolutionary origin of genomic structural variations in domestic yaks

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Genome-wide identification of bovine ADAMTS gene family and analysis of its expression profile in the inflammatory process of mammary epithelial cells

ADAM metallopeptidase with thrombospondin type 1 motif (ADAMTS) are secreted, multi-domain matrix-related zinc endopeptidases that play a role in organogenesis, assembly and degradation of extracellular matrix (ECM), cancer and inflammation. Genome-wide identification and analysis of the bovine ADAMTS gene family has not yet been carried out. In this study, 19 ADAMTS family genes were identified in Bos taurus by genome-wide bioinformatics analysis, and they were unevenly distributed on 12 chromosomes. Phylogenetic analysis shows that the Bos taurus ADAMTS are divided into eight subfamilies, with highly consistent gene structures and motifs within the same subfamily. Collinearity analysis showed that the Bos taurus ADAMTS gene family is homologous to other bovine subfamily species, and many ADAMTS genes may be derived from tandem replication and segmental replication. In addition, based on the analysis of RNA-seq data, we found the expression pattern of ADAMTS gene in different tissues. Meanwhile, we also analyzed the expression profile of ADAMTS gene in the inflammatory response of bovine mammary epithelial cells (BMECs) stimulated by LPS by qRT-PCR. The results can provide ideas for understanding the evolutionary relationship and expression pattern of ADAMTS gene in Bovidae, and clarify the theoretical basis of the function of ADAMTS in inflammation.

Whole-Genome Resequencing Highlights the Unique Characteristics of Kecai Yaks

Kecai yaks are regarded as an important genetic resource in China owing to their high fecundity and flavorful meat. However, the genetic characteristics of Kecai yaks have not been effectively characterized to date, and the relationship between Kecai yaks and other yak breeds remains to be fully characterized. In this paper, the resequencing of the Kecai yak genome is performed leading to the identification of 11,491,383 high-quality single nucleotide polymorphisms (SNPs). Through principal component, phylogenetic, and population genetic structure analyses based on these SNPs, Kecai yaks were confirmed to represent an independent population of yaks within China. In this study, marker and functional enrichment analysis of genes related to positive selection in Kecai yak was carried out, and the results show that such selection in Kecai yaks is associated with the adaptation to alpine environments and the deposition of muscle fat. Overall, these results offer a theoretical foundation for the future utilization of Kecai yak genetic resources.

Long read genome assemblies complemented by single cell RNA-sequencing reveal genetic and cellular mechanisms underlying the adaptive evolution of yak

Wild yak (Bos mutus) and domestic yak (Bos grunniens) are adapted to high altitude environment and have ecological, economic, and cultural significances on the Qinghai-Tibetan Plateau (QTP). Currently, the genetic and cellular bases underlying adaptations of yak to extreme conditions remains elusive. In the present study, we assembled two chromosome-level genomes, one each for wild yak and domestic yak, and screened structural variants (SVs) through the long-read data of yak and taurine cattle. The results revealed that 6733 genes contained high-FST SVs. 127 genes carrying special type of SVs were differentially expressed in lungs of the taurine cattle and yak. We then constructed the first single-cell gene expression atlas of yak and taurine cattle lung tissues and identified a yak-specific endothelial cell subtype. By integrating SVs and single-cell transcriptome data, we revealed that the endothelial cells expressed the highest proportion of marker genes carrying high-FST SVs in taurine cattle lungs. Furthermore, we identified pathways which were related to the medial thickness and formation of elastic fibers in yak lungs. These findings provide new insights into the high-altitude adaptation of yak and have important implications for understanding the physiological and pathological responses of large mammals and humans to hypoxia.

The genetic bases of yak adaptations to extreme conditions remains elusive. This study compares yak and cattle at a genomic and transcriptomic level, revealing a new type of endothelial cell and candidate genes related with elastic fiber formation in yak lungs that might contribute to high altitude adaptation.

The Landscape of Accessible Chromatin during Yak Adipocyte Differentiation

Although significant advancement has been made in the study of adipogenesis, knowledge about how chromatin accessibility regulates yak adipogenesis is lacking. We here described genome-wide dynamic chromatin accessibility in preadipocytes and adipocytes by using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), and thus revealed the unique characteristics of open chromatin during yak adipocyte differentiation. The chromatin accessibility of preadipocytes and adipocytes exhibited a similar genomic distribution, displaying a preferential location within the intergenic region, intron, and promoter. The pathway enrichment analysis identified that genes with differential chromatin accessibility were involved in adipogenic metabolism regulation pathways, such as the peroxisome proliferator activated receptor-γ (PPAR) signaling pathway, wingless-type MMTV integration site (Wnt) signaling pathway, and extracellular matrix-receptor (ECM–receptor) interaction. Integration of ATAC-seq and mRNA-seq revealed that genes with a high expression were associated with high levels of chromatin accessibility, especially within 1 kb upstream and downstream of the transcription start site. In addition, we identified a series of transcription factors (TFs) related to adipogenesis and created the TF regulatory network, providing the possible interactions between TFs during yak adipogenesis. This study is crucial for advancing the understanding of transcriptional regulatory mechanisms of adipogenesis and provides valuable information for understanding the adaptation of plateau species to high-altitude environments by maintaining whole body homeostasis through fat metabolism.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

The Gut Microbiota Determines the High-Altitude Adaptability of Tibetan Wild Asses (Equus kiang) in Qinghai-Tibet Plateau

It was acknowledged long ago that microorganisms have played critical roles in animal evolution. Tibetan wild asses (TWA, Equus kiang) are the only wild perissodactyls on the Qinghai-Tibet Plateau (QTP) and the first national protected animals; however, knowledge about the relationships between their gut microbiota and the host's adaptability remains poorly understood. Herein, 16S rRNA and meta-genomic sequencing approaches were employed to investigate the gut microbiota–host associations in TWA and were compared against those of the co-resident livestock of yak (Bos grunnies) and Tibetan sheep (Ovis aries). Results revealed that the gut microbiota of yak and Tibetan sheep underwent convergent evolution. By contrast, the intestinal microflora of TWA diverged in a direction enabling the host to subsist on sparse and low-quality forage. Meanwhile, high microbial diversity (Shannon and Chao1 indices), cellulolytic activity, and abundant indicator species such as Spirochaetes, Bacteroidetes, Prevotella_1, and Treponema_2 supported forage digestion and short-chain fatty acid production in the gut of TWA. Meanwhile, the enterotype identification analysis showed that TWA shifted their enterotype in response to low-quality forage for a better utilization of forage nitrogen and short-chain fatty acid production. Metagenomic analysis revealed that plant biomass degrading microbial consortia, genes, and enzymes like the cellulolytic strains (Prevotella ruminicola, Ruminococcus flavefaciens, Ruminococcus albus, Butyrivibrio fibrisolvens, and Ruminobacter amylophilus), as well as carbohydrate metabolism genes (GH43, GH3, GH31, GH5, and GH10) and enzymes (β-glucosidase, xylanase, and β-xylosidase, etc.) had a significantly higher enrichment in TWA. Our results indicate that gut microbiota can improve the adaptability of TWA through plant biomass degradation and energy maintenance by the functions of gut microbiota in the face of nutritional deficiencies and also provide a strong rationale for understanding the roles of gut microbiota in the adaptation of QTP wildlife when facing harsh feeding environments.

Mental health and social isolation under repeated mild lockdowns in Japan

The influence of repeated lockdowns on mental health and social isolation is unknown. We conducted a longitudinal study of the influence of repeated mild lockdowns during two emergency declarations in Japan, in May 2020 and February 2021. The analyses included 7893 people who participated in all online surveys. During repeated mild lockdowns, mental and physical symptoms decreased overall, while loneliness increased and social networks decreased. Subgroup analyses revealed that depression and suicidal ideation did not decrease only in the younger age group (aged 18-29 years) and that younger and middle-aged people (aged 18-49 years), women, people with a history of treatment for mental illness, and people who were socially disadvantaged in terms of income had higher levels of mental and physical symptoms at all survey times. Additionally, comprehensive extraction of the interaction structure between depression, demographic attributes, and psychosocial variables indicated that loneliness and social networks were most closely associated with depression. These results indicate that repeated lockdowns have cumulative negative effects on social isolation and loneliness and that susceptible populations, such as young people and those with high levels of loneliness, require special consideration during repeated lockdown situations.

Intestinal Microbial Diversity of Free-Range and Captive Yak in Qinghai Province

BACKGROUND

The gut microbiome is a large and complex organic assemblage with subtle and close relationships with the host. This symbiotic mechanism is important for the health and adaptability of the host to the environment. Compared with other ruminants, there are few studies on yak intestinal microbes. The study of the gut microbiota of the yak will help us better understand the correlation between the microbiota and the environmental adaptability of the host. In this study, we adapted 16S rDNA sequencing technology to investigate the diversity and composition of the intestinal microbial community in free-range yaks and captive yaks living on the Qinghai-Tibet Plateau (QTP).

RESULTS

Sequencing results showed that the intestinal microbial community diversity was significantly different between free-range yaks and captive yaks. Firmicutes and Bacteroidetes were the dominant bacteria in both free-range and captive yaks. However, there were differences between the microbes of the two analyzed feeding styles in different classification levels. Compared with the captive type, free-range yaks had a higher abundance of Ruminococcaceae, Eubacteriaceae, Desulfovibrionaceae, Elusimicrobium, and Oscillibacter, while the abundance of Succinivibrionaceae, Clostridiales, Lachnospiraceae, Prevotellaceae, Roseburia, and Barnesiella was relatively low. The feeding method may be the key factor for the formation of intestinal flora differences in yaks, while altitude did not significantly affect Qinghai yak.

CONCLUSIONS

In this study, we used 16S rDNA sequencing technology to investigate the composition of intestinal flora in free-range and captive yaks living on the QTP. The exploration of dietary factors can provide a theoretical basis for scientifically and rationally breeding yaks and provides a new direction for the development of prebiotics and microecological agents.

Structural Variants Selected during Yak Domestication Inferred from Long-Read Whole-Genome Sequencing

Structural variants (SVs) represent an important genetic resource for both natural and artificial selection. Here we present a chromosome-scale reference genome for domestic yak (Bos grunniens) that has longer contigs and scaffolds (N50 44.72 and 114.39 Mb, respectively) than reported for any other ruminant genome. We further obtained long-read resequencing data for 6 wild and 23 domestic yaks and constructed a genetic SV map of 372,220 SVs that covers the geographic range of the yaks. The majority of the SVs contains repetitive sequences and several are in or near genes. By comparing SVs in domestic and wild yaks, we identified genes that are predominantly related to the nervous system, behavior, immunity, and reproduction and may have been targeted by artificial selection during yak domestication. These findings provide new insights in the domestication of animals living at high altitude and highlight the importance of SVs in animal domestication.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus)

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau. So far, most studies on the molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on the small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) and four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. The present study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus)

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau (QTP). So far, most studies on molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity (Hd) at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) with four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. This study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

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.