Soil texture is a stronger driver of the maize rhizosphere microbiome and extracellular enzyme activities than soil depth or the presence of root hairs

Different drivers are known to shape rhizosphere microbiome assembly. How soil texture (Texture) and presence or lack of root hairs (Root Hair) of plants affect the rhizosphere microbiome assembly and soil potential extracellular enzyme activities (EEA) at defined rooting depth (Depth) is still a knowledge gap. We investigated effects of these drivers on microbial assembly in rhizosphere and on potential EEA in root-affected soil of maize. Samples were taken from three depths of root hair defective mutant rth3 and wild-type WT maize planted on loam and sand in soil columns after 22 days. Rhizosphere bacterial, archaeal, fungal and cercozoan communities were analysed by sequencing of 16S rRNA gene, ITS and 18S rRNA gene fragments. Soil potential EEA of ß-glucosidase, acid phosphatase and chitinase were estimated using fluorogenic substrates. The bacterial, archaeal and cercozoan alpha- and beta-diversities were significantly and strongly altered by Texture, followed by Depth and Root Hair. Texture and Depth had a small impact on fungal assembly, and only fungal beta-diversity was significantly affected. Significant impacts by Depth and Root Hair on beta-diversity and relative abundances at taxonomic levels of bacteria, archaea, fungi and cercozoa were dependent on Texture. Likewise, the patterns of potential EEA followed the trends of microbial communities, and the potential EEA correlated with the relative abundances of several taxa. Texture was the strongest driver of rhizosphere microbiome and of soil potential EEA, followed by Depth and Root Hair, similarly to findings in maize root architecture and plant gene expression studies.


Introduction
The rhizosphere is defined as soil influenced by plant roots, it represents a zone with complex and dynamic interactions among plant roots, soil and biota (Hartmann et al. 2008(Hartmann et al. , 2009). These interactions are modulated by physicochemical conditions such as pH, organic carbon and altered moisture. A wide variety of organisms inhabits the rhizosphere including bacteria, archaea, fungi, cercozoa, nematodes, algae, viruses and arthropods (Mendes et al. 2013). Rhizosphere microorganisms are attracted from the bulk soil to the plant roots via root exudates that consist of numerous primary and secondary metabolites serving as energy or carbon sources for their growth (Kawasaki et al. 2016;Sasse et al. 2018;Canarini et al. 2019;Cotton et al. 2019). In turn, rhizosphere microorganisms alter plant root exudations through systemic microbe-root signalling mechanisms (Korenblum et al. 2020), and contribute significantly to plant growth and health by enhancing nutrient acquisition, pathogen resistance and stress tolerance (Mendes et al. 2013;Berg et al. 2017;Mohanram and Kumar 2019). Rhizosphere inhabiting bacteria, archaea, fungi and cercozoa -hereafter termed rhizosphere microbiome -contribute to soil quality and health. Extracellular enzyme activities (EEA) produced by members of the soil microbiome accelerate the breakdown of complex organic substances for gaining energy and nutrients from soil (Nannipieri et al. 2002;Kompała-Bąba et al. 2021). Thus, they regulate the availability of nutrients for plants (Nannipieri et al. 2003;Dick 1997). Phosphomonoesterases or acid phosphatases are enzymes that catalyse the hydrolysis of organic monophosphoesters, releasing phosphate for plant and microbial uptake (Nannipieri et al. 2011). Other enzymes like β-glucosidases generally play a role in the last stage of cellulose degradation.
They hydrolyze cellobiose residues, resulting in the release of sugar monomers, which serve as a source of energy for microbial growth and activity (Gil-Sotres et al. 2005;Merino et al. 2016). N-acetyl-βd-glucosaminidases or chitinases are involved in the degradation of chitin (fungal cell wall component) and peptidoglycan (bacterial cell wall component) (Ekenler and Tabatabai 2004). Soil enzyme activities depend on the abundance and diversity of the microbiome and its metabolic activity (Zhang et al. 2017;Kompała-Bąba et al. 2021). However, our understanding of the relationship between different types of soil enzymes, their activities and corresponding microbial abundance and diversity is still limited.
Soils differing in texture have differences in pore size distribution and connectivity and rooting space of plants as observed for the two soils loam and sand used as substrates to grow maize plants . Smaller pore sizes can increase bacterial diversity in soil due to reduced connectivity of microsites (Carson et al. 2010;Hemkemeyer et al. 2018;Seaton et al. 2020), and reduced access of predators (Rutherford and Juma 1992). Some microbial taxa show preference for different particle sizes of the soil, i.e., clay and silt were shown to be more preferred substrates for specific bacteria than sand particles due Vol.: (0123456789) to their higher cation exchange capacity and nutrient content with consequences for the composition of the soil microbiome (Hemkemeyer et al. 2018;Seaton et al. 2020).
Most studies to date investigated rhizosphere microbiome assembly using composite samples of the whole root system of plants (Li et al. 2014;Silva et al. 2017;Gomes et al. 2018;Walters et al. 2018;Renoud et al. 2020;Kusstatscher et al. 2020) despite the great heterogeneity of root system morphology and architecture . Although several studies reported the microbial assembly at different scales of the plant root, e.g., at different root zones along the primary root axis (Rüger et al. 2021), at different root types (DeAngelis et al. 2009;Kawasaki et al. 2016;Pervaiz et al. 2020), root sizes (Zai et al. 2021) and root age (Wei et al. 2021), there is still lack of information regarding the microbial assembly in the rhizosphere along defined soil and rooting depth.
Root hairs were shown to play a significant role in the uptake of water and nutrients, especially in soils with low phosphate (P) contents (Leitner et al. 2010;Klamer et al. 2019). Under limited P and water shortage of soil, root hairs were known to significantly enhance P uptake compared to a root hairless mutant of barley (Ruiz et al. 2020). Especially organic acids are important rhizodeposits of root hairs that help in the acquisition of P (Pantigoso et al. 2020), and may influence the rhizosphere microbiome (Robertson-Albertyn et al. 2017). Further, the presence of root hairs also increases formation of a rhizosheath (Burak et al. 2021). The role of root hairs in microbial colonization was previously examined in barley by Robertson-Albertyn et al. (2017) who reported that root hairless barley mutants had a less complex bacterial community with lower richness and diversity than the wild type in the rhizosphere. However, the impact of root hairs on the microbial assembly remains poorly understood. The two root hairless mutant lines of barley also showed significantly distinct microbial community composition and diversity (Robertson-Albertyn et al. 2017).
Maize (Z. mays) is a plant species with a wealth of data on root traits and rhizosphere microorganisms (reviewed recently by Bonkowski et al. 2021), and thus serves as an excellent model to address specific questions of rhizosphere biology Rüger et al. 2021).
In this study, we aimed to investigate spatial microbial colonization patterns along the roots of maize grown in soil columns ( Fig. 1) in relation to 1) soil texture (Texture: loam and sand), 2) rooting depths (Depth: D1, D2 and D3) and 3) absence or presence of root hairs (Root Hair: root hair defective mutant rth3 and the corresponding wild-type WT). Further, we related soil potential EEA to the microbial abundance, diversity and dominant taxa in order to understand the mechanisms of plant-microbial interactions in the rhizosphere. Based on root gene expression  and root system architecture  analyses at D1, D2 and D3 of the maize plants, we also learned that increasing Fig. 1 22-day-old maize plant cultivated in a soil column (WT in loam, left) under a control environment and maize root system architectures from different treatments obtained from X-ray computer tomography (right). Indicated depths were used for the rhizosphere and root-affected soil sampling. L, loam and S, sand rooting depth means increasing share of young roots (≤ 7 days).
The following hypotheses were assessed: 1) loam as Texture with higher content of organic matter and sorption capacity but with lower porosity ) shows a higher microbial abundance and diversity and potential EEA (V max and K m ) than sand; 2) particularly high differences between microbial community structure are expected between uppermost and lowest rooting depths. This stratification is dependent on porosity and water-holding capacity of the Texture, which means that a bigger difference is expected in sand; and 3) the effect of Root Hair on the microbial community assembly in the rhizosphere of maize, and potential EEA in the root-affected soil are more subtle compared to Texture and Depth.

