Barley farmland harbors a highly homogeneous soil bacterial community compared to wild ecosystems in the Qinghai-Xizang Plateau

Introduction Understanding patterns and processes of microbial biogeography in soils is important for monitoring ecological responses to human activities, particularly in ecologically vulnerable areas such as the Qinghai-Xizang Plateau. Highland barley is the staple food of local people and has mainly been cultivated along the Yarlung Zangbo River valley in Xizang. Methods Here we investigated soil bacterial communities from 33 sampling sites of highland barley farmland in this region and compared them to those from wild ecosystems including alpine tundra, meadow, forest, and swamp. Additionally, the effects of environmental factors on bacterial communities, as well as the relative importance of stochastic and deterministic processes in shaping the beta diversity of soil bacterial communities in alpine ecosystems were assessed. Results In contrast to soils of wild ecosystems, these farmland samples harbored a highly homogeneous bacterial community without significant correlations with geographic, elevation, and edaphic distances. Discriminant bacterial taxa identified for farmland samples belong to Acidobacteria, with Acidobacteria Gp4 as the dominant clade. Although Acidobacteria were the most abundant members in all ecosystems, characterized bacterial taxa of meadow and forest were members of other phyla such as Proteobacteria and Verrucomicrobia. pH and organic matter were major edaphic attributes shaping these observed patterns across ecosystems. Null model analyses revealed that the deterministic assembly was dominant in bacterial communities in highland barley farmland and tundra soils, whereas stochastic assembly also contributed a large fraction to the assembly of bacterial communities in forest, meadow and swamp soils. Discussion These findings provide an insight into the consequences of human activities and agricultural intensification on taxonomic homogenization of soil bacterial communities in the Qinghai-Xizang Plateau.


Introduction
Microbial biogeography investigates the geographical variation of microorganisms across space and time, and uncovers the processes that cause this variation in microbial assemblage (Hanson et al., 2012;Xu et al., 2020;Wan et al., 2023).With the revolutionary advance in molecular phylogenetic approaches, sequencing techniques and computational methods in recent years (Kozich et al., 2013;Larin et al., 2024), extensive studies have successfully established the biogeographic patterns of bacteria across a broad spectrum of ecosystems (Nemergut et al., 2013).Now the major theme of this field has focused on identifying the mechanisms that shape these patterns (Hanson et al., 2012;Nemergut et al., 2013;Ruff et al., 2015;Cao et al., 2016;Kumar et al., 2017).
It has been acknowledged that edaphic variables, such as pH (Fierer and Jackson, 2006;Cao et al., 2016), salinity (Feng et al., 2023), and carbon availability (Fierer et al., 2007;Walsh et al., 2016), have a pivotal role in determining microbial biogeography.For example, it has been reported that the relative abundance of Acidobacteria, ubiquitous and abundant members of soil bacterial communities (Kielak et al., 2016;Delgado-Baquerizo et al., 2018), is negatively correlated with pH (Wang et al., 2023).These findings support a classic deterministic process of selection via abiotic features in determining microbial community assemblage.However, variation of microbial communities can also be affected by stochastic processes such as dispersal (Hawkins and Zeglin, 2022), evolutionary diversification (Shapiro et al., 2012), and ecological drift (Jones et al., 2023).Given the irreplaceable fundamental functions of microorganisms in biosphere, understanding both patterns and processes of microbial biogeography is not only helpful for us to investigate forces generating biodiversity, but also crucial for monitoring ecological responses to environmental changes caused by human activities.
Highland barley is the staple food for local people and an important livestock feed in the Qinghai-Xizang Plateau, which is one of the domestication centers for cultivated barley (Zeng et al., 2015).To improve the sustainable development of highland barley farmland in the ecologically vulnerable Qinghai-Xizang Plateau, it is important to know the characteristics of soil bacterial communities in the field.To this end, we firstly determined the beta diversity of soil bacterial communities in the highland barley farmland across 33 sampling sites along the Yarlung Zangbo River valley in the Qinghai-Xizang Plateau.The identified biogeographical pattern of soil bacterial communities and processes shaping it were compared to those of 33 soil samples in wild ecosystems (including tundra, swamp, meadow and forest), which are of different levels of ecological fragility.

