Peltigera frigida Lichens and Their Substrates Reduce the Influence of Forest Cover Change on Phosphate Solubilizing Bacteria

Phosphorus (P) is one of the most critical macronutrients in forest ecosystems. More than 70 years ago, some Chilean Patagonian temperate forests suffered wildfires and the subsequent afforestation with foreign tree species such as pines. Since soil P turnover is interlinked with the tree cover, this could influence soil P content and bioavailability. Next to soil microorganisms, which are key players in P transformation processes, a vital component of Patagonian temperate forest are lichens, which represent microbial hotspots for bacterial diversity. In the present study, we explored the impact of forest cover on the abundance of phosphate solubilizing bacteria (PSB) from three microenvironments of the forest floor: Peltigera frigida lichen thallus, their underlying substrates, and the forest soil without lichen cover. We expected that the abundance of PSB in the forest soil would be strongly affected by the tree cover composition since the aboveground vegetation influences the edaphic properties; but, as P. frigida has a specific bacterial community, lichens would mitigate this impact. Our study includes five sites representing a gradient in tree cover types, from a mature forest dominated by the native species Nothofagus pumilio, to native second-growth forests with a gradual increase in the presence of Pinus contorta in the last sites. In each site, we measured edaphic parameters, P fractions, and the bacterial potential to solubilize phosphate by quantifying five specific marker genes by qPCR. The results show higher soluble P, labile mineral P, and organic matter in the soils of the sites with a higher abundance of P. contorta, while most of the molecular markers were less abundant in the soils of these sites. Contrarily, the abundance of the molecular markers in lichens and substrates was less affected by the tree cover type. Therefore, the bacterial potential to solubilize phosphate is more affected by the edaphic factors and tree cover type in soils than in substrates and thalli of P. frigida lichens. Altogether, these results indicate that the microenvironments of lichens and their substrates could act as an environmental buffer reducing the influence of forest cover composition on bacteria involved in P turnover.

Introduction

Lichens are important components of the cryptogamic biota in forest ecosystems and represent interesting ecological niches for microorganisms. Beyond the classical definition of lichens, which includes a mutualistic symbiotic association between a fungus (mycobiont) and one or more photosynthetic microorganisms (photobiont/s; Nash, 2008), these organisms also harbor an associated microbiome (Cernava et al., 2015; Leiva et al., 2016, 2021; Aschenbrenner et al., 2017; Hawksworth and Grube, 2020). In recent years, lichen-associated bacteria have been described as forming a highly structured biofilm covering the lichen thallus (Grube and Berg, 2009; Bates et al., 2011; Hodkinson et al., 2012), which can play an essential role in lichens by facilitating, for example, the supply of crucial nutrients and therefore making their survival and growth possible in nutrient-poor substrates (Grimm et al., 2021). Besides, these bacterial communities associated with lichens could play an essential role in forest nutrition because of their potential to express enzymes driving carbon turnover, including cellulase, xylanase, and amylase, as well as their ability to fix nitrogen and solubilize phosphate (González et al., 2005; Liba et al., 2006; Selbmann et al., 2010; Cardinale et al., 2012; Sigurbjörnsdóttir et al., 2016; Almendras et al., 2018a,b). These bacterial communities associated with the lichen thallus are more than only an extension of those found in the substrate where the lichens grow, but an additional component of the traditional symbiosis, mainly due to the dissimilarity in abundance and the presence of specific bacterial groups throughout different lichens (Bates et al., 2011; Cardinale et al., 2012; Printzen et al., 2012; Leiva et al., 2016, 2021).

Phosphorus (P) is one of the most critical macronutrients in all ecosystems, and consequently, P availability in soils also affects the growth and performance of lichens. Microorganisms, including bacteria, play an important role by making P bioavailable (Richardson et al., 2001; Sashidhar and Podile, 2010; Chhabra et al., 2013; Sharma et al., 2013; Zhao et al., 2014). The way used by microorganisms to increase P bioavailability depends on phosphate bindings with organic and mineral molecules. The solubilization of organic phosphate involves several key enzymes, including alkaline and acid phosphatases, phytases, and phosphonatases (Rodríguez et al., 2006; Sharma et al., 2013). Mineral phosphate solubilization depends on the production of organic acids, which is essential for the dissolution of many poorly soluble mineral phosphates. The most studied organic acid is gluconic acid, produced by the membrane-bound glucose dehydrogenase, which catalyzes the direct oxidation of glucose (Goldstein, 1995). In forests, litterfall is the primary route of nutrient returns to the soil, releasing nutrients by decomposition (Osman, 2013). However, different plant species differ in the amount and quality of organic matter they add to the soil (Hooper et al., 2000; Finér et al., 2007; Šnajdr et al., 2013). Therefore, tree species may differently influence soil chemical properties, such as the pH or the relative content and chemical form of macronutrients (Binkley and Valentine, 1991; Augusto et al., 2000; Frouz et al., 2009; Iovieno et al., 2010). For example, P fluxes in coniferous and deciduous forests differ, with the concentration of P in the forest floor under coniferous trees being about three times higher than in deciduous forests, probably because the P content of coniferous foliage is larger than in deciduous one (Sohrt et al., 2017).

Peltigera frigida is one of the most conspicuous lichens in forests of the Coyhaique National Reserve in Aysén Region (Leiva et al., 2016), but its abundance is low in other parts of Southern Chile (Ramírez-Fernández et al., 2013; Zúñiga et al., 2015). It has been proposed that its adaptation could be partly due to the peculiarities of its microbiome (Leiva et al., 2021), as it could be related to soil properties due to vegetation composition and historic events suffered by the study site (Fajardo and Gundale, 2015). In particular, the main events affecting the study site were wildfires and subsequent reforestation with pines, which are reported as crucial drivers of phosphorus availability (Dickie et al., 2014; Butler et al., 2018).

Based on this information, we asked if the microenvironments associated with the terricolous cyanolichen P. frigida (thallus and substrate) reduce the influence of the tree cover type on the diversity of phosphate solubilizing bacteria (PSB). This question is driven by previous studies showing that P. frigida thalli and substrates constitute different microenvironments that select specific bacterial communities (Leiva et al., 2016, 2021), which could reduce the impact that the tree cover type exerts on the diversity of PBS in soils without lichen influence. In order to answer this question, we investigated the diversity of PSB along a gradient including sampling sites from mature native forest to second-growth forests afforested with exotic pines in the Coyhaique National Reserve (Aysén Region, Chile), considering three microenvironments related to the forest floor: (i) P. frigida thalli, (ii) their associated substrate (i.e., the soil beneath each thallus), and (iii) forest soil without lichen influence.

Materials and Methods

Selection and Characterization of Sampling Sites

The Coyhaique National Reserve (45°31′42.96”S, 72°1′51.95”W) is located in the Aysén Region of Southern Chile. In this reserve, deciduous forests with a predominance of Nothofagus pumilio constitute the typical vegetation. The native mature forests were affected by several wildfires in the region between 1930 and 1950 (Quintanilla Pérez, 2007; Fajardo and Gundale, 2015), after which natural regeneration resulted in native second-growth forests. Also, to protect the soil from erosion, plantations of various kinds of pines (mainly Pinus contorta) were part of afforestation programs organized by local authorities. However, the study area is not currently exposed to regular management practices since it is a protected area, and then the current forest community is the result of long-standing colonization.

Five ~200 m transects away ~300 m from each other (S1 to S5) were selected along a gradient in the forest cover composition present in the reserve trails (Supplementary Figure S1).

