Integration of ATAC-Seq and RNA-Seq Analysis to Identify Key Genes in the Longissimus Dorsi Muscle Development of the Tianzhu White Yak

During the postnatal stages, skeletal muscle development undergoes a series of meticulously regulated alterations in gene expression. However, limited studies have employed chromatin accessibility to unravel the underlying molecular mechanisms governing muscle development in yak species. Therefore, we conducted an analysis of both gene expression levels and chromatin accessibility to comprehensively characterize the dynamic genome-wide chromatin accessibility during muscle growth and development in the Tianzhu white yak, thereby elucidating the features of accessible chromatin regions throughout this process. Initially, we compared the differences in chromatin accessibility between two groups and observed that calves exhibited higher levels of chromatin accessibility compared to adult cattle, particularly within ±2 kb of the transcription start site (TSS). In order to investigate the correlation between alterations in chromatin accessible regions and variations in gene expression levels, we employed a combination of ATAC-seq and RNA-seq techniques, leading to the identification of 18 central transcriptional factors (TFs) and 110 key genes with significant effects. Through further analysis, we successfully identified several TFs, including Sp1, YY1, MyoG, MEF2A and MEF2C, as well as a number of candidate genes (ANKRD2, ANKRD1, BTG2 and LMOD3) which may be closely associated with muscle growth and development. Moreover, we constructed an interactive network program encompassing hub TFs and key genes related to muscle growth and development. This innovative approach provided valuable insights into the molecular mechanism underlying skeletal muscle development in the postnatal stages of Tianzhu white yaks while also establishing a solid theoretical foundation for future research on yak muscle development.

Transcriptomic analysis reveals differentially expressed genes and key immune pathways in the spleen of the yak (Bos grunniens) at different growth stage

Yak is one of the rare and unique cattle species on the Qinghai-Tibetan Plateau, which has strong adaptability to the extreme environment of the plateau. The spleens are important functional organs that enable animals to adapt to their external environment and are vital in the growth and development process. To further investigate changes in immune function during yak development, we compared the transcriptome profiles of spleen tissues among juvenile (1-day old), youth (15-months old), and prime (5-years old) yaks. Immunology of spleen development was evaluated based on histological analyses and global gene expression was examined by using RNA-sequencing (RNA-seq) technology. In this work, we found 6378 genes with significant differences between the spleen of juvenile yak and youth yak, with the largest difference between groups. There were 3144 genes with significant differences between the spleen of young yak and prime yak, with the smallest differences between groups. Further, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted for the functional annotation of these genes. GO and KEGG analysis showed that some of them were related to growth, disease, immune, and metabolism. However, the genetic mechanism underlying the adaptability of yak spleens at different ages to harsh plateau environments remains unknown. These findings are important for studying the mechanisms of spleen development in yaks of different age groups.

Differentiation of Bos grunniens and Bos taurus based on STR locus polymorphism

Differentiation of closely related biological species using molecular genetic analysis is important for breeding farm animals, creating hybrid lines, maintaining the genetic purity of breeds, lines and layering. Bos grunniens and Bos taurus differentiation based on STR locus polymorphism will help maintain the genetic isolation of these species and identify hybrid individuals. The aim of this study is to assess the differentiating potential of 15 microsatellite loci to distinguish between domestic yak (B. grunniens) bred in the Kalmak-Ashuu highland region (Kochkor district, Naryn region, Kyrgyz Republic) and cattle (B. taurus) of three breeds (Aberdeen-Angus, Holstein and Alatau) using molecular genetic analysis. The samples were genotyped at 15 microsatellite loci (ETH3, INRA023, TGLA227, TGLA126, TGLA122, SPS115, ETH225, TGLA53, BM2113, BM1824, ETH10, BM1818, CSSM66, ILSTS006 and CSRM60). Twelve of the loci were from the standard markers panel recommended by ISAG. Statistical analysis was performed using GenAlEx v.6.503, Structure v.2.3.4, PAST v.4.03, and POPHELPER v1.0.10. The analysis of the samples' subpopulation structure using the Structure v.2.3.4 and 15 STR locus genotyping showed that the accuracy of assigning a sample to B. taurus was 99.6 ± 0.4 %, whereas the accuracy of assigning a sample to B. grunniens was 99.2 ± 2.6 %. Of the 15 STRs, the greatest potential to differentiate B. grunniens and B. taurus was found in those with the maximal calculated FST values, including BM1818 (0.056), BM1824 (0.041), BM2113 (0.030), CSSM66 (0.034) and ILSTS006 (0.063). The classification accuracy of B. grunniens using only these five microsatellite loci was 98.8 ± 3.4 %, similar for B. taurus, 99.1 ± 1.2 %. The proposed approach, based on the molecular genetic analysis of 5 STR loci, can be used as an express test in Kyrgyzstan breeding and reproduction programs for B. grunniens.

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