Sampling, DNA extraction, soil physicochemical properties characterization
Highland barley farmland soils were collected from 33 sampling sites in Xizang, China, with 12 samples collected along the Nyang Qu River (an upper tributary of Yarlung Zangbo River), 12 samples from the downstream of the Yarlung Zangbo River (where the Lhasa River joined), and nine samples from the midstream of the Yarlung Zangbo River (Supplementary Table S1 and Figure 1).Non-farmland soils were collected from 30 sites across eight alpine mountains (Mt. Bangla, Mt. Gangbala, Mt. Kaluola, Mt. Luobadui, Mt. Mila, Mt. Sejila, Mt. Pading and Mt. Zhongla) in the same region (with fifteen alpine meadow sites, nine alpine tundra sites, and six high-frigid forest sites), and three sites in Aralake swamp, a dry lakebed located in the Everest Wetland Nature Reserve (Supplementary Table S1 and Figure 1).Samples were collected during the growing season when the barley starts sprouting and the alpine plants flourishing and flowering.At each site, three areas (0.5 m x 0.5 m) were probed, with dominant vegetation recorded in situ.Soil samples (0-15 cm) from three areas were randomly collected and mixed as a single sample for each sampling site.Soil samples were stored in ice bags, transported to laboratory and preserved at −80°C before DNA extraction.Total DNA of soil samples were extracted using the Fast DNA™ SPIN Kit (MP biomedicals, Cleveland, OH, United States).The quality and quantity of DNA were checked with a NanoDrop device (ND-1000, Thermo Fisher, United States).Moreover, total N, available phosphorus (P), organic matter (OM), available potassium (K), electrical conductivity (EC), pH, and ion concentrations of Ca 2+ , Mg 2+ , Na + , Cl − , K + , HCO 3 − , and SO 4 2− for air-dried soil samples were determined at the Plant Nutrient and Resource Research Institution, Beijing Academy of Agriculture and Forestry Sciences, China (Supplementary Table S1), using the methods described before (Wang et al., 2018).Non-parametric pairwise comparison (Dunn's Post-Hoc test) between the 5 alpine ecotypes showed no significant difference based on the Euclidean distances of geographic (Supplementary Figure S1A) and elevation (Supplementary Figure S1B) characteristics, however, the edaphic distances of farmland are significantly lower than that of meadow and tundra (Supplementary Figure S1C, p-values = 9.85e−04 and 0.02, respectively).Climate data such as mean annual temperature (AMT) and annual precipitation (AP) of the sampling region were acquired from DIVA-GIS1 at the 2.5 arc-minutes resolution level.AMT and AP were then predicted and cross validated with linear model using latitude, longitude and altitude as explanatory variables in R before interpolating into sampling sites (Supplementary Table S1).

Library construction for 16S rRNA gene sequencing
PCR amplifications were performed on 66 DNA samples using the 515F/806R primer set (515F: GTGYCAGCMGCCGCGGTAA and 806R: GGACTACNVGGGTWTCTAAT), which targets the V4 region of the 16S rRNA gene (Caporaso et al., 2012).The PCR products with P5/P7 adapters and double indexes were purified with double volume of AMPure XP beads, and fluorescence quantified (Qubit 2.0 Fluorometer, ThermoFisher Scientific, United States) with Qubit dsDNA HS assay kit (Invitrogen, United States).Multiple samples were pooled together in equimolar concentrations.Then 10 ng of pooled PCR products were sequenced on the Illumina MiSeq platform (paired-end reads of 300 bp; Illumina, United States) at the National Human Genome Centre of China in Shanghai, China, according to the manufacturer's manual.

Sequence processing and bacteria diversity analysis
Raw 16S rRNA gene fastaq files were processed using Mothur v.1.35.1 according to Miseq SOP (Kozich et al., 2013).Reads were firstly combined and screened to reduce sequencing errors, then aligned against the SILVA-based bacterial reference (Release 138.1).The aligned reads were further de-noised and chimeras were removed using the UCHIME algorithm (Edgar et al., 2011).Any sequences not belonging to bacterial 16S rRNA gene, such as homologs from archaea, chloroplasts, and mitochondria, were classified using the Bayesian classifier and removed before downstream analysis.The optimized sequences were then normalized to 4,052 sequences per sample (the smallest sequencing depth for tested samples) before calculating taxonomy-based diversity and phylogeny-based diversity using Mothur and FastTree, respectively, (Price et al., 2009;Schloss et al., 2009).
For taxonomy-based analysis, the quality-controlled subsample sequences were clustered into operational taxonomic units (OTUs) that share ≥97% sequence similarity.Rare OTUs constituting less than 0.001% of the total normalized filtered sequences were removed, for rare OTUs may result in inflated estimates of diversity (Bokulich et al.,  Wang et al. 10.3389/fmicb.2024.1418161Frontiers in Microbiology 04 frontiersin.org2013).Taxonomic information was then assigned for the remaining OTUs.For phylogeny-based analysis, the subsample sequences with rare OTUs removed were analyzed using FastTree software to compute an approximately-maximum-likelihood phylogenetic tree (Price et al., 2009).The phylogenetic tree of the dominant clade of bacterial orders and their relative abundance across ecosystems were visualized in EvolView (Zhang et al., 2012).Alpha diversity was analyzed using the number of OTUs observed, the inverse Simpson D (Simpson, 1949), as well as Faith's phylogenetic diversity (PD) (Helmus et al., 2007) by randomly selecting 3,000 sequences per sample 1,000 times and calculating the average (Supplementary Table S2).Beta diversity was calculated either using the taxonomy-based Bray-Curtis dissimilarity metric or the phylogeny-based weighted UniFrac metric generated from the filtered sequences (evenly sampled at 3,000 reads per sample with 1,000 iterations).The resulting distances were further visualized by using GenGIS geographical mapping analysis (Parks et al., 2009).Analysis of similarity statistics (ANOSIM) (Clarke, 1993) and homogeneity of molecular variance (HOMOVA) (Stewart and Excoffier, 1996) were used to test significant differences in composition and variation of bacterial communities.
Discriminant bacteria from phylum to family level were identified for different ecosystems using linear discriminant analysis effect size (LEfSe) (Segata et al., 2011).Taxa with significant differential abundance (Kruskal-Wallis sum-rank test, p-value <0.05; LEfSe >4) were used to generate taxonomic cladograms illustrating differences between Xizang ecosystems.Indicator OTU analysis (Dufrene and Legendre, 1997) was used to identify potential specialist for each ecosystem at the OTU level.This approach can uncover OTU, which is present in most sites of a group and has a high relative abundance within that group compared to the other groups, with statistical significance evaluated using a randomization procedure.OTUs with IV > 0.5 (indicator value) and p-value <0.05 were classified as group specialists.The high abundant indicators (OTUs representing larger than 0.05% of the total filtered and normalized sequences) identified by this method were selected for species assemblages analysis using Kendall's W coefficient of concordance (Legendre, 2005).An overall test of independence of all abundant indicator OTUs was first executed.If the null hypothesis was refused, correlated species were divided into groups and permutation tests were used to test the contribution of each OTU to the overall statistic within each group.OTUs representing larger than 0.05% of the total subsample sequences were used in the analysis to avoid inflated estimates of indicator OTUs caused by rare taxa.

Effects of environmental factors on bacterial communities
The variation in phylogenetic diversity was measured using nullmodel-based phylogenetic β-diversity metrics (βNTI) to classify community pairs into underlying drivers of deterministic and stochastic processes, as described in previous studies (Dini-Andreote et al., 2015;Jiao et al., 2022).Briefly, when βNTI is less than −2, it suggests a significantly lower phylogenetic turnover than expected under a null model of community assembly, indicating the presence of homogeneous selection, where similar species tend to coexist.Conversely, a βNTI greater than 2 indicates significantly higher phylogenetic turnover, suggesting variable selection, where dissimilar species coexist.If the absolute value of βNTI is less than 2, it indicates that stochastic processes, such as dispersal, drift, or speciation, play a dominant role in shaping the community structure.Mantel and Partial Mantel tests were performed based on weighted UniFrac distance and Euclidean distances (elevation, geographic, or edaphic distances) to study the correlation between environmental factors and soil bacterial communities.After removing highly correlated environment variables (Spearman ρ > 0.8), Hellinger transformed OTU abundance data were projected to environmental variables: geographic (G), elevation (A), edaphic attributes (E) for correlation analysis using redundancy analysis (RDA).The variance degree of bacterial taxa accounted by those three explanatory variables and their combined effects were further determined using the partitioning functions based on RDA analysis.These analyses were all carried out by using functions in the vegan package in R software (v 3.1.3and v 3.2.2) (Legendre and Legendre, 2012;Oksanen et al., 2016).Partial correlation analysis was used to control the variation of other edaphic factors when the correlation between pH/OM and relative abundance of Acidobacteria was studied.

Homogenous bacterial community in the barley field in contrast to other alpine ecosystems
From 66 soil samples of different ecosystems along the Yarlung Zangbo River valley (Figure 1A), 745,451 quality-trimmed sequences (an average of 253 bp) were generated with a reduction of 27.5% of total raw sequencing reads.In total, 5,329 OTUs were kept after removing rare OTUs and the number of observed OTUs per sample ranged from 352 to 779 (Supplementary Table S2).The highland barley farmland soils harbored the lowest variation in bacterial community structure as measured by using weighted UniFrac distance (Figure 1B and Supplementary Figure S1D) compared to those of forest (Dunn' s Post-Hoc test, p-values = 1.84e−04) and meadow (p-values = 3.72e−24), but showed no significant difference to those of tundra (p-values >0.05).HOMOVA (Table 1) showed that the amount of bacterial beta diversity within farmland samples was different from those of forest and meadow (all p-values <0.05), but not from those of tundra and swamp (all p-values >0.05).ANOSIM showed that the bacterial composition of farmland soils was different from those of wild ecosystems (all p-values <0.01;Table 1).Dissimilarity of bacterial composition between different wild ecosystems was also uncovered for 5/6 pairs (all p-values <0.05), except the meadowforest pair (p-value >0.05).
The variation patterns of the soil bacterial community within or between ecosystems were further analyzed using geospatial biodiversity mapping based on the pairwise weighted UniFrac distance (Figure 1C).Generally, independent of geographical distance between samples (Figure 1A), farmland samples formed a homogeneous cluster, which is closer to the alpine tundra cluster than the intermingled forest-meadow branches.A significant distance decay of community similarity was revealed for samples of tundra (Supplementary Figure S2, R 2 = 0.19, p-value <0.01) and meadow (R 2 = 0.12, p-value <0.001), but not for farmland (R 2 = 0.0056, p-value >0.05) and forest samples (R 2 = 0.0526, p-value Frontiers in Microbiology 05 frontiersin.org>0.05).Swamp samples were not analyzed herein due to the small sample size and their immediate vicinity (less than 10 meters' distance).In short, soil bacterial communities in the highland barley farmland were highly homogenous along the Yarlung Zangbo River valley compared to those of wild ecosystems investigated herein in the same region.

3.2
The relative importance of stochastic and deterministic processes in shaping the beta diversity of soil bacterial communities in alpine ecosystems in Xizang Null model analyses revealed that the deterministic assembly was dominant in bacterial communities in highland barley farmland (93%) and tundra (84%) soils, whereas stochastic assembly also contributed a large fraction to the assembly of bacterial communities in forest (41%), meadow (40%) and swamp (67%) soils (Figure 2).In particular, the assembly of bacterial communities in highland barley farmland and tundra soil was governed by the homogeneous selection (βNTI < −2, Figure 2), in accordance with the homogenous bacterial community in these two ecosystems contrast to other alpine ecosystems.The relationships between phylogenetic diversity and major factors were also used to infer changes in the relative influences of deterministic and stochastic assembly processes, which was assessed using the Mantel tests that correlated the weighted UniFrac distance matrices with the Euclidean distance matrices of environmental variables.The highly correlated environmental variables (Spearman ρ > 0.8) among soil physicochemical properties, climatic and topographical features (Supplementary Table S1) were removed.The final variables used were as follows: latitude, longitude, altitude, organic matter, available phosphorus, available potassium, electrical conductivity, pH, Na + , Cl − and SO 4 2− . All selected environmental features except K + and SO 4 2− exhibited significant differences between ecosystems (ANOVA, all p-values <0.01, Supplementary Table S3).
The Mantel and Partial Mantel tests revealed that the beta diversity, based on weighted UniFrac distance (Table 2), within farmland samples was not correlated with the distance matrix of geographic (latitude and longitude), elevation (altitude) or edaphic characteristics (all p-values >0.05), indicated that the deterministic assembly of bacterial communities of highland barley farmland could not be explained by measured environmental variables.By contrast, the beta diversity of tundra samples correlated with edaphic distance (ρ = 0.396, p-value <0.05) in the Mantel test, and after controlling for geographic (ρ = 0.359, p-value <0.05) and elevation distances (ρ = 0.397, p-value <0.05) in the Partial Mantel test, indicating significant effects of deterministic factors.Similarly, the beta diversity of meadow samples correlated with geographic (ρ = 0.366, p-value <0.01) and edaphic distances (ρ = 0.510, p-value <0.001) in the Mantel test.The beta diversity of forest samples showed no correlation with elevation, geographic and edaphic factors (Table 2, all p-values >0.05), despite its high level of similarity to that of meadow samples (Figure 1B).
Notably, Acidobacteria_Gp4, Acidobacteria_Gp6 and Spartobacteria dominated bacterial community composition, from class to family level.These three clades accounted for 56.9% of the total sequence in all soil samples (Figure 3B).Acidobacteria_Gp4 was 31.6 and 40.9% of the total sequences in farmland and tundra, respectively.Acidobacteria_Gp6 accounted for 26.2 and 27.2% of the total sequences in farmland and swamp, respectively.Spartobacteria was 26.8, 20.8 and 13.1% of the overall sequences in forest, alpine meadow and tundra, respectively (Figure 3B).Although Acidobacteria_Gp4 dominated in both farmland and  RDA showed that 31.6% (based on adjusted R 2 , p-value <0.001) of the variance of bacterial communities among tested samples could be explained by those selected environmental and spatial factors.Variation partitioning analysis (Figure 4A) further revealed that edaphic factors alone captured 12.7% of variance (p-value <0.001), while very limited variation could be attributed to pure effects of altitude (0.9%, p-value = 0.02) and geographic distance (2.5%, p-value = 0.316).By contrast, the combined effects of edaphic factors with either altitude (23.3%) or geographic distance (24.9%) apparently may play a greater role in shaping the observed beta diversity (Figure 4A).As shown in the RDA analysis (Figure 4B), pH and organic matter were the major environmental factors along the first axis (22.5% of total variance, p-value <0.001).The farmland cluster was separated from those of wild ecosystems, particularly meadow and forest, along these two edaphic factors.The meadow and forest samples were enriched with certain OTUs belonging to Proteobacteria and Verrucomicrobia while the farmland samples harbored more OTUs of Acidobacteria (Figure 4B).Significant correlations between some OTUs and other individual edaphic factors such as phosphorus, SO 4 2− , EC, Na + and Cl − were also observed (Figure 4B).
In addition to the strong positive correlation with pH, relative abundance of Acidobacteria showed lower levels of correlation with EC (ρ = 0.414), available phosphorus (ρ = 0.365), Na + (ρ = 0.297) and SO 4 2− (ρ = 0.288) (all p-values <0.05).When these edaphic factors were individually controlled, the partial correlation coefficient of pH and relative abundance of Acidobacteria ranged from 0.65 to 0.8 (all p-values <0.001).If the other edaphic factors were all controlled, the coefficient was 0.41 (p-value = 0.0007).This implies that the positive correlation between pH and relative abundance of Acidobacteria was robust in our test samples.By contrast, the partial correlation coefficient was −0.17 (p-value >0.05) between organic matter and Acidobacteria abundance when the other edaphic factors were all controlled, indicating a combined effect by organic matter and other edaphic attributes.Deterministic factors shaping the bacterial communities.(A) Variation partition analysis of OTUs among geographic ("G"), elevation ("A"), edaphic attributes ("E"), and their interactions ("G x A," "G x E," and "E x A").U, Unexplained variation (residual).The values in parentheses are p-values.(B) RDA of OTUs and environmental attributes.The percentage of variation explained by RDA1 and RDA2 is shown, both of the eigenvalues (RDA1 and RDA2) and all of the environmental variables are significant in permutation tests (α = 0.05).OM, organic matter; phosphorus, available phosphorus; potassium, available potassium; EC, electrical conductivity.Understanding patterns and processes of bacterial biogeography is an important way to evaluate the ecological effect caused by human activities (Geremia et al., 2016;Wu et al., 2021;Xue et al., 2021).Soil biodiversity is strongly influenced by external drivers such as climate change and nitrogen deposition but also by land-use management (Wall et al., 2015).Evidence is mounting that soil biodiversity loss and simplification of soil community composition impair multiple ecosystem functions, including plant diversity, decomposition, nutrient retention, and nutrient cycling (Wagg et al., 2014).The Qinghai-Xizang Plateau, an ecologically vulnerable region, has undergone studies investigating soil bacterial communities in individual wild ecosystems like alpine meadows and lake sediments (Xiong et al., 2012;Yang et al., 2014;Yuan et al., 2014;Wang et al., 2022;Yang et al., 2022).In the present study, we compared soil bacterial communities in the highland barley farmland to those of wild alpine ecosystems including tundra, meadow, forest and swamp along the Yarlung Zangbo River valley, which is the major region for highland barley cultivation.
Soil bacterial communities were highly homogeneous in the highland barley farmland compared to those of meadow, forest and swamp.The observed low beta diversity in farmland samples might be partially due to the relatively simple vegetation coverage by highland barley compared to that of meadow, forest and swamp.This view is supported by a report of soil microbial communities along four sites of a Qinghai-Xizang mountainous grassland (Yang et al., 2014), where vegetation diversity was revealed as one of the major factors shaping the variations of microbial communities.Although farmland samples showed a similar level of beta diversity to those from tundra, no significant relationship was found between the beta diversity of farmland samples and their corresponding geographic, elevation and edaphic distances.This is in contrast to the observations in meadow samples in this study and earlier reports on bacterial communities in lake sediments and alpine meadows in the Qinghai-Xizang Plateau (Xiong et al., 2012;Yang et al., 2014;Yuan et al., 2014).It might be reasonable to speculate that transport via human activities, such as those in the farmland, could be more efficient than dispersal via wind, water, and other hitchhiking mechanisms in wild ecosystems (Papke and Ward, 2004).Recent studies demonstrated that land conversion to agricultural land results in taxonomic and functional homogenization of soil bacteria and fungi, mainly driven by the increase in the geographic ranges of taxa in croplands (Peng et al., 2024) and disproportionate negative effects of agricultural cultivation on rare microbial taxa (Banerjee et al., 2024).These results align with our perspective, indicating that agricultural activities lead to microbial homogenization and has profound implications for the structure and function of soil ecosystems, highlighting the importance of agricultural sustainability in protecting biodiversity.
The dominant phylum was Acidobacteria in all ecosystems tested in this study, and this is in line with earlier studies of soil bacterial communities across continents suggesting Acidobacteria as ubiquitous and abundant members (Zinger et al., 2009;Yuan et al., 2014;Sun et al., 2020;Sui et al., 2022).Moreover, farmland, tundra and swamp soils characterized with higher pH and lower organic matter than meadow and forest soils had a higher abundance of Acidobacteria.This is consistent with a negative correlation between abundance of Acidobacteria and organic carbon availability, and a strong positive correlation between pH and Acidobacteria abundance.These results are inconsistent with an earlier report in North and South America soils (Jones et al., 2009) but could be partially explained at the subgroup level of Acidobacteria, since subgroups such as 1, 2, 3, 12 and 13 which have negative correlations with pH (Jones et al., 2009) were less abundant in test samples of this study than those subgroups (4, 6, 7, 16, and 17) that have positive correlations with pH.Therefore, the general conclusion of the correlation between Acidobacteria abundance and pH/organic carbon availability proposed earlier (Fierer et al., 2007;Jones et al., 2009) might be misleading.These findings suggest that more studies should be compared to improve our understanding of the ecological significance of Acidobacteria (Kielak et al., 2016), and it is important to investigate bacterial communities at a high-resolution level (Liao et al., 2023).
Agricultural intensification can reduce the diversity and density of beneficial organisms; less-intensive practices such as enhancing crop diversity by rotation and intercropping, and reducing soil tillage are beneficial to soil biodiversity and the sustainability of agriculture (Wall et al., 2015;Bender et al., 2016).In this study, a considerable similarity level of bacterial communities between farmland and tundra soils was found.Tundra is characterized by its lower vegetation coverage than meadow and forest, and represents one of the most vulnerable ecosystems on the planet (Malanson et al., 2009;Sistla and Schimel, 2013;Stark et al., 2015).Future action should be taken to increase the crop diversity of barley farmland system by using rotation and intercropping crops (Bardgett and van der Putten, 2014;Bender et al., 2016;Blaser et al., 2016).Moreover, similar edaphic conditions (such as pH and organic matter) between the highland barley farmland and tundra suggest that the nutrient management of highland barley farmland soils should be carried out.Particularly pH and organic matter representing the edaphic factors with major effects on soil bacterial community could be managed in agriculture practices (van Rijssel et al., 2022;Wang et al., 2023).While this study provides valuable insights into the bacterial biogeography of highland barley farmland and wild alpine ecosystems, it has some limitations.Firstly, our study focused on the taxonomic biogeography of bacterial communities alone, a more comprehensive analysis incorporating multiple microbial groups (e.g., bacteria, fungi and archaea) and microbial function could provide a broader perspective on the effects of land-use change.Secondly, here we focused on a relatively small geographical scale in the Qinghai-Xizang Plateau, limiting the generalizability of the findings.Broader spatial sampling across different ecoregions and elevations could provide a more comprehensive understanding of bacterial biogeography patterns in alpine ecosystems.

Conclusion
An important aim of studies of soil bacterial community is to improve our understanding of the effects of human activities on the sustainability of agricultural soils.Our study highlights the distinctiveness of soil bacterial communities in highland barley farmland compared to those in wild alpine ecosystems along the Yarlung Zangbo River valley.The homogeneity observed in farmland bacterial communities, despite their similarity to tundra samples, but with no detectable correlations with abiotic effectors, suggests that agricultural activities in highland 10. 3389/fmicb.2024.1418161Frontiers in Microbiology 11 frontiersin.orgbarley farmland may be having a profound effect on soil microbial diversity.This finding raises concerns about the potential ecological consequences of agricultural practices in this region, especially given the already fragile nature of the Qinghai-Xizang Plateau.

FIGURE 1
FIGURE 1 Beta diversity analysis of bacterial community.(A) UPGMA dendrogram based on geographic distance.Photos show plants of Xizang Highland barley, and the landscape of alpine tundra and meadow.(B) Density plot showing the distribution of weighted UniFrac distance (mean ± SD). (C) Geospatial biodiversity mapping analysis comparing bacterial communities across Xizang biomes using weighted UniFrac distance.Branches are colored based on the ecosystems they represent.Pie charts represent relative abundance of bacteria taxa at phylum level for each site.Geographic origins of soil samples from wild ecosystems are shown.

FIGURE 2
FIGURE 2The distributions of between community analog (βNTI) across alpine ecosystems.The horizontal dashed lines indicate the βNTI values of +2 and − 2. Pie plots shows the relative contribution of each ecological process to community assembly.

FIGURE 5
FIGURE 5Relationship between the relative abundance of Acidobacteria and pH or organic matter in 66 soil bacterial communities.(A) The relative abundance of each Acidobacteria subgroup in each sample (ranging from 0 to 66.0%).OM, organic matter.(B,C) The effect of OM (B) or pH (C) on the abundance of Acidobacteria.(D,E) The effect of pH on the abundance of Acidobacteria_Gp4 (D) and Acidobacteria_Gp6 (E).Samples were colored according to different ecosystems: farmland (blue), forest (green), meadow (red), swamp (purple), tundra (pink).Results from the research byJones et al. (2009) were included in our analysis (black dots).Both Spearman's rank correlation (ρ) and regression analysis (R 2 ) were used for testing the correlations, all p-values <0.001.

TABLE 1
Analysis of similarities (ANOSIM) and homogeneity of molecular variance (HOMOVA) of bacterial beta-diversity.

TABLE 2 Mantel
and partial Mantel tests for the correlation between Weighted UniFrac distances of soil bacterial communities and the Euclidean distances of all explanatory variables.Elevation, geographic and edaphic distances were analyzed herein using Spearman's rho for different ecosystems and all samples.***p-value < 0.001; **p-value < 0.01; *p-value < 0.05 (based on 10,000 permutations).NA, swamp samples are not analyzed due to the small sample size.All explanatory distances were measured as Euclidean distance using all variables standardized to have a mean of zero.