Unaltered fungal community after fire prevention treatments over widespread Mediterranean rockroses (Halimium lasianthum)

Scientific Reports volume 13, Article number: 608 (2023) Cite this article

898 Accesses

5 Altmetric

Metrics details

Mediterranean ecosystems are frequently invaded by pyrophytic scrubs such as Halimium lasianthum that colonize areas traditionally used by livestock. A diverse fungal community is associated with this kind of vegetation, playing an important ecological role in these ecosystems. However, uncontrolled expansion of these shrubs considerably increases the risk of wildfires in these stands and, hence, fire-prevention treatments are needed. To investigate the long-term effects of two different forest-fire-prevention treatments on the soil fungal community, we analyzed these communities 9 years after prescribed burning or mechanical shredding were carried out in scrubland dominated by H. lasianthum. Neither of the fire-prevention treatments had a negative long-term effect on the abundance or richness of ectomycorrhizal fungi. However, saprotrophs and lichenized fungi experienced negative effects. Soil fertility significantly affected the distribution of fungi according to their functional groups, and pH was the most influential variable in terms of the distribution of edible species. Our findings indicate that forest management practices to prevent forest fires does not negatively affect the fungal community in the long-term, but for lichens and decomposers. Moreover, prescribed burning is suggested as a more economical way of reducing the risk of wildfires without affecting the ecology of the fungal community.

Forest fire is one of the major factor in forest resource degradation and affect fungal communities, primarily through fuel accumulation on the forest floor1. The need to reduce fire risk in the forests system while minimizing adverse impacts to the resource is an important challenge for forests managers, especially in fire-prone areas. In this regard, to prevent the negative effects of fires, fuel reduction treatments (mechanical and prescribed burning) are commonly used in the forests systems2 as they are also used as part of habitat management or ecological restoration3,4. Of these technics, the mechanical mastication is attractive as a fire surrogate to reduce fuels because it is cost effective and efficient5, while the prescribed fires are usually applied at low intensity and severity and generally does not have a significant impact on the forest resources4. However, the potential adverse effects of both the mechanical and prescribed fire treatments are a concern to forests managers because information about the ecological consequences are still lacking5.

Interest in the Cistaceae family is increasing due to its relationship with the increasing frequency of wildfires and its encroachment of former agricultural land. Aromatic compounds produced by the Cistaceae promote fire6,7 and the germination of Cistaceae seed is induced by heating, which happens in a wildfire, stimulating the rapid colonization of Cistaceae in burned areas8,9. However, this type of scrubland has also been generating interest because of its relationship with fungal succession and common mycorrhizal networks10,11. One of the less-studied genera among the Cistaceae family is Halimium spp., which is widely distributed throughout the Mediterranean region in degraded forest, open sites and even dry dunes12 from Portugal to Lebanon, and is even found as far north as Belgium13,14. Halimium is even found, like its Cistus relative, play a key ecological role as a dual-mycorrhizal plant helping in the fungal succession, particularly after fire disturbances14. The life cycle of this species is approximately 20 years15, which means that this genus shows rapid regrowth, enabling long- or mid-term studies of Halimium stands to be performed after a disturbance.

Traditionally, the expansion of this kind of scrubland was controlled by cattle grazing16, which was a natural, efficient, and productive way of preventing forest fires by reducing the amount of forest fuel available. However, because fewer cattle now graze in these areas, more fuel is accumulating in these ecosystems. The lack of management of Halimium stands promotes an increased risk of wildfire, specifically crown fire which can affect neighboring forested areas due to the vertical continuity between shrubs and trees17 which would take decades to the forests stands to recover. All these factors can influence the forest resources by reducing the vegetation cover, changing plant composition and damaging fungal communities10,18,19, the latter being the focus of this study.

Fuel-reduction practices, such as prescribed burning or mechanical shredding20, need to be implemented to modify the quantity and continuity of fuels, and reduce the risk of high-intensity and severe wildfires16. On the one hand, prescribed burning causes less damage than wildfires21, and is considered a low-cost option22. In this case, methods for controlling the burn must be well established so as to minimize the damage to the forest ecosystem23. On the other hand, mechanical shredding is a safer technique than prescribed burning; however, it is more expensive24 and it is difficult to undertake on steep slopes and stony terrains25. Previous studies undertaken also with Halimium lasianthum in the north-west of the Iberian Peninsula, in the Galician region, with different objectives have compared the effect of prescribed burning and mechanical shredding treatments. However, those studies observed no differences in plant mortality between treatments, or in vegetation cover and recovery, or even in species richness, diversity, and evenness24,26,27. Although prescribed burning has been reported to remove all the litter, which creates appropriate sites for rockrose germination28 as it removes the plant material to control advance of fire close to the ground, colonization by rockroses did not affect the fungal diversity or community composition29,30. Moreover, it is suggested31 that prescribed fire fosters a diversity of mycorrhizal and saprotrophic fungi owing to the increase in pyrophytic species.

Given the importance of Halimium spp. as a host for ectomycorrhizal and arbuscular mycorrhizal fungi14,32, the goal of our study was to analyze soil fungal communities in Halimium scrubland nine years after prescribed burning or mechanical shredding were carried out (i.e., midway through the lifespan of these Halimium plants) to determine the long-term effects of fire-prevention treatments. First, we hypothesized that a high diversity of fungi would be found associated with the Halimium scrubland14,15. In addition, contrary to other long-term studies33, we secondly hypothesized that neither fire-prevention treatment (Prescribed burning and mechanical shredding) would have a negative effect on the total fungal richness29, which would support the implementation of the underused practice of prescribed burning in Southern Europe35. Furthermore, we hypothesized that none of the guilds would be affected by prescribed fire or mechanical" treatments, particularly ectomycorrhizal (ECM) fungi given that previous studies have reported that shrub recovery is successful in the short-term in this area16,26. Finally, in this sense, our fourth hypothesis was that the total and guilds fungal composition would be more affected by edaphic variables than by wildfire prevention treatments. Thus, we assessed the hypotheses through the following specific objectives: (i) to describe the soil fungal richness and diversity of Halimium scrublands in the study area; (ii) to evaluate the impacts of fire prevention treatments on soil fungal richness and on guild status and; (iii) to determine whether the soil fungal community composition was governed by wildfire prevention treatments and by the environmental factors of the study area.