The Novel Structural Variation in the GHR Gene Is Associated with Growth Traits in Yaks (Bos grunniens)

The growth hormone receptor (GHR) is a member of the cytokine/hematopoietic factor receptor superfamily, which plays an important role in the growth and development, immunity, and metabolism of animals. This study identified a 246 bp deletion variant in the intronic region of the GHR gene, and three genotypes, including type II, type ID, and type DD, were observed. Genotype analysis of structural variation (SV) was performed on 585 individuals from 14 yak breeds, and it was found that 246 bp deletion was present in each breed. The II genotype was dominant in all yak breeds except for SB yak. The association analysis of gene polymorphisms and growth traits in the ASD yak population showed that the 246 bp SV was significantly associated with body length at 6 months (p < 0.05). GHR messenger RNA (mRNA) was expressed in all the tested tissues, with significantly higher levels in the liver, muscle, and fat than in other organs. The results of transcription activity showed that the luciferase activity of the pGL4.10-DD vector was significantly higher than that of the pGL4.10-II vector (p < 0.05). Additionally, the transcription-factor binding prediction results showed that the SV in the runt-related transcription factor 1 (Runx1) transcription-factor binding site may affect the transcriptional activity of the GHR gene, regulating yak growth and development. This study showed that the novel SV of the GHR gene could be used as a candidate molecular marker for the selection of the early growth trait in ASD yak.

Expression and Localization of Fas-Associated Factor 1 in Testicular Tissues of Different Ages and Ovaries at Different Reproductive Cycle Phases of Bos grunniens

Fas-associated factor 1 (FAF1), a member of the Fas family, is involved in biological processes such as apoptosis, inflammation, cell proliferation and proteostasis. This study aimed to explore the biological role of FAF1 in testicular tissue at different ages (juveniles (1 and 2 years old), adults (3, 4, 6, and 7 years old) and old-aged animals (11 years old)) and ovaries during different reproductive cycle phases (follicular, luteal, and pregnancy phases). FAF1 mRNA, relative protein expression and protein expression localization were determined in testes and ovaries using real-time quantification, WB and immunohistochemistry (IHC), respectively. Real-time quantification of testis tissues showed that the relative expression of FAF1 mRNA in testis tissues at 3, 4 and 7 years of age was significantly higher than of those in other ages, and in ovarian tissues was significantly higher in luteal phase ovaries than those in follicular and pregnancy phase ovaries; follicular phase ovaries were the lowest. WB of testis tissues showed that the relative protein expression of FAF1 protein was significantly higher at 11 and 7 years of age; in ovarian tissue, the relative protein expression of FAF1 protein was significantly higher in follicular phase ovaries than in luteal and pregnancy phase ovaries, and lowest in luteal phase ovaries. The relative protein expression of FAF1 at 3, 4 and 7 years of age was the lowest. IHC showed that FAF1 was mainly expressed in spermatozoa, spermatocytes, spermatogonia and supporting cells; in ovarian tissue, FAF1 was expressed in ovarian germ epithelial cells, granulosa cells, cumulus cells and luteal cells. The IHC results showed that FAF1 mRNA and protein were significantly differentially expressed in testes of different ages and ovarian tissues of different reproductive cycle phases, revealing the significance of FAF1 in the regulation of male and female B. grunniens reproductive physiology. Furthermore, our results provide a basis for the further exploration of FAF1 in the reproductive physiology of B. grunniens.

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.

Comprehensive Analysis of mRNA, lncRNA, circRNA, and miRNA Expression Profiles and Their ceRNA Networks in the Longissimus Dorsi Muscle of Cattle-Yak and Yak

Cattle-yak, as the hybrid offspring of cattle (Bos taurus) and yak (Bos grunniens), demonstrates obvious heterosis in production performance. Male hybrid sterility has been focused on for a long time; however, the mRNAs and non-coding RNAs related to muscle development as well as their regulatory networks remain unclear. The phenotypic data showed that the production performance (i.e., body weight, withers height, body length, and chest girth) of cattle-yak was significantly better than that of the yak, and the economic benefits of the cattle-yak were higher under the same feeding conditions. Then, we detected the expression profiles of the longissimus dorsi muscle of cattle-yak and yak to systematically reveal the molecular basis using the high-throughput sequencing technology. Here, 7,126 mRNAs, 791 lncRNAs, and 1,057 circRNAs were identified to be differentially expressed between cattle-yaks and yaks in the longissimus dorsi muscle. These mRNAs, lncRNA targeted genes, and circRNA host genes were significantly enriched in myoblast differentiation and some signaling pathways related to muscle development (such as HIF-1 signaling pathway and PI3K-Akt signaling pathway). We constructed a competing endogenous RNA (ceRNA) network and found that some non-coding RNAs differentially expressed may be involved in the regulation of muscle traits. Taken together, this study may be used as a reference tool to provide the molecular basis for studying muscle development.

