Multi-omics reveal the gut microbiota-mediated severe foraging environment adaption of small wild ruminants in the Three-River-Source National Park, China.

The Tibetan antelope (Pantholops hodgsonii), blue sheep (Pseudois nayaur), and Tibetan sheep (Ovis aries) are the dominant small ruminants in the Three-River-Source National Park (TRSNP). However, knowledge about the association between gut microbiota and host adaptability remains poorly understood. Herein, multi-omics sequencing approaches were employed to investigate the gut microbiota-mediated forage adaption in these ruminants. The results revealed that although wild ruminants (WR) of P. hodgsoni and P. nayaur were faced with severe foraging environments with significantly low vegetation coverage and nutrition, the apparent forage digestibility of dry matter, crude protein, and acid detergent fiber was significantly higher than that of O. aries. The 16s rRNA sequencing showed that the gut microbiota in WR underwent convergent evolution, and alpha diversity in these two groups was significantly higher than that in O. aries. Moreover, indicator species, including Bacteroidetes and Firmicutes, exhibited positive relationships with apparent forage digestibility, and their relative abundances were enriched in the gut of WR. Enterotype analysis further revealed that enterotype 1 belonged to WR, and the abundance of fatty acid synthesis metabolic pathway-related enzyme genes was significantly higher than enterotype 2, represented by O. aries. Besides, the metagenomic analysis identified 14 pathogenic bacterial species, among which 10 potentially pathogenic bacteria were significantly enriched in the gut microbiota of O. aries. Furthermore, the cellulolytic strains and genes encoding cellulase and hemicellulase were significantly enriched in WR. In conclusion, our results provide new evidence of gut microbiota to facilitate wildlife adaption in severe foraging environments of the TRSNP, China.

Influences of elevational gradient on flower size and number of Gentiana lawrencei var. farreri.

Plants can adapt to environmental changes by adjusting their functional traits and biomass allocation. The size and number of flowers are functional traits related to plant reproduction. Life history theory predicts that there is a trade-off between flower size and number, and the trade-off can potentially explain the adaptability of plants. Elevation gradients in mountains provide a unique opportunity to test how plants will respond to climate change. In this study, we tried to better explain the adaptability of the alpine plant Gentiana lawrencei var. farreri in response to climate change. We measured the flower size and number, individual size, and reproductive allocation of G. lawrencei var. farreri during the flowering period along an elevation gradient from 3200 to 4000 m, and explored their relationships using linear mixed-effect models and the structural equation model. We found that with elevation increasing, individual size and flower number decreased and flower size increased, while reproductive allocation remained unchanged. Individual size positively affected flower number, but was not related to flower size; reproductive allocation positively affected flower size, but was not related to flower number; there is a clear trade-off between flower size and number. We also found that elevation decreased flower number indirectly via directly reducing individual size. In sum, this study suggests that G. lawrencei var. farreri can adapt to alpine environments by the synergies or trade-offs among individual size, reproductive allocation, flower size, and flower number. This study increases our understanding of the adaptation mechanisms of alpine plants to climate change in alpine environments.

Secondary formation of oxygenated and nitrated polycyclic aromatic compounds under stagnant weather conditions: Drivers and seasonal variation.

Severe air pollution tends to occur under stagnant weather conditions. This study focused on the occurrence and formation of PM2.5-bound polycyclic aromatic compounds (PACs) under stagnant weather conditions, in consideration of their adverse human health effect and ecological toxicity. The concentrations of PACs were higher under stagnant weather conditions than in other situations with averaged values of 46.0 ng/m3 versus 12.3-39.9 ng/m3 for total PACs. Secondary formation contributed to over half of the oxygenated and nitrated polycyclic aromatic compounds (OPAHs and NPAHs). Further analyses revealed different formation mechanisms for secondary OPAHs and NPAHs. Secondary production of OPAHs was sensitive to the variations of both temperature (T) and O3 concentration at T < 22 °C but sustained at a high level despite the fluctuation of temperature and O3 concentration at T > 22 °C. Elevated NO2 concentrations favored the formation of inorganic nitrogen-containing products over NPAHs under lower temperature and higher humidity. Stagnant weather events, accompanied by raised PAC levels occurred in all seasons, but their effects on secondary processes differed among seasons. The elevated temperature, lowered humidity, and increased NO2 level facilitated the secondary formation of OPAHs and/or NPAHs during the stagnant weather events in spring and summer. While under the temperature and humidity conditions in autumn and winter, increased NO2 levels during stagnant weather events promoted the production of secondary inorganic nitrogen-containing compounds over organic products. This study raised concern about the toxic organic pollutants in the atmosphere under stagnant weather conditions and revealed different formation mechanisms between secondary oxygenated and nitrated pollutants as well as among different seasons.