A total of 2225 fungal OTUs were classified into 13 phyla and 329 genera (Fig. 1). Although most of the OTUs were identified down to genus level, it was not possible to identify all OTUs down to species level. The soil fungal community was dominated by Ascomycota, which were mainly represented by Leotiomycetes (34% of the total OTUs identified) and Basidiomycota, which were mainly represented by Agaricomycetes (21% of the total OTUs identified) (Fig. 1A). The Krona charts showing taxonomic classification of the ectomycorrhizal trophic groups found in all the treatments is given as supplementary figures (Sup Fig. 1).

Krona charts showing fungal taxonomic classification in all treatments (prescribed, mechanical, and controls) at the phylum level (A) (name of phylum; number of OTUs; percentage); and (B) classification of fungal guilds based on the proportion of number of OTUs, including the taxa belonging to those guilds, based on Põlme et al.45.

In this study, 19 guilds were found, with saprotrophic guilds dominating the overall fungal community, including soil saprotrophs (13%) and litter saprotrophs (17%). There was also quite a large number of ECM fungal taxa (30 taxa: 3.7%) (Fig. 1B), including frequent taxa such as Pseudotomentella, Thelephora, and Sarcodon and other appreciated genera such as Boletus, Terfezia, Amanita, Russula or Rhyzpogon (Fig S1).

Species richness estimates calculated using Chao1 and ACE were higher than the observed richness values, indicating a low underestimation bias (Fig. 2). Accordingly, the observed species richness values are used from here onward to represent the soil fungal community.

Proportional species richness and richness estimation of total (A), trophic guilds (B) and fungal phyla (C,D) by extrapolating the singleton and doubleton (Chao1) values or rare species (ACE, abundance-based coverage estimator; abundance < 10 species) in soil samples collected 9 years after scrublands were subjected to prescribed burning, mechanical clearing or no treatment (Control). Different lowercase letters above bars indicate a significant difference between treatments (p < 0.05).

There were significant differences in terms of observed total fungal richness between treatments (p = 0.034) (Fig. 2A). The richness between the control and the burned plots differed significantly (p = 0.027), but not between mechanical shredding and control (p = 0.231). Moreover, Ascomycota phyla significantly differed in richness between treatments (p = 0.024), which also reflected differences in richness between the control and the burned plots (p = 0.026), but not between control and mechanical shredding (p = 0.08) and much less between burned and mechanically crushed stands (p = 0.851) (Fig. 2C). No differences between treatment plots were found for the other main phyla (Fig. 2D).

An unexpected lack of differences among treatments was found among most of the trophic guilds, including the ECM fungi (Fig. 2B). Lichenized and saprotrophic fungi showed significant differences between treatments (p < 0.01) (Table 1). However, Lichenized fungi was richer in the control plots than in the cleared (mechanical shredding) and burnt plots (p < 0.01, in both cases), but not between treated plots (p = 0.85). There was also significant difference between treatments in terms of saprotrophs fungi, with higher richness in the control plots than in the burnt plots (p < 0.001) or in the cleared plots (p = 0.03), but not between treated plots (p = 0.76).

The analysis of the total community showed a stress value of 0.161, indicating weak ties with no significant differences between treatments (p > 0.05). The distribution of the community did not show differences between treatments (Fig. 3A).

Non-metric multidimensional scaling (NMDS) of total fungal taxa (A) and trophic guilds (B) in soil samples collected 9 years after scrublands were subjected to fuel-reduction treatments. Treatments: control (dark green), burned (dark gray), cleared (brown). Ellipses represent the distribution of the fungal community related to the treatment. Isolines of the gradient of the C content (B) were also plotted on the NMDS ordinations using the ordisurf Model Gaussian function.

The analysis of guilds distribution showed weak ties with a stress for the model of 0.183 again. No significant differences were found between the treatments (p > 0.05), with pH being the unique influential environmental factor in our analysis (p = 0.026; r2 = 0.214) (Table 2). The distribution of trophic guilds (Fig. 3B) in the control was insignificantly different than the burning or mechanical shredding treatment, which tended to be more similar. The fungal communities in the control plots included more lichenized and endophytic fungal taxa, whereas ECM and arbuscular mycorrhizal fungal taxa were more associated with burnt plots. The soil fungal community in mechanically shredded plots, which was at the core of the shared community, appeared to occupy an intermediate position between burned and control plot communities.

Reports on soil fungal species from Halimium scrubland are limited. As to our knowledge, this study is the first to explain fungal communities in the Halimuum scrubland of the Edreiras mountains in Ourense province, North west Spain. As expected, Halimuum scrubland was associated with large number of fungal spies as we found a total of 2225 fungal OTUs, belonging to 329 genera and 19 different guilds, with saprotrophic composition dominating the overall fungal community. This was an indication of the importance of Halimuum scrubland in terms of diversity of soil fungal species in the study area. The present study also revealed occurrence of valuable ECM genera such as Boletus, Terfezia, Amanita, Russula and Rhizopogon among others. We also, in this study found that the total soil fungal richness seemed to be affected by prescribed burning treatment as compared to the control due to its effect on Ascomycetes, more specifically saprotrophic fungal richness, agreeing with other studies36. It is assumed that after a fire, the new ecological conditions lead to a rapid dominance of pyrophilous (Fire loving) ascomycetous species that later stabilize to form a Basidiomycota–Ascomycota community37, which is possibly due to the availability of a large amount of decomposable material following the fire38,39. The two treatments didn’t show a significant difference for the comparison of total fungal richness as expected in the second hypothesis, with the overall fungal community unaltered by the fire-prevention treatments in terms of composition at the large scale, and only insignificant altered in terms of richness.