Molecular cloning, sequence, and expression patterns of DNA damage induced transcript 3 (DDIT3) gene in female yaks (Bos grunniens)

Endoplasmic reticulum stress (ERS) plays an important role in regulating the reproductive process of female mammals, mainly involved in follicular atresia and corpus luteum regression. DNA damage induced transcript 3 (DDIT3) is a marker gene of ERS. The objectives of the present study were to clone and analyze the sequence and tissue expression characteristics of DDIT3 gene in female yaks. By reverse transcriptase-polymerase chain reaction (RT-PCR) strategy, we obtained full-length 507-bp DDIT3-cDNA, encoding for 168-aa protein. Yak DDIT3 exhibited highest and least identity with that of bison and horse, respectively. Real-time PCR analyses revealed that the expression level of DDIT3 gene in ovary was higher than that in heart, liver, kidney, spleen, lung, uterus and oviduct (p < 0.05). DDIT3 expression level in ovary and uterus during pregnancy was higher than that in follicular phase, luteal phase and fetus stage. DDIT3 was highly expressed in metaphase II oocytes and granulosa cells than that in germinal vesicle and metaphase I oocytes (p < 0.05), respectively. This is the first molecular characterization and expression patterns of DDIT3 gene in female yaks. These results indicated that the DDIT3 gene possibly plays an important role in regulating ovary function and pregnancy maintenance in yaks.

Using weighted gene co-expression network analysis (WGCNA) to identify the hub genes related to hypoxic adaptation in yak (Bos grunniens)

BACKGROUND

As a mammal living at the highest altitude in the world, the yak has strong adaptability to the harsh natural environment (such as low temperature, scarce food, especially low oxygen) of Qinghai-Tibet Plateau (QTP) after a long process of natural selection.

OBJECTIVE

Here, we used Weighted Correlation Network Analysis (WGCNA), a systematic biology method, to identify hypoxic adaptation-related modules and hub genes. The research of the adaptability of yak against hypoxia is of great significance to identify the genetic characteristics and yak breeding.

METHODS

Based on the transcriptome sequencing data (PRJNA362606), the R package DESeq2 and WGCNA were conducted to analyze differentially expressed genes (DEGs) and construct the gene co-expression network. The module hub genes were identified and characterized by the correlation of gene and trait, module membership (kME). In addition, GO and KEGG enrichment analyses were used to explore the functions of hub genes.

RESULTS

Our results revealed that 1098, 1429, and 1645 DEGs were identified in muscle, spleen, and lung, respectively. Besides, a total of 13 gene co-expression modules were detected, of which two hypoxic adaptation-related modules (saddlebrown and turquoise) were found. We identified 39 and 150 hub genes in these two modules. Functional enrichment analyses showed that 12 GO terms and 18 KEGG pathways were enriched in the saddlebrown module while 85 GO terms and 22 KEGG pathways were enriched in the turquoise module. The significant pathways related to hypoxia adaptation include FoxO signaling pathway, Thermogenesis pathway, and Retrograde endocannabinoid signaling pathway, etc. CONCLUSIONS: In this study, we obtained two hypoxia-related specific modules and identified hub genes based on the connectivity by constructing a weighted gene co-expression network. Function enrichment analysis of two modules revealed mitochondrion is the most important organelle for hypoxia adaptation. Moreover, the insulin-related pathways and thermogenic-related pathways played a major role. The results of this study provide theoretical guidance for further understanding the molecular mechanism of yak adaptation to hypoxia.

Characterization of whole mitogenome sequence of the Tongde yak (Bos grunniens)

Tongde County is located in the southeast of Qinghai Province, China, harboring rich yak genetic resources. In the present study, the complete mitochondrial genome (mitogenome) of the Tongde yak (Bos grunniens) was firstly sequenced using Illumina sequencing technique and the corresponding sequence characterization was identified. Our results showed that the mitogenome of Tongde yak is a circular molecule with 16,323 bp length consisting of 37 genes (13 protein-coding genes, 2 rRNA genes, 22 tRNA genes) and a non-coding control region (D-loop), which is consistent with most bovine species. The overall nucleotide composition was found as: A (33.72%), T (27.27%), C (25.80%), and G (13.21%), respectively, yielding a higher AT content (60.99%). The complete mitogenome sequence of Tongde yak would provide useful information for further studies on its genetic resource conservation and molecular breeding programmes in the future.