Forest cover at each site was estimated using the middle of the intercept point (Goodall, 1952) and the point-centered quarter (Mitchell, 2010) methods. The first one was used to calculate the relative coverage (RC) of each species in 2 m height, according to the formula $R{C}_i=\frac{n_i}{N_T}$, where n_i is the frequency of the species i and N_T is the number of sampling points per site (Supplementary Table S1). The second method was used to calculate various indices from the field-collected data for each site (Supplementary Tables S2–S6). In this method, a reference point was selected every 20 m and divided into four quadrants. The closest tree was described at each of these quadrants by three parameters: species, distance to the reference point (d_i), and circumference at breast height (CBH). Then, three basic values were calculated: the radius of the tree trunk (r), according to the formula $r_i=\frac{CB{H}_i}{2·\pi }$; the basal area covered by each tree with the formula $A_i=\pi ·r_i{}^2$; and the total density (TD) of each species following $T{D}_i=\frac{1}{{\left(\raisebox{1ex}{$d_i$}\!\left/ \!\raisebox{-1ex}{$Q$}\right.\right)}^2}$, where Q is the total number of sampled trees. After this, the relative density (RD) of each species was calculated with the formula $R{D}_i=\frac{n_i}{Q}$, where n_i corresponds to the total number of trees of each species. The absolute density (D) was calculated dividing the relative density by the total density according to $D_i=\frac{R{D}_i}{T{D}_i}$. Additionally, five other indices were calculated by the point-centered quarter method (Supplementary Table S7). The frequency (F) was calculated according to the formula $F_i=\frac{j_i}{k}$, where j_i is the total number of points of each species, and k is the number of sampled points. Then, the relative frequency (RF) was calculated by dividing the F value of each species by the sum for all three species, according to $R{F}_i=\frac{F_i}{\sum F}$. The coverage (C) of each species was calculated with the formula $C_i=\frac{A_i·D_i}{n_i}$, and the relative coverage (RC) by $R{C}_i=\frac{C_i}{\sum C}$. Finally, the importance value (IV) was calculated according to the formula $I{V}_i=R{F}_i+R{C}_i+R{D}_i$, where the relative density (RD) of each species was calculated with the formula $R{D}_i=\frac{n_i}{Q}$, where n_i corresponds to the total number of trees of each species and Q is the total number of sampled trees. The importance value was relativized as Relative Importance (proportion of the total importance value) as it is easier to compare.

Samples Collection

From the transect defined at each study site, five Peltigera frigida thalli and their associated substrates (Supplementary Table S8), alongside five soil samples from the surroundings (Supplementary Table S9), were collected. The lichen-associated samples were collected at least 1 m from the next closest to minimize the resampling of the same genetic individual. Besides, the soil samples were collected at least 1 m away from any lichen. The samples were placed in paper bags and transported in cooled containers. In the laboratory, the thalli were separated from the attached substrate with sterile brush and spatula. Thallus samples were stored in paper bags at room temperature, while the substrate and soil samples were sieved and stored in plastic tubes at −20°C.

Physicochemical Parameters

From each soil sample, physicochemical parameters were determined as follows: pH was measured by potentiometry using a 1:2 ratio in 1 M KCl. Moisture content (MC; Steubing et al., 2002) and total organic matter (TOM; Sadzawka et al., 2006) were calculated gravimetrically before and after desiccation and calcination, respectively. Nitrate (N-NO₃⁻), ammonium (N-NH₄⁺), and phosphorus (P Bray) concentrations were measured spectrophotometrically from 1:10 soil extracts in deionized water, 1 M KCl, and Bray 1 extract, respectively (Bray and Kurtz, 1945; Steubing et al., 2002). Phosphorus sequential fractionation was determined by following Hedley’s protocol (Hedley et al., 1982) with the modification of do Nascimento et al. (2015). The soil was sequentially extracted with distilled water, 0.5 M NaHCO₃ (pH 8), 0.1 M NaOH, and mineral phosphorus (Pm) was calorimetrically quantified (Murphy and Riley, 1962). Subsequently, all samples were digested with H₂SO₄ and K₂S₂O₈ to determine organic phosphorus (Po) as the difference between total and mineral phosphorus (Po = Pt-Pm).

DNA Extraction

DNA was extracted from thallus, substrate, and soil samples using the DNeasy PowerSoil Kit (Qiagen), with 0.25 g per soil and substrate samples, as suggested by the manufacturer, and 0.15 g per thallus sample according to pre-experiments performed (data not shown). Modifications compared to the standard protocol included an overnight incubation at 4°C after adding the precipitation solution C2 (step 8) and an increase in time from 30 s to 1 min for the centrifugation steps after bead beating (step 6), ethanol washing (step 16), and elution (step 21).

Identification of the Lichen Symbionts

From the DNA of lichen thalli, mycobionts were identified by analyzing the fungal 28S rDNA region amplified with primers LIC24R (Miadlikowska and Lutzoni, 2000) and LR7 (Vilgalys and Hester, 1990), and the ITS region amplified with primers ITS1F (Gardes and Bruns, 1993) and ITS4 (White et al., 1990). In addition, photobionts were identified using the cyanobacterial 16S rDNA region with the primers PCR1 and PCR18 (Wilmotte et al., 1993). All amplicons were sequenced with the forward primers using the Genetic Analyzer 3730XL (Applied Biosystems) using a sequencing service (Macrogen, Seoul, Korea).

Mycobiont identity was confirmed by a de novo phylogenetic analysis of the 28S and ITS sequences, carried out on the T-BAS web platform using the Peltigera reference tree (Carbone et al., 2019). In the case of the cyanobionts, the evolutionary history was inferred using the Neighbor-Joining method. For this, evolutionary distances were calculated using Kimura’s 2-parameter method, as suggested by MEGA7 (Kumar et al., 2016). The analysis involved 55 nucleotide sequences, including those from this work and reference sequences obtained from both cyanolichens and free-living cyanobacteria. All positions containing deletions were removed. Evolutionary analyses were performed in the software MEGA7 (Kumar et al., 2016). Finally, the phylogenetic tree was edited in iTOL (Letunic and Bork, 2019).

Quantification of Genes Related to Phosphate Solubilization

Using the DNA obtained from the thallus, substrate, and soil samples, quantification of the bacterial markers for phosphate solubilization, including the quinoprotein glucose dehydrogenase gene (gcd), alkaline phosphatase gene (phoD), acid phosphatase gene (phoN), phytase gene (appA), and phosphonoacetaldehyde hydrolase gene (phnX), was performed using quantitative PCR (qPCR) on an ABI 7300 Real-time PCR System (Kurth et al., 2020) and specific primers for each gene (Bergkemper et al., 2016a). These primers were reported as specific for bacteria, except phoD, which amplified genes related to some fungi but not Peltigera (Bergkemper et al., 2016a). Therefore, hereinafter we refer to phosphate solubilizing bacteria (PSB), although we cannot discard the presence of fungal genes related to P turnover. In addition, the 16S rRNA gene was also quantified as a marker for bacterial abundance (Kurth et al., 2020).

After each qPCR run, a melting curve analysis was performed to verify the amplicon specificity. The quantification of the target gene was conducted by using serial dilutions of plasmid-encoding the target genes (10⁷ to 10¹ gene copies μl⁻¹). Details on the used standards were published previously (Bergkemper et al., 2016b). All samples were diluted 1:4 to avoid inhibition during amplification caused by co-extracted humic substances based on a previous dilution test (data not shown). The efficiency of the qPCR was calculated as Efficiency (%) = [10(−1/slope) − 1]. Values were always above 85%, while the R² of the standard curve was above 0.99.

Data Analysis

The abundance data of phosphate solubilization genes were log10 transformed and subjected to the Shapiro–Wilk test. The abundance of genes and physicochemical parameters were compared using one-way ANOVA and Tukey’s multiple comparison test. Non-metric multidimensional scaling (NMDS) analysis and pairwise Adonis test were performed based on Bray-Curtis distance, while redundancy analysis (RDA) was based on Euclidean distance and a stepwise selection of environmental variables using the ordiR2step function. All statistical tests and graphics were performed in R Studio (R version 3.5.3) using the vegan (Oksanen et al., 2017) and ggplot (Wickham, 2016) packages.

Results

Characterization of the Study Sites

For this study, we selected five sampling sites representing different tree cover types. The first transect is located in a native N. pumilio mature forest (S1); the second one is spanning the transition between the native mature forest and a native second-growth forest (S2); the third one is located in the second-growth forest (S3); the fourth one is in the transition between the native second-growth forest and a second-growth forest of N. pumilio with the presence of pine specimens (S4); and finally, the last transect is the native second-growth forest with the presence of pine specimens (S5).