In line with our hypothesis, the treatments barely affected the richness of trophic guilds. However, significant differences in the richness of two types of guilds (Lichenized and saprotrophic fungi) were observed between the control plots and the both types of treated plots. Like Ascomycota, saprotrophic fungi were similarly affected by both treatments. Different studies have reported contrasting effects of fire on saprotrophs: beneficial29, determinantal36, or neutral40. This suggests that the reaction of this guild to fire is context-dependent. A previous study performed at our experimental site recorded a decrease in the depth of the soil organic layer after prescribed burning24, which possibly inhibited the subsequent development of the saprotrophic community in the short-term given that it is often found to dominate the uppermost soil layers in forest ecosystems41. However, the ability of certain fungi to form heat-resistant structures contradicts this previous assumption given that the survival of thick-walled sclerotia or spore banks might give the surviving communities a competitive advantage36. It could be suggested that the low severity of our prescribed burning had yet bigger effect than expected over the soil organic material in comparison to the control. However, not enough to be differentiated from the mechanical shredding effect, which still leaves plant remains to be decomposed, bringing the possibility to an effect more related to other environmental factors than just the organic decomposable material. In the context of a fire, the woody biomass that normally hosts saprotrophic species is reduced in postfire forests42, which ultimately reduces the habitat that these organisms need given that it is a scrub-dominated vegetation. Furthermore oak-dominated areas subjected to prescribed burning produce litter with a higher C:N ratio than scrublands, leading to a slower decomposition rate43 and an altered saprotrophic community. Short-term experiments have shown that fire negatively alters the richness and abundance of saprotrophic fungi44, suggesting that the saprotrophic part of the community structure might require more than a decade to recover to its pre-fire state45. The fungal community associated with the mechanically cleared plots differed from those reported in previous studies34,46. Changes in the richness of saprotrophic fungi may be due to changes in the micro-environmental characteristics of the site, which alter the temperature, humidity, and light exposure of the area, negatively affecting the inhabiting community47,48. These alterations can also negatively affect cyanolichens and crustose sorediate species of lichens due to the reduced canopy cover, which is crucial for their survival49. However, this argument needs empirical research to justify so that we recommend further studies related to these issues. Lichenized taxa lost in this way might take years or decades to fully reestablish50. However, the resprouting behavior of plants from buried roots gives understory plants an advantage in terms of reassembling their fungal community in postfire forest ecosystems36,51. Previous studies performed at the experimental site24,27, as well as other studies36 have reported a high survival rate of resprouting plants following a fire. Therefore, we suspect that surviving roots buried in the soil might protect mycorrhizas and fungal endophytes from fires suggesting that vegetation communities consisting of plants able to resprout after fire could be considered more suitable for prescribed burning practices than other types of vegetation. Based on the assumption that there is a close ecological relationship between the vegetational and soil fungal communities, the relatively rapid recovery of the plant community after a fire treatment, as observed in earlier studies performed at the experimental site24,26,27, might have a beneficial effect on the recovery of the soil fungal community. In general, most of the phyla and guild richness levels were not affected by either of the fire-prevention treatments, potentially giving forest managers the option of using the fire-prevention treatment that is most appropriate for the characteristics and conditions of the area20. Although neither fire-prevention treatment appeared to have any severe long-term effects on the soil fungal communities assessed in this study, the deterioration of the soil fungal community in the short-term must be considered52. In particular, C- and N-stabilizing ascomycetes fungi might suffer from repetitive burnings over short time intervals, leading to a reduced capacity for C sequestration given that these fungi become more abundant in the long-term after fire events53.

Ecological factors have been observed to have non-significant effects on the taxa community29. However, in the case of guilds, the effects of treatments showed a more noticeable influence over the guilds distribution than the taxonomic one. The only environmental variable that showed significant effect over the guild distribution was the pH (Table 2), which is in agreement with others studies36,39,54,55, which could influence indirectly through plant diversity and mold prevalence, as through the solubility of minerals as P56. Although, soil pH levels recover to the pre-fire state through time when applying prescribed burning, however it is possible that our burned stands still experience a higher pH due to the ignition of organic acids57 during post-fire colonization, influencing the distribution. The long-term effects of prescribed burning on lichens tend to last for 10 to 15 years after treatment58 as well as alterations to micro-environmental conditions that affect them negatively49, such as canopy cover lost, relative to those growing in control stands. In the case of the endophytic population, fire-prevention treatments could lead to a change in guild status (i.e., endophyte to saprobic or pathogen)59, supporting the ‘body snatcher fungi’ hypothesis, which suggest that endophytes change their status after a disturbance adopting a new kind of functional role60.

In the long term, ECM fungi tend to recover after a fire of medium or low severity37 as the remaining C and N begin to becomes more recalcitrant, which negatively affects the saprotrophic fungi, thus, enabling ECM fungi to become more dominant39. Species colonization, such in the case of Rhizopogon spp., can be fostered by heat treatments61. This could be related to heat-resistant spores62, allowing mycorrhizal fungi to assist rapidly in the reforestation63, or the creation of new niches occupied by ECM fungi64, thanks to their spore production capacity65. Arbuscular mycorrhizal fungi also appeared to be associated with this forest type, which may be related to the existence of remaining prairie areas in burned stands66 or to the high resilience of arbuscular mycorrhizal fungi when facing fire disturbances in the long-term67. Furthermore, some studies have shown an increment in the relative abundance of mycorrhizal fungi based on analyses of spores or sequences68. As the differences in terms of composition are quite discrete and the overall fungal composition did not differ significantly between treatments, either fire-prevention treatment could be used. Other ecological or socioeconomic factors could also be a key factor when selecting a fire-prevention treatment69.

This study was carried out in Halimium stands in an area where the information in relation to fire prevention treatments and soil fungi is limited. Thus, this work could provide valuable scientific insight to promote the use of fire prevention treatments for the conservation fungal resources in the study area. In this study, in accordance with our expectation, many fungi are found in association with Halimium stands, indicating the importance of these stands in terms of fungal diversity. Also, no differences were found in terms of richness between treatments, and the taxa detected under the different treatments were almost identical. However, our expectation that the fungal guilds would not be significantly affected by the fire-prevention treatments was not entirely confirmed because the treatments did have a negative effect on the richness of saprobic ascomycete fungi and lichens; however, the richness of ECM fungi was not negatively affected by either of the treatments and was even promoted compared with the richness of the control stands (although not significantly). Furthermore, the fire-prevention treatments would not have a negative impact on the fungal community was confirmed as expected, with the overall fungal community unaltered by the fire-prevention treatments in terms of composition at the large scale. The promotion of a discrete mycorrhizal characterization was observed in burned plots in the community composition analysis, suggesting that the ability of ECM fungi to support forest succession was not affected by this fire-prevention treatment. This finding supports the idea of using non-spread prescribed burning treatments in certain regions of the Mediterranean basin where these H. lasianthum-dominated stands can be found. The lower cost of prescribed burning compared with mechanical clearing is another advantage, reducing the already high costs required for wildfire prevention. Moreover, further research should be carried out to investigate how the fungal response correlates with other ecological factors involved in this ecosystem, and possibly link it to the sustainable use and management of this ecosystem.