The whole-transcriptome landscape of muscle and adipose tissues reveals the ceRNA regulation network related to intramuscular fat deposition in yak

Background

The Intramuscular fat (IMF) content in meat products, which is positively correlated with meat quality, is an important trait considered by consumers. The regulation of IMF deposition is species specific. However, the IMF-deposition-related mRNA and non-coding RNA and their regulatory network in yak (Bos grunniens) remain unknown. High-throughput sequencing technology provides a powerful approach for analyzing the association between transcriptome-related differences and specific traits in animals. Thus, the whole transcriptomes of yak muscle and adipose tissues were screened and analyzed to elucidate the IMF deposition-related genes. The muscle tissues were used for IMF content measurements.

Results

Significant differences were observed between the 0.5- and 2.5-year-old yaks. Several mRNAs, miRNAs, lncRNAs and circRNAs were generally expressed in both muscle and adipose tissues. Between the 0.5- and 2.5-year-old yaks, 149 mRNAs, 62 miRNAs, 4 lncRNAs, and 223 circRNAs were differentially expressed in muscle tissue, and 72 mRNAs, 15 miRNAs, 9 lncRNAs, and 211 circRNAs were differentially expressed in adipose tissue. KEGG annotation revelved that these differentially expressed genes were related to pathways that maintain normal biological functions of muscle and adipose tissues. Moreover, 16 mRNAs, 5 miRNAs, 3 lncRNAs, and 5 circRNAs were co-differentially expressed in both types of tissue. We suspected that these co-differentially expressed genes were involved in IMF-deposition in the yak. Additionally, LPL, ACADL, SCD, and FASN, which were previously shown to be associated with the IMF content, were identified in the competing endogenous RNA (ceRNA) regulatory network that was constructed on the basis of the IMF deposition-related genes. Three ceRNA subnetworks also revealed that TCONS-00016416 and its target SIRT1 “talk” to each other through the same miR-381-y and miR-208 response elements, whereas TCONS-00061798 and its target PRKCA, and TCONS-00084092 and its target LPL “talk” to each other through miR-122-x and miR-499-y response elements, respectively.

Conclusion

Taken together, our results reveal the potential mRNA and noncoding RNAs involved in IMF deposition in the yak, providing a useful resource for further research on IMF deposition in this animal species.

A Global Analysis of Y-STR INRA189 Polymorphism in Chinese Domestic Yak Breeds/Populations

Simple Summary

Y chromosome-specific markers have been widely used in studying the origin, migration, diversity, and population relationship in several mammalian species. So far, the investigations on yak paternal genetics were mostly based on Y-SNPs information, but very few on Y-STR markers. In this study, we comprehensively investigated the polymorphism of Y-STR INRA189 locus, identified the geographical distribution pattern of its alleles in 15 Chinese domestic yak breeds/populations, and provided basic data for yak paternal genetic analysis. Our results showed that INRA189 is an intermediate polymorphic Y-STR marker with six alleles present among 15 Chinese domestic yak breeds/populations and that the variation among six alleles was due to different numbers of repetitions of a TG dinucleotide motif.

Abstract

The objective of this study was to probe into the polymorphism of Y-STR INRA189 and identify the geographical distribution pattern of its alleles in the Chinese domestic yak gene pool. We examined the variation at INRA189 locus in 682 male yaks representing 15 breeds/populations in China. The results showed that six alleles, including five reported previously (149, 155, 157, 159, and 161 bp) and a new one identified in this study (139 bp), were detected at INRA189 locus based on genotyping analysis. The frequencies of six alleles varied among the 15 yak breeds/populations with a clear phylogeographical pattern, which revealed the paternal genetic difference among Chinese yak breeds/populations. The average polymorphism information content (PIC) among the 15 yak breeds/populations was 0.32, indicating INRA189 to be an intermediate polymorphic Y-STR marker (0.25 < PIC < 0.5) in yak. Sequence alignment revealed that the variations among six alleles at INRA189 were defined by a TG dinucleotide repeat motif, which repeated for 12, 17, 20, 21, 22, and 23 times, corresponding to the alleles of 139, 149, 155, 157, 159, and 161 bp, respectively. Therefore, we believe that the polymorphic yak Y-STR INRA189 can be used to characterize male-mediated genetic events, including paternal genetic origin, diversity, and evolution.

Bacterial communities in the solid, liquid, dorsal, and ventral epithelium fractions of yak (Bos grunniens) rumen

