Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 16 September 2025

Sec. Extreme Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1587491

Plant species-specific rhizobiome assembly in the hyper-arid Atacama Desert

J. Fortt,J. Fortt1,2J. Castro-Severyn,J. Castro-Severyn1,2A. ChoqueA. Choque1G. DonosoG. Donoso1A. StollA. Stoll3D. L. JonesD. L. Jones4C. P. SaavedraC. P. Saavedra5B. FuentesB. Fuentes1F. Remonsellez,
F. Remonsellez1,2*
  • 1Departamento de Ingeniería Química y Medio Ambiente, Universidad Católica del Norte, Antofagasta, Chile
  • 2Centro de Investigación Tecnológica del Agua y Sustentabilidad en el Desierto (CEITZASA), Universidad Católica del Norte, Antofagasta, Chile
  • 3Centro de Estudios Avanzados en Zonas Áridas (CEAZA), La Serena, Chile
  • 4School of Environmental and Natural Sciences, Bangor University, Bangor, Gwynedd, United Kingdom
  • 5Laboratorio de Microbiología Molecular, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile

The hyper-arid core of the Atacama Desert represents one of the oldest and driest regions of the world and is characterized by high aridity (precipitation <2 mm y−1), hypersaline soil conditions, extremes in temperature (−5 °C to 50 °C), intense UV irradiation and low organic matter content. Despite this, the Yungay area within the hyper-arid core is capable of supporting vegetation adapted to these extreme environmental conditions, including Distichlis spicata and Suaeda foliosa, which access deep groundwater resources. Little is known, however, about the below-ground microbial community that these plants support. To understand plant-microbe interactions in this environment, we investigated the physicochemical properties in the rhizosphere soils of D. spicata and S. foliosa. In addition, DNA was extracted from the rhizosphere soil and 16S rRNA gene amplicon sequencing performed to describe the taxonomic composition of the bacterial community. Our results revealed significant differences in the physicochemical properties between the rhizosphere soils of the two native plants. D. spicata showed higher Electrical Conductivity (EC), while S. foliosa had elevated ammonium concentrations. The microbial composition also varied between the plant species: Firmicutes (Bacillota), Proteobacteria (Pseudomonadota), Halobacteria, and Actinobacteriota (Actinomycetota) were dominant in both plant rhizosphere samples, but their relative abundance differed. In this context, Halobacteria were highly represented in the soils of D. spicata and Firmicutes (Bacillota) in those from S. foliosa. Furthermore, bacterial genera such as Enterococcus were only present in the S. foliosa rhizosphere, while Natrinema was highly represented in soil from under D. spicata (33.4%) in comparison to S. foliosa (1.5%). The microbial community of D. spicata was strongly influenced by EC, whereas that of S. foliosa correlated more with ammonium levels. These findings advance our understanding of microbial community adaptation in one of Earth’s most extreme environments and provide new insights into plant-microbe interaction in hyper-arid soils.

Introduction

The Atacama Desert, one of the driest and oldest deserts on Earth, extends from 18°S to 27°S and is characterized by extreme environmental conditions including minimal precipitation, high soil salinity, and intense solar radiation (Bull et al., 2016). This desert, particularly its hyper-arid core, experiences an aridity index (AI) of < 0.03 with mean annual rainfall below 2 mm, and temperatures fluctuating drastically between day and night (Azua-Bustos et al., 2017; Hartley et al., 2005; Lacap et al., 2011; Ritter et al., 2018). The Yungay zone, located within this hyper-arid core, presents a Mars-like surface and has remained one of the most extreme environments for millions of years (Clarke, 2006; Hartley et al., 2005; Mckay et al., 2003; Navarro-González et al., 2003). The prevailing environmental conditions create significant abiotic stress, limiting the establishment of vegetation and promoting the accumulation of salts in the soil (e.g., NaCl, KNO3, NaClO4, CaSO4), which stratify due to the absence of leaching (Azua-Bustos et al., 2017; Bull et al., 2016; Voigt et al., 2020; Wang et al., 2017). Despite the poly-extreme conditions prevailing in hyper-arid environments, many studies have reported an abundance of life, implying that biogeochemical recycling of key nutrients (e.g., C, N and P) does occur (Fuentes et al., 2022a; Fuentes et al., 2022b; Knief et al., 2020). In this regard, several studies have reported that plants and associated microbial communities in the Atacama Desert both participate in soil nutrient cycling and the rapid transport of carbon (C) and nitrogen (N), particularly in response to ephemeral inputs of moisture (Jones et al., 2018a; Jones et al., 2018b; Santander et al., 2021; Wu et al., 2021). Interestingly, in the Yungay zone it is possible to find an oasis called the “Yungay Oasis,” which supports a high density of plants and shrubs (Fletcher et al., 2012). Here it is possible to find patches of the halophytes Distichlis spicata and Suaeda foliosa (Calderon et al., 2015). These plants obtain moisture from superficial groundwater, specifically they receive water from an aqueduct that was made for a mineral processing plant that was abandoned over 80 years ago (Chávez et al., 2016; Fletcher et al., 2012; Garrido et al., 2016).

Microbial life in hyper-arid deserts like the Atacama has evolved to withstand extreme environmental pressures, playing crucial roles in nutrient cycling and ecosystem stability (Bull and Asenjo, 2013). Previous studies have shown that even in the most extreme arid soils, microbial communities facilitate the rapid turnover of C and N, especially in response to rare moisture inputs (Ewing et al., 2008; Knief et al., 2020). Several works indicate that is possible to find a great biodiversity, which includes hypolithic cyanobacteria (Warren-Rhodes et al., 2006), non-lichenized fungi (Conley et al., 2006; Fuentes et al., 2020), lichens (Rundel, 1978), cacti (Rundel et al., 1991), and even shrubs and trees (Fletcher et al., 2012), which are unevenly distributed and their development is pressured by biotic and abiotic factors, such as plant disease, biocontrol, topography and hydrology, respectively (Lin et al., 2011; Poulos and Camp, 2010). This microbial diversity, ranging from hypolithic cyanobacteria to fungi, lichens, and bacteria, thrives despite the scarcity of organic matter and water (Wu et al., 2021). In this context, the microbial processes in such environments are critical for sustaining the sparse vegetation and are particularly sensitive to any changes in moisture availability (Pointing and Belnap, 2012). Understanding the microbial diversity and its functional roles in these ecosystems therefore provides unique insights into the resilience and adaptability of life under extreme abiotic stress (Makhalanyane et al., 2015).

In desert ecosystems, the plant–soil-microbial interface (i.e., the rhizosphere) serves as a biological hotspot where plants engage in complex symbiotic relationships with diverse microbial communities (Bonfante and Anca, 2009; Dhanker et al., 2020). These rhizobiomes are essential for nutrient acquisition and stress tolerance, helping native plants from arid zones survive in environments characterized by high salinity and drought (Kim et al., 2013; Niu et al., 2018). In this context, the symbiotic relationships in the rhizosphere can modulate plant metabolism through the production of phytohormones, enhancement of osmoregulatory mechanisms, and increased photosynthetic efficiency (Dhanker et al., 2020; Ruiz-Sánchez et al., 2010). These interactions are particularly relevant in the Atacama Desert, where plants rely on microbial consortia to mitigate the harsh abiotic conditions (Neilson et al., 2012; Stotz et al., 2021). Hence, research on rhizobiomes in extreme environments such as the Atacama Desert, is of critical importance as these microbial communities play a key role in enhancing plant resilience to abiotic stressors (Bashan et al., 2014; Egamberdieva et al., 2017), and understanding these plant-bacteria interactions has the potential to uncover mechanisms of plant resilience and offer strategies for improving agricultural productivity in arid regions (Richardson et al., 2009; Vurukonda et al., 2016).

According to the unique environmental conditions of the Yungay Oasis and the known role of microbial communities in supporting plant life under extreme stress, this study aims to investigate the rhizospheric microbial communities of the main two native plants colonizing the area, namely Suaeda foliosa and Distichlis spicata. Specifically, we seek to characterize the physicochemical parameters of their rhizospheric soils and correlate these with the taxonomic composition of the bacterial communities. By analyzing the relationship between abiotic soil conditions and microbial diversity, we aim to determine how the hyper-arid soil conditions of the Yungay Oasis in the Atacama Desert modulate rhizosphere microbiomes. This study will contribute to a deeper understanding of the plant-bacteria interactions that enable survival in one of the most extreme environments on Earth, with potential applications for enhancing crop resilience in the north of Chile and similar arid regions.

Methodology

Sampling site description