Trophic niche and adaptation in highland lizards: sex has greater influences than species matching.

The plateau environments are typically arid, cool, and high altitude, posing formidable challenges to wildlife survival due to resource scarcity and harsh conditions. Unraveling ecological adaptability in severe conditions requires a deeper understanding of the niche characteristics of plateau species. Trophic niche, which is a comprehensive indicator describing the energy acquisition strategy of animals, remains relatively understudied in plateau species. Here, by combining stable isotopes and morphological data, we quantified the trophic niches of two allopatric lizard species (Phrynocephalus vlangalii and P. erythrurus) that live in the hinterland of the Qinghai-Tibetan Plateau, and explored how their trophic niches correlate with morphological and environmental factors. While both trophic niche and morphological traits were similar between species, noteworthy distinctions were observed between male and female Phrynocephalus lizards. The morphological traits associated with predation (i.e. limb length and head size) and reproduction (i.e. abdomen length), annual mean temperature, and sex played influential roles in shifting trophic niches. These results imply that sexual dimorphism may facilitate inter-sex divergence in resource utilization, leading to trophic niche variations in the highland lizards. Furthermore, extreme environmental stress can constrain interspecific divergence in morphological and trophic traits. Our findings illustrate the dynamic variations of trophic niches in highland lizards, contributing to a more comprehensive understanding of the adaptation strategies employed by lizard species in plateau environments.

An early warning signal for grassland degradation on the Qinghai-Tibetan Plateau.

Intense grazing may lead to grassland degradation on the Qinghai-Tibetan Plateau, but it is difficult to predict where this will occur and to quantify it. Based on a process-based ecosystem model, we define a productivity-based stocking rate threshold that induces extreme grassland degradation to assess whether and where the current grazing activity in the region is sustainable. We find that the current stocking rate is below the threshold in ~80% of grassland areas, but in 55% of these grasslands the stocking rate exceeds half the threshold. According to our model projections, positive effects of climate change including elevated CO2 can partly offset negative effects of grazing across nearly 70% of grasslands on the Plateau, but only in areas below the stocking rate threshold. Our analysis suggests that stocking rate that does not exceed 60% (within 50% to 70%) of the threshold may balance human demands with grassland protection in the face of climate change.

Vegetation structural shift tells environmental changes on the Tibetan Plateau over 40 years.

Structural information of grassland changes on the Tibetan Plateau is essential for understanding alterations in critical ecosystem functioning and their underlying drivers that may reflect environmental changes. However, such information at the regional scale is still lacking due to methodological limitations. Beyond remote sensing indicators only recognizing vegetation productivity, we utilized multivariate data fusion and deep learning to characterize formation-based plant community structure in alpine grasslands at the regional scale of the Tibetan Plateau for the first time and compared it with the earlier version of Vegetation Map of China for historical changes. Over the past 40 years, we revealed that (1) the proportion of alpine meadows in alpine grasslands increased from 50% to 69%, well-reflecting the warming and wetting trend; (2) dominances of Kobresia pygmaea and Stipa purpurea formations in alpine meadows and steppes were strengthened to 76% and 92%, respectively; (3) the climate factor mainly drove the distribution of Stipa purpurea formation, but not the recent distribution of Kobresia pygmaea formation that was likely shaped by human activities. Therefore, the underlying mechanisms of grassland changes over the past 40 years were considered to be formation dependent. Overall, the first exploration for structural information of plant community changes in this study not only provides a new perspective to understand drivers of grassland changes and their spatial heterogeneity at the regional scale of the Tibetan Plateau, but also innovates large-scale vegetation study paradigm.

