Water pH, not soil pH, alters bacterial community structural pattern and nitrogen cycling pathways in date palm (Phoenix dactylifera L.) roots and bulk soil under freshwater irrigation regime

Date palms are widely cultivated in arid agroecosystems, where knowledge of irrigation water effect on their soil and root-associated bacterial communities is limited. Using 16S rDNA metabarcoding, we studied soil and root-associated bacterial communities of date palms growing in United Arab Emirates. Overall, 12.18% of bacterial operational taxonomic units (OTUs) were unique to roots, whereas 41.55% were specific to soil. The bacterial diversity was lower in root and community patterns were distinct between compartments, wherein irrigation water pH was a key structuring factor in both compartments, while salinity (electrical conductivity) was important only in the soil. Co-occurrence network analysis revealed a decrease in complexity in the soil–root continuum, and specific taxa/modules also varied with water pH. We observed a higher abundance of endophyte–saprotroph (Bacillus, Streptomyces and Dongia) dual-role OTUs in both compartments, possibly involved in nutrient mobilisation and plant growth. Based on PICRUSt and trait-based analyses, we showed that these OTUs are putatively involved in the nitrogen cycle (nitrification, denitrification, and assimilatory nitrate reduction). Taken together, we suggest that irrigation water pH, not soil pH, transiently affects belowground bacterial communities and selects bacteria in specific pH ranges, which may be important for nutrient cycling (i.e. nitrogen cycle) in arid agroecosystems.