The study site was located within the Yungay area, in the hyper-arid core of the Atacama Desert in the Antofagasta region of Chile (24°03′S, 69°49′W; 948 meters above sea level; Jones et al., 2023). This study was focused on the rhizosphere soil of the native plants, namely the drought tolerant shrub Suaeda foliosa Moq. (de Froidmont, 2018) and the salt grass Distichlis spicata (L.) Greene (Pelliza et al., 2005). We selected three sampling station zones covering approximately 2,500 m2. Rhizosphere soil was collected in 50 mL falcon tubes and sterile bags from five individuals per species and three control samples corresponding to bare soil (Figure 1A). Both kinds of samples (rhizosphere and control soil) were collected from 0–20 cm depth. A total of 33 soil samples were taken specifically from the Yungay Oasis (Sampling Stations 1, 2 and 3), 30 samples from the plant rhizosphere and 3 control (unvegetated bare soil) samples. A part of each sample was stored in fresh and dry conditions at room temperature (20–25 °C) for physicochemical analysis, while the other part was stored at −80 °C for followed DNA extraction.

Figure 1
Map A shows sampling stations SST1, SST2, and SST3 in a desert area, with a control point marked. Inset map locates Antofagasta, Chile. Diagram B details rhizospheric soil study methodology: plant samples from D. spicata and S. foliosa, DNA extraction, sequencing, physiochemical analysis, and statistical validation using software and charts.

Figure 1. (A) Study area location and the three sampling stations (SST1-SST3). (B) Photographs of the plant species from the Yungay Oasis and an overview of the main methodological processes used in the study.

Soil physicochemical characterization

All samples were taken and processed as shown in Figure 1B and taken to the laboratory to conduct physicochemical, microbiological and molecular procedures. The moisture content (MC) of the soil sample was measured by drying soil at 105 °C in the oven and calculating the loss of weight after drying. Furthermore, the EC and pH were measured using standard electrodes in a soil-distilled water suspension with a solution-to-soil ratio of 1:2.5 (w/v). The soils were extracted with 0.5 M K2SO4 (1:5 w/v), centrifuged (24,000 g, 5 min), and the supernatant recovered for NO3 analysis using the vanadate colorimetric procedure of Miranda et al. (2001). The same extracts were used to determine NH4+ content using the salicylate-based colorimetric procedure of Mulvaney (1996). Dissolved organic C (DOC) and N (DON) in the extracts were measured using an Multi N/C 2100S analyzer (Analytik Jena GmbH, Jena, Germany). CaCO3 content was determined with a FOGL Benchtop Soil Calcimeter (BD Inventions, Thessaloniki, Greece). Soluble salts (Na, K, Ca) were determined in soil-distilled water extracts using a Sherwood 410 flame photometer (SciMed Ltd., Stockport, United Kingdom). The rhizosphere soil elemental composition was determined using a non-destructive S2 Picofox TXRF spectrometer (Bruker Inc., Billerica, MA) using Se as an internal standard (Towett et al., 2013). Plant-available P in soil was determined using both the Olsen P method (0.5 M NaHCO3 pH 8.5, 1:5 w/v extract) according to FAO (2021) and the acetic acid method (0.5 M CH3COOH pH 2.5, 1:5 w/v extract) according to McCray et al. (2012).

DNA extraction and sequencing

For the DNA extraction, 20 g of soil were washed and resuspended by vortexing with sterile water in a ratio of 1:2 (soil:water). Then, we filtered the supernatant using Millipore Sterivex™ Pressure Filter 0.22 μm and used the filter paper like sample for the DNA extraction. Total DNA was extracted from the samples using E. Z. N. A.® Soil DNA Kit (Omega Bio-tek, Inc.) according to the manufacturer instructions. DNA integrity, quality, and quantity were verified using 1% agarose gel electrophoresis, OD260/280 ratio, and fluorescence using a QubitTM 4.0 Fluorometer along with Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific). For metabarcoding, the total DNA samples were sent to the Environmental Sample Preparation and Sequencing Facility at the Argonne National Laboratory (Illinois, United States) for the PCR amplification of the bacterial 16S rRNA gene V4 hypervariable region [~250 bp; 515F and 806R primers; (Caporaso et al., 2011)], and the construction of paired-end (250 bp) libraries and sequencing on the MiSeq platform (Illumina Inc., San Diego, CA).

Metabarcoding data processing and analysis

Analyses were conducted in R v3.5.2 and RStudio v1.1.463 (Team, 2023; R Core Team, 2020) using the R package DADA2 v1.16.0 pipeline following https://benjjneb.github.io/dada2/tutorial.html tutorial (Callahan et al., 2016) to infer amplicon sequence variants (ASVs) present in each sample. Briefly, raw reads underwent quality control filtering with the following parameters (no ambiguous bases, Ns = 0; minimum length of 150 bp, maximum expected error rate of 2). Filtered reads were then dereplicated, denoised, and merged into paired-end sequences. After the ASV table was built, chimeras were removed and taxonomic assignment was performed using DADA2’s Ribosomal Database Project’s (RDP) naive Bayesian classifier (Wang et al., 2007) against the Silva v138 database (Quast et al., 2013). The variance-stabilizing transformation was applied to normalize data using the DESeq2 v1.28.1 (Love et al., 2014). A multi-sequence alignment was created using the R package DECIPHER v2.16.1 (Wright, 2016) to infer phylogeny with FastTree v2.1.10 (Price et al., 2009). A phyloseq-object containing the ASVs, taxonomy assignments, phylogenetic tree, and sample metadata was created using the R package Phyloseq v1.32.0 (McMurdie and Holmes, 2013). Alpha diversity indices (Shannon, Simpson, and Chao1) were calculated from this object, while beta diversity (PCoA using Bray Curtis distances) was calculated using the R package ampvis2 v2.4.5 (Andersen et al., 2018).

Statistical analysis

Differences, assumptions of normality and homogeneity of variance in soil characteristics were determined by a One-way ANOVA test from R project software. Normality of residuals was evaluated using the Shapiro–Wilk test, and homogeneity of variances was assessed using Levene’s test. When assumptions were met, ANOVA followed by a post-hoc Tukey’s HSD test was applied to identify differences among sample groups. Tukey post-hoc comparison was used when significant differences (p < 0.05) were found in each analysis. In addition, a non-parametric Kruskal–Wallis test was conducted to confirm results under a distribution-free approach, and pairwise comparisons were used to identify significant groups. To identify the most relevant physicochemical parameters, an Analysis of Principal Component (PCA) and Redundancy analysis (RDA) was performed using the ampvis2 R package (Andersen et al., 2018).

Results

Rhizosphere soil characterization

A summary of the physical and chemical properties of rhizosphere soil collected from under the two plants in the hyper-arid core of the Atacama Desert are shown in Supplementary Table 1. The measured pH in our samples was between 6.01 to 7.80, being generally lower under the plants in comparison to the unvegetated control soil. The EC was high in all samples and varied from 18.0 to 92.2 dS/m with values being significatively higher in the D. spicata rhizosphere soil in comparison to S. foliosa and control soil (Figure 2; Supplementary Table 1). Overall, the concentration of NH4+ was very low in comparison to NO3, however, NH4+ concentrations were significatively higher in the S. foliosa rhizosphere in comparison to the control and D. spicata soils (Figure 2). On average, the S. foliosa samples soils presented higher NH4+ contents, and the D. spicata rhizosphere soils presented a higher NO3 content, although some differences were apparent depending on the sampling station (Supplementary Table 1; Figure 2). Both the Olsen-P and acetic acid-P method recovered very low levels of available P, although more plant-available P was determined in the S. foliosa soil using the Olsen-P method. There were no significant differences in plant-available P determined with acetic acid or the CaCO3 and moisture content of the soils across the vegetated and control areas (Supplementary Table 1).

Figure 2
Stacked bar chart showing the relative abundance of different bacterial phyla across different treatments: Control, D. spicata, and S. foliosa. Each color represents a different phylum. Firmicutes and Proteobacteria have significant proportions in all samples. The chart includes multiple other phyla such as Bacteroidota, Actinobacteriota, and Chloroflexi, among others, each with varying abundances across samples and treatments.

Figure 2. Comparison of physicochemical parameters from S. foliosa, D. spicata and control unvegetated soil rhizosphere among the sampling stations (SST1-SST3). (A) EC, (B) Ammonium content, (C) Dissolved organic carbon and (D) Dissolved organic nitrogen. The soil samples corresponding to Sampling Station indicated in Figure 1A. Values represent means ± SEM (n = 5) and statistical differences were assessed using one-way ANOVA followed by Tukey HSD post-hoc test. (A) F = 11.64, p > 0.001. (B) F = 4.16, p = 0.013. (C) F = 4.50, p = 0.003. (D) F = 4.66, p = 0.002.

Microbial community composition

Analysis of rhizosphere soil samples revealed distinct bacterial community compositions associated with D. spicata and S. foliosa in the Yungay Oasis. The 16S rRNA amplicon sequencing identified Firmicutes (Bacillota), Proteobaceria (Pseudomonadota), Halobacteria and Actinobacteriota as the predominant phyla across all soil samples (Figure 3). Halobacteria and Proteobacteria (Pseudomonadota) were the most abundant phyla in the control soils and D. spicata rhizosphere, while Firmicutes (Bacillota) and Proteobacteria (Pseudomonadota) were most abundant in the S. foliosa rhizosphere. Of note, was the greater abundance of Halobacteria in the control soil and D. spicata rhizosphere compared to the S. foliosa rhizosphere, where Firmicutes (Bacillota) dominated, although this group was also present in D. spicata rhizosphere samples (Figure 3).