Yak (Bos grunniens) is an important and dominant livestock species in the challenging environment of the Qinghai–Tibetan Plateau. Rumen microbiota of the solid, liquid, and epithelium fractions play key roles in nutrient metabolism and contribute to host adaptation in ruminants. However, there is a little knowledge of the microbiota in these rumen fractions of yak. Therefore, we collected samples of solid, liquid, dorsal, and ventral epithelium fractions from five female yaks, then amplified bacterial 16S rRNA gene V4 regions and sequenced them using an Illumina MiSeq platform. Principal coordinates analysis detected significant differences in bacterial communities between the liquid, solid, and epithelium fractions, and between dorsal and ventral epithelium fractions. Rikenellaceae RC9, the families Lachnospiraceae and Ruminococcaceae, and Fibrobacter spp. were the abundant and enriched bacteria in solid fraction, while the genera Prevotella and Prevotellaceae UCG 003 were higher in the liquid fraction. Campylobacter spp., Comamonas spp., Desulfovibrio spp., and Solobacterium spp. were significantly higher in dorsal epithelium, while Howardella spp., Prevotellaceae UCG 001, Ruminococcaceae UCG 005, and Treponema 2 were enriched in the ventral epithelium. Comparison of predictive functional profiles among the solid, liquid, and dorsal, and ventral epithelium fractions also revealed significant differences. Microbiota in the ventral fraction of yak rumen also significantly differ from reported microbiota of cattle. In conclusion, our results improve our knowledge of the taxonomic composition and roles of yak rumen microbiota.

Variation in the Fatty Acid Synthase Gene (FASN) and Its Association with Milk Traits in Gannan Yaks

Fatty acid synthase (FASN) is an enzyme involved in the synthesis of fatty acids (FA) and plays a central role in de novo lipogenesis in mammals. This study was conducted to ascertain the relative level of expression of the FASN gene (FASN) in tissues from the yak (Bos grunniens), and to search for variation in two regions of yak FASN using polymerase chain reaction single-stranded conformational polymorphism (PCR-SSCP) analyses; it also ascertains whether that variation is associated with yak milk traits. The gene was found to be expressed in twelve tissues, with the highest expression detected in the mammary gland, followed by subcutaneous fat tissue. Two regions of the gene were analyzed in 290 Gannan yaks: A region spanning exon 24-intron 24 and a region spanning exon 34. These regions both produced two PCR-SSCP patterns, which, upon sequencing, represented different DNA sequences. This sequence variation resulted from the presence of three nucleotide substitutions: c.4296+38C/T (intron 24), c.5884A/G, and c.5903G/A, both located in exon 34. The exon 34 substitutions would result in the amino acid substitutions p.Thr1962Ala and p.Gly1968Glu if expressed. Four haplotypes spanning from the exon 24-intron 24 region to exon 34 were identified. Of these, two were common (A₁-A₂ and B₁-A₂), and two were rare (A₁-B₂ and B₁-B₂) in the yaks investigated. The presence of A₁-A₂ was associated with an increase in milk fat content (p = 0.050) and total milk solid content (p = 0.037), while diplotype A₁-A₂/B₁-A₂ had a higher milk fat content (p = 0.038) than the other diplotypes. This study suggests that further characterization of the FASN gene might provide for an improved understanding of milk traits in yaks.

Characterization of Yak Common Biofluids Metabolome by Means of Proton Nuclear Magnetic Resonance Spectroscopy

The aim of this study was to evaluate the metabolic profiles of yak (Bos grunniens) serum, feces, and urine by using proton nuclear magnetic resonance (¹H-NMR), to serve as a reference guide for the healthy yak milieu. A total of 108 metabolites, giving information about diet, protein digestion, and energy generation or gut-microbial co-metabolism, were assigned across the three biological matrices. A core metabolome of 15 metabolites was ubiquitous across all biofluids. Lactate, acetate, and creatinine could be regarded as the most abundant metabolites in the metabolome of serum, feces, and urine, respectively. Metabolic pathway analysis showed that the molecules identified could be able to give thorough information about four main metabolic pathways, namely valine, leucine, and isoleucine biosynthesis; phenylalanine, tyrosine, and tryptophan biosynthesis; glutamine and glutamate metabolism; and taurine and hypotaurine metabolism.

A global analysis of CNVs in diverse yak populations using whole-genome resequencing

Background

Genomic structural variation represents a source for genetic and phenotypic variation, which may be subject to selection during the environmental adaptation and population differentiation. Here, we described a genome-wide analysis of copy number variations (CNVs) in 16 populations of yak based on genome resequencing data and CNV-based cluster analyses of these populations.

Results

