Exploring the root‑associated microbiome of the resurrection plant Myrothamnus flabellifolia

Aims and background The resurrection plant Myrothamnus flabellifolia tolerates complete desic-cation and is a great model for studying how plants cope with extreme drought. Root-associated microbes play a major role in stress tolerance and are an attractive target for enhancing drought tolerance in staple crops. However, how these dynamics play out under the most extreme water limitation remains underex-plored. This study aimed to identify bacterial and fungal communities that tolerate extreme drought stress in the bulk soil, rhizosphere, and endosphere of M. flabellifolia . Methods High-throughput amplicon sequencing was used to characterise the microbial communities associated with M. flabellifolia . Results


Introduction
Many plants are threatened by climate change and fluctuating rainfall patterns (Lizumi and Ramankutty 2016).Although plants have evolved mechanisms to cope with periods of drought, the extreme weather patterns brought on by global change are likely to push some species beyond the limits of their tolerance.However, some unique plants that evolved in habitats with seasonal aridity have developed the ability to tolerate complete desiccation of their vegetative tissues.These species, commonly called resurrection plants, may provide critical insight into how organisms survive the extreme drought that can be used to protect vulnerable species and mitigate drought induced losses (Bewley 1979;Gaff 1989;Hilhorst and Farrant 2018).
Myrothamnus flabellifolia is one of the most iconic resurrection plants and can tolerate desiccation for 9-12 months (Farrant and Kruger 2001).Most of the research on M. flabellifolia has focused on its leaf phytochemical properties, medicinal value, and the physiological and biochemical mechanisms of vegetative desiccation tolerance (Bentley et al. 2019;Erhabor et al. 2020;Marks et al. 2022;Moore et al. 2005Moore et al. , 2006Moore et al. , 2011;;Sherwin et al. 1998).However, the belowground mechanisms associated with desiccation tolerance in M. flabellifolia remain largely unknown (Tebele et al. 2021).Processes happening in the soil and roots play a major role in sensing and responding to drought stress.The occurrence of root elongation, branching, and rhizogenesis are drought-adaptive traits in multiple models and crop species (Kato and Okami 2011).However, little is known about parallel processes in resurrection plants, and they may respond to water deficit in a very different way as a result of their unique adaptations to extreme environmental stresses.
Plant-microbe interactions are another critical aspect of underground drought response mechanisms.Plants can be considered holobionts because they are the cumulative association of multiple living species.In addition, plants can actively select microbial populations that are involved in various beneficial processes to inhabit their cells and surfaces, thereby supporting overall holobiont growth and health.For example, soybean primed with arbuscular mycorrhizal fungi (AMF) and Bradyrhizobium exhibited enhanced shoot growth, root biomass, and higher seed yield compared to non-inoculated plants under drought stress conditions (Sheteiwy et al. 2021b).The roots of nearly all plants support a vast community of microorganisms internally and on their external surfaces, which could aid during environmental stresses.In fact, under drought stress, root microbes not only protect themselves but also minimise drought induced damage to the plant host (Hartman and Tringe).The mechanisms by which microbes improve plant drought tolerance are complex but likely involve reactive oxygen species (ROS) scavenging systems.For example, the upregulation of ROS genes in Burkholderia phytofirmans associated with Solanum Tuberosum roots was observed under drought stress (Sheibani-Tezerji et al. 2015).Additionally, wheat primed with a Paenibacillus polymyxa mutant with an inactivated non-ribosomal peptide synthases greatly improved water use efficiency, water relative content, and increased the activity of antioxidant enzyme activities (Timmusk et al. 2015).The antioxidant enzymes such as superoxide dismutase, glutathione reductase, catalase, and monodehydroascorbate reductase play a significant role in scavenging ROS caused by drought and this defense mechanism enhances drought tolerance in plants.Resurrection plants have robust antioxidant systems, in which there is high activity of ascorbate, catalase, glutathione (GSH), peroxidase (POD), superoxide dismutase (SOD) and that scavenge ROS to combat drought effects (Kranner et al. 2002).How host and microbiome ROS scavenging systems interact to improve drought tolerance is not known, but differences in the ability to recruit and retain beneficial microorganisms can have major impacts on holobiont performance and survival.
Soil microbiomes have a magnitude of functions that are significant to plant health, including nutrient cycling, breakdown of soil organic matter, defense against plant pathogens, and enhancing abiotic stress tolerance (Hartman and Tringe 2019;Yadav et al. 2021).A consortium of microbes is more efficient in enhancing a crop's drought tolerance than one inoculum.A recent study showed that co-inoculation of soybean with Bacillus amyloliquefaciens and AMF alleviated drought stress by increasing the production of sugars such as sucrose, fructose, and stachyose, as well as increased protein content compared to noninoculated plants (Sheteiwy et al. 2021a).Undoubtedly, microbes regulate the physio-biochemical metabolism of the plant under water deficit conditions, which results in metabolic changes and increment of primary metabolites of primed compared to non-inoculated plants in drought stress.
A substantial amount of research has gone into studying these interactions in model and crop species, but the interaction between root and soil microbiomes in resurrection plants remains largely unknown.There is only one study that investigated rhizospheric bacteria of Ramonda serbica and Ramonda nathaliae resurrection plants using high-throughput sequencing (Đokić et al. 2010).Therefore, exploring the interaction of bacteria and fungi associated with the roots of the resurrection plant M. flabellifolia provides novel insight into how plant-microbe interactions function in the most extreme drought events.M. flabellifolia occurs in shallow rocky soil and high-temperature environments with limited rainfall in summer only.It is thus highly likely to harbour diverse microbial communities that possess desiccation tolerance traits.
The soil microbiome is typically divided into multiple compartments based on its proximity to the host plant from bulk soil to rhizosphere to endosphere (Fig. 1).The rhizosphere is the region surrounding plant roots and hosts a wide range of microorganisms, from algae, archaea, bacteria, fungi, to nematodes, protozoa and viruses (Kennedy and de Luna 2005), but the endosphere and bulk soil can also host a diversity of microbes.The most abundant clades in the rhizosphere are plant-growth-promoting rhizobacteria and mycorrhizal fungi that reside within the proximity of the root system and significantly influence plant physiological processes (Tomer et al. 2016).The microbial composition of all three compartments is influenced by multiple factors such as environmental conditions, type of habitat, contribution to plant health, and stress.In some cases, microbial colonisation of the roots leads to endocytosis and it can confer abiotic and biotic stress tolerance within the host plant (Fitzpatrick et al. 2018a;Leborgne-Castel et al. 2010).
To better understand the role of the soil microbiome in extreme drought stress, we quantified the bacterial and fungal communities in the bulk soil, Fig. 1 Soil microbiome compartments are typically categorised into three zones based on their respective proximity to the plant roots.The rhizosphere is the region around the roots into which the root hairs extend.The rhizoplane is a narrow area on the root's surface, and the endosphere is the area inside the roots.Plants actively recruit beneficial bacterial and fun-gal species into these compartments, with increasing specificity from rhizosphere to endosphere.The upper root shows the secretion of metabolites, and the bottom root shows the uptake of nutrients, water, and the process of microbes entering the root endosphere (endocytosis) rhizosphere, and endosphere of M. flabellifolia.The aims of this study were to (1) identify the microbial taxa associated with M. flabellifolia under extreme drought stress, (2) quantify the differential abundance of microbial clades and explore the ecological function of the microbiome based on previous characterisations of the specific microbiota in each compartment.