A functional metagenomics study of soil carbon and nitrogen degradation networks and limiting factors on the Tibetan plateau.

Introduction: The Three-River Source Nature Reserve is located in the core area of the Qinghai-Tibetan Plateau, with the alpine swamp, meadow and steppe as the main ecosystem types. However, the microbial communities in these alpine ecosystems, and their carbon and nitrogen degrading metabolic networks and limiting factors remain unclear. Methods: We sequenced the diversity of bacteria and fungi in alpine swamps, meadows, steppes, and their degraded and artificially restored ecosystems and analyzed soil environmental conditions. Results: The results indicated that moisture content had a greater influence on soil microbial community structure compared to degradation and restoration. Proteobacteria dominated in high moisture alpine swamps and alpine meadows, while Actinobacteria dominated in low moisture alpine steppes and artificial grasslands. A metabolic network analysis of carbon and nitrogen degradation and transformation using metagenomic sequencing revealed that plateau microorganisms lacked comprehensive and efficient enzyme systems to degrade organic carbon, nitrogen, and other biological macromolecules, so that the short-term degradation of alpine vegetation had no effect on the basic composition of soil microbial community. Correlation analysis found that nitrogen fixation was strong in meadows with high moisture content, and their key nitrogen-fixing enzymes were significantly related to Sphingomonas. Denitrification metabolism was enhanced in water-deficient habitats, and the key enzyme, nitrous oxide reductase, was significantly related to Phycicoccus and accelerated the loss of nitrogen. Furthermore, Bacillus contained a large number of amylases (GH13 and GH15) and proteases (S8, S11, S26, and M24) which may promote the efficient degradation of organic carbon and nitrogen in artificially restored grasslands. Discussion: This study illustrated the irrecoverability of meadow degradation and offered fundamental information for altering microbial communities to restore alpine ecosystems.

Herbivore assemblages affect soil microbial communities by altering root biomass and available nutrients in an alpine meadow.

Three different herbivore grazing assemblages, namely, yak grazing (YG), Tibetan sheep grazing (SG) and yak and Tibetan sheep co-grazing (MG), are practiced in alpine meadows on the Qinghai-Tibetan Plateau (QTP), but the effects of the different herbivore assemblages on soil microbes are relatively unknown. The microbial community plays an important role in the functional stability of alpine grassland ecosystems. Therefore, it is important to understand how the microbial community structure of grassland ecosystems changes under different herbivore grazing assemblages to ensure their sustainable development. To fill this gap, a field study was carried out to investigate the effects of YG, SG, and MG on plant communities, soil physico-chemical properties and microbial communities under moderate grazing intensity in alpine meadows. Grazing increased the ß-diversity of the bacteria community and decreased the ß-diversity of the fungal community. The herbivore assemblage affected the microbial community diversity, but not the plant community diversity. Total phosphorus, soil bulk density, root biomass, and plant α-diversity were correlated with both the bacterial and fungal community composition, available phosphorus and soil moisture were correlated only with the bacterial community composition, while available potassium and above-ground net primary production (ANPP) were correlated only with the fungal community composition. Soil available nitrogen, soil available phosphorus and soil bulk density were highest in SG, while ANPP was highest in MG. It was concluded that MG can improve ANPP and stabilize the soil microbial community, suggesting that MG is an effective method for sustainable use and conservation of alpine meadows on the QTP.

Long read genome assemblies complemented by single cell RNA-sequencing reveal genetic and cellular mechanisms underlying the adaptive evolution of yak.

Wild yak (Bos mutus) and domestic yak (Bos grunniens) are adapted to high altitude environment and have ecological, economic, and cultural significances on the Qinghai-Tibetan Plateau (QTP). Currently, the genetic and cellular bases underlying adaptations of yak to extreme conditions remains elusive. In the present study, we assembled two chromosome-level genomes, one each for wild yak and domestic yak, and screened structural variants (SVs) through the long-read data of yak and taurine cattle. The results revealed that 6733 genes contained high-FST SVs. 127 genes carrying special type of SVs were differentially expressed in lungs of the taurine cattle and yak. We then constructed the first single-cell gene expression atlas of yak and taurine cattle lung tissues and identified a yak-specific endothelial cell subtype. By integrating SVs and single-cell transcriptome data, we revealed that the endothelial cells expressed the highest proportion of marker genes carrying high-FST SVs in taurine cattle lungs. Furthermore, we identified pathways which were related to the medial thickness and formation of elastic fibers in yak lungs. These findings provide new insights into the high-altitude adaptation of yak and have important implications for understanding the physiological and pathological responses of large mammals and humans to hypoxia.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing.

Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.

