A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Pangenome and multi-tissue gene atlas provide new insights into the domestication and highland adaptation of yaks

BACKGROUND

The genetic diversity of yak, a key domestic animal on the Qinghai-Tibetan Plateau (QTP), is a vital resource for domestication and breeding efforts. This study presents the first yak pangenome obtained through the de novo assembly of 16 yak genomes.

RESULTS

We discovered 290 Mb of nonreference sequences and 504 new genes. Our pangenome-wide presence and absence variation (PAV) analysis revealed 5,120 PAV-related genes, highlighting a wide range of variety-specific genes and genes with varying frequencies across yak populations. Principal component analysis (PCA) based on binary gene PAV data classified yaks into three new groups: wild, domestic, and Jinchuan. Moreover, we proposed a 'two-haplotype genomic hybridization model' for understanding the hybridization patterns among breeds by integrating gene frequency, heterozygosity, and gene PAV data. A gene PAV-GWAS identified a novel gene (BosGru3G009179) that may be associated with the multirib trait in Jinchuan yaks. Furthermore, an integrated transcriptome and pangenome analysis highlighted the significant differences in the expression of core genes and the mutational burden of differentially expressed genes between yaks from high and low altitudes. Transcriptome analysis across multiple species revealed that yaks have the most unique differentially expressed mRNAs and lncRNAs (between high- and low-altitude regions), especially in the heart and lungs, when comparing high- and low-altitude adaptations.

CONCLUSIONS

The yak pangenome offers a comprehensive resource and new insights for functional genomic studies, supporting future biological research and breeding strategies.

Yak genome database: a multi-omics analysis platform

BACKGROUND

The yak (Bos grunniens) is a large ruminant species that lives in high-altitude regions and exhibits excellent adaptation to the plateau environments. To further understand the genetic characteristics and adaptive mechanisms of yak, we have developed a multi-omics database of yak including genome, transcriptome, proteome, and DNA methylation data.

DESCRIPTION