Introduction
Phoenix dactylifera L., commonly known as date palm, is a perennial tree species of significant economic and cultural importance. This species is renowned for its exceptional ability to thrive under harsh environmental conditions characterised by high temperatures, low rainfall, and high soil salinity and therefore widely cultivated in arid agroecosystems. Arid agroecosystems are distinct from other agroecosystems in terms of nutrient availability (nitrogen (N), carbon, phosphorus, and other micronutrients) due to prevailing harsh environmental conditions. The nutrients present in soil are lost in the form of gases (carbon dioxide, methane, and dinitrogen) or leach into groundwater through percolation due to agricultural practices (Nielsen and Ball, 2015). A major challenge in arid agroecosystems is the availability and quality of irrigation water because the use of unsuitable groundwater (saline and alkaline) can cause changes in soil pH and salinity, which lead to land degradation Batarseh et al., 2021). Soil pH variations can affect soil multifunctionality (potassium, calcium, magnesium, and sodium) (Wei et al., 2022) and crop health. Similarly, salinisation can affect nutrient availability by altering soil properties (flocculation and dispersion), and levels of soil organic matter (SOM) (Zhou et al., 2020), thereby making nutrients unavailable for plants. Irrigation water characteristics are also known to influence edaphic characteristics (soil organic carbon (SOC), electrical conductivity (EC) and pH), which can act as a filter for microbial selection depending on their adaptation capability in belowground compartments.
The microbiota of soil and root compartments perform various roles and work in unison to alleviate stress of their host plant. Alkaline soil pH is one of the important stressors in desert ecosystems (Guo et al., 2021) and soil bacteria from such ecosystems are adapted to alkaline pH through modifications of membrane proteins, accumulation of cationic amino acids (lysine, arginine, and histidine), increased utilisation of H + for ATP synthase activity, and activation of antiporters (Na + /H + and Na + (K + )/H + ), which are important for maintaining cytosolic pH (Padan et al., 2005). Soil bacteria are also known to produce organic acids which react with alkaline mineral ions and form neutral compounds , thereby help mitigating negative effects (i.e., nutrient availability and salinity) of soil alkalinity. Similarly, soil bacteria (i.e. Bacillus and Streptomyces) accumulate osmolytes (starch, choline, and glycine betaine-"organic-in"), salts (potassium-"salt-in"), and antioxidant enzymes (catalase and superoxide dismutase) for survival under salinity-stress conditions (Gregory and Boyd, 2021). In addition to these adaptations, bacterial members (i.e. Bacillus) also produce biofilms through the secretion of extracellular polysaccharides under drought conditions, which helps in aggregating soil particles (Tsotetsi et al., 2022), thereby improving nutrient mobilisation (iron and phosphate) and promoting plant growth in arid agroecosystems. A comprehensive understanding of soil and root-associated bacterial functionalities in arid agroecosystems is crucial for enhancing soil health and augmenting crop productivity in these challenging environments.
Bulk soil serves as the starting material for the selection of bacteria in the roots of Arabidopsis thaliana (Schlaeppi et al., 2014).
The selection of bacteria in root mainly depends on the adaptation of bacteria to stress or resource availability (root exudates and rhizodeposits), which causes variation in bacterial diversity and communities between endophytic and soil (rhizosphere/bulk) communities. The relative abundance of Actinobacteria is generally high in rhizosphere/bulk soil and is followed by Proteobacteria, Cyanobacteria, and Planctomycetes in arid ecosystems (Vaśquez-Dean et al., 2020), whereas the root shows a higher abundance of Acidobacteria, Planctomycetes, Chloroflexi, and Verrucomicrobia (Trivedi et al., 2022). The abundance of bacterial taxa and diversity are known to decrease in a gradient manner in the soil-rhizosphere-endosphere continuum . Higher soil pH can alter root-associated bacterial communities directly or spill-over effects through changes in the plant host (stunted growth, altered root morphology, and function (root exudates)) (Guo et al., 2020). Similarly, soil salinity is also known to influence bacterial communities in barley and cottoncultivated soils . In addition to edaphic characteristics, climate-related factors, such as mean annual temperature and mean annual precipitation, have also been reported as strong predictors of bacterial diversity in both local (Cherif et al., 2015) and global (Vaśquez-Dean et al., 2020) levels in arid ecosystems. However, the effect of irrigation water characteristics and associated edaphic changes on diversity and communities in the root-soil continuum is poorly studied.
Bacterial communities in belowground compartments perform various complex functions such as nutrient cycling, plant growth promotion, and biocontrol, which can be inferred comprehensively using co-occurrence network patterns. Co-occurrence indicates the simultaneous presence of two bacterial taxa in a given system but does not necessarily point to actual physical interactions. Studying microbial interactions through co-occurrence networks provides clues related to microbial diversity and community assembly. Network analysis also provides insights into bacterial traits (mutualism (positive) and competition (negative)), which are important for understanding the stability, resistance, and resilience of bacterial communities. In addition, the environmental impact (i.e. drought, salinity etc) on bacterial functions can be inferred from changes in co-occurrence patterns and keystone OTUs (de Vries et al., 2018). Knowledge of cooccurrence patterns and keystone taxa of arid agroecosystems and their changes in response to perturbations in environmental variables is not well documented, and a thorough investigation is required to better understand ecosystem functioning in arid environments.
The overall well-being of agroecosystems depends on the close interaction between the roots and the surrounding soil and their impact on biogeochemical cycles and ecosystem functions. Soil bacteria enhance SOM and N accumulation under alkaline conditions (Zhang et al., 2021), thereby improving the nutrient status of the soils. Saprotrophic soil-dwelling bacteria degrade SOM through hydrolytic enzymes and release sugars, which are fermented to organic acids and used as precursors for methanogenesis (methane and carbon dioxide) in the microanoxic niches of the soil matrix (Angel et al., 2011). Similarly, nitrogen fixation is carried out by root-associated and free-living bacteria, while ammonia oxidation and denitrification is generally carried out in the soil (Ling et al., 2022). Changes in edaphic characteristics such as soil salinity and moisture can affect nitrification and denitrification in arid soils (Meng et al., 2020). Similarly, pH also affects nitrogen cycle-related pathways, namely, nitrification (Wang et al., 2019b), denitrification (Č uhel et al., 2010), and nitrate reduction (Rafrafi et al., 2015). Therefore, belowground soil and root compartments can act as a source of climate-active gases such as methane, carbon dioxide, and nitrous oxide (Hall et al., 2012) and lead to critical nutrient loss in nutrientlimited arid soils. The effect of irrigation water characteristics on belowground bacterial communities is a crucial aspect of agroecosystems, as it has the potential to affect biogeochemical cycles and ecosystem functioning (Delgado-Baquerizo et al., 2016;Xing et al., 2022).
In this study, we aimed to elucidate the effect of irrigation water characteristics and concomitant changes in soil chemistry on bacterial diversity and communities of root and bulk soil of date palm (Phoenix dactylifera L.), which are commonly cultivated in arid agroecosystems of United Arab Emirates (UAE). We also investigated the co-occurrence network patterns of bacterial communities and prediction-based functional analyses to understand the implications of irrigation water characteristics on putative microbial interactions and nutrient cycling in soil and root compartments of date palms. Overall, unravelling the influence of irrigation water on bacterial diversity, communities, their interactions, and potential functions will further deepen our understanding of belowground bacterial communities, which has enormous impact improving plant survival, plant health and productivity in challenging arid agroecosystems.

