Temporal Dynamics of Rhizobacteria Found in Pequin Pepper, Soybean, and Orange Trees Growing in a Semi-arid Ecosystem

Harsh environmental conditions in drylands force plants and their associated microbial communities to adapt to abiotic stresses. In semi-arid environments, climatic conditions and poor agricultural management have a strong impact on plant yield and thus, enhancing soil fertility by means of beneficial microorganisms such as plant growth-promoting rhizobacteria (PGPR) has been proposed as part of sustainable agricultural management. As drylands will increase due to climate change, studying microbial community dynamics of crops under such conditions is crucial as it might favor rhizobacteria adapted to drought. While the microbiome of many native dryland crops has been characterized, the microbial community composition from non-native crops under semi-arid environmental conditions is understudied. Thus, the aim of this study was to characterize the bacterial community associated with the roots of three crops with different growth cycles, cultivated in the same semi-arid environment, to understand their microbial community composition during the season with the highest temperature in northeast Mexico. We performed high throughput sequencing of the V3-V4 region of the 16S rRNA gene from root samples of Pequin pepper, soybean and orange trees. Classified taxa were evaluated according to crop, sampling time and climatological parameters. Our findings revealed that changes in temporal dynamics of microbial communities correlate with environmental temperature. Moreover, the microbial community of pepper was more diverse and differed from that of soybean and orange. Regarding PGPR, 47.6% of the genera were shared among crops with a high relative abundance of Bacillus, but we also detected crop-specific microbial associations where Serratia was specific to orange trees and Rhodobacter to pepper. When analyzing PGPR in correlation to climatological parameters, Bacillus was found to thrive under lower precipitation rates, higher temperatures and higher evaporation rates in pepper and orange. In contrast, some PGPR commonly used in commercial biofertilizers such as Rhizobium and Azospirillum were affected by high temperatures. This study provides a better understanding of the rhizobacterial assemblies of economically relevant crops grown under a semi-arid environment.

Harsh environmental conditions in drylands force plants and their associated microbial communities to adapt to abiotic stresses. In semi-arid environments, climatic conditions and poor agricultural management have a strong impact on plant yield and thus, enhancing soil fertility by means of beneficial microorganisms such as plant growth-promoting rhizobacteria (PGPR) has been proposed as part of sustainable agricultural management. As drylands will increase due to climate change, studying microbial community dynamics of crops under such conditions is crucial as it might favor rhizobacteria adapted to drought. While the microbiome of many native dryland crops has been characterized, the microbial community composition from non-native crops under semi-arid environmental conditions is understudied. Thus, the aim of this study was to characterize the bacterial community associated with the roots of three crops with different growth cycles, cultivated in the same semi-arid environment, to understand their microbial community composition during the season with the highest temperature in northeast Mexico. We performed high throughput sequencing of the V3-V4 region of the 16S rRNA gene from root samples of Pequin pepper, soybean and orange trees. Classified taxa were evaluated according to crop, sampling time and climatological parameters. Our findings revealed that changes in temporal dynamics of microbial communities correlate with environmental temperature. Moreover, the microbial community of pepper was more diverse and differed from that of soybean and orange. Regarding PGPR, 47.6% of the genera were shared among crops with a high relative abundance of Bacillus, but we also detected crop-specific microbial associations where Serratia was specific to orange trees and Rhodobacter to pepper. When analyzing PGPR in correlation to climatological parameters, Bacillus was found to thrive under lower precipitation rates, higher temperatures and higher evaporation rates in pepper and orange. In contrast, some PGPR commonly used in commercial biofertilizers such as Rhizobium and Azospirillum were affected by high temperatures. This study provides a better understanding of the rhizobacterial assemblies of economically relevant crops grown under a semi-arid environment.