The Gut Microbiota Determines the High-Altitude Adaptability of Tibetan Wild Asses (Equus kiang) in Qinghai-Tibet Plateau.

It was acknowledged long ago that microorganisms have played critical roles in animal evolution. Tibetan wild asses (TWA, Equus kiang) are the only wild perissodactyls on the Qinghai-Tibet Plateau (QTP) and the first national protected animals; however, knowledge about the relationships between their gut microbiota and the host's adaptability remains poorly understood. Herein, 16S rRNA and meta-genomic sequencing approaches were employed to investigate the gut microbiota-host associations in TWA and were compared against those of the co-resident livestock of yak (Bos grunnies) and Tibetan sheep (Ovis aries). Results revealed that the gut microbiota of yak and Tibetan sheep underwent convergent evolution. By contrast, the intestinal microflora of TWA diverged in a direction enabling the host to subsist on sparse and low-quality forage. Meanwhile, high microbial diversity (Shannon and Chao1 indices), cellulolytic activity, and abundant indicator species such as Spirochaetes, Bacteroidetes, Prevotella_1, and Treponema_2 supported forage digestion and short-chain fatty acid production in the gut of TWA. Meanwhile, the enterotype identification analysis showed that TWA shifted their enterotype in response to low-quality forage for a better utilization of forage nitrogen and short-chain fatty acid production. Metagenomic analysis revealed that plant biomass degrading microbial consortia, genes, and enzymes like the cellulolytic strains (Prevotella ruminicola, Ruminococcus flavefaciens, Ruminococcus albus, Butyrivibrio fibrisolvens, and Ruminobacter amylophilus), as well as carbohydrate metabolism genes (GH43, GH3, GH31, GH5, and GH10) and enzymes (ß-glucosidase, xylanase, and ß-xylosidase, etc.) had a significantly higher enrichment in TWA. Our results indicate that gut microbiota can improve the adaptability of TWA through plant biomass degradation and energy maintenance by the functions of gut microbiota in the face of nutritional deficiencies and also provide a strong rationale for understanding the roles of gut microbiota in the adaptation of QTP wildlife when facing harsh feeding environments.

Biodiversity and Ecosystem Function under Simulated Gradient Warming and Grazing.

Biodiversity and ecosystem functions and their relationship with environmental response constitute a major topic of ecological research. However, the changes in and impact mechanisms of multi-dimensional biodiversity and ecosystem functions in continuously changing environmental gradients and anthropogenic activities remain poorly understood. Here, we analyze the effects of multi-gradient warming and grazing on relationships between the biodiversity of plant and soil microbial with productivity/community stability through a field experiment simulating multi-gradient warming and grazing in alpine grasslands on the Tibetan Plateau. We show the following results: (i) Plant biodiversity, soil microbial diversity and community productivity in alpine grasslands show fluctuating trends with temperature gradients, and a temperature increase below approximately 1 °C is beneficial to alpine grasslands; moderate grazing only increases the fungal diversity of the soil surface layer. (ii) The warming shifted plant biomass underground in alpine grasslands to obtain more water in response to the decrease in soil moisture caused by the temperature rise. Community stability was not affected by warming or grazing. (iii) Community stability was not significantly correlated with productivity, and environmental factors, rather than biodiversity, influenced community stability and productivity.

Cropping practices manipulate soil bacterial structure and functions on the Qinghai-Tibet Plateau.

