Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Assessment of functional and structural changes of soil fungal and oomycete communities in holm oak declined dehesas through metabarcoding analysis

Assessment of functional and structural changes of soil fungal and oomycete communities in holm... www.nature.com/scientificreports opeN Assessment of functional and structural changes of soil fungal and oomycete communities in Received: 21 December 2018 holm oak declined dehesas through Accepted: 15 March 2019 Published: xx xx xxxx metabarcoding analysis 1 1 2 Francisco J. Ruiz Gómez , Rafael M. Navarro-Cerrillo , Alejandro pérez-de-Luque , 3 4 3,4 Wolfgang oβwald , Andrea Vannini & Carmen Morales-Rodríguez Forest decline is nowadays a major challenge for ecosystem sustainability. Dehesas, which consists of savannah-like mediterranean ecosystems, are threatened by the holm oak decline in the south-west of Iberian peninsula. Phytophthora cinnamomi is considered the main agent of holm oak root rot, but little is known about the relationship between diversity of soilborne microbial community and the decline syndrome of holm oak. It would be hypothesized that the changes in the structure and functionality of the soil microbiome might influence tree health status through changes in richness and diversity of beneficial organisms such as mycorrhizal species, or fungal plant pathogens such as Fusarium spp. or Alternaria spp. t otal DNA of soil samples from declined oak dehesas was extracted and analyzed through metabarcoding techniques, to evaluate the specific composition and diversity of the fungal and oomycete communities and their relationship with the disease symptoms. the fungal community included a wide range of pathogens and abundance of ectomycorrhizal key taxa related with low defoliation degree. Phytophthora cinnamomi and Pythium spiculum did not appear among the most abundant oomycetes, nor were they related directly to defoliation levels. Moreover, a particular taxon belonging to the genus Trichoderma was strongly correlated with the scarcity of pathogenic Phytophthora spp. the diversity and composition of fungal and oomycete communities were related to the severity of the decline symptoms. the metabarcoding study of microbiome represents a powerful tool to develop biocontrol strategies for the management of the holm oak root rot. Alien-invasive species and global change are nowadays major threats for ecosystems sustainability around the world. It is possible to find examples of different die-back syndrome of woody species devastating natural or semi natural ecosystem around the world related with these species . In the last decades, different pathogenic Phytophthora spp. and Pythium spp. are emerging as important invasive pathogen species, considering the envi- 2,3 ronmental losses they cause . Some examples are Phytophthora ramorum Werres in Sudden Oak Death, ae ff ct- ing tanoak (Notholithoarpus densiflorus (Hook. & Arn.) Manos, Cannon & S.H.Oh) and several oak species in 4 5 California , Phytophthora cinnamomi Rands. in eucalyptus and Jarrah ecosystem in Australia , or Phytophthora cryptogea Pethybr. & Laff., Phytophthora cambivora (Petri) Buisman and P. cinnamomi causing ink disease in chestnuts worldwide . Among the most relevant diseases caused on European woody species by Phytophthora spp., root rot has been identified since the 1990’s ae ff cting oaks . Sessile and pedunculated oaks (Quercus petraea (Matt.), Liebl. Departamento de Ingeniería Forestal, Laboratorio de Ecofisiología de Sistemas Forestales ECSIFOR- ERSAF. Universidad de Córdoba. Campus de Rabanales, Crta. IV, km. 396, E-14071, Córdoba, Spain. Área de Genómica y Biotecnología, IFAPA, Centro Alameda del Obispo, Avda. Menéndez Pidal s/n, Apdo, 3092, 14080, Córdoba, Spain. Fachgebiet Pathologie der Waldbäume, Technische Universität München. Hans-Carl-von-Carlowitz-Platz 2, 85354, Freising, Germany. Department for Innovation in Biological, Agro-food and Forest systems (DIBAF) –University of Tuscia, Via San Camillo de Lellis snc, Viterbo, 01100, Italy. Correspondence and requests for materials should be addressed to F.J.R.G. (email: g72rugof@uco.es) Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 1 www.nature.com/scientificreports www.nature.com/scientificreports/ and Quercus robur L.) from North and Central Europe are threatened by several Phytophtora spp. , and in the Mediterranean basin, the decline of holm oak (Quercus ilex L., subsp. ballota in Spain and Portugal, and subsp. ilex in central Italy) and cork oak (Quercus suber L.) due to root rot were identified since the early 1980’s, ae ff cting 9–11 dehesa rangelands in Spain and Portugal, and wide forest areas in France and Italy . Dehesas represent ecosystems with relevant ecological and socio-economic importance covering 23% of the forested area in Spain (around 14% of the total area), accounting for over 3 mill ha in the Iberian Peninsula . In general, dehesa is a traditional land use in areas excluded from intensive agriculture, covering large extension of territory mostly characterized by a typical Mediterranean climate and soils usually shallow, acidic and deficient in nutrients . This agroecosystem is mainly composed by open woodlands of Quercus ilex and Q. suber, with an annual grassland understory, presenting canopy covers ranging from 5 to 60% . Dehesa ecosystem provides important socioeconomic benefits and many ecosystem services, such as shelter for a high degree of biodiversity or soil erosion prevention, being considered a first line barrier against desertification . Tree mortality of holm oak (Quercus ilex L.) associated with root rot is a serious threat to the sustainability 16,17 of these ecosystems. Phytophthora cinnamomi Rands. is considered the main agent of this disease , but the inu fl ence of other biotic and abiotic factors - such as the increasing frequency of extreme drought episodes and decreasing precipitation associated with global warming - is also recognized . Research efforts on holm oak decline have been focused on the interaction between the host plants and the pathogen – Q. ilex and P. cin- 19–21 namomi , and on their relationships with environmental and human factors, including socio-economic and 22,23 management strategies . Management strategies for the control of disease spreading are oen diffic ft ult to imple- ment because most of the dehesas yield is based on cattle livestock or hunting game . On the other hand, other treatments such as calcium soil amendments and phosphorous fertilization influenced tree health status but are 24,25 unable to suppress the soil inoculum . The application of sinigrin from Brassicaceae cultivars has been shown to be effective for pathogen suppression in controlled conditions experiments, but the biomass production nec- essary for an effective bio-amendment on the field is far greater than the productive ability of the dehesa soils . However, the relationship between the plant community and oomycetes and the spatial patterns of Pythium spp. and Phytophthora spp. were shown to be relevant factors influencing the dynamics of declining Mediterranean forests . Moreover, the implication of other Phytophthora spp. in the decline of Q. ilex and other Quercus spp. in 28–30 the Iberian Peninsula is recognized . In addition, little is known about the role of the soil fungal community in the decline syndrome of holm oak, or about the relationship between the latter and the presence of pathogenic species belonging to oomycetes. Die ff rent authors have highlighted the importance of the interactions among communities of soil taxa, mainly 31,32 focusing on bacterial and fungal communities . The soil fungal community interacts with other communities of biota in different ways, acting as commensal, pathogenic, antagonistic or mutualistic organisms. In the case of the plant community, its interaction is mainly associated with beneficial species such as ectomycorrhizal (EcM) and mycorrhizal (AM) fungi and other soil saprobes which improve soil fertility and structure and enhance plant defenses and health status . Ectomycorrhizal and endophytic species have been found to be related to changes 22,34–36 in the tree health status of several broad-leafed species, including Q. ilex . On the other hand, pathogenic relationships between fungal species and the plant community, including asymptomatic ones, are also related to disease expression in some cases . Oomycete communities are composed of a wide range of plant pathogens, some of which are host-specific and 38,39 others with a wide range of hosts. Their lifestyle can be either biotrophic, necrotrophic and hemibiotrophic . Phytophthora spp. communities have often been found to be the dominant in forest soils with a known history of Phytophthora-related diseases , but little is known about the relationship between the oomycete community and soil fungal or bacterial communities. Sapkota and Nicolaisen published the first study of the soil commu- nity diversity, including the taxa of fungi and oomycetes, associated with cavity spot of carrot, but in the litera- ture, there are no studies describing fungal diversity or the relationship between oomycetes and fungi in dehesas rangeland. e i Th solation of living cultures was the traditional way to assess oomycete or fungal diversity in soils, but this approach seems to be insufficient due to the time and specialized scientific knowledge required for phenotypic species identification , considering the potentially enormous number of different species – 1.5 to 3 million fun- gal species – in a high number of samples. In this sense, next-generation sequencing (NGS) techniques are an alternative methodology, able to detect the presence of a very high number of different taxa in a single sample: 40,44 successful examples include pyrosequencing and second-generation techniques, such as metabarcoding based 41,45,46 on Illumina de novo sequencing . The ability of metabarcoding techniques to characterize high-complex soil microbiome systems has represented an extraordinary improvement on the study of the plants and associated microbiota interaction. These techniques also oe ff r new opportunities for the deep assessment of microbiome ecology, given an important tool to optimize integrated control methods on management of plant diseases . The main objective of our work was to assess and describe the diversity and structure of fungal and oomy- cete communities in the holm oak dehesas rangeland ecosystems of Andalusia (southern Spain) ae ff cted by oak decline, using a metabarcoding approach, to analyse the correlations of structure and functionality of these com- munities between them, and with the defoliation level of Q. ilex. We hypothesize that the interaction between the microbiota diversity and functionality must be related with changes in tree health status, considering that fungi and oomycete are the main species involved in the holm oak decline due to root rot. To reach our main objective, the following specific objectives were addressed: i) to compare the fungal and oomycete diversity, and two defoli- ation levels in holm oak stands, in four different study zones; ii) to identify the predominant functional guilds of the fungal and oomycete communities and their relationships; and iii) to identify key taxa which might be related to holm oak decline. These results should contribute to a better understanding of the influence of soil community associations on tree health status and to the identification of important fungal and oomycete taxa related to holm oak decline in dehesas ecosystems. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 2 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 1. Mean relative OTUs abundance in 26 plots of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Abundance for fungi OTUs classified by “class” taxonomic level. (b) Abundance for oomycete OTUs classified by “genus” taxonomic level. Legends are referred to detail graph (left, “Other” pie chart) ordered by frequency. UI: Unidentified (for those OTUs which classification at the chosen taxonomic level was impossible). IS: Incertidae Sedis. The number of OTUs corresponding to each class (a) or genus (b) is showed in Supplementary Information, Table S1. Results otUs clustering and taxonomy. After bioinformatic treatment of the raw data, 1789 fungal OTUs (Operational Taxonomic Units) and 178 oomycete OTUs were confidently assigned to the genus or species level. Among the fungal taxa, Basidiomycota was the most abundant division (46.8%) followed by Ascomycota (45.1%) and Zygomycota (7.4%). e Th class Agaricomycetes was the most abundant (42.87%) and, together with Sordariomycetes, Eutriomycetes, Mortierellomyetes, and Pezizomycetes, represented the five dominant fungal classes (Fig. 1a). The class Agaricomycetes had its maximum relative abundance in the S.Nor zone (53.2%) and the lowest in Arac (39.4%). The opposite tendency was shown by Mortierellomycetes, which presented their minimum value in S.Nor (3.4%) and their maximum in Andv (11.2%) with a great value in Arac (6.8%) (Fig. 2a). No significant differences appeared for fungal classes in relative abundances between defoliation classes (Fig.  2b). Other relevant classes of Ascomycota were Dothideomycetes (5.36%) and Leotiomycetes (1.15%). Regarding oomycetes, species of the genera Phytophthora and Pythium spp. represented over 95% of the over- all clustered sequences, followed by Elongisporangium spp. (corresponding to E. undulatum), with almost 3% of the total sequence abundance, all of them corresponding to plant pathogen species (Fig. 1b). Moving eastward across the plots, the abundance of Phytophthora spp. increased and that of Pythium spp. decreased, with a sig- nificant inverse correlation between the relative abundances of the two taxa (Spearman ρ = −0.966; p < 0.001). However, none of the analysed oomycete taxa presented significant differences in abundance among the defolia- tion levels (Fig. 2c,d). The feature classification analysis of the microbiome selected 20 Oomycete OTUs and 349 Fungal OTUs as core features. Among the oomycetes, the most-abundant OTU, present in all the samples with impor- tant differences in abundance, was Phytophthora plurivora T. Jung & T.I. Burgess (oomycete OTU#1) - which accounted for 15.2% of the total frequencies, except in the case of Arac, where OTU#6 (Pythium paroecandrum Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 3 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 2. Relative abundance of oomycete and fungal taxa across study zones and mean defoliation class on the 26 studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. Vertical axis represents % of frequencies. Fungal abundance for zone (a) (%) and mean defoliation class (b) (%) is represented for taxonomic level 3 (class). Oomycete abundance for zone (c) (%) and mean defoliation class (d) (%) is represented for taxonomic level 6 (genus). Andv: Andévalo. Arac: Aracena. S.Nor: Sierra Morena & Sierra Norte. Pedr: Pedroches. UI: Unidentified taxon. Drechsler) was the most abundant, representing over 21% of the total frequencies in this zone (Table 1). The fourth-most-abundant oomycete was OTU#3, identified as Phytophthora quercina T. Jung & T.I. Burgess, whose presence was greatest in the Andv zone. Phytophthora cinnamomi (OTU#4) was the ninth-most-abundant OTU, accounting for 2.2% of the total sequences and identified in 18 of the 26 selected plots. Its presence was greatest in Pedr, where it was the third-most-abundant species (3.5% of the total frequency), and in Andv: in both cases, it showed lower abundance and frequency than P. plurivora and P. quercina. However, despite the low recovery rate on high-defoliated plots (C2, 38%), the abundance of P. cinnamomi in the high defoliated plots in which appeared was significantly high (31.8%). Two OTUs were identified as Pythium spiculum Paul. (OTU#21 and OTU#235), but they were present in only 21.2% of samples, corresponding to 10 plots, with low abundance (1.4% in total), and neither of them were classified as core features. Among the 349 fungal OTUs present in almost half of the soil samples, fungal OTUs #3 and #11, correspond- ing to Mortierella elongata Linnem., were the most abundant taxon, accounting over 6.1% of the total number of sequences (Table 1). The second most abundant OTU was OTU#2, Fusarium spp., present in all the samples ana- lyzed and representing over 3.4% of the total frequency. The third-most-abundant taxon was Russula praetervisa Sarnari (OTU#8), present and abundant in all the samples of Andv and S.Nor, but representing only 0.4% of the total abundance in Arac, and identified only in 75% of the samples of Pedr. Its abundance varied significantly with the mean defoliation (F = 3.38; p < 0.01), being higher in plots with low defoliation average (C1). Two OTUs belonging to the genus Trichoderma were classified among the 10 most-abundant species of the core biome: OTU#51, with its maximum frequency in samples from Arac, and OTU#15, representing over 1% of the total number of sequences in samples from Andv and S.Nor. Both OTUs had lower frequency rates in Pedr. Spearman’s correlation analysis showed a significant negative correlation between the abundance of Trichoderma spp. (fungal OTU#51) and that of 10 different OTUs of the oomycete core biome (Table  2), includ- ing most of the Phytophthora spp., except for OTUs #3 (P. quercina), #48 (Phytophthora sp.) and #12 (P. syringae). Regarding Pythium spp., only OTUs #6 (P. paroecandrum) and #2 (Pythium sp.) presented significant negative correlation. Despite its importance in our results, it was impossible to identify OUT#51 at species level because the well-conserved ITS1 region of Trichoderma spp. genus is unable to differentiate between species. Functional guilds. Of the OTUs confidently classified in an ecological guild (Table 3), the Soil-Plant Saprotroph guild was the most abundant in all cases (over 15%), except in S.Nor, where Plant Pathogen was the most-abundant guild (over 21%). However, when the frequencies of all the putative plant pathogenic OTUs were summed (including Plant Pathogen, Opportunistic Pathogenic Species and those OTUs classified in >1 Guild, which includes Plant Pathogen as one of those guilds), their abundance was greater than that of the rest of the guilds, in all the zones studied. Five functional guilds showed significant differences among zones, the Ectomycorrhizal (EcM) and Endophyte guilds exhibiting opposing trends, with the minimum frequency values of EcM species corresponding to the zone with maximum values for Endophyte. In the Pedroches zone (Pedr) the maximum values were found for the Endophyte, Animal Pathogen, Dung Saprotroph, and Soil-Plant Saprotroph functional guilds, with a Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 4 www.nature.com/scientificreports www.nature.com/scientificreports/ Location Average Defoliation Andevalo S. Aracena S. Norte Pedroches Low (C1) High (C2) (n = 16) (n = 10) (n = 12) (n = 14) (n = 28) (n = 24) Total (n = 52) Total % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. Total OTU ID Taxon Functional Guild (%) samples (%) samples (%) samples (%) samples (%) samples (%) samples (%) samples s__Phytophthora OTU_1 Plant Pathogen 11.1 100 17.5 100 19.7 100 7.1 100 20.9 100 7.7 100 15.2 100 plurivora s__Pythium OTU_6 Plant Pathogen 6.0 63 21.4 80 0.5 58 0.4 71 11.0 75 11.4 63 7.8 67 paroecandrum s__Phytophthora OTU_10 Plant Pathogen 1.1 69 5.8 40 11.5 75 1.8 57 7.1 67 2.7 50 6.0 62 psychrophila s__Phytophthora OTU_3 Plant Pathogen 5.5 100 2.7 90 6.1 75 4.0 86 6.0 92 2.4 81 4.7 89 quercina s__Phytophthora OTU_9 Plant Pathogen 1.5 69 2.9 50 6.7 75 0.4 100 5.2 89 0.5 50 3.5 75 citrophthora Oomycetes s__Pythium OTU_17 Plant Pathogen 5.2 56 4.1 80 1.4 75 1.3 86 3.8 78 10.9 75 3.2 73 heterothallicum s__Pythium OTU_15 Plant Pathogen 8.7 81 1.7 50 0.7 83 0.1 64 2.7 69 9.5 81 3.0 71 irregulare s__Elongis- OTU_7 porangium Plant Pathogen 3.9 100 4.6 70 0.4 83 1.3 100 2.9 92 0.8 94 2.6 90 undulatum s__Phytophthora OTU_4 Plant Pathogen 2.0 56 1.9 30 2.1 50 3.5 64 2.1 58 31.8 38 2.2 52 cinnamomi OTU_16 g__Pythium sp. Plant Pathogen 1.4 75 3.2 30 1.4 67 0.1 57 2.0 64 0.4 50 1.7 60 Total oomycete core biome frequency 46.4 65.8 50.5 20.0 63.8 78.2 49.9 s__Mortierella Soil-Plant Saprotroph OTU_3 7.5 100 4.7 100 2.8 100 1.8 100 4.7 100 3.9 100 4.4 100 elongata - Endophyte Plant Pathogen-Soil- OTU_2 g__Fusarium sp. 4.2 100 3.7 100 3.7 100 1.9 100 3.3 100 3.7 100 3.4 100 Plant Saprotroph s__Russula OTU_8 Ectomy-corrhizal 2.8 100 0.4 90 2.3 100 1.1 79 1.5 92 1.9 94 1.6 92 praetervisa s__Mortierella Soil-Plant Saprotroph OTU_11 2.8 100 1.1 80 2.0 75 6.0 100 1.4 86 1.0 100 1.3 90 elongata - Endophyte Fungi OTU_24 g__Solicoccozyma sp. Unassigned 1.0 100 0.9 100 0.8 100 1.7 100 1.1 100 1.1 100 1.1 100 OTU_9 g__Thelonectria sp. Soil-Plant Saprotroph 1.4 100 1.1 100 0.8 100 0.6 100 0.8 100 1.4 100 1.0 100 OTU_25 g__Penicillium sp. Soil-Plant Saprotroph 0.6 100 1.3 100 0.7 100 0.7 100 1.0 100 0.6 100 0.9 100 OTU_51 s__Trichoderma spp. Soil-Plant Saprotroph 0.6 100 1.1 100 0.5 92 0.7 100 0.8 100 0.7 94 0.7 98 s__Clonostachys OTU_48 Soil-Plant Saprotroph 0.6 100 1.3 100 0.2 92 0.3 100 0.6 97 0.8 100 0.7 98 rosea OTU_15 s__Trichoderma spp. Soil-Plant Saprotroph 1.1 94 0.2 100 1.2 92 0.5 93 0.9 92 0.3 100 0.7 94 Total fungal core biome frequency 22.6 15.8 15.0 15.3 16.0 15.2 15.8 Table 1. Ten most abundant oomycete and fungi OTUs found in 26 Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems, ordered by relative frequency (from core-features of each abundance matrix calculated among OTUs present in at least 50% of samples). Freq.: Number of sequences belonging at each OTU relative to the total number of sequences for the considered factor. % of samples: Percentage of samples in which the OTU was identified. significantly-lower value for EcM. Only the Animal pathogen guild showed significant die ff rences among defoli - ation levels, and no Endophyte species was identified for high defoliated ones. e S Th pearman’s correlation matrix between functional guilds (Fig.  3a) shows a significant, positive correlation of Soil-Plant Saprotrophs with Animal Pathogen abundance (ρ = 0.491, p < 0.01) and a negative one with Plant Pathogens (ρ = −0.356, p < 0.05) (Fig. 3b-1 and 2). On the other hand, EcM was negatively correlated with Soil-Plant Saprotrophs (ρ = −0.559, p < 0.01; Fig. 3b-3), Animal Pathogen (ρ = −0.492, p < 0.01) (Fig. 3b-4), Dung Saprotrophs (ρ = −0.384, p < 0.01), Opportunistic Pathogenic Species (ρ = −0.348, p < 0.05) and >1 Guild, (ρ = −0.343, p < 0.01). The Dung Saprotroph guild was positively correlated with Soil-Plant Saprotrophs (ρ = 0.332, p < 0.01) and with the species assigned to more than one guild (>1 Guild, ρ = 0.392, p > 0.01). Other significant correlations appeared between > 1 Guild OTUs, Lichenized and Opportunistic Pathogenic Species. Fungal and oomycete diversity. The Good’s Coverage indices for the selected thresholds were high for the rarefied matrices of the oomycetes (G = 0.98) and fungi (G = 0.99), indicating that over 98% of the identified OTUs were sampled in the rarefaction process and thus included in the α-diversity analysis. e S Th hannon H’ and Pielou Evenness α -diversity indices of the oomycetes exhibited similar trends according to the geographical longitude, their abundances decreasing from west to east (Fig. 4), with significant differences Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 5 www.nature.com/scientificreports www.nature.com/scientificreports/ Rank OTUs Assigned Taxon Ρ 1 OTU#1 P. plurivora −0,392** 2 OTU#6 P. paroecandrum −0,263* 3 OTU#10 P. psychrophila −0,272* 5 OTU#9 P. citrophthora −0,378** 8 OTU#7 E. undulatum −0,325* 9 OTU#4 P. cinnamomi −0,345** 12 OTU#5 P. cactorum −0,277* 13 OTU#8 P. gonapodyides −0,308* 14 OTU#2 Pythium sp. −0,395** 15 OTU#14 P. cambivora −0,379** Table 2. Significant correlations between the abundance of OTU#51 (Trichoderma spp.) and the abundance of the core biome OTUs of oomycete dataset at 26 Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. Rank: Order of the OTUs in the core features classification by abundance (1 for the most abundant, to 20 for the less abundant). ρ: Spearman’s Rho value. *Significant correlation at p < 0.05; **Significant correlation at p < 0.01. Location Defoliation sign. Functional Guild Andv Arac S.Nor Pedr Total C1 C2 t (df ) n/s n/s >1 Guild designation 3.3 ± 0.1 3.6 ± 0.2 3.8 ± 0.2 4.0 ± 0.3 3.6 ± 0.1 3.7 ± 0.1 3.6 ± 0.2 0.46 (50) b ab b a Animal Pathogen 0.9 ± 0.1 1.0 ± 0.1 0.8 ± 0.1 1.3 ± 0.1 1.0 ± 0.1** 0.9 ± 0.1 1.2 ± 0.1 −2.10 * (52) n/s n/s Arbuscular Mycorrhizal 3.4 ± 0.5 3.3 ± 0.4 5.1 ± 0.9 3.5 ± 0.9 3.8 ± 0.3 3.9 ± 0.4 3.5 ± 0.5 0.44 (50) ab c bc a n/s Dung Saprotroph 1.0 ± 0.1 0.7 ± 0.1 0.8 ± 0.1 1.1 ± 0.1 0.9 ± 0.0*** 0.9 ± 0.1 0.8 ± 0.1 0.61 (50) a a a b n/s Ectomycorrhizal 9.9 ± 0.7 9.8 ± 0.9 9.4 ± 0.4 5.5 ± 0.6 8.7 ± 0.4*** 9.0 ± 0.5 8.0 ± 0.8 1.08 (51) † b b a Endophyte 0.2 0.14 ± 0.02 0.12 ± 0.02 0.4 ± 0.1 0.22 ± 0.04** 0.22 ± 0.04 — — n/s n/s Fungal Parasite 0.24 ± 0.02 0.31 ± 0.04 0.4 ± 0.1 0.29 ± 0.04 0.30 ± 0.02 0.29 ± 0.03 0.33 ± 0.03 −0.81 (49) n/s (††) n/s Lichenized 0.3 ± 0.0 0.25 ± 0.04 0.24 ± 0.03 0.19 ± 0.04 0.25 ± 0.02 0.24 ± 0.02 0.28 ± 0.06 −0.64 (12.6) n/s n/s Oportunistic species 3.7 ± 0.2 3.4 ± 0.2 3.9 ± 0.2 4.1 ± 0.3 3.8 ± 0.1 3.8 ± 0.1 3.7 ± 0.3 0.52 (54) n/s n/s Plant Pathogen 10.0 ± 1.0 14.1 ± 5.9 21.7 ± 8.4 13.4 ± 6.1 14.5 ± 2.8 15.6 ± 3.4 11.9 ± 5.0 0.59 (50) ab ab b a n/s Soil-Plant Saprotroph 15.7 ± 0.5 15.4 ± 0.5 14.0 ± 0.5 17.1 ± 0.7 15.6 ± 0.3** 15.3 ± 0.4 16.4 ± 0.5 −1.71 (50) n/s n/s Unclassified 51.7 ± 0.7 50.5 ± 3.4 45.2 ± 4.9 51.4 ± 4.0 49.8 ± 1.7 49.0 ± 2.0 52.0 ± 3.2 −0.79 (54) Table 3. Mean ± SE relative abundance of functional guilds between study zones and defoliation levels of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. e d Th ata of functional guilds which presented significant differences (ANOVA with Bonferroni mean differences test, or Kruskal-Wallis rank sum test) were written in Italic. Andv: Andévalo. Arac: Aracena. S.Nor: Sierra Morena & Sierra Norte. Pedr: Pedroches. C1: Low average defoliation plots. C2: High average defoliation sign plots. t = t-value from the Student test with the degree of freedom and the significance level of the (df ) test.  Only one sample of Andevalo presented confident identification of endophyte organismSignificance level: *p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001; n/sp > 0.05. Values with the same letter are not significantly different (††) (p > 0.05).  Levene test p < 0.05 (variances were different between factors). between Pedr and Andv (H : 5.69, p < 0.01; H : 6.06; p < 0.01), but without differences in OTUs richness. Shannon Evenness e f Th ungal community did not show significant differences for the α -diversity indices. When these indices were compared between the mean defoliation levels for both communities, only oomycete evenness was significantly lower for the greatly-defoliated plots (χ : 4.64, DF: 1, p < 0.01); the fungal community did not show differences. Regarding β-diversity, ANOSIM showed the influence of location in the composition of the fungal commu- nity (R = 0.167; p < 0.001), whereas for the oomycete community (R = 0.078; p < 0.05) differences were only found between Pedr and Andv (R = 0.214; p < 0.01) and between Pedr and S.Nor (R = 0.125; p < 0.05). The non-linear metrics analysis (NMDS) showed a more-aggregated zonal cluster for the fungal community, in comparison with the oomycetes (Fig. 5), but also a clearer separation of the high-defoliated plots for the oomycete community. Oomycete OTU#1 (P. plurivora) had the vector which most clearly influenced the separation between defoliation classes. Discussion To the best of our knowledge, this work is the first study of soil biodiversity in holm oak dehesas ecosystems, exploring the interaction between fungal and oomycete communities and their influence in the decline symp- toms, focusing on the soil microbiome. Phytophthora spp. dominated the oomycete community, but the species usually related to Q. ilex root rot (P. cinnamomi and P. spiculum) did not appear among the most abundant, nor were they related directly to defoliation levels. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 6 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 3. Relationships between functional guilds on the 26 studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Correlation matrix between functional guilds relative abundance in each sample. *Significant correlation at p < 0.05; **significant correlation at p < 0.01. (b) Scatterplots of significant correlations between Functional Guilds: (1) Plant Pathogen/Soil-Plant Saprotroph; (2) Animal Pathogen/Soil-Plant Saprotroph; (3) Ectomycorrhizal/Soil-Plant Saprotroph; (4) Ectomycorrhizal/ Animal Pathogen. ρ: Spearman’s Rho value. p: Significance level of correlation. Figure 4. Analysis of Shannon H’ diversity, (left), OTUs richness (centre) and Pielou-e eveness vector (right) by geographical zones, for the studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems, for fungi (up) and oomycete (down). It must be considered as a limitation of our study the lack of specific data from the rhizosphere. Without this information, the discussion on the shifts of an important guild like Endophytes could not be addressed, but despite holm oak rhizosphere was not analysed, the revision of the ecological roles of the key taxa identified in this work lead us to believe that most of them should be involved in the plant-microbiota interactions, agreeing 48,49 with the results of other authors . Our results showed that the fungal community of declined dehesas includes a wide range of plant pathogens and key taxa of EcM and Trichoderma spp. which would be related to the tree health status. The metabarcoding analysis revealed that structural changes of the fungal community and the diversity level of the oomycete community were related with the mean defoliation of the holm oak dehesas ana- lyzed, agreeing with our work hypothesis. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 7 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 5. Non-metric multidimensional scaling analysis of OTUs abundance influence in defoliation levels and site location, for the 26 studied plots on Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Analysis for fungi dataset. (b) Analysis for oomycete dataset. Ellipsoid indicates normal probability of point distribution at 90% confidence level for mean defoliation classes of plots. Arrows show influence vectors for the ten most abundant OTU’s of fungi and oomycete dataset. the fungal community of dehesas has a wide range of plant pathogens and abundance of e cM. Despite the large number of fungal phytopathogens described in the studied soils, it is not possible to determine the direct influence of these species on the health status of the trees. For example, some of the fungal patho - gens identified in our samples have been found previously in dehesas as a part of the endophytic mycobiota of pasture species, such as species from the genus Colletotricum spp., Fusarium spp., Sordaria spp. or Alternaria spp. . However, the Dothideomycetes Alternaria alternata and Pyrenochaeta cava were isolated from Quercus robur roots by Kwaśna et al. , presenting significant changes in abundance according to the degree of tree defo- liation. Therefore, the increment in the abundance of some fungal pathogens, such as Dothideomycetes, in the high-defoliated plots, might be due to the pathogenicity of some species, but also to changes in the functionality of mycobiota related to the understory or grassland species. Mortierella elongata is a saprotroph often isolated in natural soils, being related their abundance with soil 52,53 improvement in nutrient quality and fertility . It plays an important role in the degradation of several com- pounds, including toxic ones , and its abundance can be considered an indicator of the resilience level of dehesas ecosystems. Of the symbiotic species, Russula was the most-abundant genus, with 76 different OTUs classified within it, including the third-most-abundant species, R. praetervisa. Ectomycorrhizal abundance and bio- diversity have been related to fine root vitality in Quercus robur forests, being considered the most-important below-ground trait related to their decline . Also, its abundance has been related to benefits for the tree health status and the establishment of Quercus spp. seedlings . In our work, we did not find a relationship between Russulaceae abundance and the individual pathogenic Phytophthora spp., but statistically-significant differences were found when R. praetervisa abundance was analysed in relation to defoliation, both results agreeing with 22,57 other studies in Q. ilex forests . In our search for plant pathogens within the core biome of the fungal community, the genus Fusarium was the most relevant. This genus includes several soil-borne plant pathogens, related to forest and agricultural plant diseases in woody, shrub and herbaceous hosts . The F. oxysporum species complex causes vascular wilts, damping-off and crown and root rot in a wide range of species, being able to produce the death of holm oak seedlings in a pathogenicity test . The role of F. eumartii in the decline of Q. robur, Q. cerris and Q. pubescens 60,61 in Italy has been recognised , and other species of this genus have been linked with diseases in other Quercus 62–64 spp. . In addition to Fusarium spp., 39 more OTUs within the fungal community core biome were assigned to plant pathogenic species, including some species known to be pathogens of oaks - such as some Peniophora spp. (among which, OTU#816 was identified in 98% of samples), Pyrenochaeta spp., Thelonectria spp., Alternaria alter- 66 51 67 nata , Sporothrix ina fl ta, Pezicula radicicola and Phomopsis spp. . The abundance of fungal pathogenic species in zones with a low H’ index for oomycetes, such as Pedr and S.Nor, coincided with a lower presence of beneficial fungal species, such as strict saprotrophs, AM or EcM. These circumstances could lead to an environment that is ideal for triggering the decline of trees. Phytophthora spp. dominate the oomycete community in declining dehesas. Phythophthora was the dominant genus among oomycetes, being the most-abundant species P. plurivora, instead of those related to 16,17,68 holm oak root rot according to the literature – P. cinnamomi and P. spiculum . Particularly, Pythium spiculum, considered a trigger of holm oak decline and mortality , was excluded from the oomycete core biome. Among the oomycete pathogens causing root rot in Quercus ilex, P. psychrophila, by its abundance, and P. quercina, by its frequency of appearance, were the most-important taxa, P. cinnamomi appearing as the ninth-most-abundant Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 8 www.nature.com/scientificreports www.nature.com/scientificreports/ oomycete, present in 52.7% of samples. However, P. cinnamomi presented a very high frequency in the samples of high-defoliation plots (C2) in which it was identified, agreeing this result with the great susceptibility that holm oak presents against this pathogen . Phytophthora plurivora is a very-common soil-borne pathogen that has been isolated from soils of many natu- ral forests in Europe. It was recorded as the most-abundant species in a survey of oak stands in Italy and in chest- nut forests of central Italy, in a study of soil Phytophthora biodiversity through metabarcoding . Its pathogenicity 71 72 73 on Quercus alba , Q. robur, Q. petraea, Q. rubra and Q. cerris is recognised. However, to our knowledge, there is no evidence of a direct pathogenic relationship with Q. ilex root rot, or of its involvement in the decline syndrome in dehesas, the work of Vettraino et al. being the only one that recorded this species (designated as P. citricola) under Q. ilex trees (subsp. ilex). The abundance of this species in our samples was significantly high in low-defoliated stands, which influenced NMDS aggregation of defoliation classes. The high recovery rate of P. plurivora and its relationship with the defoliation encourage the study of their involvement in holm oak decline. Most of the Phytophthora taxa found in our samples were also recorded in a soil analysis in the search for oomy- cete in natural Q. ilex forests of eastern Spain , where P. psychrophila and P. quercina were the most-relevant taxa, and also a qPCR survey showed P. quercina abundance in dehesas of Extremadura ; P. psychrophila and P. quercina appeared, respectively, as the third- and fourth-most-abundant oomycete species in our study. Both taxa have 74 29 been confirmed as Q. ilex root pathogens, causing root rot and tree death . However, in the work of Català et al. , no match was found for P. plurivora and P. cinnamomi. The overall specific composition of the Phythophtora spp. community in our samples was more similar to that described by Vettraino et al. , who found P. plurivora, P. quercina, P. cinnamomi and P. cactorum – in that order – as the most-frequent Phytophthora spp. isolated in stands of Q. ilex ae ff cted by oak decline. Regarding the severity of symptoms associated with the oomycete community, no significant differences in 22,57 taxa abundance were found among the defoliation categories, agreeing with the works of Corcobado et al. . Notwithstanding, defoliation and the severity of root rot should be related ; in our case, this relationship was detected for the overall structure and functionality of the oomycete community (diversity and evenness), includ- ing the influence of the fungal community (EcM abundance). Given our results, the involvement in the oak decline of all the plant pathogens present in the dehesa soils should be considered, not only that of P. spiculum and P. cinnamomi, but also other pathogenic Phytophthora species - such as P. plurivora, P. quercina and P. psychrophila - and other fungal plant pathogens such as Fusarium spp. or other phytopathogenic taxa belonging to the class Dothideomycetes. ecM and AM are key guilds related to shifting soil microbiota functionality in dehesas ecosystems. The Ectomycorrhizal (EcM) and the Arbuscular Mycorrhizal (AM) guilds, the third- and fourth-most-abundant functional guilds, respectively, had lower abundance rates overall when compared with other works in soils of dif- 45,46 ferent natural ecosystems . Quercus ilex is a species that oen p ft resents mycorrhizal association with different species of Ascomycetes and Basidiomycetes . A significant reduction in EcM and increment in plant pathogens (including Plant Pathogen and Opportunistic Species guilds) were found in Pedr, coinciding with its status as the most-high-defoliated zone, being the higher EcM abundance correlated with lower pathogen and saprobe abun- 76 57 dance, agreeing our results with other works . Corcobado et al. identified non-mycorrhizal root tips of Q. ilex as vulnerable points for pathogen invasion. Beneficial mycorrhizae, including AM and EcM, can potentially induce 77 33 higher pathogen tolerance through stimulation of the plant immune system . Therefore, either the increase of the root uptake and the effect on plant immune response might be related with lower defoliation level of trees in zones with high EcM and AM abundance. The Animal Pathogen, Soil-Plant Saprotrophs, and Dung Saprotroph guilds had their maximum values in Pedr, which showed the lowest α-diversity values for oomycetes and an increment in the abundance of Phytophthora spp., with good values for the diversity of the fungal community, but lower EcM abundance. This is another example of the greater influence of the community structure and functionality, in comparison with the specific composition. It also supports the existence of a triple relationship among the abundance of EcM and other bene- 22,57 ficial fungi, the plant pathogens and the tree health status, agreeing with other works . High diversity and changes in community structure are associated with the severity of symptoms. 45,46,78 In comparison with recent analyses of the mycobiota in soils of different ecosystems , the dehesas had a high Shannon H’ index (greater than 4.5 in all the studied zones) and a high value of evenness (exceeding 0.7). This shows a high diversity and good equilibrium of the mycobiota, but with a high proportion of plant pathogens and significant changes in the specific composition of the soil community among zones, agreeing with the results of ANOSIM. No relationship of the studied parameters with plant community changes was found, agreeing with 45,46 other authors who reported that changes in the diversity of mycobiota were decoupled from plant α-diversity . e Th values of the α -diversity indices of the oomycete community were, in general, higher than those found by Català et al. in natural forest soils of eastern Spain and were similar to (slightly higher than) those reported by Vannini et al. for chestnut forests of Italy. The lowest values of the Shannon index and evenness vector for the oomycete community occurred in Pedr, the most-defoliated zone, showing low species richness and differences in taxa abundance, agreeing with the significantly-lower values of evenness in the highly-defoliated plots. Although a direct relationship between specific oomycete taxa, or their abundance, and the severity of symptoms was not found, the β-diversity analysis showed aggregation of high-defoliated trees regarding the oomycete community. Similar results were obtained by Linaldeddu et al. when evaluating the diversity of the fungal endophytic com- munity in declining Quercus suber L. trees. The severity of defoliation, therefore, seems to be related more to the functional structure of the oomycete community – driven by the evenness - than to its specific composition, being presented the change in the specific composition of the fungal community and the reduction in competition between root pathogens as the most influential drivers related with decline symptoms. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 9 www.nature.com/scientificreports www.nature.com/scientificreports/ the abundance of Pythium spp. and Trichoderma spp. influencing the low Phytophthora spp. dominance. The abundances of Pythium and Phytophthora influenced each other, but it remains unclear how the changes were inu fl enced by location or by a natural equilibrium between the taxa. Our results show a significant decrease in oomycete α -diversity in soils of highly-defoliated stands, coinciding with an increase in the abundance of Phytophthora spp. Other authors have demonstrated that the mixture of P. cinnamomi with dif- ferent species of Pythium or other pathogenic fungi (such as Fusarium oxysporum) did not increase the severity of symptoms , supporting the hypothesis that, in soils with less diversity of oomycetes, P. cinnamomi and other pathogenic Phytophthora species have an increased ability to infect roots and ae ff ct the health status of the tree, increasing the severity of visible symptoms. Another key relationship between taxa was the significant inverse correlation between OTU#51 abundance (Trichoderma spp.) and most of the core biome taxa of oomycetes in the studied plots (Table 2). The effects of Trichoderma spp. in biocontrol strategies are well known . Trichoderma asperellum, T. hamatum, and T. virens have been isolated from holm oak roots (Q. ilex subsp. ilex) and reduced the root rot of seedlings in a controlled experiment including P. nicotianae and P. cinnamomi . Other Trichoderma spp. have also been proved to be 81–85 antagonists of diverse Phytophthora spp. . The evidence of the multiple effects of Trichoderma spp. against multiple phytopathogenic species should encourage further investigation of this relationship in the holm oak dehesas. The next steps in the search for control strategies based on the management of fungal soil biota should be the isolation of Trichoderma spp. from the studied soils to identify the fungus at species level, together with the analysis of root endosphere, which would provide useful information about the endophyte behaviour of the 86,87 present species, as it has been shown for other Trichoderma spp. . Metabarcoding is a powerful tool in the management of soilborne diseases. Biocontrol strategies are considered environmentally friendly approaches in the management of plant diseases, but the high complexity of ecological interactions between the microbiome components could lead to negative impacts on non-target organisms , and the misunderstand the mechanism of the interactions provokes the inefficiency, or even adverse effects of treatments . Thus, the study of the microbiome is an important step to develop management strategies to reduce tree decline damage through biocontrol strategies, and the knowledge about microbial community has demonstrated to have important implications in those aspects . Metabarcoding techniques could help to focus biocontrol strategies, searching for the adaptability and synergism of new products with the indigenous micro- biome, or characterizing microbiota of disease-suppressive soils and the role of key taxa in this microbiome. In this sense, further studies confirming the presence of the identified key taxa in the rhizosphere are encouraged. Additionally, it is recommended to complement metabarcoding approaches with other methods, beyond the rhiz- osphere microbiome analysis, to better understand the role of the most relevant taxa, including those taxa that are identified for the first time associated with this system. The identification of relationships between Trichoderma spp., EcM and Phytophthora spp. in soils related to holm oak decline leads to new soil management strategies in those stands ae ff cted by root rot, directly linked to the microbiome ecology. Moreover, biological control strate- gies must consider the influence of location and environmental variables in soil community , being the metabar- coding analysis an adequate technique to provide this information. Conclusion This work is a novel contribution to the study of soil biodiversity in holm oak dehesas ecosystems, analysing for the first time the fungal and oomycete communities together in dehesa soils. Our results revealed high levels of diversity in the fungal community and medium to low diversity of oomycetes, the presence of key taxa influencing holm oak health status and oomycete pathogenic abundance, and the dominance of plant pathogens and sapro- trophs in the analysed soils. The microbiome presented significant influence in the defoliation status, confirming our main hypothesis, and suggesting that both the structure of the fungal community and the level of diversity of the oomycete community were related with the tree health status. On the other hand, the high degree of soil microbiome diversity in declined dehesas revealed on this work encourages to a deeper research into the relationship between different taxa, considering a tripartite relationship (plant-pathogens-other taxa) , but also in the study of the influence of environmental factors in the soil com- munity, beyond the simple picture of a single plant-pathogen interaction. The findings about the presence of key taxa related with defoliation symptoms or pathogen abundance, must drive us to consider biocontrol strategies in the management of the holm oak root rot, through methods and applications that respect the local microbiota. Several biocontrol strategies could be explored deriving from our results, such as the application of Trichoderma spp. to control the inoculum of pathogenic oomycetes, as well as the application of mycorrhizae as a strategy to improve the health status of trees, but further research is needed to clarify how both key taxa interacts between them and with holm oak. Material and Methods study area. This study was focused on holm oak (Quercus ilex subsp. ballota) dehesas of north-west Andalusia (Fig. 6). The studied area was divided into four zones - Andévalo (Andv), Sierra de Aracena (Arac), Sierra Norte de Sevilla and S. Morena (S.Nor) and Valle de los Pedroches (Pedr) - located in the provinces of Huelva, Sevilla and Córdoba, representing the core area of the holm oak dehesa ecosystems of Andalusia (Fig. 6). This division was made based on differences in environmental factors and silvicultural characteristics (Table  4). A total of 26 plots were selected among the permanent environmental monitoring points of the Andalusian Network for Damage Monitoring in Forest Ecosystems (Red de Información Ambiental, Consejería de Medio Ambiente y Ordenación del Territorio - Junta de Andalucía, ICP Forest, EC Network level II - www.juntadeanda- lucia.es). The plots were selected based on the following criteria; (i) >85% of the trees were Q. ilex, (ii) a positive record of P. cinnamomi and (iii) the tree mortality rate. To calculate the mean plot defoliation, all the trees of the Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 10 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 6. Map of sampled plots of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. plots (24 trees × 26 plots) were assessed for defoliation following the methodology of Eichhorn et al. . The mean defoliation for each plot (%), calculated using the individual defoliation of the 24 trees per plot, ranged between 13.3% and 56.9% resulting in 14 plots classified as C1 (mean defoliation < 25%), and 12 plots classified as C2 (25% ≤ mean defoliation < 65%) (Supplementary Information, Table S2). The mean defoliation showed signif- icant differences among the study zones (χ = 11.1; DF = 3; p < 0.05), the highest value corresponding to Pedr. soil sampling. In each plot, soil samples were collected from the rhizosphere under the crown of two trees with defoliation degree similar to the mean defoliation of plot (28 samples under C1 trees and 24 samples under C2 trees, n = 52), searching for the representability of the mean plot conditions (soil, exposure and slope). For each tree selected, sub-samples were taken under the crown for two orientations (North and South) and two distances from the trunk (0.5 m and 1–2 m), making four sub-samples per tree. Then, the sub-samples of each tree were bulked and stored at 4 °C until processing. Samples were collected from November to December 2015. The soil samples were processed within the next 24 hours following their collection; soil was dried at room temperature for 24 h and well-mixed, litter and roots were separated, and 20 g of the fine fraction was collected aer p ft assing the soil through a 1-mm sieve. DNA extraction, ITS1 amplification, library generation, and sequencing. Total DNA was extracted from soil using the MoBIO Power Soil DNA Isolation kit (MoBIO Laboratories, Inc., CA, USA). Five replicates of fresh soil (0.33 g each) were taken from each soil sample. The aliquots were mixed together, and total DNA was purified and concentrated using a DNA purification kit (Wizard SV Genomic DNA Purification System; Promega Corporation, WI, USA). The DNA concentration was quantified, with a Qubit dsDNA BR Assay kit (Thermo Fisher Scientific Inc., MA, USA), prior to amplification. Additionally, two mock communities were created by pooling DNA from a pure culture of 13 different species of oomycetes and another with seven species of fungi (Supplementary Information, Table S3); one community of each taxa had the same concentration of each species (Oo= and F=) and the other had different concentrations (Oo≠ and F≠). These were used as controls and were processed in the same way as the samples, following the methodology of C. Morales Rodriguez (Unpublished). e in Th ternal transcribed spacer 1 (ITS1) was amplified in a multiplexing PCR using a set of barcoded primers: 92 93 ITS1F and ITS2 for fungi, and ITS6 and ITS7 for oomycetes . e Th primers were tagged with different barcodes to distinguish different samples. Moreover, 3 replicates of each control community were amplified with different barcoded primers of the ITS1F/ITS2 and ITS6/ITS7 primer pairs (Supplementary Information, Tables S4, S5). Negative controls and cross-reactions (Oo with ITS1F/2 and F with ITS 6/7) did not originate any PCR product. e PCR r Th eaction mixture consisted of 25 μl of Maxima Hot Start PCR Master Mix (Thermo Fisher Scientific, −1 −1 USA), 2 μM of each primer, 4 µ l of 2 mg ml BSA and 3 μl of DNA template diluted to 50 ng µ l , in a total volume of 50 μl. The thermal cycle was an initial 5-min step at 95 °C, followed by cycles (30 for fungal ITS and 35 for oomycete ITS) of 40 s of denaturation at 95 °C, 2 min of annealing (58 °C for fungal ITS amplification and 60 °C for oomycete ITS) and 1 min of extension at 72 °C, with a final extension step of 7 min at 72 °C. For each sample and the controls, three replicates of PCR for fungal ITS and five for oomycetes ITS were carried out and the products were pooled. Amplicons were purified using the MagJet NGS Cleanup and Size Selection kit (Thermo Fisher Scientific Inc., MA, USA) and the final concentration of purified PCR product was quantified using a Qubit dsDNA HS Assay kit (Thermo Fisher Scientific Inc., MA, USA). The resulting amplicons were pooled in equal amounts to give a total −1 volume of 40 µ L, with a concentration of 200 ng µ L , and were sent to Eurofins Genomics AT GmbH (Vienna, Austria), to be sequenced in an Illumina MiSeq 2 × 300 bp platform. Both fungal and oomycete libraries were Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 11 www.nature.com/scientificreports www.nature.com/scientificreports/ Andevalo Aracena Sierra Morena & S. Norte Pedroches Mean Precipitation (mm) 675.7 829.2 712.1 572 Mean Annual temperature (°C) 17.4 16.2 16.1 15.8 Mean Evapotranspiration (mm) 1445 1396 1458 1443 Max elevation (masl) 293 1046 936 943 Mean elevation (masl) 191 434 429 570 Average cover* Low Medium High Low Dominant tree age** Young to Mature Mature Young Old to very old Normal tree cover with High tree cover with Class of dehesa Cultivated dehesa Pastures shrubs and pasture shrubs and pasture Table 4. Mean values of environmental and silvicultural parameters considered for the selection of the study zones (m.a.s.l.: meters above sea level). Average values for the period 1969–2000, taken from the Andalusia Environmental Information Network - REDIAM; Consejería de Medio Ambiente, Junta de Andalucía (http:// www.juntadeandalucia.es/medioambiente/site/rediam/portada/). The average tree cover, the dominant tree age and the prevalent class of dehesa were determined as the trend of all the visited plots for each zone, following the 104 −1 −1 criteria of Costa Pérez et al. . *Low = N < 50 tree ha ; Medium = 50 < N < 80 tree ha ; High = N > 80 tree −1 ha . **Young trees = 50–90 years/20–40 cm Ø; Mature trees = 110–210 years/50–80 cm Ø; Old trees = 210– 250 years/90–110 cm Ø; Very old = > 250 years/Ø > 110 cm. included on the same Illumina run. The reads are available in the EMBL Nucleotide Sequence Database (http:// www.ebi.ac.uk/embl) under the project name “NGS_Dehesa” and accession number PRJEB28448. s equence processing, ot U selection, and abundance matrix construction. A bioinformatic pipe- 94 95 line was adapted from Bálint et al. and Sapkota & Nicolaisen . The primers were debarcoded and classified by Eurofins Genomic AT GmbH (Vienna, Austria). The sequences were joined (22 774 511 sequences) and reads with mismatches or ambiguous bases were unassigned (1 246 294 sequences). Finally, only paired sequences in which the 5′ barcode and forward primer and 3′ barcode and reverse primer were found, were used in the bioin- formatic analysis (Supplementary Information, Table S6). Singletons were excluded from the analysis. e f Th orward and reverse sequences of paired reads were trimmed, the primers were eliminated, and sequences were overlapped and paired using CLC Genomics Workbench 3.6.5 (https://www.qiagenbioinformatics.com/) (Supplementary Information, Table S7). The ITS1 region was extracted with the Fungal ITS Extractor version 2 , sequences were clustered in operational taxonomic units (OTUs) and the chimeric sequences were filtered - de-novo chimera filtering for oomycete files - using USEARCH 8.1 with a similarity level established at 98%. The consensus sequences of the OTUs were identified with the GenBank nucleotide database (nr/nt) of NCBI, using the Megablast algorithm, and the results were parsed, and the sequences retained in MEGAN 6.10.8 to assign taxonomic identifications to each OTU (Supplementary Information, Table S8). In some cases, ITS1 was una- ble to die ff rentiate between species belonging to genera which presented a high conserved ITS1 region, such as Trichoderma spp., but aer a t ft horough supervision of the taxonomy, the metabarcoding analysis showed enough accuracy to study the diversity and composition of the fungal and oomycete communities, clearly differentiating between genera in all cases. Ae ft rward, abundance tables of OTUs for oomycete and fungal files were constructed using USEARCH 8.1 with a level of identity of 98%. The amplified ITS sequences from all the species included in the internal controls were correctly identified. The mean frequency of debarcoded paired sequences did not die ff r between oomycetes and fungi ( χ = 0.727; p = 0.394), with strong correspondence between the abundance of OTUs and the estimated amount of DNA included in each control. Prior to downstream analysis of the abundance matrix, control samples and OTUs with less than 10 reads were eliminated from both the fungal and oomycete feature tables. Fungal and oomycete community analysis. Aer t ft he abundance matrix filtering, functional guilds of the identified taxa were characterized using FunGuild database (accessed on March, 2018) . Only OTUs classi- fied with a confidence level of “Probable” or “Highly Probable” were used. e Th functional guilds of the remaining OTUs, classified as “Unassigned” or “Possible”, were revised based on the literature; those that were confidently classified were added to the analysis (Supplementary Information, Table S9). The rest of the unassigned OTUs and those assigned with low confidence were placed in the “Unclassified” category. Taxa classified as soil sapro- trophs, plant saprotrophs, litter saprotrophs or wood saprotrophs were grouped in a unique Soil-Plant Saprotroph category. Those taxa that were confidently classified in more than one guild were assigned to the class > 1 Guild, except in the case of taxa that were confidently assigned to four or more guilds, including pathogen and sapro- troph guilds; these were classified in the Opportunistic Pathogenic Species class. To analyze the diversity of oomycetes and fungi, abundance matrices were processed using QIIME2 version 2017.12.1 . Rarefaction curves, Good’s coverage, Shannon H’, OTUs abundance and Pielou evenness vector and the Bray-Curtis distance matrix were calculated using the plugin Diversity 2017.12.0 . e Th most abundant and frequent OTUs, presented in at least 50% of the samples, were selected and classified as core features using the Feature-table plugin . Aer ft wards, the 10 most-abundant OTUs for fungi and oomycetes were selected among them and classified as the core biome (Supplementary Information Table S10). The relative abundance of the top 10 core features with respect to the number of paired sequences was calculated separately for each dataset Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 12 www.nature.com/scientificreports www.nature.com/scientificreports/ (oomycetes and fungi). Rarefaction curves (Supplementary Information, Fig. S1) were evaluated and Good’s Coverage was calculated to estimate the adequate depth of sampling in the abundance matrices, the abundance matrices being rarefied with an even sampling depth for each dataset, with 10 iterations to avoid bias due to uneven sequencing depth. The Shannon (H’) diversity index, OTUs diversity and Evenness vector indices were calculated to study the diversity of the fungal and oomycete communities in each sample and the Bray-Curtis dis- similarity matrix was calculated; the β-diversity was analyzed to find the differences in community composition among the study zones. statistical methods. All the statistical analyses were performed in the environment R Studio V 1.0.143 (RStudio Inc. 2015. Boston, MA, USA) running under R 3.4.0 (R Foundation for Statistical Computing, 2014. Vienna, Austria). Prior to the analysis of differences of mean defoliation between zones, and for taxon abundance among defo- liation level or study zones, the normality of variable residues was checked through the normal Q-Q graph and Shapiro-Wilk test. If the data did not satisfy the assumptions of homoscedasticity and normality, the logarithmic transformation was used and then, normality checked again. For normally-distributed variables, the differences were compared using ANOVA followed by a post hoc Bonferroni mean differences test. The Kruskal-Wallis rank sum test with the Bonferroni correction was used instead ANOVA when variables did not show normal distribu- tion. Student-t Test was used for comparisons between defoliation levels. A correlation matrix was calculated to analyse the relationships between those functional guilds present in five or more samples of all zones, thus excluding the Lichen Parasite and Endophyte guilds. Unclassified features were excluded from the correlation matrix to avoid artifacts in the analysis. The correlation analysis was used also to evaluate the relationships between the key taxa resulting from the core biome analysis. The Pearson or Spearman correlation test was used, depending on whether the two variables in each pair were distributed normally or not. Null hypotheses were rejected in all cases when p ≤ 0.05. Alpha-diversity indices were subjected to the Kruskall-Wallis rank sum test with the Bonferroni correction, to compare the diversity levels of fungal and oomycete communities in the distinct study zones and according to the mean defoliation class of the plots. The Bray-Curtis dissimilarity matrix was used in non-metric mul- tidimensional scaling analysis (NMDS, vegan package , together with Analysis of Similarity (ANOSIM, 999 permutations), to determine the influence of the study location on the composition of the fungal and oomycete communities. Other packages used for data analysis and representation were dunn.test, nortest, ggplot2, ggpurb, and devtools (https://CRAN.R-project.org/). References 1. Prospero, S. & Cleary, M. Effects of Host Variability on the Spread of Invasive Forest. Diseases. Forests 8, 80 (2017). 2. Liebhold, A. M. et al. Biological invasions in forest ecosystems. Biol. Invasions 19, 3437–3458 (2017). 3. Sena, K., Crocker, E., Vincelli, P. & Barton, C. Phytophthora cinnamomi as a driver of forest change: Implications for conservation and management. For. Ecol. Manag. 409, 799–807 (2018). 4. Rizzo, D. M. & Garbelotto, M. Sudden oak death: endangering California and Oregon forest ecosystems. Front. Ecol. Environ. 1, 197–204 (2003). 5. Burgess, T. I. et al. Current and projected global distribution of Phytophthora cinnamomi, one of the world’s worst plant pathogens. Glob. Change Biol. 23, 1661–1674 (2017). 6. Vettraino, A. M. et al. Occurrence and distribution of Phytophthora species in European chestnut stands, and their association with Ink Disease and crown decline. Eur. J. Plant Pathol. 111, 169 (2005). 7. Thomas, F. M., Blank, R. & Hartmann, G. Abiotic and biotic factors and their interactions as causes of oak decline in Central. Europe. For. Pathol. 32, 277–307 (2002). 8. Jung, T., Blaschke, H. & Osswald, W. Involvement of soilborne Phytophthora species in Central European oak decline and the effect of site factors on the disease. Plant Pathol. 49, 706–718 (2000). 9. Vettraino, A. M. et al. Occurrence of Phytophthora species in oak stands in Italy and their association with declining oak trees. For. Pathol. 32, 19–28 (2002). 10. Brasier, C. M. Phytophthora cinnamomi and oak decline in southern Europe. Environmental constraints including climate change. Ann. Sci. For. 53, 347–358 (1996). 11. Navarro-Cerrillo, R. M. et al. Los procesos de decaimiento de encinas y alcornoques. (Dirección General de Gestión del Medio Natural. Consjería de Medio Ambiente. Junta de Andalucía., 2004). 12. Díaz Esteban, M. & Pulido Díaz, F. 6310: Dehesas perennifolias de Quercus spp. in Bases ecológicas preliminares para la conservación de los tipos de hábitat de interés comunitario en España (Ministerio de Medio Ambiente y Medio Rural y Marino, 2009). 13. Marañón, T. Agro-Sylvo-Pastoral Systems in the Iberian Peninsula: Dehesas and Montados. Rangelands 10, 255–258 (1988). 14. Pulido, F. & Picardo, A. Libro Verde de la Dehesa. Documento para el debate hacia una Estrategia Ibérica de gestión. (Consejería de Medio Ambiente, Junta de Castilla y León, 2010). 15. Horta, M., Caetano, P., Medeira, C., Maia, I. & Cravador, A. Involvement of the β-cinnamomin elicitin in infection and colonisation of cork oak roots by Phytophthora cinnamomi. Eur. J. Plant Pathol. 127, 427–436 (2010). 16. Brasier, C. M., Robredo, F. & Ferraz, J. F. P. Evidence for Phytophthora cinnamomi involvement in Iberian oak decline. Plant Pathol. 42, 140–145 (1993). 17. Oßwald, W. et al. Strategies of attack and defence in woody plant–Phytophthora interactions. For. Pathol. 44, 169–190 (2014). 18. Gea-Izquierdo, G., Fernández-de-Uña, L. & Cañellas, I. Growth projections reveal local vulnerability of Mediterranean oaks with rising temperatures. For. Ecol. Manag. 305, 282–293 (2013). 19. Corcobado, T., Cubera, E., Juárez, E., Moreno, G. & Solla, A. Drought events determine performance of Quercus ilex seedlings and increase their susceptibility to Phytophthora cinnamomi. Agric. For. Meteorol. 192–193, 1–8 (2014). 20. Ruiz Gómez, F., Pérez-de-Luque, A., Sánchez-Cuesta, R., Quero, J. & Navarro Cerrillo, R. Differences in the response to acute drought and Phytophthora cinnamomi Rands Infection in Quercus ilex L. seedlings. Forests 9, 634 (2018). 21. Ruiz Gómez, F. J., Navarro-Cerrillo, R. M., Sánchez-Cuesta, R. & Pérez-de-Luque, A. Histopathology of infection and colonization of Quercus ilex fine roots by Phytophthora cinnamomi . Plant Pathol. 64, 605–616 (2015). 22. Corcobado, T., Moreno, G., Azul, A. M. & Solla, A. Seasonal variations of ectomycorrhizal communities in declining Quercus ilex forests: interactions with topography, tree health status and Phytophthora cinnamomi infections. Forestry 88, 257–266 (2015). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 13 www.nature.com/scientificreports www.nature.com/scientificreports/ 23. Duque-Lazo, J., Navarro-Cerrillo, R. M., van Gils, H. & Groen, T. A. Forecasting oak decline caused by Phytophthora cinnamomi in Andalusia: Identification of priority areas for intervention. For. Ecol. Manag. 417, 122–136 (2018). 24. Serrano, M. S., Fernández-Rebollo, P., De Vita, P. & Sánchez, M. E. Calcium mineral nutrition increases the tolerance of Quercus ilex to Phytophthora root disease ae ff cting oak rangeland ecosystems in Spain. Agrofor. Syst. 87, 173–179 (2013). 25. Navarro Cerrillo, R. M., Terán Bocero, A. I. & Sánchez, M. E. Acción preventiva y curativa del fosfonato en el control de Phytophthora cinnamomi Rands en encina y alcornoque. 32, 685–694 (2006). 26. Ríos, P. et al. Brassica -based seedmeal biofumigation to control Phytophthora cinnamomi in the Spanish “dehesa” oak trees. Phytopathol. Mediterr. 56, 392–399 (2018). 27. Gómez-Aparicio, L. et al. Spatial patterns of soil pathogens in declining Mediterranean forests: implications for tree species regeneration. New Phytol. 194, 1014–1024 (2012). 28. Corcobado, T., Cubera, E., Pérez-Sierra, A., Jung, T. & Solla, A. First report of Phytophthora gonapodyides involved in the decline of Quercus ilex in xeric conditions in Spain. New Dis. Rep. 22, 33 (2010). 29. Català, S., Berbegal, M., Pérez-Sierra, A. & Abad-Campos, P. Metabarcoding and development of new real-time specific assays reveal Phytophthora species diversity in holm oak forests in eastern Spain. Plant Pathol. 66, 115–123 (2017). 30. Mora-Sala, B., Berbegal, M. & Abad-Campos, P. The use of qPCR reveals a high frequency of Phytophthora quercina in two Spanish holm oak areas. Forests 9, 697 (2018). 31. Gdanetz, K. & Trail, F. The wheat microbiome under four management strategies, and potential for endophytes in disease protection. Phytobiomes 1, 158–168 (2017). 32. Deveau, A. et al. Bacterial–fungal interactions: ecology, mechanisms and challenges. FEMS Microbiol. Rev. (2018). https://doi. org/10.1093/femsre/fuy008. 33. Pérez-de-Luque, A. et al. The interactive effects of arbuscular mycorrhiza and plant growth-promoting rhizobacteria synergistically enhance host plant defences against pathogens. Sci. Rep. 7 (2017). 34. Aleandri, M. P. et al. Use of nursery potting mixes amended with local Trichoderma strains with multiple complementary mechanisms to control soil-borne diseases. Crop Prot. 67, 269–278 (2015). 35. Blom, J. M., Vannini, A., Vettraino, A. M., Hale, M. D. & Godbold, D. L. Ectomycorrhizal community structure in a healthy and a Phytophthora-infected chestnut (Castanea sativa Mill.) stand in central Italy. Mycorrhiza 20, 25–38 (2009). 36. Lancellotti, E. & Franceschini, A. Studies on the ectomycorrhizal community in a declining Quercus suber L. stand. Mycorrhiza 23, 533–542 (2013). 37. Malcolm, G. M., Kuldau, G. A., Gugino, B. K. & Jiménez-Gasco, M. D. M. Hidden host plant associations of soilborne fungal pathogens: an ecological perspective. Phytopathology 103, 538–544 (2013). 38. Lee, S.-J. & Rose, J. K. C. Mediation of the transition from biotrophy to necrotrophy in hemibiotrophic plant pathogens by secreted effector proteins. Plant Signal. Behav. 5, 769–772 (2010). 39. Fawke, S., Doumane, M. & Schornack, S. Oomycete interactions with plants: infection strategies and resistance principles. Microbiol. Mol. Biol. Rev. 79, 263–280 (2015). 40. Vannini, A., Bruni, N., Tomassini, A., Franceschini, S. & Vettraino, A. M. Pyrosequencing of environmental soil samples reveals biodiversity of the Phytophthora resident community in chestnut forests. FEMS Microbiol. Ecol. 85, 433–442 (2013). 41. Sapkota, R. & Nicolaisen, M. Cropping history shapes fungal, oomycete and nematode communities in arable soils and affects cavity spot in carrot. Agric. Ecosyst. Environ. 257, 120–131 (2018). 42. Martin, F. N., Abad, Z. G., Balci, Y. & Ivors, K. Identic fi ation and detection of Phytophthora: reviewing our progress, identifying our needs. Plant Dis. 96, 1080–1103 (2012). 43. Hawksworth, D. L. Global species numbers of fungi: are tropical studies and molecular approaches contributing to a more robust estimate? Biodivers. Conserv. 21, 2425–2433 (2012). 44. Tedersoo, L. et al. Global diversity and geography of soil fungi. Science 346, 1256688–1256688 (2014). 45. Veach, A. M., Stokes, C. E., Knoepp, J., Jumpponen, A. & Baird, R. Fungal communities and functional guilds shift along an elevational gradient in the Southern Appalachian Mountains. Microb. Ecol. 1–13, https://doi.org/10.1007/s00248-017-1116-6 (2017). 46. Chen, W. et al. Plant diversity is coupled with beta not alpha diversity of soil fungal communities following N enrichment in a semi-arid grassland. Soil Biol. Biochem. 116, 388–398 (2018). 47. Abdelfattah, A., Malacrinò, A., Wisniewski, M., Cacciola, S. O. & Schena, L. Metabarcoding: A powerful tool to investigate microbial communities and shape future plant protection strategies. Biol. Control 120, 1–10 (2018). 48. Coleman-Derr, D. et al. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 209, 798–811 (2016). 49. Berendsen, R. L., Pieterse, C. M. J. & Bakker, P. A. H. M. The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486 (2012). 50. Lledó, S., Santamaría, O., Rodrigo, S. & Poblaciones, M. J. Endophytic mycobiota associated with Trifolium subterraneum growing under semiarid conditions. Ann. Appl. Biol. 168, 243–254 (2016). 51. Kwaśna, H., Szewczyk, W. & Behnke-Borowczyk, J. Fungal root endophytes of Quercus robur subjected to flooding. For. Pathol. 46, 35–46 (2016). 52. Nordberg, H. et al. The genome portal of the Department of Energy Joint Genome Institute: 2014 updates. Nucleic Acids Res. 42, D26–31 (2014). 53. Tamayo-Vélez, Á. & Osorio, N. W. Soil fertility improvement by litter decomposition and inoculation with the fungus Mortierella sp. in avocado plantations of Colombia. Commun. Soil Sci. Plant Anal. 49, 139–147 (2018). 54. Li, F. et al. Application of environmental DNA metabarcoding for predicting anthropogenic pollution in rivers. Environ. Sci. Technol. 52, 11708–11719 (2018). 55. Mosca, E., Montecchio, L., Barion, G., Dal Cortivo, C. & Vamerali, T. Combined effects of thinning and decline on fine root dynamics in a Quercus robur L. forest adjoining the Italian Pre-Alps. Ann. Bot. 119, 1235–1246 (2017). 56. Dickie, I. A., Guza, R. C., Krazewski, S. E. & Reich, P. B. Shared ectomycorrhizal fungi between a herbaceous perennial (Helianthemum bicknellii) and oak (Quercus) seedlings. New Phytol. 164, 375–382 (2004). 57. Corcobado, T., Vivas, M., Moreno, G. & Solla, A. Ectomycorrhizal symbiosis in declining and non-declining Quercus ilex trees infected with or free of Phytophthora cinnamomi. For. Ecol. Manag. 324, 72–80 (2014). 58. Aoki, T., O’Donnell, K. & Geiser, D. M. Systematics of key phytopathogenic Fusarium species: Current status and future challenges. J. Gen. Plant Pathol. 80, 189–201 (2014). 59. Gallego, F. J., de Algaba, A. P. & Fernandez-Escobar, R. Etiology of oak decline in Spain. Eur. J. For. Pathol. 29, 17–27 (1999). 60. Ragazzi, A., Mugnai, L., Moricca, S., Vagniluca, S. & Dellavalle, I. Requirements and biological aspects of Fusarium eumartii and its possible role in oak decline in North-eastern Italian oak forests. For. Pathol. 23, 171–177 (1993). 61. Ragazzi, A., Moricca, S., Dellavalle, I. & Turco, E. Fusarium eumartii growth in resistant and susceptible oak species. J. Phytopathol. 149, 503–507 (2001). 62. Lynch, S. C., Zambino, P. J., Mayorquin, J. S., Wang, D. H. & Eskalen, A. Identification of new fungal pathogens of coast live oak in California. Plant Dis. 97, 1025–1036 (2013). 63. Wit, M., Sierota, Z., Oszako, T., Mirzwa-Mroz, E. & Wakulinski, W. Fusarium spp. on the above-ground organs of dying oaks - a new threat? Sylwan 159, 403–410 (2015). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 14 www.nature.com/scientificreports www.nature.com/scientificreports/ 64. Moricca, S. et al. Biotic factors ae ff cting ecosystem services in urban and peri-urban forests in Italy: The role of introduced and impending pathogens and pests. Forests 9, 65 (2018). 65. Pecoraro, L. et al. Macrofungi in Mediterranean maquis along seashore and altitudinal transects. Plant Biosyst. - Int. J. Deal. Asp. Plant Biol. 148, 367–376 (2014). 66. Kwaśna, H. & Szewczyk, W. Effects of fungi isolated from Quercus robur roots on growth of oak seedlings. Dendrobiology 75, 99–112 (2016). 67. Omari, B. E., Fleck, I., Aranda, X. & Moret, A. & MartíNadal. Effect of fungal infection on leaf gas-exchange and chlorophyll uo fl rescence in Quercus ilex. Ann. For. Sci. 58, 165–174 (2001). 68. Sánchez, M., Caetano, P., Ferraz, J. & Trapero, A. Phytophthora disease of Quercus ilex in south-western Spain. For. Pathol. 32, 5–18 (2002). 69. Serrano, M. D. et al. Phytophthora cinnamomi and Pythium spiculum as main agents of Quercus decline in southern Spain and Portugal. In IOBC/WPRS Bulletin 76, 97–100 (2012). 70. Moralejo, E., García-Muñoz, J. A. & Descals, E. Susceptibility of Iberian trees to Phytophthora ramorum and P. cinnamomi. Plant Pathol. 58, 271–283 (2009). 71. McConnell, M. E. & Balci, Y. Phytophthora cinnamomi as a contributor to white oak decline in mid-atlantic United States forests. Plant Dis. 98, 319–327 (2013). 72. Jung, T. & Burgess, T. I. Re-evaluation of Phytophthora citricola isolates from multiple woody hosts in Europe and North America reveals a new species, Phytophthora plurivora sp. nov. Persoonia 22, 95–110 (2009). 73. Jung, T., Blaschke, H. & Neumann, P. Isolation, identification and pathogenicity of Phytophthora species from declining oak stands. Eur. J. For. Pathol. 26, 253–272 (1996). 74. Pérez-Sierra, A. et al. Previously unrecorded low-temperature Phytophthora species associated with Quercus decline in a Mediterranean forest in eastern Spain. For. Pathol. 43, 331–339 (2013). 75. Coll, L., Camarero, J. J. & Aragón, J. M. D. Fine root seasonal dynamics, plasticity, and mycorrhization in 2 coexisting Mediterranean oaks with contrasting aboveground phenology. Écoscience 19, 238–245 (2012). 76. Richard, F., Moreau, P.-A., Selosse, M.-A. & Gardes, M. Diversity and fruiting patterns of ectomycorrhizal and saprobic fungi in an old-growth Mediterranean forest dominated by Quercus ilex L. Can. J. Bot. 82, 1711–1729 (2004). 77. Dawkins, K. & Esiobu, N. Arbuscular and ectomycorrhizal fungi associated with the invasive Brazilian pepper tree (Schinus terebinthifolius) and two native plants in south Florida. Front. Microbiol. 8 (2017). 78. Jimu, L., Kemler, M., Mujuru, L. & Mwenje, E. Illumina DNA metabarcoding of Eucalyptus plantation soil reveals the presence of mycorrhizal and pathogenic fungi. For. Int. J. For. Res. 91, 238–245 (2018). 79. Linaldeddu, B. T., Sirca, C., Spano, D. & Franceschini, A. Variation of endophytic cork oak-associated fungal communities in relation to plant health and water stress. For. Pathol. 41, 193–201 (2011). 80. Vinale, F. et al. Trichoderma–plant–pathogen interactions. Soil Biol. Biochem. 40, 1–10 (2008). 81. Widmer, T. L. Screening Trichoderma species for biological control activity against Phytophthora ramorum in soil. Biol. Control 79, 43–48 (2014). 82. García-Núñez, H. G. et al. Caracterización morfológica y molecular de cepas nativas de Trichoderma y su potencial de biocontrol sobre Phytophthora infestans. Rev. Mex. Fitopatol. 35, 58–79 (2017). 83. Chemeltorit, P. P., Mutaqin, K. H. & Widodo, W. Combining Trichoderma hamatum THSW13 and Pseudomonas aeruginosa BJ10–86: a synergistic chili pepper seed treatment for Phytophthora capsici infested soil. Eur. J. Plant Pathol. 147, 157–166 (2017). 84. Promwee, A., Yenjit, P., Issarakraisila, M., Intana, W. & Chamswarng, C. Efficacy of indigenous Trichoderma harzianum in controlling Phytophthora leaf fall (Phytophthora palmivora) in Thai rubber trees. J. Plant Dis. Prot. 124, 41–50 (2017). 85. Widmer, T. L. et al. Remediation of Phytophthora ramorum-infested soil with Trichoderma asperellum isolate 04–22 under ornamental nursery conditions. Biol. Control 118, 67–73 (2018). 86. Bae, H. et al. The beneficial endophyte Trichoderma hamatum isolate DIS 219b promotes growth and delays the onset of the drought response in Theobroma cacao . J. Exp. Bot. 60, 3279–3295 (2009). 87. Carrero‐Carrón, I. et al. Interactions between Trichoderma harzianum and defoliating Verticillium dahliae in resistant and susceptible wild olive clones. Plant Pathol. 67, 1758–1767 (2018). 88. Menkis, A. et al. Occurrence and impact of the root-rot biocontrol agent Phlebiopsis gigantea on soil fungal communities in Picea abies forests of northern Europe. FEMS Microbiol. Ecol. 81, 438–445 (2012). 89. Sylla, J. et al. Leaf microbiota of strawberries as ae ff cted by biological control agents. Phytopathology 103, 1001–1011 (2013). 90. Pinto, C. et al. Unravelling the diversity of grapevine microbiome. PLOS ONE 9, e85622 (2014). 91. Eichhorn, J. et al. Part IV: Visual assessment of crown condition and damaging agents. In UNECE ICP Forests Programme Co- ordinating Centre. Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. 54 (Thünen Institute of Forest Ecosystems, 2017). 92. Gao, Z., Perez-Perez, G. I., Chen, Y. & Blaser, M. J. Quantitation of major human cutaneous bacterial and fungal populations. J. Clin. Microbiol. 48, 3575–3581 (2010). 93. Cooke, D. E. L., Drenth, A., Duncan, J. M., Wagels, G. & Brasier, C. M. A molecular phylogeny of Phytophthora and related oomycetes. Fungal Genet. Biol. 30, 17–32 (2000). 94. Bálint, M., Schmidt, P.-A., Sharma, R., Thines, M. & Schmitt, I. An Illumina metabarcoding pipeline for fungi. Ecol. Evol. 4, 2642–2653 (2014). 95. Sapkota, R. & Nicolaisen, M. An improved high throughput sequencing method for studying oomycete communities. J. Microbiol. Methods 110, 33–39 (2015). 96. Nilsson, R. H. et al. An open source sowa ft re package for automated extraction of ITS1 and ITS2 from fungal ITS sequences for use in high-throughput community assays and molecular ecology. Fungal Ecol. 3, 284–287 (2010). 97. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010). 98. Huson, D. H. et al. MEGAN Community Edition - Interactive exploration and analysis of large-scale microbiome sequencing data. PLOS Comput. Biol. 12, e1004957 (2016). 99. Nguyen, N. H. et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248 (2016). 100. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010). 101. Chen, J. et al. Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinforma. Oxf. Engl. 28, 2106–2113 (2012). 102. McDonald, D. et al. e Th Biological Observation Matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. GigaScience 1, 7 (2012). 103. Oksanen, J. et al. Vegan: Community Ecology Package. (2018). 104. Costa Pérez, J. C., Martín Vicente, Á., Fernández Alés, R. & Estirado Oliet, M. Dehesas de Andalucía. Caracterización ambiental (2006). (Junta de Andalucía. Consejería de Medio Ambiente, 2006). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 15 www.nature.com/scientificreports www.nature.com/scientificreports/ Acknowledgements This work was conducted with support from the projects ESPECTRAMED (CGL2017-86161-R) and ENCINOMICA (BIO2015-64737-R), funded by the Spanish Ministry of Economy, Industry and Competitivity. e fir Th st author was a recipient of an FPU fellowship, and the work was carried out during a research stay in the TUM School of Life Sciences Weihensthephan (Freising, Germany), funded by the Spanish Ministry of Education (FPU fellowship, FPU13/00231). The authors thank the “Consejería de Medio Ambiente, Junta de Andalucía” for their support in sample collection, Rafael Sánchez de la Cuesta and Roberto Jesús Cabrera Puerto by their support in the sample processing and the staff of the Dendrochronology, Silviculture, and Climate Change Laboratory at Cordoba University, for their assistance during this research. We thank Dr. David Walker for revisions of the written English in different versions of this manuscript. Also, authors thank the support of the “Campus de Excelencia en Agroalimentación - CeiA3”. Author Contributions F.J.R.G., R.M.N.C., A.P.L. and A.V. designed the work idea and established the main hypothesis, F.J.R.G., R.M.N.C. and A.P.L. performed the sample collection campaign; F.J.R.G., A.V. and C.M.R. designed the sample treatment and the laboratory workflow, F.J.R.G. and C.M.R. performed the laboratory work and the data analysis, under the supervision of W.O.; F.J.R.G. and R.M.N.C. performed the data interpretation and wrote the manuscript with the help of all the rest of co-authors. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-41804-y. Competing Interests: The authors declare no competing interests. Publisher’s note: 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 Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2019 Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 16 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

Assessment of functional and structural changes of soil fungal and oomycete communities in holm oak declined dehesas through metabarcoding analysis

Loading next page...
 
/lp/springer-journals/assessment-of-functional-and-structural-changes-of-soil-fungal-and-NsnybCs5Kr

References (114)

Publisher
Springer Journals
Copyright
Copyright © 2019 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
DOI
10.1038/s41598-019-41804-y
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports opeN Assessment of functional and structural changes of soil fungal and oomycete communities in Received: 21 December 2018 holm oak declined dehesas through Accepted: 15 March 2019 Published: xx xx xxxx metabarcoding analysis 1 1 2 Francisco J. Ruiz Gómez , Rafael M. Navarro-Cerrillo , Alejandro pérez-de-Luque , 3 4 3,4 Wolfgang oβwald , Andrea Vannini & Carmen Morales-Rodríguez Forest decline is nowadays a major challenge for ecosystem sustainability. Dehesas, which consists of savannah-like mediterranean ecosystems, are threatened by the holm oak decline in the south-west of Iberian peninsula. Phytophthora cinnamomi is considered the main agent of holm oak root rot, but little is known about the relationship between diversity of soilborne microbial community and the decline syndrome of holm oak. It would be hypothesized that the changes in the structure and functionality of the soil microbiome might influence tree health status through changes in richness and diversity of beneficial organisms such as mycorrhizal species, or fungal plant pathogens such as Fusarium spp. or Alternaria spp. t otal DNA of soil samples from declined oak dehesas was extracted and analyzed through metabarcoding techniques, to evaluate the specific composition and diversity of the fungal and oomycete communities and their relationship with the disease symptoms. the fungal community included a wide range of pathogens and abundance of ectomycorrhizal key taxa related with low defoliation degree. Phytophthora cinnamomi and Pythium spiculum did not appear among the most abundant oomycetes, nor were they related directly to defoliation levels. Moreover, a particular taxon belonging to the genus Trichoderma was strongly correlated with the scarcity of pathogenic Phytophthora spp. the diversity and composition of fungal and oomycete communities were related to the severity of the decline symptoms. the metabarcoding study of microbiome represents a powerful tool to develop biocontrol strategies for the management of the holm oak root rot. Alien-invasive species and global change are nowadays major threats for ecosystems sustainability around the world. It is possible to find examples of different die-back syndrome of woody species devastating natural or semi natural ecosystem around the world related with these species . In the last decades, different pathogenic Phytophthora spp. and Pythium spp. are emerging as important invasive pathogen species, considering the envi- 2,3 ronmental losses they cause . Some examples are Phytophthora ramorum Werres in Sudden Oak Death, ae ff ct- ing tanoak (Notholithoarpus densiflorus (Hook. & Arn.) Manos, Cannon & S.H.Oh) and several oak species in 4 5 California , Phytophthora cinnamomi Rands. in eucalyptus and Jarrah ecosystem in Australia , or Phytophthora cryptogea Pethybr. & Laff., Phytophthora cambivora (Petri) Buisman and P. cinnamomi causing ink disease in chestnuts worldwide . Among the most relevant diseases caused on European woody species by Phytophthora spp., root rot has been identified since the 1990’s ae ff cting oaks . Sessile and pedunculated oaks (Quercus petraea (Matt.), Liebl. Departamento de Ingeniería Forestal, Laboratorio de Ecofisiología de Sistemas Forestales ECSIFOR- ERSAF. Universidad de Córdoba. Campus de Rabanales, Crta. IV, km. 396, E-14071, Córdoba, Spain. Área de Genómica y Biotecnología, IFAPA, Centro Alameda del Obispo, Avda. Menéndez Pidal s/n, Apdo, 3092, 14080, Córdoba, Spain. Fachgebiet Pathologie der Waldbäume, Technische Universität München. Hans-Carl-von-Carlowitz-Platz 2, 85354, Freising, Germany. Department for Innovation in Biological, Agro-food and Forest systems (DIBAF) –University of Tuscia, Via San Camillo de Lellis snc, Viterbo, 01100, Italy. Correspondence and requests for materials should be addressed to F.J.R.G. (email: g72rugof@uco.es) Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 1 www.nature.com/scientificreports www.nature.com/scientificreports/ and Quercus robur L.) from North and Central Europe are threatened by several Phytophtora spp. , and in the Mediterranean basin, the decline of holm oak (Quercus ilex L., subsp. ballota in Spain and Portugal, and subsp. ilex in central Italy) and cork oak (Quercus suber L.) due to root rot were identified since the early 1980’s, ae ff cting 9–11 dehesa rangelands in Spain and Portugal, and wide forest areas in France and Italy . Dehesas represent ecosystems with relevant ecological and socio-economic importance covering 23% of the forested area in Spain (around 14% of the total area), accounting for over 3 mill ha in the Iberian Peninsula . In general, dehesa is a traditional land use in areas excluded from intensive agriculture, covering large extension of territory mostly characterized by a typical Mediterranean climate and soils usually shallow, acidic and deficient in nutrients . This agroecosystem is mainly composed by open woodlands of Quercus ilex and Q. suber, with an annual grassland understory, presenting canopy covers ranging from 5 to 60% . Dehesa ecosystem provides important socioeconomic benefits and many ecosystem services, such as shelter for a high degree of biodiversity or soil erosion prevention, being considered a first line barrier against desertification . Tree mortality of holm oak (Quercus ilex L.) associated with root rot is a serious threat to the sustainability 16,17 of these ecosystems. Phytophthora cinnamomi Rands. is considered the main agent of this disease , but the inu fl ence of other biotic and abiotic factors - such as the increasing frequency of extreme drought episodes and decreasing precipitation associated with global warming - is also recognized . Research efforts on holm oak decline have been focused on the interaction between the host plants and the pathogen – Q. ilex and P. cin- 19–21 namomi , and on their relationships with environmental and human factors, including socio-economic and 22,23 management strategies . Management strategies for the control of disease spreading are oen diffic ft ult to imple- ment because most of the dehesas yield is based on cattle livestock or hunting game . On the other hand, other treatments such as calcium soil amendments and phosphorous fertilization influenced tree health status but are 24,25 unable to suppress the soil inoculum . The application of sinigrin from Brassicaceae cultivars has been shown to be effective for pathogen suppression in controlled conditions experiments, but the biomass production nec- essary for an effective bio-amendment on the field is far greater than the productive ability of the dehesa soils . However, the relationship between the plant community and oomycetes and the spatial patterns of Pythium spp. and Phytophthora spp. were shown to be relevant factors influencing the dynamics of declining Mediterranean forests . Moreover, the implication of other Phytophthora spp. in the decline of Q. ilex and other Quercus spp. in 28–30 the Iberian Peninsula is recognized . In addition, little is known about the role of the soil fungal community in the decline syndrome of holm oak, or about the relationship between the latter and the presence of pathogenic species belonging to oomycetes. Die ff rent authors have highlighted the importance of the interactions among communities of soil taxa, mainly 31,32 focusing on bacterial and fungal communities . The soil fungal community interacts with other communities of biota in different ways, acting as commensal, pathogenic, antagonistic or mutualistic organisms. In the case of the plant community, its interaction is mainly associated with beneficial species such as ectomycorrhizal (EcM) and mycorrhizal (AM) fungi and other soil saprobes which improve soil fertility and structure and enhance plant defenses and health status . Ectomycorrhizal and endophytic species have been found to be related to changes 22,34–36 in the tree health status of several broad-leafed species, including Q. ilex . On the other hand, pathogenic relationships between fungal species and the plant community, including asymptomatic ones, are also related to disease expression in some cases . Oomycete communities are composed of a wide range of plant pathogens, some of which are host-specific and 38,39 others with a wide range of hosts. Their lifestyle can be either biotrophic, necrotrophic and hemibiotrophic . Phytophthora spp. communities have often been found to be the dominant in forest soils with a known history of Phytophthora-related diseases , but little is known about the relationship between the oomycete community and soil fungal or bacterial communities. Sapkota and Nicolaisen published the first study of the soil commu- nity diversity, including the taxa of fungi and oomycetes, associated with cavity spot of carrot, but in the litera- ture, there are no studies describing fungal diversity or the relationship between oomycetes and fungi in dehesas rangeland. e i Th solation of living cultures was the traditional way to assess oomycete or fungal diversity in soils, but this approach seems to be insufficient due to the time and specialized scientific knowledge required for phenotypic species identification , considering the potentially enormous number of different species – 1.5 to 3 million fun- gal species – in a high number of samples. In this sense, next-generation sequencing (NGS) techniques are an alternative methodology, able to detect the presence of a very high number of different taxa in a single sample: 40,44 successful examples include pyrosequencing and second-generation techniques, such as metabarcoding based 41,45,46 on Illumina de novo sequencing . The ability of metabarcoding techniques to characterize high-complex soil microbiome systems has represented an extraordinary improvement on the study of the plants and associated microbiota interaction. These techniques also oe ff r new opportunities for the deep assessment of microbiome ecology, given an important tool to optimize integrated control methods on management of plant diseases . The main objective of our work was to assess and describe the diversity and structure of fungal and oomy- cete communities in the holm oak dehesas rangeland ecosystems of Andalusia (southern Spain) ae ff cted by oak decline, using a metabarcoding approach, to analyse the correlations of structure and functionality of these com- munities between them, and with the defoliation level of Q. ilex. We hypothesize that the interaction between the microbiota diversity and functionality must be related with changes in tree health status, considering that fungi and oomycete are the main species involved in the holm oak decline due to root rot. To reach our main objective, the following specific objectives were addressed: i) to compare the fungal and oomycete diversity, and two defoli- ation levels in holm oak stands, in four different study zones; ii) to identify the predominant functional guilds of the fungal and oomycete communities and their relationships; and iii) to identify key taxa which might be related to holm oak decline. These results should contribute to a better understanding of the influence of soil community associations on tree health status and to the identification of important fungal and oomycete taxa related to holm oak decline in dehesas ecosystems. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 2 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 1. Mean relative OTUs abundance in 26 plots of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Abundance for fungi OTUs classified by “class” taxonomic level. (b) Abundance for oomycete OTUs classified by “genus” taxonomic level. Legends are referred to detail graph (left, “Other” pie chart) ordered by frequency. UI: Unidentified (for those OTUs which classification at the chosen taxonomic level was impossible). IS: Incertidae Sedis. The number of OTUs corresponding to each class (a) or genus (b) is showed in Supplementary Information, Table S1. Results otUs clustering and taxonomy. After bioinformatic treatment of the raw data, 1789 fungal OTUs (Operational Taxonomic Units) and 178 oomycete OTUs were confidently assigned to the genus or species level. Among the fungal taxa, Basidiomycota was the most abundant division (46.8%) followed by Ascomycota (45.1%) and Zygomycota (7.4%). e Th class Agaricomycetes was the most abundant (42.87%) and, together with Sordariomycetes, Eutriomycetes, Mortierellomyetes, and Pezizomycetes, represented the five dominant fungal classes (Fig. 1a). The class Agaricomycetes had its maximum relative abundance in the S.Nor zone (53.2%) and the lowest in Arac (39.4%). The opposite tendency was shown by Mortierellomycetes, which presented their minimum value in S.Nor (3.4%) and their maximum in Andv (11.2%) with a great value in Arac (6.8%) (Fig. 2a). No significant differences appeared for fungal classes in relative abundances between defoliation classes (Fig.  2b). Other relevant classes of Ascomycota were Dothideomycetes (5.36%) and Leotiomycetes (1.15%). Regarding oomycetes, species of the genera Phytophthora and Pythium spp. represented over 95% of the over- all clustered sequences, followed by Elongisporangium spp. (corresponding to E. undulatum), with almost 3% of the total sequence abundance, all of them corresponding to plant pathogen species (Fig. 1b). Moving eastward across the plots, the abundance of Phytophthora spp. increased and that of Pythium spp. decreased, with a sig- nificant inverse correlation between the relative abundances of the two taxa (Spearman ρ = −0.966; p < 0.001). However, none of the analysed oomycete taxa presented significant differences in abundance among the defolia- tion levels (Fig. 2c,d). The feature classification analysis of the microbiome selected 20 Oomycete OTUs and 349 Fungal OTUs as core features. Among the oomycetes, the most-abundant OTU, present in all the samples with impor- tant differences in abundance, was Phytophthora plurivora T. Jung & T.I. Burgess (oomycete OTU#1) - which accounted for 15.2% of the total frequencies, except in the case of Arac, where OTU#6 (Pythium paroecandrum Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 3 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 2. Relative abundance of oomycete and fungal taxa across study zones and mean defoliation class on the 26 studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. Vertical axis represents % of frequencies. Fungal abundance for zone (a) (%) and mean defoliation class (b) (%) is represented for taxonomic level 3 (class). Oomycete abundance for zone (c) (%) and mean defoliation class (d) (%) is represented for taxonomic level 6 (genus). Andv: Andévalo. Arac: Aracena. S.Nor: Sierra Morena & Sierra Norte. Pedr: Pedroches. UI: Unidentified taxon. Drechsler) was the most abundant, representing over 21% of the total frequencies in this zone (Table 1). The fourth-most-abundant oomycete was OTU#3, identified as Phytophthora quercina T. Jung & T.I. Burgess, whose presence was greatest in the Andv zone. Phytophthora cinnamomi (OTU#4) was the ninth-most-abundant OTU, accounting for 2.2% of the total sequences and identified in 18 of the 26 selected plots. Its presence was greatest in Pedr, where it was the third-most-abundant species (3.5% of the total frequency), and in Andv: in both cases, it showed lower abundance and frequency than P. plurivora and P. quercina. However, despite the low recovery rate on high-defoliated plots (C2, 38%), the abundance of P. cinnamomi in the high defoliated plots in which appeared was significantly high (31.8%). Two OTUs were identified as Pythium spiculum Paul. (OTU#21 and OTU#235), but they were present in only 21.2% of samples, corresponding to 10 plots, with low abundance (1.4% in total), and neither of them were classified as core features. Among the 349 fungal OTUs present in almost half of the soil samples, fungal OTUs #3 and #11, correspond- ing to Mortierella elongata Linnem., were the most abundant taxon, accounting over 6.1% of the total number of sequences (Table 1). The second most abundant OTU was OTU#2, Fusarium spp., present in all the samples ana- lyzed and representing over 3.4% of the total frequency. The third-most-abundant taxon was Russula praetervisa Sarnari (OTU#8), present and abundant in all the samples of Andv and S.Nor, but representing only 0.4% of the total abundance in Arac, and identified only in 75% of the samples of Pedr. Its abundance varied significantly with the mean defoliation (F = 3.38; p < 0.01), being higher in plots with low defoliation average (C1). Two OTUs belonging to the genus Trichoderma were classified among the 10 most-abundant species of the core biome: OTU#51, with its maximum frequency in samples from Arac, and OTU#15, representing over 1% of the total number of sequences in samples from Andv and S.Nor. Both OTUs had lower frequency rates in Pedr. Spearman’s correlation analysis showed a significant negative correlation between the abundance of Trichoderma spp. (fungal OTU#51) and that of 10 different OTUs of the oomycete core biome (Table  2), includ- ing most of the Phytophthora spp., except for OTUs #3 (P. quercina), #48 (Phytophthora sp.) and #12 (P. syringae). Regarding Pythium spp., only OTUs #6 (P. paroecandrum) and #2 (Pythium sp.) presented significant negative correlation. Despite its importance in our results, it was impossible to identify OUT#51 at species level because the well-conserved ITS1 region of Trichoderma spp. genus is unable to differentiate between species. Functional guilds. Of the OTUs confidently classified in an ecological guild (Table 3), the Soil-Plant Saprotroph guild was the most abundant in all cases (over 15%), except in S.Nor, where Plant Pathogen was the most-abundant guild (over 21%). However, when the frequencies of all the putative plant pathogenic OTUs were summed (including Plant Pathogen, Opportunistic Pathogenic Species and those OTUs classified in >1 Guild, which includes Plant Pathogen as one of those guilds), their abundance was greater than that of the rest of the guilds, in all the zones studied. Five functional guilds showed significant differences among zones, the Ectomycorrhizal (EcM) and Endophyte guilds exhibiting opposing trends, with the minimum frequency values of EcM species corresponding to the zone with maximum values for Endophyte. In the Pedroches zone (Pedr) the maximum values were found for the Endophyte, Animal Pathogen, Dung Saprotroph, and Soil-Plant Saprotroph functional guilds, with a Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 4 www.nature.com/scientificreports www.nature.com/scientificreports/ Location Average Defoliation Andevalo S. Aracena S. Norte Pedroches Low (C1) High (C2) (n = 16) (n = 10) (n = 12) (n = 14) (n = 28) (n = 24) Total (n = 52) Total % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. % of Freq. Total OTU ID Taxon Functional Guild (%) samples (%) samples (%) samples (%) samples (%) samples (%) samples (%) samples s__Phytophthora OTU_1 Plant Pathogen 11.1 100 17.5 100 19.7 100 7.1 100 20.9 100 7.7 100 15.2 100 plurivora s__Pythium OTU_6 Plant Pathogen 6.0 63 21.4 80 0.5 58 0.4 71 11.0 75 11.4 63 7.8 67 paroecandrum s__Phytophthora OTU_10 Plant Pathogen 1.1 69 5.8 40 11.5 75 1.8 57 7.1 67 2.7 50 6.0 62 psychrophila s__Phytophthora OTU_3 Plant Pathogen 5.5 100 2.7 90 6.1 75 4.0 86 6.0 92 2.4 81 4.7 89 quercina s__Phytophthora OTU_9 Plant Pathogen 1.5 69 2.9 50 6.7 75 0.4 100 5.2 89 0.5 50 3.5 75 citrophthora Oomycetes s__Pythium OTU_17 Plant Pathogen 5.2 56 4.1 80 1.4 75 1.3 86 3.8 78 10.9 75 3.2 73 heterothallicum s__Pythium OTU_15 Plant Pathogen 8.7 81 1.7 50 0.7 83 0.1 64 2.7 69 9.5 81 3.0 71 irregulare s__Elongis- OTU_7 porangium Plant Pathogen 3.9 100 4.6 70 0.4 83 1.3 100 2.9 92 0.8 94 2.6 90 undulatum s__Phytophthora OTU_4 Plant Pathogen 2.0 56 1.9 30 2.1 50 3.5 64 2.1 58 31.8 38 2.2 52 cinnamomi OTU_16 g__Pythium sp. Plant Pathogen 1.4 75 3.2 30 1.4 67 0.1 57 2.0 64 0.4 50 1.7 60 Total oomycete core biome frequency 46.4 65.8 50.5 20.0 63.8 78.2 49.9 s__Mortierella Soil-Plant Saprotroph OTU_3 7.5 100 4.7 100 2.8 100 1.8 100 4.7 100 3.9 100 4.4 100 elongata - Endophyte Plant Pathogen-Soil- OTU_2 g__Fusarium sp. 4.2 100 3.7 100 3.7 100 1.9 100 3.3 100 3.7 100 3.4 100 Plant Saprotroph s__Russula OTU_8 Ectomy-corrhizal 2.8 100 0.4 90 2.3 100 1.1 79 1.5 92 1.9 94 1.6 92 praetervisa s__Mortierella Soil-Plant Saprotroph OTU_11 2.8 100 1.1 80 2.0 75 6.0 100 1.4 86 1.0 100 1.3 90 elongata - Endophyte Fungi OTU_24 g__Solicoccozyma sp. Unassigned 1.0 100 0.9 100 0.8 100 1.7 100 1.1 100 1.1 100 1.1 100 OTU_9 g__Thelonectria sp. Soil-Plant Saprotroph 1.4 100 1.1 100 0.8 100 0.6 100 0.8 100 1.4 100 1.0 100 OTU_25 g__Penicillium sp. Soil-Plant Saprotroph 0.6 100 1.3 100 0.7 100 0.7 100 1.0 100 0.6 100 0.9 100 OTU_51 s__Trichoderma spp. Soil-Plant Saprotroph 0.6 100 1.1 100 0.5 92 0.7 100 0.8 100 0.7 94 0.7 98 s__Clonostachys OTU_48 Soil-Plant Saprotroph 0.6 100 1.3 100 0.2 92 0.3 100 0.6 97 0.8 100 0.7 98 rosea OTU_15 s__Trichoderma spp. Soil-Plant Saprotroph 1.1 94 0.2 100 1.2 92 0.5 93 0.9 92 0.3 100 0.7 94 Total fungal core biome frequency 22.6 15.8 15.0 15.3 16.0 15.2 15.8 Table 1. Ten most abundant oomycete and fungi OTUs found in 26 Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems, ordered by relative frequency (from core-features of each abundance matrix calculated among OTUs present in at least 50% of samples). Freq.: Number of sequences belonging at each OTU relative to the total number of sequences for the considered factor. % of samples: Percentage of samples in which the OTU was identified. significantly-lower value for EcM. Only the Animal pathogen guild showed significant die ff rences among defoli - ation levels, and no Endophyte species was identified for high defoliated ones. e S Th pearman’s correlation matrix between functional guilds (Fig.  3a) shows a significant, positive correlation of Soil-Plant Saprotrophs with Animal Pathogen abundance (ρ = 0.491, p < 0.01) and a negative one with Plant Pathogens (ρ = −0.356, p < 0.05) (Fig. 3b-1 and 2). On the other hand, EcM was negatively correlated with Soil-Plant Saprotrophs (ρ = −0.559, p < 0.01; Fig. 3b-3), Animal Pathogen (ρ = −0.492, p < 0.01) (Fig. 3b-4), Dung Saprotrophs (ρ = −0.384, p < 0.01), Opportunistic Pathogenic Species (ρ = −0.348, p < 0.05) and >1 Guild, (ρ = −0.343, p < 0.01). The Dung Saprotroph guild was positively correlated with Soil-Plant Saprotrophs (ρ = 0.332, p < 0.01) and with the species assigned to more than one guild (>1 Guild, ρ = 0.392, p > 0.01). Other significant correlations appeared between > 1 Guild OTUs, Lichenized and Opportunistic Pathogenic Species. Fungal and oomycete diversity. The Good’s Coverage indices for the selected thresholds were high for the rarefied matrices of the oomycetes (G = 0.98) and fungi (G = 0.99), indicating that over 98% of the identified OTUs were sampled in the rarefaction process and thus included in the α-diversity analysis. e S Th hannon H’ and Pielou Evenness α -diversity indices of the oomycetes exhibited similar trends according to the geographical longitude, their abundances decreasing from west to east (Fig. 4), with significant differences Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 5 www.nature.com/scientificreports www.nature.com/scientificreports/ Rank OTUs Assigned Taxon Ρ 1 OTU#1 P. plurivora −0,392** 2 OTU#6 P. paroecandrum −0,263* 3 OTU#10 P. psychrophila −0,272* 5 OTU#9 P. citrophthora −0,378** 8 OTU#7 E. undulatum −0,325* 9 OTU#4 P. cinnamomi −0,345** 12 OTU#5 P. cactorum −0,277* 13 OTU#8 P. gonapodyides −0,308* 14 OTU#2 Pythium sp. −0,395** 15 OTU#14 P. cambivora −0,379** Table 2. Significant correlations between the abundance of OTU#51 (Trichoderma spp.) and the abundance of the core biome OTUs of oomycete dataset at 26 Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. Rank: Order of the OTUs in the core features classification by abundance (1 for the most abundant, to 20 for the less abundant). ρ: Spearman’s Rho value. *Significant correlation at p < 0.05; **Significant correlation at p < 0.01. Location Defoliation sign. Functional Guild Andv Arac S.Nor Pedr Total C1 C2 t (df ) n/s n/s >1 Guild designation 3.3 ± 0.1 3.6 ± 0.2 3.8 ± 0.2 4.0 ± 0.3 3.6 ± 0.1 3.7 ± 0.1 3.6 ± 0.2 0.46 (50) b ab b a Animal Pathogen 0.9 ± 0.1 1.0 ± 0.1 0.8 ± 0.1 1.3 ± 0.1 1.0 ± 0.1** 0.9 ± 0.1 1.2 ± 0.1 −2.10 * (52) n/s n/s Arbuscular Mycorrhizal 3.4 ± 0.5 3.3 ± 0.4 5.1 ± 0.9 3.5 ± 0.9 3.8 ± 0.3 3.9 ± 0.4 3.5 ± 0.5 0.44 (50) ab c bc a n/s Dung Saprotroph 1.0 ± 0.1 0.7 ± 0.1 0.8 ± 0.1 1.1 ± 0.1 0.9 ± 0.0*** 0.9 ± 0.1 0.8 ± 0.1 0.61 (50) a a a b n/s Ectomycorrhizal 9.9 ± 0.7 9.8 ± 0.9 9.4 ± 0.4 5.5 ± 0.6 8.7 ± 0.4*** 9.0 ± 0.5 8.0 ± 0.8 1.08 (51) † b b a Endophyte 0.2 0.14 ± 0.02 0.12 ± 0.02 0.4 ± 0.1 0.22 ± 0.04** 0.22 ± 0.04 — — n/s n/s Fungal Parasite 0.24 ± 0.02 0.31 ± 0.04 0.4 ± 0.1 0.29 ± 0.04 0.30 ± 0.02 0.29 ± 0.03 0.33 ± 0.03 −0.81 (49) n/s (††) n/s Lichenized 0.3 ± 0.0 0.25 ± 0.04 0.24 ± 0.03 0.19 ± 0.04 0.25 ± 0.02 0.24 ± 0.02 0.28 ± 0.06 −0.64 (12.6) n/s n/s Oportunistic species 3.7 ± 0.2 3.4 ± 0.2 3.9 ± 0.2 4.1 ± 0.3 3.8 ± 0.1 3.8 ± 0.1 3.7 ± 0.3 0.52 (54) n/s n/s Plant Pathogen 10.0 ± 1.0 14.1 ± 5.9 21.7 ± 8.4 13.4 ± 6.1 14.5 ± 2.8 15.6 ± 3.4 11.9 ± 5.0 0.59 (50) ab ab b a n/s Soil-Plant Saprotroph 15.7 ± 0.5 15.4 ± 0.5 14.0 ± 0.5 17.1 ± 0.7 15.6 ± 0.3** 15.3 ± 0.4 16.4 ± 0.5 −1.71 (50) n/s n/s Unclassified 51.7 ± 0.7 50.5 ± 3.4 45.2 ± 4.9 51.4 ± 4.0 49.8 ± 1.7 49.0 ± 2.0 52.0 ± 3.2 −0.79 (54) Table 3. Mean ± SE relative abundance of functional guilds between study zones and defoliation levels of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. e d Th ata of functional guilds which presented significant differences (ANOVA with Bonferroni mean differences test, or Kruskal-Wallis rank sum test) were written in Italic. Andv: Andévalo. Arac: Aracena. S.Nor: Sierra Morena & Sierra Norte. Pedr: Pedroches. C1: Low average defoliation plots. C2: High average defoliation sign plots. t = t-value from the Student test with the degree of freedom and the significance level of the (df ) test.  Only one sample of Andevalo presented confident identification of endophyte organismSignificance level: *p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001; n/sp > 0.05. Values with the same letter are not significantly different (††) (p > 0.05).  Levene test p < 0.05 (variances were different between factors). between Pedr and Andv (H : 5.69, p < 0.01; H : 6.06; p < 0.01), but without differences in OTUs richness. Shannon Evenness e f Th ungal community did not show significant differences for the α -diversity indices. When these indices were compared between the mean defoliation levels for both communities, only oomycete evenness was significantly lower for the greatly-defoliated plots (χ : 4.64, DF: 1, p < 0.01); the fungal community did not show differences. Regarding β-diversity, ANOSIM showed the influence of location in the composition of the fungal commu- nity (R = 0.167; p < 0.001), whereas for the oomycete community (R = 0.078; p < 0.05) differences were only found between Pedr and Andv (R = 0.214; p < 0.01) and between Pedr and S.Nor (R = 0.125; p < 0.05). The non-linear metrics analysis (NMDS) showed a more-aggregated zonal cluster for the fungal community, in comparison with the oomycetes (Fig. 5), but also a clearer separation of the high-defoliated plots for the oomycete community. Oomycete OTU#1 (P. plurivora) had the vector which most clearly influenced the separation between defoliation classes. Discussion To the best of our knowledge, this work is the first study of soil biodiversity in holm oak dehesas ecosystems, exploring the interaction between fungal and oomycete communities and their influence in the decline symp- toms, focusing on the soil microbiome. Phytophthora spp. dominated the oomycete community, but the species usually related to Q. ilex root rot (P. cinnamomi and P. spiculum) did not appear among the most abundant, nor were they related directly to defoliation levels. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 6 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 3. Relationships between functional guilds on the 26 studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Correlation matrix between functional guilds relative abundance in each sample. *Significant correlation at p < 0.05; **significant correlation at p < 0.01. (b) Scatterplots of significant correlations between Functional Guilds: (1) Plant Pathogen/Soil-Plant Saprotroph; (2) Animal Pathogen/Soil-Plant Saprotroph; (3) Ectomycorrhizal/Soil-Plant Saprotroph; (4) Ectomycorrhizal/ Animal Pathogen. ρ: Spearman’s Rho value. p: Significance level of correlation. Figure 4. Analysis of Shannon H’ diversity, (left), OTUs richness (centre) and Pielou-e eveness vector (right) by geographical zones, for the studied Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems, for fungi (up) and oomycete (down). It must be considered as a limitation of our study the lack of specific data from the rhizosphere. Without this information, the discussion on the shifts of an important guild like Endophytes could not be addressed, but despite holm oak rhizosphere was not analysed, the revision of the ecological roles of the key taxa identified in this work lead us to believe that most of them should be involved in the plant-microbiota interactions, agreeing 48,49 with the results of other authors . Our results showed that the fungal community of declined dehesas includes a wide range of plant pathogens and key taxa of EcM and Trichoderma spp. which would be related to the tree health status. The metabarcoding analysis revealed that structural changes of the fungal community and the diversity level of the oomycete community were related with the mean defoliation of the holm oak dehesas ana- lyzed, agreeing with our work hypothesis. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 7 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 5. Non-metric multidimensional scaling analysis of OTUs abundance influence in defoliation levels and site location, for the 26 studied plots on Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. (a) Analysis for fungi dataset. (b) Analysis for oomycete dataset. Ellipsoid indicates normal probability of point distribution at 90% confidence level for mean defoliation classes of plots. Arrows show influence vectors for the ten most abundant OTU’s of fungi and oomycete dataset. the fungal community of dehesas has a wide range of plant pathogens and abundance of e cM. Despite the large number of fungal phytopathogens described in the studied soils, it is not possible to determine the direct influence of these species on the health status of the trees. For example, some of the fungal patho - gens identified in our samples have been found previously in dehesas as a part of the endophytic mycobiota of pasture species, such as species from the genus Colletotricum spp., Fusarium spp., Sordaria spp. or Alternaria spp. . However, the Dothideomycetes Alternaria alternata and Pyrenochaeta cava were isolated from Quercus robur roots by Kwaśna et al. , presenting significant changes in abundance according to the degree of tree defo- liation. Therefore, the increment in the abundance of some fungal pathogens, such as Dothideomycetes, in the high-defoliated plots, might be due to the pathogenicity of some species, but also to changes in the functionality of mycobiota related to the understory or grassland species. Mortierella elongata is a saprotroph often isolated in natural soils, being related their abundance with soil 52,53 improvement in nutrient quality and fertility . It plays an important role in the degradation of several com- pounds, including toxic ones , and its abundance can be considered an indicator of the resilience level of dehesas ecosystems. Of the symbiotic species, Russula was the most-abundant genus, with 76 different OTUs classified within it, including the third-most-abundant species, R. praetervisa. Ectomycorrhizal abundance and bio- diversity have been related to fine root vitality in Quercus robur forests, being considered the most-important below-ground trait related to their decline . Also, its abundance has been related to benefits for the tree health status and the establishment of Quercus spp. seedlings . In our work, we did not find a relationship between Russulaceae abundance and the individual pathogenic Phytophthora spp., but statistically-significant differences were found when R. praetervisa abundance was analysed in relation to defoliation, both results agreeing with 22,57 other studies in Q. ilex forests . In our search for plant pathogens within the core biome of the fungal community, the genus Fusarium was the most relevant. This genus includes several soil-borne plant pathogens, related to forest and agricultural plant diseases in woody, shrub and herbaceous hosts . The F. oxysporum species complex causes vascular wilts, damping-off and crown and root rot in a wide range of species, being able to produce the death of holm oak seedlings in a pathogenicity test . The role of F. eumartii in the decline of Q. robur, Q. cerris and Q. pubescens 60,61 in Italy has been recognised , and other species of this genus have been linked with diseases in other Quercus 62–64 spp. . In addition to Fusarium spp., 39 more OTUs within the fungal community core biome were assigned to plant pathogenic species, including some species known to be pathogens of oaks - such as some Peniophora spp. (among which, OTU#816 was identified in 98% of samples), Pyrenochaeta spp., Thelonectria spp., Alternaria alter- 66 51 67 nata , Sporothrix ina fl ta, Pezicula radicicola and Phomopsis spp. . The abundance of fungal pathogenic species in zones with a low H’ index for oomycetes, such as Pedr and S.Nor, coincided with a lower presence of beneficial fungal species, such as strict saprotrophs, AM or EcM. These circumstances could lead to an environment that is ideal for triggering the decline of trees. Phytophthora spp. dominate the oomycete community in declining dehesas. Phythophthora was the dominant genus among oomycetes, being the most-abundant species P. plurivora, instead of those related to 16,17,68 holm oak root rot according to the literature – P. cinnamomi and P. spiculum . Particularly, Pythium spiculum, considered a trigger of holm oak decline and mortality , was excluded from the oomycete core biome. Among the oomycete pathogens causing root rot in Quercus ilex, P. psychrophila, by its abundance, and P. quercina, by its frequency of appearance, were the most-important taxa, P. cinnamomi appearing as the ninth-most-abundant Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 8 www.nature.com/scientificreports www.nature.com/scientificreports/ oomycete, present in 52.7% of samples. However, P. cinnamomi presented a very high frequency in the samples of high-defoliation plots (C2) in which it was identified, agreeing this result with the great susceptibility that holm oak presents against this pathogen . Phytophthora plurivora is a very-common soil-borne pathogen that has been isolated from soils of many natu- ral forests in Europe. It was recorded as the most-abundant species in a survey of oak stands in Italy and in chest- nut forests of central Italy, in a study of soil Phytophthora biodiversity through metabarcoding . Its pathogenicity 71 72 73 on Quercus alba , Q. robur, Q. petraea, Q. rubra and Q. cerris is recognised. However, to our knowledge, there is no evidence of a direct pathogenic relationship with Q. ilex root rot, or of its involvement in the decline syndrome in dehesas, the work of Vettraino et al. being the only one that recorded this species (designated as P. citricola) under Q. ilex trees (subsp. ilex). The abundance of this species in our samples was significantly high in low-defoliated stands, which influenced NMDS aggregation of defoliation classes. The high recovery rate of P. plurivora and its relationship with the defoliation encourage the study of their involvement in holm oak decline. Most of the Phytophthora taxa found in our samples were also recorded in a soil analysis in the search for oomy- cete in natural Q. ilex forests of eastern Spain , where P. psychrophila and P. quercina were the most-relevant taxa, and also a qPCR survey showed P. quercina abundance in dehesas of Extremadura ; P. psychrophila and P. quercina appeared, respectively, as the third- and fourth-most-abundant oomycete species in our study. Both taxa have 74 29 been confirmed as Q. ilex root pathogens, causing root rot and tree death . However, in the work of Català et al. , no match was found for P. plurivora and P. cinnamomi. The overall specific composition of the Phythophtora spp. community in our samples was more similar to that described by Vettraino et al. , who found P. plurivora, P. quercina, P. cinnamomi and P. cactorum – in that order – as the most-frequent Phytophthora spp. isolated in stands of Q. ilex ae ff cted by oak decline. Regarding the severity of symptoms associated with the oomycete community, no significant differences in 22,57 taxa abundance were found among the defoliation categories, agreeing with the works of Corcobado et al. . Notwithstanding, defoliation and the severity of root rot should be related ; in our case, this relationship was detected for the overall structure and functionality of the oomycete community (diversity and evenness), includ- ing the influence of the fungal community (EcM abundance). Given our results, the involvement in the oak decline of all the plant pathogens present in the dehesa soils should be considered, not only that of P. spiculum and P. cinnamomi, but also other pathogenic Phytophthora species - such as P. plurivora, P. quercina and P. psychrophila - and other fungal plant pathogens such as Fusarium spp. or other phytopathogenic taxa belonging to the class Dothideomycetes. ecM and AM are key guilds related to shifting soil microbiota functionality in dehesas ecosystems. The Ectomycorrhizal (EcM) and the Arbuscular Mycorrhizal (AM) guilds, the third- and fourth-most-abundant functional guilds, respectively, had lower abundance rates overall when compared with other works in soils of dif- 45,46 ferent natural ecosystems . Quercus ilex is a species that oen p ft resents mycorrhizal association with different species of Ascomycetes and Basidiomycetes . A significant reduction in EcM and increment in plant pathogens (including Plant Pathogen and Opportunistic Species guilds) were found in Pedr, coinciding with its status as the most-high-defoliated zone, being the higher EcM abundance correlated with lower pathogen and saprobe abun- 76 57 dance, agreeing our results with other works . Corcobado et al. identified non-mycorrhizal root tips of Q. ilex as vulnerable points for pathogen invasion. Beneficial mycorrhizae, including AM and EcM, can potentially induce 77 33 higher pathogen tolerance through stimulation of the plant immune system . Therefore, either the increase of the root uptake and the effect on plant immune response might be related with lower defoliation level of trees in zones with high EcM and AM abundance. The Animal Pathogen, Soil-Plant Saprotrophs, and Dung Saprotroph guilds had their maximum values in Pedr, which showed the lowest α-diversity values for oomycetes and an increment in the abundance of Phytophthora spp., with good values for the diversity of the fungal community, but lower EcM abundance. This is another example of the greater influence of the community structure and functionality, in comparison with the specific composition. It also supports the existence of a triple relationship among the abundance of EcM and other bene- 22,57 ficial fungi, the plant pathogens and the tree health status, agreeing with other works . High diversity and changes in community structure are associated with the severity of symptoms. 45,46,78 In comparison with recent analyses of the mycobiota in soils of different ecosystems , the dehesas had a high Shannon H’ index (greater than 4.5 in all the studied zones) and a high value of evenness (exceeding 0.7). This shows a high diversity and good equilibrium of the mycobiota, but with a high proportion of plant pathogens and significant changes in the specific composition of the soil community among zones, agreeing with the results of ANOSIM. No relationship of the studied parameters with plant community changes was found, agreeing with 45,46 other authors who reported that changes in the diversity of mycobiota were decoupled from plant α-diversity . e Th values of the α -diversity indices of the oomycete community were, in general, higher than those found by Català et al. in natural forest soils of eastern Spain and were similar to (slightly higher than) those reported by Vannini et al. for chestnut forests of Italy. The lowest values of the Shannon index and evenness vector for the oomycete community occurred in Pedr, the most-defoliated zone, showing low species richness and differences in taxa abundance, agreeing with the significantly-lower values of evenness in the highly-defoliated plots. Although a direct relationship between specific oomycete taxa, or their abundance, and the severity of symptoms was not found, the β-diversity analysis showed aggregation of high-defoliated trees regarding the oomycete community. Similar results were obtained by Linaldeddu et al. when evaluating the diversity of the fungal endophytic com- munity in declining Quercus suber L. trees. The severity of defoliation, therefore, seems to be related more to the functional structure of the oomycete community – driven by the evenness - than to its specific composition, being presented the change in the specific composition of the fungal community and the reduction in competition between root pathogens as the most influential drivers related with decline symptoms. Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 9 www.nature.com/scientificreports www.nature.com/scientificreports/ the abundance of Pythium spp. and Trichoderma spp. influencing the low Phytophthora spp. dominance. The abundances of Pythium and Phytophthora influenced each other, but it remains unclear how the changes were inu fl enced by location or by a natural equilibrium between the taxa. Our results show a significant decrease in oomycete α -diversity in soils of highly-defoliated stands, coinciding with an increase in the abundance of Phytophthora spp. Other authors have demonstrated that the mixture of P. cinnamomi with dif- ferent species of Pythium or other pathogenic fungi (such as Fusarium oxysporum) did not increase the severity of symptoms , supporting the hypothesis that, in soils with less diversity of oomycetes, P. cinnamomi and other pathogenic Phytophthora species have an increased ability to infect roots and ae ff ct the health status of the tree, increasing the severity of visible symptoms. Another key relationship between taxa was the significant inverse correlation between OTU#51 abundance (Trichoderma spp.) and most of the core biome taxa of oomycetes in the studied plots (Table 2). The effects of Trichoderma spp. in biocontrol strategies are well known . Trichoderma asperellum, T. hamatum, and T. virens have been isolated from holm oak roots (Q. ilex subsp. ilex) and reduced the root rot of seedlings in a controlled experiment including P. nicotianae and P. cinnamomi . Other Trichoderma spp. have also been proved to be 81–85 antagonists of diverse Phytophthora spp. . The evidence of the multiple effects of Trichoderma spp. against multiple phytopathogenic species should encourage further investigation of this relationship in the holm oak dehesas. The next steps in the search for control strategies based on the management of fungal soil biota should be the isolation of Trichoderma spp. from the studied soils to identify the fungus at species level, together with the analysis of root endosphere, which would provide useful information about the endophyte behaviour of the 86,87 present species, as it has been shown for other Trichoderma spp. . Metabarcoding is a powerful tool in the management of soilborne diseases. Biocontrol strategies are considered environmentally friendly approaches in the management of plant diseases, but the high complexity of ecological interactions between the microbiome components could lead to negative impacts on non-target organisms , and the misunderstand the mechanism of the interactions provokes the inefficiency, or even adverse effects of treatments . Thus, the study of the microbiome is an important step to develop management strategies to reduce tree decline damage through biocontrol strategies, and the knowledge about microbial community has demonstrated to have important implications in those aspects . Metabarcoding techniques could help to focus biocontrol strategies, searching for the adaptability and synergism of new products with the indigenous micro- biome, or characterizing microbiota of disease-suppressive soils and the role of key taxa in this microbiome. In this sense, further studies confirming the presence of the identified key taxa in the rhizosphere are encouraged. Additionally, it is recommended to complement metabarcoding approaches with other methods, beyond the rhiz- osphere microbiome analysis, to better understand the role of the most relevant taxa, including those taxa that are identified for the first time associated with this system. The identification of relationships between Trichoderma spp., EcM and Phytophthora spp. in soils related to holm oak decline leads to new soil management strategies in those stands ae ff cted by root rot, directly linked to the microbiome ecology. Moreover, biological control strate- gies must consider the influence of location and environmental variables in soil community , being the metabar- coding analysis an adequate technique to provide this information. Conclusion This work is a novel contribution to the study of soil biodiversity in holm oak dehesas ecosystems, analysing for the first time the fungal and oomycete communities together in dehesa soils. Our results revealed high levels of diversity in the fungal community and medium to low diversity of oomycetes, the presence of key taxa influencing holm oak health status and oomycete pathogenic abundance, and the dominance of plant pathogens and sapro- trophs in the analysed soils. The microbiome presented significant influence in the defoliation status, confirming our main hypothesis, and suggesting that both the structure of the fungal community and the level of diversity of the oomycete community were related with the tree health status. On the other hand, the high degree of soil microbiome diversity in declined dehesas revealed on this work encourages to a deeper research into the relationship between different taxa, considering a tripartite relationship (plant-pathogens-other taxa) , but also in the study of the influence of environmental factors in the soil com- munity, beyond the simple picture of a single plant-pathogen interaction. The findings about the presence of key taxa related with defoliation symptoms or pathogen abundance, must drive us to consider biocontrol strategies in the management of the holm oak root rot, through methods and applications that respect the local microbiota. Several biocontrol strategies could be explored deriving from our results, such as the application of Trichoderma spp. to control the inoculum of pathogenic oomycetes, as well as the application of mycorrhizae as a strategy to improve the health status of trees, but further research is needed to clarify how both key taxa interacts between them and with holm oak. Material and Methods study area. This study was focused on holm oak (Quercus ilex subsp. ballota) dehesas of north-west Andalusia (Fig. 6). The studied area was divided into four zones - Andévalo (Andv), Sierra de Aracena (Arac), Sierra Norte de Sevilla and S. Morena (S.Nor) and Valle de los Pedroches (Pedr) - located in the provinces of Huelva, Sevilla and Córdoba, representing the core area of the holm oak dehesa ecosystems of Andalusia (Fig. 6). This division was made based on differences in environmental factors and silvicultural characteristics (Table  4). A total of 26 plots were selected among the permanent environmental monitoring points of the Andalusian Network for Damage Monitoring in Forest Ecosystems (Red de Información Ambiental, Consejería de Medio Ambiente y Ordenación del Territorio - Junta de Andalucía, ICP Forest, EC Network level II - www.juntadeanda- lucia.es). The plots were selected based on the following criteria; (i) >85% of the trees were Q. ilex, (ii) a positive record of P. cinnamomi and (iii) the tree mortality rate. To calculate the mean plot defoliation, all the trees of the Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 10 www.nature.com/scientificreports www.nature.com/scientificreports/ Figure 6. Map of sampled plots of Quercus ilex dehesas of the Andalusian Network for Damage Monitoring in Forest Ecosystems. plots (24 trees × 26 plots) were assessed for defoliation following the methodology of Eichhorn et al. . The mean defoliation for each plot (%), calculated using the individual defoliation of the 24 trees per plot, ranged between 13.3% and 56.9% resulting in 14 plots classified as C1 (mean defoliation < 25%), and 12 plots classified as C2 (25% ≤ mean defoliation < 65%) (Supplementary Information, Table S2). The mean defoliation showed signif- icant differences among the study zones (χ = 11.1; DF = 3; p < 0.05), the highest value corresponding to Pedr. soil sampling. In each plot, soil samples were collected from the rhizosphere under the crown of two trees with defoliation degree similar to the mean defoliation of plot (28 samples under C1 trees and 24 samples under C2 trees, n = 52), searching for the representability of the mean plot conditions (soil, exposure and slope). For each tree selected, sub-samples were taken under the crown for two orientations (North and South) and two distances from the trunk (0.5 m and 1–2 m), making four sub-samples per tree. Then, the sub-samples of each tree were bulked and stored at 4 °C until processing. Samples were collected from November to December 2015. The soil samples were processed within the next 24 hours following their collection; soil was dried at room temperature for 24 h and well-mixed, litter and roots were separated, and 20 g of the fine fraction was collected aer p ft assing the soil through a 1-mm sieve. DNA extraction, ITS1 amplification, library generation, and sequencing. Total DNA was extracted from soil using the MoBIO Power Soil DNA Isolation kit (MoBIO Laboratories, Inc., CA, USA). Five replicates of fresh soil (0.33 g each) were taken from each soil sample. The aliquots were mixed together, and total DNA was purified and concentrated using a DNA purification kit (Wizard SV Genomic DNA Purification System; Promega Corporation, WI, USA). The DNA concentration was quantified, with a Qubit dsDNA BR Assay kit (Thermo Fisher Scientific Inc., MA, USA), prior to amplification. Additionally, two mock communities were created by pooling DNA from a pure culture of 13 different species of oomycetes and another with seven species of fungi (Supplementary Information, Table S3); one community of each taxa had the same concentration of each species (Oo= and F=) and the other had different concentrations (Oo≠ and F≠). These were used as controls and were processed in the same way as the samples, following the methodology of C. Morales Rodriguez (Unpublished). e in Th ternal transcribed spacer 1 (ITS1) was amplified in a multiplexing PCR using a set of barcoded primers: 92 93 ITS1F and ITS2 for fungi, and ITS6 and ITS7 for oomycetes . e Th primers were tagged with different barcodes to distinguish different samples. Moreover, 3 replicates of each control community were amplified with different barcoded primers of the ITS1F/ITS2 and ITS6/ITS7 primer pairs (Supplementary Information, Tables S4, S5). Negative controls and cross-reactions (Oo with ITS1F/2 and F with ITS 6/7) did not originate any PCR product. e PCR r Th eaction mixture consisted of 25 μl of Maxima Hot Start PCR Master Mix (Thermo Fisher Scientific, −1 −1 USA), 2 μM of each primer, 4 µ l of 2 mg ml BSA and 3 μl of DNA template diluted to 50 ng µ l , in a total volume of 50 μl. The thermal cycle was an initial 5-min step at 95 °C, followed by cycles (30 for fungal ITS and 35 for oomycete ITS) of 40 s of denaturation at 95 °C, 2 min of annealing (58 °C for fungal ITS amplification and 60 °C for oomycete ITS) and 1 min of extension at 72 °C, with a final extension step of 7 min at 72 °C. For each sample and the controls, three replicates of PCR for fungal ITS and five for oomycetes ITS were carried out and the products were pooled. Amplicons were purified using the MagJet NGS Cleanup and Size Selection kit (Thermo Fisher Scientific Inc., MA, USA) and the final concentration of purified PCR product was quantified using a Qubit dsDNA HS Assay kit (Thermo Fisher Scientific Inc., MA, USA). The resulting amplicons were pooled in equal amounts to give a total −1 volume of 40 µ L, with a concentration of 200 ng µ L , and were sent to Eurofins Genomics AT GmbH (Vienna, Austria), to be sequenced in an Illumina MiSeq 2 × 300 bp platform. Both fungal and oomycete libraries were Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 11 www.nature.com/scientificreports www.nature.com/scientificreports/ Andevalo Aracena Sierra Morena & S. Norte Pedroches Mean Precipitation (mm) 675.7 829.2 712.1 572 Mean Annual temperature (°C) 17.4 16.2 16.1 15.8 Mean Evapotranspiration (mm) 1445 1396 1458 1443 Max elevation (masl) 293 1046 936 943 Mean elevation (masl) 191 434 429 570 Average cover* Low Medium High Low Dominant tree age** Young to Mature Mature Young Old to very old Normal tree cover with High tree cover with Class of dehesa Cultivated dehesa Pastures shrubs and pasture shrubs and pasture Table 4. Mean values of environmental and silvicultural parameters considered for the selection of the study zones (m.a.s.l.: meters above sea level). Average values for the period 1969–2000, taken from the Andalusia Environmental Information Network - REDIAM; Consejería de Medio Ambiente, Junta de Andalucía (http:// www.juntadeandalucia.es/medioambiente/site/rediam/portada/). The average tree cover, the dominant tree age and the prevalent class of dehesa were determined as the trend of all the visited plots for each zone, following the 104 −1 −1 criteria of Costa Pérez et al. . *Low = N < 50 tree ha ; Medium = 50 < N < 80 tree ha ; High = N > 80 tree −1 ha . **Young trees = 50–90 years/20–40 cm Ø; Mature trees = 110–210 years/50–80 cm Ø; Old trees = 210– 250 years/90–110 cm Ø; Very old = > 250 years/Ø > 110 cm. included on the same Illumina run. The reads are available in the EMBL Nucleotide Sequence Database (http:// www.ebi.ac.uk/embl) under the project name “NGS_Dehesa” and accession number PRJEB28448. s equence processing, ot U selection, and abundance matrix construction. A bioinformatic pipe- 94 95 line was adapted from Bálint et al. and Sapkota & Nicolaisen . The primers were debarcoded and classified by Eurofins Genomic AT GmbH (Vienna, Austria). The sequences were joined (22 774 511 sequences) and reads with mismatches or ambiguous bases were unassigned (1 246 294 sequences). Finally, only paired sequences in which the 5′ barcode and forward primer and 3′ barcode and reverse primer were found, were used in the bioin- formatic analysis (Supplementary Information, Table S6). Singletons were excluded from the analysis. e f Th orward and reverse sequences of paired reads were trimmed, the primers were eliminated, and sequences were overlapped and paired using CLC Genomics Workbench 3.6.5 (https://www.qiagenbioinformatics.com/) (Supplementary Information, Table S7). The ITS1 region was extracted with the Fungal ITS Extractor version 2 , sequences were clustered in operational taxonomic units (OTUs) and the chimeric sequences were filtered - de-novo chimera filtering for oomycete files - using USEARCH 8.1 with a similarity level established at 98%. The consensus sequences of the OTUs were identified with the GenBank nucleotide database (nr/nt) of NCBI, using the Megablast algorithm, and the results were parsed, and the sequences retained in MEGAN 6.10.8 to assign taxonomic identifications to each OTU (Supplementary Information, Table S8). In some cases, ITS1 was una- ble to die ff rentiate between species belonging to genera which presented a high conserved ITS1 region, such as Trichoderma spp., but aer a t ft horough supervision of the taxonomy, the metabarcoding analysis showed enough accuracy to study the diversity and composition of the fungal and oomycete communities, clearly differentiating between genera in all cases. Ae ft rward, abundance tables of OTUs for oomycete and fungal files were constructed using USEARCH 8.1 with a level of identity of 98%. The amplified ITS sequences from all the species included in the internal controls were correctly identified. The mean frequency of debarcoded paired sequences did not die ff r between oomycetes and fungi ( χ = 0.727; p = 0.394), with strong correspondence between the abundance of OTUs and the estimated amount of DNA included in each control. Prior to downstream analysis of the abundance matrix, control samples and OTUs with less than 10 reads were eliminated from both the fungal and oomycete feature tables. Fungal and oomycete community analysis. Aer t ft he abundance matrix filtering, functional guilds of the identified taxa were characterized using FunGuild database (accessed on March, 2018) . Only OTUs classi- fied with a confidence level of “Probable” or “Highly Probable” were used. e Th functional guilds of the remaining OTUs, classified as “Unassigned” or “Possible”, were revised based on the literature; those that were confidently classified were added to the analysis (Supplementary Information, Table S9). The rest of the unassigned OTUs and those assigned with low confidence were placed in the “Unclassified” category. Taxa classified as soil sapro- trophs, plant saprotrophs, litter saprotrophs or wood saprotrophs were grouped in a unique Soil-Plant Saprotroph category. Those taxa that were confidently classified in more than one guild were assigned to the class > 1 Guild, except in the case of taxa that were confidently assigned to four or more guilds, including pathogen and sapro- troph guilds; these were classified in the Opportunistic Pathogenic Species class. To analyze the diversity of oomycetes and fungi, abundance matrices were processed using QIIME2 version 2017.12.1 . Rarefaction curves, Good’s coverage, Shannon H’, OTUs abundance and Pielou evenness vector and the Bray-Curtis distance matrix were calculated using the plugin Diversity 2017.12.0 . e Th most abundant and frequent OTUs, presented in at least 50% of the samples, were selected and classified as core features using the Feature-table plugin . Aer ft wards, the 10 most-abundant OTUs for fungi and oomycetes were selected among them and classified as the core biome (Supplementary Information Table S10). The relative abundance of the top 10 core features with respect to the number of paired sequences was calculated separately for each dataset Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 12 www.nature.com/scientificreports www.nature.com/scientificreports/ (oomycetes and fungi). Rarefaction curves (Supplementary Information, Fig. S1) were evaluated and Good’s Coverage was calculated to estimate the adequate depth of sampling in the abundance matrices, the abundance matrices being rarefied with an even sampling depth for each dataset, with 10 iterations to avoid bias due to uneven sequencing depth. The Shannon (H’) diversity index, OTUs diversity and Evenness vector indices were calculated to study the diversity of the fungal and oomycete communities in each sample and the Bray-Curtis dis- similarity matrix was calculated; the β-diversity was analyzed to find the differences in community composition among the study zones. statistical methods. All the statistical analyses were performed in the environment R Studio V 1.0.143 (RStudio Inc. 2015. Boston, MA, USA) running under R 3.4.0 (R Foundation for Statistical Computing, 2014. Vienna, Austria). Prior to the analysis of differences of mean defoliation between zones, and for taxon abundance among defo- liation level or study zones, the normality of variable residues was checked through the normal Q-Q graph and Shapiro-Wilk test. If the data did not satisfy the assumptions of homoscedasticity and normality, the logarithmic transformation was used and then, normality checked again. For normally-distributed variables, the differences were compared using ANOVA followed by a post hoc Bonferroni mean differences test. The Kruskal-Wallis rank sum test with the Bonferroni correction was used instead ANOVA when variables did not show normal distribu- tion. Student-t Test was used for comparisons between defoliation levels. A correlation matrix was calculated to analyse the relationships between those functional guilds present in five or more samples of all zones, thus excluding the Lichen Parasite and Endophyte guilds. Unclassified features were excluded from the correlation matrix to avoid artifacts in the analysis. The correlation analysis was used also to evaluate the relationships between the key taxa resulting from the core biome analysis. The Pearson or Spearman correlation test was used, depending on whether the two variables in each pair were distributed normally or not. Null hypotheses were rejected in all cases when p ≤ 0.05. Alpha-diversity indices were subjected to the Kruskall-Wallis rank sum test with the Bonferroni correction, to compare the diversity levels of fungal and oomycete communities in the distinct study zones and according to the mean defoliation class of the plots. The Bray-Curtis dissimilarity matrix was used in non-metric mul- tidimensional scaling analysis (NMDS, vegan package , together with Analysis of Similarity (ANOSIM, 999 permutations), to determine the influence of the study location on the composition of the fungal and oomycete communities. Other packages used for data analysis and representation were dunn.test, nortest, ggplot2, ggpurb, and devtools (https://CRAN.R-project.org/). References 1. Prospero, S. & Cleary, M. Effects of Host Variability on the Spread of Invasive Forest. Diseases. Forests 8, 80 (2017). 2. Liebhold, A. M. et al. Biological invasions in forest ecosystems. Biol. Invasions 19, 3437–3458 (2017). 3. Sena, K., Crocker, E., Vincelli, P. & Barton, C. Phytophthora cinnamomi as a driver of forest change: Implications for conservation and management. For. Ecol. Manag. 409, 799–807 (2018). 4. Rizzo, D. M. & Garbelotto, M. Sudden oak death: endangering California and Oregon forest ecosystems. Front. Ecol. Environ. 1, 197–204 (2003). 5. Burgess, T. I. et al. Current and projected global distribution of Phytophthora cinnamomi, one of the world’s worst plant pathogens. Glob. Change Biol. 23, 1661–1674 (2017). 6. Vettraino, A. M. et al. Occurrence and distribution of Phytophthora species in European chestnut stands, and their association with Ink Disease and crown decline. Eur. J. Plant Pathol. 111, 169 (2005). 7. Thomas, F. M., Blank, R. & Hartmann, G. Abiotic and biotic factors and their interactions as causes of oak decline in Central. Europe. For. Pathol. 32, 277–307 (2002). 8. Jung, T., Blaschke, H. & Osswald, W. Involvement of soilborne Phytophthora species in Central European oak decline and the effect of site factors on the disease. Plant Pathol. 49, 706–718 (2000). 9. Vettraino, A. M. et al. Occurrence of Phytophthora species in oak stands in Italy and their association with declining oak trees. For. Pathol. 32, 19–28 (2002). 10. Brasier, C. M. Phytophthora cinnamomi and oak decline in southern Europe. Environmental constraints including climate change. Ann. Sci. For. 53, 347–358 (1996). 11. Navarro-Cerrillo, R. M. et al. Los procesos de decaimiento de encinas y alcornoques. (Dirección General de Gestión del Medio Natural. Consjería de Medio Ambiente. Junta de Andalucía., 2004). 12. Díaz Esteban, M. & Pulido Díaz, F. 6310: Dehesas perennifolias de Quercus spp. in Bases ecológicas preliminares para la conservación de los tipos de hábitat de interés comunitario en España (Ministerio de Medio Ambiente y Medio Rural y Marino, 2009). 13. Marañón, T. Agro-Sylvo-Pastoral Systems in the Iberian Peninsula: Dehesas and Montados. Rangelands 10, 255–258 (1988). 14. Pulido, F. & Picardo, A. Libro Verde de la Dehesa. Documento para el debate hacia una Estrategia Ibérica de gestión. (Consejería de Medio Ambiente, Junta de Castilla y León, 2010). 15. Horta, M., Caetano, P., Medeira, C., Maia, I. & Cravador, A. Involvement of the β-cinnamomin elicitin in infection and colonisation of cork oak roots by Phytophthora cinnamomi. Eur. J. Plant Pathol. 127, 427–436 (2010). 16. Brasier, C. M., Robredo, F. & Ferraz, J. F. P. Evidence for Phytophthora cinnamomi involvement in Iberian oak decline. Plant Pathol. 42, 140–145 (1993). 17. Oßwald, W. et al. Strategies of attack and defence in woody plant–Phytophthora interactions. For. Pathol. 44, 169–190 (2014). 18. Gea-Izquierdo, G., Fernández-de-Uña, L. & Cañellas, I. Growth projections reveal local vulnerability of Mediterranean oaks with rising temperatures. For. Ecol. Manag. 305, 282–293 (2013). 19. Corcobado, T., Cubera, E., Juárez, E., Moreno, G. & Solla, A. Drought events determine performance of Quercus ilex seedlings and increase their susceptibility to Phytophthora cinnamomi. Agric. For. Meteorol. 192–193, 1–8 (2014). 20. Ruiz Gómez, F., Pérez-de-Luque, A., Sánchez-Cuesta, R., Quero, J. & Navarro Cerrillo, R. Differences in the response to acute drought and Phytophthora cinnamomi Rands Infection in Quercus ilex L. seedlings. Forests 9, 634 (2018). 21. Ruiz Gómez, F. J., Navarro-Cerrillo, R. M., Sánchez-Cuesta, R. & Pérez-de-Luque, A. Histopathology of infection and colonization of Quercus ilex fine roots by Phytophthora cinnamomi . Plant Pathol. 64, 605–616 (2015). 22. Corcobado, T., Moreno, G., Azul, A. M. & Solla, A. Seasonal variations of ectomycorrhizal communities in declining Quercus ilex forests: interactions with topography, tree health status and Phytophthora cinnamomi infections. Forestry 88, 257–266 (2015). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 13 www.nature.com/scientificreports www.nature.com/scientificreports/ 23. Duque-Lazo, J., Navarro-Cerrillo, R. M., van Gils, H. & Groen, T. A. Forecasting oak decline caused by Phytophthora cinnamomi in Andalusia: Identification of priority areas for intervention. For. Ecol. Manag. 417, 122–136 (2018). 24. Serrano, M. S., Fernández-Rebollo, P., De Vita, P. & Sánchez, M. E. Calcium mineral nutrition increases the tolerance of Quercus ilex to Phytophthora root disease ae ff cting oak rangeland ecosystems in Spain. Agrofor. Syst. 87, 173–179 (2013). 25. Navarro Cerrillo, R. M., Terán Bocero, A. I. & Sánchez, M. E. Acción preventiva y curativa del fosfonato en el control de Phytophthora cinnamomi Rands en encina y alcornoque. 32, 685–694 (2006). 26. Ríos, P. et al. Brassica -based seedmeal biofumigation to control Phytophthora cinnamomi in the Spanish “dehesa” oak trees. Phytopathol. Mediterr. 56, 392–399 (2018). 27. Gómez-Aparicio, L. et al. Spatial patterns of soil pathogens in declining Mediterranean forests: implications for tree species regeneration. New Phytol. 194, 1014–1024 (2012). 28. Corcobado, T., Cubera, E., Pérez-Sierra, A., Jung, T. & Solla, A. First report of Phytophthora gonapodyides involved in the decline of Quercus ilex in xeric conditions in Spain. New Dis. Rep. 22, 33 (2010). 29. Català, S., Berbegal, M., Pérez-Sierra, A. & Abad-Campos, P. Metabarcoding and development of new real-time specific assays reveal Phytophthora species diversity in holm oak forests in eastern Spain. Plant Pathol. 66, 115–123 (2017). 30. Mora-Sala, B., Berbegal, M. & Abad-Campos, P. The use of qPCR reveals a high frequency of Phytophthora quercina in two Spanish holm oak areas. Forests 9, 697 (2018). 31. Gdanetz, K. & Trail, F. The wheat microbiome under four management strategies, and potential for endophytes in disease protection. Phytobiomes 1, 158–168 (2017). 32. Deveau, A. et al. Bacterial–fungal interactions: ecology, mechanisms and challenges. FEMS Microbiol. Rev. (2018). https://doi. org/10.1093/femsre/fuy008. 33. Pérez-de-Luque, A. et al. The interactive effects of arbuscular mycorrhiza and plant growth-promoting rhizobacteria synergistically enhance host plant defences against pathogens. Sci. Rep. 7 (2017). 34. Aleandri, M. P. et al. Use of nursery potting mixes amended with local Trichoderma strains with multiple complementary mechanisms to control soil-borne diseases. Crop Prot. 67, 269–278 (2015). 35. Blom, J. M., Vannini, A., Vettraino, A. M., Hale, M. D. & Godbold, D. L. Ectomycorrhizal community structure in a healthy and a Phytophthora-infected chestnut (Castanea sativa Mill.) stand in central Italy. Mycorrhiza 20, 25–38 (2009). 36. Lancellotti, E. & Franceschini, A. Studies on the ectomycorrhizal community in a declining Quercus suber L. stand. Mycorrhiza 23, 533–542 (2013). 37. Malcolm, G. M., Kuldau, G. A., Gugino, B. K. & Jiménez-Gasco, M. D. M. Hidden host plant associations of soilborne fungal pathogens: an ecological perspective. Phytopathology 103, 538–544 (2013). 38. Lee, S.-J. & Rose, J. K. C. Mediation of the transition from biotrophy to necrotrophy in hemibiotrophic plant pathogens by secreted effector proteins. Plant Signal. Behav. 5, 769–772 (2010). 39. Fawke, S., Doumane, M. & Schornack, S. Oomycete interactions with plants: infection strategies and resistance principles. Microbiol. Mol. Biol. Rev. 79, 263–280 (2015). 40. Vannini, A., Bruni, N., Tomassini, A., Franceschini, S. & Vettraino, A. M. Pyrosequencing of environmental soil samples reveals biodiversity of the Phytophthora resident community in chestnut forests. FEMS Microbiol. Ecol. 85, 433–442 (2013). 41. Sapkota, R. & Nicolaisen, M. Cropping history shapes fungal, oomycete and nematode communities in arable soils and affects cavity spot in carrot. Agric. Ecosyst. Environ. 257, 120–131 (2018). 42. Martin, F. N., Abad, Z. G., Balci, Y. & Ivors, K. Identic fi ation and detection of Phytophthora: reviewing our progress, identifying our needs. Plant Dis. 96, 1080–1103 (2012). 43. Hawksworth, D. L. Global species numbers of fungi: are tropical studies and molecular approaches contributing to a more robust estimate? Biodivers. Conserv. 21, 2425–2433 (2012). 44. Tedersoo, L. et al. Global diversity and geography of soil fungi. Science 346, 1256688–1256688 (2014). 45. Veach, A. M., Stokes, C. E., Knoepp, J., Jumpponen, A. & Baird, R. Fungal communities and functional guilds shift along an elevational gradient in the Southern Appalachian Mountains. Microb. Ecol. 1–13, https://doi.org/10.1007/s00248-017-1116-6 (2017). 46. Chen, W. et al. Plant diversity is coupled with beta not alpha diversity of soil fungal communities following N enrichment in a semi-arid grassland. Soil Biol. Biochem. 116, 388–398 (2018). 47. Abdelfattah, A., Malacrinò, A., Wisniewski, M., Cacciola, S. O. & Schena, L. Metabarcoding: A powerful tool to investigate microbial communities and shape future plant protection strategies. Biol. Control 120, 1–10 (2018). 48. Coleman-Derr, D. et al. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 209, 798–811 (2016). 49. Berendsen, R. L., Pieterse, C. M. J. & Bakker, P. A. H. M. The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486 (2012). 50. Lledó, S., Santamaría, O., Rodrigo, S. & Poblaciones, M. J. Endophytic mycobiota associated with Trifolium subterraneum growing under semiarid conditions. Ann. Appl. Biol. 168, 243–254 (2016). 51. Kwaśna, H., Szewczyk, W. & Behnke-Borowczyk, J. Fungal root endophytes of Quercus robur subjected to flooding. For. Pathol. 46, 35–46 (2016). 52. Nordberg, H. et al. The genome portal of the Department of Energy Joint Genome Institute: 2014 updates. Nucleic Acids Res. 42, D26–31 (2014). 53. Tamayo-Vélez, Á. & Osorio, N. W. Soil fertility improvement by litter decomposition and inoculation with the fungus Mortierella sp. in avocado plantations of Colombia. Commun. Soil Sci. Plant Anal. 49, 139–147 (2018). 54. Li, F. et al. Application of environmental DNA metabarcoding for predicting anthropogenic pollution in rivers. Environ. Sci. Technol. 52, 11708–11719 (2018). 55. Mosca, E., Montecchio, L., Barion, G., Dal Cortivo, C. & Vamerali, T. Combined effects of thinning and decline on fine root dynamics in a Quercus robur L. forest adjoining the Italian Pre-Alps. Ann. Bot. 119, 1235–1246 (2017). 56. Dickie, I. A., Guza, R. C., Krazewski, S. E. & Reich, P. B. Shared ectomycorrhizal fungi between a herbaceous perennial (Helianthemum bicknellii) and oak (Quercus) seedlings. New Phytol. 164, 375–382 (2004). 57. Corcobado, T., Vivas, M., Moreno, G. & Solla, A. Ectomycorrhizal symbiosis in declining and non-declining Quercus ilex trees infected with or free of Phytophthora cinnamomi. For. Ecol. Manag. 324, 72–80 (2014). 58. Aoki, T., O’Donnell, K. & Geiser, D. M. Systematics of key phytopathogenic Fusarium species: Current status and future challenges. J. Gen. Plant Pathol. 80, 189–201 (2014). 59. Gallego, F. J., de Algaba, A. P. & Fernandez-Escobar, R. Etiology of oak decline in Spain. Eur. J. For. Pathol. 29, 17–27 (1999). 60. Ragazzi, A., Mugnai, L., Moricca, S., Vagniluca, S. & Dellavalle, I. Requirements and biological aspects of Fusarium eumartii and its possible role in oak decline in North-eastern Italian oak forests. For. Pathol. 23, 171–177 (1993). 61. Ragazzi, A., Moricca, S., Dellavalle, I. & Turco, E. Fusarium eumartii growth in resistant and susceptible oak species. J. Phytopathol. 149, 503–507 (2001). 62. Lynch, S. C., Zambino, P. J., Mayorquin, J. S., Wang, D. H. & Eskalen, A. Identification of new fungal pathogens of coast live oak in California. Plant Dis. 97, 1025–1036 (2013). 63. Wit, M., Sierota, Z., Oszako, T., Mirzwa-Mroz, E. & Wakulinski, W. Fusarium spp. on the above-ground organs of dying oaks - a new threat? Sylwan 159, 403–410 (2015). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 14 www.nature.com/scientificreports www.nature.com/scientificreports/ 64. Moricca, S. et al. Biotic factors ae ff cting ecosystem services in urban and peri-urban forests in Italy: The role of introduced and impending pathogens and pests. Forests 9, 65 (2018). 65. Pecoraro, L. et al. Macrofungi in Mediterranean maquis along seashore and altitudinal transects. Plant Biosyst. - Int. J. Deal. Asp. Plant Biol. 148, 367–376 (2014). 66. Kwaśna, H. & Szewczyk, W. Effects of fungi isolated from Quercus robur roots on growth of oak seedlings. Dendrobiology 75, 99–112 (2016). 67. Omari, B. E., Fleck, I., Aranda, X. & Moret, A. & MartíNadal. Effect of fungal infection on leaf gas-exchange and chlorophyll uo fl rescence in Quercus ilex. Ann. For. Sci. 58, 165–174 (2001). 68. Sánchez, M., Caetano, P., Ferraz, J. & Trapero, A. Phytophthora disease of Quercus ilex in south-western Spain. For. Pathol. 32, 5–18 (2002). 69. Serrano, M. D. et al. Phytophthora cinnamomi and Pythium spiculum as main agents of Quercus decline in southern Spain and Portugal. In IOBC/WPRS Bulletin 76, 97–100 (2012). 70. Moralejo, E., García-Muñoz, J. A. & Descals, E. Susceptibility of Iberian trees to Phytophthora ramorum and P. cinnamomi. Plant Pathol. 58, 271–283 (2009). 71. McConnell, M. E. & Balci, Y. Phytophthora cinnamomi as a contributor to white oak decline in mid-atlantic United States forests. Plant Dis. 98, 319–327 (2013). 72. Jung, T. & Burgess, T. I. Re-evaluation of Phytophthora citricola isolates from multiple woody hosts in Europe and North America reveals a new species, Phytophthora plurivora sp. nov. Persoonia 22, 95–110 (2009). 73. Jung, T., Blaschke, H. & Neumann, P. Isolation, identification and pathogenicity of Phytophthora species from declining oak stands. Eur. J. For. Pathol. 26, 253–272 (1996). 74. Pérez-Sierra, A. et al. Previously unrecorded low-temperature Phytophthora species associated with Quercus decline in a Mediterranean forest in eastern Spain. For. Pathol. 43, 331–339 (2013). 75. Coll, L., Camarero, J. J. & Aragón, J. M. D. Fine root seasonal dynamics, plasticity, and mycorrhization in 2 coexisting Mediterranean oaks with contrasting aboveground phenology. Écoscience 19, 238–245 (2012). 76. Richard, F., Moreau, P.-A., Selosse, M.-A. & Gardes, M. Diversity and fruiting patterns of ectomycorrhizal and saprobic fungi in an old-growth Mediterranean forest dominated by Quercus ilex L. Can. J. Bot. 82, 1711–1729 (2004). 77. Dawkins, K. & Esiobu, N. Arbuscular and ectomycorrhizal fungi associated with the invasive Brazilian pepper tree (Schinus terebinthifolius) and two native plants in south Florida. Front. Microbiol. 8 (2017). 78. Jimu, L., Kemler, M., Mujuru, L. & Mwenje, E. Illumina DNA metabarcoding of Eucalyptus plantation soil reveals the presence of mycorrhizal and pathogenic fungi. For. Int. J. For. Res. 91, 238–245 (2018). 79. Linaldeddu, B. T., Sirca, C., Spano, D. & Franceschini, A. Variation of endophytic cork oak-associated fungal communities in relation to plant health and water stress. For. Pathol. 41, 193–201 (2011). 80. Vinale, F. et al. Trichoderma–plant–pathogen interactions. Soil Biol. Biochem. 40, 1–10 (2008). 81. Widmer, T. L. Screening Trichoderma species for biological control activity against Phytophthora ramorum in soil. Biol. Control 79, 43–48 (2014). 82. García-Núñez, H. G. et al. Caracterización morfológica y molecular de cepas nativas de Trichoderma y su potencial de biocontrol sobre Phytophthora infestans. Rev. Mex. Fitopatol. 35, 58–79 (2017). 83. Chemeltorit, P. P., Mutaqin, K. H. & Widodo, W. Combining Trichoderma hamatum THSW13 and Pseudomonas aeruginosa BJ10–86: a synergistic chili pepper seed treatment for Phytophthora capsici infested soil. Eur. J. Plant Pathol. 147, 157–166 (2017). 84. Promwee, A., Yenjit, P., Issarakraisila, M., Intana, W. & Chamswarng, C. Efficacy of indigenous Trichoderma harzianum in controlling Phytophthora leaf fall (Phytophthora palmivora) in Thai rubber trees. J. Plant Dis. Prot. 124, 41–50 (2017). 85. Widmer, T. L. et al. Remediation of Phytophthora ramorum-infested soil with Trichoderma asperellum isolate 04–22 under ornamental nursery conditions. Biol. Control 118, 67–73 (2018). 86. Bae, H. et al. The beneficial endophyte Trichoderma hamatum isolate DIS 219b promotes growth and delays the onset of the drought response in Theobroma cacao . J. Exp. Bot. 60, 3279–3295 (2009). 87. Carrero‐Carrón, I. et al. Interactions between Trichoderma harzianum and defoliating Verticillium dahliae in resistant and susceptible wild olive clones. Plant Pathol. 67, 1758–1767 (2018). 88. Menkis, A. et al. Occurrence and impact of the root-rot biocontrol agent Phlebiopsis gigantea on soil fungal communities in Picea abies forests of northern Europe. FEMS Microbiol. Ecol. 81, 438–445 (2012). 89. Sylla, J. et al. Leaf microbiota of strawberries as ae ff cted by biological control agents. Phytopathology 103, 1001–1011 (2013). 90. Pinto, C. et al. Unravelling the diversity of grapevine microbiome. PLOS ONE 9, e85622 (2014). 91. Eichhorn, J. et al. Part IV: Visual assessment of crown condition and damaging agents. In UNECE ICP Forests Programme Co- ordinating Centre. Manual on methods and criteria for harmonized sampling, assessment, monitoring and analysis of the effects of air pollution on forests. 54 (Thünen Institute of Forest Ecosystems, 2017). 92. Gao, Z., Perez-Perez, G. I., Chen, Y. & Blaser, M. J. Quantitation of major human cutaneous bacterial and fungal populations. J. Clin. Microbiol. 48, 3575–3581 (2010). 93. Cooke, D. E. L., Drenth, A., Duncan, J. M., Wagels, G. & Brasier, C. M. A molecular phylogeny of Phytophthora and related oomycetes. Fungal Genet. Biol. 30, 17–32 (2000). 94. Bálint, M., Schmidt, P.-A., Sharma, R., Thines, M. & Schmitt, I. An Illumina metabarcoding pipeline for fungi. Ecol. Evol. 4, 2642–2653 (2014). 95. Sapkota, R. & Nicolaisen, M. An improved high throughput sequencing method for studying oomycete communities. J. Microbiol. Methods 110, 33–39 (2015). 96. Nilsson, R. H. et al. An open source sowa ft re package for automated extraction of ITS1 and ITS2 from fungal ITS sequences for use in high-throughput community assays and molecular ecology. Fungal Ecol. 3, 284–287 (2010). 97. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010). 98. Huson, D. H. et al. MEGAN Community Edition - Interactive exploration and analysis of large-scale microbiome sequencing data. PLOS Comput. Biol. 12, e1004957 (2016). 99. Nguyen, N. H. et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248 (2016). 100. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010). 101. Chen, J. et al. Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinforma. Oxf. Engl. 28, 2106–2113 (2012). 102. McDonald, D. et al. e Th Biological Observation Matrix (BIOM) format or: how I learned to stop worrying and love the ome-ome. GigaScience 1, 7 (2012). 103. Oksanen, J. et al. Vegan: Community Ecology Package. (2018). 104. Costa Pérez, J. C., Martín Vicente, Á., Fernández Alés, R. & Estirado Oliet, M. Dehesas de Andalucía. Caracterización ambiental (2006). (Junta de Andalucía. Consejería de Medio Ambiente, 2006). Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 15 www.nature.com/scientificreports www.nature.com/scientificreports/ Acknowledgements This work was conducted with support from the projects ESPECTRAMED (CGL2017-86161-R) and ENCINOMICA (BIO2015-64737-R), funded by the Spanish Ministry of Economy, Industry and Competitivity. e fir Th st author was a recipient of an FPU fellowship, and the work was carried out during a research stay in the TUM School of Life Sciences Weihensthephan (Freising, Germany), funded by the Spanish Ministry of Education (FPU fellowship, FPU13/00231). The authors thank the “Consejería de Medio Ambiente, Junta de Andalucía” for their support in sample collection, Rafael Sánchez de la Cuesta and Roberto Jesús Cabrera Puerto by their support in the sample processing and the staff of the Dendrochronology, Silviculture, and Climate Change Laboratory at Cordoba University, for their assistance during this research. We thank Dr. David Walker for revisions of the written English in different versions of this manuscript. Also, authors thank the support of the “Campus de Excelencia en Agroalimentación - CeiA3”. Author Contributions F.J.R.G., R.M.N.C., A.P.L. and A.V. designed the work idea and established the main hypothesis, F.J.R.G., R.M.N.C. and A.P.L. performed the sample collection campaign; F.J.R.G., A.V. and C.M.R. designed the sample treatment and the laboratory workflow, F.J.R.G. and C.M.R. performed the laboratory work and the data analysis, under the supervision of W.O.; F.J.R.G. and R.M.N.C. performed the data interpretation and wrote the manuscript with the help of all the rest of co-authors. Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-41804-y. Competing Interests: The authors declare no competing interests. Publisher’s note: 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 Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted 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 license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2019 Scientific Repo R ts | (2019) 9:5315 | https://doi.org/10.1038/s41598-019-41804-y 16

Journal

Scientific ReportsSpringer Journals

Published: Mar 29, 2019

There are no references for this article.