Materials and methods
Five Myrothamnus flabellifolia plants were collected from Swebeswebe Nature Preserve in the Waterberg, Limpopo Province, South Africa in the autumn of 2021.Individual M. flabellifolia plants and surrounding native soil (~2 kg) were harvested from five discrete locations (within two metres of one another) while the plants were in the state of desiccation (Fig. S1A).These plants were transported to the University of Cape Town for experimental procedures.First, each plant's root system was gently pulled out of the ground to avoid tissue damage.Loose soil attached to the roots was removed by vigorously shaking the plant, and it was classified as bulk soil.Additionally, three replicates of bulk soil were used for soil physicochemical and DNA analysis.Next, a sterilised scissor was used to cut ~5 cm roots with rhizosphere soil attached from the shoot-root joint.Five replicates of the rhizosphere were collected for soil physicochemical analysis by vigorously shaking to dislodge the rhizosphere from the roots.To prepare samples for microbial DNA extraction, five replicates of roots with rhizosphere soil attached were placed in a 50 mL conical centrifuge tube containing 15 mL of autoclaved phosphate-buffered saline (PBS).The rhizosphere soil, which was firmly attached to the root system, was removed according to Edwards et al. (2018).Briefly, tubes were vortexed at the maximum speed for 20 seconds to dislodge the rhizosphere soil from the roots (Fig. S1B).The roots were transferred using a sterile tweezer into a new 50 ml falcon tube with 25 ml PBS.The rhizosphere soil was briefly centrifuged, and the PBS buffer was discarded.Roots were washed with 20 ml of PBS before isolating the endosphere compartment and vortexed at maximum speed for 20-30 seconds.The PBS buffer was discarded, and another 20 ml of PBS was added and vortexed for 20-30 seconds.This step was repeated until there were no soil residues at the bottom of the tube.A fresh 20 ml PBS buffer was added to the root tissues, and they were cleaned by sonicating at 40 Hz for 30 seconds, and this step was repeated twice.Root tissues with no traces of soil were pulverised in liquid nitrogen using a sterile pestle and mortar.All samples were stored at -80 °C for DNA extraction and 4 °C soil physicochemical analysis, respectively.

