An Overview of Phytophthora Species Inhabiting Declining Quercus suber Stands in Sardinia (Italy)

: Cork oak forests are of immense importance in terms of economic, cultural, and ecological value in the Mediterranean regions. Since the beginning of the 20th century, these forests ecosystems have been threatened by several factors, including human intervention, climate change, wildﬁres, pathogens, and pests. Several studies have demonstrated the primary role of the oomycete Phytophthora cinnamomi Ronds in the widespread decline of cork oaks in Portugal, Spain, southern France, and Italy, although other congeneric species have also been occasionally associated. Between 2015 and 2019, independent surveys were undertaken to determine the diversity of Phytophthora species in declining cork oak stands in Sardinia (Italy). Rhizosphere soil samples were collected from 39 declining cork oak stands and baited in the laboratory with oak leaﬂets. In addition, the occurrence of Phytophthora was assayed using an in-situ baiting technique in rivers and streams located throughout ten of the surveyed oak stands. Isolates were identiﬁed by means of both morphological characters and sequence analysis of internal transcribed spacer (ITS) regions of ribosomal DNA. In total, 14 di ﬀ erent Phytophthora species were detected. Phytophthora cinnamomi was the most frequently isolated species from rhizosphere soil, followed by P. quercina , P. pseudocryptogea , and P. tyrrhenica . In contrast, P. gonapodyides turned out to be the most dominant species in stream water, followed by P. bilorbang , P. pseudocryptogea , P. lacustris , and P. plurivora . Pathogenicity of the most common Phytophthora species detected was tested using both soil infestation and log inoculation methods. This study showed the high diversity of Phytophthora species inhabiting soil and watercourses, including several previously unrecorded species potentially involved in the decline of cork oak forests.


Introduction
Cork oak (Quercus suber L.) represents an important component of the Mediterranean forests landscape, covering more than 2 million ha across southern European and northern African countries [1]. This type of forest ecosystem has great socio-economic value, providing a range of non-timber forest products, such as cork, firewood, grazing, honey, and mushrooms, playing a key role in the rural economy in less favorable regions [2,3]. In particular, cork production represents a highly sustainable non-wood product derived from forests in the western Mediterranean countries and an additional source of income for farmers. In Sardinia (Italy), which hosts more than 80% of the Italian distribution area, and checked daily for Phytophthora-like hyphae development, which was subsequently transferred to Petri dishes containing carrot-agar (CA; 16 g agar technical no.3, Oxoid Ltd., Basingstoke, UK, 200 g carrots and 1000 mL distilled water) [29].

Stream Baiting and Phytophthora Isolation
In addition to soil samples, river and stream water within ten surveyed cork oak forests ( Figure 1) were assayed using the in-situ baiting technique in spring 2018. Watercourses were chosen depending on their water flow and on the capacity to collect water from the bordering forests. They were subdivided into two different groups, mainly based on the altimetric gradient, including valley floor rivers with permanent water flow, and mountain or hill streams or water catchments into the forests, with water flow strictly correlated with the seasonal rainfall and very often drying up in summer. River baiting was made using an adaptation of the method described by Reeser et al. [30] and Hüberli et al. [31], which consisted of two squared layers of fly mesh or metallic net, sealed together, with young leaves of different plant species, such as Alnus glutinosa Mill., Arbutus unedo L., Buxus sempervirens L., Laurus nobilis L., Fraxinus ornus L., Hedera helix L., Parthenocissus quinquefolia L., Q. ilex L., Q. suber L., and Taxus baccata L. placed between the two mesh layers. Rafts were placed in-situ at the same time as their set up, and each raft contained about 20 leaves. Cork stoppers were used to float the raft over the water surface. Three rafts were placed every 200 m from each other along the watercourse. The rafts were fastened to natural restraints, such as branches and rocks, and left floating over the water for 3-4 days and then brought to the laboratory where the leaves were washed with sterile water and blotted dry on filter paper. Phytophthora isolations were made by placing small fragments cut from necrotic lesions detected on the leaf baits on SMA. Any developing colonies were sub-cultured on CA for further analyses.