Study site and sampling description
We selected seven study sites that are located across the oasis agroecosystem of Al Ain (Abu Dhabi Emirates), UAE. According to the Köppen climate classification, the climate of the study sites was classified as "Bwh", which is indicative of a tropical and subtropical desert climate (Table S1). These sites were characterized by a higher mean annual temperature (MAT) ranging from 25.3 to 27.8°C and a lower mean annual precipitation (MAP) ranging from 75 to 109 mm, which are typical for arid agroecosystems. At each study site, we have collected triplicate soil and root samples at 20-30 cm depth from date palms separated approximately 50 m away from each other in March 2020. The height and age of the date palm trees were similar across sites and they were irrigated with freshwater at a rate of 50-60 liters per tree once every three days. In total, we collected 42 samples (7 farms × 2 compartments × 3 replicates) from seven different sites. Soil and root samples were transported to the laboratory under cooling conditions and stored at −20°C until DNA isolation. In addition, we have also collected irrigation water samples from each site for chemical analyses (pH and EC). Detailed information regarding the chemical analysis (water and soil) has been provided in previous studies Shamim et al., 2022) and briefly included in Table S1. The irrigation water pH (6.96-7.99) and EC (0.32-3.22 dS.m −1 ) of the freshwater varied widely between the sampling sites. Among the study sites, the irrigation water pH was the highest in S3 and S4, while irrigation water EC was the highest in S7 ( Figure S1). Soil pH (6.74-7.74), EC (0.002-0.503 dS.m −1 ), and OM (1.7-6.08%) also fluctuated between sites. We found a weak significant correlation (Spearman's r = −0.45; P < 0.05) between water pH and water EC. Soil pH also showed a significant weak correlation (Spearman's r = −0.46; P < 0.05) with soil OM ( Figure S1).

DNA isolation and V3-V4 amplicon sequencing
Soil DNA isolation was performed according to the manufacturer's instructions using E.Z.N.A soil DNA kit. We added approximately 250 mg of soil sample to a bead-beating tube, along with a lysis buffer (SLSX and DS buffer) and a series of mechanical and chemical treatments to lyse cells and remove contaminants using silica column. Similarly, 1 g of root samples were ground with liquid nitrogen to obtain a ground biomass. This was used for lysis (AP1 buffer) with proteinase K (1.2 mg/mL), followed by binding and purification of DNA as per the DNeasy Plant Mini Kit (Qiagen, USA) protocol. The purity of the DNA was determined using Nanodrop before performing PCR. We used 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVG GGTATCTAATCC) (Herlemann et al., 2011) primers for the amplification of V3-V4 16S rRNA genes. PCR conditions and purification of barcoded amplicon libraries were based on a previous study . The amplicon libraries were sequenced using Illumina MiSeq at the Integrated Microbiome Resource (IMR) lab, Nova Scotia, Canada (https:// www.imr.bio.com) based on standard procedures. The demultiplexed sequence files from this study have been deposited in the Zenodo repository (https://zenodo.org/) (10.5281/ zenodo.60782 92 and 10.5281/zenodo.6371857).