Soil column experiment and collection of samples
To address our hypotheses, we took samples from a soil column experiment (see Ganther et al. 2021) using a three-factorial design: Texture, Depth and Root Hair.
Briefly, soil columns (25 cm height and 7 cm inner diameter) were either packed with two soils with contrasting textures loam (33.2% sand, 47.7% silt, 19.1% clay, 0.84% C org, 0.084% Nt and pH 6.21) and sand. A mixture of 16.7% loam and 83.3% quartz sand was prepared for the Texture sand (88.6% sand, 8.1% silt, 3.3% clay, 0.14% C org, 0.014% Nt and pH 6.25). Loam and sand (≤ 1 mm mesh size) were fertilized differently to achieve a similar content of plant available nutrients for plant growth based on pre-trials .
Two maize genotypes differing in root hair formation B73 wild-type (WT) and the root hair defective mutant (rth3) were used (Hochholdinger et al. 2008). Planted columns (six replicates per treatment) were grown in a growth chamber for 22 days under controlled conditions (12 h at 22 °C/12 h at 18 °C for day/ night, 65% relative humidity and 350 μmol m −2 s −1 photo-synthetically active radiation). A volumetric water content of 22 and 18% (v/v) was maintained for loam and sand, respectively.
Rhizospheres and root-affected soil were sampled at three depths (D1: 4.5-6.1, D2: 9.0-10.6 and D3: 13.5-15.1 cm from soil surface, Fig. 1). A soil slice of a thickness of approximately 1.6 cm at each depth was taken. Root segments from each depth were gently shaken to remove excessive soil before being briefly submerged, and shaken in a 15 mL centrifuge tube containing sterile 0.3% NaCl. The soil obtained from shaking off the roots is determined as root-affected soil. While the rhizosphere was obtained by centrifugation of soil resuspension at 5000 x g for 30 min at 4 °C. The rhizosphere pellet was stored at −20 °C until the total community (TC-) DNA was extracted.

TC-DNA extraction and amplification of 16S rRNA gene fragments and ITS regions for real-time PCR analysis
The TC-DNA was isolated from about 0.52 ± 0.01 g per sample using the FastDNA Spin Kit and the Geneclean Spin Kit for soil following the manufacturer's instructions (MP Biomedicals, Heidelberg, Germany).

Data analysis of amplicon sequencing
Demultiplexed sequences of the 16S rRNA gene and ITS region fragments were processed as previously described (Ganther et al. 2020;Yim et al. 2020). Briefly, raw sequence reads of the bacterial and archaeal communities were trimmed for primers using cutadapt v.2.3 (Martin 2011). Primer-trimmed sequence reads were error-corrected and merged, and amplicon sequence variants (ASVs, 100% identity) were identified using DADA2 v.1.10.0 (Callahan et al. 2016) within QIIME2 (Bolyen et al. 2018). Each ASV was given a taxonomic annotation using q2-feature classifier classify-sklearn module trained with SILVA SSU rel. 132 database (Quast et al. 2013). For the fungal community, an automated ITS pipeline, so-called PIPITS (PIPITS_PREP, PIP-ITS_FUNITS and PIPITS_PROCESS) with default parameters for the ITS2 region was followed according to Gweon et al. (2015). Read-pairs from Illumina MiSeq sequencers of ITS files were merged into a single file, followed by quality filtered for the next step, generated by the PIPITS_PREP. Then, the PIP-ITS_FUNITS identified the ITS sub-regions using HMMER3 according to Mistry et al. (2013). The PIPITS_PROCESS generated ASVs, and calculated their read abundances as well as the RDP taxonomic assignments using the UNITE fungal ITS reference data set.
Regarding the cercozoan community, sequence reads of the 18S rRNA gene fragments were processed using the customized MOTHUR pipeline v.39.5 (Schloss et al. 2009). Paired-end reads were merged, not allowing any mismatches in primer or barcode sequences, maximum two mismatches and one ambiguity in the target sequence. Assembled sequences with an overlap lower than 200 bp were removed. Merged contigs were demultiplexed, and primers and tag sequences were trimmed. Remaining reads were clustered into operational taxonomic units (OTUs) using VSEARCH (Rognes et al. 2016) according to the abundance-based greedy algorithm (agc) with a similarity threshold of 97%. Clusters represented by less than 350 reads were removed as likely to represent amplification or sequencing noise (Fiore-Donno et al. 2018). OTUs were assigned to taxa using BLAST+ (Camacho et al. 2009) with an e-value of 1e-50 and the PR2 database (Guillou et al. 2013), keeping only the best hit. Sequences were aligned with the template provided by Fiore-Donno et al. (2018), allowing gaps of maximum five nucleotides, and cleaned from chimeras using UCHIME (Edgar et al. 2011) and non-cercozoan sequences.
Sequence contingency tables showing taxonomy identifications and read abundances (ASV or OTU table) were exported for subsequent analyses. For 14 samples, below 4,000 reads were obtained after removing any reads associated to plant materials such as chloroplast and mitochondrial DNA, and these samples were excluded from the analyses. Rarefaction analyses for 16S rRNA gene fragments (bacteria and archaea), ITS regions (fungi) and 18S rRNA gene fragments (cercozoa) were performed, and they showed that the sequences covered the diversity in the analysed samples (Fig. S1).

