Air pollution, genetic factors and the risk of osteoporosis: A prospective study in the UK biobank.

Purpose: To reveal relationship between air pollution exposure and osteoporosis (OP) risk. Methods: Based on large-scale data from the UK Biobank, we evaluated the relationship between OP risk and several air pollutants. Then air pollution scores (APS) were constructed to assess the combined effects of multiple air pollutants on OP risk. Finally, we constructed a genetic risk score (GRS) based on a large genome-wide association study of femoral neck bone mineral density and assessed whether single or combined exposure to air pollutants modifies the effect of genetic risk on OP and fracture risk. Results: PM2.5, NO2, NOx, and APS were significantly associated with an increased risk of OP/fracture. OP and fracture risk raised with increasing concentrations of air pollutants: compared to the lowest APS quintile group, subjects in the highest quintile group had a hazard ratio (HR) (95% CI) estimated at 1.140 (1.072-1.213) for OP and 1.080 (1.026-1.136) for fracture. Moreover, participants with low GRS and the highest air pollutant concentration had the highest risk of OP, the HRs (95% CI) of OP were 1.706 (1.483-1.964), 1.658 (1.434-1.916), 1.696 (1.478-1.947), 1.740 (1.506-2.001) and 1.659 (1.442-1.908), respectively, for PM2.5, PM10, PM2.5-10, NO2, and NOx. Similar results were also observed for fractures. Finally, we assessed the joint effect of APS and GRS on the risk of OP. Participants with higher APS and lower GRS had a higher risk of developing OP. Similar results were observed in the joint effect of GRS and APS on fracture. Conclusions: We found that exposure to air pollution, individually or jointly, could improve the risk of developing OP and fractures, and increased the risk by interacting with genetic factors.

Genome-wide comparison of DNA methylation patterns between yak and three cattle strains and their potential association with mRNA transcription.

Yak has evolved specific adaptative mechanisms to high-altitude environment. Up to date, only a few studies reported the DNA methylation in yak. In the present study, genome-wide DNA methylome and transcriptome profiles in lung, mammary, and biceps brachii muscle tissues were compared between yak and three cattle breeds (Tibetan cattle, Sanjiang cattle, and Holstein cattle). The association between differentially expressed genes (DEGs) and differentially methylated regions (DMRs) was analyzed, and the biological functions of DEGs potentially driven by DMRs were explored by KEGG enrichment analysis. Finally, we found that yak-specific DMRs-driven DEGs were mainly involved in neuromodulation, respiration, lung development, blood pressure regulation, cardiovascular protection, energy metabolism, DNA repair, and immune functions. The higher levels of the key genes associated with these functions were observed in yak than in cattle, suggesting that DNA methylation might regulate these genes. Overall, the present study contributes basic data at the DNA methylation level to further understand the physiological metabolism in yak.

Genome-wide identification of SNPs associated with body weight in yak.

BACKGROUND: The yak is the most important livestock in the Qinghai-Tibet Plateau, and body weight directly affects the economic values of yak. Up to date, the genome-wide profiling of single-nucleotide polymorphisms (SNPs) associating with body weight has not been reported in yak. In the present study, the SNPs in 480 yaks from three breeds were analyzed using the commercial high-density (600 K) yak SNP chips. RESULTS: The results identified 12 and 4 SNPs potentially associated with body weight in male and female yaks, respectively. Among them, 9 and 2 SNPs showed significant difference in yak body weight between different genotypes at each locus in male and female yaks, respectively. Further exploration found 33 coding genes within the 100 kbp upstream or downstream to the SNP loci, which might be potentially affected by the variation of SNPs. Among them, G protein-coupled receptor kinase 4 (GRK4) might be potentially affected by the SNP AX-174555047, which has been reported to affect the functioning of two body-weight associated hormones (parathyroid hormone, PTH, and adrenomedullin, ADM). Determination of PTH and ADM levels in yak revealed positive relationship between PTH level and body weight, negative relationship between ADM level and body weight along with the variation of AX-174555047 mutation. CONCLUSIONS: These results suggested that the SNP AX-174555047 might potentially affect body weight through mediating GRK4 expression and then PTH and ADM functioning. Thus, the SNP AX-174555047 might be used as a biomarker for molecular breeding of yak. More investigations are required to validate this point.

Dynamic alterations in yak (Bos grunniens) rumen microbiome in response to seasonal variations in diet.

Rumen microorganisms play important roles in the healthy growth of yaks. This study investigated changes in yak rumen microbiome during natural grazing at the warm seasons and supplementary feeding at cold seasons. High-throughput sequencing of 16S rRNA and metagenome analysis were conducted to investigate the structures and functions of yak rumen microbial communities. The results indicated that Bacteroidetes and Firmicutes were the most abundant phyla. In addition, Bacteroidetes might play a more important role than Firmicutes during the supplementary feeding stage (spring and winter), but less during natural grazing stage (summer and autumn). KEGG analysis showed that the amino sugar and nucleotide sugar metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, starch and sucrose metabolism, and fructose and mannose metabolism were the main pathways in the microbial community, which were significantly different between seasons. The carbohydrate-active enzymes (CAZyme) annotation revealed that cellulose was an important carbon source for microorganisms in yak rumen. Glycoside hydrolases (GHs) were the most abundant class of CAZymes, followed by glycosyl transferases (GTs), which were important to digestion of oil, cellulose, and hemicellulose in food. These results contribute to the understanding of microbial components and functions in yak rumen.

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.

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.