There is an increasing awareness of the adverse environmental effects of the intensive practices used in modern crop farming, such as those that cause greenhouse gas emissions and nutrient leaching. Harnessing beneficial microbes by changing planting practices presents a promising strategy for optimizing plant growth and agricultural sustainability. However, the characteristics of soil microorganisms under different planting patterns remain uncertain. We conducted a study of soil bacterial structure and function under monoculture vs. polyculture planting regimes using 16S rRNA gene sequencing on the Qinghai-Tibet Plateau. We observed substantial variations in bacterial richness, diversity, and relative abundances of taxa between gramineous and leguminous monocultures, as well as between gramineae-legume polycultures. The number of operational taxonomic units and alpha and beta diversity were markedly higher in the leguminous monocultures than in the gramineous monocultures; conversely, network analysis revealed that the interactions among the bacterial genera in the gramineous monocultures were more complex than those in the other two planting regimes. Moreover, nitrogen fixation, soil detoxification, and productivity were increased under the gramineous monocultures; more importantly, low soil-borne diseases (e.g., animals parasitic or symbiont) also facilitated strongly suppressive effects toward soil-borne pathogens. Nevertheless, the gramineae-legume polycultures were prone to nitrate seepage contamination, and leguminous monocultures exhibited strong denitrification effects. These results revealed that the gramineous monoculture is a more promising cropping pattern on the Qinghai-Tibetan Plateau. Understanding the bacterial distribution patterns and interactions of crop-sensitive microbes presents a basis for developing microbial management strategies for smart farming.

Warming has a minor effect on surface soil organic carbon in alpine meadow ecosystems on the Qinghai-Tibetan Plateau.

The alpine meadow ecosystem on the Qinghai-Tibetan Plateau (QTP) is very sensitive to warming and plays a key role in regulating global carbon (C) cycling. However, how warming affects the soil organic carbon (SOC) pool and related C inputs and outputs in alpine meadow ecosystems on the QTP remains unclear. Here, we combined two field experiments and a meta-analysis on field experiments to synthesize the responses of the SOC pool and related C cycling processes to warming in alpine meadow ecosystems on the QTP. We found that the SOC content of surface soil (0-10 cm) showed a minor response to warming, but plant respiration was accelerated by warming. In addition, the warming effect on SOC was not correlated with experimental and environmental variables, such as the method, magnitude and duration of warming, initial SOC content, mean annual temperature, and mean annual precipitation. We conclude that the surface SOC content is resistant to climate warming in alpine meadow ecosystems on the QTP.

Evenness is important in assessing progress towards sustainable development goals.

Sustainable development goals (SDGs) emphasize a holistic achievement instead of cherry-picking a few. However, no assessment has quantitatively considered the evenness among all 17 goals. Here, we propose a systematic method, which first integrates both the evenness and the overall status of all goals, to distinguish the ideal development pathways from the uneven ones and then revisit the development trajectory in China from 2000 to 2015. Our results suggest that, despite the remarkable progress, a bottleneck has occurred in China since 2013 due to the stagnant developments in some SDGs. However, many far-reaching policies in China have been targeting these deficiencies since then, providing a perspective on how a country approaches sustainable development by promoting evenness among all SDGs. Our results also indicate that regions with the slowest progress are the developed provinces, owing to the persistent uneven status of all goals. Our study demonstrates the importance of adopting evenness in assessing and guiding sustainable development.

Ambient climate determines the directional trend of community stability under warming and grazing.

Changes in ecological processes over time in ambient treatments are often larger than the responses to manipulative treatments in climate change experiments. However, the impacts of human-driven environmental changes on the stability of natural grasslands have been typically assessed by comparing differences between manipulative plots and reference plots. Little is known about whether or how ambient climate regulates the effects of manipulative treatments and their underlying mechanisms. We collected two datasets, one a 36-year long-term observational dataset from 1983 to 2018, and the other a 10-year manipulative asymmetric warming and grazing experiment using infrared heaters with moderate grazing from 2006 to 2015 in an alpine meadow on the Tibetan Plateau. The 36-year observational dataset shows that there was a nonlinear response of community stability to ambient temperature with a positive relationship between them due to an increase in ambient temperature in the first 25 years and then a decrease in ambient temperature thereafter. Warming and grazing decreased community stability with experiment duration through an increase in legume cover and a decrease in species asynchrony, which was due to the decreasing background temperature through time during the 10-year experiment period. Moreover, the temperature sensitivity of community stability was higher under the ambient treatment than under the manipulative treatments. Therefore, our results suggested that ambient climate may control the directional trend of community stability while manipulative treatments may determine the temperature sensitivity of the response of community stability to climate relative to the ambient treatment. Our study emphasizes the importance of the context dependency of the response of community stability to human-driven environmental changes.