Soil potential extracellular enzyme activity (potential EEA)
Soil potential EEA (V max and K m ) analyses in the root-affected soil were performed using fluorogenically labelled substrates (Marx et al. 2005;German et al. 2011) based on 4-methylumbelliferone -MUF: 4-MUF-D-glucoside for β-glucosidase (BG), 4-MUF-phosphate for acid phosphatase (AP) and 4-MUF-N-Acetyl-ß-D-glucosaminide for chitinase (NG). Resuspensions of root-affected soil (50 mL) were prepared using low energy sonication (40 J s −1 output energy) for 2 min. Thereafter, 50 μL of soil suspension, 100 μL substrate solution (2.5, 5, 10, 20, 50, 100 μM) and 50 μL of buffer (MES) were transferred into a 96-well microplate (Tian et al. 2020). Fluorescence was measured at 30, 60 and 120 min at 360 nm excitation and 465 nm emission wave-length and at a slit width of 35 nm with a plate reader (TECAN Infinite F200 Pro). Calibration curves were included in every series of enzyme measurements. Enzyme activities were expressed as MUF release in nM g −1 dry soil h −1 . Michaelis-Menten equation was used to determine enzyme kinetic parameters V max , and K m (Eq. 1) Where v is the rate of enzyme-mediated reaction, S is the substrate concentration and K m is an affinity constant equal to the substrate concentration at half of the maximum reaction rate V max .