Figure 3
Heatmap depicting relative abundance of various microbial taxa across different conditions: Control, *D. spicata*, and *S. foliosa*. Abundance ranges from 0 to over 40 percent, with color coding from yellow to dark blue indicating increasing abundance. The table includes genera such as *Halobactero*, *Firmicutes*, *Proteobacteria*, among others, with varying abundance values across different sample sites (SST1, SST2, SST3). A color gradient legend on the right indicates abundance percentages.

Figure 3. Bacterial community composition at phylum level in the hyper-arid core of the Atacama Desert. Samples include rhizosphere soil from D. spicata and S. foliosa collected at three sampling stations (SST1–SST3), alongside those collected from control unvegetated soil. Bars represent the mean relative abundance from five replicates per plant species rhizosphere and three replicates for the unvegetated control soils.

The 40 most abundant bacterial genera were identified across the three sampling stations to characterize differences between the rhizobiomes of both plant species (Figure 4). In accordance with the results shown above, most of these taxa belong to the Proteobacteria (Pseudomonadota), Firmicutes (Bacillota), Halobacteria and Actinobacteria phyla. However, several abundant genera showed substantial variation across sampling stations, with some being highly abundant at certain sampling stations while rare or absent at others, namely: Natrinema, Buchnera, Enterococcus and Salibacterium, supporting a high variability in these rhizosphere soil communities. In this context, the variability could be attributed to the extreme conditions of the environment, where microhabitats may differ significantly even within a small area such as native plant rhizosphere. Hence, we hypothesized the native plants in these extreme environments might exhibit low specificity in their microbial associations, potentially forming symbiotic relationships with a broad range of microorganisms that have basic functional requirements to colonize these environments. The results indicate that between D. spicata and S. foliosa rhizosphere, major differences existed in the relative abundance of the genus Natrinema (33.4% vs. 1.5% in SST3), Buchnera (3% vs. 42.3% in SST3) and Enterococcus (0% in SST1, SST2 and SST3 vs. 23.2, 7.9 and 10.9% for SST1, SST2 and SST3, respectively). Key genera in the control soil and D. spicata rhizosphere included: Natrinema and Haloterrigena with 31.8 and 33.3% for Natrinema and 22.3 and 8.8% for Haloterrigena in the control and SST3 samples, respectively. Furthermore, in the S. foliosa rhizosphere, we found a dominance of the Firmicutes (Bacillota) phylum, specifically for Buchnera (42.3%) genus in SST3 and Enterococcus genus with 23.2, 7.9 and 10.9% for SST1, SST2 and SST3, respectively. In addition, other Firmicutes (Bacillota) genera such as Paenisporosarcina (7.8%) and Paenibacillus (7.5%) were also abundant at SST1. SST 3 tended to be more heterogeneous in terms of bacterial taxa abundance and, the control soil was similar to D. spicata rhizosphere in taxonomic composition with the presence of genus Natrinema, Haloterrigena, Scopulibacillus and Nocardioles (Figure 4).

Figure 4
Principal Coordinates Analysis (PCoA) plot displaying soil samples. The x-axis represents PCo1 with 15.5% variation, and the y-axis shows PCo2 with 12.8% variation. Green and blue polygons represent different treatments,

Figure 4. Taxonomic composition and genus abundance of the rhizosphere bacterial community of two native plant species (D. spicata, S. foliosa) sampled in three sampling stations (SST1-SST3) in the hyper-arid core of the Atacama Desert. Color gradient represents the relative abundance of the top 40 bacterial genera. Figure represents the mean relative abundance from five replicates per plant species rhizosphere and three replicates for the unvegetated control soils.

Beta diversity analysis revealed distinct bacterial community compositions in the rhizospheres of the two plant species (Figure 5). The bacterial communities from control soils clustered with those of D. spicata, while communities associated with S. foliosa formed a separate group (Figure 5).

Figure 5
Biplot of redundancy analysis showing two axes: RDA1 (10.7%/70.2%) and RDA2 (4.5%/29.8%). Green and blue points represent soil samples of *D. spicata* and *S. foliosa*. Environmental variables include pH, temperature, moisture, NH4, EC, and BR, with arrows indicating their direction and influence.

Figure 5. Principal component analysis (PCA) of the rhizosphere bacterial community composition of soil samples collected from under two native plant species (D. spicata, S. foliosa) alongside control unvegetated soil in the hyper-arid core of the Atacama Desert. The percentage of variance explained is shown on each axis.

We also tested to what extent the measured soil properties could help explain the observed segregation in the microbial communities for the two plant species. Our results showed that the variation was related to several physicochemical parameters, including NH4+ content, pH and EC. Higher levels of Br, pH and EC were associated with microbial communities from D. spicata rhizosphere according to the RDA1 axis, and NH4+ content were clearly associated with microbial communities from S. foliosa (Figure 6; Supplementary Table 2).

Figure 6
Four bar graphs labeled A to D compare different treatments on soil characteristics. Graph A shows Electrical Conductivity (EC) with higher measures in D. spicata. Graph B displays ammonium content with higher values for S. foliosa. Graph C shows Dissolved Organic Carbon (DOC) dominated by S. foliosa. Graph D indicates Dissolved Organic Nitrogen (DON) with a peak in D. spicata SST2. Control, D. spicata, and S. foliosa are represented by different colors. Error bars indicate variability.

Figure 6. Redundancy analysis (RDA) of rhizosphere bacterial community composition collected from under two native plant species (D. spicata, S. foliosa) alongside control unvegetated soil in the hyper-arid core of the Atacama Desert. The RDA was based on Hellinger-transformed amplicon sequence variant (ASV) relative abundances. Points represent individual soil samples (colored by plant species for rhizosphere soil, gray for the control soil) with polygons grouping samples by plant species. Arrows indicate the explanatory power of physicochemical parameters in describing the observed variation in bacterial community composition. Percentages on axes show explained variance in the unconstrained and constrained analyses.

Discussion

Environmental conditions in the hyper-arid of the Atacama Desert

In this study, we characterized the rhizosphere microbial community associated with D. spicata and S. foliosa which are located in the hyper-arid core of the Atacama Desert, one of the most extreme environments on Earth. The climatic factors, mainly low precipitation, high day-night temperature regimes and UV irradiation, limit biological productivity and activity, chemical reactions and weathering process (Plaza et al., 2018). Further, these hyper-arid soils possess a low soil organic C content, low microbial activity, and low rates of mineral weathering (Ewing et al., 2006; Knief et al., 2020; Mörchen et al., 2019). However, in this study, we show that it is possible to find an abundance of microorganisms even in hyper-arid soils, specifically in the rhizosphere from native plant species such as D. spicata and S. foliosa. The differences in soil physicochemical parameters in the rhizosphere soil in comparison to the unvegetated control site, alongside shifts in microbial community composition, strongly suggest that the plants modulate their environment, as described in other sites of the Atacama Desert associated to flowering desert phenomena in the north of Chile (Araya et al., 2020). Our three replicate sample stations (SST1, SST2 and SST3) showed significant variations in NH4+, NO3, Olsen-P, pH, EC, and soluble C (Supplementary Table 1), while CaCO3 and moisture content showed no difference (Supplementary Table 1). According to the study area, the wet layers below the surface of the hyper arid core of the Atacama have been reported (Azua-Bustos et al., 2020). It should be noted that we sampled the rhizospheric soil (20–25 cm depth) and were not able to sample the deep soil layers close to the water table. Although the origin of water is outside the scope of this study, the deep-water sources could be related to the unusual rain events as have been indicated by Azua-Bustos et al. (2020). Precipitation events could lead to a reload of a deep groundwater system present in the Aguas Blancas basin (Dirección General de Aguas, 2017). It has a low recharge rate, limited to very sporadic and rare rainfall events (Herrera and Custodio, 2014). Previous reports have shown that the hyper-arid core of the Atacama Desert experiences large day-night variations in air temperature and relative humidity (Azua-Bustos et al., 2015). When combined with low amounts of rainfall, soil moisture, soil surface compaction, high salinity and nutrient imbalance it makes it a highly challenging environment for plants to survive (Marquet et al., 1998). Despite these environmental conditions, we can find dominant plants in the Yungay oasis such as the salty grass D. spicata, which has a water-efficient C4 metabolism, is hyper-salt tolerant can survive at >0.5 M NaCl and, can withstand temperatures of up to 57 °C (Golden et al., 1995; Lazarus et al., 2011). The drought-tolerant shrub, S. foliosa, with C3 metabolism represents the other dominant higher plant species in the Yungay Oasis (Akhani et al., 1997; Nakano and Asada, 1981) and, has previously been used for the phytoremediation of salt-affected soils (e.g., high Valley of Cochabamba in Bolivia; de Froidmont, 2018).

Microbial communities associated with native plants in the hyper-arid core of the Atacama Desert

