GbyE: an integrated tool for genome widely association study and genome selection based on genetic by environmental interaction.

BACKGROUND: The growth and development of organism were dependent on the effect of genetic, environment, and their interaction. In recent decades, lots of candidate additive genetic markers and genes had been detected by using genome-widely association study (GWAS). However, restricted to computing power and practical tool, the interactive effect of markers and genes were not revealed clearly. And utilization of these interactive markers is difficult in the breeding and prediction, such as genome selection (GS). RESULTS: Through the Power-FDR curve, the GbyE algorithm can detect more significant genetic loci at different levels of genetic correlation and heritability, especially at low heritability levels. The additive effect of GbyE exhibits high significance on certain chromosomes, while the interactive effect detects more significant sites on other chromosomes, which were not detected in the first two parts. In prediction accuracy testing, in most cases of heritability and genetic correlation, the majority of prediction accuracy of GbyE is significantly higher than that of the mean method, regardless of whether the rrBLUP model or BGLR model is used for statistics. The GbyE algorithm improves the prediction accuracy of the three Bayesian models BRR, BayesA, and BayesLASSO using information from genetic by environmental interaction (G × E) and increases the prediction accuracy by 9.4%, 9.1%, and 11%, respectively, relative to the Mean value method. The GbyE algorithm is significantly superior to the mean method in the absence of a single environment, regardless of the combination of heritability and genetic correlation, especially in the case of high genetic correlation and heritability. CONCLUSIONS: Therefore, this study constructed a new genotype design model program (GbyE) for GWAS and GS using Kronecker product. which was able to clearly estimate the additive and interactive effects separately. The results showed that GbyE can provide higher statistical power for the GWAS and more prediction accuracy of the GS models. In addition, GbyE gives varying degrees of improvement of prediction accuracy in three Bayesian models (BRR, BayesA, and BayesCpi). Whatever the phenotype were missed in the single environment or multiple environments, the GbyE also makes better prediction for inference population set. This study helps us understand the interactive relationship between genomic and environment in the complex traits. The GbyE source code is available at the GitHub website ( https://github.com/liu-xinrui/GbyE ).

Whole-genome resequencing of major populations revealed domestication-related genes in yaks.

BACKGROUND: The yak is a symbol of the Qinghai-Tibet Plateau and provides important basic resources for human life on the plateau. Domestic yaks have been subjected to strong artificial selection and environmental pressures over the long-term. Understanding the molecular mechanisms of phenotypic differences in yak populations can reveal key functional genes involved in the domestication process and improve genetic breeding. MATERIAL AND METHOD: Here, we re-sequenced 80 yaks (Maiwa, Yushu, and Huanhu populations) to identify single-nucleotide polymorphisms (SNPs) as genetic variants. After filtering and quality control, remaining SNPs were kept to identify the genome-wide regions of selective sweeps associated with domestic traits. The four methods (π, XPEHH, iHS, and XP-nSL) were used to detect the population genetic separation. RESULTS: By comparing the differences in the population stratification, linkage disequilibrium decay rate, and characteristic selective sweep signals, we identified 203 putative selective regions of domestic traits, 45 of which were mapped to 27 known genes. They were clustered into 4 major GO biological process terms. All known genes were associated with seven major domestication traits, such as dwarfism (ANKRD28), milk (HECW1, HECW2, and OSBPL2), meat (SPATA5 and GRHL2), fertility (BTBD11 and ARFIP1), adaptation (NCKAP5, ANTXR1, LAMA5, OSBPL2, AOC2, and RYR2), growth (GRHL2, GRID2, SMARCAL1, and EPHB2), and the immune system (INPP5D and ADCYAP1R1). CONCLUSIONS: We provided there is an obvious genetic different among domestic progress in these three yak populations. Our findings improve the understanding of the major genetic switches and domestic processes among yak populations.

miR-129 Regulates Yak Intramuscular Preadipocyte Proliferation and Differentiation through the PI3K/AKT Pathway.

miR-129 plays a crucial role in regulating various cellular processes, including adipogenesis; however, its downstream molecular mechanisms remain unclear. In this study, we demonstrated that miR-129 promotes yak adipogenesis in vitro via the PI3K/AKT pathway. Overexpression and interference of miR-129 in yak intramuscular preadipocytes (YIMAs) enhanced and inhibited cell differentiation, respectively, with corresponding changes in cell proliferation. Further investigation revealed that miR-129 enhances AKT and p-AKT activity in the AKT pathway without affecting cell apoptosis, and a specific inhibitor (LY294002) was used to confirm that miR-129 regulates YIMAs proliferation and differentiation through the PI3K/AKT pathway. Our findings suggest that miR-129 promotes yak adipogenesis by enhancing PI3K/AKT pathway activity. This study provides the foundation to precisely elucidate the molecular mechanism of miR-129 in YIMAs adipogenesis and develop advanced miRNA-based strategies to improve meat nutrition and obesity-related ailments in beef production.

Single-cell transcriptome analysis and in vitro differentiation of testicular cells reveal novel insights into male sterility of the interspecific hybrid cattle-yak.

BACKGROUND: Interspecific hybridization plays vital roles in enriching animal diversity, while male hybrid sterility (MHS) of the offspring commonly suffered from spermatogenic arrest constitutes the postzygotic reproductive isolation. Cattle-yak, the hybrid offspring of cattle (Bos taurus) and yak (Bos grunniens) can serve as an ideal MHS animal model. Although meiotic arrest was found to contribute to MHS of cattle-yak, yet the cellular characteristics and developmental potentials of male germline cell in pubertal cattle-yak remain to be systematically investigated. RESULTS: Single-cell RNA-seq analysis of germline and niche cell types in pubertal testis of cattle-yak and yak indicated that dynamic gene expression of developmental germ cells was terminated at late primary spermatocyte (meiotic arrest) and abnormal components of niche cell in pubertal cattle-yak. Further in vitro proliferation and differentially expressed gene (DEG) analysis of specific type of cells revealed that undifferentiated spermatogonia of cattle-yak exhibited defects in viability and proliferation/differentiation potentials. CONCLUSION: Comparative scRNA-seq and in vitro proliferation analysis of testicular cells indicated that not only meiotic arrest contributed to MHS of cattle-yak. Spermatogenic arrest of cattle-yak may originate from the differentiation stage of undifferentiated spermatogonia and niche cells of cattle-yak may provide an adverse microenvironment for spermatogenesis.

Transcriptomic analysis of yak longissimus dorsi muscle identifies genes associated with tenderness.

Meat tenderness is an important sensory index when consumers choose meat products, which determines the value of meat products and consumers' buying intentions. Yak meat is rich in nutrition and unique in flavor, which is favored by consumers. However, its meat has the deficiencies of low tenderness and poor taste, which has a negative impact on the value of its meat products and customer satisfaction. To identify the genes affecting the yak meat tenderness, we used RNA-seq to analyze the longissimus dorsi muscle of yaks with different tenderness, screened a total of 1120 differentially expressed genes (DEGs). Meanwhile, 23 pathways were significantly enriched. By further analysis, we identified eight genes related to yak meat tenderness (WNT5A, ARID5B, SERPINE1 KLHL40, RUNX1, MAFF, RFX7 and ARID5A). Notably, SERPINE1 was involved in the significant enrichment pathways of 'complement and coagulation cascade pathway', 'HIF-1 signaling pathway' and 'AGE-RAGE signaling pathway in diabetic complications' which can regulate meat tenderness. This implies that SERPINE1 may play an important regulatory role among them. The DEGs associated with yak meat quality screened in this work will be helpful to identify potential biomarkers related to meat tenderness.

Epigenomics Analysis of the Suppression Role of SIRT1 via H3K9 Deacetylation in Preadipocyte Differentiation.

Sirtuin 1 (SIRT1) overexpression significantly inhibits lipid deposition during yak intramuscular preadipocyte (YIMA) differentiation; however, the regulatory mechanism remains unknown. We elucidated the role of SIRT1 in YIMA differentiation using lentivirus-mediated downregulation technology and conducted mRNA-seq and ChIP-seq assays using H3K9ac antibodies after SIRT1 overexpression in order to reveal SIRT1 targets during YIMA adipogenesis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed in order to identify the functional annotation of common genes. In addition, a potential target of SIRT1 was selected to verify its effects on the differentiation and proliferation of YIMAs. SIRT1 interfered with lipid deposition and promoted YIMA differentiation. In total, 143,518 specific peaks were identified after SIRT1 overexpression, where genes associated with downregulation peaks were enriched in transcription, gene expression, lipid-related processes, and classical lipid-related pathways. The H3K9ac signal in the whole genome promoter region (2 kb upstream and downstream of the transcription start site (TSS)) was weakened, and the peaks were distributed across all gene functional regions. Genes that lost signals in their TSS region or gene body region were enriched in both biological processes and pathways associated with lipogenesis. The ChIP-seq results revealed 714 common differential genes in mRNA-seq, which were enriched in "MAPK signaling", "lipid and atherosclerosis", "mTOR signaling", and "FoxO signaling" pathways. A total of 445 genes were downregulated in both their H3K9ac signals and mRNA expression, and one of their most significantly enriched pathways was FoxO signaling. Nine genes (FBP2, FPGT, HSD17B11, KCNJ15, MAP3K20, SLC5A3, TRIM23, ZCCHC10, and ZMYM1) lost the H3K9ac signal in their TSS regions and had low mRNA expression, and three genes (KCNJ15, TGM3, and TRIM54) had low expression but lost their H3K9ac signal in the gene body region. The interference of TRIM23 significantly inhibited fat deposition during preadipocyte differentiation and promoted cell proliferation by increasing S-phase cell numbers. The present study provides new insights into the molecular mechanism of intramuscular fat content deposition and the epigenetic role of SIRT1 in adipocyte differentiation.

Label-free measurement of wall shear stress in the brain venule and arteriole using dual-wavelength third-harmonic-generation line-scanning imaging.

Wall shear stress (WSS) is of fundamental physiological and pathological significance. Current measurement technologies suffer from poor spatial resolution or cannot measure instantaneous values in a label-free manner. Here we demonstrate dual-wavelength third-harmonic-generation (THG) line-scanning imaging, for instantaneous wall shear rate and WSS measurement in vivo. We used the soliton self-frequency shift to generate dual-wavelength femtosecond pulses. Simultaneous acquisition of dual-wavelength THG line-scanning signals extract blood flow velocities at adjacent radial positions for instantaneous wall shear rate and WSS measurement. Our results show the oscillating behavior of WSS in brain venules and arterioles at micron spatial resolution in a label-free manner.

Dysregulated genes in undifferentiated spermatogonia and Sertoli cells are associated with the spermatogenic arrest in cattleyak.

Hybrid male sterility (HMS) is a reproductive isolation mechanism limiting the formation of fertile offspring through interspecific fertilization. Cattleyak is the interspecific hybrid presenting significant heterosis in several economic traits, but HMS restricted its wide reproduction in cettleyak breeding. In this study, we detected the specifically expressed genes of a variety of cells (undifferentiated spermatogonia, primary spermatocytes, secondary spermatocytes, haploid spermatids, sperm, Sertoli cells, Leydig cells, and macrophages) in the testis of yak and cattleyak, and found that the spermatogenesis of cattleyak might be blocked at meiosis I, and the expression of niche factors (NR5A1, GATA4, STAR, CYP11A1, CD68, TNF, and CX3CR1) in undifferentiated spermatogonia niche was abnormal. Then we isolated the undifferentiated spermatogonia and Sertoli cells from yak and cattleyak by enzyme digestion, and detected the specific genes in the two bovid testicular cells as well as the proliferation capacity of the undifferentiated spermatogonia. These results indicated that weak proliferation ability and scarce number of undifferentiated spermatogonia and abnormal gene expressions in Sertoli cells may contribute to male sterility of cattleyak.

Codon Usage Bias and Cluster Analysis of the MMP-2 and MMP-9 Genes in Seven Mammals.

Matrix metalloproteinase (MMP)-2 and MMP-9 are a family of Zn2+ and Ca2+-dependent gelatinase MMPs that regulate muscle development and disease treatment, and they are highly conservative during biological evolution. Despite increasing knowledge of MMP genes, their evolutionary mechanism for functional adaption remains unclear. Moreover, analysis of codon usage bias (CUB) is reliable to understand evolutionary associations. However, the distribution of CUB of MMP-2 and MMP-9 genes in mammals has not been revealed clearly. Multiple analytical software was used to study the genetic evolution, phylogeny, and codon usage pattern of these two genes in seven species of mammals. Results showed that the MMP-2 and MMP-9 genes have CUB. By comparing the content of synonymous codon bases amongst seven mammals, we found that MMP-2 and MMP-9 were low-expression genes in mammals with high codon conservation, and their third codon preferred the G/C base. RSCU analysis revealed that these two genes preferred codons encoding delicious amino acids. Analysing what factors influence CUB showed that the third base distributors of these two genes were C/A and C/T, and GC3S had a wide distribution range on the ENC plot reference curve under no selection or mutational pressure. Thus, mutational pressure is an important factor in CUB. This study revealed the usage characteristics of the MMP-2 and MMP-9 gene codons in different mammals and provided basic data for further study towards enhancing meat flavour, treating muscle disease, and optimizing codons.

Angiogenesis-Inhibitory Piperidine Alkaloids from the Leaves of Microcos paniculata.

Eight new 2,6-disubstituted piperidin-3-ol alkaloids (1-8), featuring a C10 unsaturated alkyl side chain, together with three previously reported analogues (9-11) were isolated from the leaves of medicinal plant Microcos paniculata. Their structures and absolute configurations were elucidated unambiguously by means of 1D and 2D NMR spectroscopic data analysis, modified Mosher's method, Snatzke's method, and quantum chemical electronic circular dichroism (ECD) calculations, as well as single-crystal X-ray crystallography. The isolates were evaluated for their antiangiogenic effects on human umbilical vein endothelial cells (HUVECs). Compound 2 displayed an inhibitory effect on tube formation of HUVECs in a concentration-dependent manner.

Comparative transcriptome analysis of winter yaks in plateau and plain.

The yak is an important source for the people living and ecological environment in the Qinghai-Tibet Plateau. In every winter, many domestic yaks will lose bodyweight or dead under cold and food scarcity. Moving the plateau yaks to farm in the plain is a useful approach to reduce their environmental stress and gain more production. In this study, we measured growth, slaughter and beef quality traits every month and sequenced mRNA expression levels of muscles of two groups yaks living in plateau and plain respectively. We found there is significant difference (p-value <0.01) in the third (60 days), fourth (90 days), fifth (120 days) and sixth (150 days) weights between subpopulations in the plateau and plain. We identified 540 different expressed genes (DEGs) including 123 known genes and 417 unknown genes. Using the weighted correlation network analysis (WGCNA) to build a co-express network, the modules were strong relative to weight traits. The findings highlighted the underlying way and a relative network to yield a new view about gene expression between the yaks living plateau and plain.

Adipose-Secreted Exosomes and Their Pathophysiologic Effects on Skeletal Muscle.

Due to its prominent secretory activity, adipose tissue (AT) is now considered a major player in the crosstalk between organs, especially with skeletal muscle. In which, exosomes are effective carriers for the intercellular material transfer of a wide range of molecules that can influence a series of physiological and pathological processes in recipient cells. Considering their underlying roles, the regulatory mechanisms of adipose-secreted exosomes and their cellular crosstalk with skeletal muscle have received great attention in the field. In this review, we describe what is currently known of adipose-secreted exosomes, as well as their applications in skeletal muscle pathophysiology.

Differences in proteomic profiles between yak and three cattle strains provide insights into molecular mechanisms underlying high-altitude adaptation.

Yaks display unique properties of the lung and heart, enabling their adaptation to high-altitude environments, but the underlying molecular mechanisms are still largely unknown. In the present study, the proteome differences in lung and heart tissues were compared between yak (Bos grunniens) and three cattle strains (Bos taurus, Holstein, Sanjiang and Tibetan cattle) using the sequential window acquisition of all theoretical mass spectra/data-independent acquisition (SWATH/DIA) proteomic method. In total, 51,755 peptides and 7215 proteins were identified. In the lung tissue, there were 162, 310 and 118 differential abundance proteins (DAPs) in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. In the heart tissue, there were 71, 57 and 78 DAPs in Tibetan, Holstein and Sanjiang cattle compared to yak respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DAPs were enriched for the retinol metabolism and toll-like receptor categories in lung tissue. The changes in these two pathways may regulate hypoxia-induced factor and immune function in yaks. Moreover, DAPs in heart tissues were enriched for cardiac muscle contraction, Huntington's disease, chemical carcinogenesis and drug metabolism-cytochrome P450. Further exploration indicated that yaks may alter cardiac function through regulation of type 2 ryanodine receptor (RyR2) and Ca2+ -release channels. The present results are useful to further develop an understanding of the mechanisms underlying adaptation of animals to high-altitude conditions.

Testis transcriptome profiling identified lncRNAs involved in spermatogenic arrest of cattleyak.

Cattleyaks are the crossbred offspring between cattle and yaks, exhibiting the prominent adaptability to the harsh environment as yaks and much higher growth performances than yaks around Qinghai-Tibet plateau. Unfortunately, cattleyak cannot be effectively used in yak breeding due to its male infertility resulted from spermatogenic arrest. In this study, we performed RNA sequencing (RNA-seq) and bioinformatics analysis to determine the expression profiles of long noncoding RNA (lncRNA) from cattleyak and yak testis. A total of 604 differentially expressed (DE) lncRNAs (135 upregulated and 469 downregulated) were identified in cattleyak with respect to yak. Through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, we identified several DE lncRNAs regulating the mitotic cell cycle processes by targeting the genes significantly associated with the mitotic cell cycle checkpoint and DNA damage checkpoint term and also significantly involved in p53 signaling pathway, mismatch repair and homologous recombination pathway (P < 0.05). The reverse transcription PCR (RT-PCR) and quantitative Real-Time PCR (qRT-PCR) analysis of the randomly selected fourteen DE lncRNAs and the seven target genes validated the RNA-seq data and their true expressions during spermatogenesis in vivo. Molecular cloning and sequencing indicated that the testis lncRNAs NONBTAT012170 and NONBTAT010258 presented higher similarity among different cattleyak and yak individuals. The downregulation of these target genes in cattleyak contributed to the abnormal DNA replication and spermatogenic arrest during the S phase of mitotic cell cycle. This study provided a novel insight into lncRNA expression profile changes associated with spermatogenic arrest of cattleyak.

Whole-genome resequencing provides insights into the evolution and divergence of the native domestic yaks of the Qinghai-Tibet Plateau.

BACKGROUND: On the Qinghai-Tibet Plateau, known as the roof ridge of the world, the yak is a precious cattle species that has been indispensable to the human beings living in this high-altitude area. However, the origin of domestication, dispersal route, and the divergence of domestic yaks from different areas are poorly understood. RESULTS: Here, we resequenced the genome of 91 domestic yak individuals from 31 populations and 1 wild yaks throughout China. Using a population genomics approach, we observed considerable genetic variation. Phylogenetic analysis suggested that the earliest domestications of yak occurred in the south-eastern QTP, followed by dispersal to the west QTP and northeast to SiChuang, Gansu, and Qinghai by two routes. Interestingly, we also found potential associations between the distribution of some breeds and historical trade routes such as the Silk Road and Tang-Tibet Ancient Road. Selective analysis identified 11 genes showing differentiation between domesticated and wild yaks and the potentially positively selected genes in each group were identified and compared among domesticated groups. We also detected an unbalanced pattern of introgression among domestic yak, wild yak, and Tibetan cattle. CONCLUSIONS: Our research revealed population genetic evidence for three groups of domestic yaks. In addition to providing genomic evidence for the domestication history of yaks, we identified potential selected genes and introgression, which provide a theoretical basis and resources for the selective breeding of superior characters and high-quality yak.

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.

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.

Signature of high altitude adaptation in the gluteus proteome of the yak.

Yak is the unique Bovidae species in the Qinghai-Tibetan Plateau. A previous proteomic study has compared the yak muscle tissue to one cattle strain using the isobaric tags for relative and absolute quantification approach. In this study, to further investigate the molecular mechanisms underlying yak adaptation, the proteomic profiles of gluteus were compared between yak and one moderate-altitude cattle strain (Tibetan cattle) and two low-altitude cattle strains (Holstein and Sanjiang cattle) using a label-free quantitative method. The comparisons identified 20, 364, 143 upregulated proteins and 4, 6, 37 downregulated proteins in yak, compared with Tibetan, Holstein, and Sanjiang cattle, respectively. Protein-protein interaction analysis indicated that these differentially expressed proteins were mainly related to "oxidative phosphorylation" and "electron transport chain." Further analysis revealed that NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 11, NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 4, cytochrome C oxidase subunit 6A2, mitochondrial and cytochrome c oxidase subunit NDUFA4 were all increased in the yak, suggesting that yak might increase mitochondrial capacity to sustain metabolic rates under high altitude conditions, which might be a long-term adaptive mechanism underlying adaptation to high altitude environments. Yak increased the level of thioredoxin reductase 2 to protect themselves from oxidative damages. Moreover, the increased expression levels of phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit alpha isoform and caveolin-1 in yak suggested that yaks promoted glucose uptake for adaptation to high altitude. These results provided more information to better understand the molecular mechanisms underlying yak adaption.

Detection and integrated analysis of lncRNA and mRNA relevant to plateau adaptation of Yak.

Known as the 'ship of the plateau', through thousands of years evolution and cruelty environments selection containing low oxygen and strong ultraviolet radiation, yaks have adapted plateau environments and supplied important goods and materials for the people in the Qinghai-Tibet Plateau. This study aimed to identify differentially expressed (DE) genes and novel long non-coding RNAs (lncRNAs) of yaks for the Plateau adaptation and their underlying co-expression and regulatory network. We carried out RNA-seq analysis for cerebral and cerebellar tissue specimens of Bos taurus, Bos grunniens × Bos Taurus and B. grunniens. Furthermore, 12,072 pseudo lncRNAs were predicted using three software. In total, 4,257 significant DE transcripts were identified using the Ballgown R package (p < .01), of which 1,021 were protein-coding genes, 14 were known lncRNAs, and 661 were novel lncRNAs. Using WGCNA, a co-expression network of DE mRNAs and lncRNAs comprising 5 modules was generated to determine functional associations clusters. This study reveals a valuable sub-network comprising 8 hub genes, one known lncRNA and 5 novel lncRNAs in the major module. These hub genes are associated with blood pressure regulation, generation of reactive oxygen species and metabolism. The analysis of co-expressed genes thus provides a basis for the regulatory mechanisms in PA in Yaks and for the detection of additional genes between cross-breed and parent populations.

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.