The study site is located in the Edreiras mountains (42° 8′ 02″ N–7° 26′ 17″ W, 1330 m a.s.l.), in Ourense province (NW Spain). It experiences Mediterranean climate conditions with an average rainfall of 1100 mm year−1 and a mean annual temperature of 10 °C. The soils are characterized as Alumi-umbric Regosols (FAO, 1998) developed on schists, and the sites, on average, have a 10% slope. The experimental area is mainly covered by heather, the main species being Erica australis L., together with Pterospartum tridentatum (L.) Wilk. and Halimium lasianthum (Lam.) Spach. Cattle and roe deer are very frequent in this area and are two important agents of ecosystem changes. In the past, prescribed burning was commonly motivated by pastoral considerations. For further information about this study area, see previous studies in this study area that investigated shrub recovery16 and resprouting24, seedling emergence and microbial responses70 to fuel reduction treatments.

Plots were established in the same local area to reduce the influence that differences in altitude71, aspect72, slope73, and climatic conditions74,75 may have on different fungal communities. The experiment involved the analysis of two different fire-prevention treatments to obtain different plot conditions, i.e., plots that had either been burned by prescribed burning or cleared by mechanical shredding in 2010, as well as control plots in which vegetation had been growing naturally since 2003 when the last wildfire affected the study area. A total of four experimental blocks were established. Each block was composed of nine rectangular plots, with 50 × 16 m2 plots (800 m2) per treatment (i.e., burned, cleared, and control plots), to reduce the influence of variations in vegetation and stoniness. A total of 36 plots were surveyed.

The following parameters were used to choose the meteorological window according to Fernandez et al.2: a temperature of 20 °C, a relative humidity range of 40 to 60%, and a low wind speed of 2 m/s. For the high dead fuels, the required moisture content was between 10 and 15%. For the organic layer, it was > 80%. The amount of soil moisture was almost at field capacity. With a strip width of 10 m, plots were burned using a strip head fire. The fires were captured on videotape, and the flame length was calculated by comparing it to stakes that were placed at various points across the plot and were known to be that length.

A steel-track tractor was used to masticate the above-ground biomass into a patchy layer of small diameter woody debris that was about 5 cm thick and stayed on the soil surface. In order to achieve treatment homogeneity, the operator made deliberate passes through the vegetation. Visual inspections revealed that the root systems and beneath soil were unaltered. The shrub clearance procedure involved manually cutting shrub with a trimmer from the base of the plants and removing it from the parcels. In our case, there was a greater than 95% clearance of vegetation.

After removing superficial stones and roots, 15 soil core samples were collected per plot between May and June 2019. The samples were mixed to form one pooled soil sample per plot to preserve the heterogeneity of the terrain, i.e., 36 pooled soil samples in total. The samples were placed in labeled plastic bags, transported to the laboratory by using an ice bag, and then processed within 24 h. Samples were air-dried and sieved through a 1-mm mesh sieve before analyzing them to remove rocky and oversized soil material. Each sample was subjected to genomic DNA analysis, which was performed using 0.25 g of soil which was stored at 4 °C for until subjected to DNA extraction. Soil physicochemical analysis was performed using two 20 g samples of each soil sample. Soil pH was determined by water-based suspension at 1:2.5. Also, to determine dry matter (%) following UNE-ISO 11465 rule and total phosphorus (P) using the Olsen methodology. For total carbon (C) and total nitrogen (N) contents (%) the Dumas methodology was used (Table 3).

A PowerSoil™ DNA Isolation Kit (MoBio Laboratories Inc., Carlsbad, CA, USA) was used to extract DNA from 0.25 g of soil per sample. PCR reactions of each sample were carried out in triplicate to minimize PCR biases. PCR reactions were performed in 20 μL reaction volumes containing 11.22 μL of Modified Quantization (MQ) water, 1.60 μL of DNA template, 2.00 μL of 10 × buffer, 1.40 μL of MgCl2 (50 mM), 1.60 μL of dNTPs (10 mM), 0.50 μL of bovine serum albumin (2%), 0.80 μL of reverse and forward primers (10 μM), and 0.08 μL of Platinum Taq polymerase (Invitrogen, Carlsbad, CA, USA). The following PCR conditions were used: an initial denaturation step at 94 °C for 3 min; followed by 35 cycles of 94 °C for 45 s, 50 °C for 1 min, and 72 °C for 1.5 min; and finally, one cycle of 72 °C for 10 min. PCR was used to amplify the ITS2 region (ca. 250 bp) of the nuclear ribosomal DNA repeat using primers fITS776 and ITS477, as described by Geml71. Sample-specific Multiplex Identification DNA-tags were used to label the ITS4 primer. A negative control comprising MQ water instead of DNA was included on each set of PCR replicates, which underwent PCR under the same experimental conditions and was shown to be amplicon-free on a gel. Sequencing was performed using an Illumina MiSeq platform (BaseClear BV, Leiden, the Netherlands).

Above-ground vegetation was surveyed by establishing a 50-m long vegetation transect in each plot. Plant cover and height were measured using a line interception methodology, as described by Kent 78.