The Yak Genome Database ( http://yakgenomics.com/ ) integrates the research results of genome, transcriptome, proteome, and DNA methylation, and provides an integrated platform for researchers to share and exchange omics data. The database contains 26,518 genes, 62 transcriptomes, 144,309 proteome spectra, and 22,478 methylation sites of yak. The genome module provides access to yak genome sequences, gene annotations and variant information. The transcriptome module offers transcriptome data from various tissues of yak and cattle strains at different developmental stages. The proteome module presents protein profiles from diverse yak organs. Additionally, the DNA methylation module shows the DNA methylation information at each base of the whole genome. Functions of data downloading and browsing, functional gene exploration, and experimental practice were available for the database.

CONCLUSION

This comprehensive database provides a valuable resource for further investigations on development, molecular mechanisms underlying high-altitude adaptation, and molecular breeding of yak.

Genome Variation Map of Domestic Qinghai-Tibet Plateau Yaks by SLAF-Seq Reveals Genetic Footprint during Artificial Selection

The yak (Bos grunniens) was domesticated in the high-altitude QTP. Research about their genetic diversity and population structure is limited. In this study, we resequenced the genome of 494 domestic yaks using Specific-Locus Amplified Fragment Sequencing (SLAF-seq). The survey was conducted on six populations sampled from isolated locations in China in order to analyze their structure and genetic diversity. These six domestic populations were clearly grouped into two independent clusters, with Jinchuan, Changtai, and Jiulong showing a tight genetic relationship with the wild yak. Nerve development pathways were enriched with GO enrichment analysis of 334 domesticated genes. Major genomic regions associated with the differentiation of domestic yaks were detected. These findings provide preliminary information on the yak genome variability, useful to understand the genomic characteristics of different populations in QTP.

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 reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in Qinghai, China

The Qinghai Province of China is located in the northeast region of the Qinghai-Tibetan Plateau (QTP) and carries abundant yak genetic resources. Previous investigations of archaeological records, mitochondrial DNA, and Y chromosomal markers have suggested that Qinghai was the major center of yak domestication. In the present study, we examined the genomic diversity, differentiation, and selection signatures of 113 Qinghai yak, including 42 newly sequenced Qinghai yak and 71 publicly available individuals, from nine yak breeds/populations (wild, Datong, Huanhu, Xueduo, Yushu, Qilian, Geermu, Tongde, and Huzhu white) using high-depth whole-genome resequencing data. We observed that most of Qinghai yak breeds/populations have abundant genomic diversity based on four genomic parameters (nucleotide diversity, inbreeding coefficients, linkage disequilibrium decay, and runs of homozygosity). Population genetic structure analysis showed that Qinghai yak have two lineages with two ancestral origins and that nine yak breeds/populations are clustered into three distinct groups of wild yak, Geermu yak, and seven other domestic yak breeds/populations. F _ST values showed moderate genetic differentiation between wild yak, Geermu yak, and the other Qinghai yak breeds/populations. Positive selection signals were detected in candidate genes associated with disease resistance (CDK2AP2, PLEC, and CYB5B), heat stress (NFAT5, HSF1, and SLC25A48), pigmentation (MCAM, RNF26, and BOP1), vision (C1QTNF5, MFRP, and TAX1BP3), milk quality (OPLAH and GRINA), neurodevelopment (SUSD4, INSYN1, and PPP1CA), and meat quality (ZRANB1), using the integrated PI, composite likelihood ratio (CLR), and F _ST methods. These findings offer new insights into the genetic mechanisms underlying target traits in yak and provide important information for understanding the genomic characteristics of yak breeds/populations in Qinghai.

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.

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.

Molecular Regulation of Yak Preadipocyte Differentiation and Proliferation by LncFAM200B and ceRNA Regulatory Network Analysis

The positive regulatory role of lncFAM200B in differentiation and lipid deposition in yak intramuscular preadipocytes has been demonstrated in our previous study. However, the regulatory mechanisms remain unclear. In this study, we aimed to produce complete mRNA and microRNA (miRNA) profiles after adenovirus-mediated lncFAM200B overexpression in yak preadipocytes using high-throughput sequencing. We constructed a competing endogenous RNA (ceRNA) network with lncFAM200B as the core and identified the functions of the selected target miRNA during cell proliferation and differentiation. We obtained 118 differentially expressed genes (DEGs) after lncFAM200B overexpression, 76 of which were up-regulated, including Notch signaling members NOTCH3, DTX3L, and HES4, and 42 DEGs were down-regulated, including genes related to the cell cycle (CCNA2, BUB1, CDC20, TOP2A, and KIF20A). Additionally, many ubiquitin-mediated proteolysis pathway members were also significantly up-regulated (BUA7, PML, TRIM21, and TRIM25). MiRNA sequencing showed that 13 miRNAs were significantly up-regulated, and 12 miRNAs were down-regulated. Among them, 29 targets of 10 differentially expressed miRNAs (DEMs) were differentially expressed, including miR-152-FBXO33, miR-6529a-TRIM21, miR-148c-NOTCH3, and the miR-6529b-HES4 axis. We further verified that overexpression and inhibition of miR-6529a can inhibit and promote, respectively, the proliferation and differentiation of preadipocytes. Taken together, our study not only revealed the regulatory network of lncFAM200B during yak preadipocytes differentiation but also laid a foundation for elucidating the cause for lower intramuscular fat content in yaks at the molecular level.

Genome-Wide Association Study of Body Weight Trait in Yaks

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

Proteomics discovery of protein biomarkers linked to yak meat tenderness as determined by label-free mass spectrometry

Tenderness is one of the most important qualities in meat. A proteomic approach is a suitable way to ensure meat tenderness. Thirty-six tenderloin samples from yak were classified as exhibiting high (n = 12) or low (n = 12) tenderness and were evaluated using label-free proteomics for the identification of the proteins and pathways most influential in tenderness variability. Between the two groups, proteomic changes were mainly caused by 33 differentially expressed proteins as displayed in reference patterns in heat maps. The expression of ENO2, SUCLG2, ETFDH, PGM1, TNNT3, TNNT1, HSDL2, GPI, ALAD, and COL1A1 proteins was very different between yak meats with high and low tenderness, and therefore, they are candidate biomarkers of yak meat tenderness. Furthermore, bioinformatics analyses revealed that the identified proteins are related to pentose phosphate, glycolysis, the citrate cycle, fatty acid metabolism, and the calcium signaling pathway.

Genome-Wide DNA Methylation and Hydroxymethylation Changes Revealed Epigenetic Regulation of Neuromodulation and Myelination in Yak Hypothalamus

Both 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are important epigenetic modifications in neurodevelopment. However, there is little research examining the genome-wide patterns of 5mC and 5hmC in brain regions of animals under natural high-altitude conditions. We used oxidative reduced representation bisulfite sequencing (oxRRBS) to determine the 5mC and 5hmC sites in the brain, brainstem, cerebellum, and hypothalamus of yak and cattle. We reported the first map of genome-wide DNA methylation and hydroxymethylation in the brain, brainstem, cerebellum, and hypothalamus of yak (living at high altitudes) and cattle. Overall, we found striking differences in 5mC and 5hmC between the hypothalamus and other brain regions in both yak and cattle. Genome-wide profiling revealed that 5mC level decreased and 5hmC level increased in the hypothalamus than in other regions. Furthermore, we identified differentially methylated regions (DMRs) and differentially hydroxymethylated regions (DhMRs), most of which overlapped with each other. Interestingly, transcriptome results for these brain regions also showed distinctive gene levels in the hypothalamus. Finally, differentially expressed genes (DEGs) regulated by DMRs and DhMRs may play important roles in neuromodulation and myelination. Overall, our results suggested that mediation of 5mC and 5hmC on epigenetic regulation may broadly impact the development of hypothalamus and its biological functions.

Whole-Transcriptome Analysis of Yak and Cattle Heart Tissues Reveals Regulatory Pathways Associated With High-Altitude Adaptation

BACKGROUND

The yak (Bos grunniens) is an important livestock species that can survive the extremely cold, harsh, and oxygen-poor conditions of the Qinghai-Tibetan Plateau and provide meat, milk, and transportation for the Tibetans living there. However, the regulatory network that drive this hypoxic adaptation remain elusive.

RESULTS

The heart tissues from LeiRoqi (LWQY) yak and their related cattle (Bos Taurus) breeds, which are two native cattle breeds located in high altitude (HAC) and low altitude (LAC) regions, respectively, were collected for RNA sequencing. A total of 178 co-differentially expressed protein-coding transcripts (co-DETs) were discovered in each of the LAC-vs-LWQY and LAC-vs-HAC comparison groups, including NFATC2, NFATC1, ENPP2, ACSL4, BAD, and many other genes whose functions were reported to be associated with the immune-system, endocrine-system, and lipid metabolism. Two and 230 lncRNA transcripts were differentially expressed in the LAC-vs-LWQY and LAC-vs-HAC comparisons' respectively, but no lncRNA transcripts that were co-differentially expressed. Among the 58 miRNAs that were co-differentially expressed, 18 were up-regulated and 40 were down-regulated. In addition, 640 (501 up-regulated and 139 down-regulated) and 152 (152 up-regulated and one down-regulated) circRNAs showed differential expression in LAC-vs-LWQY and LAC-vs-HAC comparison groups, respectively, and 53 up-regulated co-differentially expressed circRNAs were shared. Multiple co-DETs, which are the targets of miRNAs/lncRNAs, are significantly enriched in high-altitude adaptation related processes, such as, T cell receptor signaling, VEGF signaling, and cAMP signaling. A competing endogenous RNA (ceRNA) network was constructed by integrating the competing relationships among co-differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs. Furthermore, the hypoxic adaptation related ceRNA network was constructed, and the six mRNAs (MAPKAPK3, PXN, NFATC2, ATP7A, DIAPH1, and F2R), the eight miRNAs (including miR-195), and 15 circRNAs (including novel-circ-017096 and novel-circ-018073) are proposed as novel and promising candidates for regulation of hypoxic adaptation in the heart.

CONCLUSION

In conclusion, the data recorded in the present study provides new insights into the molecular network of high-altitude adaptation along with more detailed information of protein-coding transcripts and non-coding transcripts involved in this physiological process, the detailed mechanisms behind how these transcripts "crosstalk" with each other during the plateau adaptation are worthy of future research efforts.

Methylome and transcriptome profiles in three yak tissues revealed that DNA methylation and the transcription factor ZGPAT co-regulate milk production

Background

Domestic yaks play an indispensable role in sustaining the livelihood of Tibetans and other ethnic groups on the Qinghai-Tibetan Plateau (QTP), by providing milk and meat. They have evolved numerous physiological adaptations to high-altitude environment, including strong blood oxygen transportation capabilities and high metabolism. The roles of DNA methylation and gene expression in milk production and high-altitudes adaptation need further exploration.

Results

We performed genome-wide DNA methylome and transcriptome analyses of breast, lung, and biceps brachii muscle tissues from yaks of different ages. We identified 432,350 differentially methylated regions (DMRs) across the age groups within each tissue. The post-mature breast tissue had considerably more differentially methylated regions (155,957) than that from the three younger age groups. Hypomethylated genes with high expression levels might regulate milk production by influencing protein processing in the endoplasmic reticulum. According to weighted gene correlation network analysis, the “hub” gene ZGPAT was highly expressed in the post-mature breast tissue, indicating that it potentially regulates the transcription of 280 genes that influence protein synthesis, processing, and secretion. The tissue network analysis indicated that high expression of HIF1A regulates energy metabolism in the lung.

Conclusions

This study provides a basis for understanding the epigenetic mechanisms underlying milk production in yaks, and the results offer insight to breeding programs aimed at improving milk production.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s12864-020-07151-3.

A chromosome-scale reference genome and genome-wide genetic variations elucidate adaptation in yak

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu Kush ranges of the Himalayas. The yak genome was sequenced in 2012, but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyse it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied towards high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein‐coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here, we identified breed‐specific sequences and genes by whole‐genome alignment, which may facilitate yak breed identification.