INTRODUCTION
Drylands represent around 41% of the total Earth's land surface, being mostly composed by semi-arid regions (Huang et al., 2016) which are characterized by adverse environmental conditions that urge plants and other organisms to adapt in order to alleviate abiotic stress (Soussi et al., 2016). In drylands, plants are prone to experience high temperature, drought, high salinity and osmotic stress. Surviving these climatic conditions requires plants to adapt and acclimate, overall impacting crop yield (Vimal et al., 2017) and consequently, affecting the global production of cereals, oil crops, fruits and vegetables, and experiencing hundreds of billion-dollar losses (Dwivedi et al., 2018). In this environment, plant's growth is limited due to low nutrient availability caused by decreased primary productivity, high nutrient immobilization and loss by soil erosion (Cui et al., 2019).
Additionally, interactions of climate change and human activities have a strong impact on semi-arid locations. In recent years the intensive practices of modern agriculture have brought awareness of their environmental impact, for instance in greenhouse gas emissions and nutrient leaching. Besides, climate change has contributed to an expansion of the drylands all over the world with semi-arid regions having the largest growth (Huang et al., 2016). More than ever, finding sustainable alternatives to fertilize and protect plants against phytopathogens while increasing food production is extremely important (Ferreira et al., 2019). Therefore, practices that aim at enhancing soil fertility and at the same time maintain crop yield by supporting diverse microbiota from the rhizosphere have been proposed (Hartman et al., 2018).
The rhizosphere, which is rich in nutrients from root exudates, represents a stable and protected interface with the adjacent soil. Its microbiota is composed of bacteria and other microorganisms that collectively regulate ecosystem processes and play key roles in nutrient cycling, i.e., by decomposing organic matter as well as transforming and fixing soil nutrients (Hartman et al., 2018;Singh et al., 2019). Understanding rhizosphere communities is not trivial, as they are affected by interactions between soil types, plant species and their growth stages (Smalla et al., 2001;Xu et al., 2009;Philippot et al., 2013;Yuan et al., 2015;Qiao et al., 2017). In addition, an increasing amount of research has shown that in drylands the composition of bacterial assemblies is favored by drought adaptation to alleviate plant stress (Soussi et al., 2016). The root microbiome in these conditions is especially sensitive to disturbances. This aspect has been considered key to make drylands particularly prone to desertification caused by poor management and climatic events (Taketani et al., 2017).
Insights on the impact of soil moisture on wheat rhizosphere microbiome in a semi-arid ecosystem have identified significant differences in specific taxa classified as PGPR (Mavrodi et al., 2018). However, current studies that evaluate plant root microbiome on drylands are mainly limited to native crops Fonseca-García et al., 2016;Dastogeer et al., 2018;Khan et al., 2020) and do not focus on introduced species that may require growth enhancers for their production. Hence, to produce biofertilizers specific for crops cultivated in drylands it is necessary to conduct studies of microbiota dynamics in field conditions for each crop.
The aim of this study is to characterize the rhizosphere bacterial community of three crops with different cycles growing under the same field conditions on a semi-arid ecosystem. A perennial crop of Pequin pepper (Capsicum annuum var. glabriusculum) and annual crops of soybean (Glycine max L.) and orange trees (Citrus sinensis) were analyzed to identify crop specific rhizobacteria and the effects of climatic conditions on microbial communities. Identified taxa were evaluated according to sampling time and crop type. The PGPR community was analyzed in correlation with climatological parameters characteristic of drylands. To our knowledge, this is the first study that evaluates the microbiome dynamics from non-native plants with different growth cycles in a semi-arid environment.