We used cutadapt (Martin, 2013) to trim off poor-quality ends in both directions (3′ and 5′) using a quality criteria value of q = 15. The next step was to join both sequences of each sample using USEARCH v.10.0.24079 and cutadapt, with a minimum sequence length of 200 bp. Primers (ITS4 reverse and fITS7 forward) were trimmed and sequences with expected errors of > 1 were removed. Then, sequences were combined into a single sample and read count numbers were recorded to generate an operational taxonomic unit (OTU) showing the number of times that an OTU was detected in each sample. We assigned sequences to taxonomic groups based on their similarity to the curated UNITE database (version v.8.0 released on November 18th, 2018), which contains identified and unidentified sequences assigned to species hypothesis (SH) groups defined based on dynamic sequence similarity thresholds80. Excluding OTUs with < 70% similarity or with < 200 bp pairwise alignment length to a fungal sequence, a total of 2225 OTUs were obtained, representing a total of 292,098 quality-filtered sequences. Taxonomic studies generally only consider OTUs with > 95% identity; however, given that we were analyzing fungal communities from an ecological point of view, all the OTUs detected were considered to be of interest in this study.

Assignment of ecological guilds was performed using the PlutoF web workbench (https://plutof.ut.ee)81. OTUs with > 90% similarity to a fungal SH group with known ecological function were assigned to functional groups according to Põlme82. Some additional categories were assigned to the ectomycorrhizal (ECM) fungi following criteria used by Geml73, Põlme et al.45 and data published by Agerer83, Tedersoo and Smith84, and the DEEMY database (http://deemy.de).

The statistical significance of variable groups under different treatments was calculated using a Monte Carlo permutation test (499 permutations). Soil variables, fungal proportional richness, and relative abundance were compared across treatments using Linear Mixed Effects models (LME, p ≤ 0.05), which were developed using the package Nlme85, where site was defined as random and fire-prevention treatments were defined as fixed factors, and using Tukey’s HSD test. The R software environment was used for all these statistical analyses (version 3.5.3; R Development Core Team 2019). To avoid a downward-biased estimate of the species richness of the observed assemblage, two estimates, the Chao1 estimator and the abundance-based coverage estimator (ACE), were computed to extrapolate the singleton and doubleton values of rare species (abundance < 10 species), respectively, using the ‘estimateR’ function in the ‘vegan’ package86,87,88. The relationship between soil fungal composition and the edaphic variables was visualized using non-metric multidimensional scaling (NMDS) based on a Hellinger-transformed OTU and environmental scaled data matrix using the metaMDS function from the ‘vegan’ package. The effects of treatments were analyzed using a permutational multivariate ANOVA (PerMANOVA) based on 999 permutations using the adonis function in the ‘vegan’ package89. The correlation of NMDS axes scores with explanatory variables was assessed using the envfit function in R. The correlations of the environmental variables with the composition of the soil fungal community were determined by the Mantel Test using Bray–Curtis distance on the total species matrix and Euclidean distance on environmental parameters.

Provisional Submission number GenBank: SUB11719157. We state that we dispose all the data necessary to be required for proving the rigorous study here provided.

Cairney, J. W. G. & Bastias, B. A. Influences of fire on forest soil fungal communities. Can. J. For. Res. 37, 207–215 (2007).

Article Google Scholar

Fernández, C., Vega, J. A. & Fonturbel, T. The effects of fuel reduction treatments on runoff, infiltration and erosion in two shrubland areas in the north of Spain. J. Environ. Manage. 105, 96–102 (2012).

Article Google Scholar

Reazin, C., Morris, S., Smith, J. E., Cowan, A. D. & Jumpponen, A. Fires of differing intensities rapidly select distinct soil fungal communities in a Northwest US ponderosa pine forest ecosystem. For. Ecol. Manage. 377, 118–127 (2016).

Article Google Scholar

Durán-Manual, F. et al. Prescribed burning in Pinus cubensis-dominated tropical natural forests: A myco-friendly fire-prevention tool. For. Syst. 31, e012 (2022).

Google Scholar

Busse, M. D., Hubbert, K. R., Fiddler, G. O., Shestak, C. J. & Powers, R. F. Lethal soil temperatures during burning of masticated forest residues. Int. J. Wildl. Fire 14, 267–276 (2005).

Article Google Scholar

Frazão, D. F. et al. Cistus ladanifer (Cistaceae): A natural resource in Mediterranean-type ecosystems. Planta 247, 289–300 (2018).

Article Google Scholar

Keeley, J. E., Bond, W. J., Bradstock, R. A., Pausas, J. G. & Rundel, P. W. Fire in mediterranean ecosystems. Fire Medit. Ecosyst. https://doi.org/10.1017/cbo9781139033091 (2011).

Article Google Scholar

Louro, R., Peixe, A. & Santos-silva, C. New insights on Cistus salviifolius L. micropropagation. J. Bot. Sci. 6, 10–14 (2017).

CAS Google Scholar

Valbuena, L., Tarrega, R. & Luis, E. Influence of heat on seed germination of Cistus laurifolius and Cistus ladanifer. J. Wildl. Fire 2, 15–20 (1992).

Article Google Scholar

Martín-Pinto, P., Vaquerizo, H., Peñalver, F., Olaizola, J. & Oria-De-Rueda, J. A. Early effects of a wildfire on the diversity and production of fungal communities in Mediterranean vegetation types dominated by Cistus ladanifer and Pinus pinaster in Spain. For. Ecol. Manage. 225, 296–305 (2006).

Article Google Scholar

Comandini, O., Contu, M. & Rinaldi, A. C. An overview of Cistus ectomycorrhizal fungi. Mycorrhiza 16, 381–395 (2006).

Article CAS Google Scholar

Zuzunegui, M. et al. Growth response of Halimium halimifolium at four sites with different soil water availability regimes in two contrasted hydrological cycles. Plant Soil 247, 271–281 (2002).

Article Google Scholar

Civeyrel, L. et al. Molecular systematics, character evolution, and pollen morphology of Cistus and Halimium (Cistaceae). Plant Syst. Evol. 295, 23–54 (2011).

Article Google Scholar

Leonardi, M., Furtado, A. N. M., Comandini, O., Geml, J. & Rinaldi, A. C. Halimium as an ectomycorrhizal symbiont: New records and an appreciation of known fungal diversity. Mycol. Prog. 19, 1495–1509 (2020).

Article Google Scholar

Oria-De-Rueda, J. A., Martín-Pinto, P. & Olaizola, J. Bolete productivity of cistaceous scrublands in northwestern Spain. Econ. Bot. 62, 323–330 (2008).

Article Google Scholar

Fernández, C., Vega, J. A. & Fonturbel, T. Does shrub recovery differ after prescribed burning, clearing and mastication in a Spanish heathland?. Plant Ecol. 216, 429–437 (2015).

Article Google Scholar

Ponte, E. D., Costafreda-Aumedes, S. & Vega-Garcia, C. Lessons learned from arson wildfire incidence in reforestations and natural stands in Spain. Forests 10, 1–18 (2019).

Article Google Scholar

Franco-Manchón, I., Salo, K., Oria-de-Rueda, J. A., Bonet, J. A. & Martín-Pinto, P. Are wildfires a threat to fungi in European Pinus forests? A case study of boreal and Mediterranean forests. Forests 10, 309 (2019).

Article Google Scholar

Mediavilla, O., Oria-de-Rueda, J. A. & Martin-Pinto, P. Changes in sporocarp production and vegetation following wildfire in a Mediterranean Forest Ecosystem dominated by Pinus nigra in Northern Spain. For. Ecol. Manage. 331, 85–92 (2014).

Article Google Scholar

Tomao, A., Antonio Bonet, J., Castaño, C. & De-Miguel, S. How does forest management affect fungal diversity and community composition? Current knowledge and future perspectives for the conservation of forest fungi. For. Ecol. Manage. 457, 117678 (2020).

Article Google Scholar

Espinosa, J., Rodríguez de Rivera, O., Madrigal, J., Guijarro, M. & Hernando, C. Predicting potential cambium damage and fire resistance in Pinus nigra Arn. ssp. salzmannii. For. Ecol. Manage. 474, 118372 (2020).

Article Google Scholar

Potts, J. B. & Stephens, S. L. Invasive and native plant responses to shrubland fuel reduction: Comparing prescribed fire, mastication and treatment season. Biol. Conserv. 142, 1657–1664 (2009).

Article Google Scholar

Agee, J. K. & Skinner, C. N. Basic principles of forest fuel reduction treatments. For. Ecol. Manage. 211, 83–96 (2005).

Article Google Scholar

Fernández, C., Vega, J. A. & Fonturbel, T. Fuel reduction at a Spanish heathland by prescribed fire and mechanical shredding: Effects on seedling emergence. J. Environ. Manage. 129, 621–627 (2013).

Article Google Scholar

Huggett, R. J., Abt, K. L. & Shepperd, W. Efficacy of mechanical fuel treatments for reducing wildfire hazard. For. Policy Econ. 10, 408–414 (2008).

Article Google Scholar

Fernández, C. & Vega, J. A. Shrub recovery after fuel reduction treatments and a subsequent fire in a Spanish heathland. Plant Ecol. 215, 1233–1243 (2014).

Article Google Scholar

Fernández, C., Vega, J. A. & Fonturbel, T. Does fire severity influence shrub resprouting after spring prescribed burning?. Acta Oecologica 48, 30–36 (2013).

Article ADS Google Scholar

Ellsworth, J. W., Harrington, R. A. & Fownes, J. H. Seedling emergence, growth, and allocation of Oriental bittersweet: Effects of seed input, seed bank, and forest floor litter. For. Ecol. Manage. 190, 255–264 (2004).

Article Google Scholar

Castaño, C. et al. Resistance of the soil fungal communities to medium-intensity fire prevention treatments in a Mediterranean scrubland. For. Ecol. Manage. 472, 118217 (2020).

Article Google Scholar

Anderson, I. C., Bastias, B. A., Genney, D. R., Parkin, P. I. & Cairney, J. W. G. Basidiomycete fungal communities in Australian sclerophyll forest soil are altered by repeated prescribed burning. Mycol. Res. 111, 482–486 (2007).

Article CAS Google Scholar

Hernández-Rodríguez, M. et al. Soil fungal community composition in a Mediterranean shrubland is primarily shaped by history of major disturbance, less so by current fire fuel reduction treatments. Unpublished (2015).

Oria de Rueda, J. A., Martín-Pinto, P. & Olaizola, J. Boletus edulis PRODUCTION IN XEROPHILIC AND PIROPHITIC SCHRUBS OF Cistus ladanifer AND Halimium lasianthum IN WESTERN SPAIN. in IV International Workshop on Edible Mycorrhizal Mushrooms (2005).

Hart, B. T. N., Smith, J. E., Luoma, D. L. & Hatten, J. A. Recovery of ectomycorrhizal fungus communities fifteen years after fuels reduction treatments in ponderosa pine forests of the Blue Mountains. Oregon. For. Ecol. Manage. 422, 11–22 (2018).

Article Google Scholar

Hernández-Rodríguez, M., Oria-de-Rueda, J. A., Pando, V. & Martín-Pinto, P. Impact of fuel reduction treatments on fungal sporocarp production and diversity associated with Cistus ladanifer L. ecosystems. For. Ecol. Manage. 353, 10–20 (2015).

Article Google Scholar

Fernandes, P. M. Scientific support to prescribed underburning in southern Europe: What do we know?. Sci. Total Environ. 630, 340–348 (2018).

Article ADS CAS Google Scholar

Day, N. J. et al. Wildfire severity reduces richness and alters composition of soil fungal communities in boreal forests of western Canada. Glob. Chang. Biol. 25, 2310–2324 (2019).

Article ADS Google Scholar

Salo, K., Domisch, T. & Kouki, J. Forest wildfire and 12 years of post-disturbance succession of saprotrophic macrofungi (Basidiomycota, Ascomycota). For. Ecol. Manage. 451, 117454 (2019).

Article Google Scholar

Zakaria, A. J. & Boddy, L. Mycelial foraging by Resinicium bicolor: Interactive effects of resource quantity, quality and soil composition. FEMS Microbiol. Ecol. 40, 135–142 (2002).

Article CAS Google Scholar

Hul, S. et al. Fungal community shifts in structure and function across a boreal forest fire chronosequence. Appl. Environ. Microbiol. 81, 7869–7880 (2015).

Article ADS Google Scholar

Vázquez-Veloso, A. et al. Prescribed burning in spring or autumn did not affect the soil fungal community in Mediterranean Pinus nigra natural forests. For. Ecol. Manage. 512, 120161 (2022).

Article Google Scholar

Lindahl, B. D. et al. Spatial separation of litter decomposition and mycorrhizal nitrogen uptake in a boreal forest. New Phytol. 173, 611–620 (2007).

Article CAS Google Scholar

Salomón, R., Rodríguez-Calcerrada, J., González-Doncel, I., Gil, L. & Valbuena-Carabaña, M. On the general failure of coppice conversion into high forest in Quercus pyrenaica stands: A genetic and physiological approach. Folia Geobot. 52, 101–112 (2017).

Article Google Scholar

Williams, R. J., Hallgren, S. W. & Wilson, G. W. T. Frequency of prescribed burning in an upland oak forest determines soil and litter properties and alters the soil microbial community. For. Ecol. Manage. 265, 241–247 (2012).

Article Google Scholar

Semenova-Nelsen, T. A., Platt, W. J., Patterson, T. R., Huffman, J. & Sikes, B. A. Frequent fire reorganizes fungal communities and slows decomposition across a heterogeneous pine savanna landscape. New Phytol. 224, 916–927 (2019).

Article Google Scholar

Oliver, A. K., Callaham, M. A. & Jumpponen, A. Soil fungal communities respond compositionally to recurring frequent prescribed burning in a managed southeastern US forest ecosystem. For. Ecol. Manage. 345, 1–9 (2015).

Article Google Scholar

Sanz-Benito, I., Mediavilla, O., Casas, A., Oria-de-Rueda, J. A. & Martín-Pinto, P. Effects of fuel reduction treatments on the sporocarp production and richness of a Quercus/Cistus mixed system. For. Ecol. Manage. 503, 119798 (2022).

Article Google Scholar

Santos-Silva, C., Gonçalves, A. & Louro, R. Canopy cover influence on macrofungal richness and sporocarp production in montado ecosystems. Agrofor. Syst. 82, 149–159 (2011).

Article Google Scholar

Lin, W. R. et al. The impacts of thinning on the fruiting of saprophytic fungi in Cryptomeria japonica plantations in central Taiwan. For. Ecol. Manage. 336, 183–193 (2015).

Article Google Scholar

Aragón, G., López, R. & Martínez, I. Effects of Mediterranean dehesa management on epiphytic lichens. Sci. Total Environ. 409, 116–122 (2010).

Article ADS Google Scholar

Hämäläinen, A., Kouki, J. & Lohmus, P. The value of retained Scots pines and their dead wood legacies for lichen diversity in clear-cut forests: The effects of retention level and prescribed burning. For. Ecol. Manage. 324, 89–100 (2014).

Article Google Scholar

Schimmel, J. & Granstrom, A. Fire severity and vegetation response in the boreal Swedish. Ecol. Soc. Am. 77, 1436–1450 (1996).

Google Scholar

Hinojosa, M. B., Albert-Belda, E., Gómez-Muñoz, B. & Moreno, J. M. High fire frequency reduces soil fertility underneath woody plant canopies of Mediterranean ecosystems. Sci. Total Environ. 752, 141877 (2021).

Article ADS CAS Google Scholar

Clemmensen, K. E. et al. Carbon sequestration is related to mycorrhizal fungal community shifts during long-term succession in boreal forests. New Phytol. 205, 1525–1536 (2015).

Article CAS Google Scholar

Tedersoo, L. et al. Disentangling global soil fungal diversity. Science 346, 1052–1053 (2014).

Article Google Scholar

Adamo, I. et al. Sampling forest soils to describe fungal diversity and composition. Which is the optimal sampling size in Mediterranean pure and mixed pine oak forests?. Fungal Biol. https://doi.org/10.1016/j.funbio.2021.01.005 (2021).

Article Google Scholar

Tedersoo, L. et al. Regional-scale in-depth analysis of soil fungal diversity reveals strong pH and plant species effects in northern Europe. Front. Microbiol. 11, 1953 (2020).

Article Google Scholar

Peay, K., Garbelotto, M. & Bruns, T. Evidence of dispersal limitation in soil microorganisms: Isolation reduces species richness on mycorrhizal tree islands. Ecology 91, 3631–3640 (2010).

Article Google Scholar

Koivula, M. & Vanha-Majamaa, I. Experimental evidence on biodiversity impacts of variable retention forestry, prescribed burning, and deadwood manipulation in Fennoscandia. Ecol. Process. 9, 1–22 (2020).

Article Google Scholar

Fox, S. et al. Fire as a driver of fungal diversity—A synthesis of current knowledge. Mycologia 00, 1–27 (2022).

Google Scholar

Raudabaugh, D. B. et al. Where are they hiding? Testing the body snatchers hypothesis in pyrophilous fungi. Fungal Ecol. 43, 100870 (2020).

Article Google Scholar

Izzo, A., Canright, M. & Bruns, T. D. The effects of heat treatments on ectomycorrhizal resistant propagules and their ability to colonize bioassay seedlings. Mycol. Res. 110, 196–202 (2006).

Article Google Scholar

Kipfer, T., Moser, B., Egli, S., Wohlgemuth, T. & Ghazoul, J. Ectomycorrhiza succession patterns in Pinus sylvestris forests after stand-replacing fire in the Central Alps. Oecologia 167, 219–228 (2011).

Article ADS Google Scholar

Glassman, S. I., Levine, C. R., Dirocco, A. M., Battles, J. J. & Bruns, T. D. Ectomycorrhizal fungal spore bank recovery after a severe forest fire: Some like it hot. ISME J. 10, 1228–1239 (2016).

Article Google Scholar

Buscardo, E. et al. Impact of wildfire return interval on the ectomycorrhizal resistant propagules communities of a Mediterranean open forest. Fungal Biol. 114, 628–636 (2010).

Article Google Scholar

Pringle, A., Vellinga, E. & Peay, K. The shape of fungal ecology: Does spore morphology give clues to a species’ niche?. Fungal Ecol. 17, 213–216 (2015).

Article Google Scholar

Zhang, K., Cheng, X., Shu, X., Liu, Y. & Zhang, Q. Linking soil bacterial and fungal communities to vegetation succession following agricultural abandonment. Plant Soil 431, 19–36 (2018).

Article ADS CAS Google Scholar

Xiang, X. et al. Arbuscular mycorrhizal fungal communities show low resistance and high resilience to wildfire disturbance. Plant Soil 397, 347–356 (2015).

Article CAS Google Scholar

Dove, N. C., Klingeman, D. M., Carrell, A. A., Cregger, M. A. & Schadt, C. W. Fire alters plant microbiome assembly patterns: Integrating the plant and soil microbial response to disturbance. New Phytol. 230, 2433–2446 (2021).

Article CAS Google Scholar

Fernandes, P. M. Fire-smart management of forest landscapes in the Mediterranean basin under global change. Landsc. Urban Plan. 110, 175–182 (2013).

Article Google Scholar

Fontúrbel, M. T., Fernández, C. & Vega, J. A. Prescribed burning versus mechanical treatments as shrubland management options in NW Spain: Mid-term soil microbial response. Appl. Soil Ecol. 107, 334–346 (2016).

Article Google Scholar

Geml, J. et al. Large-scale fungal diversity assessment in the Andean Yungas forests reveals strong community turnover among forest types along an altitudinal gradient. Mol. Ecol. 23, 2452–2472 (2014).

Article CAS Google Scholar

Chu, H. et al. Effects of slope aspects on soil bacterial and arbuscular fungal communities in a boreal forest in China. Pedosphere 26, 226–234 (2016).

Article Google Scholar

Geml, J. Soil fungal communities reflect aspect-driven environmental structuring and vegetation types in a Pannonian forest landscape. Fungal Ecol. 39, 63–79 (2019).

Article Google Scholar

Castaño, C. et al. Soil microclimate changes affect soil fungal communities in a Mediterranean pine forest. New Phytol. 220, 1211–1221 (2018).

Article Google Scholar

Collado, E. et al. Mushroom productivity trends in relation to tree growth and climate across different European forest biomes. Sci. Total Environ. 689, 602–615 (2019).

Article ADS CAS Google Scholar

Ihrmark, K., Bödeker, I. & Cruz-Martinez, K. New primers to amplify the fungal ITS2 region—Evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 82, 666–677 (2012).

Article CAS Google Scholar

White, T., Bruns, S., Lee, S. & Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications (eds Innis, M. A. et al.) 315–322 (Academic Press, 1990).

Google Scholar

Kent, M. Vegetation Description and Data Analysis: A Practical Approach (Wiley, 2011).

Google Scholar

Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010).

