Hirokazu Toju
Hiroko Kurokawa
Tanaka Kenta- 1Center for Ecological Research, Kyoto University, Kyoto, Japan
- 2Precursory Research for Embryonic Science and Technology, Japan Science and Technology Agency, Kawaguchi, Japan
- 3Forestry and Forest Products Research Institute, Tsukuba, Japan
- 4Sugadaira Research Station, Mountain Science Center, University of Tsukuba, Nagano, Japan
In terrestrial ecosystems, plants interact with diverse taxonomic groups of bacteria and fungi in the phyllosphere and rhizosphere. Although recent studies based on high-throughput DNA sequencing have drastically increased our understanding of plant-associated microbiomes, we still have limited knowledge of how plant species in a species-rich community differ in their leaf and root microbiome compositions. In a cool-temperate semi-natural grassland in Japan, we compared leaf- and root-associated microbiomes across 137 plant species belonging to 33 plant orders. Based on the whole-microbiome inventory data, we analyzed how sampling season as well as the taxonomy, nativeness (native or alien), lifeform (herbaceous or woody), and mycorrhizal type of host plants could contribute to variation in microbiome compositions among co-occurring plant species. The data also allowed us to explore prokaryote and fungal lineages showing preferences for specific host characteristics. The list of microbial taxa showing significant host preferences involved those potentially having some impacts on survival, growth, or environmental resistance of host plants. Overall, this study provides a platform for understanding how plant and microbial communities are linked with each other at the ecosystem level.
Introduction
Plants interact with various taxonomic groups of microbes both in the phyllosphere and rhizosphere (van der Heijden et al., 1998; Berendsen et al., 2012; Bai et al., 2015; Peay et al., 2016). Diverse bacteria and yeasts, for example, are present on leaf surfaces, involved in underappreciated metabolic pathways (Mercier and Lindow, 2000; Delmotte et al., 2009; Hacquard et al., 2015). In addition to those epiphytes, a number of bacteria and filamentous fungi are known to inhabit leaf tissue (Ding and Melcher, 2016; Hamonts et al., 2018), playing pivotal roles in resistance of host plants against biotic and abiotic environmental stresses (Schardl and Phillips, 1997; Arnold et al., 2003; Hardoim et al., 2015; Santoyo et al., 2016). In root systems, mycorrhizal fungi provide plants with soil phosphorus and/or nitrogen, fueling hosts’ growth (Parniske, 2008; Smith and Read, 2008; Tedersoo et al., 2010). Likewise, some endophytic fungal taxa have been known to enhance nutritional conditions of host plants (Newsham, 2011; Hiruma et al., 2016; Almario et al., 2017). Moreover, rhizosphere/endophytic bacteria and fungi associated with roots can increase disease resistance of host plants, possibly by stimulating host immune systems (Ramamoorthy et al., 2001; Weston et al., 2012; Pieterse et al., 2014; Hacquard et al., 2017) or by suppressing soil pathogens directly (Compant et al., 2005; Gao et al., 2010; Durán et al., 2018; Kwak et al., 2018). Thus, understanding of the compositions of plant microbiomes is a prerequisite for understanding the physiology and ecology of plants in terrestrial ecosystems (van der Heijden et al., 2008; Schlaeppi and Bulgarelli, 2015; Toju et al., 2018a).
While exploration of plant microbiomes has been accelerated since the emergence of high-throughput DNA sequencing (Öpik et al., 2009; Lundberg et al., 2012; Bai et al., 2015), we still have limited knowledge of how diverse plant species co-occurring in a grassland or forest ecosystem can differ in their microbiome compositions (Toju et al., 2016a). Moreover, most plant microbiome studies target only bacteria or fungi [but see (Agler et al., 2016)) in either above- or below-ground systems [but see (Bai et al., 2015; Wagner et al., 2016)], precluding comprehensive understanding of microbiome compositions. Given that bacteria and fungi can interact with each other within hosts (Frey-Klett et al., 2007; Hoffman and Arnold, 2010) and that above- and below-ground ecological processes can be interlinked (Bever et al., 2010; Mangan et al., 2010; van der Putten et al., 2013), the targets of plant microbiome studies need to be expanded toward a better understanding of the processes that organize plant and microbial communities in the wild. Studies comparing microbiome compositions across tens (or more) of plant species co-occurring in natural ecosystems (Toju et al., 2014, 2018b), in particular, will allow us to examine what kinds of host properties can contribute to the organization of leaf- and root-associated microbial communities.
In this study, we sampled leaves and roots of 137 plant species representing 111 genera, 55 families, and 33 orders in a cool-temperate grassland in Japan, thereby performing a high-throughput sequencing analysis of both prokaryote and fungal communities associated with plants. The sample set of diverse plant species allowed us to examine what host properties can contribute to variation in leaf and root microbiome compositions in an ecosystem. Furthermore, we statistically tested how each prokaryote or fungal genus showed preferences for seasons as well as preferences for nativeness (native or alien), lifeform (herbaceous or woody), and mycorrhizal type (ectomycorrhizal, arbuscular mycorrhizal, non-mycorrhizal, or variable mycorrhizal) of host plants. Overall, this study, for the first time, shows how more than 100 plant species in a single ecosystem can differ in their leaf and root microbiome compositions depending on their characteristics. The statistical results on plant–microbe associations shed light on underappreciated diversity of host–symbiont associations in grasslands, providing fundamental information for conserving and restoring terrestrial ecosystems.
Materials and Methods
Sampling
Fieldwork was conducted in Sugadaira Research Station, Mountain Science Center, University of Tsukuba, Sugadaira, Ueda, Nagano Prefecture, Japan (36.524°N; 138.349°E; 1340 m asl) (Figure 1). In Sugadaira Research Station, 6 ha of a semi-natural grassland has been maintained by mowing plants in autumn and thereby preventing the community succession to a forest. Thus, woody plant species that occurred in the grassland are shrubs or saplings of tall trees colonized from surrounding forests. In total, 200 plant species have been observed from the grassland, including some alien species non-native to the Japanese Archipelago [plant species uncited in Ohashi et al. (2016)].
Figure 1. Grassland of Sugadaira Research Station.
In the grassland, both native and alien plant species were sampled to reveal the compositions of prokaryote and fungal communities associated with leaves and roots through summer and autumn (July 19–20, August 16–18, and September 7–8) in 2017. We targeted only non-reproductive plant individuals that had neither flower buds, flowers, nor fruits so that plant physiology and chemistry would not be affected by reproduction [sensu Obeso (2002)]. We tried to sample as many plant species as possible within the sampling days in each month. Note that root systems of multiple plant species were tangled with each other at the study site due to the dominance of perennial plants. Therefore, we sampled 1–8 liters of soil (A horizon) including root systems for each target plant individuals and quickly washed the root system in a nearby laboratory to carefully trace root tips directly connected to above-ground tissue of the target plant. A 1-cm2 disk of a mature leaf and a 2-cm segment of a terminal root were collected from each plant sample and preserved at -20°C until DNA extraction. After the sampling, remaining plant organs of rare plant species were replaced at the original sampling positions. In total, 289 plant individuals representing 138 plant species (112 genera, 55 families, 33 orders) were collected (Supplementary Data S1). The sampling points were scattered across the 6-ha semi-natural grassland of Sugadaira Research Station: sampling in August and September was conducted in the spatial positions that had not been disturbed by sampling in previous months. The fieldwork and research were permitted by Sugadaira Research Station, Mountain Science Center, University of Tsukuba.
DNA Extraction, PCR, and Sequencing
Each leaf or root sample was surface-sterilized by immersing it in × 1/100 NaClO (Nacalai Tesque) for 1 min and it was subsequently washed in ethanol twice. DNA was extracted with a cetyltrimethylammonium bromide (CTAB) method after pulverizing the roots with 4 mm zirconium balls at 25 Hz for 3 min using a TissueLyser II (Qiagen) (Toju et al., 2013).
For each of the leaf and root samples, the 16S rRNA V4 region of the prokaryotes and the internal transcribed spacer 1 (ITS1) region of fungi were amplified. The PCR of the 16S rRNA region was performed with the forward primer 515f (Caporaso et al., 2011) fused with 3–6-mer Ns for improved Illumina sequencing quality (Lundberg et al., 2013) and the forward Illumina sequencing primer (5′- TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG- [3–6-mer Ns] – [515f] -3′) and the reverse primer 806rB (Apprill et al., 2015) fused with 3–6-mer Ns and the reverse sequencing primer (5′- GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA G [3–6-mer Ns] - [806rB] -3′) (0.2 μM each). To prevent the amplification of mitochondrial and chloroplast 16S rRNA sequences of host plants, specific peptide nucleic acids [mPNA and pPNA; Lundberg et al. (2013)] (0.7 μM each) were added to the reaction mix of KOD FX Neo (Toyobo). The temperature profile of the PCR was 94°C for 2 min, followed by 35 cycles at 98°C (denaturation) for 10 s, 78°C (annealing of PNA) for 10 s, 60°C (annealing of primers) for 30 s, and 68°C (extension) for 50 s, and a final extension at 68°C for 5 min. To prevent generation of chimeric sequences, the ramp rate through the thermal cycles was set to 1°C/sec (Stevens et al., 2013). Illumina sequencing adaptors were then added to respective samples in the supplemental PCR using the forward fusion primers consisting of the P5 Illumina adaptor, 8-mer indexes for sample identification (Hamady et al., 2008) and a partial sequence of the sequencing primer (5′- AAT GAT ACG GCG ACC ACC GAG ATC TAC AC - [8-mer index] - TCG TCG GCA GCG TC -3′) and the reverse fusion primers consisting of the P7 adaptor, 8-mer indexes, and a partial sequence of the sequencing primer (5′- CAA GCA GAA GAC GGC ATA CGA GAT - [8-mer index] - GTC TCG TGG GCT CGG -3′). KOD FX Neo was used with a temperature profile of 94°C for 2 min, followed by 8 cycles at 98°C for 10 s, 55°C for 30 s, and 68°C for 50 s (ramp rate = 1°C/s), and a final extension at 68°C for 5 min. The PCR amplicons of the samples were then pooled after a purification/equalization process with the AMPureXP Kit (Beckman Coulter). Primer dimers, which were shorter than 200 bp, were removed from the pooled library by supplemental purification with AMpureXP: the ratio of AMPureXP reagent to the pooled library was set to 0.6 (v/v) in this process.
The PCR of fungal ITS1 region was performed with the forward primer ITS1F_KYO1 (Toju et al., 2012) fused with 3–6-mer Ns for improved Illumina sequencing quality (Lundberg et al., 2013) and the forward Illumina sequencing primer (5′- TCG TCG GCA GCG TCA GAT GTG TAT AAG AGA CAG- [3–6-mer Ns] – [ITS1F_KYO1] -3′) and the reverse primer ITS2_KYO2 (Toju et al., 2012) fused with 3–6-mer Ns and the reverse sequencing primer (5′- GTC TCG TGG GCT CGG AGA TGT GTA TAA GAG ACA G [3–6-mer Ns] - [ITS2_KYO2] -3′). The buffer and polymerase system of KOD FX Neo was used with a temperature profile of 94°C for 2 min, followed by 35 cycles at 98°C for 10 s, 58°C for 30 s, and 68°C for 50 s, and a final extension at 68°C for 5 min. Illumina sequencing adaptors and 8-mer index sequences were then added in the second PCR as described above. The amplicons were purified and pooled as described above.
The sequencing libraries of the prokaryote 16S and fungal ITS regions were processed in an Illumina MiSeq sequencer (run center: KYOTO-HE; 15% PhiX spike-in). Because the quality of forward sequences is generally higher than that of reverse sequences in Illumina sequencing, we optimized the MiSeq run setting in order to use only forward sequences. Specifically, the run length was set 271 forward (R1) and 31 reverse (R4) cycles to enhance forward sequencing data: the reverse sequences were used only for discriminating between 16S and ITS1 sequences based on the sequences of primer positions.
Bioinformatics
The raw sequencing data were converted into FASTQ files using the program bcl2fastq 1.8.4 distributed by Illumina. The output FASTQ files were demultiplexed with the program Claident v0.2. 2018.05.29 (Tanabe and Toju, 2013; Tanabe, 2018), by which sequencing reads whose 8-mer index positions included nucleotides with low ( < 30) quality scores were removed. Only forward sequences were used in the following analyses after removing low-quality 3′-ends using Claident. Noisy reads (Tanabe, 2018) were subsequently discarded and then denoised dataset consisting of 2,973,811 16S and 2,774,197 ITS1 reads were obtained. The sequencing data were deposited to DNA Data Bank of Japan (DDBJ1) (Accession No.: DRA007062).
For each dataset of 16S and ITS1 regions, filtered reads were clustered with a cut-off sequencing similarity of 97% using the program VSEARCH (Rognes et al., 2014) as implemented in Claident. The operational taxonomic units (OTUs) representing less than 10 sequencing reads were subsequently discarded (Supplementary Data S2). The molecular identification of the remaining OTUs was performed based on the combination of the query-centric auto-k-nearest neighbor (QCauto) method (Tanabe and Toju, 2013) and the lowest common ancestor (LCA) algorithm (Huson et al., 2007) as implemented in Claident (Supplementary Data S2). Note that taxonomic identification results based on the combination of the QCauto search and the LCA taxonomic assignment are comparable to, or sometimes more accurate than, those with alternative approaches (Tanabe and Toju, 2013; Toju et al., 2016a; Toju et al., 2016b).
For each combination of target region (16S or ITS1) and sample type (root or soil), we obtained a sample × OTU matrix, in which a cell entry depicted the number of sequencing reads of an OTU in a sample (Supplementary Data S3). The cell entries whose read counts represented less than 0.1% of the total read count of each sample were removed to minimize effects of PCR/sequencing errors (Peay et al., 2015). The filtered matrix was then rarefied to 500 reads per sample using the “rrarefy” function of the vegan 2.5-2 package (Oksanen et al., 2012) of R 3.5.1 (R-Core-Team, 2018). Samples with less than 500 reads were discarded in this process: the numbers of OTUs in the rarefied sample × OTU matrices were 1,470, 5,638, 1,537, and 3,367 for leaf prokaryote, root prokaryote, leaf fungal, and root fungal datasets, respectively (Supplementary Data S4). For each dataset, we also obtained order- and genus-level matrices, which represented order- and genus-level taxonomic compositions of microbes (prokaryotes or fungi), respectively (Supplementary Data S5).
Prokaryote and Fungal Diversity
Relationships between the number of sequencing reads and that of detected OTUs were examined for respective data matrices (leaf prokaryote, root prokaryote, leaf fungal, and root fungal datasets) with the “rarecurve” function of the R vegan package. Likewise, relationships between the number of samples and that of prokaryote/fungal orders or genera were examined with the vegan “specaccum” function. The order-level taxonomic compositions of leaf prokaryotes, root prokaryotes, leaf fungi, and root fungi were visualized in bar graphs for respective plant orders.
Factors Contributing to Microbiome Compositions
For each dataset (leaf prokaryote, root prokaryote, leaf fungal, or root fungal dataset), factors contributing to microbial community compositions were examined with the permutational analysis of variance [PERMANOVA; Anderson (2001)] using the vegan “adonis” function (10,000 permutations). Sampling month (July, August, or September) and four variables representing host plant properties were included as explanatory variables. Specifically, order-level plant taxonomy, plant nativeness (native or alien) [based on the list of plant species native to the Japanese Archipelago Ohashi et al. (2016)], plant lifeform (herbaceous or woody), and plant mycorrhizal type [ectomycorrhizal (EcM), arbuscular mycorrhizal (AM), non-mycorrhizal (NM), or variable mycorrhizal (NM-AM)] (Brundrett, 2009) were included as variables representing host properties. In each model, a matrix representing order- or genus-level taxonomic compositions of prokaryotes/fungi was used as the input response matrix. The “Bray-Curtis” metric of β-diversity was used in the PERMANOVA analyses.
Preferences of Each Prokaryote/Fungal Genus
To explore prokaryote/fungal genera that preferentially occurred on plant samples with specific properties, a series of linear regression analyses were performed. For each prokaryote/fungal genus that appeared in 30 or more samples in a genus-level matrix (the leaf prokaryote, root prokaryote, leaf fungal, or root fungal genus-level matrix), a linear regression model of relative abundance (values in a genus-level matrix) was constructed by incorporating a sample property as an explanatory variable. In each model, sampling month (July, August, or September), plant nativeness (native or alien), plant lifeform (herbaceous or woody), or plant mycorrhizal type was examined. By z-transforming response variables (i.e., zero-mean and unit-variance), a standardized coefficient was obtained for each combination of a genus and a sample property. Because most plant orders included a few plant species in our datasets, the regression analysis was not applied to plant taxonomy.
Results
Prokaryote and Fungal Diversity
After a series of quality filtering and rarefaction procedures, 41.1 (SD = 22.1), 143.4 (SD = 37.9), 54.5 (SD = 18.8), and 46.0 (SD = 22.5) OTUs per sample, on average, were detected, from the leaf prokaryote, root prokaryote, leaf fungal, and root fungal datasets, respectively (Supplementary Figure S1). The number of the samples from which 500 or more reads of prokaryote/fungal sequences were obtained varied among datasets (Figure 2): high proportions of host 16S rRNA and ITS sequences swamped leaf datasets in various taxonomic groups of plants, presumably due to higher concentrations of host organelle/nuclear DNA in leaves than in roots. In total, microbiome data of any of the leaf prokaryote, root prokaryote, leaf fungal, or root fungal community were obtained from 284 plant individuals representing 137 plant species, 111 genera, 55 families, and 33 orders. The numbers of prokaryote orders and genera were higher in root samples than in leaf samples, while those of fungal orders and genera showed opposite patterns (Figure 2).
Figure 2. Relationship between the number of leaf/root samples and that of prokaryote/fungal taxa observed. (A) Number of prokaryote orders. In total, sequencing data were successfully obtained from 188 leaf and 270 root samples. Blue and red curves represent leaf and root samples, respectively. (B) Number of fungal orders. In total, sequencing data were successfully obtained from 213 leaf and 261 root samples. (C) Number of prokaryote genera (188 leaf and 270 root samples). (D) Number of fungal genera (213 leaf and 261 root samples).
The leaf prokaryote communities of the examined plants were dominated by the order Rhizobiales, while diverse bacterial taxa constituted the root prokaryote communities (Figures 3A,B). In the leaf fungal communities, the order Capnodiales were the most abundant, while root fungal community compositions varied considerably among host plant orders (Figures 3C,D).
Figure 3. Prokaryote and fungal community compositions. (A) Order-level compositions of prokaryotes in leaf samples. Mean proportions of sequencing reads are shown for each plant order. In total, sequencing data were successfully obtained from 188 leaf samples. (B) Order-level compositions of prokaryotes in root samples (270 root samples). (C) Order-level compositions of fungi in leaf samples (213 leaf samples). Mean proportions of sequencing reads are shown for each plant order. (D) Order-level compositions of fungi in root samples (261 root samples).
Factors Contributing to Microbiome Compositions
In the PERMANOVA, sampling month had significant effects on the leaf prokaryote, root prokaryote, and leaf fungal community compositions but not on the root fungal community structure (Table 1 and Supplementary Table S1). Meanwhile, order-level host taxonomy influenced the root prokaryote, leaf fungal, and root fungal community compositions but not the leaf prokaryote community structure (Table 1). The nativeness of host plants (native or alien) had significant impacts on the root prokaryote and the root fungal (genus-level) community compositions (Table 1). The analysis also showed that host plant lifeform (herbaceous or woody) had significant effects on the leaf fungal community structure (Table 1). In contrast, mycorrhizal type of host plants had significant effects on none of the communities examined in the PERMANOVA (Table 1).
Table 1. Factors contributing to variation in genus-level community compositions of bacteria and fungi.
Preferences of Each Prokaryote/Fungal Genus
In the linear regression analyses, the relative abundances of a bacterial genus and eight fungal genera changed through the sampling months in the leaf sample data (Figures 4A,B). For example, the fungal genera Leucosporidium, Taphrina, and Dioszegia in the leaf fungal community appeared preferentially in July, while the bacterial genus Amnibacterium preferentially occurred in September (Figure 4). On the other hand, no bacterial/fungal genera showed preferences for sampling months in the root microbiome data (Figures 4C,D). Regarding the nativeness of hosts, bacteria in the genera Rhodococcus, Mucilaginibacter, Deinococcus, and Pseudomonas in the leaf microbiome data showed statistically significant preferences for alien plant species, while two leaf-associated bacterial genera, Actinoallomurus and Singulisphaera, showed preferences for native plant species (Figures 5A,B). No fungal genus showed significant preferences for the nativeness of host plants (Figures 5C,D). Although mycorrhizal type of host plants did not have significant effects in the community-level statistical analysis (Table 1), nine bacterial and two fungal genera showed preferences for host mycorrhizal type (Figure 6). For example, the bacterial genera Nocardioides and Pseudonocardia preferentially occurred in leaves of non-mycorrhizal plants (Figure 6A). Likewise, three bacterial genera, Ferrimicrobium, Mycobacterium, and Nocardioides, preferentially appeared in the roots of non-mycorrhizal plants (Figure 6B). Meanwhile, Flavisolibacter and Phenylobacterium showed preferences for ectomycorrhizal plant roots, while Rubrobacter and Gemmata displayed preferences for variable mycorrhizal (NM-AM) plants (Figure 6B). In the fungal community, Nigrospora in the leaf data and Phialocephala in the root data occurred preferentially in ectomycorrhizal plant samples (Figures 6C,D). Among the prokaryote and fungal genera examined, none showed statistically significant preferences for the lifeform of host plants (Supplementary Figure S2).
Figure 4. Effects of sampling month on the relative abundances of respective prokaryote and fungal genera. For each prokaryote/fungal genus that appeared in 30 or more samples, a linear model of z-standardized (zero-mean and unit-variance) relative abundance was constructed by incorporating sampling month as an explanatory variable. A standardized coefficient was then obtained for each genus × month combination. The P-values were converted to false discovery rates (FDRs) by taking into account the number of the genera and months examined. (A) Leaf prokaryotes. (B) Root prokaryotes. (C) Leaf fungi. (D) Root fungi.
Figure 5. Effects of sampling host plant nativeness on the relative abundances of respective prokaryote and fungal genera. For each prokaryote/fungal genus that appeared in 30 or more samples, a linear model of z-standardized (zero-mean and unit-variance) relative abundance was constructed by incorporating host plant nativeness (native/alien) as an explanatory variable. (A) Leaf prokaryotes. (B) Root prokaryotes. (C) Leaf fungi. (D) Root fungi.
Figure 6. Effects of host-plant mycorrhizal type on the relative abundances of respective prokaryote and fungal genera. For each prokaryote/fungal genus that appeared in 30 or more samples, a linear model of z-standardized (zero-mean and unit-variance) relative abundance was constructed by incorporating mycorrhizal type of host plants as an explanatory variable. EcM, ectomycorrhizal; AM, arbuscular mycorrhizal; NM, non-mycorrhizal; NM-AM, variable mycorrhizal (i.e., non-mycorrhizal or arbuscular mycorrhizal). (A) Leaf prokaryotes. (B) Root prokaryotes. (C) Leaf fungi. (D) Root fungi.
Discussion
Based on a high-throughput sequencing dataset, we herein compared leaf and root microbiome compositions across co-occurring plant species in an ecosystem. By targeting a plant-species-rich grassland in the cool-temperate climate, we compared leaf- and root-associated microbial communities across 33 plant orders (Figure 3) and then performed a series of statistical analyses on factors that may influence community compositions of plant-associated microbes (Figures 4–6 and Table 1). Hereafter, we discuss potential contributions of the factors examined, focusing on preferences of each microbial taxon for host characteristics.
An interesting finding of this study is that, while the compositions of leaf prokaryote, root prokaryote, and leaf fungal communities changed through the sampling months [cf. (Davison et al., 2012; Copeland et al., 2015)], root fungal community compositions did not significantly shift during the period (Table 1). This pattern possibly represents difference in basic environmental features between above- and below-ground systems and/or difference in phenological patterns between prokaryote and fungal communities. For example, above-ground biotic/abiotic environments may be more dynamic than below-ground environments, resulting in rapid turnover of microbial communities. Moreover, above-ground parts of plants are more likely to be accessed by wind-dispersed spores and inocula than below-ground parts (Berg et al., 2014): hence, above-ground microbiome processes may be susceptible to continual immigration. In addition to potential contrasting features of above- vs. below-ground systems, difference in basic ecology between bacteria and fungi may have contributed to the varied phenological patterns. While mycorrhizal and endophytic fungi usually persist on/around host root systems in the form of hyphal networks (Lian et al., 2006; Smith and Read, 2008; Hiruma et al., 2016; Almario et al., 2017), bacterial communities may consist mainly of opportunistic symbionts [sensu (Hardoim et al., 2008)], which undergo rapid population growth under favorable environmental conditions and are subsequently replaced by others. Year-round comparative studies on leaf and root microbiomes are awaited for gaining more comprehensive understandings of microbiome dynamics.
Among the microbial communities examined, both root-associated prokaryote and fungal communities significantly varied between native and alien plant species (Table 1). Meanwhile, a series of regression analyses targeting respective microbial genera indicated that some bacterial genera occurred preferentially in the leaves or roots of sampled plants (Figures 5A,B). Among the genera showing preferences for alien plant leaves, Rhodococcus is known to involve plant pathogens (Vereecke et al., 2000), while Mucilaginibacter includes species with xylan- and pectin-degrading abilities (Pankratov et al., 2007) (Figure 5A). Singulisphaera, which showed preferences for native plant roots (Figure 5B), is a genus of bacteria reported from cold, acidic environments (Kulichevskaya et al., 2008). The analysis also showed that some genera in the phylum Actinobacteria (Rhodococcus and Actinoallomurus) showed preferences for native or alien plant species. Given that many actinomycete bacteria produce chemicals suppressing other microbes (Qin et al., 2011; Bérdy, 2012), their ecological roles in ecosystems are of particular interest. Overall, these results suggest that various taxonomic groups of bacteria are preferentially associated with native or alien plant species, potentially affecting invasiveness of alien plants both positively and negatively (Mitchell and Power, 2003; Moora et al., 2011; Majewska et al., 2015).
At the whole community level, mycorrhizal types of host plants did not have significant effects on plant microbiome compositions (Table 1), unlike previous studies comparing root-associated microbial communities between co-occurring plants with different mycorrhizal types (Toju et al., 2016a; Toju and Sato, 2018). Nonetheless, the analyses for respective microbial taxa highlighted diverse bacterial and fungal genera showing statistically significant preferences for host mycorrhizal types (Figure 6). Among the bacteria showing preferences for non-mycorrhizal plants, Ferrimicrobium includes species adapted to low pH conditions (Johnson et al., 2009), while the genus Nocardioides consists of bacteria described as endophytes (Qin et al., 2012; Han et al., 2013). Nigrospora, which occurred preferentially in ectomycorrhizal plant leaves (Figure 6C), is known to include fungi producing chemicals with antiviral and antifungal functions (Kim et al., 2001; He et al., 2012). We also found that possibly endophytic fungi in the genus Phialocephala (Fernando and Currah, 1996; Grünig et al., 2008) showed preferences for ectomycorrhizal plants. Thus, the list of microbes preferentially associated with plants with specific mycorrhizal types (Figure 6) sheds new light on potential diversity of bacteria and non-mycorrhizal fungi, whose physiological and ecosystem-scale functions remain to be investigated (Hardoim et al., 2015; Peay et al., 2016).
Although the data collected in this study provide fundamental information of microbial diversity in a grassland ecosystem, the statistical results should be interpreted with caution. For example, it should be acknowledged that the small number of samples per plant species may have affected the comparison of microbiome compositions among plant taxa (Figure 3). The identification of plant roots is time-consuming especially in species-rich grasslands consisting mainly of perennial plants with tangled root systems, limiting the throughput of sampling. Therefore, for more comprehensive profiling of plant microbiomes, we may need to increase the throughput of plant species identification based on molecular taxonomic assignment (i.e., DNA barcoding) of host plants (Hollingsworth et al., 2009; Toju et al., 2013). Another potential pitfall is that the presence of unidentified bacteria and fungi in the dataset may have biased the statistical analyses. Although databases of microbes have been continually updated (Abarenkov et al., 2010), there remain many bacterial and fungal lineages whose taxonomy has not yet been fixed. In particular, below-ground microbiomes are known to involve a number of poorly investigated taxa, whose physiological and ecological functions remain to be uncovered (Buée et al., 2009; Fierer, 2017). Thus, with more reference microbial databases (Langille et al., 2013; Nguyen et al., 2016), we will be able to examine whether the patterns found in the present analysis hold after assigning unidentified OTUs to right categories.
We herein revealed how diverse bacterial and fungal taxa were associated with leaves and roots of the 138 plant species co-occurring in a cool-temperate grassland, focusing on potential contributions of host plant characteristics on microbiome compositions. Although recent ecological studies have highlighted possible feedbacks between plant and microbial community dynamics (Bever et al., 2010; Mangan et al., 2010; van der Putten et al., 2013), we still have limited knowledge of the processes by which species-rich plant communities are maintained by phyllosphere and rhizosphere microbiomes. Accumulating comprehensive inventory data of microbiomes associated with whole plant communities is a prerequisite for advancing our understanding of ecosystem-scale processes. Case studies in various types of terrestrial ecosystems in diverse climatic regions will allow us to elucidate how plant species with different mycorrhizal types often coexist in natural ecosystems (Booth, 2004; Kadowaki et al., 2018) or why some ecosystems are resistant against alien plants, while others are heavily disturbed by invasive species (Mitchell and Power, 2003; Reinhart and Callaway, 2006).
Author Contributions
HT, HK, and TK designed the work, performed the fieldwork, and wrote the manuscript. HT conducted molecular experiments and analyzed the data.
Funding
This work was financially supported by JST PRESTO (JPMJPR16Q6) to HT, by the JSPS Grant-in-Aid for Scientific Research (B) Grant number JP17H03736 to HK, and by the promotion cost for functional enhancement of the Mountain Science Center to TK.
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.
Acknowledgments
We thank Qing-Wei Wang, Yuka Yamamoto, and Rie Watanabe for their help in fieldwork and Sarasa Amma and Hiroki Kawai for their support in molecular experiments. We are also grateful to reviewers for their constructive comments on the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00241/full#supplementary-material
Footnotes
References
Abarenkov, K., Henrik Nilsson, R., Larsson, K. H., Alexander, I. J., Eberhardt, U., Erland, S., et al. (2010). The UNITE database for molecular identification of fungi–recent updates and future perspectives. New Phytol. 186, 281–285. doi: 10.1111/j.1469-8137.2009.03160.x
Agler, M. T., Ruhe, J., Kroll, S., Morhenn, C., Kim, S.-T., Weigel, D., et al. (2016). Microbial hub taxa link host and abiotic factors to plant microbiome variation. PLoS Biol. 14:e1002352. doi: 10.1371/journal.pbio.1002352
Almario, J., Jeena, G., Wunder, J., Langen, G., Zuccaro, A., Coupland, G., et al. (2017). Root-associated fungal microbiota of nonmycorrhizal Arabis alpina and its contribution to plant phosphorus nutrition. Proc. Natl. Acad. Sci. U.S.A. 114, E9403–E9412. doi: 10.1073/pnas.1710455114
Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. doi: 10.1111/j.1442-9993.2001.01070.pp.x
Apprill, A., Mcnally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 75, 129–137. doi: 10.3354/ame01753
Arnold, A. E., Mejía, L. C., Kyllo, D., Rojas, E. I., Maynard, Z., Robbins, N., et al. (2003). Fungal endophytes limit pathogen damage in a tropical tree. Proc. Natl. Acad. Sci. U.S.A. 100, 15649–15654. doi: 10.1073/pnas.2533483100
Bai, Y., Müller, D. B., Srinivas, G., Garrido-Oter, R., Potthoff, E., Rott, M., et al. (2015). Functional overlap of the Arabidopsis leaf and root microbiota. Nature 528, 364–369. doi: 10.1038/nature16192
Bérdy, J. (2012). Thoughts and facts about antibiotics: where we are now and where we are heading. J. Antibiot. 65, 385–395. doi: 10.1038/ja.2012.27
Berendsen, R. L., Pieterse, C. M., and Bakker, P. A. (2012). The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. doi: 10.1016/j.tplants.2012.04.001
Berg, G., Grube, M., Schloter, M., and Smalla, K. (2014). Unraveling the plant microbiome: looking back and future perspectives. Front. Microbiol. 5:148. doi: 10.3389/fmicb.2014.00148
Bever, J. D., Dickie, I. A., Facelli, E., Facelli, J. M., Klironomos, J., Moora, M., et al. (2010). Rooting theories of plant community ecology in microbial interactions. Trends Ecol. Evol. 25, 468–478. doi: 10.1016/j.tree.2010.05.004
Booth, M. G. (2004). Mycorrhizal networks mediate overstorey-understorey competition in a temperate forest. Ecol. Lett. 7, 538–546. doi: 10.1111/j.1461-0248.2004.00605.x
Brundrett, M. C. (2009). Mycorrhizal associations and other means of nutrition of vascular plants: understanding the global diversity of host plants by resolving conflicting information and developing reliable means of diagnosis. Plant Soil 320, 37–77. doi: 10.1007/s11104-008-9877-9
Buée, M., Reich, M., Murat, C., Morin, E., Nilsson, R., Uroz, S., et al. (2009). 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol. 184, 449–456. doi: 10.1111/j.1469-8137.2009.03003.x
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Compant, S., Duffy, B., Nowak, J., Clément, C., and Barka, E. A. (2005). Use of plant growth-promoting bacteria for biocontrol of plant diseases: principles, mechanisms of action, and future prospects. Appl. Environ. Microbiol. 71, 4951–4959. doi: 10.1128/AEM.71.9.4951-4959.2005
Copeland, J. K., Yuan, L., Layeghifard, M., Wang, P. W., and Guttman, D. S. (2015). Seasonal community succession of the phyllosphere microbiome. Mol. Plant Microbe Interact. 28, 274–285. doi: 10.1094/MPMI-10-14-0331-FI
Davison, J., Öpik, M., Zobel, M., Vasar, M., Metsis, M., and Moora, M. (2012). Communities of arbuscular mycorrhizal fungi detected in forest soil are spatially heterogeneous but do not vary throughout the growing season. PLoS One 7:e41938. doi: 10.1371/journal.pone.0041938
Delmotte, N., Knief, C., Chaffron, S., Innerebner, G., Roschitzki, B., Schlapbach, R., et al. (2009). Community proteogenomics reveals insights into the physiology of phyllosphere bacteria. Proc. Natl. Acad. Sci. U.S.A. 106, 16428–16433. doi: 10.1073/pnas.0905240106
Ding, T., and Melcher, U. (2016). Influences of plant species, season and location on leaf endophytic bacterial communities of non-cultivated plants. PLoS One 11:e0150895. doi: 10.1371/journal.pone.0150895
Durán, P., Thiergart, T., Garrido-Oter, R., Agler, M., Kemen, E., Schulze-Lefert, P., et al. (2018). Microbial interkingdom interactions in roots promote Arabidopsis survival. Cell 175, 973–983. doi: 10.1016/j.cell.2018.10.020
Fernando, A. A., and Currah, R. S. (1996). A comparative study of the effects of the root endophytes Leptodontidium orchidicola and Phialocephala fortinii (Fungi Imperfecti) on the growth of some subalpine plants in culture. Can. J. Bot. 74, 1071–1078. doi: 10.1139/b96-131
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Frey-Klett, P., Garbaye, J. A., and Tarkka, M. (2007). The mycorrhiza helper bacteria revisited. New Phytol. 176, 22–36. doi: 10.1111/j.1469-8137.2007.02191.x
Gao, F.-K., Dai, C.-C., and Liu, X.-Z. (2010). Mechanisms of fungal endophytes in plant protection against pathogens. Afr. J. Microbiol. Res. 4, 1346–1351.
Grünig, C. R., Queloz, V., Sieber, T. N., and Holdenrieder, O. (2008). Dark septate endophytes (DSE) of the Phialocephala fortinii s.l. –Acephala applanata species complex in tree roots: classification, population biology, and ecology. Botany 86, 1355–1369. doi: 10.1139/B08-108
Hacquard, S., Garrido-Oter, R., González, A., Spaepen, S., Ackermann, G., Lebeis, S., et al. (2015). Microbiota and host nutrition across plant and animal kingdoms. Cell Host Microbe 17, 603–616. doi: 10.1016/j.chom.2015.04.009
Hacquard, S., Spaepen, S., Garrido-Oter, R., and Schulze-Lefert, P. (2017). Interplay between innate immunity and the plant microbiota. Ann. Rev. Phytopathol. 55, 565–589. doi: 10.1146/annurev-phyto-080516-035623
Hamady, M., Walker, J. J., Harris, J. K., Gold, N. J., and Knight, R. (2008). Error-correcting barcoded primers for pyrosequencing hundreds of samples in multiplex. Nat. Methods 5, 235–237. doi: 10.1038/nmeth.1184
Hamonts, K., Trivedi, P., Garg, A., Janitz, C., Grinyer, J., Holford, P., et al. (2018). Field study reveals core plant microbiota and relative importance of their drivers. Environ. Microbiol. 20, 124–140. doi: 10.1111/1462-2920.14031
Han, J.-H., Kim, T.-S., Joung, Y., Kim, M. N., Shin, K.-S., Bae, T., et al. (2013). Nocardioides endophyticus sp. nov. and Nocardioides conyzicola sp. nov., isolated from herbaceous plant roots. Int. J. Syst. Evol. Microbiol. 63, 4730–4734. doi: 10.1099/ijs.0.054619-0
Hardoim, P. R., Van Overbeek, L. S., Berg, G., Pirttilä, A. M., Compant, S., Campisano, A., et al. (2015). The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiol. Mol. Biol. Rev. 79, 293–320. doi: 10.1128/MMBR.00050-14
Hardoim, P. R., Van Overbeek, L. S., and Van Elsas, J. D. (2008). Properties of bacterial endophytes and their proposed role in plant growth. Trends Microbiol. 16, 463–471. doi: 10.1016/j.tim.2008.07.008
He, J. W., Chen, G. D., Gao, H., Yang, F., Li, X. X., Peng, T., et al. (2012). Heptaketides with antiviral activity from three endolichenic fungal strains Nigrospora sp., Alternaria sp. and Phialophora sp. Fitoterapia 83, 1087–1091. doi: 10.1016/j.fitote.2012.05.002
Hiruma, K., Gerlach, N., Sacristán, S., Nakano, R. T., Hacquard, S., Kracher, B., et al. (2016). Root endophyte Colletotrichum tofieldiae confers plant fitness benefits that are phosphate status dependent. Cell 165, 464–474. doi: 10.1016/j.cell.2016.02.028
Hoffman, M. T., and Arnold, A. E. (2010). Diverse bacteria inhabit living hyphae of phylogenetically diverse fungal endophytes. Appl. Environ. Microbiol. 76, 4063–4075. doi: 10.1128/AEM.02928-09
Hollingsworth, P. M., Forrest, L. L., Spouge, J. L., Hajibabaei, M., Ratnasingham, S., Van Der Bank, M., et al. (2009). A DNA barcode for land plants. Proc. Natl. Acad. Sci. U.S.A. 106, 12794–12797. doi: 10.1073/pnas.0905845106
Huson, D. H., Auch, A. F., Qi, J., and Schuster, S. C. (2007). MEGAN analysis of metagenomic data. Genome Res. 17, 377–386. doi: 10.1101/gr.5969107
Johnson, D. B., Bacelar-Nicolau, P., Okibe, N., Thomas, A., and Hallberg, K. B. (2009). Ferrimicrobium acidiphilum gen. nov., sp. nov. and Ferrithrix thermotolerans gen. nov., sp. nov.: heterotrophic, iron-oxidizing, extremely acidophilic actinobacteria. Int. J. Syst. Evol. Microbiol. 59, 1082–1089. doi: 10.1099/ijs.0.65409-0
Kadowaki, K., Yamamoto, S., Sato, H., Tanabe, A. S., Hidaka, A., and Toju, H. (2018). Mycorrhizal fungi mediate the direction and strength of plant–soil feedbacks differently between arbuscular mycorrhizal and ectomycorrhizal communities. Commun. Biol. 1:196. doi: 10.1038/s42003-018-0201-9
Kim, J. C., Choi, G. J., Park, J. H., Kim, H. T., and Cho, K. Y. (2001). Activity against plant pathogenic fungi of phomalactone isolated from Nigrospora sphaerica. Pest Manag. Sci. 57, 554–559. doi: 10.1002/ps.318
Kulichevskaya, I. S., Ivanova, A. O., Baulina, O. I., Bodelier, P. L., Damste, J. S. S., and Dedysh, S. N. (2008). Singulisphaera acidiphila gen. nov., sp. nov., a non-filamentous, Isosphaera-like planctomycete from acidic northern wetlands. Int. J. Syst. Evol. Microbiol. 58, 1186–1193. doi: 10.1099/ijs.0.65593-0
Kwak, M.-J., Kong, H. G., Choi, K., Kwon, S.-K., Song, J. Y., Lee, J., et al. (2018). Rhizosphere microbiome structure alters to enable wilt resistance in tomato. Nat. Biotechnol. 36, 1100–1109. doi: 10.1038/nbt.4232
Langille, M. G., Zaneveld, J., Caporaso, J. G., Mcdonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Lian, C., Narimatsu, M., Nara, K., and Hogetsu, T. (2006). Tricholoma matsutake in a natural Pinus densiflora forest: correspondence between above-and below-ground genets, association with multiple host trees and alteration of existing ectomycorrhizal communities. New Phytol. 171, 825–836. doi: 10.1111/j.1469-8137.2006.01801.x
Lundberg, D. S., Lebeis, S. L., Paredes, S. H., Yourstone, S., Gehring, J., Malfatti, S., et al. (2012). Defining the core Arabidopsis thaliana root microbiome. Nature 488, 86–90. doi: 10.1038/nature11237
Lundberg, D. S., Yourstone, S., Mieczkowski, P., Jones, C. D., and Dangl, J. L. (2013). Practical innovations for high-throughput amplicon sequencing. Nat. Methods 10, 999–1002. doi: 10.1038/nmeth.2634
Majewska, M. L., Błaszkowski, J., Nobis, M., Rola, K., Nobis, A., Łakomiec, D., et al. (2015). Root-inhabiting fungi in alien plant species in relation to invasion status and soil chemical properties. Symbiosis 65, 101–115. doi: 10.1007/s13199-015-0324-4
Mangan, S. A., Schnitzer, S. A., Herre, E. A., Mack, K. M., Valencia, M. C., Sanchez, E. I., et al. (2010). Negative plant-soil feedback predicts tree-species relative abundance in a tropical forest. Nature 466, 752–755. doi: 10.1038/nature09273
Mercier, J., and Lindow, S. (2000). Role of leaf surface sugars in colonization of plants by bacterial epiphytes. Appl. Environ. Microbiol. 66, 369–374. doi: 10.1128/AEM.66.1.369-374.2000
Mitchell, C. E., and Power, A. G. (2003). Release of invasive plants from fungal and viral pathogens. Nature 421, 625–627. doi: 10.1038/nature01317
Moora, M., Berger, S., Davison, J., Öpik, M., Bommarco, R., Bruelheide, H., et al. (2011). Alien plants associate with widespread generalist arbuscular mycorrhizal fungal taxa: evidence from a continental-scale study using massively parallel 454 sequencing. J. Biogeogr. 38, 1305–1317. doi: 10.1111/j.1365-2699.2011.02478.x
Newsham, K. K. (2011). A meta-analysis of plant responses to dark septate root endophytes. New Phytol. 190, 783–793. doi: 10.1111/j.1469-8137.2010.03611.x
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Obeso, J. R. (2002). The costs of reproduction in plants. New Phytol. 155, 321–348. doi: 10.1046/j.1469-8137.2002.00477.x
Ohashi, H., Kadota, Y., Murata, J., Yonekura, K., and Kihara, H. (2016). Wild Flowers of Japan. Tokyo: Heibonsha.
Oksanen, J., Blanachet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’hara, R. B., et al. (2012). Vegan: Community Ecology Package. R Package Version 2.0-3. Available at: http://CRAN.R-project.org/package=vegan
Öpik, M., Metsis, M., Daniell, T., Zobel, M., and Moora, M. (2009). Large-scale parallel 454 sequencing reveals host ecological group specificity of arbuscular mycorrhizal fungi in a boreonemoral forest. New Phytol. 184, 424–437. doi: 10.1111/j.1469-8137.2009.02920.x
Pankratov, T. A., Tindall, B. J., Liesack, W., and Dedysh, S. N. (2007). Mucilaginibacter paludis gen. nov., sp. nov. and Mucilaginibacter gracilis sp. nov., pectin-, xylan-and laminarin-degrading members of the family Sphingobacteriaceae from acidic Sphagnum peat bog. Int. J. Syst. Evol. Microbiol. 57, 2349–2354. doi: 10.1099/ijs.0.65100-0
Parniske, M. (2008). Arbuscular mycorrhiza: the mother of plant root endosymbioses. Nat. Rev. Microbiol. 6, 763–775. doi: 10.1038/nrmicro1987
Peay, K. G., Kennedy, P. G., and Talbot, J. M. (2016). Dimensions of biodiversity in the Earth mycobiome. Nat. Rev. Microbiol. 14, 434–447. doi: 10.1038/nrmicro.2016.59
Peay, K. G., Russo, S. E., Mcguire, K. L., Lim, Z., Chan, J. P., Tan, S., et al. (2015). Lack of host specificity leads to independent assortment of dipterocarps and ectomycorrhizal fungi across a soil fertility gradient. Ecol. Lett. 18, 807–816. doi: 10.1111/ele.12459
Pieterse, C. M., Zamioudis, C., Berendsen, R. L., Weller, D. M., Van Wees, S. C., and Bakker, P. A. (2014). Induced systemic resistance by beneficial microbes. Ann. Rev. Phytopathol. 52, 347–375. doi: 10.1146/annurev-phyto-082712-102340
Qin, S., Xing, K., Jiang, J.-H., Xu, L.-H., and Li, W.-J. (2011). Biodiversity, bioactive natural products and biotechnological potential of plant-associated endophytic actinobacteria. Appl. Microbiol. Biotechnol. 89, 457–473. doi: 10.1007/s00253-010-2923-6
Qin, S., Yuan, B., Zhang, Y. J., Bian, G. K., Tamura, T., Sun, B. Z., et al. (2012). Nocardioides panzhihuaensis sp. nov., a novel endophytic actinomycete isolated from medicinal plant Jatropha curcas L. Antonie Van Leeuwenhoek 102, 353–360. doi: 10.1007/s10482-012-9745-8
Ramamoorthy, V., Viswanathan, R., Raguchander, T., Prakasam, V., and Samiyappan, R. (2001). Induction of systemic resistance by plant growth promoting rhizobacteria in crop plants against pests and diseases. Crop Prot. 20, 1–11. doi: 10.1016/S0261-2194(00)00056-9
R-Core-Team (2018). R 3.5.1: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Reinhart, K. O., and Callaway, R. M. (2006). Soil biota and invasive plants. New Phytol. 170, 445–457. doi: 10.1111/j.1469-8137.2006.01715.x
Rognes, T., Mahé, F., Flouri, T., Quince, C., and Nichols, B. (2014). Vsearch: Program. Available at: https://github.com/torognes/vsearch
Santoyo, G., Moreno-Hagelsieb, G., Orozco-Mosqueda Mdel, C., and Glick, B. R. (2016). Plant growth-promoting bacterial endophytes. Microbiol. Res. 183, 92–99. doi: 10.1016/j.micres.2015.11.008
Schardl, C. L., and Phillips, T. D. (1997). Protective grass endophytes: where are they from and where are they going? Plant Dis. 81, 430–438. doi: 10.1094/PDIS.1997.81.5.430
Schlaeppi, K., and Bulgarelli, D. (2015). The plant microbiome at work. Mol. Plant Microbe Interact. 28, 212–217. doi: 10.1094/MPMI-10-14-0334-FI
Stevens, J. L., Jackson, R. L., and Olson, J. B. (2013). Slowing PCR ramp speed reduces chimera formation from environmental samples. J. Microbiol. Methods 93, 203–205. doi: 10.1016/j.mimet.2013.03.013
Tanabe, A. S. (2018). Claident v0.2.2018.05.29, A Software Distributed by Author. Available at: http://www.fifthdimension.jp/
Tanabe, A. S., and Toju, H. (2013). Two new computational methods for universal DNA barcoding: a benchmark using barcode sequences of bacteria, archaea, animals, fungi, and land plants. PLoS One 8:e76910. doi: 10.1371/journal.pone.0076910
Tedersoo, L., May, T. W., and Smith, M. E. (2010). Ectomycorrhizal lifestyle in fungi: global diversity, distribution, and evolution of phylogenetic lineages. Mycorrhiza 20, 217–263. doi: 10.1007/s00572-009-0274-x
Toju, H., Guimarães, P. R. Jr., Olesen, J. M., and Thompson, J. N. (2014). Assembly of complex plant–fungus networks. Nat. Commun. 5:5273. doi: 10.1038/ncomms6273
Toju, H., Peay, K. G., Yamamichi, M., Narisawa, K., Hiruma, K., Naito, K., et al. (2018a). Core microbiomes for sustainable agroecosystems. Nat. Plants 4, 247–257. doi: 10.1038/s41477-018-0139-4
Toju, H., Tanabe, A. S., and Sato, H. (2018b). Network hubs in root-associated fungal metacommunities. Microbiome 6:116. doi: 10.1186/s40168-018-0497-1
Toju, H., and Sato, H. (2018). Root-associated fungi shared between arbuscular mycorrhizal and ectomycorrhizal conifers in a temperate forest. Front. Microbiol. 9:433. doi: 10.3389/fmicb.2018.00433
Toju, H., Tanabe, A., and Ishii, H. (2016a). Ericaceous plant–fungus network in a harsh alpine–subalpine environment. Mol. Ecol. 25, 3242–3257. doi: 10.1111/mec.13680
Toju, H., Yamamoto, S., Tanabe, A. S., Hayakawa, T., and Ishii, H. S. (2016b). Network modules and hubs in plant-root fungal biome. J. R. Soc. Interface 13:20151097. doi: 10.1098/rsif.2015.1097
Toju, H., Tanabe, A. S., Yamamoto, S., and Sato, H. (2012). High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. PLoS One 7:e40863. doi: 10.1371/journal.pone.0040863
Toju, H., Yamamoto, S., Sato, H., Tanabe, A. S., Gilbert, G. S., and Kadowaki, K. (2013). Community composition of root-associated fungi in a Quercus-dominated temperate forest: “codominance” of mycorrhizal and root-endophytic fungi. Ecol. Evol. 3, 1281–1293. doi: 10.1002/ece3.546
van der Heijden, M. G., Bardgett, R. D., and Van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
van der Heijden, M. G., Klironomos, J. N., Ursic, M., Moutoglis, P., Streitwolf-Engel, R., Boller, T., et al. (1998). Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. Nature 396, 69–72. doi: 10.1038/23932
van der Putten, W. H., Bardgett, R. D., Bever, J. D., Bezemer, T. M., Casper, B. B., Fukami, T., et al. (2013). Plant–soil feedbacks: the past, the present and future challenges. J. Ecol. 101, 265–276. doi: 10.1111/1365-2745.12054
Vereecke, D., Burssens, S., Simón-Mateo, C., Inzé, D., Van Montagu, M., Goethals, K., et al. (2000). The Rhodococcus fascians-plant interaction: morphological traits and biotechnological applications. Planta 210, 241–251. doi: 10.1007/PL00008131
Wagner, M. R., Lundberg, D. S., Tijana, G., Tringe, S. G., Dangl, J. L., and Mitchell-Olds, T. (2016). Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. Nat. Commun. 7:12151. doi: 10.1038/ncomms12151
Weston, D. J., Pelletier, D. A., Morrell-Falvey, J. L., Tschaplinski, T. J., Jawdy, S. S., Lu, T.-Y., et al. (2012). Pseudomonas fluorescens induces strain-dependent and strain-independent host plant responses in defense networks, primary metabolism, photosynthesis, and fitness. Mol. Plant Microbe Interact. 25, 765–778. doi: 10.1094/MPMI-09-11-0253
Keywords: biodiversity, endophytes, host specificity, mycorrhizal fungi, plant-associated microbiomes, plant-soil feedbacks
Citation: Toju H, Kurokawa H and Kenta T (2019) Factors Influencing Leaf- and Root-Associated Communities of Bacteria and Fungi Across 33 Plant Orders in a Grassland. Front. Microbiol. 10:241. doi: 10.3389/fmicb.2019.00241
Received: 21 October 2018; Accepted: 29 January 2019;
Published: 19 February 2019.
Edited by:
Devin Coleman-Derr, United States Department of Agriculture, United StatesReviewed by:
Derek Persoh, Ruhr-Universität Bochum, GermanyCollin M. Timm, Johns Hopkins University, United States
Copyright © 2019 Toju, Kurokawa and Kenta. 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) and the copyright owner(s) 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: Hirokazu Toju, dG9qdS5oaXJva2F6dS40Y0BreW90by11LmFjLmpw