The first method used to estimate the forest cover, the middle of the intercept point method, calculates the relative coverage of the canopy of a tree species at a certain height (in this case, 2 m). The sum of this variable can reach values greater than 1 since the canopies of the trees overlap. According to this method, N. pumilio had the highest relative tree canopy coverage in all five sites. The sites S2 and S3, with second-growth forests, showed a less relative coverage of the tree canopy of N. dombeyi. Finally, the presence of P. contorta displaced the relative coverage of N. dombeyi in sites S4 and S5 (Table 1). The second method, the point-centered quarter method, calculates relative importance as the sum of the relative density, relative frequency, and relative volume for each species, being a way of measuring tree dominance. According to this method, N. pumilio is the most dominant tree species in all the forest sites, but it decreases in S4 with a concomitant increase of N. dombeyi, while P. contorta increased its dominance from S4 to S5 (Table 1). In summary, from S1 to S5, there is a decrease in the coverage of N. pumilio and an increase in that of P. contorta. N. dombeyi was also present but with no apparent coverage gradient among the sites (Table 1).

TABLE 1

Table 1. Forest cover of the arboreal species present in the different sites (Nothofagus pumilio, Nothofagus dombeyi, and Pinus contorta).

The physicochemical characteristics of the soils of the five sites are summarized in Table 2. Total organic matter, P Bray, soluble Pm (Pm-H₂O), labile Pm (Pm-HCO₃⁻), and labile Po (Po-HCO₃⁻) were significantly higher in S5 than in S1. Notably, soluble Pm was about 15-fold higher, while P-Bray was 3-fold higher in S1 than S5. Conversely, moisture content and N-NO₃⁻ tend to significantly decrease from S1 to S5. Finally, pH was higher in S2 and S4 than in S3 and S5, while N-NH₄ and Po-NaOH were similar in all sites.

TABLE 2

Table 2. Physicochemical parameters of the soil samples from the different sites.

Abundance of Phosphate Solubilization Marker Genes

Five specimens of P. frigida thalli were selected from each of the study sites. The identification of the mycobionts was confirmed by a phylogenetic analysis based on the 28S rRNA gene (Supplementary Figure S2) and the ITS hypervariable region (ITS-HR; Miadlikowska et al., 2003; Supplementary Figure S3). Further, all specimens selected were associated with the same photobiont haplotype namely Nostoc sp. C01 (Zúñiga et al., 2015) according to a 16S rRNA gene-based phylogenetic analysis (Supplementary Figure S4).

From the samples of the three different microenvironments (i.e., P. frigida thalli, substrates, and the surrounding soils) collected from the five sites (i.e., S1 to S5), we measured by qPCR the relative abundance of one gene that codes for an enzyme that triggers mineral phosphate solubilization (gcd), four genes that codes for enzymes involved in organic phosphate solubilization (phoD, phoN, appA, phnX), and the 16S rRNA gene as a bacterial marker (Figure 1). The quinoprotein glucose dehydrogenase gene (gcd; Figure 1A) was the most abundant marker in all soil samples, and unexpectedly its abundance was ~40,000 times higher in soils than in substrates and ~ 3,000 times higher in soils than in thalli in all the study sites. However, there were no significant differences in the abundance of this marker gene comparing the same microenvironment among the sampling sites. The abundance of the alkaline phosphatase marker gene (phoD; Figure 1B) also did not show differences when comparing the sampling sites, and although in the microenvironments of the extreme sites (i.e., S1 and S5) its abundance did not show significant differences, in the soil samples of the intermediate sites its abundance was significantly lower than in substrates and thalli. On the contrary, the abundance of the acid phosphatase marker gene (phoN; Figure 1C) was significantly less abundant in the extreme sites (i.e., S1 and S5) than in the other sites, independent of the microenvironment. The abundance of the phytase marker gene (appA; Figure 1D) and the phosphonoacetaldehyde hydrolase marker gene (phnX; Figure 1E) had the highest variation in soils throughout the sites, but appA showed the highest abundance in thalli and phnX in substrates. Finally, the 16S rRNA marker gene (Figure 1F) abundance almost did not change among sites, except for the soil samples from S5 that presented the highest abundance for this microenvironment. However, the abundance of the bacterial marker gene was significantly higher in thalli than in soils and substrates from the intermediate sites (i.e., S2 to S4).

FIGURE 1

Figure 1. Abundance of phosphate solubilization markers per site. Each value is the mean ± SD of five replicates. Different lowercase and capital letters represent significant difference (p < 0.05; one-way ANOVA) between microenvironments and sites, respectively. (A) Quinoprotein glucose dehydrogenase gene (gcd). (B) Alkaline phosphatase gene (phoD). (C) Acid phosphatase gene (phoN). (D) Phosphonatase gene (appA). (E) Phosphonoacetaldehyde hydrolase gene (phnX). (F) 16S rRNA gene as marker for bacterial abundance.

Relation of Explanatory Variables With the Diversity of Phosphate Solubilizing Marker Genes

NMDS multivariate analysis of the diversity of phosphate solubilization marker genes shows a clear separation of soil samples from those of lichen-associated microenvironments (i.e., substrates and thalli) in all the sampling sites (Figure 2A; Table 3). On the contrary, the separation of substrate and thallus samples, although statistically significative, is partially overlapped, being more evident in the intermediate than in the extreme sampling sites (i.e., S1 and S5). When comparing the diversity in the different sampling sites (Figure 2B; Table 4), there are two well-defined groups, the first one including S2, S3 and S4, and the second one including S1 and S5. Notoriously, the 95% ellipse encompassing the points of S1 is the largest one and encloses the 95% ellipse of S5, which is noticeably smaller.

FIGURE 2

Figure 2. Non-metric multidimensional scaling (NMDS) multivariate analysis of the diversity of phosphate solubilization marker genes based on Bray-Curtis distance. Ellipses of 95% confidence are included for enclosing samples from each microenvironment (A; i.e., thallus, substrate, and soil) and sites (B; i.e., S1 to S5).

TABLE 3

Table 3. PERMANOVA (adonis function in R) to evaluate the effect of the microenvironment on the distribution of the variance of each sample type.

TABLE 4

Table 4. PERMANOVA (adonis function in R) to evaluate the effect of the site on the distribution of the variance of each sample type.

To determine the effect of the edaphic and tree cover parameters on the abundance of phosphate solubilizing genes, we performed redundancy analyses (RDA) and ordistep for selecting the variables that significantly contribute to the model. The RDA analysis of soil samples, which explained 60.1% of the total variance, corroborates that sites S2, S3, and S4 grouped closely together and that S1 is the most diverse of all sites, in contrast to S5, which is the less diverse (Supplementary Figure S5A). The variables that significantly contributed to this model were the relative abundance of tree canopies of N. pumilio, N dombeyi, and P. contorta, pH, and the labile fraction of inorganic mineral phosphorus (Pm-HCO₃⁻; Table 5). On the other hand, the variance explained by the RDA analyses for substrates (Supplementary Figure S5B) and thalli (Supplementary Figure S5C) are not robust enough (36.7 and 47.5% of the total variance explained, respectively) to appropriately assert patterns of phosphate solubilizing bacteria diversity. In these cases, the only significant explanatory variables were the relative abundance of N. pumilio tree canopies and the labile fraction of mineral phosphorus (Pm-HCO₃⁻), respectively (Table 5).

TABLE 5

Table 5. Stepwise selection of significant environmental factors explaining the abundance variation of phosphate solubilization and mineralization markers.

Discussion