Morphological and Molecular Identification
The isolates obtained from soil samples and stream water were first grouped based on their colony growth patterns after 5-7 days at 20 • C in the dark on CA. In addition, morphological features of sporangia, oogonia, antheridia, chlamydospores, hyphal swellings, and aggregations were examined under the Leitz Diaplan compound microscope (Leitz, Wetzlar, Germany) and compared with species descriptions in the literature [21,[32][33][34][35][36][37][38]. A subset of representative morphotypes (88 isolates) was selected for molecular analyses, which consisted of DNA extraction, amplification, sequencing, and analysis of sequences of the entire region of the internal transcribed spacers (ITS1 and ITS2) and the 5.8 S gene of the rDNA. DNA was extracted from mycelium fragments, using the extraction kit InstaGene™ Matrix (BioRad Laboratories, Hercules, CA, USA). The amplification of the ITS region was carried out with a thermocycler (Hybaid PCR Express), using the forward primers ITS1 or ITS6 and the reverse primer ITS4 [39,40]. A total volume of 50 µL, consisting of 18.2 µL of water, 5 µL of BSA, 5 µL of dNTPs, 5 µL of both ITS6 and ITS4 primers, 10 µL of the buffer, 0.3 µL of Go Taq polymerase, and 1.5 µL of DNA from each morphotype, was used for standard PCR (polymerase chain reaction). The cycle used for the amplification of the ITS genes regions was as follows: initial denaturation of 1 min at 94 • C, followed by 35 cycles of 1 min at 94 • C, 1 min at 56 • C, 1 min at 72 • C, followed by a final elongation phase of 7 min at 72 • C. The PCR products were purified using the EUROGOLD gel extraction kit (EuroClone S.p.A., Pero, Italy). After quantification, purified PCR amplicons and the sequencing primers were sent to BMR Genomics sequencing service (https://www.bmr-genomics.it). DNA sequence chromatograms were viewed and edited using BioEdit v. 5.0.6 software [41]. Heterozygous sites observed were labeled according to the IUPAC coding system. Isolates were assigned to a species when sequence identities were above a 99% cut-off with respect to those of ex-type isolates or key isolates. All sequences were deposited at GenBank (http://www.ncbi.nlm.nih.gov/), and the accession numbers are given in Table 1.

Pathogenicity Tests
Pathogenicity of the most frequent Phytophthora species isolated was assayed using the soil infestation method described by Jung et al. [27], with some modification as reported by Scanu et al. [29]. In particular, two isolates of P. cinnamomi, P. gonapodyides, P. pseudocryptogea, P. psychrophila, P. quercina, and P. tyrrhenica were grown in individual 500 mL Erlenmeyer flasks containing an autoclaved mixture of 250 mL of vermiculite and 150 mL of Lolium italicum seeds thoroughly moistened with 100 mL of carrot juice (200 mL/L carrot juice, 3 g/L CaCO 3 , and 800 mL/L distilled water). Flasks were incubated at 20 • C for 1 month, then 20 mL of inoculum was collected and inserted inside the soil of 2-year-old cork oak seedlings (provided by the Regional Agency Fo.Re.S.T.A.S.). The substrate in the controls received a sterile mixture of vermiculite/seeds-carrot juice at the same ratio. To stimulate the production of sporangia and pathogen spread and infection via zoospores, pots were flooded immediately after inoculation for 48 h, and flooding was repeated at three-week intervals by immerging pots in 10 L buckets just to 1 cm above the soil surface. There were eight replicates per isolate and controls. After 5 months of incubation at 20 • C (±2 • C), 70% relative humidity with a 12/12 h photoperiod, seedlings were visually assessed for symptoms, and the mortality rate was recorded; then each plant was removed from the pot, and the root system gently washed under tap water. Single roots were cut off at the collar, and after scanning, the total root length of all the plant root systems was measured using the APS Assess 2.0 software (The American Phytopathological Society, St. Paul, MN, USA). The remaining soil was baited following the method described above to determine whether the pathogen was still viable. Re-isolations were also made directly from necrotic roots using SMA. Pathogenicity of the above isolates (except for those of P. tyrrhenica) was further tested using freshly cut logs of cork oak following the method described by Brasier and Kirk [42]. Four logs (1.4 m long and 20 cm in diameter) were cut from stems of living cork oak trees 24 h before the experiment, and the cut ends were sealed with a liquid waterproofing membrane. In each log, three bands were marked around the log circumference, 30 cm apart from each other, with 5 inoculation points per band, about 15 cm apart. After sterilizing the bark with 70% ethanol, a 7 mm diameter hole was punched through the bark to the wood surface with a steel cork borer. The same-sized plug was taken from the edge of an actively growing colony on CA and used as inoculum by inserting into the hole replacing the bark plug. Three control inoculation points per log were inoculated with a sterile CA plug and covered with the removed piece of bark. Moist cotton wool was placed over the wounds, covered with a 5 × 5 cm piece of aluminum foil, and sealed with an adhesive PVC tape. There were four replicates per isolate. Inoculated logs were covered individually in loose polythene sleeves (sealed at both ends) and incubated at 20 • C (±2 • C) in an air-conditioned laboratory and checked weekly for the appearance of symptoms. After 45 days, the experiment was finished, and logs were destructively sampled by removing the periderm with a drawknife to expose the phloem. Each lesion's outline was then recorded on tracing paper and scanned on an Epson Perfection V30 photo scanner, and the lesion area calculated using APS Assess software, as described by Scanu and Webber [43]. Re-isolation of all the inoculated Phytophthora species onto SMA was attempted from the lesion margins.
Statistical analyses for both pathogenicity tests were performed using XLSTAT software (Addinsoft). Data were first checked for normality and then subjected to analysis of variance (ANOVA). Statistical differences among mean values of root lengths and lesion areas were determined using Fisher's protected least significant difference (LSD) test. Differences with p < 0.05 were considered significant.

Symptomatology
A wide range of symptoms of decline was observed on cork oak trees across all the investigated sites. These included rapid dieback of the crown in both mature ( Figure 2a) and young oak trees (Figure 2b), which was frequently observed in early autumn, especially after a long summer and drought conditions. In the case of afforestation sites (10 to 20-year-old), the infections could reach epidemic levels and cause extensive mortality of oak trees. Other symptoms included shoot dieback and increased transparency of the whole crown, leaf chlorosis, and abundant proliferation of epicormic shoots on stems and branches (Figure 2c,d). At the collar level, trees showed necrotic bark lesions frequently associated with black exudation and very often girdling the stem (Figure 2e). In the root system of declining oak trees, an extensive loss of both lateral small woody roots and fine roots and callusing or open cankers on suberized roots were observed.

Soilborne Phytophthora Species
Phytophthora species were recovered from 68.5% of the 295 soil samples tested. The highest level of soil infestation was detected in the afforestation sites (80.4%), while in natural forests, the percentage of positive trees was 61.7%. In total, 224 isolates were obtained from rhizosphere soil samples collected from around symptomatic trees in declining cork oak stands (Table 1). All isolates conformed morphologically to previously known Phytophthora species. ITS sequence analysis of the isolates confirmed the morphological identification of all Phytophthora species. BLAST searches in GenBank showed 99-100% similarity with reference sequences, including those of ex-type cultures or representative isolates (Table 1). In total, eight Phytophthora species belonging to five (clade 3, 6, 7, 8, and 12) out of the twelve known phylogenetic clades were isolated, including P. cinnamomi, P. gonapodyides, P. pseudocryptogea, P. psychrophila, P. quercina, P. syringae, P. tyrrhenica, and P. ×cambivora (Table 1 and Figure 3a,c). Phytophthora cinnamomi from clade 7c was the most frequent species isolated from both natural forests and afforestation stands. It was detected from almost all investigated afforestation sites (from 30 out of 39 sites) with an infection rate of 80.2%, while its incidence was markedly lower in natural cork oak forests (45.7% of 122 investigated trees). At one afforestation site (QS35), P. cinnamomi was recovered from 27 out of 30 cork oak trees sampled (Table 1). It was the only species recovered in 10 investigated stands. Similarly, the second most common species, P. quercina, from clade 12 occurred in both afforestation and natural stands, with an infestation rate ranging from 8.1% to 30.2%, respectively. It was the only species isolated from declining trees in sites QS6, QS13, and QS25. All P. quercina isolates had identical ITS sequences; however, a certain phenotypic variation among the isolates was observed. Both P. pseudocryptogea (clade 8a) and P. tyrrhenica (clade 7a) were recovered from eight declining cork oak stands, with an infestation rate of

Soilborne Phytophthora Species
Phytophthora species were recovered from 68.5% of the 295 soil samples tested. The highest level of soil infestation was detected in the afforestation sites (80.4%), while in natural forests, the percentage of positive trees was 61.7%. In total, 224 isolates were obtained from rhizosphere soil samples collected from around symptomatic trees in declining cork oak stands (Table 1). All isolates conformed morphologically to previously known Phytophthora species. ITS sequence analysis of the isolates confirmed the morphological identification of all Phytophthora species. BLAST searches in GenBank showed 99-100% similarity with reference sequences, including those of ex-type cultures or representative isolates (Table 1). In total, eight Phytophthora species belonging to five (clade 3, 6, 7, 8, and 12) out of the twelve known phylogenetic clades were isolated, including P. cinnamomi, P. gonapodyides, P. pseudocryptogea, P. psychrophila, P. quercina, P. syringae, P. tyrrhenica, and P. ×cambivora (Table 1 and Figure 3a,c). Phytophthora cinnamomi from clade 7c was the most frequent species isolated from both natural forests and afforestation stands. It was detected from almost all investigated afforestation sites (from 30 out of 39 sites) with an infection rate of 80.2%, while its incidence was markedly lower in natural cork oak forests (45.7% of 122 investigated trees). At one afforestation site (QS35), P. cinnamomi was recovered from 27 out of 30 cork oak trees sampled (Table 1). It was the only species recovered in 10 investigated stands. Similarly, the second most common species, P. quercina, from clade 12 occurred in both afforestation and natural stands, with an infestation rate ranging from 8.1% to 30.2%, respectively. It was the only species isolated from declining trees in sites QS6, QS13, and QS25. All P. quercina isolates had identical ITS sequences; however, a certain phenotypic variation among the isolates was observed. Both P. pseudocryptogea (clade 8a) and P. tyrrhenica (clade 7a) were recovered from eight declining cork oak stands, with an infestation rate of around 10%. The ITS sequences of all P. pseudocryptogea isolates from rhizosphere soil matched the ex-holotype isolate (GenBank no. KP288376). However, they had a unique polymorphism at position 56 (C instead of Y) and were heterozygous at position 650 (Y instead of T). Phytophthora gonapodyides from clade 6b occurred only in natural contexts, isolated from 12 declining trees in seven sites, and always in association with other Phytophthora species. Among the less frequently isolated species, P. psychrophila (clade 3) and P. ×cambivora (clade 7a) were detected only from natural forests at a very low infestation rate (2.5% and 1%, respectively), while P. syringae was exclusively found at one afforestation site from two symptomatic trees. In 18 soil samples from declining trees and 12 sites, multiple Phytophthora species were detected. Phytophthora cinnamomi was isolated along with P. quercina in four samples and along with both P. quercina and P. tyrrhenica (two samples) or P. pseudocryptogea (one sample). In four cases, P. tyrrhenica was isolated together with P. cinnamomi and, in one case, with P. quercina.

Waterborne Phytophthora Species
In total, 115 Phytophthora isolates were detected in all watercourses monitored through ten selected declining cork oak stands. Based on morphological analyses and molecular identification,

Waterborne Phytophthora Species
In total, 115 Phytophthora isolates were detected in all watercourses monitored through ten selected declining cork oak stands. Based on morphological analyses and molecular identification, these isolates belonged to five phylogenetic clades (clade 2, 6, 7, 8, and 9) corresponding to nine formally known Phytophthora species, including Phytophthora bilorbang, Phytophthora chlamydospora, P. gonapodyides, P. hydropathica, P. lacustris, P. plurivora, P. pseudocryptogea, P. ×cambivora, and P. ×multiformis (Table 2 and Figure 3b,c). Overall, more than 50% of the isolates obtained were identified as P. gonapodyides, which was recovered from all watercourses surveyed. Interestingly, eight isolates detected from stands QS1, QS8, and QS9 were heterozygous at position 106 (R instead of G). Another isolate from stand QS1 (GON3) shared the heterozygous position 106 and had unique polymorphisms at positions 106 (A instead of R or G) and 517 (T instead of G), respectively. Two isolates, PH255 (QS10) and PH267 (QS9), also differed from the other P. gonapodyides isolates by having unique polymorphisms at positions 145 (T instead of C) and 517 (T instead of G), respectively.
Phytophthora bilorbang was the second most widespread species isolated from five investigated streams, four of which only flow seasonally. All isolates from stand QS21 were heterozygous at position 106 (Y instead of T). Phytophthora pseudocryptogea and P. plurivora were isolated only from permanent water bodies, while P. lacustris occurred from both permanent and intermittent watercourses ( Table 2). The ITS sequences of all P. lacustris isolates differed from the ex-type culture (GenBank no. AF266793) having a heterozygous site at position 783 (S instead of C). In addition, one isolate was heterozygous at position 458 (Y instead of C). Almost all isolates identified as P. pseudocryptogea from river water had identical ITS sequences than those isolates obtained from rhizosphere soil, differing from the ex-holotype by 2 bp. Moreover, one isolate (PH269) from stand Q12 differed from the ex-holotype isolate of P. pseudocryptogea and the other isolates obtained in this study by 4-5 bp at positions 56, 601, 650, 728, and 733. Phytophthora chlamydospora and P. ×cambivora were recovered from two streams, while P. ×multiformis and P. hydropathica were exclusively isolated from QS10 and QS13, respectively. The two P. chlamydospora isolates from stand QS9 were heterozygous at position 666 (Y instead of T). All isolates identified as P. hydropathica differed from the ex-type culture (GenBank no. EU583793) by having two heterozygous sites at positions 413 (S instead of C) and 665 (K instead of T) and by a unique polymorphism at position 628 (G instead of C). The ITS sequences of P. ×cambivora often generated overlapping ITS sequences starting at position 396 in both directions. This was a consistent pattern observed in all isolates of P. ×cambivora obtained from river water. The non-overlapping sequences up the indel positions were identical to that of the neotype culture of P. ×cambivora (GenBank no. KU899179).
Looking at the diversity of Phytophthora species across rivers, the two geographically close rivers at sites QS9 and QS10 hosted the highest number of Phytophthora species, followed by a stream in stand QS12 with five species detected. Only two species were isolated from water bodies at stands QS1-3, QS13, and QS21, and they were all from clade 6, except for a clade 9 species at site QS13 detected together with P. gonapodyides. a BIL = P. bilorbang, CAM = P. ×cambivora, CHL = P. chlamydospora, GON = P. gonapodyides, HYD = P. hydropathica, LAC = P. lacustris, MUL = P. ×multiformis, PLU = P. plurivora, PSE = P. pseudocryptogea.

Pathogenicity Test
The soil infestation experiment showed that all Phytophthora species tested were able to cause a significant reduction of the root system in 2-year-old cork oak seedlings (Figure 4a). The mean root length was significantly higher in control seedlings (p < 0.05) than in seedlings infected with Phytophthora isolates. Phytophthora cinnamomi was the most aggressive species causing a root length reduction near to 65% compared to the control seedlings, followed by P. pseudocryptogea (56.2%), P. tyrrhenica (44.6%), and P. quercina (40.2%). Phytophthora gonapodyides and P. psychrophila caused a root length reduction below 35%. Phytophthora cinnamomi was the only species associated with extensive lesions on the mother root, with lesions in some cases reaching 15 mm in length. Apart from P. psychrophila, all the other Phytophthora species were re-isolated from limited necrotic lesions on taproot. No symptoms of pathogen infection could be seen on the roots of control seedlings.
length was significantly higher in control seedlings (p < 0.05) than in seedlings infected with Phytophthora isolates. Phytophthora cinnamomi was the most aggressive species causing a root length reduction near to 65% compared to the control seedlings, followed by P. pseudocryptogea (56.2%), P. tyrrhenica (44.6%), and P. quercina (40.2%). Phytophthora gonapodyides and P. psychrophila caused a root length reduction below 35%. Phytophthora cinnamomi was the only species associated with extensive lesions on the mother root, with lesions in some cases reaching 15 mm in length. Apart from P. psychrophila, all the other Phytophthora species were re-isolated from limited necrotic lesions on taproot. No symptoms of pathogen infection could be seen on the roots of control seedlings.
In the log inoculation tests, lesions in the phloem tissue of cork oak caused by P. cinnamomi were significantly larger (p < 0.0001) compared with the negative controls, with a mean necrosis area of 80 cm 2 (Figure 4b). Phytophthora pseudocryptogea and P. gonapodyides also showed considerable aggressiveness on inoculated logs. The mean lesion area formed by P. cinnamomi was approximately two to three times larger (p < 0.0001) than that developed following inoculation with P. pseudocryptogea and P. gonapodyides, respectively (Figure 4b). Phytophthora psychrophila and P. quercina were not able to colonize phloem tissues producing lesions that did not differ significantly from the negative controls (p > 0.05). Apart from P. psychrophila and P. quercina, all Phytophthora species were readily re-isolated from the necrotic lesions. In contrast, the controls developed only limited discoloration around the inoculation point and never yielded any Phytophthora.  In the log inoculation tests, lesions in the phloem tissue of cork oak caused by P. cinnamomi were significantly larger (p < 0.0001) compared with the negative controls, with a mean necrosis area of 80 cm 2 (Figure 4b). Phytophthora pseudocryptogea and P. gonapodyides also showed considerable aggressiveness on inoculated logs. The mean lesion area formed by P. cinnamomi was approximately two to three times larger (p < 0.0001) than that developed following inoculation with P. pseudocryptogea and P. gonapodyides, respectively (Figure 4b). Phytophthora psychrophila and P. quercina were not able to colonize phloem tissues producing lesions that did not differ significantly from the negative controls (p > 0.05). Apart from P. psychrophila and P. quercina, all Phytophthora species were readily re-isolated from the necrotic lesions. In contrast, the controls developed only limited discoloration around the inoculation point and never yielded any Phytophthora.

Discussion
The extensive surveys made over four years across declining cork oak stands in Sardinia, together with morphological and ITS sequences analyses, have revealed the occurrence of 14 Phytophthora taxa from seven of the 12 known phylogenetic clades [21]. These include species common in Mediterranean oak soil, such as P. cinnamomi, P. gonapodyides, and P. quercina, and the less widespread species P. pseudocryptogea, P. psychrophila, P. syringae, P. tyrrhenica, and P. ×cambivora [12,17,20,44,45]. The detection of nine Phytophthora species from stream and river water represents the first attempt to look at the diversity of aquatic species in such forest ecosystems. Apart from P. gonapodyides and P. pseudocryptogea, all the other species identified from watercourses, P. bilorbang, P. chlamydospora, P. hydropathica, P. lacustris, P. plurivora, P. ×cambivora, and P. ×multiformis, were never reported in cork oak ecosystems.
Most of the previous surveys on Mediterranean oak decline have focused on the association of P. cinnamomi and Q. ilex (ssp. ballota and rotundifolia) [16,17,[46][47][48][49], while cork oak has been less studied as it appears to be more tolerant to the disease due to its defense response mechanisms to Phytophthora infection [14,15,50]. Our study represents the first extensive survey on the distribution of P. cinnamomi and other congeneric species in declining cork oak stands.
Among the Phytophthora species detected from soil samples, P. cinnamomi was the most common species encountered. Listed as one of the 100 worst invasive alien species, P. cinnamomi is considered one of the most devastating plant pathogens worldwide [51,52]. It was first associated with the severe dieback and mortality of Mediterranean oaks, including both cork and holm oak, in the Iberian peninsula by Brasier in 1992 [53]; and since then, P. cinnamomi has been reported across the Mediterranean basin, and this is well documented for European countries [14][15][16][17]44,54]. In this study, the pathogen was detected only from rhizosphere soil; however, it was occasionally isolated from bleeding lesions on the stem (data not shown), as reported by Robin et al. [15]. Although cork oak has been shown to be less susceptible than holm oak [46,55], the high ability of P. cinnamomi to colonize phloem tissues, as exhibited in the log inoculation trial in this study and previously [56], together with its widespread occurrence across Sardinian stands, could suggest less tolerance of the Sardinian cork oak population [57]. This hypothesis, however, needs further investigation.
Phytophthora quercina is the second most prevalent species from soil samples. This oak-specific pathogen has been previously reported in central and southern Europe, causing a chronic decline in Quercus cerris, Quercus faginea, Q. ilex, Quercus petrea, Quercus pubescens, and Quercus robur [20,44,58,59]. Although P. quercina was recorded from two cork oak plantations in Spain [25], this appears to be the first widespread occurrence in natural cork oak forests. Interestingly, two distinct phenotypes were observed amongst the isolates detected from cork oak trees, supporting the hypothesis that P. quercina originated from Europe [21]. Multigene sequencing and phylogenetic analyses are currently underway to investigate the genetic population structure of a large number of P. quercina isolates from different oaks and various geographic provenances (Scanu and Jung, unpublished).
Two other slow-growing and homothallic species are detected at low frequency, and these are identified as P. psychorphila and P. syringae. Together with P. quercina, both species were previously reported from Mediterranean oaks in Spain, and their pathogenicity was demonstrated on both Q. ilex and Q. faginea [20]. Pathogenicity tests in the present study showed both species were not able to invade the inner bark of cut logs of cork oak; however, this did not correlate with root susceptibility. Previous results obtained by Perez-Sierra et al. [20] suggest these species are well adapted to the Mediterranean climate and may act as fine root nibblers, the incidence of which varies depending on the occurrence of extreme climatic events, such as recurrent drought and wet seasons [60]. Both P. psychrophila and P. quercina are already reported in Sardinia from declining Q. ilex trees, while P. syringae is associated with dieback and mortality of Juniperus phoenicea on Caprera Island [28]. Due to their low maximum temperature for growth (around 25 • C), typical of Phytophthora species from cool temperate regions, a potential seasonal activity, as suggested for P. cinnamomi and other cool-temperature pathogens, may occur [61][62][63][64]. Of note, P. syringae was detected only from a new plantation, suggesting its possible new introduction through infected plant material [25,65]. This could also happen for the other Phytophthora species considering the massive afforestation effort in Sardinia between 1990 and 2010 through EU programs, like EEC Regulation 2080/92 [4,25].
The finding of the recently described P. tyrrhenica confirmed its original description from declining cork oak trees [21]. Apart from one site (QS39), it was detected at nine stands and always coexisted with other congeneric species. So far, this cryptic species has been recovered from cork and holm oak trees in Sardinia and Sicily (Italy), respectively, and it is considered endemic to the Mediterranean basin [21,66]. Similarly, P. pseudocryptogea from clade 8 was the only species obtained from QS5 and QS38; otherwise, it co-occurred with P. cinnamomi or other species. It is noteworthy that all the previous isolates identified as Phytophthora cryptogea in Sardinia from forest trees, including oaks and Pinus radiate, are indeed P. pseudocryptogea [28,37,48,67]. Phytophthora pseudocryptogea was recently reported as one of the most widespread species in riparian thermo-Mediterranean forest stands and from five rivers in Sicily [66]. The clade 6 species P. gonapodyides had a more scattered distribution (seven stands and 12 trees), and its finding on cork oak confirmed results obtained by Jung et al. [66]. Finally, for the first time, we recorded the hybrid species P. ×cambivora on Q. suber. In Mediterranean regions, this species is frequently associated with "chestnut ink disease" [68], as well as from other Fagaceae [66,69] and Pinaceae trees [70]. In a recent survey in Sicily, P. ×cambivora was isolated from three different Mediterranean oaks, namely Q. cerris, Q. ilex, and Q. pubescens [66].
The detection of nine Phytophthora species from five phylogenetic clades in 10 rivers within declining cork oak stands was unexpected since similar diversity rates are often reported in more diverse forest ecosystems and with higher sampling rates from surveys in Australia, Europe, the USA, and South Africa [31,32,66]. The Phytophthora assemblage from watercourses was different from that detected from soil samples at the same sites, which was consistent with previous surveys [38,66,71,72]. Only three species were shared between terrestrial and aquatic environments, and these were P. pseudocryptogea, P. gonapodyides, and P. ×cambivora. Interestingly, the most common Phytophthora species isolated from rhizosphere soil, P. cinnamomi, was never detected from streams running through declining cork oak stands. As reported by previous similar studies, clade 6 species were the most common inhabitants of streams and rivers, highlighting their specific lifestyle to aquatic environments [73,74]. Phytophthora gonapodyides occurred in all investigated sites, followed by P. bilorbang, P. lacustris, and P. chlamydospora isolated from five, four, and two watercourses, respectively. Both P. bilorbang and P. gonapodyides have been previously reported in Sardinia on Caprera Island, detected from both rhizosphere soil beneath declining Mediterranean maquis vegetation and holm oak trees [28]. In Italy, P. bilorbang has also been reported from riparian ecosystems of Alnus glutinosa, in Sardinia, and on 15-year-old olive trees in Calabria [75,76], while P. gonapodyides is generally encountered in Mediterranean forest ecosystems, including holm oak [20,44,45]. In the pathogenicity tests, P. gonapodyides caused significant lesions in cork oak logs, confirming its ability to colonize bark and xylem tissues in both artificial and natural infections [77], behaving as a weak pathogen able to survive as a saprophyte on twigs and leaves playing a role in the breakdown of trees debris [78,79]. The other two clade 6 species, P. lacustris and P. chlamydospora, are common species previously reported in Italy from water bodies in forest ecosystems [36,66,80]. Phytophthora plurivora and P. pseudocryptogea are cosmopolitan pathogens with a broad host range [34,37,81]. The DNA of both species was recently detected from holm oak stands across different regions in Spain [23,24], and the pathogenicity tests showed they are amongst the most aggressive species on inoculated holm oak seedlings [82]. The detection of the clade 9 P. hydropathica is not surprising since it was recently detected from river water in Sicily and again using metabarcoding in Spain [24,66]. The origin of P. hydropathica is unknown but considering its low frequency (only one river) and occurrence from ornamental plants in a nursery in Italy [83], a recent introduction into wild environments is most likely.
Of note is the detection of P. ×cambivora and P. ×multiformis, two stable hybrid species from clade 7a that have evolved elsewhere. While P. ×cambivora has been previously reported from Mediterranean oaks in Italy [44], and Q. suber soil in the present study, the isolation of P. ×multiformis most likely occurred due to the presence of Alnus glutinosa trees along the river where the hybrid was isolated. The occurrence of multiple heterozygous sites in the ITS sequences of some isolates, including P. chlamydospora, P. gonapodyides, P. hydropathica, and P. lacustris, together with mixed unreadable ITS sequences generated for some of these isolates could indicate their possible hybrid nature. However, since the ITS region is not a particularly useful locus for studying interspecific hybrids due to the presence of the highly variable non-coding regions ITS1 and ITS2 [83], further molecular analyses, such as cloning, sequencing of other nuclear and mitochondrial genes, estimation of nuclear DNA content by flow cytometry, as well as morphological characterization of the isolates are required [38,84,85].

Conclusions
This unexpected high diversity of Phytophthora species in cork oak stands in such a small geographic area underlines how limited our current knowledge of oomycete diversity in Mediterranean oak forests is [86]. Recent molecular studies in Spanish oak forests have revealed the presence of several Phytophthora species besides P. cinnamomi [12,22,24,87], most of which were detected also in this survey. Although the high-throughput amplicon pyrosequencing of environmental DNA represents a very useful tool for assessing Phytophthora diversity in environmental samples, very little is known on the biological status of the detected microorganisms; therefore, specific baiting technique and metagenomic approaches should be carried out in parallel [72,88].
Pathogenicity tests and results obtained on the susceptibility of Mediterranean oaks to Phytophthora taxa from this and other prior studies [15,46,55,56,62] suggest these previously unrecorded species may play a relevant role in the aetiology of cork oak decline either acting as fine root "nibblers" [20,33,58] or shaping recruitment patterns due to their negative effects on seedling establishment [89,90]. Future studies will be based around understanding the ecological role of all Phytophthora species recovered in such ecosystems as well as their possible interactions with hosts and a changing environment that could promote the establishment of invasive Phytophthora in cork oak forests. Funding: This work was in part financially supported by the Sardinian Regional Government (Convenzione fra la Regione Autonoma Della Sardegna, Assessorato della Difesa dell'Ambiente e l'AGRIS Dipartimento della Ricerca per il Sughero e la Silvicoltura per "L'attuazione del programma biennale di lotta ai lepidotteri defogliatori della sughera"). This study was also supported by the "fondo di Ateneo per la ricerca 2019", an internal funding provided by the University of Sassari.