Site Selection, Sampling, and DNA Extraction
Pequin pepper, soybean and orange trees were sampled from the Tecnologico de Monterrey experimental field site in Hualahuises,Nuevo Leon,Mexico (24 • 52.50' N,99 • 37.26' W). This site has been described by the Mexican Environmental and Natural Resources Secretary as a semi-arid location based on Penman's aridity index. The sampling was done in May, July and September 2017. The sites planted with soybean and pepper were in close proximity (10 m apart), whereas the orange trees and young peppers were planted ∼50 m from these crops (Figure 1). Two composite soil samples were taken each at 0-30 cm depth and were analyzed by the ISO certificated company Fertilab (Guanajuato, Mexico). The soils are of a medium clay-loam texture with pH 8.1-8.5 and only slight differences were found among the sites regarding soil characteristics (Supplementary Table 1). Lateral roots were collected at a depth of 10 cm from the ground and at least 10 cm of root-length. For Pequin pepper, plants with different ages, young (2 months old) and mature (1 year old), were sampled; for the orange trees, peripheral roots under the tree canopy were taken from approximately a radial length of 2 m, details from each sampling  Table 1. The collected roots were kept in Whirl-Pak bags at 4 • C for transportation. Roots were cleaned using a sieve, washed with tap water and 96% ethanol, clean roots were cut in 1 cm long pieces and stored at −20 • C.
In total, 33 samples (four biological replicates) were selected for DNA extraction performed by the FastDNA Spin Kit for Soil, according to the standard protocol, with the lysing matrix type A, but with one extra ¼ inch ceramic bead to assure complete rupture of the roots (Senés-Guerrero et al., 2014). DNA was stored at −20 • C until sequencing library preparation.

PCR, Library Preparation, and Sequencing
Amplification of the variable regions V3-V4 (∼550 bp, including indexes and adapters) from the 16S rRNA gene was performed using the primers Bakt_341F (5 ′ -CCTACGGGNGGCWGCAG-3 ′ ) and Bakt_805R (5 ′ -GACTACHVGGGTATCTAATCC-3 ′ ), following the Illumina protocol for 16S Metagenomic Sequencing Library Preparation. Each PCR was performed in triplicate per sample and visualized in a 2% agarose gel. Nextera v2 indexes (Illumina, San Diego, CA, USA) were attached to pooled PCR products according to the afore mentioned protocol. Amplicon clean-up was done using AMPure XP  Relative read abundance of genera present in Pequin pepper (young and mature shown separately), orange and soybean. "Un" stands for unclassified to the genus level. Genera represented by < 1% were grouped as "Others".
Frontiers in Sustainable Food Systems | www.frontiersin.org quantification in a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA). Samples were pooled at a concentration of 6.3 pM and loaded together with 30% Phix control into an Illumina R MiSeq sequencer using the MiSeq Reagent Kit v3 (Illumina, San Diego, CA, USA) for paired-end sequencing (2 × 300 bp) carried out at the sequencing facilities of Tecnologico de Monterrey, Campus Monterrey.

Bioinformatic and Statistical Analyses
Raw reads were used as input for QIIME2 version 2019.4 (Bolyen et al., 2019) and a standard pipeline for 16S rDNA amplicon bioinformatics analysis was followed. Briefly, forward and reverse reads were joined and filtered according to Q score prior to denoising with Deblur, obtaining a table with samples grouped in Amplicon Sequence Variants (ASVs). To assign taxonomy, a classifier was constructed using 99% similarity and the GreenGenes database version 13.4. The classifier was trained using the primers and the length of the samples. Taxa bar plots were generated to assign the corresponding taxonomy to the ASV table. Files of phyla and genera were downloaded from view.qiime2.org in CSV format to continue further analysis.
Results were analyzed using RStudio (version 1.1.463; R Core Team, 2019). Rarefaction curves were constructed per sample. Diversity and richness indexes were computed with the vegan package (Oksanen et al., 2019) rarefied to the median read abundance, the samples that had less reads than the median were kept rarefied to the maximum. To analyze temporal dynamics, data was normalized with the DESeq2 package and used as input for perMANOVA, principal coordinate analysis (PCoA) and non-metric multidimensional scaling (NMDS) representing the distribution with 95% interval ellipses and eliminating the genera that were present in less than 10% of the samples. Statistical non-parametric analyses were done using the vegan package and plotted with the package ggplot2 (Wickham, 2009). Unclassified genera were denoted by "Un" and the previous taxonomic level to which they were classified. For relative read abundance, taxa bar plots and presence-absence plots were generated where pepper samples were grouped in young and mature plants. Genera with relative read abundance <1% were classified as "Others." Genera known to promote plant growth were selected as PGPR genera. Additionally, a CIRCOS plot of the PGPR genera was constructed using the package circlize (Gu et al., 2014), relating the presence of PGPR genera to each crop.
To understand the effects of climatic conditions, Mexico's Water National Commission datasets were used to construct boxplots comparing rainfall, evaporation rate, maximal and minimal temperature of the year from which sampling was done (2017) to historic data from 1998 to 2018 (CONAGUA Gobierno de México, 2020). The climatic variables were measured at a climatological station located at 2 km from the experimental field. These climatological parameters were used to analyze their correlation with PGPR genera using Pearson's correlation with the corrplot package (Wei and Simko, 2017). Historic climatic data was also analyzed from 1998 to 2018 to classify the sampling site as semi-arid.