Article CAS Google Scholar

Kõljalg, U. et al. Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277 (2013).

Article Google Scholar

Abarenkov, K. et al. Plutof-a web based workbench for ecological and taxonomic research, with an online implementation for fungal its sequences. Evol. Bioinforma. 2010, 189–196 (2010).

Google Scholar

Põlme, S. et al. FungalTraits: A user-friendly traits database of fungi and fungus-like stramenopiles. Fungal Divers. 105, 1–16 (2020).

Article Google Scholar

Agerer, R. Fungal relationships and structural identity of their ectomycorrhizae. Mycol. Prog. 5, 67–107 (2006).

Article Google Scholar

Tedersoo, L. & Smith, M. E. Lineages of ectomycorrhizal fungi revisited: Foraging strategies and novel lineages revealed by sequences from belowground. Fungal Biol. Rev. 27, 83–99 (2013).

Article Google Scholar

Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & Team, R. C. Nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 3.1–128. http://CRAN.R-project.org/package=nlme (2016).

Chao, A. & Chiu, C. Species richness: Estimation and comparison. Wiley StatsRef https://doi.org/10.1002/9781118445112.stat03432.pub2 (2016).

Article Google Scholar

Chiu, C. H., Wang, Y. T., Walther, B. A. & Chao, A. An improved nonparametric lower bound of species richness via a modified good-turing frequency formula. Biometrics 70, 671–682 (2014).