In total, we identified 51,461 CNV events and defined 3174 copy number variation regions (CNVRs) that covered 163.8 Mb (6.2%) of yak genome with more “loss” events than both “gain” and “both” events, and we confirmed 31 CNVRs in 36 selected yaks using quantitative PCR. Of the total 163.8 Mb CNVR coverage, a 10.8 Mb region of high-confidence CNVRs directly overlapped with the 52.9 Mb of segmental duplications, and we confirmed their uneven distributions across chromosomes. Furthermore, functional annotation indicated that the CNVR-harbored genes have a considerable variety of molecular functions, including immune response, glucose metabolism, and sensory perception. Notably, some of the identified CNVR-harbored genes associated with adaptation to hypoxia (e.g., DCC, MRPS28, GSTCD, MOGAT2, DEXI, CIITA, and SMYD1). Additionally, cluster analysis, based on either individuals or populations, showed that the CNV clustering was divided into two origins, indicating that some yak CNVs are likely to arisen independently in different populations and contribute to population difference.

Conclusions

Collectively, the results of the present study advanced our understanding of CNV as an important type of genomic structural variation in yak, and provide a useful genomic resource to facilitate further research on yak evolution and breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5451-5) contains supplementary material, which is available to authorized users.

Identification and diversity of Y-chromosome haplotypes in Qinghai yak populations

The aim of the present study was to perform a preliminary analysis of the characterization and diversity of Y-chromosome haplotypes/haplogroups in yak of Qinghai Province, China. A total of 322 male yaks from nine populations belonging to three officially recognized breeds (Gaoyuan, Huanhu and Datong) were sampled. Animals were genotyped using six previously reported Y-SNPs present in the SRY, USP9Y, UTY, AMELY and OFD1Y genes and four new Y-SNPs in the OFD1Y gene (g.569A>C, g.578A>C, g.608G>T and g.653G>C) identified in this study. Seven Y-chromosome haplotypes (H1-H7) were identified according to the combination of the 10 Y-SNPs. H1, H2 and H6 were the most common and shared haplotypes across all yak populations/breeds. Private haplotypes H3 and H7 were detected in the Datong breed; H4 in Guoleimude, Qumalai, Qilian, Tianjun and Ganglong populations; and H5 in Qumalai of Gaoyuan breed. Haplotype clustering and network analyses inferred two haplogroups, Y1 and Y2, indicating two divergent lineages of paternal origins of Qinghai yak. The analysis of molecular variance showed a significant difference among individuals (P < 0.0001) with more than 93% of the total genetic variation present within populations, suggesting a weak genetic structure among Qinghai yak populations. The overall Y-haplotype diversity was 0.538 ± 0.028, showing a relatively high diversity in Qinghai yak. The Gaoyuan and Datong breeds had similar haplotype diversities (0.547 ± 0.030 and 0.553 ± 0.083, respectively), which were higher than that of the Huanhu breed (0.441 ± 0.098). Our results support the conservation and sustainable use of unique yak genetic resources in Qinghai.

Evaluation of foot and mouth vaccination for yak (Bos grunniens) in Pakistan

In northern Pakistan, many farming communities rely on domestic yak (Bos grunniens) as a principle source of income. A 2006 participatory disease surveillance report from this region indicated that foot-and-mouth disease (FMD) is the most prevalent annual disease of yak. Our objectives of this study were to determine exposure levels of yak to FMD virus; implement a vaccination program based on current, regional FMD virus serotypes and subtypes; and quantify immune responses following vaccination. Blood samples were used to determine pre-vaccination exposure of animals to FMD virus by antibody presence to non-structural proteins of FMD virus using a 3-ABC trapping indirect ELISA. Vaccine used consisted of FMD serotypes 'O' (PanAsia-2), 'A' (Iran-05), and 'Asia-1' (Shamir), but changed later during the study to match newly circulating viruses in the country ('O'-PanAsia-2; 'A'-Turk-06 and Asia-1-Sindh-08). Three hundred sixty-three blood samples were tested from selected villages to determine pre-vaccination FMD virus exposure in yak with an average of 37.7%. Immune responses from initial vaccination and booster dose 30 days later showed clear protective levels (as mean percent inhibition) of antibodies against structural proteins of serotypes 'O,' 'A,' and 'Asia-1.' These responses remained above threshold positive level even at day 210 following initial vaccination. Results of sero-surveillance and anecdotal information of repeated FMD outbreaks demonstrate the persistence of FMD virus of yak in northern Pakistan. Laboratory results and field observations clearly indicated that yak can be protected against FMD with a good quality vaccine with FMD serotype(s) matching current, regionally circulating FMD virus.