Physicochemical soil analysis
Soil samples were prepared in microwave-assisted acid digestion (Mars 6 Microwave Digester) using concentrated boric, hydrofluoric, and nitic acid (Al-Harahsheh et al. 2009).Before the elemental analysis, the method was validated using certified reference materials (CRMs), and calibration curves for each analyte were constructed.The targeted elements of, calcium (Ca), copper (Cu), iron (Fe), potassium (K), magnesium, manganese (Mn), sodium (Na), phosphorus (P), and zinc (Zn) were analysed using Varian ES 730 inductively coupled plasma-optical emission spectroscopy (Agilent Technologies, Inc).The organic carbon (OC) and organic matter (OM) measurements/analysis of the soil were performed based on the established Walkley Black protocol (Walker-Black 1934).The soil pH, moisture content, exchangeable acidity, and aluminium ions were also determined according to (Robertson et al. 1999).We used a one-way analysis of variance (ANOVA) with Turkey's pairwise comparison test to test for differences in physicochemical soil analysis between bulk and rhizosphere soil compartments.

DNA extraction, PCR amplification, library preparation, and sequencing
Total DNA was extracted from the bulk soil, rhizosphere, and endosphere samples using DNeasy PowerSoil (QIAGEN, 12888-50) kit following the manufacturer's protocol.The purity and quantity of the extracted DNA were assessed using gel electrophoresis and a Nanodrop ND-2000 UV-Vis spectrophotometer (Thermo Fisher Scientific, USA).The V3-V4 hypervariable region of the 16S rRNA gene from each sample was amplified using the 341F (5'-TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG-3′) and the 505R (5-GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACAG-3′) primers.The mitochondrial peptide nucleic acid (mPNA, 5-GGC AAG TGT TCT TCGGA-3) and plastid peptide nucleic acid (pPNA, 5-GGC TCA ACC CTG GACAG-3) clamps (PNA Bio, Newbury Park, CA, USA) were included to minimise the amplification of host plant DNA in the PCR reaction (Fitzpatrick et al. 2018b).Each 25 μL PCR reaction consisted of the following ingredients: 12.5 μL of 2x KAPA HiFi HotStart mix, 1.25 μL of 0.25 μM mPNA, 1.25 μL of 0.25 μM pPNA, 5 μL of 1 μM forward primer, 5 μL of 1 μM forward reverse, and 2.5 μL microbial genomic DNA.The reaction mixture was placed in the thermocycler for PCR program: denature at 95 °C for 3 min, followed by 25 cycles of 95 °C for 30 sec denaturation, 75 °C for 10 sec PNAs annealing, 55 °C for 30 sec primer annealing, 72 °C for one min elongation, and 72 °C for five min final extension before holding at 4 °C.The ITS1 and ITS2 regions were amplified using the NSA3 (5'-AAA CTC TGT CGT GCT GGG GATA-3′) and NLC2 (5'-GAG CTG CAT TCC CAA ACA ACTC-3′) primers.These Dikaryomycota-specific primers were developed to reduce the amplification of host plant DNA (Martin and Rygiewicz 2005).All reactions, including negative (sterile H 2 O) and positive (ZymoBIOMICS standard) controls, were performed in duplicates.The PCR amplicons of both bacterial and fungal species were purified using Agencourt AMPure XP beads (Beckman Coulter) to remove free primers and primer-dimers molecules.The PCR amplicons were analysed using an Agilent D1000 ScreenTape.Libraries were generated using the Illumina 16S rRNA and fungal metagenomic protocol.For library preparations, the 16S rRNA and ITS1 regions were amplified using a pool of nested primers modified to include the binding site for the sequencing primer and regions for index binding.Similar PCR conditions were used for the index PCR.The ITS1 amplicon varied in length relative to taxa.Library quantification, normalisation, and pooling were performed prior to Illumina sequencing.The 16S V3-V4 libraries were sequenced on an Illumina MiSeq sequencing instrument (Illumina, Inc) using a MiSeq reagent v2 kit (500 cycles) with 2 × 250 bp paired-end reads.The ITS1 libraries were sequenced using a MiSeq reagent nano v2 kit (500 cycles) with 2 × 250 bp paired-end reads.

