- 1Laboratory of Soil Microbial Ecology and Biogeochemistry (LEMiBiS), Institute of Agri-Food, Animal and Environmental Sciences (ICA3), Universidad de O’Higgins, San Fernando, Chile
- 2Center of Applied Ecology and Sustainability (CAPES), Santiago, Chile
- 3Sustainable Agroecosystems, Department of Environmental Systems Science, ETH Zurich, Zurich, Switzerland
- 4Institute of Biology, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile
Forest fires alter soil microbial communities that are essential to support ecosystem recovery following land burning. These alterations have different responses according to soil abiotic pre- and post-fire conditions and fire severity, among others, and tend to decrease along vegetation recovery over time. Thus, understanding the effects of fires on microbial soil communities is critical to evaluate ecosystem resilience and restoration strategies in fire-prone ecosystems. We studied the state of community-level physiological profiles (CLPPs) and the prokaryotic community structure of rhizosphere and bulk soils from two fire-affected sclerophyll forests (one surveyed 17 months and the other 33 months after fire occurrence) in the Mediterranean climate zone of central Chile. Increases in catabolic activity (by average well color development of CLPPs), especially in the rhizosphere as compared with the bulk soil, were observed in the most recently affected site only. Legacy of land burning was still clearly shaping soil prokaryote community structure, as shown by quantitative PCR (qPCR) and Illumina MiSeq sequencing of the V4 region of the 16S rRNA gene, particularly in the most recent fire-affected site. The qPCR copy numbers and alpha diversity indexes (Shannon and Pielou’s evenness) of sequencing data decreased in burned soils at both locations. Beta diversity analyses showed dissimilarity of prokaryote communities at both study sites according to fire occurrence, and NO3– was the common variable explaining community changes for both of them. Acidobacteria and Rokubacteria phyla significantly decreased in burned soils at both locations, while Firmicutes and Actinobacteria increased. These findings provide a better understanding of the resilience of soil prokaryote communities and their physiological conditions in Mediterranean forests of central Chile following different time periods after fire, conditions that likely influence the ecological processes taking place during recovery of fire-affected ecosystems.
Introduction
Fires are key events regulating the structure and function of ecosystems (Yang et al., 2020); they have been important drivers of vegetation evolution and adaptation in most of the Mediterranean biomes across the globe. However, in the case of the Mediterranean ecosystem of central Chile, as opposed to the other four biomes of this type across the globe, land burning has not been considered an evolutionary force to select for fire-prone vegetation (Rundel et al., 2016, 2018). Fires have also represented an important pressure inducing land degradation worldwide, especially in ecosystems of warm and dry summers (Esteves et al., 2012), which are highly susceptible to the effects of climate change (Fletcher and Zielhofer, 2013). This is the case of Mediterranean ecosystems, which, despite occupying only about 2.2% of the land surface worldwide (Rundel et al., 2018), are considered biodiversity hotspots constantly threatened by anthropogenic disturbances (García-Vega and Newbold, 2020). This is of particular concern in Chile, where fire events are expected to increase due to extended megadroughts synchronized with human interventions (Fuentes-Castillo et al., 2019).
In forests, land burning not only exert an individual impact on each of their biotic and abiotic components but also on the relationship among these constituents, compromising the functionality of the whole ecosystem (Brown and Smith, 2000). Land burning can alter forest belowground conditions that are essential to support aboveground life (Certini, 2005). Several studies have reported a greater sensitivity of soil biological properties as compared with physicochemical properties to fire occurrence (Hart et al., 2005; Mataix-Solera et al., 2009). For microbial indicators, phylogenetic community structure has been shown to evidence greater sensitivity to land burning than those targeting soil ecosystem functions such as microbial biomass, respiration, and specific enzyme activities related to C, N, and P cycles (Pérez-Valera et al., 2019). The impact on soil microbial status is of particular importance since microbe-mediated soil processes are one of the most important drivers of ecosystem recovery in fire-affected lands (Neary et al., 1999). Forest fires can threaten soil organisms due to their different sensitivity to soil heating (Hart et al., 2005; Mataix-Solera et al., 2009; Rodríguez et al., 2018). In addition, indirect effects of fires occur by altering abiotic factors known to shape microbial community structure, such as soil pH (Lauber et al., 2009; Rousk et al., 2010), soil organic matter (SOM) and soil organic carbon (SOC) (Eilers et al., 2010), soil oxygen (e.g., aeration) (Fierer, 2017), and nutritional conditions (Lauber et al., 2008), all of which are affected by fires (Hart et al., 2005). Direct and indirect forest fire effects on soil microbiota can reduce the total number of microbial species (Certini, 2005), increase the proportion of bacteria over fungi (Fultz et al., 2016), and change the abundance of particular taxonomic groups (Pressler et al., 2019). Moreover, several studies have shown negative effects of land burning on soil microbial activity and functional diversity (Mataix-Solera et al., 2009; Wang et al., 2016; Alcañiz et al., 2018; Sadeghifar et al., 2020).
Shifts in prokaryote community composition following fires have been reported in Mediterranean ecosystems, which have been related to local environmental conditions prior to land burning (Goberna et al., 2012). Soil prokaryote communities play key soil functions, including soil weathering, primary production, and organic matter decomposition, all of which are essential in the ecological functioning of soils and susceptible to fires (Pérez-Valera et al., 2019). For example, increases in Proteobacteria and Firmicutes, and a decrease in Acidobacteria in Mediterranean burned soils after 2 and 3 years of fire have been previously reported (Rodríguez et al., 2018). Proteobacteria are recognized by their ability to respond to labile C sources particularly in the rhizosphere, while Acidobacteria have the ability to degrade cellulose and lignin (Lagos et al., 2015). Therefore, changes in soil bacteria composition induced by fires and alterations of soil abiotic properties will have a direct effect of soil ecosystem functions (Pérez-Valera et al., 2019). Although several studies have reported effects of forest fires on microbial composition at phylum level, the evaluation of such effects on Mediterranean ecosystems of central Chile is rare, especially when considering potential relationships with soil functionality and physicochemical properties after fire.
The degree of alterations in soil microbial community structure and functional diversity after fires greatly depends on soil physicochemical characteristics prior and after land burning, fire severity, vegetation type, and other site-specific conditions (Certini, 2005; Hart et al., 2005; Mataix-Solera et al., 2009). The effect of forest fires on soil microbial conditions also depends on time after fires, as immediate to long-lasting effects have been reported (Mataix-Solera et al., 2009; Singh et al., 2017). Fire perturbance to belowground conditions is usually expected to last until vegetation recovers over ecological successions (Hart et al., 2005). Soil microbial activity following fires in Mediterranean forests have been shown to increase within the first year of fire occurrence and then reach values comparable with unburned conditions after 32 months of land burning (Bárcenas-Moreno et al., 2011). Other studies have evidenced different responses of microbial activity indicators following fires, which mainly remain altered within 1 year after fires to later resemble unburned conditions over a 3-year time frame (Sadeghifar et al., 2020). In boreal forest ecosystems, bacterial community composition and diversity have also been shown to be affected following 1 year of land burning and then recovered after 11 years (Xiang et al., 2014). As the literature suggests, ecosystems affected by fires take varying time spans to recover after land burning. Nonetheless, it is known that changes in soil microbial conditions and ecosystem functioning rapidly occur within the earliest stages following fires; thus, studies focusing on post-fire dynamics in the scale of months to years are of particular interest to understand recovery in relation to plant succession and changes in soil properties (Knelman et al., 2017; Pérez-Valera et al., 2019). This is particularly true for understudied environments that are less represented in the literature where local evidence is required, as it is the case for fire-affected soils of the Mediterranean ecosystem of central Chile.
The present work aimed to evaluate the state of community-level physiological profile and prokaryotic community structure of soil compartments (bulk and rhizosphere soils) at two fire-affected sclerophyll forests in the Mediterranean zone of central Chile. These sites had similar arboreous plant composition but contrasting land location and time elapsed after fire (17 months and the other 33 months after fire occurrence). We hypothesized that fire occurrence affects community-level physiological profiles and prokaryote community structure to a different extent depending on the time elapsed after fires, with the least recent fire affected soils closely resembling unburned conditions. To accomplish this, we used microbial community-level physiological profiles to evaluate carbon source utilization patterns, molecular approaches such as qPCR, and high-throughput sequencing of 16S rRNA genes to assess abundance and structure of the prokaryotic community, and classical soil physicochemical analyses to describe microbial habitats.
Materials and Methods
Research Sites
This study was conducted within the Mediterranean climate zone of Central Chile, specifically in the O’Higgins administrative region, at two locations affected by one forest fire event, either in January of 2017 or May of 2018 (Figure 1). The first research site (affected by a low-severity fire in January of 2017; García-Carmona et al., 2021) is in the Pumanque commune (34°35′44.99″ S; 71°42′17.349″ W) toward the east front of the coastal mountain range at an elevation of 100 m. At this place, a mean monthly precipitation of 36.7 mm is registered for the period 2000–2017 (Center of Climate and Resilience Research, n.d.). Soils are classified as Aquic Dystric Xerochrepts (Inceptisol), which originated from alluvium-colluvial parent materials, exhibited mainly loam to loamy sand textures, are stratified with 50–120 cm of depth, and have a mean slope of 15% (CIREN, 1996). The second site [affected by a low-severity fire in January of 2018, according to visual inspections recommended by Keeley (2009)] is in the Requínoa commune (34°14′19.314″ W, 70°41′26.52″ S) on the Piedmont of the west side of the Andes Mountain ridges at an elevation of 1,000 m. Here, a mean monthly precipitation of 25.0 mm was registered for the period 2000–2017 (Center of Climate and Resilience Research, n.d.). Soils at this place are classified as Mollic Haploxeralfs (Alfisol); they were developed on volcanic material, and exhibit loamy sand textures and shallow depths (20–80 cm depth), with a mean slope of 17% (CIREN, 1996). Vegetation at both sites is dominated by typical woody native species of the sclerophyll forest of Central Chile, including the tree species Peumus boldus, Lithraea caustica, and Quillaja saponaria and shrubs such as Trevoa trinervis (Supplementary Figure 1A). In addition, the herbaceous vegetation is composed by approximately 18 species in Requínoa and 30 species in Pumanque, of which six are common for both sites [Erodium sp., Loasa triloba, Bromus hordeaceus, Hypochaeris radicata, Anthriscus caucalis, Oxalis micrantha (Supplementary Figure 1B)].
Figure 1. Sampling sites consisting of sclerophyll forests located within the Mediterranean climate zone of Central Chile, in the O’Higgins administrative region. White dots represent the study locations, which were affected by forest fires in January of 2017 (Pumanque) and May of 2018 (Requínoa). Image created using QGIS 3.12.3.
Experimental Design and Soil Sampling
In October of 2019, study sites of ∼0.7 and ∼1 ha in Pumanque and Requínoa, respectively, were chosen for this study. Within these, burned and unburned (undisturbed) areas were selected (Supplementary Figure 2), which represents the main variation source (treatment) in this study. Thus, Pumanque was surveyed 33 months following land burning while Requínoa was evaluated 17 months after fire occurrence. At each burned and unburned areas, three experimental plots of 20 × 20 m were marked (total = 12), and these were separated from the center of each other by 40–60 m, where the closest plots between burned and unburned areas were approximately 10 m apart from their border (Supplementary Figure 2). Inside each experimental plot, a 5 × 5 m grid was projected to obtain nine equidistant points from which soil samples were obtained (Supplementary Figure 3). At each sampling point, a radial buffer area of 1 m was selected and shovel excavated (up to rooting deep, 5–10 cm) to aseptically obtain (1) rhizosphere soil samples by hand shaking and collecting soil attached to roots of all herbaceous species described for each site (Supplementary Figure 1), and (2) bulk soil samples from soil not attached to plant roots. Rhizosphere and bulk soil samples (soil compartments) from each of the nine sampling points per plot were pooled to obtain 12 composite samples for bulk soil (∼2 kg each) and 12 composite samples for rhizosphere (∼300 g each) in total. Thus, within each study location (i.e., Pumanque and Requínoa), the studied variation sources were (1) treatment (burned and unburned soils) and (2) soil compartment (bulk soil and rhizosphere). For every composite soil sample, an aliquot was kept at 4°C for gravimetric water content and physiological profiling, a second portion was stored at −80°C until DNA extraction and further molecular analyses, and the rest was air-dried at room temperature for physicochemical analyses.
Physicochemical Analysis
Soil physicochemical properties, except for gravimetric water content (GWC), were measured only for bulk soil samples due to the small amount of rhizosphere samples obtained. The GWC was determined by oven-drying 10 g of samples at 105°C for 48 h. Soil pH and electrical conductivity (EC) were determined in a 1:2.5 and a 1:5 (w/v) water extract, respectively (Sadzawka et al., 2006). Nitrate (NO3–) and ammonium (NH4+) contents were determined by a distillation-titration method using a 2M potassium chloride (KCl) extraction solution (Sadzawka et al., 2006). Organic matter (OM) was assessed by acid oxidation and colorimetric measurement (Sadzawka et al., 2006). Total carbon (C) and nitrogen (N) contents were determined by the Dumas dry combustion method using a LECO TRUSPEC analyzer (Leco Corporation, MI, United States). Total elements (P, K, Fe, Mg, S, Mo, Na, Ca) were obtained from acid digestion and later determined by inductively coupled plasma atomic emission spectrometry (ICP-AES) (Agilent Technologies, Victoria, Australia) (U.S. EPA, 2007). Soil bulk density (BD) was determined by the clod method (Baver et al., 1973); clay, silt, and sand contents (soil texture) were determined by the soil hydrometer method (Bouyoucos, 1962); and aggregate stability (AS) was determined by wet sieving (Kemper and Rosenau, 1986), using a wet sieving apparatus (Eijkelkamp Soil & Water, Giesbeek, Netherlands).
Microbial Community-Level Physiological Profile (CLPP)
Carbon source utilization patterns of soil microbial communities, also known as community-level physiological profiles (CLPPs) (Garland and Mills, 1991), were assessed using Biolog EcoPlates (Biolog Inc., Hayward, CA, United States) for all samples (i.e., including both bulk soil and rhizosphere samples). Every plate included the following 31 different carbon sources: (1) carbohydrates: α-D-lactose, i-erythritol, D-xylose, β-methyl-D-glucoside, D-cellobiose, D-mannitol, N-acetyl-D-glucosamine; (2) carboxylic acids: 2-hydroxy-benzoic acid, 4-hydroxy-benzoic acid, D-glucosaminic acid, pyruvic acid methyl ester, α-ketobutyric acid, D-galactonic acid γ-lactone, D-malic acid, itaconic acid, γ-hydroxy-butyric acid, D-galacturonic acid; (3) amino acids: L-threonine, L-phenylalanine, glycyl-L-glutamic acid, L-asparagine, L-arginine, L-serine; (4) polymers: α-cyclodextrin, glycogen, tween 80, tween 40; (5) miscellaneous: D,L-α-glycerol-phosphate, glucose-1-phosphate; and (6) amines/amides: phenylethylamine, putrescine. In addition, a blank well was considered in triplicates. Each Biolog EcoPlate well was inoculated with a 150-μl aliquot of a 10–3 soil dilution obtained from 5 g of soil suspended in 45 ml of sterile saline solution (0.85% NaCl). The plates were incubated at 25°C and absorbance at 590 nm was determined after 24, 48, 72, 96, 120, 144, and 168 h using the Infinite 200 PRO NanoQuant Microplate Reader (Tecan Group Ltd., Männedorf, Switzerland). At each time interval, optical density (OD) values were corrected by subtracting the control (blank well) values from each plate well. The OD values at 120 h for each C-substrate by sample were averaged to determine microbial catabolic activity, this time corresponded to the mid-exponential growth phase of the average well color development (AWCD-CLPP) (Garland, 1996; Sofo and Ricciuti, 2019).
Soil DNA Extraction and Molecular Analyses
Soil DNA was isolated from 0.25 g of fresh soil from all samples using the DNeasy PowerSoil DNA isolation kit (QIAGEN, Valencia, CA, United States) following the manufacturer’s instructions. Prokaryotes (bacteria and archaea) were surveyed by quantitative PCR (qPCR) assay and bar-coded amplicon sequencing using the primers 515F (5′ GTG YCA GCM GCC GCG GTA A 3′) (Parada et al., 2016) and 806R (5′ GGA CTA CNV GGG TWT CTA AT 3′) (Apprill et al., 2015) complementary to the V4 region of the 16S rRNA genes. For qPCR, a TaqMan probe (16S probe 5′ TGT AGC RGT GAA ATK CGT AG 3′) was designed for the conserved sequence region amplified by 515F and 806R, by which the dye 6-carboxy-fluorescein (6-FAM) and the Black Hole Quencher (BHQ1) were attached to their 5′ and 3′ ends, respectively. The qPCRs were carried out on a StepOne Plus Real Time PCR System (Thermo Fisher Scientific, United States), in a final volume of 20 μl, containing 10 μl of NZY qPCR Probe Master Mix ROX plus (NZYTech), 0.9 μM of the probe, 0.4 μM of the amplification primers, 2 μl of template DNA, and ultrapure water up to 20 μl. The reaction mixture was incubated as follows: an initial incubation at 95°C for 10 min, followed by 40 cycles of denaturation at 95°C for 15 s, hybridization at 50°C for 1 min, extension at 60°C for 1 min, and a final extension step at 60°C for 30 s. The intensity of the fluorescence emitted by the probe at each cycle of the PCR reaction was registered by using the software StepOne (Applied Biosystems, Foster City, CA, United States), which allowed the estimation of the quantification cycle (Cq). Standard for qPCR was prepared from a known unburned sample (ranging from 5 × 1010 to 5 × 106 16S copy numbers). This standard was obtained from the amplification of the 16S rRNA region complementary to the 515F/806R primer set, which was carried out on a Bio-Rad T100 Thermal Cycler (Bio-Rad, United States), in a final volume of 25 μl, containing 2.5 μl of template DNA, 0.5 μM of the primers, 12.5 μl of Supreme NZYTaq 2 × Green Master Mix (NZTech), and ultrapure water up to 25 μl. The PCR reaction mixture was incubated as follows: an initial denaturation at 95°C for 5 min, followed by 25 cycles of 95°C for 30 s, 50°C for 45 s, 72°C for 45 s, and a final extension step at 72°C for 10 min. The PCR product was purified and used to generate the standard curve in the qPCR experiment. The standard curve was performed by a five-fold dilution series of 16S copy numbers from the qPCR standard by plotting the Cq values against the total 16S rRNA gene copy numbers (Jiang, 2018). Efficiency and R2 values for the 16S quantification were 84.65 and 0.99%, respectively.
Bar-coded amplicon sequencing was performed on an Illumina MiSeq PE300 platform (Illumina, San Diego, CA, United States) at the AllGenetics & Biology SL laboratory service facility (A Coruña, Spain). PCRs were carried out in a final volume of 25 μl, containing 2.5 μl of template DNA, 0.5 μM of the primers, 12.5 μl of Supreme NZYTaq 2 × Green Master Mix (NZYTech, Lisbon, Portugal), and ultrapure water up to 25 μl. The reaction mixture was incubated as follows: an initial denaturation at 95°C for 5 min, followed by 25 cycles of 95°C for 30 s, 50°C for 45 s, 72°C for 45 s, and a final extension step at 72°C for 10 min. The oligonucleotide indices that are required for multiplexing different libraries in the same sequencing pool were attached in a second PCR round with identical conditions but only 5 cycles and 60°C as the annealing temperature. No template controls were included to check for cross-contamination during library preparation. The libraries were run on 2% agarose gels stained with GreenSafe (NZYTech, Lisbon, Portugal) and imaged under UV light to verify the library size. Libraries were purified using the Mag-Bind RXNPure Plus magnetic beads (Omega Bio-tek, Norcross, GA, United States), following the instructions provided by the manufacturer. Then, they were pooled in equimolar amounts according to the quantification data provided by the Qubit dsDNA HS Assay (Thermo Fisher Scientific, Waltham, MA, United States).
Processing and Analysis of Sequencing Data
Sequencing data were processed with a customized pipeline largely based on “VSEARCH” (Rognes et al., 2016). Illumina MiSeq reads can be found in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under accession number PRJNA784510. PhiX reads were removed with Bowtie2 (Langmead and Salzberg, 2012) and primer sequences were trimmed with Cutadapt (Martin, 2011). Paired-end reads were merged by using the fastq_mergepairs algorithm in “VSEARCH” and quality filtered using a maximum expected error of one using the fastq_filter algorithm in “VSEARCH” (Edgar and Flyvbjerg, 2015). Merged and quality filtered sequences were dereplicated using derep_fulllength algorithm in “VSEARCH” and the UNOISE algorithm implemented in “VSEARCH” (Edgar, 2016) was used to retrieve amplicon sequence variants (ASVs). The uchime algorithm implemented in “VSEARCH” (Edgar et al., 2011) was applied to detect and remove chimeric ASVs. Remaining ASVs were verified by using Metaxa2 to evaluate the presence of non-target sequences for the 16S rRNA gene (Bengtsson-Palme et al., 2015). The final ASV table was obtained by mapping the high-quality sequences against the verified ASV centroids using the usearch_global algorithm implemented in “VSEARCH”. Taxonomic classification was performed by using the Naïve Bayesian classifier Sintax implemented in “VSEARCH” (Edgar, 2016) to classify the sequences against the non-redundant SILVA v132 database (Quast et al., 2012). ASVs assigned to mitochondria, chloroplasts, and eukaryotes were removed. Rare ASVs with an abundance < 0.005%, sparse ASVs that occur in less than four samples (which represent 1/6 of the sampling effort), and ASV outliers for which the greatest abundance was more than 0.15 times in a sample compared with the second most abundant value across the other samples were removed from the ASV table. Finally, the filtered ASV table was normalized by iterative subsampling with 100 iterations using the rrarefy function implemented in the “vegan” package in R (Oksanen et al., 2020).
Statistical Analysis
Statistical analysis were performed in R statistics version 3.5.1 (R Core Team, 2021). Treatments (burned vs. unburned) and soil compartments (bulk soil and rhizosphere) were considered as the variation sources for statistical analyses, while research location (i.e., sites) was not considered due to the natural differences in biophysical characteristics described previously. The analysis for physicochemical properties of bulk soil samples were performed following normality and homoscedasticity confirmation. To assess differences within treatments (i.e., burned and unburned conditions), a Student’s t-test was applied for each of the variables analyzed. For CLPP data sets, a two-way ANOVA considering treatment and soil compartment (i.e., bulk and rhizosphere soil) was performed after normality and homoscedasticity were confirmed. These datasets included the AWCD-CLPP values and the diversity of substrate utilization, determined by the Shannon diversity index (Shannon-CLPP) by: H′ = −Σ pi × ln pi (Shannon, 1948), where pi is the ratio of the respiration rate of every single C-substrate to the sum of all substrates, both obtained from OD values at 120 h of incubation. A heatmap with hierarchical clustering (based on Euclidean distance) was performed for the 31 different carbon sources, considering the average of each C-substrates per treatment and compartment, with the function heatmap from the package “stats” (R Core Team, 2021). The rarefied sequencing data were used to determine α-diversity indexes (ASVs, as observed richness, Shannon diversity index, and Pielou’s evenness) with the diversity function in “vegan.” A Spearman correlation analysis between soil physicochemical, AWCD-CLPP, Shannon-CLPP, and microbial α-diversity indexes was applied with the function cor from the package “corrplot” (Wei and Simko, 2021). To assess differences within treatments and soil compartment for each of these indexes, a two-way ANOVA was conducted after normality and homoscedasticity were confirmed. Differences in β-diversity were assessed by calculating Bray–Curtis dissimilarities based on the iteratively rarefied ASV abundance tables with a PERMANOVA and PERMDISP by using the functions adonis2 and betadisper from the package “vegan”, respectively. This allowed the detection of factors with significant effects on prokaryotic community structure. Thus, a canonical analysis of principal coordinates (CAP, Anderson and Willis, 2003) was applied with the function CAPdiscrim from the package “BiodiversityR” (Kindt and Coe, 2005). A distance-based redundancy analysis (dbRDA, Legendre and Andersson, 1999) was used to constrain β-diversity patterns by the bulk soil physicochemical properties using the function dbrda and envfit from the package “vegan.” In addition, relative abundance data of ASVs was Z-transformed to show relative change between burned and unburned soils. A univariate PERMANOVA (Anderson, 2001) and PERMDISP (Anderson, 2001) were applied on Z-transformed relative abundance of ASVs at the phylum level and also at the order taxa and genus taxa levels to evaluate treatment effects on each soil compartment with the adonis2 and betadisper function from the package “vegan” (Oksanen et al., 2020).
Results
Soil Abiotic Properties
Statistically significant differences between burned and unburned bulk soils were observed for NH4+ and NO3– in Pumanque (taken 33 months after fire occurrence) and for P, Fe, Na, and AS in Requínoa (taken 17 months after fire occurrence) (Supplementary Table 1). In Pumanque, NH4+ contents almost doubled in burned soils, and NO3– contents decreased by approximately two thirds compared with unburned soils (Supplementary Table 1). In Requínoa, total P and Fe contents dropped by 43 and 14% in burned soils, respectively, and Na increased by almost 70% in burned soils. In terms of physical properties, AS was significantly lower in burned soils in Requínoa. Although non-significant, NH4+ and NO3– contents showed opposite results in Requínoa compared with Pumanque. In addition, also not significant, a decreasing trend in OM, C, and N and an increasing trend in soil pH and BD were observed for burned soils at both sites.
Microbial Community-Level Physiological Profile (CLPP)
Microbial catabolic activity and diversity of substrate utilization [here assessed by AWCD-CLPP and Shannon diversity index (Shannon-CLPP), respectively] responded differently to forest fires at both sites. In Pumanque (assessed 33 months after fire), there were no effects of burning on AWCD-CLPP or Shannon-CLPP (Table 1). In addition, no differences in carbon utilization were observed between bulk and rhizosphere samples at this site (Figure 2), although higher values for AWCD and Shannon were observed for rhizosphere samples in both burned and unburned conditions (Figure 2). On the contrary, in Requínoa (assessed burned 17 months after fire), soils showed significantly higher catabolic activity, with AWCD values in the burned soils being 37% higher than in the unburned soils (Table 1 and Figure 2). Significant differences were also observed between soil compartment (Table 1), with rhizosphere samples showing higher AWCD values than bulk soil samples in both burned and unburned conditions (Figure 2). Shannon diversity of substrate utilization was also significantly higher in the rhizosphere when compared with the bulk soil, whereas burning had no effect (Table 1 and Figure 2).
Table 1. F-values from two-way ANOVA for biological activity and α-diversity of prokaryotic communities, permutational multivariate analysis (PERMANOVA), and permutation test for homogeneity of multivariate dispersion (PERMDISP) for prokaryotic community (β-diversity) (the table is split by study locations).
Figure 2. Bar plots of the average well color development (AWCD) (mean ± SE, n = 3) and Shannon diversity index based on CLPP (Shannon-CLPP) (mean ± SE, n = 3) of microbial communities for each study location contrasting treatments (burned vs. unburned) and soil compartment (bulk soil vs. rhizosphere soil). Different uppercase letters show significant (p < 0.05) differences by treatment. Different lowercase letters show significant (p < 0.05) differences by compartment within each treatment. Interactions between treatment and compartment were not significant (p > 0.05).
Clustered heatmaps for the average (n = 3) of each C-substrate utilization showed different catabolic capacity for soil microbial communities in both study locations (Figure 3), where utilization decreased for recalcitrant substrates (e.g., polymeric compounds and carboxylic acids) compared with labile C sources (e.g., carbohydrates). Soil microbial communities from Requínoa showed higher C utilization capacity than those from Pumanque, with unclear trends regarding the type of C, while C sources in the form of carbohydrates and carboxylic acids were generally more used by soil microorganisms in Pumanque. A closed look to C source consumption indicated that L-asparagine (ASP) was the C source most consumed in soils from Pumanque, especially rhizosphere of burned soils. Similarly, ASP was the C source most consumed in soils from Requínoa particularly in rhizosphere of burned soils, in addition to galacturonic acid (GALACTU) and N-acetyl-D-glucosamine (NAGA), which were specific for this site.
Figure 3. Clustered heatmaps based on the average (n = 3) of carbon source utilization patterns of soil microbial communities, also referred as community-level physiological profile (CLPP) for (A) Pumanque and (B) Requínoa, by treatment (burned vs. unburned) and soil compartment (R = rhizosphere; BS = bulk soil). Row dendrograms show the Euclidean distance between treatments, while column dendrograms show the distance between the utilization of C sources. Yellow boxes represent higher microbial C utilization. The following 31 different carbon sources were used: (1) carbohydrates: α-D-lactose (LAC), i-erythritol (ERY), D-xylose (XYL), β-methyl-D-glucoside (BMDG), D-cellobiose (CELL), D-mannitol (MAN), N-acetyl-D-glucosamine (NAGA); (2) carboxylic acids: 2-hydroxy-benzoic acid (2HBA), 4-hydroxy-benzoic acid (4HBA), D-glucosaminic acid (GLUCA), pyruvic acid methyl ester (PAME), α-keto-butyric acid (KETO), D-galactonic acid γ-lactone (GALAC), D-malic acid (MAL), itaconic acid (ITA), γ-hydroxy-butyric acid (GHBA), D-galacturonic acid (GALACTU); (3) amino acids: L-threonine (THRE), L-phenylalanine (PHE), glycyl-L-glutamic acid (GLYC), L-asparagine (ASP), L-arginine (ARG), L-serine (SER); (4) polymers: α-cyclodextrin (CYCLO), glycogen (GLYC), Tween 80 (T80), Tween 40 (T40); (5) miscellaneous: D,L-α-glycerol-phosphate (GLYPHO), glucose-1-phosphate (GLUC); and (6) amines/amides: phenylethylamine (PHE), putrescine (PUT).
Prokaryotic Abundance, Alpha and Beta Diversity, and Taxonomy Assignments
DNA-based analyses were more sensitive to detect the effects of fire than the substrate utilization profiles (Table 1). Unburned soils in both Pumanque and Requínoa showed significantly higher 16S rRNA gene copy numbers when compared with the burned soils (Table 1 and Figure 4). No differences in gene copy numbers were observed between bulk and rhizosphere soils (Table 1 and Figure 4).
Figure 4. Bar plots of abundance (mean ± SE, n = 3) and alpha diversity (mean ± SE, n = 3) of prokaryotic communities for each study locations showing effect of treatment (burned vs. unburned) and soil compartment (bulk soil vs. rhizosphere soil). Different uppercase letters show significant (p < 0.05) differences by treatment. Different lowercase letters show significant (p < 0.05) differences by compartment. Interactions between treatment and compartment were not significant (p < 0.05).
A total of 3,782 curated prokaryote ASVs were obtained across all samples. Statistically significant differences in observed richness were detected only between burned and unburned samples in Requínoa (Table 1), where ASVs in burned soils decreased by 12% (2,687 ± 143.7) when compared with unburned soils (3,041 ± 229.1) (Figure 4). No differences in richness were observed between soil compartments. Shannon diversity (Shannon-16S) and Pielou’s evenness (Evenness-16S) were higher in unburned soils at both study locations, with alterations due to soil compartment only in Pumanque (Table 1 and Figure 4).
Beta diversity analyses showed differences in soil prokaryotic community structure according to fire occurrence at both study locations (Table 1 and Figure 5, and Supplementary Figure 6). Moreover, in Pumanque, communities also separated apart according to their compartment (bulk soil and rhizosphere), and larger dissimilarities were observed between rhizosphere of burned and unburned soils when compared with the bulk soils (Figure 5A). In Requínoa, prokaryotic community structures also showed dissimilarities between burned and unburned soils (Figure 5B), while differences between compartment were not conclusive due to low CAP reclassification success. On the other hand, Bray–Curtis distances between burned and unburned soils from Requínoa were approximately two times greater than Pumanque. Distance-based redundancy analysis (dbRDA) for bulk soil data sets showed discrimination of prokaryote communities by forest fire associated to some soil physicochemical properties. The most parsimonious and significant dbRDA models showed that soil pH and NO3– content explained prokaryotic community changes in Pumanque (Figure 5C), while contents of NO3–, and total contents of Fe, N, Ca, and NO3– were important explanatory variables of prokaryote community changes in Requínoa (Figure 5D).
Figure 5. Beta diversity for prokaryotic communities from burned and unburned bulk and rhizosphere soil samples. Canonical analysis of principal coordinates (CAP) based on Bray–Curtis dissimilarities constrained by burned and unburned as well as bulk and rhizosphere soils from Pumanque (A) and Requínoa (B). The CAP reclassification (in percent) for each treatment is shown next to each cluster. Distance-based redundancy analysis (dbRDA) based on Bray–Curtis dissimilarities constrained by soil physicochemical properties (pH, NO3–, and total contents of Fe, N, and Ca) in burned and unburned bulk soils from Pumanque (C) and Requínoa (D). Soil physicochemical properties (explanatory variables) are shown by solid eigenvectors. PERMANOVA was performed by each site independently (Table 1).
Taxonomy assignments resulted in a total of twenty-two and two classified bacterial and archaeal phyla, respectively. The most abundant phyla (>70%) across samples were Proteobacteria, Actinobacteria, and Acidobacteria (Supplementary Figure 5A). A total of sixteen and eight classified prokaryotic phyla were affected by fire occurrence in Pumanque and Requínoa, respectively. In Pumanque, forest fire negatively affected the relative abundance of Thaumarchaeota, Nitrospirae, Entotheonellaeota, Elusimicrobia, Acidobacteria, and Rokubacteria in both soil compartments, except for Rokubacteria with differences in rhizosphere samples only (Figure 6). Contrarily, relative abundance of the candidate division FBP, Gemmatimonadetes, Planctomycetes, Fibrobacteres, Hydrogenedentes, Armatimonadetes, Cyanobacteria, Firmicutes, Actinobacteria, and Patescibacteria increased in burned soils as compared with unburned soils, with similar responses in bulk soil and rhizosphere. However, relative abundance of Actinobacteria and Patescibacteria only increased in rhizosphere and bulk soil samples, respectively. In Requínoa, forest fire negatively affected the relative abundance of Rokubacteria, WS2, Planctomycetes, Acidobacteria, and Verrucomicrobia, which showed similar responses in bulk soil and rhizosphere. In contrast, Actinobacteria, Firmicutes, and Bacteroidetes increased in burned soils, with similar responses in both soil compartments.
Figure 6. Relative change (z-transformed) of prokaryote abundance, at the phylum level, between unburned (blue triangle) and burned (red spot) soils for Pumanque and Requínoa. The panel to the left shows the relative change in abundance from the overall mean ± SE (considering bulk soil and rhizosphere). The middle and right panels show the relative change in abundance for bulk and rhizosphere soils, respectively. Bulk soil and rhizosphere panels show the relative change of phylum between burned and unburned soil for each soil compartments. The significance of the PERMANOVA test is shown at left side of each relative change in the overall panel. Significance codes are based on p-values as follows: “***” 0.001, “**” 0.01, “*” 0.05, ns = not significant. PERMANOVA was not applied on bulk soil and rhizosphere due to scarce data.
From the classified phyla significantly affected by fires, there were also 33 and 26 orders showing significant changes due to fires in Pumanque and Requínoa, respectively (Supplementary Figure 7). In Pumanque, the relative abundance of 16 orders decreased due to fire occurrence, which mainly belonged to Actinobacteria (Supplementary Figure 7). Nitrososphaerales (Thaumarchaeota) was also among the taxa most affected by fire at this site, in both soil compartments. Instead, 17 orders increased in burned soils from Pumanque, which mainly belonged to Actinobacteria, Gemmatimonadetes, Planctomycetes, and Acidobacteria (Supplementary Figure 7). Among them, Hydrogenedentiales (Hydrogenedentes) showed the greatest increase in relative abundance due to fire occurrence. In Requínoa, 17 orders were negatively affected by fire and most of them belonged to Actinobacteria, Acidobacteria, and Planctomycetes (Supplementary Figure 7). Only ten orders increased in relative abundance after fire, which mainly belonged to Actinobacteria (Supplementary Figure 7). Fires reduced the relative abundance of Pyrinomonadales (Acidobacteria), Rokubacteriales (Rokubacteria), Solirubrobacterales (Actinobacteria), and Pseudonocardiales (Actinobacteria) at both locations. Contrarily, Bacillales (Firmicutes) was the only order that consistently increased in burned soils at both study locations (Supplementary Figure 7).
Discussion
Soil physicochemical properties did not show strong differences according to fire occurrence, especially in Pumanque where only clear alterations of NO3– and NH4+ were found in burned soils. In fire-affected soils, NH4+ is known to form as a result of organic matter mineralization induced by land burning, which is later transformed to NO3– by nitrification (Certini, 2005). However, this latter transformation could be differently altered after fire, as some findings have reported a decrease in nitrifying microorganisms after fires (Acea and Carballas, 1996), while others have reported higher nitrification rates occurring 12 months after fire (Covington and Sackett, 1992). Thus, lower nitrification rates in Pumanque, and the opposite trend for Requínoa, could partly explain differences in NO3– and NH4+ after forest fire, which is similar to those reported by Ball et al. (2010) who found higher nitrification rates in recent burned soil (12 years old) compared with a site without fire for 75 years. In addition, NH4+ and NO3– contents also respond to their different mobility in soils, where NO3– is prone to leaching over time while NH4+ could be retained in burned soils due to adsorption to organic matter (OM) or mineral surfaces (Mroz et al., 1980). A combination of these biotic and abiotic processes would likely relate to the negative correlations observed between NO3– and NH4+ contents in soils from Pumanque and Requínoa, where in the former NH4+ increased and NO3– decreased 33 months after land burning, and in the latter NH4+ decreased and NO3– increased after 17 months of forest fire occurrence. On the other hand, lower P and Fe contents in burned soils from Requínoa might be related to losses of Fe/P bound to organic matter (OM) in fire-affected soils (Norouzi and Ramezanpour, 2013), which might have migrated from surface to lower horizons (Mosugu et al., 1999). In addition, both study locations showed a slight increase in pH in burned soils, which was expected due to the OH– and cation release from the denaturation of OM occurring after fires (Certini, 2005; Alcañiz et al., 2018).
In burned soils, microbial activity can be promoted by the combustion of plant residues that increase C and N input to microorganisms (Díaz-Raviña et al., 1996; Liu et al., 2000), especially in the rhizosphere that is considered as a hotspot for microbial activity in soils (Prashar et al., 2014). Thus, negative effects of forest fires on plant cover and diversity indirectly affect soil microorganisms associated with roots (Certini, 2005; Hart et al., 2005). In our study, forest fires differently affected microbial catabolic activity assessed by Biolog EcoPlates. Influence of land burning was mostly evident in Requínoa, where forest fires increased the microbial catabolic activity, especially in the rhizosphere. Different responses of microbial metabolic fingerprints based on Biolog EcoPlates have been related to vegetation coverage after fires in Mediterranean pine forest (Moya et al., 2021). For instance, microbial catabolic activity has been shown to increase a few months after forest fires, especially under soils with high vegetation cover, while C substrate utilization patterns have also been observed to increase in burned soils, but with no relation to vegetation (Moya et al., 2021). However, opposite results have been reported in desert grassland ecosystems, where fires were shown to reduce soil microbial substrate utilization with variations observed across seasons (Liu et al., 2000). In our study, although herbaceous species richness in Pumanque (assessed 33 months after land burning) was greater than in Requínoa (assessed after 17 months following fires), the latter showed significant increase of substrate usage in burned soils. Thus, conditions in burned soils in Pumanque might closely resemble biological functionality of unburned soils at this site. Despite dissimilar responses of substrate utilization observed in soils from Pumanque and Requínoa after land burning, values were always higher for rhizosphere samples at both sites, supporting the influence of vegetation on this biological indicator (Zhang et al., 2020). Higher values associated to the rhizosphere using CLPP approaches (assessed either by Biolog EcoPlates or MicroResp) in other Mediterranean regions (de Armas-Ricard et al., 2016), and different ecosystems such as boreal coniferous forests (Hernesmaa et al., 2005) or agricultural fields (Sneha et al., 2021), have been attributed to root exudates that serve as C sources for microbial communities (Grayston et al., 1997; Rodríguez-Loinaz et al., 2008).
Carbon substrate utilization capacity has been shown to differ following land burning (Barreiro et al., 2015; Moya et al., 2021). Immediately after experimental fires, C sources in the form of carboxylic acids, amino acids, carbohydrates, and phenolic acids have been shown to greatly increase, while over a year all C source utilization tends to resemble unburned conditions (Barreiro et al., 2015). In our study, no clear patterns were observed for soils from Requínoa; however, soils from Pumanque showed higher utilization of carbohydrates and carboxylic acids, especially in rhizosphere samples of the burned soils. This can be likely due to the recovery of C labile forms reported to occur over a year after fire occurrence (Martín et al., 2009), which in the case of Pumanque could have been promoted by the recovery of herbaceous vegetation. The microbial L-asparagine (ASP) utilization increased in rhizosphere samples in burned soils from Pumanque and Requínoa (as well as other carbohydrates and organic acids in the case of the latter), which can be associated with the release of ammonia by microbial activity of asparaginases (Kandeler et al., 2011). However, as previously mentioned, NH4+ and NO3– contents showed opposite results at each site, which can likely relate to transformation of ammonia to nitrate differently affected at both sites. In the case of Pumanque, the consumption of ASP followed by low nitrification rates, as observed in other ecosystems few years after fire (Cobo-Díaz et al., 2015), can relate to higher ammonia detected at this site, while in Requínoa, ASP utilization followed by rapid microbial uptake and conversion to nitrate (Covington and Sackett, 1992; Hart et al., 2005; Switzer et al., 2012) likely relates to higher nitrate contents registered at this site. These findings support the hypothesis proposed by Yang et al. (2020), who stated that forest fires could stimulate functional microbial groups related to N cycling such as ureC, amoA (Ball et al., 2010), and nifH (Sun et al., 2016) genes.
In our study, the assessments of soil prokaryotic communities, based on molecular ribosomal markers, evidenced a greater sensitivity of this method as compared with CLPP to elucidate changes in response to the occurrence of fires. This might reflect different sensitivity of soil genetic and functional diversity after forest fires under the studied conditions. Although metabolic activity and the diversity of C substrate utilization were not affected by forest fires in Pumanque, burned soils at this site showed lower prokaryote abundance and α-diversity indexes, and clear dissimilarities in community composition between burned and unburned samples, as well as rhizosphere and bulk soils. These findings suggest the presence of functionally redundant microbial groups at this site, which coincides with other studies showing recovery of soil functionality but not bacterial community structure after fires (Pérez-Valera et al., 2019). Similarly, burned soils from Requínoa, which showed higher metabolic activity compared with unburned soils, had lower prokaryote 16S rRNA gene copy numbers and diversity, and clear distinct community structure. These results agree with several studies, including those performed in other Mediterranean regions, that have reported decreases in soil microbial abundance (e.g., direct counts, qPCR) and diversity (e.g., based on DNA data) due to negative effect of heating (Mataix-Solera et al., 2009; Alcañiz et al., 2018; Pressler et al., 2019; Marín and Rojas, 2020). Low prokaryote abundance and diversity in burned soils can be likely associated with initial detrimental effects of forest fires on soil microorganisms followed by slow recovery of genomic pools, which can be partly explained by persistent changes in soil physicochemical properties mentioned previously. Long-term effects of forest fires on soil microbial communities can also be associated with post-fire climate conditions, substrate availability, plant recolonization, changes in soil physicochemical properties, and time after fires, among others (Hart et al., 2005; Mataix-Solera et al., 2009; Switzer et al., 2012; Singh et al., 2017; Li et al., 2019). Based on our results, low 16S rRNA gene copy numbers and α-diversity in burned soils still evidence the legacy of fires despite the time elapsed after fires. This observation supports findings from other ecosystems across the globe, which suggest recovery of soil microbial communities beyond 10 years following fire occurrence (Xiang et al., 2014; Pressler et al., 2019), as opposed to those reporting conditions comparable with unburned counterpart close to 3 years after fires (Bárcenas-Moreno et al., 2011; Sadeghifar et al., 2020). In the case of our study, similar to Mediterranean ecosystems located at other latitudes, legacy of fires in soil microbial community structure and microbial activity is still detectable within an approximate 1- to 3-year period, which is particularly critical in the Chilean ecosystem with no fire-prone vegetation (Rundel et al., 2018, 2016).
Prokaryote communities were dissimilar according to fire occurrence at both study locations. Dissimilarities of prokaryote communities in Pumanque were also observed between bulk and rhizosphere soil samples. Thus, such differences probably resulted from greater vegetation richness that increased the input of organic matter with different quality to bulk soil in addition to C input from rhizosphere, by which these differences were not attributed to forest fires. On the other hand, greater distances between microbial community structures were found between burned and unburned soils from Requínoa compared with Pumanque, which might reflect the greater number prokaryotic genera that significantly changed after fire. In Pumanque, pH and the available NO3– content significantly explained variability in the community structure, while available NO3– and total Fe, N, and Ca contents were associated with prokaryote community changes in Requínoa. Soil pH and N has been shown to be important factors to determine the prokaryote community composition (Mandakovic et al., 2018; Li et al., 2019; Yang et al., 2020), which we confirm in this study. For instance, Li et al. (2019) reported that soil pH, total N, and available NH4+ contents were the main drivers of changes in microbial communities affected by fires. Similar to our study, Pérez-Valera et al. (2020) reported that soil nitrogen content and pH were strong predictors of soil microbial functions after forest fire in a Mediterranean region in Spain. Effects of soil pH on bacterial communities can be associated with (1) physiological pressure on soil bacteria, reducing the growth of some taxa unable to survive in some pH ranges, and (2) pH effects on soil properties, such as substrate availability (Mandakovic et al., 2018), which may drive changes in bacterial community composition (Lauber et al., 2009). On the other hand, NO3– could be associated with biological functioning described previously since it is immediately leached if not taken up (Li et al., 2019). This also suggests differences in composition of prokaryotic groups involved in N cycling. In this sense, decreasing total Fe content in burned soils from Requínoa might also result in lower abundance of prokaryotes able to synthesize Fe-dependent nitrogenases (Raymond et al., 2004).
In our study, some prokaryote taxa were consistently affected by fires at both sites, likely due to different sensitivity to forest fires reported for several microbial groups (Mataix-Solera et al., 2009). Firmicutes and Actinobacteria significantly increased in burned soils from both locations as also reported by Cobo-Díaz et al. (2015) in a fire-affected forest of the Mediterranean Basin, while Acidobacteria and Rokubacteria were reduced by fires at both research sites. Indeed, the former two phyla have been found to increase following fires due to their resistance to heat (Pérez-Valera et al., 2019). Moreover, Rodríguez et al. (2018) reported an increase in Firmicutes and a decrease in Acidobacteria in Spanish Mediterranean burned soils after 2 and 3 years of fire. Similarly, Li et al. (2019) showed higher abundance of Proteobacteria and Actinobacteria, and lower abundance of Acidobacteria, Verrucomicrobia, and Chloroflexi after wildfires. In our study, Firmicutes found to increase in burned soils at both locations belonged to the order of Bacillales, particularly to the genera Tumebacillus and Bacillus in the case of Pumanque (Supplementary Figure 9) and to Tumebacillus, Conhella, and Paenibacillus in the case of Requínoa (Supplementary Figure 9). Spore- or endospore-forming bacteria within all these genera have been previously described (Khianngam et al., 2010; Kim and Kim, 2016; Puri et al., 2016). This might explain the increase of relative abundance of Firmicutes in burned soils in our study, as such resistance structures allow to better cope with higher temperature and also promote proliferation of spore- or endospore-forming microbial groups due to spore germination enhanced by fire (Dworkin, 2006; Pérez-Valera et al., 2019). In the case of members of the Acidobacteria, which decreased in burned soils from Pumanque and Requínoa, these belonged to the Pyrinomonadales and DS-100 orders in Pumanque and to the Pyrinomonadales, Acidobacteriales, and Blastocatellales and 11–24 orders in Requínoa. Previous studies have shown that different groups of Acidobacteria have variable responses to changes in soil pH, nutritional conditions, and C source availability (Eichorst et al., 2011; Kielak et al., 2016). It has been shown that Acidobacteria subdivisions 1, 2, 3, 12, 13, and 15 negatively correlate with soil pH, while subdivisions 4, 6, 7, 10, 11, 11, 16, 17, 18, 22, and 25 show an opposite behavior (Kishimoto et al., 1991; Janssen, 2006; Sait et al., 2006; Jones et al., 2009; Kielak et al., 2016); even subdivision 6 can either increase or decrease according to soil pH (Kielak et al., 2016). Thus, increases in soil pH registered in our study as a result of land burning can likely explain the decrease of Acidobacteria susceptible to this condition (Li et al., 2020; Sun et al., 2021). Interestingly, Thaumarchaeota was the archaeal phylum most negatively affected by forest fires in Pumanque, particularly of the order Nitrososphaerales (Supplementary Figure 8). Species of this order can oxidize ammonia (Tourna et al., 2011; Brochier-Armanet et al., 2012; Bustamante et al., 2012), similar to species of Nitrospirae phylum such as those within the Nitrospirales order (Vijayan et al., 2021), which also decreased in abundance after forest fire in Pumanque. This in part might be related to potentially lower rates of nitrification in burned soils from Pumanque.
Our results indicate that forest fires have differently altered prokaryote community-level physiological profiles and community structure in the study sites located in the Mediterranean climate zone of central Chile, providing local evidence that support and complement previous findings in other Mediterranean ecosystems across the globe. Legacy of land burning was still clearly shaping soil prokaryotic community structure rather than community-level physiological profiles, particularly in the most recent fire-affected site. However, we acknowledge that these findings cannot be completely disentangled from site effects (e.g., soil type, vegetation, and climatic conditions); thus, resilience of soil prokaryotic communities can also be explained by these factors. Nonetheless, changes in the relative abundance of coinciding taxonomic groups were observed in fire-affected soils at both locations, despite environmental difference between sites. Moreover, variability between community structures was moderated by changes in soil physicochemical properties known to be affected by fires and possibly by post-fire substrate availability. In an ecosystem that is more frequently affected by fires and whose vegetation is not well adapted to such events, as is the case for Mediterranean ecosystems of central Chile, our results allow for a better understanding of the state of soil prokaryotes and their physiological conditions after different time periods following fire, and suggest that biotic belowground responses to fire occurrence at this particular Mediterranean ecosystem are comparable with other Mediterranean biomes across the globe. These findings highlight the importance of including such soil microbiological assessments in combination to aboveground conditions to better understand forest ecosystem resilience and restoration processes under local conditions where natural recovery is hampered.
Data Availability Statement
The dataset generated for this study can be found in the raw sequencing data deposited in the Sequence Read Archive (SRA) of NCBI under accession number PRJNA784510.
Author Contributions
HA, CY, and CR designed the research and performed the laboratory work. HA, CY, TG-C, and MH performed the data analysis. HA and CR wrote the manuscript and designed the tables and figures. CY, TG-C, and MH made revisions of the manuscript. All authors revised the article and approved the final version.
Funding
This study was funded by the National Agency for Research and Development (ANID), through the Fund for Scientific and Technological Development (FONDECYT) Project No. 11180869. TG-C acknowledges the support of the Marie Skłodowska-Curie Actions ROOTPHENOBIOME 698 839235 (H2020-EU.1.3.2.).
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.
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.2022.824813/full#supplementary-material
References
Acea, M. J., and Carballas, T. (1996). Changes in physiological groups of microorganisms in soil following wildfire. FEMS Microbiol. Ecol. 20, 33–39. doi: 10.1111/j.1574-6941.1996.tb00302.x
Alcañiz, M., Outeiro, L., Francos, M., and Úbeda, X. (2018). Effects of prescribed fires on soil properties: a review. Sci. Total Environ. 613–614, 944–957. doi: 10.1016/j.scitotenv.2017.09.144
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Anderson, M. J., and Willis, T. J. (2003). Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology 84, 511–525. doi: 10.1890/0012-9658(2003)084[0511:caopca]2.0.co;2
Apprill, A., McNally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 75, 129–137. doi: 10.3354/ame01753
Ball, P. N., MacKenzie, M. D., DeLuca, T. H., and Montana, W. E. H. (2010). Wildfire and charcoal enhance nitrification and ammonium-oxidizing bacterial abundance in dry montane forest soils. J. Environ. Qual. 39, 1243–1253. doi: 10.2134/jeq2009.0082
Bárcenas-Moreno, G., García-Orenes, F., Mataix-Solera, J., Mataix-Beneyto, J., and Bååth, E. (2011). Soil microbial recolonisation after a fire in a Mediterranean forest. Biol. Fertil. Soils 47, 261–272. doi: 10.1007/s00374-010-0532-2
Barreiro, A., Fontúrbel, M. T., Lombao, A., Martín, A., Vega, J. A., Fernández, C., et al. (2015). Using phospholipid fatty acid and community level physiological profiling techniques to characterize soil microbial communities following an experimental fire and different stabilization treatments. CATENA 135, 419–429. doi: 10.1016/j.catena.2014.07.011
Baver, L. D., Gardner, W. H., and Gardner, W. R. (1973). Física de suelo. México: Centro Regional de Ayuda Técnica.
Bengtsson-Palme, J., Hartmann, M., Eriksson, K. M., Pal, C., Thorell, K., Larsson, D. G. J., et al. (2015). Metaxa 2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data. Mol. Ecol. Resour. 15, 1403–1414. doi: 10.1111/1755-0998.12399
Bouyoucos, G. J. (1962). Hidrometer method improved for marking partied size analyses of soil. Agron. J. 54, 464–465. doi: 10.2134/agronj1962.00021962005400050028x
Brochier-Armanet, C., Gribaldo, S., and Forterre, P. (2012). Spotlight on the Thaumarchaeota. ISME J. 6, 227–230. doi: 10.1038/ismej.2011.145
Brown, J. K., and Smith, J. K. (2000). Wildland Fire in Ecosystems Effects of Fire on Flora. Ogden, UT: U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station.
Bustamante, M., Verdejo, V., Zúñiga, C., Espinosa, F., Orlando, J., and Carú, M. (2012). Comparison of water availability effect on ammonia-oxidizing bacteria and archaea in microcosms of a Chilean semiarid soil. Front. Microbiol. 3:282. doi: 10.3389/fmicb.2012.00282
Center of Climate and Resilience Research (n.d.). Datos de Precipitación [WWW Document]. Available online at: https://www.cr2.cl/datos-de-precipitacion/?cp_Precipitacion=2.Web (accessed August 30, 2021).
Certini, G. (2005). Effects of fire on properties of forest soils: a review. Oecologia 143, 1–10. doi: 10.1007/s00442-004-1788-8
CIREN (1996). Estudio Agrológico. Descripciones de Suelos Materiales y Símbolos. Publicación 114. VI Región. Santiago: CIREN.
Cobo-Díaz, J. F., Fernández-González, A. J., Villadas, P. J., Robles, A. B., Toro, N., and Fernández-López, M. (2015). Metagenomic assessment of the potential microbial nitrogen pathways in the rhizosphere of a mediterranean forest after a wildfire. Microb. Ecol. 69, 895–904. doi: 10.1007/s00248-015-0586-7
Covington, W. W., and Sackett, S. S. (1992). Soil mineral nitrogen changes following prescribed burning in ponderosa pine. For. Ecol. Manag. 54, 175–191. doi: 10.1016/0378-1127(92)90011-W
de Armas-Ricard, M., Orlando, J., Bustamante, R., and Carú, M. (2016). Microbial communities of bulk and Eschscholzia californica rhizosphere soils at two altitudes in Central Chile. J. Soil Sci. Plant Nutr. 16, 1–13. doi: 10.4067/S0718-95162016005000001
Díaz-Raviña, M., Prieto, A., and Bååth, E. (1996). Bacterial activity in a forest soil after soil heating and organic amendments measured by the thymidine and leucine incorporation techniques. Soil Biol. Biochem. 28, 419–426. doi: 10.1016/0038-0717(95)00156-5
Dworkin, M. (2006). “Prokaryotic life cycles,” in The Prokaryotes, eds M. Dworkin, S. Falkow, E. Rosenberg, K. H. Schleifer, and E. Stackebrandt (New York, NY: Springer New York), 140–166. doi: 10.1007/0-387-30742-7_6
Edgar, R. C. (2016). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. BioRxiv [Preprint] doi: 10.1101/081257
Edgar, R. C., and Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 31, 3476–3482. doi: 10.1093/bioinformatics/btv401
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Eichorst, S. A., Kuske, C. R., and Schmidt, T. M. (2011). Influence of plant polymers on the distribution and cultivation of bacteria in the phylum acidobacteria. Appl. Environ. Microbiol. 77, 586–596. doi: 10.1128/AEM.01080-10
Eilers, K. G., Lauber, C. L., Knight, R., and Fierer, N. (2010). Shifts in bacterial community structure associated with inputs of low molecular weight carbon compounds to soil. Soil Biol. Biochem. 42, 896–903. doi: 10.1016/j.soilbio.2010.02.003
Esteves, T. C. J., Kirkby, M. J., Shakesby, R. A., Ferreira, A. J. D., Soares, J. A. A., Irvine, B. J., et al. (2012). Mitigating land degradation caused by wildfire: application of the PESERA model to fire-affected sites in central Portugal. Geoderma 191, 40–50. doi: 10.1016/j.geoderma.2012.01.001
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fletcher, W. J., and Zielhofer, C. (2013). Fragility of western mediterranean landscapes during holocene rapid climate changes. CATENA 103, 16–29. doi: 10.1016/j.catena.2011.05.001
Fuentes-Castillo, T., Scherson, R. A., Marquet, P. A., Fajardo, J., Corcoran, D., Román, M. J., et al. (2019). Modelling the current and future biodiversity distribution in the Chilean Mediterranean hotspot. The role of protected areas network in a warmer future. Divers. Distrib. 25, 1897–1909. doi: 10.1111/ddi.12988
Fultz, L. M., Moore-Kucera, J., Dathe, J., Davinic, M., Perry, G., Wester, D., et al. (2016). Forest wildfire and grassland prescribed fire effects on soil biogeochemical processes and microbial communities: two case studies in the semi-arid Southwest. Appl. Soil Ecol. 99, 118–128. doi: 10.1016/j.apsoil.2015.10.023
García-Carmona, M., Marín, C., García-Orenes, F., and Rojas, C. (2021). Contrasting organic amendments induce different short-term responses in soil abiotic and biotic properties in a fire-affected native Mediterranean forest in Chile. J. Soil Sci. Plant Nutr. 21, 2105–2114. doi: 10.1007/s42729-021-00506-z
García-Vega, D., and Newbold, T. (2020). Assessing the effects of land use on biodiversity in the world’s drylands and Mediterranean environments. Biodivers. Conserv. 29, 393–408. doi: 10.1007/s10531-019-01888-4
Garland, J. L. (1996). Analytical approaches to the characterization of samples of microbial communities using patterns of potential C source utilization. Soil Biol. Biochem. 28, 213–221. doi: 10.1016/0038-0717(95)00112-3
Garland, J. L., and Mills, A. L. (1991). Classification and characterization of heterotrophic microbial communities on the basis of patterns of community-level sole-carbon-source utilization. Appl. Environ. Microbiol. 57, 2351–2359. doi: 10.1128/aem.57.8.2351-2359.1991
Goberna, M., García, C., Insam, H., Hernández, M. T., and Verdú, M. (2012). Burning fire-prone mediterranean shrublands: immediate changes in soil microbial community structure and ecosystem functions. Microb. Ecol. 64, 242–255. doi: 10.1007/s00248-011-9995-4
Grayston, S. J., Vaughan, D., and Jones, D. (1997). Rhizosphere carbon flow in trees, in comparison with annual plants: the importance of root exudation and its impact on microbial activity and nutrient availability. Appl. Soil Ecol. 5, 29–56. doi: 10.1016/S0929-1393(96)00126-6
Hart, S. C., DeLuca, T. H., Newman, G. S., MacKenzie, M. D., and Boyle, S. I. (2005). Post-fire vegetative dynamics as drivers of microbial community structure and function in forest soils. For. Ecol. Manag. 220, 166–184. doi: 10.1016/j.foreco.2005.08.012
Hernesmaa, A., Bjorklof, K., Kiikkila, O., Fritze, H., Haahtela, K., and Romantschuk, M. (2005). Structure and function of microbial communities in the rhizosphere of Scots pine after tree-felling. Soil Biol. Biochem. 37, 777–785. doi: 10.1016/j.soilbio.2004.10.010
Janssen, P. H. (2006). Identifying the Dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72, 1719–1728. doi: 10.1128/AEM.72.3.1719-1728.2006
Jiang, W. (2018). A protocol for quantizing total bacterial 16S rDNA in plasma as a marker of microbial translocation in vivo. Cell. Mol. Immunol. 15, 937–939. doi: 10.1038/cmi.2018.3
Jones, R. T., Robeson, M. S., Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J. 3, 442–453. doi: 10.1038/ismej.2008.127
Kandeler, E., Poll, C., Frankenberger, W. T., and Tabatabai, M. A. (2011). “Nitrogen cycle enzymes,” in Methods of Soil Enzymology, ed. R. P. Dick (Madison, WI: Soil Science Society of America), 211–245. doi: 10.2136/sssabookser9.c10
Keeley, J. E. (2009). Fire intensity, fire severity and burn severity: a brief review and suggested usage. Int. J. Wildland Fire 18:116. doi: 10.1071/WF07049
Kemper, W. D., and Rosenau, R. C. (1986). Aggregate stability and size distribution. Methods Soil Anal. 9, 425–442. doi: 10.2136/sssabookser5.1.2ed.c17
Khianngam, S., Tanasupawat, S., Akaracharanya, A., Kim, K. K., Lee, K. C., and Lee, J.-S. (2010). Cohnella thailandensis sp. nov., a xylanolytic bacterium from Thai soil. Int. J. Syst. Evol. Microbiol. 60, 2284–2287. doi: 10.1099/ijs.0.015859-0
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A., and Kuramae, E. E. (2016). The ecology of Acidobacteria: moving beyond genes and genomes. Front. Microbiol. 7:744. doi: 10.3389/fmicb.2016.00744
Kim, J.-H., and Kim, W. (2016). Tumebacillus soli sp. nov., isolated from non-rhizosphere soil. Int. J. Syst. Evol. Microbiol. 66, 2192–2197. doi: 10.1099/ijsem.0.001009
Kindt, R., and Coe, R. (2005). Tree Diversity Analysis; A Manual and Software for Common Statistical Methods for Ecological and Biodiversity Studies. Nairobi: World Agroforestry Centre (ICRAF).
Kishimoto, N., Kosako, Y., and Tano, T. (1991). Acidobacterium capsulatum gen. nov., sp. nov.: an acidophilic chemoorganotrophic bacterium containing menaquinone from acidic mineral environment. Curr. Microbiol. 22, 1–7. doi: 10.1007/BF02106205
Knelman, J., Graham, E., Ferrenberg, S., Lecoeuvre, A., Labrado, A., Darcy, J., et al. (2017). Rapid shifts in soil nutrients and decomposition enzyme activity in early succession following forest fire. Forests 8:347. doi: 10.3390/f8090347
Lagos, M. L., Maruyama, F., Nannipieri, P., Mora, M. L., Ogram, A., and Jorquera, M. A. (2015). Current overview on the study of bacteria in the rhizosphere by modern molecular techniques: a mini-review. J. Soil Sci. Plant Nutr. 15, 504–523. doi: 10.4067/S0718-95162015005000042
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. doi: 10.1128/AEM.00335-09
Lauber, C. L., Strickland, M. S., Bradford, M. A., and Fierer, N. (2008). The influence of soil properties on the structure of bacterial and fungal communities across land-use types. Soil Biol. Biochem. 40, 2407–2415. doi: 10.1016/j.soilbio.2008.05.021
Legendre, P., and Andersson, M. J. (1999). Distance-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 69, 1–24. doi: 10.1890/0012-9615(1999)069[0001:dbratm]2.0.co;2
Li, J., Li, C., Kou, Y., Yao, M., He, Z., and Li, X. (2020). Distinct mechanisms shape soil bacterial and fungal co-occurrence networks in a mountain ecosystem. FEMS Microbiol. Ecol. 96:fiaa030. doi: 10.1093/femsec/fiaa030
Li, W., Niu, S., Liu, X., and Wang, J. (2019). Short-term response of the soil bacterial community to differing wildfire severity in Pinus tabulaeformis stands. Sci. Rep. 9:1148. doi: 10.1038/s41598-019-38541-7
Liu, X., Lindemann, W. C., Whitford, W. G., and Steiner, R. L. (2000). Microbial diversity anti activity of disturbed soil in the northern chihuahuan desert. Biol. Fertil. Soils 32, 243–249. doi: 10.1007/s003740000242
Mandakovic, D., Rojas, C., Maldonado, J., Latorre, M., Travisany, D., Delage, E., et al. (2018). Structure and co-occurrence patterns in microbial communities under acute environmental stress reveal ecological factors fostering resilience. Sci. Rep. 8:5875. doi: 10.1038/s41598-018-23931-0
Marín, C., and Rojas, C. (2020). Early responses of soil health indicators to organic amendments and plant establishment at a fire-affected sclerophyllous forest. Agro Sur 48, 17–27. doi: 10.4206/agrosur.2020.v48n2-03
Martín, A., Díaz-Raviña, M., and Carballas, T. (2009). Evolution of composition and content of soil carbohydrates following forest wildfires. Biol. Fertil. Soils 45, 511–520. doi: 10.1007/s00374-009-0363-1
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17:10. doi: 10.14806/ej.17.1.200
Mataix-Solera, J., Guerrero, C., García-Orenes, F., Bárcenas, G., and Torres, M. (2009). “Forest fire effects on soil microbiology,” in Fire Effects on Soils and Restoration Strategies, eds A. Cerdà and P. Robichaud (Enfield: Science Publishers), 133–175. doi: 10.1201/9781439843338-c5
Mosugu, M. E., Chude, V. O., Esu, I. E., Kparmwang, T., and Malgwi, W. B. (1999). Contents and profile distribution of three forms of free iron oxides in three ultisois and an alfisol in Nigeria. Commun. Soil Sci. Plant Anal. 30, 1013–1024. doi: 10.1080/00103629909370264
Moya, D., Fonturbel, M. T., Lucas-Borja, M. E., Peña, E., Alfaro-Sanchez, R., Plaza-Álvarez, P. A., et al. (2021). Burning season and vegetation coverage influenced the community-level physiological profile of Mediterranean mixed-mesogean pine forest soils. J. Environ. Manag. 277:111405. doi: 10.1016/j.jenvman.2020.111405
Mroz, G. D., Jurgensen, M. F., Harvey, A. E., and Larsen, M. J. (1980). Effects of fire on nitrogen in forest floor horizons. Soil Sci. Soc. Am. J. 44, 395–400. doi: 10.2136/sssaj1980.03615995004400020038x
Neary, D. G., Klopatek, C. C., DeBano, L. F., and Ffolliott, P. F. (1999). Fire effects on belowground sustainability: a review and synthesis. For. Ecol. Manag. 122, 51–71. doi: 10.1016/S0378-1127(99)00032-8
Norouzi, M., and Ramezanpour, H. (2013). Effect of fire on chemical forms of iron and manganese in forest soils of Iran. Environ. Forensics 14, 169–177. doi: 10.1080/15275922.2013.781077
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020). vegan: Community Ecology Package.
Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023
Pérez-Valera, E., Goberna, M., and Verdú, M. (2019). Fire modulates ecosystem functioning through the phylogenetic structure of soil bacterial communities. Soil Biol. Biochem. 129, 80–89. doi: 10.1016/j.soilbio.2018.11.007
Pérez-Valera, E., Verdú, M., Navarro-Cano, J. A., and Goberna, M. (2020). Soil microbiome drives the recovery of ecosystem functions after fire. Soil Biol. Biochem. 149:107948. doi: 10.1016/j.soilbio.2020.107948
Prashar, P., Kapoor, N., and Sachdeva, S. (2014). Rhizosphere: its structure, bacterial diversity and significance. Rev. Environ. Sci. Biotechnol. 13, 63–77. doi: 10.1007/s11157-013-9317-z
Pressler, Y., Moore, J. C., and Cotrufo, M. F. (2019). Belowground community responses to fire: meta-analysis reveals contrasting responses of soil microorganisms and mesofauna. Oikos 128, 309–327. doi: 10.1111/oik.05738
Puri, A., Padda, K. P., and Chanway, C. P. (2016). Evidence of nitrogen fixation and growth promotion in canola (Brassica napus L.) by an endophytic diazotroph Paenibacillus polymyxa P2b-2R. Biol. Fertil. Soils 52, 119–125. doi: 10.1007/s00374-015-1051-y
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). 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
R Core Team (2021). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Raymond, J., Siefert, J. L., Staples, C. R., and Blankenship, R. E. (2004). The natural history of nitrogen fixation. Mol. Biol. Evol. 21, 541–554. doi: 10.1093/molbev/msh047
Rodríguez, J., González-Pérez, J. A., Turmero, A., Hernández, M., Ball, A. S., González-Vila, F. J., et al. (2018). Physico-chemical and microbial perturbations of Andalusian pine forest soils following a wildfire. Sci. Total Environ. 634, 650–660. doi: 10.1016/j.scitotenv.2018.04.028
Rodríguez-Loinaz, G., Onaindia, M., Amezaga, I., Mijangos, I., and Garbisu, C. (2008). Relationship between vegetation diversity and soil functional diversity in native mixed-oak forests. Soil Biol. Biochem. 40, 49–60. doi: 10.1016/j.soilbio.2007.04.015
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Rousk, J., Bååth, E., Brookes, P. C., Lauber, C. L., Lozupone, C., Caporaso, J. G., et al. (2010). Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351. doi: 10.1038/ismej.2010.58
Rundel, P. W., Arroyo, M. T. K., Cowling, R. M., Keeley, J. E., Lamont, B. B., Pausas, J. G., et al. (2018). Fire and plant diversification in mediterranean-climate regions. Front. Plant Sci. 9:851. doi: 10.3389/fpls.2018.00851
Rundel, P. W., Arroyo, M. T. K., Cowling, R. M., Keeley, J. E., Lamont, B. B., and Vargas, P. (2016). Mediterranean biomes: evolution of their vegetation, floras, and climate. Annu. Rev. Ecol. Evol. Syst. 47, 383–407. doi: 10.1146/annurev-ecolsys-121415-032330
Sadeghifar, M., Agha, A. B. A., and Pourreza, M. (2020). Comparing soil microbial eco-physiological and enzymatic response to fire in the semi-arid Zagros woodlands. Appl. Soil Ecol. 147:103366. doi: 10.1016/j.apsoil.2019.103366
Sadzawka, A., Carrasco, R., Mora, M., Flores, H., and Neaman, A. (2006).. Métodos de Análisis Recomendados Para los Suelos de Chile. Revisión 2006. Serie Actas INIA N° 34. Santiago: Instituto de Investigaciones Agropecuarias (INIA).
Sait, M., Davis, K. E. R., and Janssen, P. H. (2006). Effect of pH on isolation and distribution of members of subdivision 1 of the phylum Acidobacteria occurring in soil. Appl. Environ. Microbiol. 72, 1852–1857. doi: 10.1128/AEM.72.3.1852-1857.2006
Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Singh, A. K., Kushwaha, M., Rai, A., and Singh, N. (2017). Changes in soil microbial response across year following a wildfire in tropical dry forest. For. Ecol. Manag. 391, 458–468. doi: 10.1016/j.foreco.2017.02.042
Sneha, G. R., Swarnalakshmi, K., Sharma, M., Reddy, K., Bhoumik, A., Suman, A., et al. (2021). Soil type influence nutrient availability, microbial metabolic diversity, eubacterial and diazotroph abundance in chickpea rhizosphere. World J. Microbiol. Biotechnol. 37:167. doi: 10.1007/s11274-021-03132-0
Sofo, A., and Ricciuti, P. (2019). A standardized method for estimating the functional diversity of soil bacterial community by Biolog® EcoPlatesTM assay—the case study of a sustainable olive orchard. Appl. Sci. 9:4035. doi: 10.3390/app9194035
Sun, H., Santalahti, M., Pumpanen, J., Köster, K., Berninger, F., Raffaello, T., et al. (2016). Bacterial community structure and function shift across a northern boreal forest fire chronosequence. Sci. Rep. 6:32411. doi: 10.1038/srep32411
Sun, M., Chen, B., Wang, H., Wang, N., Ma, T., Cui, Y., et al. (2021). Microbial interactions and roles in soil fertility in seasonal freeze-thaw periods under different straw returning strategies. Agriculture 11:779. doi: 10.3390/agriculture11080779
Switzer, J. M., Hope, G. D., Grayston, S. J., and Prescott, C. E. (2012). Changes in soil chemical and biological properties after thinning and prescribed fire for ecosystem restoration in a Rocky Mountain Douglas-fir forest. For. Ecol. Manag. 275, 1–13. doi: 10.1016/j.foreco.2012.02.025
Tourna, M., Stieglmeier, M., Spang, A., Konneke, M., Schintlmeister, A., Urich, T., et al. (2011). Nitrososphaera viennensis, an ammonia oxidizing archaeon from soil. Proc. Natl. Acad. Sci. 108, 8420–8425. doi: 10.1073/pnas.1013488108
U.S. EPA (2007). Method 3051A: Microwave Assisted Acid Digestion of Sediments, Sludges, Soils, and Oils. Revision 1. Washington, DC: U.S. EPA.
Vijayan, A., Vattiringal Jayadradhan, R. K., Pillai, D., Prasannan Geetha, P., Joseph, V., and Isaac Sarojini, B. S. (2021). Nitrospira as versatile nitrifiers: taxonomy, ecophysiology, genome characteristics, growth, and metabolic diversity. J. Basic Microbiol. 61, 88–109. doi: 10.1002/jobm.202000485
Wang, C., Wang, G., Wang, Y., Rafique, R., Ma, L., Hu, L., et al. (2016). Fire alters vegetation and soil microbial community in alpine meadow. Land Degrad. Dev. 27, 1379–1390. doi: 10.1002/ldr.2367
Wei, T., and Simko, V. (2021). R Package “corrplot”: Visualization of a Correlation Matrix (Version 0.90).
Xiang, X., Shi, Y., Yang, J., Kong, J., Lin, X., Zhang, H., et al. (2014). Rapid recovery of soil bacterial communities after wildfire in a Chinese boreal forest. Sci. Rep. 4:3829. doi: 10.1038/srep03829
Yang, S., Zheng, Q., Yang, Y., Yuan, M., Ma, X., Chiariello, N. R., et al. (2020). Fire affects the taxonomic and functional composition of soil microbial communities, with cascading effects on grassland ecosystem functioning. Glob. Chang. Biol. 26, 431–442. doi: 10.1111/gcb.14852
Keywords: rhizosphere, bacteria, Biolog EcoPlates, ecosystem recovery, wildfires
Citation: Aponte H, Galindo-Castañeda T, Yáñez C, Hartmann M and Rojas C (2022) Microbial Community-Level Physiological Profiles and Genetic Prokaryotic Structure of Burned Soils Under Mediterranean Sclerophyll Forests in Central Chile. Front. Microbiol. 13:824813. doi: 10.3389/fmicb.2022.824813
Received: 29 November 2021; Accepted: 24 February 2022;
Published: 28 April 2022.
Edited by:
Yoanna Eissler, Universidad de Valparaíso, ChileReviewed by:
Julieta Orlando, University of Chile, ChileGergely Krett, Eötvös Loránd University, Hungary
Tibor Szili-Kovács, Hungarian Academy of Sciences, Hungary
Copyright © 2022 Aponte, Galindo-Castañeda, Yáñez, Hartmann and Rojas. 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: Claudia Rojas, claudia.rojas@uoh.cl