Statistical analysis
Amplicon sequencing data, ASV or OTU species richness and the Shannon diversity index (Shannon et al. 1948) were evaluated using rarefied reads of the 16S rRNA gene fragments (4121), ITS regions (4993) and relative abundance reads of the 18S rRNA gene fragments. Microbial abundances (copy numbers of 16S rRNA gene and ITS region fragments) and soil potential EEA were fitted into bacterial, archaeal, fungal and cercozoan beta-diversities as environmental variables applying Redundancy Analysis (RDA) using scaling = 1 (vegan package, Oksanen et al. 2020 The 30 ASVs or OTUs with highest relative abundances (bacteria, archaea, fungi and cercozoa) were used to generate heatmaps (presenting the average relative abundance ASVs or OTUs per sample). Further, the Spearman's rank correlation coefficient analysis was performed between the relative abundances of those 30 ASVs or OTUs and the soil potential EEA (BG, AP or NG), applying p value correction based on Benjamini and Hochberg (1995).
The Levene-Test was applied to check the variance homogeneity of the data. The three-way (Texture x Depth x Root Hair) ANOVA was performed to test effects of Texture, Depth and Root Hair on microbial assembly and potential EEA in the rhizosphere and root-affected soil of maize, respectively at p < 0.05. The two-way ANOVA (Depth x Root Hair) was applied per Texture (loam or sand). When the two-way ANOVA indicated significant differences between treatments, then multiple comparisons were followed in loam or sand using Tukey-Test, applying p value correction (at p < 0.05) based on Benjamini and Hochberg (1995).

Results
Microbial abundance and diversity in the rhizosphere of maize affected by Texture, Depth and Root Hair Copy numbers of 16S rRNA gene determined by qPCR were more than three orders of magnitude higher than the copies of ITS fragments (Table S1). Three-way ANOVA of 16S rRNA gene copy numbers revealed a significant effect of Texture and Depth, while the ITS copy numbers did not differ significantly (Tables S1 & S2). The effects of interactions of Texture:Depth and Texture:Root Hair on 16S rRNA gene and ITS fragment copy numbers were significant, but they differed between loam and sand (Table S2).
After sequence processing, the reads clustered into 28,590 bacterial and archaeal, 2694 fungal ASVs and 409 cercozoan OTUs.
Texture significantly affected bacterial and archaeal ASV richness and Shannon indices, with a higher diversity in loam than in sand as shown in the box-plots and by three-way ANOVA (Table 1; Fig. 2A).
Depth showed a significant effect on bacterial and archaeal Shannon indices, dependent on Texture (a significant interaction of Texture:Depth, Table 1). The effect of Root Hair on ASV richness and Shannon indices was dependent on Depth (a significant interaction of Depth:Root Hair, Table 1). Pairwise comparisons revealed that Shannon index in the rth3 rhizosphere was lower at D3 compared to D1 in both loam and sand ( Fig. 2A). This pattern was not observed for the WT rhizosphere. Further, significantly less ASVs and a lower Shannon index were recorded in loam at D3 in rth3 compared to WT rhizosphere ( Fig. 2A).
The RDA showed that bacterial, archaeal, fungal and cercozoan beta-diversities differed between Texture (Fig. 3A, B & C) and Depth (Figs. S2, S3 & S4), and the significance was confirmed by the PERMANOVA ( Table 2). The significant effect of Depth and Root Hair on bacterial, archaeal, fungal and cercozoan beta-diversity was dependent on Texture (significant interactions of Texture:Depth and Texture:Root Hair, Table 2).
Depth explained a higher variation in bacterial, archaeal, fungal and cercozoan beta-diversities in sand than in loam (PERMANOVA, Table 3). For bacteria and archaea, the greatest differences were observed between D1 vs. D3 for both rth3 and WT grown in loam or in sand (Table 4). Significant differences for fungal beta-diversity between D1 vs. D2 and D1 vs. D3 were only detected for rth3 grown in sand (Table 4). Further, for WT grown in loam significant differences of the fungal beta-diversity were found between D1 vs. D3. For cercozoa, the greatest differences in beta-diversity were shown between D2 vs. D3 in the rhizosphere of WT (in loam) and rth3 (in loam or sand) (Table 4).
Overall, the PERMANOVA tests revealed the Texture as the strongest driver, followed by Depth and Root Hair on bacterial, archaeal, fungal and cercozoan assembly (Table 2).
Bacterial and archaeal ASVs were affiliated to 35 phyla but 25 of them had a relative abundance below 1%, and thus they were grouped together and assigned as "Other". The phylum Proteobacteria was most dominant in sand, while in loam the dominant phylum was Firmicutes (Fig. 4A). Acidobacteria and Thaumarchaeota were far more abundant in loam than in sand (Fig. 4A).
Changes in relative abundances of bacteria and archaea at the phylum level along examined depths were also highly modulated by Texture. For instance, the relative abundance of Proteobacteria was significantly lower at D1 than D3 in the rhizosphere of both WT and rth3 in sand, and the opposite was observed in loam (Fig. 4A). The effects of Depth on Thaumarchaeota were only observed in the rhizosphere of maize grown in loam (Fig. 4A). Overall, Root Hair had a subtle effect on the relative abundances at the phylum level (Fig. 4A).
For fungi, all ASVs were affiliated to eight phyla (four of them were grouped together, and renamed as "Other" due to a low relative abundance ≤1%). The fungal phyla Ascomycota and Basidiomycota were dominant in loam or sand (Fig. 4B). Effects of Depth on fungal relative abundances at phylum level were only observed for Mortierellomycota that displayed significantly lower relative abundances at D3 than D1 for WT grown in loam (Fig. 4B). Overall, no effect of Root Hair on fungal relative abundances was detected at the phylum level in the rhizosphere of maize grown in loam or sand (Fig. 4B).
Cercozoan OTUs were affiliated to 19 orders, 11 of them were grouped together and re-named "Other" as each represented less than 1% of the total number of reads. The orders Glissomonadida and Cercomonadida dominated in the rhizosphere of maize grown in loam or sand and at all three depths (Fig. 4C). In loam, no clear effect of Depth or Root Hair was detected, while in sand an increase in the relative abundance of the order Glissomonadida was found with increasing depth, but only in the rth3 rhizosphere. Cryomonadida were mainly detected in sand, and their relative abundances were higher at D1 than at D2 or D3 in WT and rth3 rhizosphere (Fig. 4C).
To increase the level of resolution effects of Texture, Depth and Root Hair on the rhizosphere microbiome were also analysed at the ASV or OTU level.
The 30 bacterial or archaeal ASVs with the highest relative abundance can be depicted from Fig. 5A. ASVs belonging to archaeal family Nitrososphaeraceae (ASV_2 &_10) and the bacterial genera Terrimonas (ASV_19) and Bacillus (ASV_11, _12 & _13) were dominant in loam, while the ASVs affiliated to bacterial genera Dyella (ASV_41 & _45), Massilia (ASV_1, _4 & _14) and Streptomyces (ASV_30) were more abundant in sand (Fig. 5A). The relative abundance of Bacillus (ASV_13) was higher at D3 than D1 for the rth3 grown in loam (at p < 0.05). In sand, the relative abundance of Massilia (ASV_1 & _14) was significantly higher at D3 than D1 for the WT or rth3 rhizosphere. In contrast, the relative abundance of ASV_30 belonging to Streptomyces was significantly higher at D1 than D2 or D3 for the WT or rth3 rhizosphere in sand. Root Hair affected the distribution of several bacterial ASVs in the rhizosphere dependent on Texture and Depth. For instance, ASVs affiliated to bacterial Dyella (ASV_41 & _45) in the rhizosphere of rth3 grown in sand had significantly higher relative abundance at D1 than in the WT rhizosphere (Fig. 5A). It is important to note that not all ASVs affiliated to a certain genus or family showed identical Depth and Root Hair-dependent relative abundances.
The majority of the 30 most abundant fungal ASVs did not show Texture, Depth and Root Hair-dependent differences in their relative abundances (Fig. 6A). But ASV_5 affiliated to the genus Trichophaea and the ASV_7 (unclassified Ascomycota) were dominant in loam and in sand, respectively.
This similar pattern was also revealed for the 30 most abundant cercozoan OTUs (Fig. 7A). Sequences belonging to cercozoan families Allapsidae (OTU_15) and Sandonidae (OTU_8) were dominant in loam, while two other unclassified cercozoan OTUs again of the families Allapsidae (OTU_4) and Sandonidae (OTU_11) were more abundant in sand (Fig. 7A). The relative abundance of OTUs of the families Sandonidae (OTU_11) and Allapsida (OTU_4) was highest at D2 and D3 rhizosphere, respectively, for both WT or rth3 grown in sand. In loam, the OTU assigned to the Allapsidae (OTU_15) showed the lowest relative abundance at D2 for both WT and rth3 rhizosphere.
Soil potential extracellular enzyme activities (EEA) affected by Texture, Depth and Root Hair In general, the soil potential EEA (V max / V) were approximately 2.5-7 times significantly higher in loam than sand (Table 5), and the significant effect of Texture was confirmed by three-way ANOVA ( Table S3). Effects of Depth on the soil potential EEA were Texture-dependent, and a significant interaction of Texture:Depth was indicated (Table S3). The β-glucosidase activity (BG_V) was not affected by Depth in loam or in sand. Acid phosphatase (AP_V) showed lower activity at D3 compared to D1 or D2 for rth3 and WT in loam, but the opposite pattern was observed in sand (Table 5). Chitinase (NG_V) had a higher enzyme activity at D1 than D2 for rth3 (loam) and at D1 than D2 and D3 for WT (loam). Root Hair affected only the BG_V at D1 in sand as WT had a more than two-fold activity higher than rth3.
The soil potential EEA affinity (K m / K) BG_K showed different enzyme affinities at D2 vs. D3 in loam of WT (Table 5). AP_K at D3 in loam and at D2 in sand were significantly higher than other depths of the rth3. NG_K indicated a similar enzyme affinity at all depths. Overall, enzyme kinetics also revealed Texture as the strongest driver, followed by Depth and Root Hair.
Data of soil potential EEA (V max and K m ) were incorporated into the microbial beta-diversity in RDA plots to examine their linkage to the bacteria and archaea (Fig. 3A), fungi (Fig. 3B) and cercozoa (Fig. 3C) of both textures and to increase the resolution per Texture (Figs. S2, S3 & S4). The fitted vectors showed that bacterial and archaeal communities had significantly high activities of BG_V and AP_V in loam, while in sand significantly high enzyme affinity BG_K was observed (Fig. 3A). Spearman's rank correlation coefficients were calculated to discover potential interactions between the relative abundances of the 30 dominant bacterial and archaeal ASVs and the soil potential EEA (Fig. 5B). We found that several bacterial ASVs affiliated to Flavisolibacter (ASV_32) in loam and Paenibacillus (ASV_47) and Massilia (ASV_1, _4 & _14) in sand were significantly and positively correlated to AP_V. The fungal communities had higher AP_V, BG_V, NG_V and NG_K in loam than in sand but differences were not significant (Fig. 3B). Spearman's rank correlation coefficients indicated several ASVs were positively correlated to AP_V, BG_V and NG_V, but they were not significant (Fig. 6B).
The cercozoan communities had significantly higher NG_V in loam than in sand (Fig. 3C). The relative abundance of OTUs affiliated to Allapsidae_OTU004 and Eocercomonas_OTU006 was significantly and positively correlated to AP_V, while the Paracerco-monas_OTU007 was correlated to activity and affinity of NG_V and NG_K, respectively in sand (Fig. 7B).
Overall, the copy numbers of 16S rRNA gene and ITS region fragments were strongly linked with AP_V ( Fig. 3A, B & 3C).
The RDA plot including the soil potential EEA and microbial communities per Texture (Figs. S2, S3 & S4) revealed that the bacterial, archaeal, fungal and cercozoan communities had higher AP_V at D2/D3 than D1 in sand. While in loam, the NG_V was higher at D1 for bacterial, archaeal and cercozoan communities.

Discussion
Effects of Texture on the maize rhizosphere microbiome Substrates with different textures sand and loam used for growing maize in the present work differed per se in pore size distribution, and maize root-induced changes in porosity in the rhizosphere . Due to higher organic matter, silt and clay contents, loam has higher sorption sites than sand causing lower nutrient mobility and leaching substances in loam . The dilution of loam with 83.3% quartz sand (for Texture sand) indeed resulted in a decrease in abundance (copy numbers of 16S rRNA gene fragments) and alpha-diversity (observed ASVs/OTUs and Shannon indices) of bacteria, archaea and cercozoa (Tables 1, S1 & S2; Fig. 2). Higher pore connectivity of sand was recently already reported in different studies to decrease bacterial diversity in soil due to easier migration of microbial cells between pores (Carson et al. 2010;Hemkemeyer et al. 2018;Seaton et al. 2020).
In the present study, changes in cercozoan richness and diversity displayed a similar pattern as those observed for bacteria and archaea (Table 1; Fig. 2). This is likely due to feeding traits of cercozoa. Most cercozoan taxa mainly feed on bacteria , and these bacterivorous cercozoan taxa affect bacteria and archaea in soil due to their specificity in feeding patterns (Kreuzer et al. 2006;Rosenberg et al. 2009;Flues et al. 2017;Henkes et al. 2018). Subsequently, cercozoa are bottom up controlled by bacteria, as the cercozoan community assembly is affected by bacterial defense mechanisms (Jousset 2012), and their growth depends on the availability of their major food source. In addition, the differences in Texture of the two soils might have changed these predator-prey interactions. Loam with a higher content of fine soil pores, might have restricted the access of protists to their bacterial prey (Rutherford and Juma 1992). Table 2 Global PERMANOVA analysis (R 2 value) to reveal effects of Texture, Depth and Root Hair on bacterial and archaeal (16S rRNA gene fragments), fungal (ITS regions) and cercozoan (18S rRNA gene fragments) beta-diversities in the rhizosphere of maize R 2 + * significant differences at p < 0.05 (in bold). n = 6, except for L_rth3_D3, S_WT_D1 and S_rth3_D3 (n = 4) and S_WT_ D3 (n = 5) for bacteria and archaea; S_WT_D2 and S_WT_D3 (n = 5), S_rth3_D3 (n = 4) and S_rth3_D2 (n = 3) for fungi. For cercozoa, n = 6, except for S_WT_D3 (n = 5) and S_rth3_D3 (n = 3). L, loam and S, sand  Table 3 PERMANOVA analysis (R 2 value) to reveal effects of Texture, Depth and Root Hair on bacterial and archaeal (16S rRNA gene fragments), fungal (ITS regions) and cercozoan (18S rRNA gene fragments) beta-diversities in the rhizosphere of maize -in loam or sand R 2 + * significant differences at p < 0.05 (in bold). n = 6, except for L_rth3_D3, S_WT_D1 and S_rth3_D3 (n = 4) and S_WT_D3 (n = 5) for bacteria and archaea; S_WT_D2 and S_WT_D3 (n = 5), S_rth3_D3 (n = 4) and S_rth3_D2 (n = 3) for fungi. For cercozoa, n = 6, except for S_WT_D3 (n = 5) and S_rth3_D3 (n = 3). L, loam and S, sand  Table 4 PERMANOVA analysis (R 2 value) to reveal effects of Depth on bacterial and archaeal (16S rRNA gene fragments), fungal (ITS regions) and cercozoan (18S rRNA gene fragments) beta-diversities in the rhizosphere of maize R 2 + * significant differences at p < 0.05 (in bold). n = 6, except for L_rth3_D3, S_WT_D1 and S_rth3_D3 (n = 4) and S_WT_D3 (n = 5) for bacteria and archaea; S_WT_D2 and S_WT_D3 (n = 5), S_rth3_D3 (n = 4) and S_rth3_D2 (n = 3) for fungi. For cercozoa, n = 6, except for S_WT_D3 (n = 5) and S_rth3_D3 (n = 3). L, loam and S, sand Compared to bacteria, archaea and cercozoa, Texture revealed less effects on the fungal abundance (copy numbers of ITS regions, Tables S1 & S2), alpha- (Table 1) and beta-diversity (Table 2) as previously observed for other soils (Hartmann et al. 2014;Yim et al. 2017;Seaton et al. 2020). Fungal hyphae can easily bridge pore spaces (Ritz and Young 2004), while single cell organisms depend on the connectivity of the water film for dispersal.
Findings using maize roots obtained from the three depth layers of the soil column experiments of the present study, Ganther et al. (2021) found a higher expression level of genes for aquaporins involved in passive water transport, genes related to plant immunity function, e.g., gene related to pathogenesisrelated protein 5 (PR5), chitinase, ethylene-intensive 3 transcription factor and genes related to secondary metabolite production in sand compared to loam. Genes coding for aquaporin (Marulanda et al. 2010), chitinases (Shoresh and Harman 2010), PR5 (Anisimova et al. 2021) and metabolites (Cotton et al. 2019;Murphy et al. 2021) that were differentially expressed between loam and sand, were previously shown to be involved in plant microbe interactions, and might be influenced by the Texture-dependent microbiome composition.

Effects of Depth on the maize rhizosphere microbiome
In an identical experiment, we also learned that share of young roots increased with increasing rooting depth ; Fig. 1). Young or fine roots play active roles in respiration, transport and absorption of water and nutrients, and they also have much more carbon release as exudates for microbial cells compared to the older roots (Nikolova et al. 2020). Thus, they might favor r-strategists (King et al. 2021;Wei et al. 2021) as observed in the present study for Massilia and Bacillus. The assembly of bacterial, archaeal and cercozoan communities in different root zones, i.e., distance to root tips of the primary root of nine-day-old WT maize plants grown in loam, was previously investigated by Rüger et al. (2021). They reported that along the primary root axis, the alpha-diversity was higher at root tips compared to older root regions. In our study, a mixture of roots per depth layer was used for obtaining the rhizosphere microbial pellet, but as the proportion of young roots increased with depth ), we observed a slightly increased alpha-diversity of bacteria and archaea in the rhizosphere of 22-day-old WT maize grown in loam at the D3 compared to D1 ( Fig. 2A), confirming the findings by Rüger et al. (2021).
The soil water content showed a depth gradient due to gravity in the soil column . As plants developed, differences in root length density with depth evolved, and might have introduced additional gradients in nutrient distribution. Therefore, differences in water, O 2 availability and nutrient resources contributed to differences in microbial abundance and diversity at the three depth levels analysed as already reported in earlier studies (Schlüter et al. 2019;Schimel et al. 1999;Schlatter et al. 2018;Li et al. 2020). Further, microbe-microbe interactions in terms of competition and facilitation of niche occupancy (Sasse et al. 2018;Cotton et al. 2019), top-down control by microbial predators, i.e., protists (Bonkowski et al. 2021) -these factors also cause heterogeneity along rooting depth, and thus shape the rhizosphere microbiome.
Larger differences in bacterial and archaeal betadiversity were observed between D1 vs. D3 than between D1 vs. D2 or D2 vs. D3 (Table 4), and similar findings were reported for plant root gene expression analysis using roots from the same experiment and depths ). Depth affected also the structure of acdS gene carrying bacteria in TC-DNAs of the same rhizosphere samples investigated in the present study (Gebauer et al. 2021).
Depth-related effects on fungal gene copy numbers and alpha-diversity were not observed in the present study (at 4.5-15.1 cm depth) likely due to relatively short soil columns being used.