Bioinformatics analysis of 16S rRNA and ITS data
The resulting fastq files were processed using R for mac OS X GUI version 4.1.2(Urbanek 2022).Raw sequencing paired-end reads were imported to RStudio for pre-processing steps.The quality profile of the reads was inspected using the DADA2 version 1.22.0(Callahan et al. 2016), and the raw sequence reads with poor quality average scores (< 30) were discarded.The bacterial and fungal reads were filtered and trimmed using DADA2 to eliminate primer and adaptor sequences.In addition, error rates for each consensus quality score were evaluated.Samples were processed independently after sharing information between samples using pool 'pseudo'.The denoised forward and reverse reads longer than ten bp were merged into a multiple sequence alignment using DECIPHER package, and amplicon sequence variances (ASVs) were obtained.Chimeric sequences were identified and removed using DADA2.Taxonomic annotation was performed using the most updated and extensive SILVA and UNITE databases for bacteria and fungi, respectively (Nilsson et al. 2019;Quast et al. 2012).
The bacterial and fungal alpha-diversity were calculated based on the pairwise Wilcoxon test metric at the ASV level and was measured using phyloseq, picante, and stat R packages.The beta-diversity of bacterial and fungal communities was assessed by computing Bray-Curtis distance across different microbial taxa of three compartments.The ordination was performed using non-metric multidimensional scaling (NMDS).The differences in microbial composition between compartments were transformed with total sum scaling (Nilsson et al. 2019) from 16S rRNA and ITS data.A non-parametric paired samples Wilcoxon test was performed to evaluate the alphadiversity and taxonomic differences among various compartments.The microbiome dissimilarity distance metrics were also analysed by the multivariate permutation analysis of variance (PERMANOVA), together with the nonparametric statistical Adonis test (vegan package) to establish the homogeneity of dispersion among three compartments.The PERMANOVA test was performed using 999 permutations.A Dirichlet-Multinomial distribution test was used for relative abundance using the HMP package and this test is analogous to a paired sample t-test.However, it evaluates whether taxa frequencies observed in all groups of metagenomic samples are equal.Differentially abundant ASVs were identified using the DESeq2 package and a pseudo count was added to bypass the errors in the logarithm transformation.Taxa on the heatmap were selected from the annotation table fitted to a negative binomial generalised linear model (nbFLM) using analysis of deviance (alpha = 0.05) at the genus level.

Soil physicochemical analysis
There were significant physiochemical and elemental differences between bulk and rhizosphere soil (Table 1).The bulk soil had significantly lower moisture content than the rhizosphere soil (P < 0.01), and exchangeable hydrogen ions were significantly (P < 0.05) lower in the rhizosphere than the bulk soil.Thus the bulk soil was more acidic compared to the rhizosphere.In general, with the exception of Fe and Mn, the rhizosphere had higher concentrations of elements tested (Table 1).Amplicon 16S rRNA and ITS sequencing data A total of 21,871 bacterial ASVs and 1332 fungal ASVs were retained after filtering across all compartments.The Venn diagrams of bacteria and fungi were constructed by creating a presence/absence matrix summarising the total number of ASVs in all samples for each compartment.The ASVs that were present in the bulk soil, rhizosphere, and endosphere were compared.We identified unique and common taxa across the different soil compartments (Fig. 2A and  B).The majority of bacterial and fungal taxa were specific to the different soil compartments.A total of 5584 ASVs of bacterial taxa were shared across all three compartments (Fig. 2A).The bulk soil, rhizosphere, and endosphere each harboured 3215, 5833, and 2899 unique ASVs, respectively.The rhizosphere had the richest diversity of bacterial taxa compared to other compartments and the rhizosphere and bulk soil shared 2751 ASVs.Similarly, with the fungal communities, the rhizosphere had greater unique ASVs (699) compared to other compartments (Fig. 2B).The shared ASVs across all compartments, and between the rhizosphere and endosphere were the same (22).However, the fungal bulk and rhizosphere soil shared higher ASVs (132) relative to other shared ASVs between compartments.Furthermore, the bulk soil and endosphere showed lower overlap in bacterial (275) and fungal (4) ASVs.There was a high overlap of bacterial and fungal taxa between bulk soil and rhizosphere, and the endosphere had the least unique ASVs.