Rhizobacteria Community Composition of Pequin Pepper, Soybean, and Orange Trees
As a first approach to describe temporal dynamics of rhizobacteria of Pequin pepper, soybean and orange trees, we characterized their root endophytic bacterial community at three sampling times in a single field site under semi-arid environmental conditions. Taxonomic classification was assigned at two distinct levels, phylum and genus (Figure 2). Regarding the most abundant phyla, crop microbial composition was highly similar among samples (Figure 2A), where the most abundant phyla were Proteobacteria ranging from 41.1 to 43.6% relative read abundance and Bacteroidetes from 37.2 to 40.7%. Among the tested crops, orange was the only one that possessed Acidobacteria and Chloroflexi in an abundance higher than 1%. At the genus level, the Venn diagram showed a total of 547 genera. Five hundred thirty-four genera were identified in Pequin pepper, 291 in orange tree and 257 in soybean roots ( Figure 2B). Differences in shared genera among crops revealed that the Pequin pepper had 209 specific genera, absent in the other crops, while soybean and orange had only 2 and 10 specific genera, respectively. The three crops shared 209 genera among them. Young and mature Pequin pepper samples were grouped as they did not significantly differ in microbial community composition (Supplementary Figure 2,  Supplementary Table 2). They shared 45.2% of the total amount of genera that represented 55.8 and 57.1% of the relative abundance of all genera, for young and mature plants, respectively. Soybean samples had a lower abundance of unassigned genera of 17.6% compared to pepper that ranged from 43.9 to 44.5% and soybean with 46.1% ( Figure 2C). In addition, orange trees and soybean possessed fewer genera with a predominance of Bacillus (48.1%) in orange and Niastella (28.3%) in soybean.
Additionally, to compare crop bacterial assemblies, as well as the two ages of pepper plants, diversity and richness at the genus level was calculated using the Shannon diversity index and Chao1 richness index. Pequin pepper showed a higher diversity and richness of the microbial root endosphere community compared to soybean but not to orange (P < 0.05) ( Figure 3A); there were no significant differences between young and mature Pequin pepper ( Figure 3B).

Crop Specificity and Temporal Dynamics of Rhizobacteria
Microbial community differences among crops, sampling times and the interaction between these parameters were shown to be statistically significant when analyzing results at the genus level (Table 2). Furthermore, these dissimilarities in microbial communities were visualized by performing NMDS and PCoAs analyses of crop and sampling time (Figure 4) and sampling time for each crop (Figure 5). Both methods yielded similar results, where rhizobacteria community composition of Pequin pepper was significantly different from that of orange and soybean ( Figure 4A).
Regarding temporal dynamics, May samples were significantly different from July and September (Figure 4B). In the same way, differences between sampling times of individual crops resulted in May microbial communities being more dissimilar for Pequin pepper and soybean (Figures 5A,C). However, in the case of the orange trees all sampling times were similar ( Figure 5B). Concerning climatological parameters during sampling, May exhibited in average lower minimum temperature (20.0 ± 3.2 • C) compared to July and September (Supplementary Figure 3). July exhibited the highest average maximum evaporation rate (5.9 ± 2.0 mm; daily average) and lowest rainfall (46.1 mm; accumulated for the month), whereas September had the highest average rainfall (180.9 mm; accumulated for the month) and a maximum temperature (35.2 ± 3.1 • C) very similar to May (34.4 ± 3.1 • C) and July (36.7 ± 2.9 • C).

