Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Oceanic climate changes threaten the sustainability of Asia's water tower.

Water resources sustainability in High Mountain Asia (HMA) surrounding the Tibetan Plateau (TP)-known as Asia's water tower-has triggered widespread concerns because HMA protects millions of people against water stress1,2. However, the mechanisms behind the heterogeneous trends observed in terrestrial water storage (TWS) over the TP remain poorly understood. Here we use a Lagrangian particle dispersion model and satellite observations to attribute about 1 Gt of monthly TWS decline in the southern TP during 2003-2016 to westerlies-carried deficit in precipitation minus evaporation (PME) from the southeast North Atlantic. We further show that HMA blocks the propagation of PME deficit into the central TP, causing a monthly TWS increase by about 0.5 Gt. Furthermore, warming-induced snow and glacial melt as well as drying-induced TWS depletion in HMA weaken the blocking of HMA's mountains, causing persistent northward expansion of the TP's TWS deficit since 2009. Future projections under two emissions scenarios verified by satellite observations during 2020-2021 indicate that, by the end of the twenty-first century, up to 84% (for scenario SSP245) and 97% (for scenario SSP585) of the TP could be afflicted by TWS deficits. Our findings indicate a trajectory towards unsustainable water systems in HMA that could exacerbate downstream water stress.

Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

Changes in above- versus belowground biomass distribution in permafrost regions in response to climate warming.

Permafrost regions contain approximately half of the carbon stored in land ecosystems and have warmed at least twice as much as any other biome. This warming has influenced vegetation activity, leading to changes in plant composition, physiology, and biomass storage in aboveground and belowground components, ultimately impacting ecosystem carbon balance. Yet, little is known about the causes and magnitude of long-term changes in the above- to belowground biomass ratio of plants (η). Here, we analyzed η values using 3,013 plots and 26,337 species-specific measurements across eight sites on the Tibetan Plateau from 1995 to 2021. Our analysis revealed distinct temporal trends in η for three vegetation types: a 17% increase in alpine wetlands, and a decrease of 26% and 48% in alpine meadows and alpine steppes, respectively. These trends were primarily driven by temperature-induced growth preferences rather than shifts in plant species composition. Our findings indicate that in wetter ecosystems, climate warming promotes aboveground plant growth, while in drier ecosystems, such as alpine meadows and alpine steppes, plants allocate more biomass belowground. Furthermore, we observed a threefold strengthening of the warming effect on η over the past 27 y. Soil moisture was found to modulate the sensitivity of η to soil temperature in alpine meadows and alpine steppes, but not in alpine wetlands. Our results contribute to a better understanding of the processes driving the response of biomass distribution to climate warming, which is crucial for predicting the future carbon trajectory of permafrost ecosystems and climate feedback.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep.

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Effectiveness of a comprehensive package based on electronic medication monitors at improving treatment outcomes among tuberculosis patients in Tibet: a multicentre randomised controlled trial.

BACKGROUND: WHO recommends that electronic medication monitors, a form of digital adherence technology, be used as a complement to directly observed treatment (DOT) for tuberculosis, as DOT is inconvenient and costly. However, existing evidence about the effectiveness of these monitors is inconclusive. Therefore, we evaluated the effectiveness of a comprehensive package based on electronic medication monitors among patients with tuberculosis in Tibet Autonomous Region (hereafter Tibet), China. METHODS: This multicentre, randomised controlled trial recruited patients from six counties in Shigatse, Tibet. Eligible participants had drug-susceptible tuberculosis and were aged 15 years or older when starting standard tuberculosis treatment. Tuberculosis doctors recruited patients from the public tuberculosis dispensary in each county and the study statistician randomly assigned them to the intervention or control group based on the predetermined randomised allocation sequence. Intervention patients received an electronic medication monitor box. The box included audio medication-adherence reminders and recorded box-opening data, which were transmitted to a cloud-based server and were accessible to health-care providers to allow remote adherence monitoring. A linked smartphone app enabled text, audio, and video communication between patients and health-care providers. Patients were also provided with a free data plan. Patients selected a treatment supporter (often a family member) who was trained to support patients with using the electronic medication monitor and app. Patients in the control group received usual care plus a deactivated electronic medication monitor, which only recorded and transmitted box-opening data that was not made available to health-care providers. The control group also had no access to the app or trained treatment supporters. The primary outcome was a binary indicator of poor monthly adherence, defined as missing 20% or more of planned doses in the treatment month, measured using electronic medication monitor opening data, and verified by counting used medication blister packages during consultations. We recorded other secondary treatment outcomes based on national tuberculosis reporting data. We analysed the primary outcome based on the intention-to-treat population. This trial is registered at ISRCTN, 52132803. FINDINGS: Between Nov 17, 2018, and April 5, 2021, 278 patients were enrolled into the study. 143 patients were randomly assigned to the intervention group and 135 patients to the control group. Follow-up ended when the final patient completed treatment on Oct 4, 2021. In the intervention group, 87 (10%) of the 854 treatment months showed poor adherence compared with 290 (37%) of the 795 months in the control group. The corresponding adjusted risk difference for the intervention versus control was -29·2 percentage points (95% CI -35·3 to -22·2; p<0·0001). Five of the six secondary treatment outcomes also showed clear improvements, including treatment success, which was found for 133 (94%) of the 142 individuals in the intervention arm and 98 (73%) of the 134 individuals in the control arm, with an adjusted risk difference of 21 percentage points (95% CI 12·4-29·4); p<0·0001. INTERPRETATION: The interventions were effective at improving tuberculosis treatment adherence and outcomes, and the trial suggests that a comprehensive package involving electronic medication monitors might positively affect tuberculosis programmes in high-burden and low-resource settings. FUNDING: TB REACH.