Seroprevalence and Risk Factors of Fascioliasis in Yaks, Bos grunniens, from Three Counties of Gansu Province, China

The aim of this study was to determine the seroprevalence and risk factors of fascioliasis in yaks, Bos grunniens, from 3 counties of Gansu Province in China. A total of 1,584 serum samples, including 974 samples from white yaks from Tianzhu, 464 from black yaks from Maqu, and 146 from black yaks from Luqu County, were collected and analyzed using ELISA to detect IgG antibodies against Fasciola hepatica. The overall F. hepatica seroprevalence was 28.7% (454/1,584), with 29.2% in white yaks (284/974) and 27.9% in black yaks (170/610). The seroprevalence of F. hepatica in yaks from Tianzhu, Luqu, and Maqu was 29.2%, 22.6%, and 29.5%, respectively. Female yaks (30.9%) had higher F. hepatica seroprevalence than male yaks (23.4%). Also, F. hepatica seroprevalence varied by different age group from 24.1% to 33.8%. Further, the seroprevalence ranged from 21.8% to 39.1% over different seasons. Interestingly, the season and age of yaks were associated with F. hepatica infection in yaks in the investigated areas. These findings provided a basis for further studies on this disease in yaks from 3 counties of Gansu Province in northwestern China, which may ultimately support the development of effective control strategies of fascioliasis in these areas.

Protein and mRNA expression of follicle-stimulating hormone receptor and luteinizing hormone receptor during the oestrus in the yak (Bos grunniens)

Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) have a central role in follicle growth, maturation and oestrus, but no clear pathway in the seasonal oestrus of yak (Bos grunniens) has been found. To better understand the role of FSH and LH in seasonal oestrus in the yak, six yaks were slaughtered while in oestrus, and the pineal gland, hypothalamus, pituitary gland, and gonads were collected. Using real-time PCR and immunohistochemical assays, we determined the mRNA and protein expression of the FSH and LH receptors (FSHR and LHR) in these organs. The analysis showed that the FSHR mRNA expression level was higher in the pituitary gland tissue compared with LHR (p < .01) during oestrus. By contrast, there was low expression of FSHR and LHR mRNA in the pineal gland and hypothalamus. FSHR mRNA expression was higher than that of LHR (p < .05) in the ovary, whereas LHR mRNA expression was higher than that of FSHR (p < .01) in the uterus. FSHR and LHR proteins were located in the pinealocyte, synaptic ribbon and synaptic spherules of the pineal gland and that FSH and LH interact via nerve fibres. In the hypothalamus, FSHR and LHR proteins were located in the magnocellular neurons and parvocellular neurons. FSHR and LHR proteins were localized in acidophilic cells and basophilic cells in the pituitary gland, and in surface epithelium, stromal cell and gland epithelium in the uterus. In the ovary, FSHR and LHR protein were present in the ovarian follicle. Thus, we concluded that FSHR and LHR are located in the pineal gland, hypothalamus, pituitary and gonad during oestrus in the yak. However, FSHR was mainly expressed in the pituitary gland and ovaries, whereas LHR was mainly expressed in the pituitary gland and uterus.

Prevalence of Circulating Antibodies to Bovine Herpesvirus 1 in Yaks (Bos grunniens) on the Qinghai-Tibetan Plateau, China

Bovine Herpesvirus 1 (BoHV-1) causes infections with many clinical signs, including rhinotracheitis, encephalitis, and genital lesions. The virus occurs worldwide in bovines, and in recent years, it has been reported in yaks (Bos grunniens) inhabiting the Tibetan Plateau in China. However, there is little epidemiologic data describing BoHV-1 infections in China's yak herds. We conducted a cross-sectional study on the Qinghai-Tibetan Plateau (QTP) in China July 2011-July 2012 to estimate the prevalence of BoHV-1 antibody in yak herds. We collected 1,840 serum samples from yaks on the QTP, in Tibet (988 yaks), Qinghai (475 yaks), and Sichuan (377 yaks) Provinces. Using an enzyme-linked immunosorbent assay, we found that 381 (38.6%) of the Tibetan samples, 212 (44.6%) of the Qinghai samples, and 105 (27.9%) of the Sichuan samples had detectable antibodies to BoHV-1. Given that this high prevalence of infection in yaks could result in heavy economic losses, we suggest that an effective management program, including vaccination and strategies for infection control, be developed.

The complete mitochondrial genome of the Qinghai Plateau yak Bos grunniens (Cetartiodactyla: Bovidae)