The microbial community in the rhizosphere soil was dominated by Firmicutes (Bacillota), Proteobacteria (Pseudomonadota), Halobacteria and Actinobacteria (Figure 3). However, the microbial community composition may vary in different situations, for example, if we compared the soil surface with subsoil layers, or the interaction of other soil fractions with the rhizosphere is analyzed (Castro-Severyn et al., 2024; Crits-Christoph et al., 2016; Fuentes et al., 2022a; Warren-Rhodes et al., 2019). Our results indicate that the taxonomic composition of D. spicata rhizospheric soil exhibited more variability compared to S. foliosa, which accounts for a more stable or selected community. In this context, some taxa were found to be very abundant at some sampling sites while being completely absent or having very low representation at others (Figure 4). Similar patterns have been reported in arid environments such as Atacama and North China Desert, which showed high variability in their soil microbial communities (Vásquez-Dean et al., 2020). Additionally, studies in the Namib Desert have reported differences in the microbial communities associated with the roots of Tetraena simplex, Tetraena stapffi and Stipagrostis sp. (Lopes et al., 2024). This may be related to the plant’s need to generate interaction with more specialized organisms capable of tolerating the high salinity or pH changes induced in the soil, which are caused by the plant’s physiologic processes to eliminate salts and thus thrive under desert conditions (Pelliza et al., 2005; Conticello et al., 2002). Under salt stress, certain adapted crops show enhanced root performance (Vaseva et al., 2021). This adaptation is supported by alterations in plant processes (e.g., Na+/K+ transporter selectivity) as well as microbial mechanisms such as the production of compounds such as indole acetic acid or the enhancement of antioxidant enzyme activity, which can help the plant tolerate stress conditions (Fortt et al., 2022; Khalid and Aftab, 2020; Zhao et al., 2021). Our analysis revealed dominance by the Proteobacteria (Pseudomonadota), Firmicutes (Bacillota), Halobacteria and Actinobacteria phyla. At the genus level, significant differences between the D. spicata and S. foliosa rhizosphere communities were observed for Natrinema, Buchnera, Enterococcus and Salibacterium (Figure 4). Similar patterns of microbial distribution have been observed in other extreme environments, where genera like Natrinema and Haloterrigena are known to dominate hypersaline soils due to their ability to thrive in high-salinity environments, like saline lakes or salt pans (Crognale et al., 2013; Jiang et al., 2012; Mirete et al., 2015). These genera have been consistently associated with saline environments, suggesting that their presence in the D. spicata rhizosphere could be directly related to the high salt concentrations observed in these soils (Jones et al., 2023; Remonsellez et al., 2018). Furthermore, Buchnera and Enterococcus, typically associated with gut symbiosis, have also been found in plant environments, where their presence may be explained by their role in N cycling or plant growth promotion in stressful conditions (Bartelme et al., 2018; Zhang et al., 2020). Moreover, the dominance of Firmicutes (Bacillota) in S. foliosa rhizosphere, particularly Enterococcus and Buchnera, could be an adaptive response to the ammonium-rich environment, as seen in other studies where these genera have thrived in nutrient-stressed ecosystems related to phosphate limitation (Knief et al., 2020; Neilson et al., 2017). We ascribe the intrinsic variability observed in microbial composition between sampling sites to microenvironmental differences, as extreme conditions in the Atacama Desert are typically highly heterogeneous even within small areas (Crits-Christoph et al., 2016). Also, the assembly mechanisms and ecosystem functions of hypolithic communities from the Qaidam Desert in China, which occupy highly specialized microenvironments beneath translucent rocks in arid environments, are not yet fully understood (Lai et al., 2024). As well as a study from Chihuahuan Desert in Mexico related to plant-bacteria interaction, specifically the promote facilitation from nurse plant relationship have been reported that facilitation is enhanced by plant species that provide a more homogeneous microenvironment (Sánchez-Martín et al., 2024).

Other studies demonstrated that Actinobacteria and Chloroflexi dominate soil microbial communities in the hyper-arid margin of the Atacama Desert with vegetable cover (Neilson et al., 2012, 2017). Furthermore, most recent studies showed the dominance of different groups of Actinobacteria, Proteobacteria and Chloroflexi in Atacama Desert soils (Knief et al., 2020), while Connon et al. (2007) detected Actinobacteria, Proteobacteria (Pseudomonadota) and Firmicutes (Bacillota) in Yungay soil (0–20 cm depth), similar to our findings. Our analysis revealed a strong correlation between EC and the abundance of Halobacteria in the D. spicata rhizosphere, reflecting the relationship between soil salinity and halophilic bacterial communities. Furthermore, the higher ammonium content in the S. foliosa rhizosphere could favor Firmicutes phyla, although they are not commonly classed as ammoniaphilic and, the high Br content may be related to salt excretion by D. spicata and the geochemical composition of the surrounding soil (Morris et al., 2019; Figure 6). Soil moisture is essential for bacterial diversity in desert soils (Bottos et al., 2020), but in this work the moisture does not present statistical differences between the rhizosphere of two plant species, perhaps because in 2022 the Yungay area experienced a slight episode of precipitation and fog (Weather and Climate, 2022). Castro-Severyn et al. (2024) demonstrated that rhizosphere soil enrichment cultures from S. foliosa and D. spicata were dominated by Klebsiella, followed by Brevibacillus, reflecting the selection patterns at the phylum level Proteobacteria (Pseudomonadota) and Firmicutes (Bacillota), respectively for both plant species. In this context, different species of the Klebsiella genus have been repeatedly detected in many soil types and environmental samples, including soils from the Yungay area (Ekwanzala et al., 2019; Thomas et al., 2006). The soil moisture could be ancillary to other factors in shaping bacterial diversity and, the clay-rich soils can act as a possible “water reservoir” and help shape microbial life conditions (Fuentes et al., 2022a). Moreover, the salinity of the soil, reflected by the EC, was high in soil from the D. spicata rhizosphere (Figure 2A), and it seems to shape bacterial communities due to intense selective pressure, as few bacteria are capable of growing over large gradients of salt concentrations (Figures 3, 5). We measured high EC values (over 64 to 92 dS/m) in D. spicata rhizosphere, while in S. foliosa rhizosphere, the EC values fluctuated between 18 and 24 dS/m (Supplementary Table 1). In this regard, we observed a particular accumulation of halophilic taxa in the D. spicata rhizosphere; specifically, in samples from Sampling Station 3. This halophytic trait is reflected in the high contents of NaCl in the leaves, the salt crystals presence on the leaf surface in D. spicata and the presence of high amounts of salt in the phyllosphere soil (Jones et al., 2023). In the same work, Jones et al. (2023) presented a comparison between vegetated and unvegetated soils, suggesting that D. spicata is effective at removing salt from the soil and translocating it to above-ground component (e.g., much lower levels of Na were seen in the roots relative to the shoots and in the subsoil relative to the phyllosphere soil). These observations are consistent with the foliar Na excretion through salt glands in this plant, as reported in other studies (Hasanuzzaman et al., 2014; Semenova et al., 2010). Furthermore, it is interesting to note that some of the bacterial families present in the soil profile are known to include halophilic taxa or have been shown to be present in other saline environments, such as Proteobacteria (Pseudomonadota), Actinobacteria, and Firmicutes (Bacillota) which have been detected in different saline environments (Bartelme et al., 2018; Crognale et al., 2013; Jiang et al., 2012; Mirete et al., 2015; Remonsellez et al., 2018; Zhang et al., 2020).

Ecology importance of S. foliosa and D. spicata microbial communities

To date, only one study has reported NH3 emissions occurring during the microbial processing of organic-N in the Atacama Desert (Jones et al., 2023), this information can explain the high levels of nitrogen compounds (TN and NO3) in our samples from D. spicata. Furthermore, different studies have reported an association between microbial richness and diversity in the Atacama Desert hyper-arid soils with water availability, relative humidity and TC content (Crits-Christoph et al., 2016; Fuentes et al., 2022a; Knief et al., 2020; Neilson et al., 2017; Schulze-Makuch et al., 2018). Furthermore, it would be important to highlight the functional potential of the microbial community associated with D. spicata and S. foliosa. In this context, there could be related physiological functions such as salt excretion in D. spicata plants (Yuan et al., 2016; Morris et al., 2019) with microbial functions such as osmolytes biosynthesis. Moreover, the high ammonium content in the S. foliosa rhizosphere could be an ammonium accumulation under stress conditions (Rütting et al., 2011). In this respect, other studies have reported microbial functions related to nitrate and nitrite ammonification in plants from arid environments (Kuypers et al., 2018). Despite this knowledge, to date there has been no metagenomic data relating to microbial community functions and plant-bacteria interactions from soils of the hyper-arid core of the Atacama Desert. Therefore, the search for life in extreme environments should consider bacterial, archaea and fungi communities functioning to better understand how native plants thrive under the combined abiotic stresses operating in the Atacama Desert.

Conclusion