Bioinformatics analyses
We analysed 907,507 raw sequence reads using the Divisive Amplicon Denoising Algorithm 2 v1.12 DADA2 R package (Callahan et al., 2016, accessed on 20 December 2021). We performed primer removal, trimming, error learning, denoising, merging, amplicon sequence variant generation, and removal of chimeric reads using appropriate R packages, as described by Loganathachetti et al. (2022). An additional clustering step was performed at 97% sequence similarity using vsearch v2.15 (Rognes et al., 2016) and a total of 6,279 Operational Taxonomic Units (OTUs) were obtained. We performed a series of steps, namely, clustering, singleton, and chimera removal, using vsearch v2.15, which finally yielded 5,504 OTUs (425,521 reads) after the removal of 1,245 OTUs (12,523 reads). The taxonomic classification of OTUs (5,504 OTUs: 425,521 reads) was carried out by comparing them with the Silva database v138.1 using assignTaxonomy function of the DADA2 R package (Quast et al., 2012). We removed 15 archaeal (45 reads), 11 chloroplast (94,833 reads), 5 mitochondrial (20201 reads), 1 non-bacterial (2 reads), 2 unclassified kingdom level OTUs (36 reads), and 1091 OTUs containing less than five reads (6,852 reads) from the count table. After completing the aforementioned processes, the OTU table consisted of 3487 OTUs (310,404 reads) from 38 samples (read range: 1,147-20,298). We normalised the OTU table based on the sample with the lowest sequence count (1,147), which was used for diversity and community analyses using the rrarefy function of the Vegan R package (Oksanen et al., 2020).
We predicted putative functional pathways of the OTUs using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) for both the soil and root samples. Initially, OTUs were assigned to KEGG orthologues (KO), and then MetaCyc pathway abundance was calculated based on KO. The abundance percentage of putative assimilatory nitrate reduction to ammonium (ANRA) OTUs in the soil and root samples was calculated by comparing the genus information obtained from nirA and narG databases. These databases were prepared after downloading the bacterial taxa information against "nirA" and "narG" search in the Genbank, NCBI (https://www.ncbi.nlm.nih.gov/genbank/, accessed on 29 November 2022) followed by manual curation (Tables S2, S3).

Statistical analyses
We used R v4.0.3 (R Core Development Team, 2020) for all statistical analyses unless stated otherwise. The OTU abundance table was arcsine-transformed to minimise heterogeneity. The environmental variables (water pH, water EC, soil pH, soil EC, and soil OM) were Z-transformed to Z-scores (0-1) to increase homogeneity. We applied the analyses of variance (ANOVA) test and Tukey's HSD post-hoc test to statistically analyse the significant differences in bacterial diversity (richness, diversity, and evenness) between soil and root compartments using the Agricolae R package. We performed Non-Metric Multi-Dimensional Scaling (NMDS) analysis to identify the differences in the community structure between soil and root compartments and Permutational Analysis of Variance (PERMANOVA) for identifying environmental variables important for community structuring using Vegan R package. A forward selection approach was applied for the environmental variables, and the model with the lowest Akaike information criterion was considered the final model. Briefly, we initially tested the significance of individual factors, and the final model was prepared by ranking them based on their R 2 values. NMDS and PERMANOVA analyses were performed separately for the soil and root compartments, respectively. In functional prediction analyses, the significantly different MetaCyc pathways between the soil and root compartments were identified by t-test using Prism 6 software and visualised using the Pheatmap R package.
Co-occurrence network analysis was performed separately for OTUs detected in soil and root using the iNAP workflow (Feng et al., 2022, accessed on 15 September 2022), which was based on a previous method by Deng et al. (2012). First, we performed majority selection which filtered OTUs based on occurrence in at least 80% of samples. The filtered OTU table was subjected to SparrCC correlation and an adjacency matrix was generated (Friedman and Alm, 2012). The OTUs in the adjacency matrix that passed the correlation threshold (<−0.6 to >0.6) significantly (P < 0.05) were used for calculating network centralities (within module: Zi and between module: Pi) using "Greedy modularity optimization" method (Newman, 2004) using the iNAP workflow. We utilized Gephi v0.9 to calculate topological characteristics, including average path length, average degree, and modularity, and to visualise the cooccurrence network patterns. The average degree is defined as the number of connections in relation to the total number of nodes, and the average path length points are the mean distances of all pairs of nodes. Higher modularity indicates denser interactions of nodes within modules and fewer interactions among the nodes of other modules. The nodes in the co-occurrence networks were categorised as module hubs (Zi > 2.5, Pi ≤ 0.62), network hubs (Zi > 2.5, Pi > 0.62), connector hubs (Zi ≤ 2.5, Pi > 0.62), and peripheral OTUs (Zi ≤ 2.5, Pi ≤ 0.62) (Olesen et al., 2007). Except for peripheral OTUs, all types of nodes (connector-, module-, and network hubs) are called as keystone OTUs (Deng et al., 2012).

Bacterial diversity and community structure of soil and root compartments
Of the total OTUs, 41.55% were found only in the soil, 12.18% were found only in the roots, and 46.27% were recorded in both soil and root samples ( Figure S2). Bacterial richness, Shannon diversity index, and evenness were significantly higher in the soil than in the roots ( Figures 1A-C). The rarefaction curves appeared to reach a plateau for both the soil and root compartments, indicating a sufficient sequencing depth ( Figure S2). We found distinct bacterial community structures in the soil and roots in the NMDS ordination space (Figures 1D-F) (R 2 = 0.1635; P = 0.0001). At the individual compartment level, PERMANOVA analyses showed that irrigation water pH was a major factor influencing bacterial community structure in both the soil (R 2 = 0.2157; P < 0.0001) and root (R 2 = 0.12375; P = 0.0018) (Table 1), whereas soil EC was a significant additional factor in the soil affecting bacterial community structure (R 2 = 0.0935; P = 0.0359).

Co-occurrence network patterns and topological characteristics of soil and root compartments
The co-occurrence networks were generated from 436 soil and 129 root associated bacterial OTUs. These networks were used for the calculation of topological characteristics (average degree, network diameter, average path length, modularity, modules, nodes, co-occurrence, and co-exclusion) and identification of keystone OTUs (connector, module and network hubs), which are shown in Figure 2 and Table S4. The complexity of the cooccurrence patterns (higher average degree, higher nodes, lower Effects of compartment, irrigation water pH, irrigation water EC (electrical conductivity dS m −1 ), soil pH, soil EC (electrical conductivity dS m −1 ), and soil OM were tested on overall (soil vs root), soil-and root-associated bacterial communities. P-values were obtained after 9, 999 permutations and boldface indicates statistical significance (P < 0.05). * Indicates significance P < 0.05, while *** indicates P < 0.0005.  Bacterial diversity and community patterns of soil and root compartments. The boxplots show (A) richness, (B) Shannon's diversity index, and (C) evenness tested using ANOVA followed by Tukey's HSD post-hoc test. Statistical inference is highlighted within the plot (P < 0.05). The box spans the interquartile range (IQR; first quartile to the third) with the median indicated by a dark horizontal line, the whiskers show the 1.5xIQR. The bacterial community structure in the NMDS ordination for (D) soil and root compartments, (E) soil compartment, and (F) root compartment. The arrows point in the direction of maximum increase of individual vector variables and had significant effects (P < 0.05) on the NMDS ordination. The S1-S7 points site numbers in (E, F). modularity, and lower average path length) was relatively higher in the soil than in the roots. The soil samples consisted of four modules (SM1-4) while the root contained four modules (RM1-3). Soil modules (SM3 and SM4) were positively correlated to irrigation water pH (P < 0.05). Similarly, specific root modules also showed positive correlation to irrigation water pH (RM2 and RM3) and soil pH (RM2) respectively (P < 0.05) (Figure 2), while RM1 module was negatively correlated to soil OM (P < 0.05). The correlation and the significance of the module-trait (water pH, water EC, soil pH, soil EC, and soil OM) relationships of both soil and root are shown in Figure 2. List of keystone OTUs detected in both compartments are presented in Table S4.
Soil compartment was enriched with Microvirga, Ammoniphilus, and Gaiella, whereas root showed enhancement of Streptomyces, Mycobacterium, Devosia, Allo_Neo_Para_Rhizobium, and Steroidobacter, while Bacillus and Pedomicrobium were highly abundant in both compartments (Table 2; Figure 3C). We performed heatmap clustering to identify taxa that differed significantly between soil and root compartments. Bacterial co-occurrence and co-exclusion network patterns and module-trait relationships of soil and root compartments. Co-occurrence and coexclusion network patterns of (A) soil, (B) root, (C) module-trait analysis of root, (D) topological characteristics, and (E) module-trait analysis of soil. The edges indicating co-occurrence are coloured in green, while co-exclusion is coded in red. The colour of the nodes in the networks and the first column of the module-trait relationship indicate respective modules. The direction of the correlation (scale 1 to −1) and statistical significance (P-value) of module-trait analysis of soil and root modules is given in the tiles. Connector hubs (CH), module hubs (MH) and network hubs (NH) are labelled within the networks.
Of the 426 MetaCyc pathways, 84 were found to be significantly (P < 0.05) different between the root and soil. Biogeochemical pathways related to anaerobic degradation, fermentation, methanogenesis, and nitrifier denitrification-related pathways were enriched in the soil (Figure 4). Similarly, amino acid synthesis-and nitrate reduction-related pathways were enriched in the roots. The percentage of taxa capable of assimilatory nitrate reduction was higher in the root than in the soil ( Figure S3).

Discussion
We find a lower bacterial diversity in root and distinct bacterial communities between compartments (soil and root) of date palms across sampling sites, in which irrigation water pH was an important structure factor for both compartments, while soil EC is an additional factor for soil bacterial communities. The complexity and stability of the network was higher in soil compared to root, and the abundance of OTUs in specific network modules which are known to participate in nitrogen cycle (nitrification, denitrification, and assimilatory nitrate reduction) increases according to water pH in both compartments. We demonstrate that water pH, and not soil pH, plays a major role in altering bacterial communities and modulating specific modules/taxa involved in the nitrogen cycling in arid agroecosystem.

Soil and root-associated bacterial community structure based on irrigation water pH
We attribute lower bacterial diversity (richness, Shannon diversity index and evenness) and selection of only 12.8% of total OTUs in the root to the "root effect", since date palm roots act as "gatekeepers" by selecting specific bacteria from the surrounding soil (Mosqueira et al., 2019) through possible modulations in root exudates and rhizodeposits (Bulgarelli et al., 2012). In this study, we found distinct bacterial communities between soil and root compartments, which were shaped mainly by irrigation water pH. The relative abundance of bacterial taxa in the soil and root compartments of date palm. Relative abundance percentage of bacteria at the (A) phylum, (B) order, and (C) genus levels. * Indicates significantly higher operational taxonomic units (OTUs; P < 0.05) in the marked compartment (root/soil).
Soil pH is known to structure bacterial communities directly by selecting specific bacteria depending on pH preference or indirectly through changes in nutrient availability (Lammel et al., 2018), especially SOC by altering its dissolution . Soil pH has also been reported to structure bacterial communities distinctly according to gradients of agroecosystems (Rousk et al., 2010) and arid (Zeng et al., 2019) ecosystems. It is also known as a global predictor for bacterial diversity (Bahram et al., 2018) and communities. Therefore, we expect an irrigation water pHdependent filtering effect in roots (Lopes et al., 2021;Xiong et al., 2021), which possibly amplifies host selection pressure and transiently select specific OTUs (i.e. Gaiella-soil and Streptomyces-root), thereby distinctly structure the bacterial communities between compartments. In addition, the higher percolation rate of water through soil (aridsol) in the study regions (Abdelfattah and Shahid, 2007) may also play an important role in transient filtering in arid agroecosystems, while the water pH variations may be attributed to the possible increase in carbonate and bicarbonate (Batarseh et al., 2021) levels. In this study, apart from water pH, soil bacterial communities were also structured based on soil EC, wherein salinity negatively affects bacterial growth through plasmolysis and selectively allows colonisation of salinity-tolerant OTUs (i.e. Bacillus), thereby shaping bacterial communities (O'Brien et al., 2019;Zhang et al., 2019).

Soil bacterial co-occurrence patterns are more complex and affected by water chemistry and soil OM
Higher resource availability (i.e. root exudates, rhizodeposits, and organic acids) in roots possibly alter the bacterial co-occurrence pattern from a highly complex (higher average degree, higher nodes, lower modularity, and lower average path length) to a less complex pattern in the soil-root continuum. This, along with higher bacterial diversity in soil, indicates higher network stability compared to roots in this study since complex networks were reported to be more resistant to changes in environmental conditions compared to less-complex simple networks (Santolini and Barabaśi, 2018;Wang et al., 2019a). We observed relatively dense connections within modules than between modules (modularity) in the root than in the soil, which indicates niche partitioning of root-associated bacterial communities of date palm due to selection pressure exerted by the plant host (Qian et al., 2019). In addition, we found a higher percentage of co-occurrence interactions among OTUs in root, pointing a possible enhancement of mutualistic interactions useful for survival and growth under arid conditions (Mosqueira et al., 2019). We detected a higher abundance of putative dual-role root-associated bacteria (Hardoim et al., 2008), Bacillus (biofilm) (Khan et al., 2021), Streptomyces (biocontrol) (Saeed et al., 2017) and Dongia (heat stress), whose modules (RM2 and RM3) correlated positively to irrigation water pH, indicating their versatile role in plant homeostasis and growth promotion. Similarly, we found enrichment of specific OTUs, Hyphomicrobium (denitrification) (Martineau et al., 2015), Romboutsia (acetogenesis), Novibacillus (biocontrol), Gracilibacillus (biocontrol), Pedomicrobium (iron cycle), and Gemmatimonadaceae_unclassified (nitrous oxide reduction) (Oshiki et al., 2022) and their respective modules (SM3 and SM4) to water pH in soil, which perhaps help in regulating soil health and nutrient cycling in arid soil. We also found a module (SM2), which correlates negatively with soil OM concentrations and positively with soil pH, which indicates the preference of bacteria (i.e. Microvirga) for oligotrophic conditions (Zi et al., 2022) and neutral-alkaline pH (Li et al., 2020). Detailed prediction-based functional analysis of OTUs is useful for understanding the irrigation water characteristics associated changes in nutrient cycling of agroecosystems.

Modulation of nitrogen cycle related OTUs and putative pathways in soil and root compartments
The enrichment of aerobic denitrification (nitrifier denitrification)-related MetaCyc pathways along with a higher abundance of Nitrosomonadaceae OTUs indicate that aerobic denitrification is predominantly carried out in the soil. In this study, the enrichment of specific OTUs (i.e. Nitrosomonadaceae) in relation to water pH increase indicate that water pH is an important factor for denitrification in soil instead of soil pH (Č uhel et al., 2010). We also observed a higher percentage of root-associated bacterial OTUs (Bacillus, Streptomyces, Mycobacterium, Actinoplanes and Hyphomicrobium), which are putatively involved in assimilatory nitrate reduction to ammonium (ANRA) (Fahrbach et al., 2008;Xu et al., 2017;Zhang et al., 2021). In addition, the enhancement of ANRA MetaCyc pathways suggest the Heatmap showing putative MetaCyc pathway abundance of bacterial OTUs. The intensity of the colours in the heatmap indicates the relatively higher/lower abundance (scale of −4 to 4) of the MetaCyc pathways (P < 0.05) in the soil and root compartments calculated using PICRUSt2.
importance of root-associated bacteria in preventing nitrogen loss, in which nitrate is converted to ammonium and incorporated into biomass through glutamine biosynthesis (Bernard and Habash, 2009). This is further supported by the presence of OTUs in roots putatively containing nitrite reductase (nirA) and nitrate reductase (narG) at higher abundance.

Conclusion
We show a distinct variation in bacterial diversity and community structure of belowground compartments (root and soil) and changes in water and soil chemistry were the important factors affecting the community. Soil bacterial communities were complex, connected and stable, while root-associated bacterial communities were relatively less complex and stable. Rootassociated bacterial communities show a higher amount of mutualistic interactions compared to soil indicating their critical role in plant homeostasis. The perturbations in bacterial communities and modules to irrigation water pH provide insights into the possible adaptation mechanisms (endophytic-saprotrophic dual-role) of bacterial OTUs (Bacillus, Streptomyces, Mycobacterium, Dongia and Gaiella) in the belowground compartments of arid agroecosystems. Irrigation water pH may modulate specific modules/OTUs (Bacillus, Streptomyces, Mycobacterium, Actinoplanes and Hyphomicrobium) and possibly increase the abundance of aerobic denitrification (soil) and ANRA pathways (root). Taken together, our findings suggest that irrigation water pH, rather than soil pH, alters bacterial diversity and communities in both soil and root compartments by allowing the growth of selected bacteria in specific pH ranges, which may be important for nutrient cycling (i.e. the nitrogen cycle) in arid agroecosystems.

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 below: https://zenodo.org/record/ 6078292#.Y73QrezP2HE, 10.5281/zenodo.6078292; https:// z e n o d o . o r g / r e c o r d / 6 3 7 1 8 5 7 # . Y 7 3 Q w e z P 2 H E , 10.5281/zenodo.6371857.

Author contributions
DSL: Performed bioinformatic, statistical analyses and drafted the article; SM: Conceived the project idea, performed fieldwork and funding acquisition. Both authors contributed to the article and approved the submitted version.

Funding
This research work is supported by Asian University Alliance (AUA) grant, G00003654 and #12R133.