We assessed the abundance of phosphate solubilizing bacteria (PSB) in three microenvironments related to the forest floor: (i) Peltigera frigida thalli, (ii) their associated substrate (i.e., the thallus subjacent soil), and (iii) the forest soil. These habitats were assessed in five sites with different tree cover in the Coyhaique National Reserve of the Chilean Patagonia characterized by temperate forests of Nothofagus species (Donoso, 1993; Veblen et al., 1996; Fajardo and González, 2009). We selected five sites representing: a native mature forest (S1), a native second-growth forest (S3), a native second-growth forest with the presence of pine specimens (S5), and the transition zones between these vegetative formations (S2 and S4, respectively). According to the methods used to calculate the forest coverage, N. pumilio was the dominating tree species with the greatest relative importance and coverage at all sites, and P. contorta was increasing its abundance in the last two sites (S4 and S5) with a decrease of N. dombeyi.

The Abundance of PSB in Forest Soils Is Affected by the Tree Cover Type

Changes in the soil P content are influenced by multiple and diverse factors. One of them are wildfires, which in Patagonian forests has been reported to have long-term effects on soil properties, such as decreases in organic C and N and increases in pH and extractable P (Alauzis et al., 2004). Another factor is the influence of trees on soil chemical properties through their root exudates and the loss of their leaves in deciduous species (such as N. pumilio), contributing to changes in pH and the relative content and availability of macronutrients (Binkley and Valentine, 1991; Augusto et al., 2000; Frouz et al., 2009; Iovieno et al., 2010). Additionally, enzymatic activities of phosphorus solubilizing microorganisms could also affect the soil P content, which not only responds to the soil P availability but also to the balance with other nutrients, mainly nitrogen (Sorkau et al., 2018). In fact, Margalef et al. (2017) showed that the availability of nitrogen in soil is the main factor explaining phosphatase activity in temperate climates.

The quantification of the PSB marker genes in the soil samples evidenced two noticeable results. The first one was the high abundance of gcd, which codes an enzyme that allows the direct oxidation of glucose to gluconic acid, one of the main mechanisms of extracellular Pm solubilization studied in bacteria (Sharma et al., 2013; Sindhu et al., 2014; An and Moe, 2016). Therefore, the solubilization through gluconic acid of Pm from mineral complexes found in soils would be important in the sampling sites. There are other organic acids such as malic, lactic, citric and oxalic acids, which can also solubilize Pm (Sharma et al., 2013), but they have been less studied in this regard because the enzymes associated with their production play different intracellular roles and could be not relevant in the P cycle (Grafe et al., 2018). The second one was the pattern of abundance of the acid phosphatase marker gene (phoN), which was lower in the extreme sampling sites (i.e., S1 and S5) and more similar to each other than to the abundance in the intermediate sites. This could be related, in part, to the higher values of available nitrogen (ammonium plus nitrate) measured in the extreme sampling sites regarding the transition sampling sites (i.e., S2 and S4). In addition, the intermediate sites (S2-S3-S4) correspond to sectors of young trees and, therefore, they have a more active growth than those in the mature forest (S1) and then a greater demand for phosphorus that could explain the higher abundance of this marker gene. The measured pH in the soil samples (5.3 to 5.6), suggests that the activity of the enzyme coded by phoN could be more relevant than the activity of the alkaline phosphatase (encoded by phoD), as has been shown for other forest soils (Staddon et al., 1998; Herbien and Neal, 2008).

On the other hand, the growth dynamics of the tree species present in the sampling sites is different. P. contorta is a fast-growing exotic conifer, which would demand higher amounts of nutrients, such as P, compared to the native species N. pumilio (Fajardo and Gundale, 2015). Our results are in accordance with a comparative study about the P flux in coniferous and deciduous forests of the northern hemisphere, where the concentration of P in the forest floor under coniferous trees was about three times higher than in deciduous forests (Sohrt et al., 2017). Since the P content of coniferous foliage is larger than in deciduous one (Sohrt et al., 2017) and litterfall plays a critical role in nutrient cycling (Attiwill and Adams, 1993), the tree cover could affect the abundance of phosphate solubilizing microorganisms in soils. Indeed, in our study, the tree cover type strongly affected the abundance of PSB in the forest soil, being the relative coverage of all three tree canopies identified as a significant driver. The most noteworthy result regarding the effect of the tree cover influence on the PSB marker genes is the reduced diversity detected in the site with the highest presence of pines (S5) compared with the site of native mature forest (S1), although in the former the number of bacteria is higher.

Other significant explanatory variables were the pH, reported as one of the main factors controlling bacterial communities in soils (Fierer and Jackson, 2006), and labile Pm, which could be the primary source of P for PSB (Sims and Pierzynski, 2005). pH was higher in the transition zones (S2 and S4) than in the native second-growth forests independently of the presence of pines (S3 and S5), but we found a significant increasing gradient of total organic matter and various phosphorus fractions from the native mature forest (S1) to the native second-growth forest with the presence of pine specimens (S5). Although the phosphorus content values in our samples were in the range of soils in southern Chilean forests (Redel et al., 2008), differences were still observed between the sites.

It is important to mention that other factors, such as rhizospheric exudates and the association with mycorrhizae, could affect soil P dynamics (Chen et al., 2008). Interestingly, all the tree species of our study associate with ectomycorrhizal fungi (Bradbury et al., 1998; Godoy and Marín, 2019), which are highly efficient in P turnover by combining a high capacity to explore soils, the formation of polyphosphates, and the activity of acid phosphatase and phytase enzymes (Chen et al., 2008). However, they are affected by forest anthropogenic disturbances such as fires and afforestation (Marín et al., 2017), which could explain our observations that some bacterial markers related to the solubilization of organic phosphate (e.g., phoN) were lower in the soils of the sites where ectomycorrhizae would be contributing more importantly to the P turnover (i.e., S1 and S5).

Our results reinforce the concept that microorganisms are susceptible to changes in environmental conditions and management practices; thus, they are good indicators for assessing changes in soil quality or soil recovery from disturbances (Quintanilla Pérez, 2007), such as the afforestation with exotic trees after wildfires. A deep understanding of microbial communities in forest habitats with different forest cover is then essential for predicting the response of forest ecosystems to changes in environmental conditions (Lladó et al., 2017).

Lichen Microhabitats Shape Their Associated Bacteria, Reducing the Impact of the Tree Cover on the Abundance of PSB

Thalli and substrates could be considered microenvironments with a high level of nutrient turnover (N, P or C), as biological soil crusts (Kurth et al., 2021) or the rhizosphere (Kuzyakov and Blagodatskaya, 2015). According to our results, the solubilization of Po from stable molecules such as phytate (through the enzyme coded by appA) or phosphonate (through the enzyme coded by phnX) seems to be more important in these microenvironments, respectively, than the solubilization of Pm, which would be more relevant in soils due to the higher abundance of the gcd marker gene. Phytate is accumulated by plants (Balaban et al., 2016), and to a lesser extent by microorganisms (Turner et al., 2021), as a P storage compound. The capacity to hydrolyse phytate has been reported as widespread among lichens, which potentially contributes to P capture from atmospheric deposits and plant leachates (Higgins and Crittenden, 2015). Phosphonates, on the other hand, can be found in membranes or excretions of microorganisms (Horsman and Zechel, 2017; Calcott et al., 2018) and groups of phosphonate biosynthetic genes have been reported in lichenized Nostoc (Gagunashvili and Andrésson, 2018). Therefore, considering the close contact between the thallus and the substrate in foliose lichens, as P. frigida, that the secondary metabolites produced by the lichen would diffuse to its substrate (Leiva et al., 2016), and that the main photobiont of Peltigera lichens were detected in both thalli and substrates (Zúñiga et al., 2017), the lichen-related microenvironments could be a rich medium in this type of compounds.

