Martina Cerri1*
Rumakanta Sapkota2
Andrea Coppi3
Valentina Ferri1Bruno Foggi3Daniela Gigante4Lorenzo Lastrucci3Roberta Selvaggi4Roberto Venanzoni4
Mogens Nicolaisen2Francesco Ferranti1
Lara Reale1- 1Department of Agricultural, Food and Environmental Sciences, University of Perugia, Perugia, Italy
- 2Department of Agroecology, Aarhus University, Aarhus, Denmark
- 3Department of Biology, University of Florence, Florence, Italy
- 4Department of Chemistry, Biology and Biotechnology, University of Perugia, Perugia, Italy
Phragmites australis (Cav.) Trin. ex Steud. die-back is a widely-studied phenomenon that was first discovered in northern Europe and that, until recently, was almost unknown in the Mediterranean basin. It has been described as a complex syndrome affecting reed populations leading to their retreat and decline, with significant impacts on valuable ecosystem services. Among the factors that cause the decline, soil-living microorganisms can be crucial. The aims of this study were to analyze the diversity of oomycetes communities associated with reed stands, and to understand whether they could play a key role in the decline. Variations in the structure of oomycetes communities were studied by metabarcoding of the internal transcribed spacer (ITS) 1 region of ribosomal DNA, from the sediments of five Italian freshwater ecosystems. They were chosen to cover a large variability in terms of surface area, water depth, microclimate, and presence of documented reed retreat. From 96 samples collected from reed roots, rhizosphere, and bulk soil, we assembled 207661 ITS1 reads into 523 OTUs. We demonstrated that oomycete communities were structured by several factors, among which the most important was die-back occurrence. Our study also indicates that Pythiogeton spp. could be potentially involved in the development of die-back. The role of heavy metals in the soil was also explored, and cadmium concentration was shown to affect oomycetes distribution. This study represents a significant step forward for the characterization of microbial communities associated with reed die-back syndrome and helps to gain knowledge of the complexity of these important wet ecosystems.
Introduction
Phragmites australis is one of the most widespread angiosperms, with a cosmopolitan distribution range. It forms extensive stands in many types of aquatic habitats, especially in lakes and river shores, marshes, coastal brackish swamps, and lagoons (Engloner, 2009). The importance of this plant is widely acknowledged: it is used in phytoremediation (Guo et al., 2014), it protects the shoreline from bank erosion, and serves as a food resource or protection for arthropods, birds and mammals (Ostendorp et al., 2003). Even if this species sometimes shows an expansive behavior (Bertness et al., 2002; Foggi et al., 2011), since the 1950s an irreversible retreat of the population known as reed die-back syndrome (RDBS) has been observed from several wetlands of northern and central Europe (for a review, see van der Putten, 1997). In Italy, it was first detected in a brackish lagoon (Fogli et al., 2002), and recently in freshwater lakes (Gigante et al., 2011, 2014; Lastrucci et al., 2016a). Typical symptoms of RDBS include reduced plant height, weaker culms, abnormal rhizomes, formation of clumps (Armstrong et al., 1996; Fogli et al., 2002), and eventually retreat of the population (van der Putten, 1997). Some of the possible abiotic and biotic causes of reed decline has been widely investigated (Brix, 1999; Fogli et al., 2002). Among the abiotic factors, prolonged submersion of the plant appears to have dramatic consequences; Lastrucci et al. (2016a) demonstrated the existence of a relationship between reed die-back and permanent-artificially induced flooding conditions. Several macromorphological traits associated with RDBS (such as clumping habit, reduced culm height or diameter, high rate of dead buds, see Gigante et al., 2011, 2014; Lastrucci et al., 2016a for details) were significantly more often observed in permanently flooded stands than in the emerged stands. Furthermore, even if P. australis can tolerate high concentrations of heavy metals (Bonanno and Lo Giudice, 2010), some studies have highlighted their negative effect on plant fitness, suggesting a role in the RDBS (Gigante et al., 2014; Lastrucci et al., 2016a). Cd and Co, for example, seemed to be related to some traits associated with the die-back symptoms in Montepulciano (Lastrucci et al., 2016a).
The understanding of the possible biotic causes of RDBS is yet limited; the role of microorganisms living in the soil is still an issue of concern and deserves to be investigated, since plants are metaorganisms establishing close symbiotic relationships with their associated microorganisms (Mendes et al., 2013; Coats and Rumpho, 2014; Vandenkoornhuyse et al., 2015; Berg et al., 2016). To date, only a few studies have been carried out about the interaction between P. australis and bacteria or fungi (Wirsel et al., 2002; Neubert et al., 2006; Clay et al., 2016; Soares et al., 2016), and even less studies have explored the relationship between P. australis and communities of oomycetes in natural ecosystems (Nechwatal et al., 2008). Several oomycete pathogens have been shown to be involved in root disease complexes in many Poaceae (Toda et al., 2015; Grijalba et al., 2017). A previous study conducted in Lake Constance (Germany) suggested an association of the oomycete Pythium phragmitis with reed decline (Nechwatal et al., 2005). Therefore, we decided to investigate the involvement of oomycetes in RDBS in our study area.
Phylogenetic studies have shown that oomycetes are a diverse group of fungus-like eukaryotic microorganisms closely related to diatoms and sea weeds, and probably linked to the marine environment during evolution (Thines, 2014). They have colonized almost all ecosystems, from semi-arid regions (Mirzaee et al., 2013) to Arctic and Antarctic (Bridge et al., 2008) and include both saprophytes and pathogens of plants, insects, crustaceans, fish, vertebrate animals, and various microorganisms.
Culture-based techniques have allowed isolated microbes to be studied in detail, but they fail to represent the full microbial diversity inhabiting the natural environment, limiting the analysis to those that grow under laboratory conditions (Rondon et al., 2000; Fierer et al., 2007; Mendes et al., 2013). Molecular techniques, such as metabarcoding, overcomes these difficulties and allows the identification of microbes in situ (Ravin et al., 2015). Metabarcoding studies of oomycetes are limited in number and in efficacy, but recently, an improved strategy developed by Sapkota and Nicolaisen (2015) and based on a new annealing temperature, was demonstrated to successfully uncover the complexity of these microorganisms in soils.
The present study represents the first investigation of the diversity of reed-associated oomycetes; additionally, it is the first contribute to assess whether declining reed-beds harbor unique communities of oomycetes that differ from those of the healthy stands, and to understand the factors that influence their distribution. We evaluated the effect of the location, the die-back occurrence and submersion as driving forces in determining community composition. Henderson et al. (2017) suggested that heavy metals can inhibit the growth of oomycetes growth, being particularly toxic at zoospores stage, and with toxicity in order of Ag > Cu > Ni > Co > Zn; therefore, the effect of heavy metals on reed-associated oomycetes communities was investigated.
Materials and Methods
Site Description
Five wetlands from Central Italy were considered in this study: (1) Colfiorito Marsh, in Umbria Region (43°01′23″ N, 12°52′36″ E) at 756 m a.s.l.; (2) Lake Trasimeno, in Umbria Region (43°08′05.5″ N, 12°06′04.6″ E) at an average altitude of 257 m a.s.l.; (3) Lake Vico in Lazio Region (42°18′58.40″ N, 12°10′5.89″ E), 507 m a.s.l.; (4) Lake Chiusi in Tuscany Region (43°03′22.11″ N, 11°57′55.79″ E) at 252 m a.s.l.; (5) the wetland of “Le Morette” (43°48′30.38″ N, 10°48′20.14″ E) at Padule of Fucecchio Marsh, in Tuscany, at 13 m a.s.l. (Figure 1). These five ecosystems were chosen to cover a large variability in terms of: (i) surface area, from large lakes such as Trasimeno and Vico, to the small wetland of Le Morette and Colfiorito marshes; (ii) water depth, from the deep waters of Vico, to the shallow waters of Le Morette or Colfiorito marshes; (iii) microclimate, from temperate Apennine areas such as Colfiorito, to the submediterranean climate of Vico or Chiusi lakes; (iv) previous reed retreat documented by aerial images (Gigante and Venanzoni, 2012; Gigante et al., 2014).
FIGURE 1. Map showing the distribution of Italian lakes selected for sampling.
Soil and Roots Sampling
The sampling design was based on the results of former investigations in the same areas (Lastrucci et al., 2016b). Two types of stands could be identified for sampling: healthy stands (group N on Supplementary Table S1) and RDBS stands (group Y on Supplementary Table S1), the latter showing a number of RDBS diagnostic traits and, above all, the clumping habit. This is a clustered growth of the reed stems deriving from an uncontrolled outgrowth of dormant buds caused by the breaking of apical dominance, and represents a widely acknowledged and easily detectable trait of reed die-back (Armstrong et al., 1996; van der Putten, 1997; Clevering, 1998; Dinka and Szeglet, 2001; Gigante et al., 2011). Sampling was carried out during two different periods: in February 2015, bulk and rhizosphere soils were sampled from three healthy stands and three RDBS stands from four wetlands (Colfiorito, Trasimeno, Vico, and Chiusi) (48 samples); in September 2015, two healthy stands and two RDBS stands were sampled from all five wetlands, collecting, where possible, bulk, rhizosphere, and roots (48 samples). Bulk soil (40 g for molecular analysis and 300 g for chemical analysis) was collected at an approximate depth of 30 cm, 50 cm away from the plant, using a soil core sampler. Belowground part of the stem with roots were collected; roots were cut and thoroughly washed in distilled water for several times, while the soil attached to the roots was considered as rhizosphere. Soil and root samples were stored at -20°C until DNA extraction. For a summary of the samples used in this study, see Supplementary Table S1.
DNA Extraction, DNA Amplification, and Pyrosequencing
Two different procedures were used for DNA extraction, one for soil samples (bulk soil and rhizosphere) and the second one for root samples. Soil samples weighing 40 g were mixed, freeze dried for 48 h, and ground; from these samples, 250 mg were used for total soil DNA extraction using the PowerLyzerTM PowerSoil® DNA Isolation Kit (Mo Bio Laboratories, Carlsbad, CA, United States) according to the manufacturer’s instructions, except that samples were homogenized in a Geno/Grinder 2000 (SPEX CertiPrep, Metuchen, NJ, United States) at 1500 rpm for three times, 30 s each. For root samples, 100 mg of fresh material were ground in a mortar with liquid nitrogen; DNA was extracted using the DNeasy Plant Mini kit (Quiagen GmbH, Hilden, Germany) according to the manufacturer’s instructions.
The internal transcribed spacer (ITS) 1 region of ribosomal DNA was used as a marker. Each sample was barcoded with 10-nucleotide multiplex identifier (MID) primer tags before pooling PCR products (Supplementary Table S1). Amplification strategy was as described by Sapkota and Nicolaisen (2015): the PCR reaction mixture consisted of 1× PCR reaction buffer, 1.5 mM of MgCl2, 0.2 mM of dNTPs, 1 μM of each primer, and 1 U of GoTaq Flexi polymerase (Promega Corporation, Madison, WI, United States) in a total volume of 25 μl containing 24 μl of reaction mixture and 1 μl of template. Tagged PCR amplicons from the 96 samples were pooled in equimolar amounts, electrophoresed on a 1.5% agarose gel and the visible smear of 200–400 base pairs was excised and purified from the gel using QIAquick Gel Extraction Kit (Quiagen, Hilden, Germany). The pooled samples were sequenced by Eurofins MWG (Ebersberg, Germany) on a GS Junior 454 Sequencer (Roche Diagnostics) using titanium chemistry. The raw data were deposited on the Sequence Read Archive website1 under the SRA number SRP111776.
Sequence Analysis and Statistics
Sequences were analyzed using the bioinformatics pipeline available in QIIME v. 1.8 (Caporaso et al., 2010). De-multiplexing and quality filtering was carried out using split_libraries.py command using quality score of 50, discarding reads less than 150, homopolymer length of 10 nucleotides, maximum number of ambiguous bases of 6 and other default settings. The ITS1 region was extracted using ITSx extractor version 1.0.6 (Bengtsson-Palme et al., 2013). ITS1 reads were then clustered using the pick _otus.py script at 98% similarity level (Vettraino et al., 2012) using usearch61 with a minimum cluster size of two, thus excluding singletons (Edgar, 2010). For taxonomic assignment, a reference database described earlier (Sapkota and Nicolaisen, 2015) was used. In addition, for species identification, a representative sequence from each of the most abundant operational taxonomic unit (OTUs) with at least 100 reads in total was further subjected to Basic Local Alignment Search Tool (BLAST) matches at NCBI. At least the top five BLAST hits were evaluated: OTUs with identity and coverage ≥ 95% were assigned at species level, or at genus or family level if lower. Only reads assigned to oomycetes were retained in the OTU table.
The OTU table was used for diversity based calculations. Data visualization and statistical analysis were carried out in R (R Core Team, 2015). A species accumulation curve was generated using the specaccum function available in the ‘vegan’ package (Oksanen et al., 2015). Samples representing less than 400 reads were removed before performing diversity based calculations. Species richness and evenness were calculated by rarifying the OTU table at an even sampling depth of 400 reads, while the OTU table was transformed to relative abundance before beta diversity based calculations. The species evenness was evaluated by Pielou’s index and the species richness was estimated as the number of OTUs. β-diversity distance matrices were obtained using Bray–Curtis dissimilarity and visualized using principal coordinates analysis (PCoA) plot. Permutational multivariate analysis of variance (Anderson, 2001) with a permutation number of 999 was carried out to test the null hypothesis of no differences among a priori defined groups using the adonis function in ‘vegan’ package.
Covariance biplots were calculated using composition data analysis on centered log-ratio transformed OTU table (Gloor and Reid, 2016), removing samples less than 400 reads.
The OTU table was also used to perform the Dufrene-Legendre indicator species analysis in R, using the labdsv package (Roberts, 2012) to test the OTUs associated with die-back. OTUs with significant values (P < 0.05) and an indicator value > 0.25 were set as cut-off values to define indicator species (Dufrene and Legendre, 1997). The difference in relative abundance of OTUs was determined by non-parametric statistical method using Kruskal–Wallis test for ANOVA and Wilcoxon rank sum tests for pairwise comparison. The network analysis was made in QIIME using OTU table (with the command make_otu_network.py) and the resulting network files were visualized in cytoscape2. Samples and OTUs were subjected to network analysis via the spring-embedded algorithm using eweights, where samples represent larger nodes than OTUs. All OTUs are connected to the sample to which they belong with a line (edge). Shared OTUs are thus connected to several samples, and occupy the center of network, whereas OTUs that only occur in one or few samples occupy the periphery of the network. Samples are thus clustered based on shared OTU’s and the degree to which samples cluster is based on the number of OTUs shared between samples.
Phylogenetic Analysis and Assessment of Taxa Associated to the RDBS
Using the indicator species analysis, 12 OTUs were found to have significantly higher abundance in die-back samples in the total dataset; they were aligned together with the ITS1 sequences of the accessions of Pythiogeton available in GenBank. ITS1 sequences of the two most abundant Pythium species in our dataset were used as outgroup. The Neighbor-joining tree was assembled with MEGA 7.0 (Kumar et al., 2015) and tested with a bootstrap of 1000 replicates, to ascertain the reliability of any given branch pattern.
Determination of Heavy Metals in Soil and Plant Tissues
Concentrations of heavy metals (Cd, Pb, Zn, Cr, Ni, Cu, Al) in soil samples and in plant tissues were determined by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES, Ultima 2, HORIBA Scientific) equipped with an ultrasonic nebulizer (U-5000AT, CETAC Technologies). The nebulizer gas flowrate was kept at 1 L min-1 and the plasma operated at a radiofrequency power of 1000 W. Supplementary Table S2 shows the wavelength used. Instrumental detection limits were in the range of 0.1–1.5 mg kg-1. Commercially produced (ICP multi-element standard solution IV CertiPUR®, VWR Merck Chemicals and Reagents) standard solutions (1000 mg L-1) in nitric acid were used to prepare appropriate elemental calibration standards. An acid digestion was performed on the soil (Ministero dell’Ambiente e della Tutela del Territorio [MATT] and Agenzia per la Protezione dell’Ambiente e per i Servizi Tecnici [APAT], 2005) and on the plant tissue samples (Du Laing et al., 2003) prior the analytical determination. The soils were air-dried, disaggregated using a mortar and pestle to pass through a 2 mm mesh sieve, dried at 105°C for 24 h and digested as follows: 8 ml of ultrapure nitric acid (Millipore Suprapur®, 65%) and 2 ml of ultrapure solution of hydrogen peroxide (Millipore Suprapur®, 30%) were added to 0.200 g of soil samples and digested in a Mars Microwave Oven, working at a power of 1000 W. Microwave digestion consisted of two steps: 130°C (200 psi) for 1 min, 180°C (300 psi) for 10 min. The mixture was cooled, filtered (Whatman Grade No. 42, particle retention 2.5 mm) and diluted with ultrapure water to 50 ml. Plant samples were preliminarily dissected into roots and rhizomes and dried at 70°C for 48 h. Acid digestion of tissues (0.4 g) was carried by the same procedure used for soils, but ultrapure water was added to reach the volume of 25 ml. The measure of the heavy metal was based on two replicates of the overall analytical procedure for each sample of soil and plant tissue examined.
Results
Data Structure and Taxonomic Composition
In total, 235065 reads were obtained from 454-pyrosequencing. After quality control, ITS extraction and taxonomic assignments, a total of 207661 reads were clustered into 523 OTUs assigned to oomycetes, excluding singletons. Two samples failed to yield any reads; thus, reads were obtained from 94 samples in total. The number of reads per sample revealed high variation (min, mean, max: 4, 2464, 13461).
The accumulation curve showed that the number of OTUs hardly approached a plateau phase (Figure 2), indicating a higher oomycetes diversity in the area than was covered in this study and that diversity would increase with increasing sample number. In our dataset, Pythiales dominated the reads (66%), followed by Lagenidiales (25%), Saprolegniales (6.5%) and Peronosporales (2.4%), while at genus level we detected the presence of three dominating taxa: Pythiogeton (30.9%), Pythium (29.9%), and Lagena (22%) (Figure 3).
FIGURE 2. Species accumulation curve for the 96 samples, obtained with the function specaccum of the ‘vegan’ R package; it shows the relationship between observed OTUs and number of samples.
FIGURE 3. Boxplot showing the relative abundance of the most abundant taxa of oomycetes communities at order, family, and genus levels. Box borders represent the first and third quartiles, and central lines represent medians.
Oomycete Diversity Analysis
The incidence (the number of samples in which a certain OTU occurred) was checked for the 523 OTUs identified in this study: 166 OTUs were found in only one sample, while only one OTU (denovo0; Lagena) was found in more than 50% of the samples (see Supplementary Table S3). Lagena sp. (denovo0) had a significantly higher relative abundance in root samples in comparison to rhizosphere (Kruskal–Wallis test, P = 0.023) and bulk soil samples (Kruskal–Wallis test, P = 0.017) (Figure 4). The species richness was compared in relation to different sites, plant compartments (root, rhizosphere, or bulk soil), temporary or permanent flooding (i.e., absence or presence of die-back). Interestingly, none of the studied factors significantly influenced OTUs richness. The species evenness estimated by Pielou’s index, is significantly influenced by the compartment (df = 2; P = 0.017). Pairwise Wilcoxon rank sum tests revealed that the evenness is significantly different between root and bulk soil samples (P = 0.018).
FIGURE 4. Lagena sp. distribution in bulk soil, rhizosphere, and root samples. Boxplot show the 25% and 75% quartiles of the data, the central line represents the median, bars extend to 95% confidence limits.
Factor Shaping Oomycete Communities
The effects of the geographical location, the compartment (bulk soil, rhizosphere soil, or root sample), die-back occurrence, and the available chemical data were subjected to permutation multivariate analysis of variance to evaluate factors that affected the structure (beta diversity) of oomycetes communities. The factors that were found to be significant in explaining the diversity were die-back occurrence (R2 = 0.049, P = 0.001), compartment (R2 = 0.043, P = 0.009), and the interaction between lake and die-back occurrence (R2 = 0.073, P = 0.011) (Table 1). Among the heavy metals, only cadmium concentration (mg⋅Kg-1) is significantly related to oomycetes distribution (R2 = 0.035, P = 0.002) (Table 2).
TABLE 1. Results from permutational multivariate analysis of variance to partition the variance in oomycete communities for the combined data. R_B_RT indicates the compartment: R, rhizosphere; B, bulk soil; RT, roots.
TABLE 2. Results from permutational multivariate analysis of variance to test the effect of heavy metals in the soils on oomycetes communities.
In support of this, a PCoA plot revealed clustering based on the plant health status, demonstrating the importance of die-back in explaining the diversity of oomycetes (Figure 5). The clustering of samples based on die-back was supported by covariance biplots (Supplementary Figure S1). The influence of die-back in oomycetes communities was also visible via network based analysis (Figure 6). The samples with higher number of shared OTUs were closer and clustered in the center of the network. The majority of die-back samples clustered in the lower core of the network. The network analysis also showed that a large part of OTUs were shared between healthy and RDBS samples.
FIGURE 5. Principal coordinates analysis (PCoA) plot based on Bray–Curtis matrix. Samples are colored based on RDBS occurrence.
FIGURE 6. Oomycetes OTU networks in healthy and die-back samples. Larger nodes indicate samples, whereas smaller nodes represent OTUs. All OTUs are connected to the sample to which they belong with a line (edge) making shared OTUs connected to more than one sample. Healthy samples and OTUs as well as lines are indicated with blue, and die-back samples with red. Samples with shared OTUs are located in the center of the plot, while samples with unique OTUs are located in the periphery.
Difference between Healthy and RDBS Samples at OTU Level
Indicator species analysis of the total dataset revealed 12 OTUs that were associated with RDBS (Table 3).
TABLE 3. List of indicator species that were significantly different calculated using the Dufrene and Legendre analysis in R, considering the total dataset, only roots samples, and only rhizosphere samples.
Those 12 OTUs belonged to the genus Pythiogeton. The relative abundance of Pythiogeton was significantly higher in die-back samples in comparison to healthy (P = 0.0002) (Supplementary Figure S2). A phylogenetic tree (Figure 7) shows that eight of the OTUs (denovo594, denovo1113, denovo208, denovo8, denovo3, denovo245, denovo302, denovo246) are similar to Pythiogeton microzoosporum, although the clades do not show a high bootstrap support. Denovo1 and denovo284 are similar (bootstrap = 98) and close to the clade containing Pythiogeton oblongilobum (bootstrap = 63). Denovo34 and denovo152 cluster together (bootstrap = 100) and they are similar to Pythiogeton abundance.
FIGURE 7. Neighbor-joining phylogenetic tree based on ITS1 region showing the diversity of Pythiogeton species in our root, rhizosphere, and bulk soil samples. Sequences from this study are shown with denovo numbers, whereas reference sequences are listed with the NCBI accession numbers. The two most abundant Pythium species were used as outgroups. Numbers at branching points are bootstrap percentages based on 1000 replications. Only values > 50% are shown.
Discussion
Die-back syndrome is a complex phenomenon affecting mature stands of P. australis. It represents a serious threat to reed-dominated ecosystems, given the importance of the plant in protecting the shoreline from bank erosion and in providing food resource or protection for arthropods, birds and mammals (Ostendorp et al., 2003). Even though different abiotic factors (for example, permanent flooding of the stand and presence of heavy metals) have been suggested as possible causes of this syndrome (Gigante et al., 2014; Lastrucci et al., 2016a), a clear understanding of this phenomenon is still missing. The role of biotic factor is also poorly understood; in particular, the interaction between plant and microorganisms is not completely explored yet. To this aim, the present study analyzed the distribution of oomycetes communities in five freshwater ecosystems and investigated a possible correlation with the health status of P. australis.
As reported by Wielgoss et al. (2009), the oomycete P. phragmitis was considered to be often associated with reed die-back in Lake Constance, Germany, especially if co-occurring with permanent flooding condition of the plant. Interestingly, we did not detect P. phragmitis, in any of our 96 samples. We attempted a qPCR approach, using the primers described by Wielgoss et al. (2009) for the detection of P. phragmitis, however, with no success (data not shown). The reason for the absence of P. phragmitis in samples from Italy could be the different bioclimatic conditions between Germany and Italy.
In our study, instead, indicator species analysis revealed 12 OTUs which showed significant differences in abundance between stands where P. australis is affected by die-back syndrome and stands in good conservation status. All the OTUs identified as indicator species in RDBS belong to the genus Pythiogeton. Pythiogeton is a poorly studied oomycetes genus containing 16 known species (MycoBank3). It is difficult to study because of its reduced capacity to grow in culture; indeed, it was believed to be a saprophyte until recent investigations proved its pathogenicity, also toward Poaceae members (Doan et al., 2013; Huang et al., 2013). As already demonstrated (Natvig and Gleason, 1983; Huang et al., 2013), Pythiogeton species are more virulent in anaerobic condition, and this could be the reason why they were detected in die-back stands. In fact, reed die-back, as showed in several former studies (Gigante et al., 2011, 2014; Lastrucci et al., 2016a) and even confirmed in the presently considered areas (Lastrucci et al., 2016b), mostly affects stands with prolonged or even permanent conditions of submersions. Indeed, in the study areas the stands affected by RDBS (group Y) are permanently flooded, with water persistence even in the driest period (i.e., the end of summer), while those in a good health status (group N) are only temporarily flooded, emerging at least at the end of summer. A phylogenetic tree was constructed to understand the relationship between the known species of Pythiogeton and the OTUs assembled in this study. Eight of them are closely related to P. microzoosporum, two OTUs were similar to the branch with P. oblongilobum, and two were close to P. abundance. P. oblongilobum was isolated by Huang et al. (2013) in the stem debris of the water bamboo (Zizania latifolia). P. abundance was firstly isolated from decaying leaves of common cattail (Drechsler, 1932), and later by Huang et al. (2013). In previous investigations (Xu et al., 2012; Nelson and Karp, 2013; Sapkota and Nicolaisen, 2015) Pythiogeton was not detected at all, or only at low levels; instead, in our work, it was the most abundant genus, together with Pyhtium. One reason can be in the different approaches used: Nelson and Karp (2013) used the cloning, Xu et al. (2012) used different primer pairs to include both fungi and oomycetes, while we used primers specific for oomycetes. Another obvious reason can be in the different considered ecosystem, as Sapkota and Nicolaisen (2015) worked on carrot fields in Denmark, while we work on freshwater ecosystems. At last, the co-occurrence of Pythiogeton with RDBS could also be the reason of the findings.
In the PCoA plot, most of the samples from healthy stands clustered together, while RDBS samples were spread in the plot. Nevertheless, the network analysis indicated that a large part of the OTUs are shared between the two categories, healthy or die-back affected samples, and the main differences between them were the differences in Pythiogeton spp. distribution. These results supported a possible relation between die-back occurrence and oomycetes communities composition.
Interestingly, we found that OTUs richness is not statistically different among the considered lakes, nor between compartments or health status of the plant. Species evenness, on the other hand, was significantly affected by the compartment (root, rhizosphere, bulk soil). This means that species are not evenly distributed and their relative abundance change based on the sample type. Lagena sp. was indeed more widely present in the roots samples, compared to rhizosphere and bulk soil; it is an obligate intercellular parasite of plant roots (Barr and Dèsaulniers, 1990), observed also in wild grasses (Blackwell, 2011; Bakker et al., 2017), but little investigation has been done to elucidate its taxonomy and pathogenicity because of difficulties in cultivation. In this study, hardly any differences in the relative abundance of Lagena were detected between healthy and declining stands, indicating that it cannot be related to RDBS development in our samples.
In this study, we also took the chemistry of the soils into account, as previous investigations highlighted a possible role of heavy metals in oomycetes growth (Henderson et al., 2017). Among the heavy metals analyzed, the concentration of cadmium was the only one that showed a significant relationship with oomycete distribution. The toxicity of cadmium for soil microbial biomass growth and metabolic activity is well-known (Xu et al., 2013) and Cd2+ concentrations has been demonstrated to affect soil microbial diversity (Chien et al., 2008; Ribeiro et al., 2016). Several studies have also demonstrated the effect of Cadmium on P. australis fitness (Ederli et al., 2004) and in RDBS (Lastrucci et al., 2016a).
To summarize, RDBS remains a very intricate phenomenon, which is the result of the interaction between abiotic and biotic factors. Stress conditions determined by flooding and heavy metals could in fact be worsened by oomycetes presence. This is the first metabarcoding study of oomycete diversity in reed beds ecosystems. Our approach allowed to highlight results that would have been missed with cultivation-based studies, and demonstrated the presence of Pythiogeton as associated to declining reed stands, which occur to be permanently flooded. Hopefully, these new insights will contribute to unravel the diversity of these poorly studied organisms living in reed-dominated ecosystems.
Author Contributions
MC, MN, LR, and FF designed the work; MC, VF, DG, AC, and LL collected the samples; MC, RSe, and RSa acquired and analyzed the data; MC drafted the manuscript; RSe, MN, DG, LR, AC, LL, FF, RV, BF, RSa critically revised the article. All the authors approved the version of the manuscript to be published.
Funding
Partial funding of this research was provided by the Italian Ministry of Education, University, and Research (Grant RBFR13P7PR).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgment
The authors would like to thank S. E. Rasmussen, B. R. Eleuteri, L. Picchiarelli, M. Muzzatti, and M. Chiappini for the kind support.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01550/full#supplementary-material
Footnotes
References
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46.
Armstrong, J., Armstrong, W., and Van der Putten, W. (1996). Phragmites die-back: bud and root death blockages within the aeration and vascular systems. New Phytol. 133, 339–414. doi: 10.1111/j.1469-8137.1996.tb01907.x
Bakker, M. G., Moorman, T. B., and Kaspar, T. C. (2017). Isolation of cultivation-resistant oomycetes, first detected as amplicon sequences, from roots of herbicide-terminated winter rye. Phytobiomes J. 1, 24–35. doi: 10.1094/PBIOMES-10-16-0011-R
Barr, D. J. S., and Dèsaulniers, N. L. (1990). The life cycle of Lagena radicicola, an oomycetous parasite of wheat roots. Can. J. Bot. 68, 813–824. doi: 10.1139/b90-108
Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., et al. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol. Evol. 4, 914–919. doi: 10.1111/2041-210X.12073
Berg, G., Rybakova, D., Grube, M., and Köberl, M. (2016). The plant microbiome explored: implications for experimental botany. J. Exp. Bot. 67, 995–1002. doi: 10.1093/jxb/erv466
Bertness, M. D., Ewanchuk, P. J., and Silliman, B. R. (2002). Anthropogenic modification of New England salt marsh landscapes. Proc. Natl. Acad. Sci. U.S.A. 99, 1395–1398. doi: 10.1073/pnas.022447299
Blackwell, W. H. (2011). The genus Lagena (Stramenopila: Oomycota), taxonomic history and nomenclature. Phytologia 93, 157–167.
Bonanno, G., and Lo Giudice, R. (2010). Heavy metal bioaccumulation by the organs of Phragmites australis (common reed) and their potential use as contamination indicators. Ecol. Indic. 10, 639–645. doi: 10.1016/j.ecolind.2009.11.002
Bridge, P. D., Newsham, K. K., and Denton, G. J. (2008). Snow mould caused by a Pythium sp.: a potential vascular plant pathogen in the maritime Antarctic. Plant Pathol. 57, 1066–1072. doi: 10.1111/j.1365-3059.2008.01868.x
Brix, H. (1999). The European research project on reed die-back and progression (EUREED). Limnologica 29, 5–10. doi: 10.1016/S0075-9511(99)80033-4
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chien, C., Kuo, Y., Chen, C., Hung, C., Yeh, C., and Yeh, W. (2008). Microbial diversity of soil bacteria in agricultural field contaminated with heavy metals. J. Environ. Sci. 20, 359–363. doi: 10.1016/S1001-0742(08)60056-X
Clay, K., Shearin, Z. R. C., Bourke, K. A., Bickford, W. A., and Kowalski, K. P. (2016). Diversity of fungal endophytes in non-native Phragmites australis in the Great Lakes. Biol. Invasions 18, 2703–2716. doi: 10.1007/s10530-016-1137-y
Clevering, O. A. (1998). Effects of litter accumulation and water table on morphology and productivity of Phragmites australis. Wetl. Ecol. Manag. 5, 275–287. doi: 10.1023/A:1008233912279
Coats, V. C., and Rumpho, M. E. (2014). The rhizosphere microbiota of plant invaders: an overview of recent advances in the microbiomics of invasive plants. Front. Microbiol. 5:368. doi: 10.3389/fmicb.2014.00368
Dinka, M., and Szeglet, P. (2001). Some characteristics of reed (Phragmites australis/Cav./Trin ex Steudel) that indicate different health between vigorous and die-back stands. Verh. Internat. Ver. Limnol. 27, 3364–3369.
Doan, H. K., Davis, R. M., Sartori, F. F., and Marcum, D. B. (2013). First report of a Pythiogeton sp. causing root and basal stalk rot of wild rice in California. Plant Dis. 98, 851. doi: 10.1094/PDIS-10-13-1098-PDN
Drechsler, C. (1932). A specie of Pyhtiogeton isolated from decaying leaf-sheaths of the common cat-tail. J. Wash. Acad. Sci. 22, 421–449.
Du Laing, G., Tack, F. M. G., and Verloo, M. G. (2003). Performance of selected destruction methods for the determination of heavy metals in reed plants (Phragmites australis). Anal. Chim. Acta 497, 191–198. doi: 10.1016/j.aca.2003.08.044
Dufrene, M., and Legendre, P. (1997). Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol. Monogr. 67, 345–366. doi: 10.2307/2963459
Ederli, L., Reale, L., Ferranti, F., and Pasqualini, S. (2004). Responses induced by high concentration of cadmium in Phragmites australis roots. Physiol. Plant. 121, 66–74. doi: 10.1111/j.0031-9317.2004.00295.x
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Engloner, A. I. (2009). Structure, growth dynamics and biomass of reed (Phragmites australis) – A review. Flora Morphol. Distrib. Funct. Ecol. Plants 204, 331–346. doi: 10.1016/j.flora.2008.05.001
Fierer, N., Breitbart, M., Nulton, J., Salamon, P., Lozupone, C., Jones, R., et al. (2007). Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of bacteria, archaea, fungi, and viruses in soil. Appl. Environ. Microbiol. 73, 7059–7066. doi: 10.1128/AEM.00358-07
Foggi, B., Lastrucci, L., Viciani, D., Brunialti, G., and Benesperi, R. (2011). Long-term monitoring of an invasion process: the case of an isolated small wetland on a Mediterranean Island. Biologia 66, 638–644. doi: 10.2478/s11756-011-0057-7
Fogli, S., Marchesini, R., and Gerdol, R. (2002). Reed (Phragmites australis) decline in a brackish wetland in Italy. Mar. Environ. Res. 53, 465–479. doi: 10.1016/S0141-1136(02)00091-0
Gigante, D., Angiolini, C., Landucci, F., Maneli, F., Nisi, B., Vaselli, O., et al. (2014). New occurrence of reed bed decline in southern Europe: do permanent flooding and chemical parameters play a role? C. R. Biol. 337, 487–498. doi: 10.1016/j.crvi.2014.05.005
Gigante, D., and Venanzoni, R. (2012). “Il declino della popolazione di Phragmites australis al Lago Trasimeno,” in Tutela Ambientale del lago Trasimeno, ed. A. Martinelli (Umbria: Libri/A.R.P.A), 109–120.
Gigante, D., Venanzoni, R., and Zuccarello, V. (2011). Reed die-back in southern Europe? A case study from Central Italy. C. R. Biol. 334, 327–336. doi: 10.1016/j.crvi.2011.02.004
Gloor, G. B., and Reid, G. (2016). Compositional analysis: a valid approach to analyze microbiome high-throughput sequencing data. Can. J. Microbiol. 62, 692–703. doi: 10.1139/cjm-2015-0821
Grijalba, P. E., Palmucci, H. E., and Guillin, E. (2017). Identification and characterization of Pythium graminicola, causal agent of kikuyu yellows in Argentina. Trop. Plant Pathol. 42, 284–290. doi: 10.1007/s40858-017-0149-1
Guo, L., Ott, D. W., and Cutright, T. J. (2014). Accumulation and histological location of heavy metals in Phragmites australis grown in acid mine drainage contaminated soil with or without citric acid. Environ. Exp. Bot. 105, 46–54. doi: 10.1016/j.envexpbot.2014.04.010
Henderson, L., Lilje, E., Robinson, K., Gleason, F. H., and Lilje, O. (2017) “Effects of toxic metals on Chytrids, fugal-like organisms, and higher fungi” in The Fungal Community: Its Organization and Role in the Ecosystem, 4th Edn, eds J. Daton and J. F. White (Boca Raton, FL: CRC Press).
Huang, J. H., Chen, C. Y., Lin, Y. S., Ann, P. J., Huang, H. C., and Chung, W. H. (2013). Six new species of Pythiogeton in Taiwan, with an account of the molecular phylogeny of this genus. Mycoscience 54, 130–147. doi: 10.1016/j.myc.2012.09.007
Kumar, S., Stecher, G., and Tamura, K. (2015). MEGA7: molecular evolutionary genetics analysis version 7.0. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Lastrucci, L., Gigante, D., Vaselli, O., Nisi, B., Viciani, D., Reale, L., et al. (2016a). Sediment chemistry and flooding exposure: a fatal cocktail for Phragmites australis in the Mediterranean basin? Ann. Limnol. 52, 365–377. doi: 10.1051/limn/2016023
Lastrucci, L., Lazzaro, L., Coppi, A., Reale, L., Ferranti, F., Venanzoni, R., et al. (2016b). Indagini Cartografiche e Macro-Morfologiche per L’individuazione e lo Studio del Declino di Phragmites australis in Italia Centrale. Abetone: 50° Congresso SISV.
Mendes, R., Garbeva, P., and Raaijmakers, J. M. (2013). The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663. doi: 10.1111/1574-6976.12028
Ministero dell’Ambiente e della Tutela del Territorio [MATT] and Agenzia per la Protezione dell’Ambiente e per i Servizi Tecnici [APAT] (2005). Progetto Nazionale di Monitoraggio Acque Superficiali (IRSA, CNR), Gli Ecosistemi e i Sedimenti: Caratterizzazione dei Sedimenti, 2005. Rome: Ministero dell’Ambiente e della Tutela del Territorio.
Mirzaee, M. R., Ploch, S., Runge, F., Telle, S., Nigrelli, L., and Thines, M. (2013). A new presumably widespread species of Albugo parasitic to Strigosella spp. (Brassicaceae). Mycol. Prog. 12, 45–52. doi: 10.1007/s11557-012-0811-4
Natvig, D. O., and Gleason, F. H. (1983). Oxygen uptake by obligately- fermentative aquatic fungi: absence of a cyanide-sensitive component. Arch. Microbiol. 134, 5–8. doi: 10.1007/BF00429397
Nechwatal, J., Wielgoss, A., and Mendgen, K. (2005). Pythium phragmitis sp. nov., a new species close to P. arrhenomanes as a pathogen of common reed (Phragmites australis). Mycol. Res. 109, 1337–1346. doi: 10.1017/S0953756205003990
Nechwatal, J., Wielgoss, A., and Mendgen, K. (2008). Flooding events and rising water temperatures increase the significance of the reed pathogen Pythium phragmitis as a contributing factor in the decline of Phragmites australis. Hydrobiologia 613, 109–115. doi: 10.1007/s10750-008-9476-z
Nelson, E. B., and Karp, M. A. (2013). Soil pathogen communities associated with native and non- native Phragmites australis populations in freshwater wetlands. Ecol. Evol. 3, 5254–5267. doi: 10.1002/ece3.900
Neubert, K., Mendgen, K., Brinkmann, H., and Wirsel, S. G. R. (2006). Only a few fungal species dominate highly diverse mycofloras associated with the common reed. Appl. Environ. Microbiol. 72, 1118–1128. doi: 10.1128/AEM.72.2.1118-1128.2006
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., et al. (2015). vegan: Community Ecology Package. R Package Version 2.2-1. Available at: http://CRAN.R-project.org/package=vegan
Ostendorp, W., Dienst, M., and Schmieder, K. (2003). Disturbance and rehabilitation of lakeside Phragmites reeds following an extreme flood in Lake Constance (Germany). Hydrobiologia 50, 687–695. doi: 10.1023/B:HYDR.0000008622.60094.6d
R Core Team (2015). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Ravin, N. V., Mardanov, A. V., and Skryabin, K. G. (2015). Metagenomics as a tool for the investigation of uncultured microorganisms. Russ. J. Genet. 51, 431–439. doi: 10.1134/S1022795415050063
Ribeiro, T. C., Weiblen, C., de Avila Botton, S., Pereira, D. I. B., de Jesus, F. P. K., Verdi, C. M., et al. (2016). In vitro susceptibility of the oomycete Pythium insidiosum to metallic compounds containing cadmium, lead, copper, manganese or zinc. Med. Mycol. doi: 10.1093/mmy/myw115 [Epub ahead of print].
Roberts, D. W. (2012). labdsv: Ordination and Multivariate Analysis for Ecology. Package Version 1.6-1. Available at: http://cran.r-project.org/web/packages/labdsv/labdsv.pdf
Rondon, M. R., August, P. R., Bettermann, A. D., Brady, S. F., Grossman, T. H., Liles, M. R., et al. (2000). Cloning the soil metagenome: a strategy for accessing the genetic and functional diversity of uncultured microorganisms. Appl. Environ. Microbiol. 66, 2541–2547. doi: 10.1128/AEM.66.6.2541-2547.2000
Sapkota, R., and Nicolaisen, M. (2015). An improved high throughput sequencing method for studying oomycete communities. J. Microbiol. Methods 110, 33–39. doi: 10.1016/j.mimet.2015.01.013
Soares, M. A., Li, H., Kowalski, K. P., Bergen, M., Torres, M. S., and White, J. F. (2016). Evaluation of the functional roles of fungal endophytes of Phragmites australis from high saline and low saline habitats. Biol. Invasions 18, 2689–2702. doi: 10.1007/s10530-016-1160-z
Thines, M. (2014). Phylogeny and evolution of plant pathogenic oomycetes-a global overview. Eur. J. Plant Pathol. 138, 431–447. doi: 10.1007/s10658-013-0366-5
Toda, T., Iwasa, A., Shinichi, F., and Furuya, H. (2015). Widespread occurrence of Pythium arrhenomanes pathogenic to rice seedlings around Japanese rice fields. Plant Dis. 99, 1823–1831. doi: 10.1094/PDIS-01-15-0124-RE
van der Putten, W. H. (1997). Die-back of Phragmites australis in European wetlands: an overview of the European research programme on reed die-back and progression (1993–1994). Aquat. Bot. 59, 263–275. doi: 10.1016/S0304-3770(97)00060-0
Vandenkoornhuyse, P., Quaiser, A., Duhamel, M., Le Van, A., and Dufresne, A. (2015). The importance of the microbiome of the plant holobiont. New Phytol. 206, 1196–1206. doi: 10.1111/nph.13312
Vettraino, A. M., Bonats, P., Tomassini, A., Bruni, N., and Vannini, A. (2012). Pyrosequencing as a tool for detection of Phytophthora species: error rate and risk of false MOTUs. Lett. Appl. Microbiol. 55, 390–396. doi: 10.1111/j.1472-765x.2012.03310.x
Wielgoss, A., Nechwatal, J., Bogs, C., and Mendgen, K. (2009). Host plant development, water level and water parameters shape Phragmites australis-associated oomycete communities and determine reed pathogen dynamics in a large lake. FEMS Microbiol. Ecol. 69, 255–265. doi: 10.1111/j.1574-6941.2009.00701.x
Wirsel, S. G. R., Leibinger, W., Ernst, M., and Mendgen, K. (2002). Genetic diversity of fungi closely associated with common reed. New Phytol. 149, 589–598. doi: 10.1046/j.1469-8137.2001.00038.x
Xu, L., Ravnskov, S., Larsen, J., Nilsson, R. H., and Nicolaisen, M. (2012). Soil fungal community structure along a soil health gradient in pea fields examined using deep amplicon sequencing. Soil Biol. Biochem. 46, 26–32. doi: 10.1016/j.soilbio.2011.11.010
Keywords: common reed, freshwater ecosystem, metabarcoding, microbial ecology, Phragmites australis, rhizosphere, wetlands
Citation: Cerri M, Sapkota R, Coppi A, Ferri V, Foggi B, Gigante D, Lastrucci L, Selvaggi R, Venanzoni R, Nicolaisen M, Ferranti F and Reale L (2017) Oomycete Communities Associated with Reed Die-Back Syndrome. Front. Plant Sci. 8:1550. doi: 10.3389/fpls.2017.01550
Received: 17 May 2017; Accepted: 24 August 2017;
Published: 07 September 2017.
Edited by:
Davide Bulgarelli, University of Dundee, United KingdomReviewed by:
Ida Karlsson, Swedish University of Agricultural Sciences, SwedenMatthew G. Bakker, Agricultural Research Service (USDA), United States
Copyright © 2017 Cerri, Sapkota, Coppi, Ferri, Foggi, Gigante, Lastrucci, Selvaggi, Venanzoni, Nicolaisen, Ferranti and Reale. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martina Cerri, Y2VycmkubWFydGluYUBnbWFpbC5jb20=