Alpha and Beta-diversity
To understand the community diversity of M. flabellifolia's root microbiome, we characterised the alpha (within the community) and beta (across the communities) diversity of the bulk soil, rhizosphere, and endosphere compartments.Rarefaction curves at the ASV level show the sequencing depth of both bacterial and fungal community diversity (Fig. S2A and S2B).Alpha-diversity, both in terms of species richness using Shannon's index and species evenness using Simpson was measured.In bacterial communities, species richness was significantly higher in bulk and rhizosphere soil compared to the endosphere (Wilcoxon test, P = 0.05) (Fig. 2C).Similarly, fungal richness was significantly higher (P = 0.02) in bulk and rhizosphere soil than in the endosphere (Fig. 2D).
There was no significant difference between bulk and rhizosphere soil for bacteria (P = 0.39) or fungi (P = 0.25).Moreover, the bulk and rhizosphere soil exhibited a high degree of evenness in their bacterial and fungal communities, whereas there was more variability in the endosphere (Fig. 2C-D).The bulk and rhizosphere soil had a greater diversity of bacteria and fungi compared to the endosphere compared.Next, we calculated beta-diversity-a measure of variation in species diversity between compartments (Fig. 3A and B).The PERMANOVA indicated significant variation in bacterial and fungal communities across compartments.The separation between the groups was statistically significant (PERMANOVA, P = 0.01), providing strong evidence for a difference in bacterial composition between soil compartments and endosphere.(Fig. 3A).There was an overlap in bacterial ASVs between the rhizosphere and bulk soil compartments, which was further observed in hierarchical clustering.The fungal community showed substantial variation in species diversity between compartments with a high statistical difference (PERMANOVA, P = 0.001) in variation between three compartments (Fig. 3B).Moreover, hierarchical cluster analysis further confirmed the clustering of compartments with Bray-Curtis dissimilarity distances ranging from 0.1 to 1 distance across compartments for both bacterial (Fig. 3C) and fungal ASVs (Fig. 3D).However, the fungal dissimilarity distance between the compartments was one, indicating that observations were well clustered and had distinct taxa compared to bacterial hierarchy.Taken together, the endosphere microbiome was more distinct compared to the bulk and rhizosphere soil for both bacteria and fungi clustering.

Relative and differential abundance
To better understand the possible functional impact of the root microbiome, relative abundances of different taxa were evaluated (Fig. 4).Not only were different taxa present in different compartments, but the abundance of these bacterial and fungal taxa differed significantly.The top six most abundant bacterial phyla across all compartments were Acidobacteriota, Actinobacteriota, Chloroflexota, Planctomycetota, Pseudomonadota, and WPS-2 (Fig. 4A).The relative abundance of Chloroflexota, Crenarchaeota, Cyanophyta, Gemmatimonadota, RCP2-54, and WPS-2 was more enriched in bulk soil (Fig. 4A).In contrast, the rhizosphere had a high abundance of Bdellovibrionota, Dependentiae, Bacillota, GAL-15, Spirochaetota, and Verrucomicrobiota.Furthermore, we observed an enrichment of Acidobacteriota, Actinobacteriota, Elusimicrobiota, Myxococcota, Patescibacteria, Planctomycetota, and Pseudomonadota in the endosphere.Interestingly, the relative abundance of Chloroflexota decreased from bulk-rhizosphere soil to endosphere, whereas Actinobacteriota increased from bulk-rhizosphere to endosphere.Interestingly there was no significant difference in the relative abundance of bacterial taxa between the rhizosphere and endosphere, but we found that the bacterial microbiome composition was significantly different between the soil compartments (bulk and rhizosphere soil) and endosphere (P = 0.00054).The most dominant fungal phyla found across all compartments were Ascomycota and Basidiomycota (Fig. 4B), but the fungal microbiome differed significantly between the three compartments (P = 0.01).The rhizosphere harboured a more diverse fungal microbiome than the endosphere and bulk soil, and the significantly enriched phyla were Chytridiomycota, Mortierellomycota, Mucoromycota, Rozellomycota, and Zoopagomycota.The endosphere and rhizosphere had a similar composition of the dominant bacterial phyla, while fungal composition across compartments was dominated by Ascomycota.

Environmental factors influencing the microbial composition
To explore the effects of edaphic factors on the microbial composition of bulk and rhizosphere soil of M. flabellifolia, we performed a redundancy analysis (RDA) to relate elemental soil variables to microbiome composition for both bacteria (Fig. 5A and C) and fungi (Fig. 5B and D).Bulk soil samples clustered together more than the rhizosphere samples, indicating a higher variability within the rhizosphere compared to the bulk soil compartment.The soil physicochemical factors pH, moisture, organic carbon and matter, and soil elements positively correlated with a high abundance of bacterial and fungal species in the rhizosphere compartment.Furthermore, the soil elements Ca, K, Mg, Na, P, and Zn were associated with the enrichment of bacterial genera Gordonia, Taonella, and Rhodobacter (Fig. S3A), whereas Ca, Fe, and Mn were significantly associated with the enrichment of fungal genera Spegazzinia, Dothideales and Vishniacozyma in both bulk and rhizosphere soil (Fig. S3B).There was a strong correlation between the rhizosphere microbiome.and soil physicochemical factors.