The Qinghai Plateau yak Bos grunniens (Cetartiodactyla: Bovidae) is an important primitive local breed in the Qinghai-Tibetan Plateau and adjacent regions. In this study, its complete mitochondrial genome sequence has been assembled and characterized from high-throughput Illumina sequencing data. This genome is 16 322 bp in length, and contains 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a non-coding D-loop or control region. The nucleotide composition is asymmetric (33.73% A, 25.79% C, 13.19% G, and 27.29% T) with an overall A + T content of 61.02%. The gene arrangement and the composition are similar to most other vertebrates. These data would contribute to our better understanding its population genetics and evolutionary history.

Molecular Characterization and Expression Analysis of Insulin-like Growth Factor-1 and Insulin-like Growth Factor Binding Protein-1 Genes in Qinghai-Tibet Plateau Bos grunniens and Lowland Bos taurus

Insulin-like growth factor-1 (IGF-1) and insulin-like growth factor binding protein-1 (IGFBP-1) play a pivotal role in regulating cellular hypoxic response. In this study, we cloned and characterized the genes encoding IGF-1 and IGFBP-1 to improve the current knowledge on their roles in highland Bos grunniens (Yak). We also compared their expression levels in the liver and kidney tissues between yaks and lowland cattle. We obtained full-length 465 bp IGF-1 and 792 bp IGFBP-1, encoding 154 amino acids (AA) IGF-1, and 263 AA IGFBP-1 protein, respectively using reverse transcriptase-polyerase chain reaction (RT-PCR) technology. Analysis of their corresponding amino acid sequences showed a high identity between B. grunniens and lowland mammals. Moreover, the two genes were proved to be widely distributed in the examined tissues through expression pattern analysis. Real-time PCR results revealed that IGF-1 expression was higher in the liver and kidney tissues in B. grunniens than in Bos taurus (p<0.05). The IGFBP-1 gene was expressed at a higher level in the liver (p<0.05) of B. taurus than B. grunniens, but it has a similar expression level in the kidneys of the two species. These results indicated that upregulated IGF-1 and downregulated IGFBP-1 are associated with hypoxia adaptive response in B. grunniens.

Complete mitochondrial genome of the Yakow (Bos primigenius taurus × Bos grunniens) in China

The Chinese Yakow, Bos primigenius taurus × Bos grunniens, is a large and commercially important hybrid in family Bovidae. We first determined and annotated its complete mitochondrial genome. The mitogenome is 16,322 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.7%, T: 27.3%, C: 25.8% and G: 13.2%. The complete mitogenome of Yakow (B. p. taurus × B. grunniens) could provide an important data to further explore the taxonomic status of Yakow (B. p. taurus × B. grunniens) in B. grunniens and enrich the genetic information or evolutionary history of the Bovidae.

Complete mitochondrial genome of the a rare subspecies of genus Bos, Tianzhu white yak from Tibetan area in China

The population of domestic yak, Tianzhu white yak, from Tibetan area in China is considered as a rare Bos grunniens species. We first determined and annotated its complete mitochondrial genome. The mitogenome is 16,319 bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes and a control region. As in other mammals, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes, which are encoded on the light strand. Its overall base composition is A: 33.7%, T: 27.2%, C: 25.8% and G: 13.2%. The complete mitogenome of the new subspecies of Bos grunniens could provide an important data to further explore the taxonomic status of the subspecies.

Maternal phylogeny of a newly-found yak population in china

The Jinchuan yak is a new yak population identified in Sichuan, China. This population has a special anatomical characteristic: an additional pair of ribs compared with other yak breeds. The genetic structure of this population is unknown. In the present study, we investigated the maternal phylogeny of this special yak population using the mitochondrial DNA variation. A total of 23 Jinchuan yaks were sequenced for a 823-bp fragment of D-loop control region and three individuals were sequenced for the whole mtDNA genome with a length of 16,371-bp. To compare with the data from other yaks, we extracted sequence data from Genebank, including D-loop of 398 yaks (from 12 breeds) and 55 wild yaks, and whole mitochondrial genomes of 53 yaks (from 12 breeds) and 21 wild yaks. A total of 127 haplotypes were defined, based on the D-loop data. Thirteen haplotypes were defined from 23 mtDNA D-loop sequences of Jinchuan yaks, six of which were shared only by Jinchuan, and one was shared by Jinchuan and wild yaks. The Jinquan yaks were found to carry clades A and B from lineage I and clade C of lineage II, respectively. It was also suggested that the Jinchuan population has no distinct different phylogenetic relationship in maternal inheritance with other breeds of yak. The highly haplotype diversity of the Pali breed, Jinchuan population, Maiwa breed and Jiulong breed suggested that the yak was first domesticated from wild yaks in the middle Himalayan region and the northern Hengduan Mountains. The special anatomic characteristic that we found in the Jinchuan population needs further studies based on nuclear data.