The hyper-arid core of the Atacama Desert represents one of the driest ecosystems on Earth, while the prevailing environmental conditions make it one of the most extreme places for life to establish. In this study, we characterized the rhizosphere bacterial communities and soil properties associated with two native plant species (D. spicata and S. foliosa) in the Yungay Oasis. In this context, it is important to remark that to date no studies related to plant-bacteria interaction in this area from the hyper arid core of the Atacama Desert. Our analyses revealed distinct microbial community structures correlating with specific soil parameters, particularly EC and ammonium content. This study advances our fundamental understanding of plant-microbe interactions in extreme environments by identifying bacterial communities capable of thriving in the hyper-arid core of the Atacama Desert and their community structure depends on the soil parameters and plant interaction in the rhizospheric soil environment. Future research should focus on functional studies by metagenomic approach to elucidate the specific roles of these microbial communities in plant-bacteria interactions within these extreme soils.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found at: https://www.ncbi.nlm.nih.gov/bioproject/, accession number PRJNA971922.

Author contributions

JF: Conceptualization, Writing – original draft, Writing – review & editing, Data curation, Formal analysis, Investigation, Methodology, Software. JC-S: Conceptualization, Formal analysis, Investigation, Supervision, Validation, Writing – review & editing. AC: Formal analysis, Investigation, Methodology, Writing – review & editing. GD: Investigation, Methodology, Writing – review & editing. AS: Funding acquisition, Supervision, Writing – review & editing. DJ: Formal analysis, Investigation, Methodology, Supervision, Writing – review & editing. CS: Formal analysis, Funding acquisition, Writing – review & editing. BF: Writing – review & editing, Formal analysis, Funding acquisition. FR: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was sponsored by ANID (Agencia Nacional de Investigación y Desarrollo de Chile) grants; Regular FONDECYT 1220902; 2022 Estrategia en Sequía FSEQ210029 to Castro-Severyn JC-S and FR; 2023 ANID National Doctoral Fellowship N° 21230757 to JF, Atracción e capital humano avanzado del extranjero: Modalidad Estadías Cortas MEC80190012 to JF and BF Regular FONDECYT 1210633 and 2023 Anillo ATE220007 to CPS.

Acknowledgments

We thank Emily Cooledge for help with the laboratory work. The funders had no role in study design, data collection, or analysis.

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.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1587491/full#supplementary-material

References

Akhani, H., Trimborn, P., and Ziegler, H. (1997). Photosynthetic pathways in Chenopodiaceae from Africa, Asia and Europe with their ecological, phytogeographical and taxonomical importance. Plant Syst. Evol. 206, 187–206.

Google Scholar

Andersen, K. S., Kirkegaard, R. H., Karst, S. M., and Albertsen, M. (2018). Ampvis2: an R package to analyse and visualise 16S rRNA amplicon data. BioRxiv, 299537. doi: 10.1101/299537

Crossref Full Text | Google Scholar

Araya, J. P., González, M., Cardinale, M., Schnell, S., and Stoll, A. (2020). Microbiome dynamics associated with the Atacama flowering desert. Front. Microbiol. 10:3160. doi: 10.3389/fmicb.2019.03160

PubMed Abstract | Crossref Full Text | Google Scholar

Azua-Bustos, A., Caro-Lara, L., and Vicuña, R. (2015). Discovery and microbial content of the driest site of the hyperarid Atacama Desert, Chile. Environ. Microbiol. Rep. 7, 388–394. doi: 10.1111/1758-2229.12261

PubMed Abstract | Crossref Full Text | Google Scholar

Azua-Bustos, A., González-Silva, C., and Corsini, G. (2017). The hyperarid core of the Atacama Desert, an extremely dry and carbon deprived habitat of potential interest for the field of carbon science. Front. Microbiol. 8:933. doi: 10.3389/fmicb.2017.00993

PubMed Abstract | Crossref Full Text | Google Scholar

Azua-Bustos, A., Fairén, A. G., Silva, C. G., Carrizo, D., Fernández-Martínez, M. Á., Arenas-Fajardo, C., et al. (2020). Inhabited subsurface wet smectites in the hyperarid core of the Atacama Desert as an analog for the search for life on Mars. Scientific Reports. 10:19183. doi: 10.1038/s41598-020-76302-z

Crossref Full Text | Google Scholar

Bartelme, R. P., Oyserman, B. O., Blom, J. E., Sepulveda-Villet, O. J., and Newton, R. J. (2018). Stripping away the soil: plant growth promoting microbiology opportunities in aquaponics. Front. Microbiol. 9:8. doi: 10.3389/fmicb.2018.00008

PubMed Abstract | Crossref Full Text | Google Scholar

Bashan, Y., de- Bashan, L. E., Prabhu, S. R., and Hernandez, J. P. (2014). Advances in plant growth-promoting bacterial inoculant technology: formulations and practical perspectives (1998-2013). Plant Soil 378, 1–33. doi: 10.1007/s11104-013-1956-x

Crossref Full Text | Google Scholar

Bonfante, P., and Anca, I. A. (2009). Plants, mycorrhizal fungi, and bacteria: a network of interactions. Ann. Rev. Microbiol. 63, 363–383. doi: 10.1146/annurev.micro.091208.073504

Crossref Full Text | Google Scholar

Bottos, E. M., Laughlin, D. C., Herbold, C. W., Lee, C. K., McDonald, I. R., and Cary, S. C. (2020). Abiotic factors influence patterns of bacterial diversity and community composition in the dry valleys of Antarctica. FEMS Microbiol. Ecol. 96:fiaa042. doi: 10.1093/FEMSEC/FIAA042

PubMed Abstract | Crossref Full Text | Google Scholar

Bull, A. T., and Asenjo, J. A. (2013). Microbiology of hyper-arid environments: recent insights from the Atacama Desert, Chile. Anton. Leeuw. Int. J. Gen. Mol. Microbiol. 103, 1173–1179. doi: 10.1007/s10482-013-9911-7

PubMed Abstract | Crossref Full Text | Google Scholar

Bull, A. T., Asenjo, J. A., Goodfellow, M., and Gómez-Silva, B. (2016). The Atacama Desert: technical resources and the growing importance of novel microbial diversity. Ann. Rev. Microbiol. 70, 215–234. doi: 10.1146/annurev-micro-102215-095236

PubMed Abstract | Crossref Full Text | Google Scholar

Calderon, G., Garrido, M., and Acevedo, E. (2015). Prosopis tamarugo Phil.: a native tree from the Atacama Desert groundwater table depth thresholds for conservation. Rev. Chil. Hist. Nat. 88, 1–8. doi: 10.1186/s40693-015-0048-0

Crossref Full Text | Google Scholar

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | Crossref Full Text | Google Scholar

Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 108, 4516–4522. doi: 10.1073/pnas.1000080107

PubMed Abstract | Crossref Full Text | Google Scholar

Castro-Severyn, J., Fortt, J., Sierralta, M., Alegria, P., Donoso, G., Choque, A., et al. (2024). Rhizospheric bacteria from the Atacama Desert hyper-arid core: cultured community dynamics and plant growth promotion. Microbiol. Spectrum 12:e0005624. doi: 10.1128/spectrum.00056-24

PubMed Abstract | Crossref Full Text | Google Scholar

Chávez, R. O., Clevers, J. G. P. W., Decuyper, M., de Bruin, S., and Herold, M. (2016). 50 years of water extraction in the Pampa del Tamarugal basin: can Prosopis tamarugo trees survive in the hyper-arid Atacama Desert (northern Chile)? J. Arid Environ. 124, 292–303. doi: 10.1016/j.jaridenv.2015.09.007

Crossref Full Text | Google Scholar

Clarke, J. D. A. (2006). Antiquity of aridity in the Chilean Atacama Desert. Geomorphology 73, 101–114. doi: 10.1016/j.geomorph.2005.06.008

PubMed Abstract | Crossref Full Text | Google Scholar

Conley, C. A., Ishkhanova, G., McKay, C. P., and Cullings, K. (2006). A preliminary survey of non-lichenized fungi cultured from the hyperarid Atacama Desert of Chile. Astrobiology 6, 521–526. doi: 10.1089/ast.2006.6.521

PubMed Abstract | Crossref Full Text | Google Scholar

Connon, S. A., Lester, E. D., Shafaat, H. S., Obenhuber, D. C., and Ponce, A. (2007). Bacterial diversity in hyperarid Atacama Desert soils. J. Geophys. Res. Biogeosci. 112, 1–9. doi: 10.1029/2006JG000311

Crossref Full Text | Google Scholar

Conticello, L., Cerazo, B., and Bustamante, A. (2002). Dinámica de comunidades hidrófilas asociadas a canales de riego en el Alto Valle del Río Negro (Argentina). Gayana. Bot. 59, 13–20. doi: 10.4067/S0717-66432002000100003

Crossref Full Text | Google Scholar

Crits-Christoph, A., Gelsinger, D. R., Ma, B., Wierzchos, J., Ravel, J., Davila, A., et al. (2016). Functional interactions of archaea, bacteria and viruses in a hypersaline endolithic community. Environmental microbiology. 18, 2064–2077. doi: 10.1111/1462-2920.13259