As an essential element, P plays a relevant role as a limiting macronutrient in the growth and activity of lichens. For example, it has been observed that a single immersion in phosphate solution can double the annual growth of the lichen Lobaria pulmonaria (McCune and Caldwell, 2009), as well as that fertilization with P, could increase nitrogen fixation in the cyanolichens Peltigera aphtosa and Peltigera polydactyla (Weiss et al., 2005). Considering the relevance of P in the development of lichens, different studies have been focused on the possible role that the lichen microbiota would be playing in functions related to phosphate solubilization processes. Indeed, several publications report potential PSB in lichen thalli (Bates et al., 2012; Grube et al., 2015; Sigurbjörnsdóttir et al., 2015). Thus, despite the lower abundance of organic phosphate solubilizers in comparison to mineral phosphate solubilizers, the former could play an essential role in recycling the old parts of the thallus (Grube et al., 2015), releasing phosphate that members of the lichen symbiosis could use. It is important to highlight that phosphatase and phytase activity have been reported in several species of lichens (Whiteford et al., 2010; Higgins and Crittenden, 2015), attributed to lichen enzymes being produced by the mycobiont. However, the primers chosen for this study were designed to amplify bacterial genes (Bergkemper et al., 2016a), which would indicate that the lichen microbiome is at least contributing to P turnover. However, since the lichens of this study are terricolous (i.e., their substrate is soil), the high abundance of the gcd gene, as the molecular marker of mineral phosphate solubilizers, is in accordance with analyses of metagenomes of soils where this is the most abundant gene of a set of 40 genes that encode enzymes related to the P cycle (Bergkemper et al., 2016b; Grafe et al., 2018).

Previous studies have found that thalli and substrates constitute microenvironments that differentially structure the associated microbial communities (Grube and Berg, 2009; Bates et al., 2011; Hodkinson et al., 2012; Ramírez-Fernández et al., 2014; Leiva et al., 2016). In particular, a recent study has shown that the bacterial community associated with P. frigida thalli is more than a mere extension of the microbiota of their substrates, with several bacterial being highly abundant in lichens but almost absent in substrates (Leiva et al., 2021). In our study, some markers of P turnover were more abundant in thalli, which could be related to the fact that part of bacteria related to this function would be recruited from the substrates and would later be enriched in the thalli. This proposition is related to Cardinale et al. (2012), who suggests that those microbial groups with specific functions within the microbiome would be acquired from the new environment colonized by lichen propagules because they are better adapted to local conditions.

In contrast to soils, PSB in substrates and thalli were only correlated with N. pumilio cover and labile Pm, respectively, and RDA analysis did not explain much of the variance. This supports our proposal that the microenvironments of lichens and their substrates act as an environmental buffer reducing the influence of forest cover composition on bacteria involved in P turnover, and that other factors driving the association with lichens could be of importance, like microbiome specificity.

Conclusion

Although the phosphorus turnover in forests is a complex process involving not only different organisms, such as plants, fungi, and bacteria, but also their interactions, altogether our results show that the diversity of genes related to phosphate solubilizing bacteria in forest soils is strongly affected by the tree cover. Conversely, lichens could act as an environmental buffer reducing the influence of forest cover composition on their associated bacteria involved in phosphorus turnover.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author Contributions

CeM, MC, and JO: conceptualization. MS, MC, and JO: funding acquisition. CeM, DL, CaM, and JO: methodology. CeM, DL, and MG: formal analysis. JO: project administration. MS and JO: supervision. All authors revised and approved the final version of the manuscript.

Funding

This research was funded by ANID—National Fund for Scientific and Technological Development (FONDECYT) —1181510 and ANID—Millennium Science Initiative Program—ICN2021_002.

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.

Acknowledgments

The authors want to sincerely thank J. L. Parraguez for his valuable fieldwork assistance. Permits for sampling at the Coyhaique National Reserve were obtained from CONAF (Corporación Nacional Forestal). CeM acknowledges the support provided by the Vicerrectoría de Asuntos Académicos of the Universidad de Chile to finance her internship at the Helmholtz Zentrum München.

Supplementary Material

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

References

Alauzis, M. V., Mazzarino, M. J., Raffaele, E., and Roselli, L. (2004). Wildfires in NW Patagonia: long-term effects on a Nothofagus forest soil. For. Ecol. Manag. 192, 131–142. doi: 10.1016/j.foreco.2003.11.014

CrossRef Full Text | Google Scholar

Almendras, K., García, J., Carú, M., and Orlando, J. (2018a). Nitrogen-fixing bacteria associated with Peltigera cyanolichens and Cladonia chlorolichens. Molecules 23:3077. doi: 10.3390/molecules23123077

PubMed Abstract | CrossRef Full Text | Google Scholar

Almendras, K., Leiva, D., Carú, M., and Orlando, J. (2018b). Carbon consumption patterns of microbial communities associated with Peltigera lichens from a Chilean temperate forest. Molecules 23:2746. doi: 10.3390/molecules23112746

CrossRef Full Text | Google Scholar

An, R., and Moe, L. A. (2016). Regulation of pyrroloquinoline quinone-dependent glucose dehydrogenase activity in the model rhizosphere-dwelling bacterium Pseudomonas putida KT2440. Appl. Environ. Microbiol. 82, 4955–4964. doi: 10.1128/AEM.00813-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Aschenbrenner, I. A., Cernava, T., Erlacher, A., Berg, G., and Grube, M. (2017). Differential sharing and distinct co-occurrence networks among spatially close bacterial microbiota of bark, mosses and lichens. Mol. Ecol. 26, 2826–2838. doi: 10.1111/mec.14070

PubMed Abstract | CrossRef Full Text | Google Scholar

Attiwill, P. M., and Adams, M. A. (1993). Nutrient cycling in forests. New Phytol. 124, 561–582. doi: 10.1111/j.1469-8137.1993.tb03847.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Augusto, L., Ranger, J., Ponette, Q., and Rapp, M. (2000). Relationships between forest tree species, stand production and stand nutrient amount. Ann. For. Sci. 57, 313–324. doi: 10.1051/forest:2000122

CrossRef Full Text | Google Scholar

Balaban, N. P., Suleimanova, A. D., Valeeva, L. R., Chastukhina, I. B., Rudakova, N. L., Sharipova, M. R., et al. (2016). Microbial phytases and phytate: exploring opportunities for sustainable phosphorus management in agriculture. Am J Mol Biol 7, 11–29. doi: 10.4236/AJMB.2017.71002

CrossRef Full Text | Google Scholar

Bates, S. T., Berg-Lyons, D., Lauber, C. L., Walters, W. A., Knight, R., and Fierer, N. (2012). A preliminary survey of lichen associated eukaryotes using pyrosequencing. Lichenologist 44, 137–146. doi: 10.1017/S0024282911000648

CrossRef Full Text | Google Scholar

Bates, S. T., Cropsey, G. W. G., Caporaso, J. G., Knight, R., and Fierer, N. (2011). Bacterial communities associated with the lichen symbiosis. Appl. Environ. Microbiol. 77, 1309–1314. doi: 10.1128/AEM.02257-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergkemper, F., Kublik, S., Lang, F., Krüger, J., Vestergaard, G., Schloter, M., et al. (2016a). Novel oligonucleotide primers reveal a high diversity of microbes which drive phosphorous turnover in soil. J. Microbiol. Methods 125, 91–97. doi: 10.1016/j.mimet.2016.04.011

CrossRef Full Text | Google Scholar

Bergkemper, F., Schöler, A., Engel, M., Lang, F., Krüger, J., Schloter, M., et al. (2016b). Phosphorus depletion in forest soils shapes bacterial communities towards phosphorus recycling systems. Environ. Microbiol. 18, 1988–2000. doi: 10.1111/1462-2920.13188

PubMed Abstract | CrossRef Full Text | Google Scholar

Binkley, D., and Valentine, D. (1991). Fifty-year biogeochemical effects of green ash, white pine, and Norway spruce in a replicated experiment. For. Ecol. Manag. 40, 13–25. doi: 10.1016/0378-1127(91)90088-D

CrossRef Full Text | Google Scholar

Bradbury, S. M., Danielson, R. M., and Visser, S. (1998). Ectomycorrhizas of regenerating stands of lodgepole pine (Pinus contorta). Can. J. Bot. 76, 218–227. doi: 10.1139/b97-172

CrossRef Full Text | Google Scholar