Effects of Root Hair on the maize rhizosphere microbiome
Lack of Root Hair had little effects on microbial assembly compared to Texture and Depth in the present study, but was significant for bacteria, archaea and cercozoa (Table 2). Only in loam Root Hair had also an effect on fungal beta-diversity. This suggests that there might be a link to the findings by Lippold et al. (2021) that total P uptake was significantly lower for rth3 than WT in loam.
The small effect of Root Hair was in line with root gene expression ) and with bacterial community carrying acdS gene (Gebauer et al. 2021) observed for the same experiment.
Root hairs of the WT used in this study were relatively short, of a length of about 0.24 mm Phalempin et al. 2021). This resulted in small differences between the WT and rth3 in the rhizosphere zone and in their role in water and nutrient uptake. This might also explain the rather small effects on rhizosphere microbiome observed. Further, the subtle effects of Root Hair in microbial assembly might be attributed to fertilization of the soils and to well irrigation .
In a previous work, the root hairless rth2 maize accumulated less biomass and P than WT plants, but only under water stress (Klamer et al. 2019), suggesting that testing the genotypes under different P and soil moisture levels could reveal stronger effects of Root Hair on the rhizosphere microbiome. Further, the effects of Root Hair on the rhizosphere microbiome might be larger for other plant species showing longer root hairs, i.e., for barley and growing under water deficiency or drought (Marin et al. 2021).