Crossref Full Text | Google Scholar

Crognale, S., Máthé, I., Cardone, V., Stazi, S. R., and Ráduly, B. (2013). Halobacterial community analysis of Mierlei saline lake in Transylvania (Romania). Geomicrobiol J. 30, 801–812. doi: 10.1080/01490451.2013.774073

Crossref Full Text | Google Scholar

de Froidmont, C. (2018). Contribution to the identification of solution for the remediation of salt - affected soils of the high valley of Cochabamba in Bolivia: use of gypsum, Sulphur and Suaeda Foliosa (master thesis, Gembloux agro-bio tech, University of Liege.

Google Scholar

Dhanker, R., Chaudhary, S., Kumari, A., Kumar, R., and Goyal, S. (2020). Symbiotic signaling: insights from arbuscular mycorrhizal symbiosis. In: Varma, A., Tripathi, S., Prasad, R. (eds.). Plant Microbe Symbiosis. Springer, Cham. 75–103. doi: 10.1007/978-3-030-36248-5_5

Crossref Full Text | Google Scholar

Dirección General de Aguas (2017). Balance de Gestión Integral. Santiago, Chile, Ministerio de Obras Públicas. Avaliable online at: https://dga.mop.gob.cl/uploads/sites/13/2024/08/BGI_DGA_2017.pdf

Google Scholar

Egamberdieva, D., Wirth, S. J., Shurigin, V. V., Hashem, A., and Abd Allah, E. F. (2017). Endophytic bacteria improve plant growth, symbiotic performance of chickpea (Cicer arietinum L.) and induce suppression of root rot caused by Fusarium solani under salt stress. Front. Microbiol. 8:1887. doi: 10.3389/fmicb.2017.01887

PubMed Abstract | Crossref Full Text | Google Scholar

Ekwanzala, M. D., Dewar, J. B., Kamika, I., and Momba, M. N. B. (2019). Tracking the environmental dissemination of carbapenem-resistant Klebsiella pneumoniae using whole genome sequencing. Sci. Total Environ. 691, 80–92. doi: 10.1016/j.scitotenv.2019.06.533

PubMed Abstract | Crossref Full Text | Google Scholar

Ewing, S. A., Sutter, B., Owen, J., Nishiizumi, K., Sharp, W., Cliff, S. S., et al. (2006). A threshold in soil formation at earth’s arid-hyperarid transition. Geochim. Cosmochim. Acta 70, 5293–5322. doi: 10.1016/j.gca.2006.08.020

Crossref Full Text | Google Scholar

Ewing, S. A., Yang, W., DePaolo, D. J., Michalski, G., Kendall, C., Stewart, B. W., et al. (2008). Non-biological fractionation of stable ca isotopes in soils of the Atacama Desert, Chile. Geochim. Cosmochim. Acta 72, 1096–1110. doi: 10.1016/j.gca.2007.10.029

Crossref Full Text | Google Scholar

FAO (2021). El estado mundial de la agricultura y la alimentación. Food & Agriculture Organization of the United Nations, Rome.

Google Scholar

Fletcher, L. E., Valdivia-Silva, J. E., Perez-Montaño, S., Condori-Apaza, R. M., Conley, C. A., and McKay, C. P. (2012). Variability of organic material in surface horizons of the hyper-arid Mars-like soils of the Atacama Desert. Adv. Space Res. 49, 271–279. doi: 10.1016/j.asr.2011.10.001

Crossref Full Text | Google Scholar

Fortt, J., González, M., Morales, P., Araya, N., Remonsellez, F., Coba De La Peña, T., et al. (2022). Bacterial modulation of the plant ethylene signaling pathway improves tolerance to salt stress in lettuce (Lactuca sativa L.). Front. Sustain. Food Syst. 6:768250. doi: 10.3389/fsufs.2022.768250

Crossref Full Text | Google Scholar

Fuentes, B., Choque, A., Gómez, F., Alarcón, J., Castro-Nallar, E., Arenas, F., et al. (2022a). Influence of physical-chemical soil parameters on microbiota composition and diversity in a deep hyperarid core of the Atacama Desert. Front. Microbiol. 12:794743. doi: 10.3389/fmicb.2021.794743

PubMed Abstract | Crossref Full Text | Google Scholar

Fuentes, B., Gómez, F., Valdez, C., Videla, A., Castro-Severyn, J., Barahona, S., et al. (2022b). Effects of altitude on soil properties in coastal fog ecosystems in Morro Moreno National Park, Antofagasta, Chile. Eur. J. Soil Sci. 73:e13217. doi: 10.1111/ejss.13217

Crossref Full Text | Google Scholar

Fuentes, A., Herrera, H., Charles, T. C., and Arriagada, C. (2020). Fungal and bacterial microbiome associated with the rhizosphere of native plants from the Atacama Desert. Microorganisms 8:209. doi: 10.3390/microorganisms8020209

PubMed Abstract | Crossref Full Text | Google Scholar

Garrido, M., Silva, P., and Acevedo, E. (2016). Water relations and foliar isotopic composition of Prosopis tamarugo phil., an endemic tree of the Atacama Desert growing at three levels of water table depth. Front. Plant Sci. 7:375. doi: 10.3389/fpls.2016.00375

PubMed Abstract | Crossref Full Text | Google Scholar

Golden, A. M., Baldwin, J. G., and Mundo-Ocampo, M. (1995). Description of Tylenchorhynchus thermophilus n. sp. (nematoda: tylenchina) from saltgrass in Death Valley, California. J. Nematol. 27:312.

Google Scholar

Hartley, A. J., Chong, G., Houston, J., and Mather, A. E. (2005). 150 million years of climatic stability: evidence from the Atacama Desert, northern Chile. J. Geol. Soc. 162, 421–424. doi: 10.1144/0016-764904-071

Crossref Full Text | Google Scholar

Hasanuzzaman, M., Nahar, K., Alam, M. M., Bhowmik, P. C., Hossain, M. A., Rahman, M. M., et al. (2014). Potential use of halophytes to remediate saline soils. Biomed. Res. Int. 2014:589341. doi: 10.1155/2014/589341

Crossref Full Text | Google Scholar

Herrera, C., and Custodio, E. (2014). Origen de las aguas de pequeños manantiales de la Costa del norte de Chile, en las cercanías de Antofagasta. Andean Geol. 41, 314–341. doi: 10.5027/andgeoV41n2-a03

Crossref Full Text | Google Scholar

Jiang, F., Cao, S. J., Li, Z. H., Fan, H., Li, H. F., Liu, W. J., et al. (2012). Salisediminibacterium halotolerans gen. Nov., sp. nov., a halophilic bacterium from soda lake sediment. Int. J. Syst. Evol. Microbiol. 62, 2127–2132. doi: 10.1099/ijs.0.034488-0

PubMed Abstract | Crossref Full Text | Google Scholar

Jones, D. L., Fuentes, B., Arenas-Díaz, F., Remonsellez, F., van Hall, R., Atkinson, B. S., et al. (2023). Life at the extreme: plant-driven hotspots of soil nutrient cycling in the hyper-arid core of the Atacama Desert. Soil Biol. Biochem. 184:109128. doi: 10.1016/j.soilbio.2023.109128

Crossref Full Text | Google Scholar

Jones, D. L., Magthab, E. A., Gleeson, D. B., Hill, P. W., Sánchez-Rodríguez, A. R., Roberts, P., et al. (2018a). Microbial competition for nitrogen and carbon is as intense in the subsoil as in the topsoil. Soil Biol. Biochem. 117, 72–82. doi: 10.1016/j.soilbio.2017.10.024

Crossref Full Text | Google Scholar

Jones, D. L., Olivera-Ardid, S., Klumpp, E., Knief, C., Hill, P. W., Lehndorff, E., et al. (2018b). Moisture activation and carbon use efficiency of soil microbial communities along an aridity gradient in the Atacama Desert. Soil Biol. Biochem. 117, 68–71. doi: 10.1016/j.soilbio.2017.10.026

Crossref Full Text | Google Scholar

Khalid, A., and Aftab, F. (2020). Effect of exogenous application of IAA and GA3 on growth, protein content, and antioxidant enzymes of Solanum tuberosum L. grown in vitro under salt stress. In Vitro Cellular Develop. Biol. Plant 56, 377–389. doi: 10.1007/s11627-019-10047-x

Crossref Full Text | Google Scholar

Kim, Y., Glick, B., Bashan, Y., and Ryu, C. (2013). “Enhancement of plant drought tolerance by microbes” in Plant responses to drought stress. ed. R. Aroca (Berlin, Heidelberg: Springer), 383–413.

Google Scholar

Knief, C., Bol, R., Amelung, W., Kusch, S., Frindte, K., Eckmeier, E., et al. (2020). Tracing elevational changes in microbial life and organic carbon sources in soils of the Atacama Desert. Glob. Planet. Chang. 184:103078. doi: 10.1016/j.gloplacha.2019.103078

Crossref Full Text | Google Scholar

Kuypers, M. M., Marchant, H. K., and Kartal, B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol, 16, 263–276. doi: 10.1038/nrmicro.2018.9

Crossref Full Text | Google Scholar

Lacap, D. C., Warren-Rhodes, K. A., McKay, C. P., and Pointing, S. B. (2011). Cyanobacteria and chloroflexi-dominated hypolithic colonization of quartz at the hyper-arid core of the Atacama Desert, Chile. Extremophiles 15, 31–38. doi: 10.1007/s00792-010-0334-3

PubMed Abstract | Crossref Full Text | Google Scholar

Lai, Z., Liu, Z., Zhao, Y., Qin, S., Zhang, W., Lang, T., et al. (2024). Distinct microbial communities under different rock-associated microhabitats in the Qaidam Desert. Environ. Res. 250:118462. doi: 10.1016/j.envres.2024.118462

PubMed Abstract | Crossref Full Text | Google Scholar

Lazarus, B. E., Richards, J. H., Gordon, P. E., Oki, L. R., and Barnes, C. S. (2011). Plasticity tradeoffs in salt tolerance mechanisms among desert Distichlis spicata genotypes. Funct. Plant Biol. 38, 187–198. doi: 10.1071/FP10192

PubMed Abstract | Crossref Full Text | Google Scholar

Lin, Y. C., Chang, L. W., Yang, K. C., Wang, H. H., and Sun, I. F. (2011). Point patterns of tree distribution determined by habitat heterogeneity and dispersal limitation. Oecologia 165, 175–184. doi: 10.1007/s00442-010-1718-x

PubMed Abstract | Crossref Full Text | Google Scholar

Lopes, T., Santos, J., Matos, D., Sá, C., Pina, D., Pinto, R., et al. (2024). Soil bacteria from the Namib Desert: insights into plant growth promotion and osmotolerance in a hyper-arid environment. Land 13:1678. doi: 10.3390/land13101678

Crossref Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

Makhalanyane, T. P., Valverde, A., Gunnigle, E., Frossard, A., Ramond, J. B., and Cowan, D. A. (2015). Microbial ecology of hot desert edaphic systems. FEMS Microbiol. Rev. 39, 203–221. doi: 10.1093/femsre/fuu011

PubMed Abstract | Crossref Full Text | Google Scholar

Marquet, P. A., Bozinovic, F., Bradshaw, G., Cornelius, C., Gonzalez, H., Gutierrez, J., et al. (1998). Los ecosistemas del desierto de Atacama y área andina adyacente en el norte de Chile. Rev. Chil. Hist. Nat. 71, 593–617.

Google Scholar

McCray, J. M., Wright, A. L., Luo, Y., and Ji, S. (2012). Soil phosphorus forms related to extractable phosphorus in the Everglades agricultural area. Soil Sci. 177, 31–38. doi: 10.1097/SS.0b013e31823782da

Crossref Full Text | Google Scholar

Mckay, C., Friedmann, E., Gómez-Silva, B., Cáceres-Villanueva, L., and Andersen, D. (2003). Temperature and moisture conditions for life in the extreme arid region of the Atacama Desert: four years of observations including el Niño of 1997-1998. Astrobiology 3, 393–406. doi: 10.1089/15311070376901646

Crossref Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217

PubMed Abstract | Crossref Full Text | Google Scholar

Miranda, K. M., Espey, M. G., and Wink, D. A. (2001). A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide 5, 62–71. doi: 10.1006/niox.2000.0319

PubMed Abstract | Crossref Full Text | Google Scholar

Mirete, S., Mora-Ruiz, M. R., Lamprecht-Grandío, M., de Figueras, C. G., Rosselló-Móra, R., and González-Pastor, J. E. (2015). Salt resistance genes revealed by functional metagenomics from brines and moderate-salinity rhizosphere within a hypersaline environment. Front. Microbiol. 6:1121. doi: 10.3389/fmicb.2015.01121

PubMed Abstract | Crossref Full Text | Google Scholar

Mörchen, R., Lehndorff, E., Diaz, F. A., Moradi, G., Bol, R., Fuentes, B., et al. (2019). Carbon accrual in the Atacama Desert. Glob. Planet. Chang. 181:102993. doi: 10.1016/j.gloplacha.2019.102993

Crossref Full Text | Google Scholar

Morris, L., Yun, K., Rutter, A., and Zeeb, B. A. (2019). Characterization of excreted salt from the Recretohalophytes Distichlis spicata and Spartina pectinata. J. Environ. Qual. 48, 1775–1780. doi: 10.2134/jeq2019.03.0102

Crossref Full Text | Google Scholar

Mulvaney, R. L. (1996). “Nitrogen inorganic forms” in Methods of soil analysis: Part 3 chemical methods. eds. D. L. Sparks, A. L. Page, P. A. Helmke, R. H. Loeppert, and P. N. Soltanpour (Madison, WI: American Society of Agronomy and Soil Science Society of America), 1123–1184.

Google Scholar

Nakano, Y., and Asada, K. (1981). Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts. Plant Cell Physiol. 22, 867–880. doi: 10.1093/oxfordjournals.pcp.a076232

Crossref Full Text | Google Scholar

Navarro-González, R., Rainey, F., Molina, P., Bagaley, D., Hollen, B., de la Rosa, J., et al. (2003). Mars-like soils in the Atacama Desert, Chile, and the dry limit of microbial life. Science 302, 1018–1021. doi: 10.1126/science.1089143

Crossref Full Text | Google Scholar

Neilson, J. W., Califf, K., Cardona, C., Copeland, A., van Treuren, W., Josephson, K. L., et al. (2017). Significant impacts of increasing aridity on the arid soil microbiome. MSystems 2, 10–1128. doi: 10.1128/mSystems.00195-16

PubMed Abstract | Crossref Full Text | Google Scholar

Neilson, J. W., Quade, J., Ortiz, M., Nelson, W. M., Legatzki, A., Tian, F., et al. (2012). Life at the hyperarid margin: novel bacterial diversity in arid soils of the Atacama Desert, Chile. Extremophiles 16, 553–566. doi: 10.1007/s00792-012-0454-z

PubMed Abstract | Crossref Full Text | Google Scholar

Niu, X., Song, L., Xiao, Y., and Ge, W. (2018). Drought-tolerant plant growth-promoting rhizobacteria associated with foxtail millet in a semi-arid and their potential in alleviating drought stress. Front. Microbiol. 8:2580. doi: 10.3389/fmicb.2017.02580

Crossref Full Text | Google Scholar

Pelliza, A., Borrelli, L., and Bonvissuto, G. (2005). El Pasto Salado (Distichlis Spp.) En La Patagonia: Una Forrajera adaptada a la aridez y a la salinidad. Rev. Cient. Agropecu. 9, 119–131.

Google Scholar

Plaza, C., Zaccone, C., Sawicka, K., Méndez, A. M., Tarquis, A., Gascó, G., et al. (2018). Soil resources and element stocks in drylands to face global issues. Sci. Rep. 8:13788. doi: 10.1038/s41598-018-32229-0

PubMed Abstract | Crossref Full Text | Google Scholar

Pointing, S. B., and Belnap, J. (2012). Microbial colonization and controls in dryland systems. Nat. Rev. Microbiol. 10, 551–562. doi: 10.1038/nrmicro2831

PubMed Abstract | Crossref Full Text | Google Scholar

Poulos, H. M., and Camp, A. E. (2010). Topographic influences on vegetation mosaics and tree diversity in the Chihuahuan desert borderlands. Ecology 91, 1140–1151. doi: 10.1890/08-1808.1

PubMed Abstract | Crossref Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). Fasttree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | Crossref Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team, R. (2020). R: A language and environment for statistical computing [Manual].