Bray, R. H., and Kurtz, L. T. (1945). Determination of total, organic, and available forms of phosphorus in soils. Soil Sci. 59, 39–46. doi: 10.1097/00010694-194501000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Butler, O. M., Elser, J. J., Lewis, T., Mackey, B., and Chen, C. (2018). The phosphorus-rich signature of fire in the soil–plant system: a global meta-analysis. Ecol. Lett. 21, 335–344. doi: 10.1111/ele.12896

PubMed Abstract | CrossRef Full Text | Google Scholar

Calcott, M. J., Ackerley, D. F., Knight, A., Keyzers, R. A., and Owen, J. G. (2018). Secondary metabolism in the lichen symbiosis. Chem. Soc. Rev. 47, 1730–1760. doi: 10.1039/C7CS00431A

PubMed Abstract | CrossRef Full Text | Google Scholar

Carbone, I., White, J. B., Miadlikowska, J., Arnold, A. E., Miller, M. A., Magain, N., et al. (2019). T-BAS version 2.1: tree-based alignment selector toolkit for evolutionary placement of DNA sequences and viewing alignments and specimen metadata on curated and custom trees. Microbio. Resour. Announce. 8, e00328–e00319. doi: 10.1128/MRA.00328-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardinale, M., Grube, M., Castro, J. V., Müller, H., and Berg, G. (2012). Bacterial taxa associated with the lung lichen Lobaria pulmonaria are differentially shaped by geography and habitat. FEMS Microbiol. Lett. 329, 111–115. doi: 10.1111/j.1574-6968.2012.02508.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cernava, T., Müller, H., Aschenbrenner, I. A., Grube, M., and Berg, G. (2015). Analyzing the antagonistic potential of the lichen microbiome against pathogens by bridging metagenomic with culture studies. Front. Microbiol. 6:620. doi: 10.3389/fmicb.2015.00620

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. R., Condron, L. M., and Xu, Z. H. (2008). Impacts of grassland afforestation with coniferous trees on soil phosphorus dynamics and associated microbial processes: A review. For. Ecol. Manag. 255, 396–409. doi: 10.1016/j.foreco.2007.10.040

CrossRef Full Text | Google Scholar

Chhabra, S., Brazil, D., Morrissey, J., Burke, J. I., O’Gara, F., Dowling, N., et al. (2013). Characterization of mineral phosphate solubilization traits from a barley rhizosphere soil functional metagenome. Microbiology 2, 717–724. doi: 10.1002/mbo3.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickie, I. A., St John, M. G., Yeates, G. W., Morse, C. W., Bonner, K. I., Orwin, K., et al. (2014). Belowground legacies of Pinus contorta invasion and removal result in multiple mechanisms of invasional meltdown. AoB Plants 6, 1–15. doi: 10.1093/aobpla/plu056

PubMed Abstract | CrossRef Full Text | Google Scholar

do Nascimento, C. A. C., Pagliari, P. H., Schmitt, D., He, Z., and Waldrip, H. (2015). Phosphorus concentrations in sequentially fractionated soil samples as affected by digestion methods. Sci. Rep. 5:17967. doi: 10.1038/srep17967

PubMed Abstract | CrossRef Full Text | Google Scholar

Donoso, C. (1993). Bosques templados de Chile y Argentina. Variación, estructura y dinámica. First Edn. Chile: Editorial Universitaria.

Google Scholar

Fajardo, A., and González, M. E. (2009). Replacement patterns and species coexistence in an Andean Araucaria-Nothofagus forest. J. Veg. Sci. 20, 1176–1190. doi: 10.1111/j.1654-1103.2009.01117.x

CrossRef Full Text | Google Scholar

Fajardo, A., and Gundale, M. J. (2015). Combined effects of anthropogenic fires and land-use change on soil properties and processes in Patagonia, Chile. For. Ecol. Manag. 357, 60–67. doi: 10.1016/j.foreco.2015.08.012

CrossRef Full Text | Google Scholar

Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 103, 626–631. doi: 10.1073/pnas.0507535103

PubMed Abstract | CrossRef Full Text | Google Scholar

Finér, L., Helmisaari, H.-S., Lõhmus, K., Majdi, H., Brunner, I., Børja, I., et al. (2007). Variation in fine root biomass of three European tree species: beech (Fagus sylvatica L.), Norway spruce (Picea abies L. karst.), and scots pine (Pinus sylvestris L.). Plant Biosyst. 141, 394–405. doi: 10.1080/11263500701625897

CrossRef Full Text | Google Scholar

Frouz, J., Pižl, V., Cienciala, E., and Kalčík, J. (2009). Carbon storage in post-mining forest soil, the role of tree biomass and soil bioturbation. Biogeochemistry 94, 111–121. doi: 10.1007/s10533-009-9313-0

CrossRef Full Text | Google Scholar

Gagunashvili, A. N., and Andrésson, Ó. S. (2018). Distinctive characters of Nostoc genomes in cyanolichens. BMC Genomics 19, 1–18. doi: 10.1186/S12864-018-4743-5

CrossRef Full Text | Google Scholar

Gardes, M., and Bruns, T. D. (1993). ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts. Mol. Ecol. 2, 113–118. doi: 10.1111/j.1365-294X.1993.tb00005.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Godoy, R., and Marín, C. (2019). “Mycorrhizal studies in temperate rainforests of southern Chile,” in Mycorrhizal Fungi in South America. eds. M. Pagano and M. Lugo (Cham: Springer), 315–341.

Google Scholar

Goldstein, A. H. (1995). Recent progress in understanding the molecular genetics and biochemistry of calcium phosphate solubilization by gram negative Bacteria. Biol. Agric. Hortic. 12, 185–193. doi: 10.1080/01448765.1995.9754736

CrossRef Full Text | Google Scholar

González, I., Ayuso-Sacido, A., Anderson, A., and Genilloud, O. (2005). Actinomycetes isolated from lichens: evaluation of their diversity and detection of biosynthetic gene sequences. FEMS Microbiol. Ecol. 54, 401–415. doi: 10.1016/j.femsec.2005.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodall, D. W. (1952). Some considerations in the use of point quadrats for the analysis of vegetation. Australian J. Biol. Sci. 5, 1–41.

Google Scholar

Grafe, M., Goers, M., von Tucher, S., Baum, C., Zimmer, D., Leinweber, P., et al. (2018). Bacterial potentials for uptake, solubilization and mineralization of extracellular phosphorus in agricultural soils are highly stable under different fertilization regimes. Environ. Microbiol. Rep. 10, 320–327. doi: 10.1111/1758-2229.12651

PubMed Abstract | CrossRef Full Text | Google Scholar

Grimm, M., Grube, M., Schiefelbein, U., Zühlke, D., and Bernhardt, J. (2021). The lichens’ microbiota, still a mystery? Front. Microbiol. 12:623839., doi: 10.3389/fmicb.2021.623839

PubMed Abstract | CrossRef Full Text | Google Scholar

Grube, M., and Berg, G. (2009). Microbial consortia of bacteria and fungi with focus on the lichen symbiosis. Fungal Biology Rev. 23, 72–85. doi: 10.1016/j.fbr.2009.10.001

CrossRef Full Text | Google Scholar

Grube, M., Cernava, T., Soh, J., Fuchs, S., Aschenbrenner, I., Lassek, C., et al. (2015). Exploring functional contexts of symbiotic sustain within lichen-associated bacteria by comparative omics. ISME J. 9, 412–424. doi: 10.1038/ismej.2014.138

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawksworth, D. L., and Grube, M. (2020). Lichens redefined as complex ecosystems. New Phytol. 227, 1281–1283. doi: 10.1111/nph.16630

PubMed Abstract | CrossRef Full Text | Google Scholar

Hedley, M. J., Stewart, J. W. B., and Chauhan, B. S. (1982). Changes in inorganic and organic soil phosphorus fractions induced by cultivation practices and by laboratory incubations. Soil Sci. Soc. Am. J. 46, 970–976. doi: 10.2136/sssaj1982.03615995004600050017x

CrossRef Full Text | Google Scholar

Herbien, S. A., and Neal, J. L. (2008). Soil pH and phosphatase activity. Commun. Soil Sci. Plant Anal. 21, 439–456. doi: 10.1080/00103629009368244