Relative abundance of microbial taxa affected by Texture, Depth and Root Hair
The dominant phyla in loam were Firmicutes followed by Proteobacteria (Fig. 4) as previously also reported by Ganther et al. (2020), for the same Texture and WT maize, using the same soil column experimental set-up to investigate the effects of X-ray computer tomography on soil bacterial communities. In contrast to loam, the relative abundance of Proteobacteria was dominant in the rhizosphere of WT and rth3 grown in sand. In particular, the genus Massilia was remarkably dominant in the Texture sand that was obtained by mixing loam with quartz sand. Texture-dependent changes in root exudation patterns or the soil physicochemical characteristics might have facilitated the successful rhizosphere assembly of the three ASVs affiliated to Massilia. Acidobacteria and also Thaumarchaeota were mainly detected in loam indicating a preference of these phyla to the soil conditions in loam, e.g., particle size preferences as recently reported by Hemkemeyer et al. (2018). Wei et al. (2021) revealed that most of the enriched taxa at younger roots or at a lower depth were r-strategists, they proliferated, and responded quickly to nutrients available. Indeed, our findings revealed significantly higher relative abundances of ASVs affiliated to the genera Bacillus (in loam) and Massilia (in sand) in the rhizosphere of roots sampled from D3 compared to D1. Isolates affiliated to these genera are also known as r-strategists (Ofek et al. 2012;Wei et al. 2021). Further, strains from both genera were reported to have catalase activity (Yuan et al. 2021), and thus their increased abundance in the maize rhizosphere might reduce H 2 O 2 concentration, and thus reduce H 2 O 2 stress. Root gene expression analyses of maize showed that stress-or defense-related genes such as peroxidases were highly expressed at lower than upper depths .
Slight decreases in oxygen availability, which coincide with the increase in water content with depth, might explain differences in relative abundances of ASVs affiliated to Bacillus. Bacillus isolates can adapt to a lower oxygen availability (Hartmann et al. 2014) that might support the present finding that their relative abundances were enriched at D3 in loam (Fig. 5A).
The fungal relative abundances were not affected by Depth and Root Hair at the phylum level and most of the 30 ASVs observed in the present work (Figs. 4B & 6A). This was likely due to our observation using a short soil column, and was under growth chamber conditions.
ASVs affiliated to the genus Dyella that were detected in significantly higher relative abundances in the rhizosphere of rth3 than WT grown in sand Fig. 4 Effects of Texture, Depth and Root Hair on relative abundances of A) bacterial and archaeal (16S rRNA gene fragments), B) fungal (ITS regions) and C) cercozoan (18S rRNA gene fragments) communities at phylum or order levels in the rhizosphere of maize (> 1%). Letters indicate significant differences of each bacterial, archaeal, fungal and cercozoan phylum or order level in loam or sand, Tukey test applying "BH" p value correction, p < 0.05 (Benjamini and Hochberg 1995). n = 6, except for L_rth3_D3, S_WT_D1 and S_rth3_D3 (n = 4) and S_WT_D3 (n = 5) for bacteria and archaea; S_WT_D2 and S_WT_D3 (n = 5), S_rth3_D3 (n = 4) and S_rth3_D2 (n = 3) for fungi. n = 6 for cercozoa, except for S_WT_D3 (n = 5) and S_rth3_D3 (n = 3). L, loam and S, sand ◂ (Fig. 5A), are known to be involved in nitrogen fixation (Swarnalakshmi et al. 2020). Based on plant roots and rhizosphere microbiome feedback, strains belonging to Dyella might have been recruited to assist the rth3 plants for N-uptake in sand. Interestingly, root genes coding for nitrate transport were shown to be less expressed in the rth3 compared to WT .
Linkage between bacteria, archaea, fungi, cercozoa and soil potential extracellular enzyme activities (EEA) Both plant roots and microbial cells are able to release enzymes into soil (Cabugao et al. 2017) but we assume that the fraction of plant-originated enzymes decreases with the distance from the root, and as we used root-affected soil for the EEA analyses, the detected potential EEA were assumed to be of microbial origin.
In line with the bacterial and archaeal abundance and alpha-diversity, higher potential EEA (V max ) were recorded in loam than in sand (Tables 5 & S1). This is explained by a greater soil organic carbon content in loam Feng et al. 2019;Ren et al. 2018).
The K m /K in soil studies is generally termed as an apparent affinity constant, which is surely affected to a certain extent by soil physical properties. The twoto-three-fold higher K m for acid phosphatase (AP_K) in sand compared to loam (Table 5), indicated the presence of different enzyme systems with a lower affinity resulting in a decline in the overall enzyme function under substrate limitation in loam German et al. 2012). Further, as discussed above, due to the higher sorption capacity of loam, the affinity of AP_K was strongly reduced compared to sand (Table 5).
Variations in K m values at different depths were more pronounced for the AP, followed by NG (chitinase) and BG (ß-glucosidase). They revealed changes in functional traits of microorganisms, i.e., Spearman's rank correlation coefficient (r) to soil potential microbial extracellular enzyme activities (B). Different letters indicate significant differences between treatments of each ASV in loam or sand, Tukey test applying "BH" p value correction, p < 0.05 (Benjamini and Hochberg 1995). BG, ß-glucosidase; AP, acid phosphatase; NG, chitinase; V, V max and K, K m ; Blank, no significant correlation and * significant correlation. n = 6, except for L_rth3_D3, S_WT_D1 and S_rth3_D3 (n = 4) and S_WT_D3 (n = 5). L, loam and S, sand changes in the metabolic activity and microbial community composition (Blagodatskaya et al. 2021).
The enzyme activity and affinity were fitted as environmental variables of the bacterial, archaeal, fungal and cercozoan communities (Fig. 3). Higher activities of BG_V, AP_V and NG_V were shown in loam than in sand demonstrating different abilities of dominant bacteria, archaea, fungi or cercozoa to produce or release the enzymes investigated. The correlation of AP_V with multiple ASVs or OTUs indicated functional redundancy and possible mutualistic interactions in P acquisition and in cellulose degradation within the communities (Banerjee et al. 2016). Positive correlation of AP_V with ASVs affiliated to Massilia and the cercozoa Allapsidae and Eocercomonas in sand possibly indicate different specialized taxa in production of specific enzymes. The bacterial, archaeal and cercozoan communities showed significantly higher AP_V at lower depths in sand A higher enzyme affinity in sand (Fig. 3) to the decomposed compounds indicated different taxa produced distinctly different enzyme systems (Fontaine and Barot 2005;Blagodatskaya et al. 2009;Blagodatskaya and Kuzyakov 2013). The stronger linkage of enzyme parameters with bacteria and archaea than fungi observed in the present study might be due to the lower abundance of fungi but possibly indicated a broader functional diversity within bacterial and archaeal communities that are able to decompose plant and microbial residues as well as to hydrolyze organic P compounds, while fungi invested more resources in P acquisition (Smith et al. 2011;Chiba et al. 2021). Although the potential EAA were analysed in soil taken at a larger distance from root surface, we did find the correlation to the microbiome in the rhizosphere (Figs. 5, 6 & 7).