Article MathSciNet MATH Google Scholar

Oksanen, J. et al. Vegan: Community Ecology Package. R package version 2.4–2. https://CRAN.R-project.org/package=vegan. (2017).

Oksanen, J., Blanchet, F., Kindt, R. & Al, E. vegan: Community Ecology Package. R package version 2.3–0. (2015).

Download references

This work was supported by resources provided by the Spanish Ministery of Science and innovation. We would also like to express our gratitude to the people involved in the fieldwork of this study.

This research work was partially supported by the project MYCOINFOR (PID2019-105188RB-I00), which was funded by the Spanish Ministry for Science and Innovation. Ignacio Sanz-Benito’s contract was funded by the National Program for Youth Guarantee by Junta Castilla y León, partially funded by Fondo Social Europeo.

Sustainable Forest Management Research Institute UVa-INIA, Avenida Madrid, S/N, 34004, Palencia, Spain

Pablo Martín-Pinto, María Santos, Juan Andrés Oria-de-Rueda, Aitor Vázquez-Veloso, Olaya Mediavilla & Ignacio Sanz-Benito

Centro de Investigación Forestal-Lourizán, Xunta de Galicia, P.O. Box. 127, 36080, Pontevedra, Spain

Cristina Fernández & Teresa Fontúrbel