CrossRef Full Text | Google Scholar

Higgins, N. F., and Crittenden, P. D. (2015). Phytase activity in lichens. New Phytol. 208, 544–554. doi: 10.1111/nph.13454

PubMed Abstract | CrossRef Full Text | Google Scholar

Hodkinson, B. P., Gottel, N. R., Schadt, C. W., and Lutzoni, F. (2012). Photoautotrophic symbiont and geography are major factors affecting highly structured and diverse bacterial communities in the lichen microbiome. Environ. Microbiol. 14, 147–161. doi: 10.1111/j.1462-2920.2011.02560.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hooper, D. U., Bignell, D. E., Brown, V. K., Brussard, L., Dangerfield, J. M., Wall, D. H., et al. (2000). Interactions between aboveground and belowground biodiversity in terrestrial ecosystems: patterns, mechanisms, and feedbacks. Bioscience 50, 1049–1061. doi: 10.1641/0006-3568(2000)050[1049:IBAABB]2.0.CO;2

CrossRef Full Text | Google Scholar

Horsman, G. P., and Zechel, D. L. (2017). Phosphonate biochemistry. Chem. Rev. 117, 5704–5783. doi: 10.1021/acs.chemrev.6b00536

PubMed Abstract | CrossRef Full Text | Google Scholar

Iovieno, P., Alfani, A., and Bååth, E. (2010). Soil microbial community structure and biomass as affected by Pinus pinea plantation in two Mediterranean areas. Appl. Soil Ecol. 45, 56–63. doi: 10.1016/j.apsoil.2010.02.001

CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurth, J. K., Albrecht, M., Karsten, U., Glaser, K., Schloter, M., and Schulz, S. (2020). Correlation of the abundance of bacteria catalyzing phosphorus and nitrogen turnover in biological soil crusts of temperate forests of Germany. Biol. Fertil. Soils 57, 179–192. doi: 10.1007/s00374-020-01515-3

CrossRef Full Text | Google Scholar

Kurth, J. K., Albrecht, M., Karsten, U., Glaser, K., Schloter, M., and Schulz, S. (2021). Correlation of the abundance of bacteria catalyzing phosphorus and nitrogen turnover in biological soil crusts of temperate forests of Germany. Biol. Fertil. Soils 57, 179–192. doi: 10.1007/s00374-020-01515-3

CrossRef Full Text | Google Scholar

Kuzyakov, Y., and Blagodatskaya, E. (2015). Microbial hotspots and hot moments in soil: concept & review. Soil Biol. Biochem. 83, 184–199. doi: 10.1016/j.soilbio.2015.01.025

CrossRef Full Text | Google Scholar

Leiva, D., Clavero-León, C., Carú, M., and Orlando, J. (2016). Intrinsic factors of Peltigera lichens influence the structure of the associated soil bacterial microbiota. FEMS Microbiol. Ecol. 92, 1–9. doi: 10.1093/femsec/iw178

CrossRef Full Text | Google Scholar

Leiva, D., Fernández-Mendoza, F., Acevedo, J., Carú, M., Grube, M., and Orlando, J. (2021). The bacterial community of the foliose macro-lichen Peltigera frigida is more than a mere extension of the microbiota of the subjacent substrate. Microb. Ecol. 81, 965–976. doi: 10.1007/s00248-020-01662-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Letunic, I., and Bork, P. (2019). Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256–W259. doi: 10.1093/nar/gkz239

PubMed Abstract | CrossRef Full Text | Google Scholar

Liba, C. M., Ferrara, F. I. S., Manfio, G. P., Fantinatti-Garboggini, F., Albuquerque, R. C., Pavan, C., et al. (2006). Nitrogen-fixing chemo-organotrophic bacteria isolated from cyanobacteria-deprived lichens and their ability to solubilize phosphate and to release amino acids and phytohormones. J. Appl. Microbiol. 101, 1076–1086. doi: 10.1111/j.1365-2672.2006.03010.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lladó, S., López-Mondéjar, R., and Baldrian, P. (2017). Forest soil bacteria: diversity, involvement in ecosystem processes, and response to global change. Microbiol. Mol. Biol. Rev. 81, 1–27. doi: 10.1128/mmbr.00063-16

CrossRef Full Text | Google Scholar

Margalef, O., Sardans, J., Fernández-Martínez, M., Molowny-Horas, R., Janssens, I. A., Ciais, P., et al. (2017). Global patterns of phosphatase activity in natural soils. Sci. Rep. 7, 1–13. doi: 10.1038/s41598-017-01418-8

CrossRef Full Text | Google Scholar

Marín, C., Godoy, R., Valenzuela, E., Schloter, M., Wubet, T., Boy, J., et al. (2017). Functional land-use change effects on soil fungal communities in Chilean temperate rainforests. J. Soil Sci. Plant Nutr. 17, 985–1002. doi: 10.4067/S0718-95162017000400011

CrossRef Full Text | Google Scholar

McCune, B., and Caldwell, B. A. (2009). Treatment doubles single phosphorus growth of lichen transplants cyanobacterial. Ecol. Soc. Am. 90, 567–570.

Google Scholar

Miadlikowska, J., and Lutzoni, F. (2000). Phylogenetic revision of the genus Peltigera (lichen-forming Ascomycota) based on morphological, chemical, and large subunit nuclear ribosomal DNA data. Int. J. Plant Sci. 161, 925–958. doi: 10.1086/317568

CrossRef Full Text | Google Scholar

Miadlikowska, J., Lutzoni, F., Goward, T., Zoller, S., and Posada, D. (2003). New approach to an old problem: incorporating signal from gap-rich regions of ITS and rDNA large subunit into phylogenetic analyses to resolve the Peltigera canina species complex. Mycologia 95, 1181–1203. doi: 10.2307/3761919

CrossRef Full Text | Google Scholar

Mitchell, K. (2010). Quantitative analysis by the point-centered quarter method. arXiv preprint. Available at: arXiv:1010.3303 (Accessed September 10, 2019).

Google Scholar

Murphy, J., and Riley, J. P. (1962). A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27, 31–36. doi: 10.1016/S0003-2670(00)88444-5

CrossRef Full Text | Google Scholar

Nash, T. H. (2008). Lichen Biology. 2nd ed. United Kingdom: Cambridge University Press, doi: 10.1017/CBO9780511790478

CrossRef Full Text | Google Scholar

Oksanen, J., Kindt, R., Legendre, P., O’Hara, B., Simpson, G. L., Solymos, P. M., et al. (2017). Community Ecology Package. R package version 2.4–2. Community ecology package 2.5–6, 190. Available at: https://CRAN.R-project.org/package=vegan

Google Scholar

Osman, K. T. (ed.) (2013). “Chemical Properties of Forest Soils” in Forest Soils (Berlin: Springer), 45–61.

Google Scholar

Printzen, C., Fernández-Mendoza, F., Muggia, L., Berg, G., and Grube, M. (2012). Alphaproteobacterial communities in geographically distant populations of the lichen Cetraria aculeata. FEMS Microbiol. Ecol. 82, 316–325. doi: 10.1111/j.1574-6941.2012.01358.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Quintanilla Pérez, V. (2007). Perturbaciones a la vegetación nativa por grandes fuegos de 50 años atrás, en bosques Nordpatagónicos. Caso de estudio en Chile meridional. Anales de geografía de la Universidad Complutense.

Google Scholar

Ramírez-Fernández, L., Zúñiga, C., Caru, M., and Orlando, J. (2014). Environmental context shapes the bacterial community structure associated to Peltigera cyanolichens growing in Tierra. World J. Microbiol. Biotechnol. 30, 1141–1144. doi: 10.1007/s11274-013-1533-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramírez-Fernández, L., Zúñiga, C., Méndez, M. A., Carú, M., and Orlando, J. (2013). Genetic diversity of terricolous Peltigera cyanolichen communities in different conservation states of native forest from southern Chile. Int. Microbiol. 16, 243–252. doi: 10.2436/20.1501.01.200