Plant Growth Promoting Rhizobacteria Community by Crop and Sampling Time
Among the 547 identified genera in this study, 21 (3.84%) represented genera known to play roles as PGPR. Of these, 11 PGPR genera were found in all crops. Bacillus, Pseudomonas, and Agrobacterium were highly abundant in all crops, with variations in abundance depending on sampling time (Figure 6). Soybean exhibited the most stable PGPR community, where microbial component proportions slightly changed with time ( Figure 6D). Contrary, Pequin pepper and orange tree showed strong changes in abundance of various genera such as Rhizobium, Bacillus, Pseudomonas and Sinorhizobium (Figures 6A-C). There are similarities among crops such as the presence of Agrobacterium and Bacillus in all samples at all sampling times. However, in contrast to global genera distribution, PGPR genera showed differences between young and mature Pequin peppers through sampling times (Figures 6A,B). For instance, the most abundant genus in May is Rhizobium (66.5%) in young Pequin pepper, whereas in mature Pequin pepper is Sinorhizobium (70.8%). Bacillus was found in larger quantity in July for all samples of Pequin pepper and orange, in the latter representing up to 99.1% of the PGPR community ( Figure 6C). This tendency, however, is not observed in soybean samples, where Bacillus has the highest relative abundance in May with 8.21%. For soybean, there was not a high dominance of a single genus at any sampling point, where Agrobacterium, Flavobacterium, and Sinorhizobium presented high abundances ( Figure 6D). The number of samples in which each PGPR genus was identified revealed that a few genera were only observed in specific crops, e.g., Serratia in orange, Azospirillum in young Pequin pepper and Rhodobacter in young and mature Pequin pepper ( Figure 7A). Consequently, we found few PGPR genera specific to one crop and sampling time, such as Serratia, Azospirillum and Alicyclobacillus (Figure 7B).
Since PGPR genera showed differential distributions among sampling times, their correlation with climatological parameters characteristic of drylands (i.e., rainfall, evaporation rate and temperature) was tested using Pearson's correlation (Figure 8). Overall, young Pequin pepper had the highest amount of PGPR genera influenced by climatological parameters, Mycobacterium, Arthrobacter, Mesorhizobium, Rhizobium, Azospirillum, and Burkholderia were significantly affected by an increase in temperature. This trend was similar in mature Pequin pepper for Mycobacterium, Rhizobium and, specific to this crop, in Sinorhizobium. Also, in soybean Arthrobacter presence was affected by high temperatures. For orange trees, Bacillus was significantly correlated to all climatological parameters, affected by high rainfall and increasing its abundance as evaporation rate and temperature increased.