Transcriptome analysis identified long non-coding RNAs involved in the adaption of yak to high-altitude environments.

The mechanisms underlying yak adaptation to high-altitude environments have been investigated using various methods, but no report has focused on long non-coding RNA (lncRNA). In the present study, lncRNAs were screened from the gluteus transcriptomes of yak and their transcriptional levels were compared with those in Sanjiang cattle, Holstein cattle and Tibetan cattle. The potential target genes of the differentially expressed lncRNAs between species/strains were predicted using cis and trans models. Based on cis-regulated target genes, no KEGG pathway was significantly enriched. Based on trans-regulated target genes, 11 KEGG pathways in relation to energy metabolism and three KEGG pathways associated with muscle contraction were significantly enriched. Compared with cattle strains, transcriptional levels of acyl-CoA dehydrogenase, acyl-CoA-binding protein, 3-hydroxyacyl-CoA dehydrogenase were relatively higher and those of glyceraldehyde 3-phosphate dehydrogenase, phosphoglycerate mutase 1, pyruvate kinase and lactate/malate dehydrogenase were relatively lower in yak, suggesting that yaks activated fatty acid oxidation but inhibited glucose oxidation and glycolysis. Besides, NADH dehydrogenase and ATP synthase showed lower transcriptional levels in yak than in cattle, which might protect muscle tissues from deterioration caused by reactive oxygen species (ROS). Compared with cattle strains, the higher transcriptional level of glyoxalase in yak might contribute to dicarbonyl stress resistance. Voltage-dependent calcium channel/calcium release channel showed a lower level in yak than in cattle strains, which could reduce the Ca2+ influx and subsequently decrease the risk of hypertension. However, levels of EF-hand and myosin were higher in yak than in cattle strains, which might enhance the negative effects of reduced Ca2+ on muscle contraction. Overall, the present study identified lncRNAs and proposed their potential regulatory functions in yak.

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.

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.

Comparisons of lung and gluteus transcriptome profiles between yaks at different ages.

The yak, Bos grunniens, is the only large mammal in the Qinghai-Tibet Plateau and has been bred to provide meat, milk, and transportation. Previous studies indicate that the immune system contributes to the yak's adaptation to high-altitude environments. In order to further investigate changes in immune function during yak development, we compared the transcriptome profiles of gluteus and lung tissues among yaks at 6, 30, 60, and 90 months of age. Analyses of significantly differentially expressed genes (DEGs) in lung tissues revealed that immune function was more activated at 6-months and less activated at 90-months than in the 30 and 60-month-old animals. DEG exploration in gluteal tissues revealed that immune functions were more highly activated at both 6 and 90-months, compared with 30 and 60-months. Immune system activation in the muscle and lung tissues of 30-month-old yaks may increase their resistance to infections, while decreased may be due to aging. Furthermore, the higher immune activation status in the gluteal tissues in 90-month-old yaks could be due to muscle injury and subsequent regeneration, which is supported by the fact that 5 unigenes related with muscle injury and 3 related to muscle regeneration displayed greater expression levels at 90-months than at 30 and 60-months. Overall, the present study highlights the important role of the immune system in yak development, which will facilitate future investigations.

Transcriptome profiles revealed the mechanisms underlying the adaptation of yak to high-altitude environments.

The yak is a valuable species in the Qinghai-Tibet Plateau of China. Nevertheless, the molecular mechanisms underlying its adaptation to high-altitude environments remain largely unknown. In the present study, comparative transcriptome sequencing was performed for lung and gluteus tissues from two species of low-altitude cattle (Sanjiang and Holstein cattle), Tibetan cattle (living at a moderate altitude), and yak (living at a high altitude) and the differentially expressed genes were validated using real-time quantitative PCR. The results showed that CD36 antigen was up-regulated and CD59 antigen was down-regulated in yak in comparison to the other animals, which might promote the development of red blood cells and inhibit the development of lymphocytes in yak. In addition, thrombospondin type 1, coagulation factor 5/8, and fibronectin were all down-regulated, but serpin and alpha 2-macroglobulin (A2M) were up-regulated. These differences would inhibit blood coagulation, thus reducing the risk of pulmonary edema. The expression levels of the calcium-release, potassium, and transient receptor potential channels decreased in yak, minimizing membrane depolarization and the harmful effects of pulmonary edema. Eleven KEGG pathways associated with innate immunity were more activated in yak and Tibetan cattle than in other cattle strains, which should reduce their risk of infection and disease. These changes together might facilitate the adaptation of yak and Tibetan cattle to live in high-altitude habitats.