PubMed Abstract | CrossRef Full Text | Google Scholar

Redel, Y., Rubio, R., Godoy, R., and Borie, F. (2008). Phosphorus fractions and phosphatase activity in an Andisol under different forest ecosystems. Geoderma 145, 216–221. doi: 10.1016/j.geoderma.2008.03.007

CrossRef Full Text | Google Scholar

Richardson, A. E., Hadobas, P. A., Hayes, J. E., O’hara, C. P., and Simpson, R. J. (2001). Utilization of phosphorus by pasture plants supplied with myo-inositol hexaphosphate is enhanced by the presence of soil micro-organisms. Plant Soil 229, 47–56. doi: 10.1023/A:1004871704173

CrossRef Full Text | Google Scholar

Rodríguez, H., Fraga, R., Gonzalez, T., and Bashan, Y. (2006). Genetics of phosphate solubilization and its potential applications for improving plant growth-promoting bacteria. Plant Soil 287, 15–21. doi: 10.1007/s11104-006-9056-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadzawka, R. A., Carrasco, R. M., Grez, Z. R., Mora, G. M., Flores, P. H., and Neaman, A. (2006). Métodos de análisis recomendados para los suelos de Chile. Santiago: Instituto de Investigaciones Agropecuarias.

Google Scholar

Sashidhar, B., and Podile, A. R. (2010). Mineral phosphate solubilization by rhizosphere bacteria and scope for manipulation of the direct oxidation pathway involving glucose dehydrogenase. J. Appl. Microbiol. 109, 1–12. doi: 10.1111/j.1365-2672.2009.04654.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Selbmann, L., Zucconi, L., Ruisi, S., Grube, M., Cardinale, M., and Onofri, S. (2010). Culturable bacteria associated with Antarctic lichens: affiliation and psychrotolerance. Polar Biol. 33, 71–83. doi: 10.1007/s00300-009-0686-2

CrossRef Full Text | Google Scholar

Sharma, S. B., Sayyed, R. Z., Trivedi, M. H., and Gobi, T. a (2013). Phosphate solubilizing microbes: sustainable approach for managing phosphorus deficiency in agricultural soils. Springerplus 2:587. doi: 10.1186/2193-1801-2-587

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigurbjörnsdóttir, M. A., Andrésson, Ó. S., and Vilhelmsson, O. (2015). Analysis of the Peltigera membranacea metagenome indicates that lichen-associated bacteria are involved in phosphate solubilization. Microbiology (NY) 161, 989–996. doi: 10.1099/mic.0.000069

CrossRef Full Text | Google Scholar

Sigurbjörnsdóttir, M. A., Andrésson, Ó. S., and Vilhelmsson, O. (2016). Nutrient scavenging activity and antagonistic factors of non-photobiont lichen-associated bacteria: a review. World J. Microbiol. Biotechnol. 32, 1–11. doi: 10.1007/s11274-016-2019-2

CrossRef Full Text | Google Scholar

Sims, T. J., and Pierzynski, G. M. (2005). Chemistry of phosphorus in soils. Soil Sci. Soc. Am. J. 677, 151–192. doi: 10.2136/sssabookser8.c2

CrossRef Full Text | Google Scholar

Sindhu, S. S., Phour, M., Choudhary, S. R., and Chaudhary, D. (2014). “Phosphorus cycling: prospects of using rhizosphere microorganisms for improving phosphorus nutrition of plants,” in Geomicrobiology and Biogeochemistry. eds. N. Parmar and A. Singh (Berlin: Springer), 199–237.

Google Scholar

Šnajdr, J., Dobiášová, P., Urbanová, M., Petránková, M., Cajthaml, T., Frouz, J., et al. (2013). Dominant trees affect microbial community composition and activity in post-mining afforested soils. Soil Biol. Biochem. 56, 105–115. doi: 10.1016/j.soilbio.2012.05.004

CrossRef Full Text | Google Scholar

Sohrt, J., Lang, F., and Weiler, M. (2017). Quantifying components of the phosphorus cycle in temperate forests. Wiley interdisciplinary reviews. WaterSA 4:e1243. doi: 10.1002/wat2.1243

CrossRef Full Text | Google Scholar

Sorkau, E., Boch, S., Boeddinghaus, R. S., Bonkowski, M., Fischer, M., Kandeler, E., et al. (2018). The role of soil chemical properties, land use and plant diversity for microbial phosphorus in forest and grassland soils. J. Plant Nutr. Soil Sci. 181, 185–197. doi: 10.1002/jpln.201700082

CrossRef Full Text | Google Scholar

Staddon, W. J., Duchesne, L. C., and Trevors, J. T. (1998). Acid phosphatase, alkaline phosphatase and arylsulfatase activities in soils from a jack pine (Pinus banksiana lamb.) ecosystem after clear-cutting, prescribed burning, and scarification. Biol. Fertil. Soils 27, 1–4. doi: 10.1007/S003740050390

CrossRef Full Text | Google Scholar

Steubing, L., Godoy, R., and Alberdi, M. (2002). Métodos en Ecología Vegetal. 1° Edición. Santiago de Chile: Editorial Universitaria.

Google Scholar

Turner, D., Kropinski, A. M., and Adriaenssens, E. M. (2021). A roadmap for genome-based phage taxonomy. Viruses 13, 1–10. doi: 10.3390/v13030506

CrossRef Full Text | Google Scholar

Veblen, T. T., Donoso, C., Kitzberger, T., and Rebertus, A. J. (1996). Ecology of southern Chilean and Argentinean Nothofagus forests. The ecolog. biogeograp. Nothofagus forests. 10, 93–353.

Google Scholar

Vilgalys, R., and Hester, M. (1990). Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. J. Bacteriol. 172, 4238–4246. doi: 10.1128/jb.172.8.4238-4246.1990

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiss, M., Hobbie, S. E., and Gettel, G. M. (2005). Contrasting Responses of nitrogen-fixation in Arctic lichens to experimental and ambient nitrogen and phosphorus availability. Arct. Antarct. Alp. Res. 37, 396–401. doi: 10.1657/1523-0430(2005)037[0396:CRONIA]2.0.CO;2

CrossRef Full Text | Google Scholar

White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR protocols: a guide to methods and applic. 18, 315–322. doi: 10.1016/B978-0-12-372180-8.50042-1

CrossRef Full Text | Google Scholar

Whiteford, E., Crittenden, P. D., Hogan, E., Minnullina, G., Smith, R. I., and Crittenden, P. D. (2010). Effects of nitrogen enrichment on phosphatase activity and nitrogen: phosphorus relationships in Cladonia portentosa. New Phytol. 186, 911–925. doi: 10.1111/j.1469-8137.2010.03222.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag.

Google Scholar

Wilmotte, A., Van der Auwera, G., and De Wachter, R. (1993). Structure of the 16 S ribosomal RNA of the thermophilic cyanobacterium Chlorogloeopsis HTF ('Mastigocladus laminosus HTF’) strain PCC7518, and phylogenetic analysis. FEBS Lett. 317, 96–100. doi: 10.1016/0014-5793(93)81499-P

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, K., Penttinen, P., Zhang, X., Ao, X., Liu, M., Yu, X., et al. (2014). Maize rhizosphere in Sichuan, China, hosts plant growth promoting Burkholderia cepacia with phosphate solubilizing and antifungal abilities. Microbiol. Res. 169, 76–82. doi: 10.1016/j.micres.2013.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zúñiga, C., Leiva, D., Carú, M., and Orlando, J. (2017). Substrates of Peltigera lichens as a potential source of cyanobionts. Microb. Ecol. 74, 561–569. doi: 10.1007/s00248-017-0969-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zúñiga, C., Leiva, D., Ramírez-Fernández, L., Carú, M., Yahr, R., and Orlando, J. (2015). Phylogenetic diversity of Peltigera cyanolichens and their photobionts in southern Chile and Antarctica. Microb. environ./JSME 30, 172–179. doi: 10.1264/jsme2.ME14156

PubMed Abstract | CrossRef Full Text | Google Scholar