University for Sustainable Development Eberswalde, Schickler Street 5, 16225, Eberswalde, Germany

Tim Stadler

Department of Vegetal Production and Natural Resources, University of Valladolid, Avenida Madrid, S/N, 34004, Palencia, Spain

Pablo Martín-Pinto

You can also search for this author in PubMed Google Scholar

Conceptualization: C.F., M.S., T.F., P.M.P.; Methodology: O.M., J.A. O.R., C.F., M.S., T.F., P.M.P.; Formal analysis and investigation: I.S.B., O.M., A.V.V., P.M.P.; Writing—original draft preparation: I.S.B., T.S., A.V.V., P.M.P.; Writing—review and editing: I.S.B., T.S., A.V.V., P.M.P., O.M.; Supervision: I.S.B., T.S., A.V.V., P.M.P., J.A.O.R.

Correspondence to Pablo Martín-Pinto.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Martín-Pinto, P., Fernández, C., Santos, M. et al. Unaltered fungal community after fire prevention treatments over widespread Mediterranean rockroses (Halimium lasianthum). Sci Rep 13, 608 (2023). https://doi.org/10.1038/s41598-023-27945-1

Download citation

Received: 20 June 2022

Accepted: 10 January 2023

Published: 12 January 2023

DOI: https://doi.org/10.1038/s41598-023-27945-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.