DISCUSSION
For semi-arid ecosystems, it has been shown that bacterial assemblies of adapted plants have a reduced diversity and are strongly influenced by climatological parameters (Taketani et al., 2017). In places with soil aridity and drought, rhizosphere microbiota is driven by plant drought stress alleviation and/or adaptation to environmental conditions (Cherni et al., 2019; see Jansson and Hofmockel, 2019 for a review). For example, Acidobacteria was shown to decrease in abundance as aridity increases, while Chloroflexi and Proteobacteria showed the opposite behavior (Maestre et al., 2015). Here, we also observed this pattern except for Chloroflexi, which we only found in an abundance of around 1% in orange tree roots and which is known for multiple adaptation mechanisms to harsh environments (Lacap et al., 2011). Moreover, Proteobacteria and Firmicutes were reported as the most abundant phyla under drought stress for orange trees (Cherni et al., 2019) but we observed a high abundance of Bacteroidetes, instead of Firmicutes, and Proteobacteria in all three crops. Contrary to the climatic conditions in our study, it was reported that both Firmicutes and Bacteroidetes increased in abundance under heavy rain and low temperatures (Št'ovíček et al., 2017). Also, a higher abundance in Bacteroidetes has been shown to be the result of rich organic root exudates composed of sugar, amino and organic acids (Sánchez-Cañizares et al., 2017).
For pepper plants (Capsicum annuum L.) cultivated in arid ecosystems, it was found that the root endosphere is colonized by a restricted community of bacteria compared to the rhizosphere and adjacent soil (Marasco et al., 2012), but no information is available on whether pepper plant-age affects the microbial communities. Interestingly, no statistical difference was found when analyzing young and mature Pequin pepper root bacteria at the genus level, which is in contrast to some other plants (e.g., rice) where endophytic microbiota communities are characteristic of younger plants when exposed to drought stress (Edwards et al., 2018). In soybean, the rhizosphere community has been reported as a subset of the bulk soil bacterial assembly with some microbial groups being associated with the rhizosphere, such as bacteria belonging to the phyla Chloroflexi, Proteobacteria, Acidobacteria, and Bacteroidetes. These associations are related to plant needs of nutrient uptake and growth enhancement, impacting the diversity of the root microbiome (Mendes et al., 2014).
The dynamics of bacterial assemblies can be related to climatic stressors and plant development (Cao et al., 2019). Bacterial communities in drylands are affected significantly by soil aridity index (evaporation rate/rainfall), total water holding capacity, total organic carbon, total nitrogen and pH . Here, results indicate that significant differences are due to crop specificity and climatic conditions, most pronounced FIGURE 8 | Plant growth promoting rhizobacteria (PGPR) correlated with climatological parameters. Correlation matrix was constructed using Pearson lineal correlation method. "*": significant correlation (P ≤ 0.05). Red and blue color represents positive and negative correlation, respectively. +1 means a perfect positive correlation and −1 means a perfect negative correlation. Rainfall, Evaporation, Tmax and Tmin are the climatic variables measured at the climatological station located at 2 km from the experimental field during 2017.
in May. Regarding climatic conditions, even though July had the harshest environmental circumstances due to the lowest average rainfall and highest temperature and evaporation rate, we hypothesize that significant differences in May root microbiome could have been the result of low or nearly lacking rainfall in the previous 4 months (Supplementary Figure 3). In contrast to July, higher rainfall and a lower evaporation rate occur in September. In leguminous trees, a higher diversity was observed during rainy season compared to the dry season, even though the taxa abundance distribution remained consistent (Taketani et al., 2017). Indeed, microbial response to rain is shown to be dependent not only on its intensity but also on ambient temperature (Št'ovíček et al., 2017). This may indicate that changes in microbial communities may be observed after September and October, the only 2 months where rainfall increases.
With respect to crop specificity, pepper is reported to be more sensitive to water stress, requiring a specific bacteria consortium to reduce water deprivation such as high abundances of Bacillus spp. and Klebsiella spp., known to improve drought stress (Marasco et al., 2012). For orange trees, Bacillus was the most abundant genus at this time, whereas Bacillus has generally been reported as a dominant genus of the endophytic communities in Citrus (Wu and Srivastava, 2012). In the case of soybean, sampling times corresponded to plant developmental stages (i.e., vegetative, flowering and fruiting, and senescent). It has been previously recognized that for soybean the variation in rhizosphere bacterial communities is more influenced by plant growth than by environmental factors (Sugiyama, 2019), which may be reflected by the young soybean microbial composition differing from that of mature or senescent plants, in this study.
In comparison, the orange tree root microbiome remained more static, which is in accordance to previous reports indicating that tree-microbiome interactions differ from annual plants, since trees need decades to grow and strongly depend on the nutrients from the soil (Colin et al., 2017).
Adaptation to the harsh environment of drylands leads to changes in root bacterial assemblies, increasing the presence of PGPR (Dimkpa et al., 2009). When analyzing the correlation of climatological parameters and PGPR genera, Bacillus was found to increase its relative abundance under high temperature and evaporation rates as well as under low rainfall. In this matter, it has been previously reported that Bacillus species enhance resistance to drought stress in various plants (Dimkpa et al., 2009) by inducing systemic resistance and creating a primed physiological state that may cause increased tolerance against abiotic stress (Chakraborty et al., 2006). PGPR strains development for field application is often hampered by the difficulty of bacteria survival under harsh environmental conditions, including temperature (Kaushal and Wani, 2016). Indeed, Arthrobacter, Azospirillum, Burkholderia, Rhizobium, Mesorhizobium, and Sinorhizobium were shown to be affected by high temperature (Huang et al., 2017). Burkholderia spp. have been reported to reduce evapotranspiration under drought stress and improve survival of maize (Huang et al., 2017). Azospirillum, Rhizobium, Mesorhizobium, and Sinorhizobium species have nitrogen fixing activity, form symbiosis with plants and promote their growth (Antoun and Kloepper, 2001;Zarea, 2017). While Azotobacter spp. also can fix nitrogen in a non-symbiotic way, other plant growth promoting traits like phytohormone production and phosphorous solubilization are responsible of promoting plant growth (Antoun and Kloepper, 2001;El_Komy et al., 2020). Arthrobacter spp. have caused mild biotic stress in pepper plants while increasing plant biomass and induced osmotic stress relief in a controlled environment (Sziderics et al., 2007). PGPR genera susceptibility to high temperature is problematic since drylands are extremely prone to present heat stress. Since young Pequin pepper plants showed more PGPR genera affected by high temperatures, our results indicate that they presented PGPR more susceptible to climatic conditions than mature plants. Consistently, it has been reported that the root microbiota community presents high dynamics during the vegetative plant growth phase and stabilizes in composition during the remaining plant life cycle (Edwards et al., 2018); once well-established, the root associated bacteria assembly is structurally more robust (Dombrowski et al., 2017).
In conclusion, pepper root bacterial community dynamics in drylands may be explained by environmental stresses such as long periods of drought, short seasons with rainfall and high temperature and evaporation rates. While these environmental conditions can possibly shape bacterial population in pepper (perennial), for soybean, which is an annual plant, growth stages may be more significant than climatic conditions. For orange trees, the microbial assemblies were more robust. Overall, the presence of Bacillus was correlated to a decrease in rainfall and increase of temperature and evaporation. Our results revealed significant preferential crop-microbe associations, as well as different profiles of PGPR genera affected by rainfall, evaporation rate and temperature. Additionally, several drought stress resistance mediating bacteria may be susceptible to high temperature. Together, these findings suggest that dynamics of root bacterial assemblies are influenced by crop type and semi-arid climatic conditions. Moreover, they provide a better understanding of the suitability of commonly used PGPR genera to promote plant growth on drylands. Further efforts are needed to evaluate potential PGPR candidates for biofertilizers under these conditions.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih. gov/bioproject/PRJNA657454, PRJNA657454.

AUTHOR CONTRIBUTIONS
AMD-G performed bioinformatic and statistical analyses. AS and CS-G conceived the research project. AP, JIF-R, and CS-G collected field samples, constructed the library for Illumina sequencing, and sequenced the samples. MSG-H performed statistical analyses. All authors wrote and revised the manuscript.

ACKNOWLEDGMENTS
We are grateful to Cristina Chuck and Juan Valiente for giving us access to the plants in the field site and to Diego Diaz-Vazquez for his collaboration with geographic localization and data research.