Fig. 6
Effects of Texture, Depth and Root Hair on 30 fungal (ITS regions) ASVs with highest relative abundances in the rhizosphere of maize (A) and their Spearman's rank correlation coefficient (r) to soil potential microbial extracellular enzyme activities (B). Different letters indicate significant differences between treatments of each ASV in loam or sand, Tukey test applying "BH" p value correction, p < 0.05 (Benjamini and Hochberg 1995). BG, ß-glucosidase; AP, acid phosphatase; NG, chitinase; V, V max and K, K m ; Blank, no significant correlation and * significant correlation. n = 6, except for S_WT_D2 and S_WT_D3 (n = 5), S_rth3_D3 (n = 4) and S_rth3_D2 (n = 3). L, loam and S, sand Effects of Texture, Depth and Root Hair on 30 cercozoan (18S rRNA gene fragments) OTUs with highest relative abundances in the rhizosphere of maize (A) and their Spearman's rank correlation coefficient (r) to soil potential microbial extracellular enzyme activities (B). Different letters indicate significant differences of each OUT in loam or sand, Tukey test applying "BH" p value correction, p < 0.05 (Benjamini and Hochberg 1995). BG, ß-glucosidase; AP, acid phosphatase; NG, chitinase; V, V max and K, K m ; Blank, no significant correlation and * significant correlation. n = 6, except for S_WT_D3 (n = 5) and S_rth3_D3 (n = 3). L, loam and S, sand Table 5 Soil potential microbial extracellular enzyme activities affected by Texture, Depth and Root Hair (nM g −1 h −1 for V/V max and μM for K/ K m ) The data were presented as mean ± standard error of mean. Different letters indicate significant differences between treatments in loam or sand, Tukey test applying "BH" p value correction, p < 0.05 (Benjamini and Hochberg 1995). BG, ß-glucosidase; AP, acid phosphatase; NG, chitinase; V, V max and K, K m . n = 6, except for L_rth3_D3 (BG, AP and NG), L_WT_D3 (AP), S_rth3_D1 (BG and AP) and S_WT_D1/D2 (AP), n = 5; L_WT_D3 (NG), S_rth3_D2 (BG and AP) and S_WT_D3 (NG), n = 4; and S_rth3_D3 (BG) and S_WT_D3 (BG), n = 3 and n.a., not analysed due to lack of sample materials 23.85 ± 9.9 a 278.90 ± 13.8 b 23.45 ± 3.3 ab 9.69 ± 0.5 a 6.69 ± 1.7 a