Phylogenetic evidence for Sino-Tibetan origin in northern China in the Late Neolithic.

The study of language origin and divergence is important for understanding the history of human populations and their cultures. The Sino-Tibetan language family is the second largest in the world after Indo-European, and there is a long-running debate about its phylogeny and the time depth of its original divergence1. Here we perform a Bayesian phylogenetic analysis to examine two competing hypotheses of the origin of the Sino-Tibetan language family: the 'northern-origin hypothesis' and the 'southwestern-origin hypothesis'. The northern-origin hypothesis states that the initial expansion of Sino-Tibetan languages occurred approximately 4,000-6,000 years before present (BP; taken as AD 1950) in the Yellow River basin of northern China2-4, and that this expansion is associated with the development of the Yangshao and/or Majiayao Neolithic cultures. The southwestern-origin hypothesis states that an early expansion of Sino-Tibetan languages occurred before 9,000 years BP from a region in southwest Sichuan province in China5 or in northeast India6, where a high diversity of Tibeto-Burman languages exists today. Consistent with the northern-origin hypothesis, our Bayesian phylogenetic analysis of 109 languages with 949 lexical root-meanings produced an estimated time depth for the divergence of Sino-Tibetan languages of approximately 4,200-7,800 years BP, with an average value of approximately 5,900 years BP. In addition, the phylogeny supported a dichotomy between Sinitic and Tibeto-Burman languages. Our results are compatible with the archaeological records, and with the farming and language dispersal hypothesis7 of agricultural expansion in China. Our findings provide a linguistic foothold for further interdisciplinary studies of prehistoric human activity in East Asia.

A late Middle Pleistocene Denisovan mandible from the Tibetan Plateau.

Denisovans are members of a hominin group who are currently only known directly from fragmentary fossils, the genomes of which have been studied from a single site, Denisova Cave1-3 in Siberia. They are also known indirectly from their genetic legacy through gene flow into several low-altitude East Asian populations4,5 and high-altitude modern Tibetans6. The lack of morphologically informative Denisovan fossils hinders our ability to connect geographically and temporally dispersed fossil hominins from Asia and to understand in a coherent manner their relation to recent Asian populations. This includes understanding the genetic adaptation of humans to the high-altitude Tibetan Plateau7,8, which was inherited from the Denisovans. Here we report a Denisovan mandible, identified by ancient protein analysis9,10, found on the Tibetan Plateau in Baishiya Karst Cave, Xiahe, Gansu, China. We determine the mandible to be at least 160 thousand years old through U-series dating of an adhering carbonate matrix. The Xiahe specimen provides direct evidence of the Denisovans outside the Altai Mountains and its analysis unique insights into Denisovan mandibular and dental morphology. Our results indicate that archaic hominins occupied the Tibetan Plateau in the Middle Pleistocene epoch and successfully adapted to high-altitude hypoxic environments long before the regional arrival of modern Homo sapiens.

Seismic constraints and geodynamic implications of differential lithosphere-asthenosphere flow revealed in East Asia.