Warming homogenizes apparent temperature sensitivity of ecosystem respiration.

Warming-induced carbon loss through terrestrial ecosystem respiration (Re) is likely getting stronger in high latitudes and cold regions because of the more rapid warming and higher temperature sensitivity of Re (Q 10). However, it is not known whether the spatial relationship between Q 10 and temperature also holds temporally under a future warmer climate. Here, we analyzed apparent Q 10 values derived from multiyear observations at 74 FLUXNET sites spanning diverse climates and biomes. We found warming-induced decline in Q 10 is stronger at colder regions than other locations, which is consistent with a meta-analysis of 54 field warming experiments across the globe. We predict future warming will shrink the global variability of Q 10 values to an average of 1.44 across the globe under a high emission trajectory (RCP 8.5) by the end of the century. Therefore, warming-induced carbon loss may be less than previously assumed because of Q 10 homogenization in a warming world.

Domestication Shapes the Community Structure and Functional Metagenomic Content of the Yak Fecal Microbiota.

Domestication is a key factor of genetic variation; however, the mechanism by which domestication alters gut microbiota is poorly understood. Here, to explore the variation in the structure, function, rapidly evolved genes (REGs), and enzyme profiles of cellulase and hemicellulose in fecal microbiota, we studied the fecal microbiota in wild, half-blood, and domestic yaks based on 16S rDNA sequencing, shotgun-metagenomic sequencing, and the measurement of short-chain-fatty-acids (SCFAs) concentration. Results indicated that wild and half-blood yaks harbored an increased abundance of the phylum Firmicutes and reduced abundance of the genus Akkermansia, which are both associated with efficient energy harvesting. The gut microbial diversity decreased in domestic yaks. The results of the shotgun-metagenomic sequencing showed that the wild yak harbored an increased abundance of microbial pathways that play crucial roles in digestion and growth of the host, whereas the domestic yak harbored an increased abundance of methane-metabolism-related pathways. Wild yaks had enriched amounts of REGs in energy and carbohydrate metabolism pathways, and possessed a significantly increased abundance of cellulases and endohemicellulases in the glycoside hydrolase family compared to domestic yaks. The concentrations of acetic, propionic, n-butyric, i-butyric, n-valeric, and i-valeric acid were highest in wild yaks. Our study displayed the domestic effect on the phenotype of composition, function in gut microbiota, and SCFAs associated with gut microbiota, which had a closely association with the growth performance of the livestock. These findings may enlighten the researchers to construct more links between economic characteristics and gut microbiota, and develop new commercial strains in livestock based on the biotechnology of gut microbiota.

Environment and host species identity shape gut microbiota diversity in sympatric herbivorous mammals.

The previous studies have reported that the mammalian gut microbiota is a physiological consequence; nonetheless, the factors influencing its composition and function remain unclear. In this study, to evaluate the contributions of the host and environment to the gut microbiota, we conducted a sequencing analysis of 16S rDNA and shotgun metagenomic DNA from plateau pikas and yaks, two sympatric herbivorous mammals, and further compared the sequences in summer and winter. The results revealed that both pikas and yaks harboured considerably more distinct communities between summer and winter. We detected the over-representation of Verrucomicrobia and Proteobacteria in pikas, and Archaea and Bacteroidetes in yaks. Firmicutes and Actinobacteria, associated with energy-efficient acquisition, significantly enriched in winter. The diversity of the microbial community was determined by the interactive effects between the host and season. Metagenomic analysis revealed that methane-metabolism-related pathway of yaks was significantly enriched in summer, while some pathogenic pathways were more abundant in pikas. Both pikas and yaks had a higher capacity for lipid degradation in winter. Pika and yak shared more OTUs when food shortage occurred in winter, and this caused a convergence in gut microbial composition and function. From winter to summer, the network module number increased from one to five in pikas, which was different in yaks. Our study demonstrates that the host is a dominant factor in shaping the microbial communities and that seasonality promotes divergence or convergence based on dietary quality across host species identity.