Discussion
Plant-microbe interactions play an important role in plant growth, health, and productivity.Despite the Fig. 5 Redundancy analysis of soil physicochemical parameters of bulk and rhizosphere soil and their influence on bacterial (A and C) and fungal communities (B and D) Plant Soil importance of plant growth-promoting microbes in plants, particularly in the rhizosphere, microbial species associated with resurrection plants have not been explored (Tebele et al. 2021).Our study provides one of the first high-resolution characterisations of the root microbiome of a resurrection plant.We described the bacterial and fungal microbiomes associated with the resurrection plant M. flabellifolia under air-dry conditions (≤ 10% RWC) and investigated how the microbial community composition varied across compartments from bulk soil to the rhizosphere, to the endosphere.While this provided only a snapshot of communities present under extreme conditions, precluding a full understanding of progressive changes in community dynamics during dehydration, this study does allow speculation on the role of these communities and their contribution to the unique resilience of M. flabellifolia.Our analyses revealed substantial differences in the microbiome of the bulk soil, rhizosphere, and endosphere compartments.Species richness and microbial dynamism were observed in the bulk soil and, to a higher degree, in the rhizosphere compartment.In contrast, the endosphere microbiome was notably reduced in such properties, with only a few select species being present, suggesting that strong filtering occurs at the soil-to-root interface.These findings suggest that regardless of extreme drought conditions, M. flabellifolia harbors a diverse microbial community and may actively regulate microbial assemblages in different compartments.
The most dominant bacterial phyla across all compartments were Acidobacteriota, Actinobacteriota, Chloroflexota, Planctomycetota, Pseudomonadota, and WPS-2, whereas the most dominant fungal phyla were Ascomycota, Basidiomycota, Mortierellomycota and Mucoromycota.These results are parallel to previous studies that also reported the presence of Actinobacteriota, Chloroflexota, and Planctomycetota in the resurrection plant Ramonda as well as in chickpea, sorghum, and other plant hosts (Ahkami et al. 2017;Đokić et al. 2010;Xu et al. 2018;Yadav et al. 2021).However, there are notable differences between the findings on Ramonda spp and the current study.Most importantly, the presence of Flavobacterium, Mucilaginibacter, and Nocardioides in M. flabellifolia but absent in the Ramonda species' rhizosphere, suggests that these microbes could be environmentally or geographically determined (M.flabellifolia occurs in southern Africa whereas Ramonda spp.occurs in Eurasia).In desiccation sensitive species, such as sorghum, Xu et al. (2018) have a low abundance of Acidobacteriota in roots under drought stress.In contrast, our study shows a high enrichment of Acidobacteriota in the endosphere of M. flabellifolia suggesting a possible role for these taxa in coping with extreme environmental stresses.The abundance of these distinctive microbial strains in M. flabellifolia demonstrates that extreme environmental stress may restructure the plant microbiome.Drought stress undoubtedly introduces physiological, metabolic, and genetic responses while also impacting nutrient availability in both microbiome and host plants.Bulk soil samples were clustered together more strongly compared to the rhizosphere compartment (Fig. 5).There was variability in the rhizosphere samples, suggesting that the rhizosphere of each plant hosts a slightly different microbial composition and this could be due to the landscape topography of the natural environment.
We identified differences in soil physicochemical factors that were associated with changing microbiome composition across soil compartments.The bioavailability of micronutrients such as Ca, Cu, Mg, Mn, Na, and Zn in the rhizosphere soil of M. flabellifolia under drought stress seems to impact the relative microbiome abundance and hints at an intriguing plant-microbe-soil interaction.Previous studies have conclusively shown that arbuscular mycorrhizal enhances plant growth by increasing root access to immobile mineral ions and binding to heavy metals and transporting them from roots to shoot tissues (Srinivasagam et al. 2013).These findings showed that rhizospheric fungal species have synergistic interaction with bacterial species in enhancing nutrient uptake and heavy metals resistance in M. flabellifolia.The content of soil moisture and exchangeable H + were significantly correlated with bacterial and fungal communities in the rhizosphere compartment.The exchangeable H + is the predominant factor that influences soil pH (Zheng et al. 2020), therefore, indicating that rhizosphere soil is less acidic relative to bulk soil and this enhances soil nutrients.Taken together, the rhizosphere soil had higher moisture content and micronutrients and it was less acidic compared to bulk soil, which could attract diverse and beneficial microbial communities.
The enrichment of particular bacteria and fungi under extreme drought could indicate a role for Vol.: (0123456789) those microbes in supporting drought tolerance through close linkage with the host plant.We identified numerous taxa that were uniquely enriched in different compartments of M. flabellifolia.The enriched taxa could point towards a fundamental role of microbes in enhancing the drought tolerance of M. flabellifolia.The M. flabellifolia microbiome is extremely complex consisting of over 900 unique bacterial and fungal taxa in each compartment.To summarize this complexity, we highlight selected taxa from each compartment and discuss their possible role in mitigating drought stress.
This study found a significant enrichment of monoderm (Actinobacteriota, Chloroflexota, Bacillota) and diderm (Acidobacteriota, Bacteriodota, and Pseudomonadota) lineages in rhizosphere and endosphere compared to the bulk soil, except for Chloroflexota, which showed high enrichment in bulk soil.The bulk soil compartment of the M. flabellifolia microbiome hosted an abundance of bacterial taxa in Coleofascicus, Flavisolibacter, and Leptolyngbya.These bacterial genera are known to be resistant to drought stress and can improve the physical and biological conditions of rhizosheath soil (Liu et al. 2021;Moreira et al. 2021).The finding of Cyanophyta in bulk soil is significant, as these species are normally restricted to the upper surface of soil crusts (due to their photosynthetic ability).Such species have been shown to play a significant role in soil stabilisation and the addition of organic carbon (Gao et al. 2020b) and their role in the bulk soil could be significant.Interestingly, our study revealed that the bulk soil compartment hosts unique, beneficial bacterial and fungal taxa, and this might be due to the ecological niche of M. flabellifolia.
The significantly enriched fungal taxa in bulk soil were Epicoccum dendrobii, Metarhizium lepidiotae, and Mycenastrum.These enriched fungal genera also have multiple functions related to resilience.For instance, Epicoccum dendrobii is an antifungal agent that enters the plant tissue through stomatal cells and secretes compounds that inhibit anthracnose lesion development caused by pathogens (Bian et al. 2021).Metarhizium spp. is an endophyte that proliferates propagule levels in the rhizosphere (Steinwender et al. 2015).Remarkably, Metarhizium spp was found in the bulk soil of M. flabellifolia, instead of the endosphere compared to previous studies.These findings indicate that both bacteria and fungi have synergic effects in M. flabellifolia, which is in accordance with the study by Sheteiwy et al. (2021b) reported the co-inoculation of soybean seedlings with AMF and Bradyrhizobium improved root length and dry weight of the plant as compared to the controls.The microbial diversity in the bulk soil depends on the plant species and habitat, and the bulk soil compartment of M. flabellifolia hosted a high relative abundance of Chloroflexota and Cyanophyta compared to other compartments.
A plethora of bacterial and fungal taxa were significantly enriched in the rhizosphere compared to other compartments.The bacterial genera include Flavobacterium, Nocardioides, and Streptomyces.These bacterial genera produce indoleacetic acid (IAA) which improves xylem and phloem formation in roots and further increases the formation of lateral and adventitious roots to mitigate drought effects (Khan et al. 2014).This suggests that these species play a vital role during drought stress in M. flabellifolia because a large root system plays a significant role in sourcing water in the deeper soil profiles.A recent study reported the upregulation of the TaIAA15-1A gene is enhanced by exogenous abscisic acid (ABA) treatment and overexpression of TaIAA15-1A in Brachypodium improved drought tolerance through the regulation of antioxidant pathways and ABA signalling, which further reduces ROS and decreased malondialdehyde content (Su et al. 2023).Moreover, high amounts of iron in the rhizosphere soil might be associated with the abundance of Acidicapsa species (Fig. S3A), which produces IAA (Kielak et al. 2016).Interestingly maize primed with Pseudomonas species enhanced the upregulation of protein kinases (SnRKs), which are positive regulators of ABA receptors (SkZ et al. 2018).This may suggest that numerous diverse protein kinases involved in the desiccation tolerance of M. flabellifolia (Ma et al. 2015) might be enhanced by Pseudomonatota to regulate extreme drought stress.Water limitation and nutrient starvation in plants stimulate the transcription of multiple drought-responsive genes that ameliorate drought stress effects.
The enriched fungal taxa in the rhizosphere included Actinomucor, Aspergillus, and Penicillium spp.Actinomucor sp. has been shown to produce abundant quantities of proteases, lipases, and amylases that hydrolyse the components of sorghum seeds (Gao et al. 2020a).Therefore, the presence of hydrolytic-enzymeproducing species under drought conditions in the Vol:.( 1234567890) rhizosphere suggests that they might be involved in the rearrangement of metabolism through the selective degradation of short-lived proteins (Vaseva et al. 2012).A recent study by Sheteiwy et al. (2021b) reported high antioxidant activity and upregulation of catalase and peroxidase gene expression in the soybean seedling under drought stress was enhanced by Bradyrhizobium and mycorrhizal co-inoculant.The presence of trehalose in M. flabellifolia's leaves has been related to fungal species, and trehalose is a desiccation-induced osmoticum in resurrection plants (Farrant 2000;Moore et al. 2011).Although there was variation in the rhizosphere microbiome, this compartment hosted diverse microbial communities that may confer drought tolerance and promote plant growth to the host.
The endosphere of M. flabellifolia was significantly less diverse in microbial taxa than the other two compartments, suggesting that stringent filtering occurs at the soil-root interface.Bacterial genera significantly enriched in the endosphere were Mucilaginibacter, TMa7, and unclassified Pirellulaceae, and fungal taxa were Cladosporium and Penicillium.The Mucilaginibacter and Cladosporium genera are known to produce exopolysaccharides which form biofilms and adhere to root surfaces (Fan et al. 2018;Kielak et al. 2016;Mahapatra and Banerjee 2013).Exopolysaccharides also act as defense mechanisms against drought by sustaining a hydrated microenvironment and reducing water loss (Morcillo and Manzanera 2021).As the rhizosphere soil moisture was significantly greater than bulk soil, it is possible that the presence of exopolysaccharide-producing microbes may slow water loss and alleviate drought stress.These results provide further support for the hypothesis that root-associated microorganisms confer drought tolerance to their host plant.In addition, bacterial (Mucilaginibacter) and fungal genera (Penicillium and Aspergillus) are known to be resistant to heavy metals and remove or absorb Zn 2+ , Cd 2+, and Cu 2+ (Anahid et al. 2011;Fan et al. 2018).This suggests that these species do not only tolerate drought stress but also may play a role in minimising heavy metal toxicity in M. flabellifolia.