During the last 50 Ma, the East Asian continent has been a zone of massive continental collision and lithospheric deformation. While the consequences of this for Asian surface and lithospheric deformation have been intensively studied over the past 4 decades, the relationships between lithospheric deformation and underlying asthenospheric flow have been more difficult to constrain. Here we present a high resolution 3-D azimuthal anisotropy model for the northeastern Tibetan Plateau and its eastward continuation based on surface-wave tomography and shear-wave splitting measurements. This model shows that eastward lateral flow of asthenosphere beneath the northeastern Tibetan Plateau is being blocked by thick Ordos and Sichuan cratonic keels. The damming effect of these keels induces flow to first rotate around the Ordos keel and then transition into strong east-west flow beneath the thinner lithosphere that forms the lithospheric suture between the two cratonic keels. We further find that asthenosphere flow directions can differ from those of overlying lithosphere, with the asthenosphere neither being passively dragged by overlying lithosphere, nor being able to drag the overlying plate to mimic its subsurface flow. Finally, the region of eastward-channeled asthenospheric flow from Tibet underlies a belt of stronger intracontinental deformation in eastern China.

Warming-induced tipping points of Arctic and alpine shrub recruitment.

Shrub recruitment, a key component of vegetation dynamics beyond forests, is a highly sensitive indicator of climate and environmental change. Warming-induced tipping points in Arctic and alpine treeless ecosystems are, however, little understood. Here, we compare two long-term recruitment datasets of 2,770 shrubs from coastal East Greenland and from the Tibetan Plateau against atmospheric circulation patterns between 1871 and 2010 Common Era. Increasing rates of shrub recruitment since 1871 reached critical tipping points in the 1930s and 1960s on the Tibetan Plateau and in East Greenland, respectively. A recent decline in shrub recruitment in both datasets was likely related to warmer and drier climates, with a stronger May to July El Niño Southern Oscillation over the Tibetan Plateau and a stronger June to July Atlantic Multidecadal Oscillation over Greenland. Exceeding the thermal optimum of shrub recruitment, the recent warming trend may cause soil moisture deficit. Our findings suggest that changes in atmospheric circulation explain regional climate dynamics and associated response patterns in Arctic and alpine shrub communities, knowledge that should be considered to protect vulnerable high-elevation and high-latitude ecosystems from the cascading effects of anthropogenic warming.

Use of δ18Oatm in dating a Tibetan ice core record of Holocene/Late Glacial climate.

Ice cores from the northwestern Tibetan Plateau (NWTP) contain long records of regional climate variability, but refrozen meltwater and dust in these cores has hampered development of robust timescales. Here, we introduce an approach to dating the ice via the isotopic composition of atmospheric O2 in air bubbles (δ18Oatm), along with annual layer counting and radiocarbon dating. We provide a robust chronology for water isotope records (δ18Oice and d-excess) from three ice cores from the Guliya ice cap in the NWTP. The measurement of δ18Oatm, although common in polar ice core timescales, has rarely been used on ice cores from low-latitude, high-altitude glaciers due to (1) low air pressure, (2) the common presence of refrozen melt that adds dissolved gases and reduces the amount of air available for analysis, and (3) the respiratory consumption of molecular oxygen (O2) by micro-organisms in the ice, which fractionates the δ18O of O2 from the atmospheric value. Here, we make corrections for melt and respiration to address these complications. The resulting records of water isotopes from the Guliya ice cores reveal climatic variations over the last 15,000 y, the timings of which correspond to those observed in independently dated lake and speleothem records and confirm that the Guliya ice cap existed before the Holocene. The millennial-scale drivers of δ18Oice are complex and temporally variable; however, Guliya δ18Oice values since the mid-20th century are the highest since the beginning of the Holocene and have increased with regional air temperature.

A Tibetan ice core covering the past 1,300 years radiometrically dated with 39Ar.

Ice cores from alpine glaciers are unique archives of past global and regional climate conditions. However, recovering climate records from these ice cores is often hindered by the lack of a reliable chronology, especially in the age range of 100 to 500 anni (a) for which radiometric dating has not been available so far. We report on radiometric 39Ar dating of an ice core from the Tibetan Plateau and the construction of a chronology covering the past 1,300 a using the obtained 39Ar ages. This is made possible by advances in the analysis of 39Ar using the laser-based detection method atom trap trace analysis, resulting in a twofold increase in the upper age limit of 39Ar dating. By measuring the anthropogenic 85Kr along with 39Ar we quantify and correct modern air contamination, thus removing a major systematic uncertainty of 39Ar dating. Moreover, the 85Kr data for the top part of the ice core provide information on firn processes, including the age difference between the ice and its enclosed gas. This first application of 39Ar and 85Kr to an ice core facilitates further ice cores from nonpolar glaciers to be used for recovering climate records of the Common Era, a period including pronounced anomalies such as the Little Ice Age and the Medieval Warm Period.