Google Scholar

Remonsellez, F., Castro-Severyn, J., Pardo-Esté, C., Aguilar, P., Fortt, J., Salinas, C., et al. (2018). Characterization and salt response in recurrent halotolerant Exiguobacterium sp. SH31 isolated from sediments of Salar de Huasco, chilean altiplano. Front. Microbiol. 9:228. doi: 10.3389/fmicb.2018.02228

PubMed Abstract | Crossref Full Text | Google Scholar

Richardson, A. E., Barea, J. M., McNeill, A. M., and Prigent-Combaret, C. (2009). Acquisition of phosphorus and nitrogen in the rhizosphere and plant growth promotion by microorganisms. Plant Soil 321, 305–339. doi: 10.1007/s11104-009-9895-2

Crossref Full Text | Google Scholar

Ritter, B., Binnie, S. A., Stuart, F. M., Wennrich, V., and Dunai, T. J. (2018). Evidence for multiple Plio-Pleistocene lake episodes in the hyperarid Atacama Desert. Quat. Geochronol. 44, 1–12. doi: 10.1016/j.quageo.2017.11.002

Crossref Full Text | Google Scholar

Ruiz-Sánchez, M., Aroca, R., Muñoz, Y., Polón, R., and Ruiz-Lozano, J. M. (2010). The arbuscular mycorrhizal symbiosis enhances the photosynthetic efficiency and the antioxidative response of rice plants subjected to drought stress. J. Plant Physiol. 167, 862–869. doi: 10.1016/j.jplph.2010.01.018