Conclusion
Our research on the relationship between Texture, Depth and Root Hair has three general implications. First, we have shown that Texture was the strongest driver of the rhizosphere microbial assembly and of potential EEA for both WT and rth3 plants. Second, we have demonstrated that Depth was another driver of the rhizosphere microbiome, suggesting that the abiotic environment may differ between the different layers of the column. Third, the small impact of Root Hair on the rhizosphere microbiome and potential EEA in root-affected soil of maize raises the question of the importance of the root hair length and the low plasticity of root hair defective mutant related to limited P availability as reported by Lippold et al. (2021). Overall, our hypotheses were only partly confirmed which might be due to the resolution level of amplicon sequence analyses, but also due to the experimental design of the column experiment.
Most excitingly, the results of the present study as well as the previously published data on the plant side, e.g., plant gene expression and root system architecture from same and identical column experiments Lippold et al. 2021) showed the same drivers, and highlighted the close linkage between the plant and their rhizosphere microbiome and potential EEA.
Acknowledgements This project was carried out in the framework of the priority program 2089 "Rhizosphere spatiotemporal organization -a key to rhizosphere functions" funded by DFG (Deutsche Forschungsgemeinschaft), project numbers 403637238, 403640293, 403641192, 403664478 and 403635931. Seeds of the maize WT and rth3 were provided by Caroline Marcon and Frank Hochholdinger (University of Bonn). We also would like to acknowledge all participants involved in the sampling at UFZ Halle, Dr. Doreen Babin and Ilse-Marie Jungkurth for reading the MS.
Author contributions MG and MT designed and carried out experiments. ZI, LR, MG and BY collected the samples. BY prepared rhizosphere DNA and amplicon sequencing analyses of 16S rRNA gene fragments and ITS regions. LM prepared amplicon libraries of 16S rRNA gene and ITS region fragments, sequenced and produced ASV Tables. LR performed amplicon-sequencing of cercozoa. ZI analysed soil microbial extracellular enzyme activities. AHB contributed data analyses. BY, ZI, LR, MT, DV, MB, EB and KS wrote the manuscript with contributions from all authors.
Funding Open Access funding enabled and organized by Projekt DEAL. This research was conducted within the research program "Rhizosphere Spatiotemporal Organization -a Key to Rhizosphere Functions" of the German Science Foundation, funded by the Deutsche Forschungsgemeinschaft (project numbers 403637238, 403640293, 403641192, 403664478 and 403635931). AHB was supported by the German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, funded by the Deutsche Forschungsgemeinschaft (DFG-FZT 118, 202548816).
Data availability All raw sequences for bacteria, archaea and fungi were deposited at NCBI within the sequence read archive accession PRJNA677863. Raw sequences of cercozoa are available at the ENA under the accession number PRJEB49274.

Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.