Genomic insights into zokors' phylogeny and speciation in China.

The phylogeny and speciation of subterranean zokors in China are unclear, as previous studies on morphology and limited molecular markers have generated conflicting results. This study unraveled the complex evolutionary history of eight zokor species in China based on de novo assembly at chromosome level and whole-genome sequencing of 23 populations. We found extensive phylogenetic discordances between nuclear and mitochondrial phylogenies, and different coalescent phylogenies, which could be explained by introgression and incomplete lineage sorting (ILS). The recent Qinghai-Tibet Plateau uplift (∼3.60 million y ago; Mya) drove Eospalax to speciate into clade A and clade B (∼3.22 Mya), and discordant phylogenies in this node were mainly attributed to introgression rather than ILS. Clade A rapidly diverged into three lineages due to geographical isolation and glaciation, while glaciation and C4 plant expansion contributed to the speciation of clade B. ILS contributed to the discordances of two rapidly radiated nodes rather than introgression. The effective population sizes (Ne's) of all the species of Eospalax were affected by three glaciations. Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of all the species pairs. Positively selected genes putatively related to specific inhabitation adaptations were identified, such as heart development, neurogenesis, DNA repair, and immune response. Climate, geological tectonism, and C4 vegetation shaped the adaptation and speciation of zokors in China.

Genetic adaptation of skin pigmentation in highland Tibetans.

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes.

Tibetan cashmere goats are not only served as a valuable model for studying adaptation to hypoxia and high-altitude conditions but also playing a pivotal role in bolstering local economies through the provision of premium quality cashmere yarn. In this study, we performed an integration and network analysis of metabolomic, transcriptomic and proteomic to elucidate the role of differentially expressed genes, important metabolites, and relevant cellular and metabolic pathways between the fine (average 12.04 ± 0.03 µm of mean fiber diameter) and coarse cashmere (average 14.88 ± 0.05 µm of mean fber diameter) producing by Tibetan cashmere goats. We identified a distinction of 56 and 71 differential metabolites (DMs) between the F and C cashmere groups under positive and negative ion modes, respectively. The KEGG pathway enrichment analysis of these DMs highlighted numerous pathways predominantly involved in amino acid and protein metabolism, as indicated by the finding that the most impactful pathway was the mammalian target of rapamycin (mTOR) signalling pathway. In the F group, we identified a distinctive metabolic profile where amino acid metabolites including serine, histidine, asparagine, glutamic acid, arginine, valine, aspartic acid, tyrosine, and methionine were upregulated, while lysine, isoleucine, glutamine, tryptophan, and threonine were downregulated. The regulatory network and gene co-expression network revealed crucial genes, metabolites, and metabolic pathways. The integrative omics analysis revealed a high enrichment of several pathways, notably encompassing protein digestion and absorption, sphingolipid signalling, and the synaptic vesicle cycle. Within the sphere of our integrative analysis, DNMT3B was identified as a paramount gene, intricately associated with significant proteins such as HMCN1, CPB2, GNG12, and LRP1. Our present study delineated the molecular underpinnings governing the variations in cashmere characteristics by conducting comprehensive analyses across metabolomic, transcriptomic, and proteomic dimensions. This research provided newly insights into the mechanisms regulating cashmere traits and facilitated the advancement of selective breeding programs aimed at cultivating high-quality superfine Tibetan cashmere goats.

Polymorphic transposable elements provide new insights on high-altitude adaptation in the Tibetan Plateau.

Several studies demonstrated that populations living in the Tibetan plateau are genetically and physiologically adapted to high-altitude conditions, showing genomic signatures ascribable to the action of natural selection. However, so far most of them relied solely on inferences drawn from the analysis of coding variants and point mutations. To fill this gap, we focused on the possible role of polymorphic transposable elements in influencing the adaptation of Tibetan and Sherpa highlanders. To do so, we compared high-altitude and middle/low-lander individuals of East Asian ancestry by performing in silico analyses and differentiation tests on 118 modern and ancient samples. We detected several transposable elements associated with high altitude, which map genes involved in cardiovascular, hematological, chem-dependent and respiratory conditions, suggesting that metabolic and signaling pathways taking part in these functions are disproportionately impacted by the effect of environmental stressors in high-altitude individuals. To our knowledge, our study is the first hinting to a possible role of transposable elements in the adaptation of Tibetan and Sherpa highlanders.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing.

BACKGROUND: Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method. RESULTS: The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism. CONCLUSION: Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

High genetic diversity of the himalayan marmot relative to plague outbreaks in the Qinghai-Tibet Plateau, China.

Plague, as an ancient zoonotic disease caused by Yersinia pestis, has brought great disasters. The natural plague focus of Marmota himalayana in the Qinghai-Tibet Plateau is the largest, which has been constantly active and the leading source of human plague in China for decades. Understanding the population genetics of M. himalayana and relating that information to the biogeographic distribution of Yersinia pestis and plague outbreaks are greatly beneficial for the knowledge of plague spillover and arecrucial for pandemic prevention. In the present research, we assessed the population genetics of M. himalayana. We carried out a comparative study of plague outbreaks and the population genetics of M. himalayana on the Qinghai-Tibet Plateau. We found that M. himalayana populations are divided into two main clusters located in the south and north of the Qinghai-Tibet Plateau. Fourteen DFR genomovars of Y. pestis were found and exhibited a significant region-specific distribution. Additionally, the increased genetic diversity of plague hosts is positively associated with human plague outbreaks. This insight gained can improve our understanding of biodiversity for pathogen spillover and provide municipally directed targets for One Health surveillance development, which will be an informative next step toward increased monitoring of M. himalayana dynamics.

Metabolomics and transcriptomics of embryonic livers reveal hypoxia adaptation of Tibetan chickens.

BACKGROUND: Exploring the hypoxia adaptation mechanism of Tibetan chicken is of great significance for revealing the survival law of Tibetan chicken and plateau animal husbandry production. To investigate the hypoxia adaptation of Tibetan chickens (TBCs), an integrative metabolomic-transcriptomic analysis of the liver on day 18 of embryonic development was performed. Dwarf laying chickens (DLCs), a lowland breed, were used as a control. RESULTS: A total of 1,908 metabolites were identified in both TBCs and DLCs. Energy metabolism and amino acid metabolism related differentially regulated metabolites (DRMs) were significantly enriched under hypoxia. Important metabolic pathways including the TCA cycle and arginine and proline metabolism were screened; PCK1, SUCLA2, and CPS1 were found to be altered under hypoxic conditions. In addition, integrated analysis suggested potential differences in mitochondrial function, which may play a crucial role in the study of chicken oxygen adaptation. CONCLUSIONS: These results suggest that hypoxia changed the gene expression and metabolic patterns of embryonic liver of TBCs compared to DLCs. Our study provides a basis for uncovering the molecular regulation mechanisms of hypoxia adaptation in TBCs with the potential application of hypoxia adaptation research for other animals living on the Qinghai-Tibet plateau, and may even contribute to the study of diseases caused by hypoxia.

Comparative transcriptomic analysis delineates adaptation strategies of Rana kukunoris toward cold stress on the Qinghai-Tibet Plateau.

BACKGROUND: Cold hardiness is fundamental for amphibians to survive during the extremely cold winter on the Qinghai-Tibet plateau. Exploring the gene regulation mechanism of freezing-tolerant Rana kukunoris could help us to understand how the frogs survive in winter. RESULTS: Transcriptome of liver and muscle of R. kukunoris collected in hibernation and spring were assisted by single molecule real-time (SMRT) sequencing technology. A total of 10,062 unigenes of R. kukunoris were obtained, and 9,924 coding sequences (CDS) were successfully annotated. Our examination of the mRNA response to whole body freezing and recover in the frogs revealed key genes concerning underlying antifreeze proteins and cryoprotectants (glucose and urea). Functional pathway analyses revealed differential regulated pathways of ribosome, energy supply, and protein metabolism which displayed a freeze-induced response and damage recover. Genes related to energy supply in the muscle of winter frogs were up-regulated compared with the muscle of spring frogs. The liver of hibernating frogs maintained modest levels of protein synthesis in the winter. In contrast, the liver underwent intensive high levels of protein synthesis and lipid catabolism to produce substantial quantity of fresh proteins and energy in spring. Differences between hibernation and spring were smaller than that between tissues, yet the physiological traits of hibernation were nevertheless passed down to active state in spring. CONCLUSIONS: Based on our comparative transcriptomic analyses, we revealed the likely adaptive mechanisms of R. kukunoris. Ultimately, our study expands genetic resources for the freezing-tolerant frogs.