PubMed Abstract | Crossref Full Text | Google Scholar

Rundel, P. W. (1978). The ecological role of secondary lichen substances. Biochem. Syst. Ecol. 6, 157–170.

Google Scholar

Rundel, P. W., Dillon, M. O., Palma, B., Mooney, H. A., and Gulmon, S. L. (1991). The phytogeography and ecology of the coastal Atacama and Peruvian deserts. Aliso 13, 1–49.

Google Scholar

Rütting, T., Boeckx, P., Müller, C., and Klemedtsson, L. (2011). Assessment of the importance of dissimilatory nitrate reduction to ammonium for the terrestrial nitrogen cycle. Biogeosciences. 8, 1779–1791. doi: 10.5194/bg-8-1779-2011

Crossref Full Text | Google Scholar

Sánchez-Martín, R., Montesinos-Navarro, A., Ochoterena, H., Pisanty, I., Rodríguez-Sánchez, M., Verdú, M., et al. (2024). Homogeneous microenvironmental conditions under nurses promote facilitation. Funct. Ecol. 38, 350–362. doi: 10.1111/1365-2435.14486

Crossref Full Text | Google Scholar

Santander, C., García, S., Moreira, J., Aponte, H., Araneda, P., Olave, J., et al. (2021). Arbuscular mycorrhizal fungal abundance in elevation belts of the hyperarid Atacama Desert. Fungal Ecol. 51:101060. doi: 10.1016/j.funeco.2021.101060

Crossref Full Text | Google Scholar

Schulze-Makuch, D., Wagner, D., Kounaves, S. P., Mangelsdorf, K., Devine, K. G., De Vera, J. P., et al. (2018). Transitory microbial habitat in the hyperarid Atacama Desert. Proc. Natl. Acad. Sci. USA 115, 2670–2675. doi: 10.1073/pnas.1714341115

PubMed Abstract | Crossref Full Text | Google Scholar

Semenova, G. A., Fomina, I. R., and Biel, K. Y. (2010). Structural features of the salt glands of the leaf of Distichlis spicata “Yensen 4a” (Poaceae). Protoplasma 240, 75–82. doi: 10.1007/s00709-009-0092-1

Crossref Full Text | Google Scholar

Stotz, G. C., Salgado-Luarte, C., Vigil, A. T., De La Cruz, H. J., Pastén-Marambio, V., and Gianoli, E. (2021). Habitat-islands in the coastal Atacama Desert: loss of functional redundancy, but not of functional diversity, with decreased precipitation. Ann. Bot. 127, 669–680. doi: 10.1093/aob/mcaa206

PubMed Abstract | Crossref Full Text | Google Scholar

Team, R. (2023). RStudio: Integrated Development for R version 4.2. 3 [Software]. RStudio Team.

Google Scholar

Thomas, G. H., Southworth, T., León-Kempis, M. R., and Leech, A. (2006). Novel ligands for the extracellular solute receptors of two bacterial TRAP transporters. Microbiology 152, 187–198. doi: 10.1099/mic.0.28334-0

PubMed Abstract | Crossref Full Text | Google Scholar

Towett, E. K., Shepherd, K. D., and Cadisch, G. (2013). Quantification of total element concentrations in soils using total X-ray fluorescence spectroscopy (TXRF). Sci. Total Environ. 463-464, 374–388. doi: 10.1016/j.scitotenv.2013.05.068

PubMed Abstract | Crossref Full Text | Google Scholar

Vaseva, I. I., Simova-Stoilova, L., Kirova, E., Mishev, K., Depaepe, T., Van Der Straeten, D., et al. (2021). Ethylene signaling in salt-stressed Arabidopsis thaliana ein2-1 and ctr1-1 mutants – a dissection of molecular mechanisms involved in acclimation. Plant Physiol. Biochem. 167, 999–1010. doi: 10.1016/j.plaphy.2021.09.029

PubMed Abstract | Crossref Full Text | Google Scholar

Vásquez-Dean, J., Maza, F., Morel, I., Pulgar, R., and González, M. (2020). Microbial communities from arid environments on a global scale. A systematic review. Biol. Res. 53:29. doi: 10.1186/s40659-020-00296-1

PubMed Abstract | Crossref Full Text | Google Scholar

Voigt, C., Klipsch, S., Herwartz, D., Chong, G., and Staubwasser, M. (2020). The spatial distribution of soluble salts in the surface soil of the Atacama Desert and their relationship to hyperaridity. Glob. Planet. Chang. 184:103077. doi: 10.1016/j.gloplacha.2019.103077

Crossref Full Text | Google Scholar

Vurukonda, S. S. K. P., Vardharajula, S., Shrivastava, M., and SkZ, A. (2016). Enhancement of drought stress tolerance in crops by plant growth promoting rhizobacteria. Microbiol. Res. 184, 13–24. doi: 10.1016/j.micres.2015.12.003

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Chen, J., Zhou, S., Wang, X., Chen, Y., Lin, X., et al. (2017). 16S rRNA gene high-throughput sequencing reveals shift in nitrogen conversion related microorganisms in a CANON system in response to salt stress. Chem. Eng. J. 317, 512–521. doi: 10.1016/j.cej.2017.02.096

Crossref Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07

PubMed Abstract | Crossref Full Text | Google Scholar

Warren-Rhodes, K. A., Lee, K. C., Archer, S. D. J., Cabrol, N., Ng-Boyle, L., Wettergreen, D., et al. (2019). Subsurface microbial habitats in an extreme desert Mars-analog environment. Front. Microbiol. 10:69. doi: 10.3389/fmicb.2019.00069

PubMed Abstract | Crossref Full Text | Google Scholar

Warren-Rhodes, K. A., Rhodes, K. L., Pointing, S. B., Ewing, S. A., Lacap, D. C., Gómez-Silva, B., et al. (2006). Hypolithic cyanobacteria, dry limit of photosynthesis, and microbial ecology in the hyperarid Atacama Desert. Microb. Ecol. 52, 389–398. doi: 10.1007/s00248-006-9055-7

PubMed Abstract | Crossref Full Text | Google Scholar

Weather and Climate. (2022). Historical weather in Antofagasta, Chile. Available online at: https://www.weatherandclimate.eu/85442_3 (Accessed October, 2022).

Google Scholar

Wright, E. S. (2016). Using DECIPHER v2.0 to analyze big biological sequence data in R. R J. 8, 352–359. doi: 10.32614/rj-2016-025

Crossref Full Text | Google Scholar

Wu, Q., Ren, H., Bisseling, T., Chang, S. X., Wang, Z., Li, Y., et al. (2021). Long-term warming and nitrogen addition have contrasting effects on ecosystem carbon exchange in a desert steppe. Environ. Sci. Technol. 55, 7256–7265. doi: 10.1021/acs.est.0c06526

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, F., Leng, B., and Wang, B. (2016). Progress in studying salt secretion from the salt glands in recretohalophytes: how do plants secrete salt?. Frontiers in Plant Science. 7:977. doi: 10.3389/fpls.2016.00977

Crossref Full Text | Google Scholar

Zhang, X., Ward, B. B., and Sigman, D. M. (2020). Global nitrogen cycle: critical enzymes, organisms, and processes for nitrogen budgets and dynamics. Chem. Rev. 120, 5308–5351. doi: 10.1021/acs.chemrev.9b00613

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, S., Zhang, Q., Liu, M., Zhou, H., Ma, C., and Wang, P. (2021). Regulation of plant responses to salt stress. Int. J. Mol. Sci. 22:4609. doi: 10.3390/ijms22094609

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: rhizosphere, Suaeda foliosa , Distichlis spicata , microbial community, hyper-arid core

Citation: Fortt J, ​Castro-Severyn J, Choque A, Donoso G, Stoll A, Jones DL, Saavedra CP, Fuentes B and Remonsellez F (2025) Plant species-specific rhizobiome assembly in the hyper-arid Atacama Desert. Front. Microbiol. 16:1587491. doi: 10.3389/fmicb.2025.1587491

Received: 04 March 2025; Accepted: 12 August 2025;
Published: 16 September 2025.

Edited by:

Ondrej Uhlik, University of Chemistry and Technology in Prague, Czechia

Reviewed by:

Sakineh Abbasi, Institut National de recherche pour l’agriculture, l’alimentation et l’environnement (INRAE), France
Olli H. Tuovinen, The Ohio State University, United States
Stacey J. Doherty, Cold Regions Research and Engineering Laboratory, United States

Copyright © 2025 Fortt, Castro-Severyn, Choque, Donoso, Stoll, Jones, Saavedra, Fuentes and Remonsellez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: F. Remonsellez, ZnJlbW9uc2VAdWNuLmNs

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.