Conclusion
Myrothamnus flabellifolia hosts a diverse microbiome that may assist the plant in surviving extreme drought stress.Understanding the microbial composition and the distribution across three compartments and characterising the core taxa in the bulk soil, rhizosphere, and endosphere provides a great opportunity to identify the microbes involved in enhancing M. flabellifolia drought tolerance.Taken together, this study demonstrated that the most dominant bacterial phyla across all compartments were Acidobacteriota, Actinobacteriota, Chloroflexota, Planctomycetota, Pseudomonadota, and WPS-2, and fungal phyla were Ascomycota, Basidiomycota, Mortierellomycota and Mucoromycota.Additionally, the rhizosphere hosted a high level of bacterial and fungal diversity compared to other compartments.Future studies should consider investigating M. flabellifolia under rehydration and dehydration conditions using meta-transcriptomic and metabolomic approaches to explore the expressed microbial genes and responsive metabolites under extreme drought conditions.The identified microbial communities associated with M. flabellifolia roots could have the potential to improve plant resilience and well-being.
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/.

Fig. 2 Fig. 3
Fig.2The Venn diagram shows the number of A) bacterial and B) fungal amplicon sequence variants (ASVs) that are unique and shared across the three soil compartments.Alphadiversity of C) bacterial and D) fungal communities in the three compartments.A Shannon diversity index was used to measure species richness (number of different ASVs in the sample).Simpson diversity index measures the evenness of the microbial community (abundance of ASVs relative to each other) within each compartment

Fig. 4
Fig. 4Relative abundance of A) Bacterial B) Fungal microbial community composition of bulk soil, rhizosphere, and endosphere compartments at phylum level during extreme drought stress.Significant differences (P < 0.01) in the Dirichlet-multinomial distribution of bacterial and fungal phyla between bulk and rhizosphere soil, and bulk soil and endo-

Table 1
Physicochemical factors and abundance of elements (μg/g) in the bulk and rhizosphere soil of M. flabellifolia NB: Mean ± standard error (n = 5).The asterisk on the same column indicates significantly different values (*P < 0.05, **P < 0.01) based on Tukey's pairwise significantly different test.OC Organic carbon, OM organic matter, Exch acidity exchangeable acidity, Exch Al exchangeable aluminium ion, Exch H exchangeable hydrogen, Ca calcium, Cu copper, Fe iron, K potassium, Mn magnesium, manganese, Na sodium, P phosphorus, Zn zinc