Structural variability and niche differentiation in the rhizosphere and endosphere bacterial microbiome of field-grown poplar trees
© The Author(s). 2017
Received: 11 October 2016
Accepted: 3 February 2017
Published: 23 February 2017
The plant microbiome represents one of the key determinants of plant health and productivity by providing a plethora of functional capacities such as access to low-abundance nutrients, suppression of phytopathogens, and resistance to biotic and/or abiotic stressors. However, a robust understanding of the structural composition of the bacterial microbiome present in different plant microenvironments and especially the relationship between below-ground and above-ground communities has remained elusive. In this work, we addressed hypotheses regarding microbiome niche differentiation and structural stability of the bacterial communities within different ecological plant niches.
We sampled the rhizosphere soil, root, stem, and leaf endosphere of field-grown poplar trees (Populus tremula × Populus alba) and applied 16S rRNA amplicon pyrosequencing to unravel the bacterial communities associated with the different plant habitats.
We found that the structural variability of rhizosphere microbiomes in field-grown poplar trees (P. tremula × P. alba) is much lower than that of the endosphere microbiomes. Furthermore, our data not only confirm microbiome niche differentiation reports at the rhizosphere soil–root interface but also clearly show additional fine-tuning and adaptation of the endosphere microbiome in the stem and leaf compartment. Each plant compartment represents an unique ecological niche for the bacterial communities. Finally, we identified the core bacterial microbiome associated with the different ecological niches of Populus.
Understanding the complex host–microbe interactions of Populus could provide the basis for the exploitation of the eukaryote–prokaryote associations in phytoremediation applications, sustainable crop production (bio-energy efficiency), and/or the production of secondary metabolites.
KeywordsPopulus tremula × Populus alba Bacterial microbiome Rhizosphere Endosphere Microbiome niche differentiation 16S rRNA amplicon pyrosequencing
Inter-organismal associations between eukaryotic and prokaryotic organisms are one of the most studied research areas in (micro)biology in recent years. The massive interest in this topic is reflected by numerous studies ranging from the human microbiome [1, 2] and host–genotype associations therein  and gut microfauna of insects [4–6] to microbiota associated with plants [7–17]. In fact, most eukaryotes maintain close mutualistic relationships with microorganisms that are, in most cases, linked to their nutrient acquisition and thereby crucial for their performance and survival [18, 19]. Furthermore, the associated prokaryotic communities may play important roles in the regulation of the eukaryote immune system [20–23].
Plant–microbe interactions are of specific interest, not only to get a better understanding of their role during plant growth and development but also to allow exploitation of their relationships in phytoremediation applications, sustainable crop production, and the production of secondary metabolites [24–26]. The plant microbiome, often referred to as the host’s second or extended genome, comprises diverse microbial classes, including bacteria and archaea, fungi, oomycetes, and viruses. In its entirety, the plant microbiome represents one of the key determinants of plant health and productivity by providing a plethora of functional capacities [27–30]. More specifically, bacterial microbiota may improve nutrient bioavailability and transport from the soil as well as increase host tolerance to biotic (and abiotic stresses), promote stress resistance, and influence crop yield and quality. In return, the host plant delivers habitation and a constant supply of energy and carbon sources to the microbiota [29, 31]. Virtually all tissues of a plant host bacterial communities: at the soil–root interface (rhizosphere/rhizoplane), inside the plants tissues (root, stem, and leaf endosphere), and at the air–plant interface (phyllosphere environment). To a lesser extent, we can also distinguish the bacterial colonization of the anthosphere (flower) , the spermosphere (seeds) [33, 34], and the carposphere (fruit) . All these microenvironments provide specific biotic and abiotic conditions for the residing bacterial communities.
Within plant–bacteria research, most attention has been dedicated to niche differentiation of bacterial communities at the rhizosphere soil–root interface [12, 14–16, 36–38]. For example, Gottel et al. compared the bacterial (and fungal) microbiota of mature poplar (Populus deltoides) trees using 16S ribosomal RNA (rRNA) gene pyrosequencing and revealed highly different root endophytic bacterial communities as compared to the rhizosphere soil . Bulgarelli et al.  and Lundberg et al.  obtained qualitatively similar results in a study on the bacterial root microbiota of Arabidopsis. In contrast to the knowledge concerning the differentiation of the bacterial microbiome at the rhizosphere–endosphere barrier, a robust understanding of the structural composition of the bacterial microbiome present in different plant microenvironments and especially the relationship between below-ground and above-ground communities in field conditions has remained elusive. Recently, Coleman-Derr et al., Fonseca-Garcia et al., and Tardif et al. observed significant plant compartment effects respectively in the microbiome of Agave species, cacti, and willow [39–41]. Alternatively, Ottesen et al. reported significant differentiation of the epiphytic microbiome across different plant organs of tomato plants . Other studies have focused on the leaf and root microbiomes [43, 44].
Here, we evaluate microbiome niche differentiation of bacterial communities associated with the rhizosphere soil and the root, stem, and leaf endosphere of field-grown wild-type poplar trees (Populus tremula × Populus alba) using 16S rRNA pyrosequencing. Populus is widely considered as the model of choice to study the biology of woody perennials and also provides an ideal model to explore the large variety of plant–microbe interactions [8, 9, 15, 36, 45–47]. Hybrid poplars are among the fastest growing trees and provide high economic flexibility with end-use applications such as the production of biofuels, pulp, and paper and other bio-based products such as chemicals and adhesives . Furthermore, poplar trees can be grown on marginal land (land not suitable for food production) thereby evading the food versus fuel debate [49–51]. Sequencing of the poplar genome along with the availability of large natural provenances and breeding pedigrees, and the first successful use of gene editing have also opened biotechnological possibilities to accelerate breeding and genetic engineering [52–58]. In the present study, we focussed on two main questions: (i) How variable are bacterial communities associated with different field-grown trees within the same study site? (ii) Do bacterial communities in the endosphere differentiate among the plant niches, and how do they relate to the rhizosphere communities?
Quality metrics of pyrosequencing analysis
Quality metrics of pyrosequencing analysis
A. Total number of reads and read length before and after quality checking and trimming
Total # of raw reads before QC
Average read length before QC
405 ± 96
Total # of assigned reads after QC
Average read length after QC
207 ± 4
B. Assigned reads
Average # of reads
5058 ± 615
5311 ± 643
2761 ± 1174
3034 ± 960
26.09 ± 0.01a
5.01 ± 0.55b
2.60 ± 0.35b
2.21 ± 0.65b
Normalization to 2000 reads per sample
C. Non-target rRNA (%)
0.01 ± 0.01
0.44 ± 0.17
0.03 ± 0.02
0.03 ± 0.01
D. Unclassified reads
34.07 ± 1.10a
4.74 ± 0.32b
19.05 ± 4.32b
3.59 ± 1.03b
Furthermore, we determined the co-amplification of non-target 16S rRNA (archaeal, chloroplast, and mitochondrial sequences) and the number of singletons identified within each plant compartment (%), as well as the number of reads that could not be unambiguously classified at the phylum level (Table 1). We found a distinct plant compartment effect in the retrieval of singletons. Significantly more singletons were obtained from the rhizosphere soil as compared to all other plant compartments (F (3, 44) = 454.7, P < 0.001) (Table 1). Under our optimized PCR conditions , no mitochondrial 16S rRNA sequences were co-amplified from any of the plant compartments. Minute fractions of chloroplast/plastidal 16S rRNA sequences were co-amplified from root, stem, and leaf samples (ranging from 0.01 to 0.44% of the normalized reads). Finally, in the rhizosphere, we identified a small portion of reads, which were assigned to the taxonomic domain Archaea (0.03%). In the rhizosphere soil, a large fraction of reads could not be unambiguously classified at the phylum level (34.07%). In the plant compartments, we were able to classify the majority of reads and only a relatively small proportion of reads remained unclassified (ranging from 3.59 to 19.05%). Unclassified reads at the phylum level were removed from the dataset for further analysis (Table 1).
Alpha rarefaction curves and alpha diversity
Alpha diversity, the microbial diversity within each sample, was analyzed based on the OTU richness, the inverse Simpson diversity index, and Pielou’s evenness (Fig. 2). To control for differences in sampling effort across plant compartments, we rarefied each sample to 2000 sequences per sample before calculating the diversity indices. OTU richness was highly dependent on plant compartment (P < 0.05) with high richness values for rhizosphere soil (848.9 ± 7.9) and consistently decreased richness estimates in the root samples (243.7 ± 9.6) and stem samples (126.7 ± 11.9). OTU richness indices of the leaf samples (118.3 ± 17.2) were comparable with those of the stem samples. For diversity and evenness estimates, we found a clear separation between the rhizosphere soil samples and endosphere samples (P < 0.05). Higher diversity and evenness measures were observed for the rhizosphere soil samples as compared to the samples of the endosphere plant compartments. In contrast, all endosphere compartments revealed highly comparable diversity and evenness estimates. Furthermore, to control for bias in the used community estimators, alternative estimators were calculated which resulted in highly similar conclusions (see Additional file 1).
Analysis of similarity (ANOSIM)
Rhizosphere soil vs root
Rhizosphere soil vs stem
Rhizosphere soil vs leaf
Root vs stem
Root vs leaf
Stem vs leaf
Top members of the bacterial microbiome within each plant compartment
Indicator species analysis
OTU (Genus or higher)
Relative abundance (%)
Finally, to provide a complete overview of the OTU distribution within the plant compartments, we calculated the proportion of OTUs uniquely identified in each specific plant compartment as well as the OTUs shared by the different plant compartments (Additional file 7). The proportion of OTUs shared by all plant compartments was 16.4%. Approximately 26% of all OTUs were exclusively found in the rhizosphere samples compared to the root samples (7.1%), the stem samples (3.2%), and the leaf samples (5.5%). Around 6% of the total OTUs were exclusively observed in the endosphere compartments. Finally, we clearly observed a higher overlap in OTUs from the rhizosphere soil samples and the root samples (11.7%) as compared to rhizosphere soil and the stem samples (1.8%) and the rhizosphere soil and leaf samples (2.8%).
Quality of the pyrosequencing analysis
We used an optimized PCR approach to reduce co-amplification of chloroplast and mitochondrial 16S rRNA . In many studies, the high homology between bacterial 16S rRNA genes, chloroplast 16S rRNA genes, and plant nuclear and mitochondrial 16S rRNA genes [60, 61] and moreover the high abundance of chloroplast 16S rRNA genes in these environments led to undesired co-amplification of non-target sequences [12, 15, 16, 36, 44]. Our optimized PCR approach resulted in very low co-amplification of these sequences and high retrieval of bacterial 16S rRNA sequences (Table 1). The highest retrieval of chloroplast 16S rRNA sequences was observed in the stem samples, corroborating results from our primer optimization  and reinforcing our view that the balance between the amount of endophytic bacterial DNA (bacterial cell density) and chloroplast DNA seems to play a more important role than the absolute chloroplast concentration. Finally, we also considered the number of singletons (sequences only found once in the dataset) obtained from each plant compartment (Table 1). Remarkably, we found high levels of singletons in rhizosphere soil and a decreasing number of singletons in other plant compartments. Singletons have been shown to comprise up to and beyond 60% of taxa in some surveys [62, 63] and are generally considered as being problematic since they represent inherently unreplicated data . Most singletons arise from DNA sequencing errors (insertions, deletions, low-quality reads, inadequate clustering and formation of chimeric sequences) [65–68] creating false sequences and artificially inflating diversity estimates [69–71]. In our experimental setup, sequencing error (and potential creation of erroneous (singleton) sequences) is expected to be similar for all plant compartments with possibly a minor impact of carry-over contaminants in the rhizosphere soil samples, which could potentially increase PCR error. A confounding factor in this respect could be the use of different DNA extractions kits for the rhizosphere samples and the endosphere samples. Previously, we focused on extracting DNA from all the studied plant compartments (rhizosphere soil, root, stem, and leaf samples). , but we were unable to extract high-quality DNA (and quantity) from all four plant compartments using the same DNA extraction kit. To ensure high-quality and quantity DNA from all studied plant compartments and reduce bias from low DNA retrieval, we opted for a different DNA extraction kit for the rhizosphere soil samples and the endosphere samples. Nonetheless, a certain amount of bias may have been introduced in the results as a consequence of differences in the lysis efficiency of different DNA extraction kits [72, 73]. Notwithstanding these elements, our results indicate that the high discrepancy in the number of singletons between the plant compartments could in fact be attributable to more genuine rare (singleton) OTUs in the rhizosphere soil (Table 1). Indeed, the rhizosphere soil is renowned for its vast microbial diversity [74, 75]. For further analysis, we chose a conservative approach and treated all singletons as potentially erroneous and removed them from the data sets [68, 76]. However, the involvement of this rare biosphere in community dynamics and their ecological roles are largely unknown, but they could contribute to community stability by enabling fast responses to altering environmental conditions .
(i) How variable are bacterial communities associated with different field-grown trees?
We observed remarkably dissimilar shapes of the OTU rarefaction curves when comparing rhizosphere soil and endosphere samples (Fig. 1). Rhizosphere soil samples displayed uniform rarefaction curves (Fig. 1a) whereas the variation in the shape of the rarefaction curves from the endophytic samples was much higher, especially for the stem and leaf samples (Fig. 1b–d). High variability of endophytic OTU richness, as depicted by the rarefaction curves, could possibly be caused by sporadic and non-uniform colonization of the roots and aerial plant compartments of Populus . Gottel et al. attributed part of the variation to their inability to sequence the bacterial endophytic community deeply and uniformly enough because of the high co-amplification of organellar 16S rRNA (67,000 chloroplast and 65,000 mitochondrial sequences) . However, our data exhibit roughly the same pattern without the co-amplification of non-target DNA (Table 1) and with high Good’s coverage estimates (Fig. 1). Therefore, our data suggest considerable variation in endophytic colonization as a major reason for the high variability in the rarefaction curves. Indeed, rhizosphere/rhizoplane colonization is primarily driven by (a) the deposition of large amounts of carbon (e.g., root exudates, mucilage by the root caps, etc.) by plants (rhizodeposition) and (b) the relatively simple or inelaborate chemo-attraction of the bacteria (and other microorganisms) to the root exudates [78–81]. Although, since root exudates and mucilage-derived nutrients attract a myriad of organisms to the rhizosphere environment, plant-associated bacteria have to be highly competitive to successfully colonize the root zone . In contrast to rhizosphere/rhizoplane colonization, endophytic competence (i.e., ability to successfully colonize the host plant) can require specific traits (e.g., expression of genes involved in chemotaxis, the formation of flagella and pilli, the production of cell-wall degrading enzymes, etc.) and intricate interplay between rhizospheric soil-borne bacteria and the host plants innate immune system [12, 20, 25, 30, 82].
Furthermore, we also clearly observed more variation in the bacterial community structures in the endosphere as compared to the rhizosphere communities, especially in stem and leaf samples (Fig. 3 and Additional file 3). As mentioned previously, a possible confounding factor in the interpretation of these results is the use of different DNA extraction kits for the rhizosphere and endosphere samples. Nonetheless, the within group variation, as depicted by ANOSIM analysis (Additional file 3), of rhizosphere soil bacterial assemblages is very low. The soil biome is one of the richest microbial ecosystems on Earth with an estimated bacterial diversity of >2000 species within 0.5 g of soil [74, 75, 83]. Furthermore, the root exudation process is heterogeneous in space and time [84, 85]. Despite these factors, the formation of distinctive rhizosphere bacterial communities mediated by rhizodeposition (and chemo-attraction to photoassimilates) seems to be a very consistent and stable process across different poplar individuals. In contrast, variation within endophytic communities is much higher (Fig. 3 and Additional file 3). As mentioned previously, endophytic colonization and formation of stable communities appears to be a more variable process, as suggested by our results from the alpha rarefaction curves (Fig. 1), from the PCA analyses (Fig. 3a), the relative abundance of bacterial phyla (Fig. 4) and the ANOSIM results of the bacterial community structures (Additional file 3). Crucial factors underlining this variability are the nature of endophytic colonization and competence (e.g., bacterial motility, ability to produce cell-wall degrading enzymes) [25, 82], interplay with the host plants innate immune system . and acute fluctuations in abiotic conditions (temperature, humidity, access to nutrients, etc.) which differ from the buffered fluctuations in the rhizosphere [31, 86]. However, in contrast, OTU richness and OTU diversity (Fig. 2) were markably higher in the rhizosphere soil as compared to the endosphere samples.
(ii) Do bacterial communities present in the endosphere differentiate within the plant niches, and how do they relate to the rhizosphere communities?
To control for differences in sampling effort across plant compartments, we rarefied each sample to 2000 sequences per sample, although rarefying and using linear models of abundance have been scrutinized recently by McMurdie and Holmes . Initially, we estimated alpha diversity focussing on OTU richness, evenness and diversity. We found that richness estimates were highly dependent on plant compartment with rhizosphere soil, root and stem compartments clearly differentiated from each other by decreasing OTU richness (Fig. 2). These results are in concordance with the general views of endophytic colonization. Rhizodeposition and root exudation by the host plant in the root zone fuels chemo-attraction and colonization of the rhizosphere soil and rhizoplane, thereby leading to the formation of distinctive, highly rich, and diverse rhizosphere microbiomes [78–81]. After rhizoplane colonization, adaptation to an endophytic lifestyle is dependent on the ability of the soil-borne bacteria to pass (actively or passively) the endodermis and pericycle, reach the xylem vessels, and finally lead to systemic colonization of the plant [25, 82]. Systemic plant colonization by certain bacterial species is re-enforced by the proportion of OTUs shared by all the plant compartments (16.4%, Additional file 7). The rhizosphere soil–root interface acts as a selective barrier, and endophytic competence/colonization is limited to specific bacterial species. The great loss of diversity and evenness (Fig. 2a–c) from rhizosphere soil to endophytic compartments supports this view and indicates that only a limited number of bacteria can adapt to an endophytic lifestyle (loss of diversity) (Fig. 2c) and these bacterial strains will therefore dominate endophytic assemblages (loss of evenness) (Fig. 2b).
To compare the bacterial community structures present in the plant compartments, we clustered all samples using principal component analysis (PCA) and hierarchical clustering (Bray–Curtis dissimilarities) (Fig. 3). At the phylum level and OTU level, all samples strongly clustered according to plant compartment (P < 0.01) and rendered microbiota significantly dissimilar from each other (Fig. 3 and Table 2) (see Additional file 2). Again to put the results in a broader context, the caveat of using different DNA extraction kits for the rhizosphere samples and the endosphere samples may have introduced a certain amount of bias in these results. However, previously, we observed the same niche differentiation for the cultivable bacteria of poplar trees in the same field study . Niche differentiation between rhizosphere and root endophyte microbiome has also been described for mature poplar trees growing in natural ecosystems (P. deltoides) [15, 36], for Arabidopsis thaliana [12, 16] and other plant species [10, 25, 37]. Recently, Bulgarelli et al.  proposed a two-step selection model for root microbiota differentiation from the rhizosphere where rhizodeposition and host genotype-dependent fine-tuning converge to select specific endophytic assemblages. Bulgarelli et al. argue that substrate-driven selection in the rhizosphere is expected to persist in the endosphere . Indeed, our data suggest additional fine-tuning and niche differentiation of microbiota in the aerial plant organs (both at the phylum and OTU level), with the stem and leaf bacterial assemblages being remarkably dissimilar from the root and rhizosphere (Fig. 3 and Additional file 2) (Table 2). This in agreement with the studies of (a) Coleman-Derr et al.  and Fonseca-Garcia et al. , who revealed that the composition of bacterial communities in plants native to semi-arid and arid ecosystems (Agave species and cacti) were primarily determined by the plant compartment and (b) Tardif et al., who reported significant plant compartment effects in the willow microbiome . Each of the plant microenvironments or ecological niches (rhizosphere soil, root, stem, and leaf) provide relevant biotic and abiotic gradients such as availability of soluble organic compounds [31, 88, 89]. This is further highlighted by the existence of specific proportions of OTUs, which were exclusively found in different plant compartments (e.g., 25.7% unique OTUs in the rhizosphere soil samples (Additional file 7)). The distribution of all identified OTUs across the different plant compartments (Additional file 7) also highlights several other aspects: (a) the inability of a large number of OTUs to colonize the plant (25.7% of all OTUs), (b) the existence of obligate endophytes which are only observed in the endosphere compartments (5.9% of all OTUs) and are strictly dependent on their host plant for survival , (c) the existence of facultative endophytes which may exist inside (endosphere) and outside the host plant (rhizosphere soil) , and (d) although most endophytic bacteria colonizing the host plant originate from the rhizosphere soil , some may originate elsewhere (e.g., colonization of the phyllosphere via aerosols and subsequently the leaf endosphere ) as evidenced by the proportion of OTUs uniquely identified in the leaf samples (5.5%).
Drivers of microbiome niche differentiation
At the phylum level, Actinobacteria and Proteobacteria (mostly Alpha- and Betaproteobacteria) and to a lesser extent Bacteroidetes, Firmicutes, and Acidobacteria dominated the rhizobacterial assemblages. The ratio between Proteobacteria and Acidobacteria in rhizosphere bacterial communities has previously been shown to be an indicator of soil nutrient-content where Proteobacteria were linked to nutrient-rich soils and Acidobacteria to nutrient-poor soils [36, 91, 92]. Similarly to studies in Arabidopsis [12, 16], rice , and poplar [15, 36], the relative abundance of Acidobacteria and Actinobacteria decreased from the rhizosphere soil to the root microbiota whereas the relative abundance of Proteobacteria (mostly Alpha) increased in the root endosphere. Across different unrelated plant host species, the host-associated bacterial microbiota in the rhizosphere and root endosphere are consistently enriched with members belonging to the phylum Proteobacteria [12, 14–17, 36, 37, 93–95]. From our results, we can conclude that also the stem and leaf microbiota are dominated by Proteobacteria albeit with different OTU level members, mostly belonging to the Alpha- and Gammaproteobacteria (Figs. 4 and 5 and Additional files 4 and 5). The large overlap in key community members of endophytic bacterial assemblages across different plant host species demonstrates that endophytic competence (efficient colonization) and dealing with niche-specific plant settings (nutrient availability/variability, oxygen levels, etc.) is reserved for a minority of bacterial phyla. Enrichment and depletion of specific bacteria within the plant-associated microbiome are not passive processes but rather depend on active selection of microbial consortia by the plant host and/or opportunistic colonization of the available ecological niches by certain bacteria [14, 19, 31]. A remarkable phylum, candidate division TM7 (recently renamed phylum Candidatus Saccharibacteria), which has only been described from 16S rRNA gene sequence and genome data [96, 97], showed highly variable colonization capacities in the root and stem endosphere (Fig. 4 and Additional files 4 and 5). Phylum Candidatus Saccharibacteria is a highly ubiquitous phylum found in soils, sediments, wastewater, animals, and plant microbiomes [9, 12, 15, 97]. Furthermore, Shakya et al. also reported high variability in the relative abundance of phylum Candidatus Saccharibacteria (albeit in the rhizosphere microbiome of poplar)  possibly suggesting high sensitivity of these members to mild variations in abiotic and/or biotic stressors, strict nutritional requirements, variable responses of the plant’s innate immune system, strong influence of microbe-microbe interactions, or possible interactions with the plant host genotype.
Finally, at the OTU level (genus or higher), rhizosphere soil communities were dominated primarily by Rhizobiales, Actinomycetales, Burkholderiales, Arthrobacter, and Variovorax which were characteristically isolated from rhizosphere soil samples [10, 12, 36, 98]. Root endophytic assemblages were dominated by Rhizobiales, Rhizobium, and candidate division TM7 (with high variability). Dominant members of the stem samples are Pseudomonas, candidate division TM7, Methylobacterium, and Deinococcus. Finally, leaf samples mainly contained of Pseudomonas, Sphingomonas, and Methylobacterium. All of the above mentioned OTUs, which have been isolated from a variety of plant samples, may provide beneficial effects on plant health and growth [90, 99–102]. Remarkable in the stem (16%) and leaf endosphere (27%) is the efficient colonization capacity of Pseudomonas (Fig. 5 and Additional file 5). Niche-specific adaptation of Pseudomonas putida has previously been described by Wu et al. . We previously observed the same enrichment of Pseudomonas in the stem and leaf samples irrespective of the 16S rRNA primer pair used . Since aerosol samples were found to harbor abundant Pseudomonas and Sphingomonas sequences , enrichment of these bacteria in the leaf endosphere may occur via dual origins, colonization of the rhizosphere and/or leaf stomatal colonization. Furthermore, Sphingomonas and Methylobacterium, both abundantly present in the leaf endosphere, were shown to harbor specific adaptation strategies such as TonB-dependent receptors to survive in the phyllosphere environment [100, 104, 105].
We proved that the structural variability of rhizosphere microbiomes in field-grown poplar trees (P. tremula x P. alba) is much lower than that of the endosphere microbiomes. The formation of rhizosphere bacterial communities appears to be a more stable and controlled process whereas endophytic colonization of the roots, stems, and leaves is highly variable. Furthermore, our data not only confirm microbiome niche differentiation reports at the rhizosphere soil–root interface but also clearly show additional fine-tuning and adaptation of the endosphere microbiome in the stem and leaf compartment. Each plant compartment represents an unique ecological niche for the bacterial communities. Future studies which include the analysis of specific host genotype effects (clones, genetically modified genotypes, etc.) could provide more insight into the plasticity or responsiveness of the bacterial communities to specific changes in the host plant. Finally, we identified a core bacterial microbiome associated with the different ecological niches of Populus. This could provide the basis for more detailed (isolation) studies of the identified abundant OTUs and gain further insight into the complex host–microbe interactions of Populus.
Field trial and sampling
A poplar field trial located in Ghent, Belgium, was selected to obtain samples for this study. This field trial was established in April 2009 and contains female poplar clones (P. tremula × P. alba cv. “717-1B4”). Poplars were micropropagated in vitro, and ramets were grown in soil in the greenhouse for 9 months. Thereafter, the stems were cut 10 cm above soil level, and plants were maintained for an additional 10 days in the greenhouse. Finally, coppiced trees were transferred to the field in May 2009. The trees were planted in a density of 15,000 trees per hectare with an inter-plant distance of 0.75 m [54, 106]. Poplar trees were sampled in October 2012. At the time of sampling, the height of the trees was on average approximately 3.5–4.5 m. Fifteen individual trees were sampled for the rhizosphere soil and root samples, and 11 trees were sampled for the stem and leaf samples. Collected samples included rhizosphere soil, roots, stems, and leaves. The root samples were collected at a depth of 5–10 cm below ground level. The rhizosphere soil was strictly defined as soil particles adhering to the roots. For the stem and leaf samples, one complete branch of each of the 11 poplar individuals was collected. Sampled branches were directly connected to the central trunk and had on average a circumference of approximately 4–7 cm and a height of approximately 80–140 cm To standardize and maximize reproducibility of the stem samples, several small stem “cores” with bark (5–7 cores; 1 cm each) were collected from each branch from the base to the top of the offshoot to represent the stem compartment. For the leaf samples, all leaves from the sampled offshoot were collected to represent the leaf compartment.
Processing of samples
The samples were processed as described by Beckers et al. . Briefly, the root samples were depleted from soil particles by shaking on a platform (20 min, 120 rpm). The soil particles directly dislodged from roots represented the “rhizosphere soil” compartment. Subsequently “root,” “stem,” and “leaf” compartments were cleared from epiphytic bacteria by sequential washing (surface sterilization) with (a) sterile Millipore water (30 s), (b) 70% (v/v) ethanol (2 min), (c) sodium hypochlorite solution (2.5% active Cl− with 0.1% Tween 80) (5 min), and (d) 70% (v/v) ethanol (30 s) and finalized by rinsing the samples five times with sterile Millipore water. The plant samples were portioned into small fragments using a sterile scalpel and were subsequently macerated in a sterile phosphate saline buffer (PBS; 130 mM NaCl, 7 mM Na2HPO4, 3 mM NaH2PO4, pH 7.4) using a Polytron PR1200 mixer (Kinematica A6). Sterilization and homogenization of the plant samples were performed under aseptic conditions in a laminar airflow. Finally, quadruplicate aliquots of each sample (1.5 ml) of the homogenized plant material (root, stem, or leaf) were stored for all poplar individuals at −80 °C until DNA was extracted.
To minimize DNA extraction bias, DNA was extracted in quadruplicate from the rhizosphere soil, root, stem, and leaf samples [107, 108]. Approximately 250 mg of rhizosphere soil was used for each individual DNA extraction. DNA was extracted using the Power Soil DNA Isolation Kit following the protocol provided by the manufacturer (MoBio, Carlsbad, CA, USA). For the plant tissues, aliquots of homogenized plant material (1.5 ml) were first centrifuged (13,400 rpm, 30 min.) to collect all cells. Supernatants were discarded and DNA extractions were performed on pelleted plant material. DNA was extracted from plant samples using the Invisorb Spin Plant Mini Kit according to the manufacturer’s protocol (Stratec Biomedical AG, Birkenfeld, Germany).
PCR amplification and 454 pyrosequencing
Quadruplicate DNA samples from all compartments were individually amplified using a Techne TC-5000 thermocycler (Bibby Scientific Limited, Staffordshire, UK). Based on previous optimization experiments with 16S rRNA primer pairs , we selected primer 799F (5′-AACMGGATTAGATACCCKG-3′), with three mismatches with the poplar chloroplast 16S rRNA, and primer 1391R (5′-GACGGGCGGTGWGTRCA-3′). Furthermore, we included negative controls to evaluate the presence of contaminating sequences in reagents, which were checked using gel-electrophoresis (1.5% agarose gel, 90 V, 30 min.). A first round of PCR amplification was conducted using these primers without the Roche 454 pyrosequencing adaptors and sample-specific barcodes. Each 25 μl PCR reaction contained approximately 10 ng of DNA and was carried out using the FastStart High Fidelity PCR System (Roche Applied Science, Mannheim, Germany). Each reaction contained 2.75 μl FastStart 10× reaction buffer, 1.8 mM MgCl2, 0.2 mM dNTP mix, 0.4 μM of each primer, and 2 U FastStart HiFi polymerase. Cycling conditions included initial denaturation at 94 °C for 3 min, followed by 35 cycles of denaturation at 94 °C for 1 min, annealing at 53 °C for 1 min, and extension at 72 °C during 1 min; a final extension phase was performed at 72 °C during 10 min. PCR amplicon pools were cleared from residual primers and primer dimers by separating the PCR products on a 1.5% agarose gel (90 V, 30 min.), excising the bacterial product (amplicon length = 592 bp) and extracting the DNA from the gel slices using the QIAQuick gel extraction kit (Qiagen Benelux N.V., Venlo, The Netherlands). Mitochondrial by-products (1000 bp) were eliminated via this gel-purification. Following the first round of PCR amplification and gel-purification of the PCR products, a second round of PCR amplification was performed with primer 967F (5′CAACGCGAAGAACCTTACC-3′)-1391R(5′-GACGGGCGGTGWGTRCA-3′) to reduce the amplicon length (424 bp) to a more suitable length for 454 pyrosequencing. The forward primer (967F) was fused to the Roche 454 pyrosequencing adaptor A including a sample-specific 10-bp barcode (multiplex identifiers, MIDs). The reverse primer (1391R) was fused to adaptor B (Roche Applied Science, Mannheim, Germany). PCR cycling conditions were identical as described above, except for the number of PCR cycles that was lowered to 25.
Subsequently, quadruplicate PCR amplicon pools from the corresponding samples were grouped together resulting in 15 samples (rhizosphere soil and root) and 11 samples (stem and leaf) per plant compartment (15 biological replicates × 2 plant compartments + 11 biological replicates × 2 plant compartments = total of 52 samples). PCR amplicon pools were purified to remove PCR primers and primer dimers using the QIAquick PCR purification kit (Qiagen Benelux B.V., Venlo, the Netherlands). Following purification, the quality of the amplicon pools was evaluated using an Agilent 2100 Bioanalyzer system (Agilent Technologies, Diegem, Belgium) according to the manufacturer’s protocol. Finally, purified amplicon libraries were quantified with the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA) and a Fluostar Omega plate reader (BMG Labtech, Ortenberg, Germany) and pooled in equimolar concentrations. Rhizosphere samples (15) and root samples (15) were each separately pooled into two amplicon libraries. Stem and leaf samples were grouped into an additional library consisting of 22 samples (11 stem and 11 leaf samples). Each amplicon library (total of 3) was sequenced on one eighth of a Picotiter Plate on a Roche Genome Sequencer FLX+ using Titanium chemistry (Roche Applied Science, Mannheim, Germany) by Macrogen (Seoul, Korea).
Sequencing generated three individual standard flowgram format (SFF) files, which were analyzed separately using the software package mothur (version 1.33.2) following the standard operating protocol outlined in https://www.mothur.org/wiki/454_SOP . Briefly, the sequencing error was reduced by denoising (shhh.flows, Mothur implementation of Amplicon Noise algorithm) and quality trimming, which removed reads shorter than 200 bases, reads with homopolymers longer than 8 bases, and reads containing ambiguous bases. Unique sequences were identified, while archiving the abundance data of the unique sequences, and aligned using align.seqs with the SILVA reference alignment (Release 119) . Within the unique reads, chimeric sequences were identified using the Uchime tool  followed by their removal from the dataset. Unique sequences were classified using the mothur implementation of the Bayesian classifier where an 80% bootstrap cut-off value was used for assigning taxonomic classifications. Abundance data of sequences matching “Chloroplast” and “Mitochondria” were archived and these sequences were removed from the data sets. Subsequently, pairwise distances were calculated between all remaining unique sequences and a distance matrix was created. Average neighbor clustering was performed and, using a 0.03 OTU definition (97% sequence similarity cut-off level), a majority consensus taxonomy was obtained for each OTU. To minimize the impact of sequencing artefacts, we removed singletons from the datasets . Subsequently, rarefaction curves were assembled, and Good’s coverage scores were calculated in mothur based on 10,000 iterations. To calculate diversity indices (richness, diversity, evenness) while controlling for the sampling effort, each sample was rarefied to 2000 sequences. OTU richness, corresponding to the number of observed OTUs per sample (sobs), inverse Simpson diversity indices , and Pielou’s evenness indices  were calculated in mothur based on 10,000 iterations. To exclude bias in the community richness, evenness, and diversity estimators, we included several alternative estimators (Additional file 1).
Statistical analyses were performed in R 2.15.1 (The R Foundation for Statistical Computing, Vienna, Austria) . Normal distributions of the data were checked with the Shapiro–Wilk test and homoscedasticity of variances was analyzed using either Bartlett’s or the Fligner–Killeen test. Significant differences in the variance of parameters were evaluated, depending on the distribution of the estimated parameters, either with ANOVA or the Kruskal–Wallis rank sum test. Post hoc comparisons were conducted by either the Tukey’s honest significant differences tests or pairwise Wilcoxon rank sum tests. ANOVA was used to test the effect of the plant compartment (rhizosphere soil, root, stem, leaf) on the read abundances. Hierarchical clustering (based on Bray–Curtis dissimilarities) and principal component analyses (PCA) were performed in and displayed with PRIMER 7 . To statistically support the visual clustering of the bacterial communities in the PCA analyses, the different plant compartments were compared using permutation-based hypothesis tests: tests of the multivariate null hypotheses of no differences among a priori defined groups were examined using ANOSIM (an analog of univariate ANOVA) with the Spearman rank correlation method in PRIMER 7. Indicator species analysis was performed using the multipat function of the indicspecies package in R (version 1.7.1) . P values were corrected for multiple comparisons using the false discovery rate (FDR) with the Benjamini–Hochberg method. Taxonomic dendrogram (Fig. 5) was generated with one representative sequence of each OTU using Unipro UGENE and displayed with the use of iTOL (Interactive Tree Of Life) .
This work was funded by the Fund for Scientific Research Flanders (FWO-Vlaanderen) through the research project number G032912N, a Ph.D. grant for Michiel Op De Beeck, and a post-doc grant for Nele Weyens. Furthermore. This work has been financially supported by the UHasselt Methusalem project 08M03VGRJ.
Availability of data and materials
The standard flowgram format (SFF) files were deposited in the NCBI Sequence Read Archive (SRA) under the Bioproject number PRJNA341696.
BB, MODB, NW, JV, and WB planned and designed the research and experiments. BB performed the experiments and analyzed the data. BB wrote the manuscript. MODB, NW, JV, and WB reviewed and edited the manuscript. JV and WB acquired the funds for the study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464:59–65.PubMedPubMed CentralView ArticleGoogle Scholar
- Human Microbiome Project Consortium (HMPC). Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.View ArticleGoogle Scholar
- Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome. Nat Rev Microbiol. 2011;9:279–90.PubMedView ArticleGoogle Scholar
- Dillon RJ, Dillon VM. The gut bacteria of insects: nonpathogenic interactions. Annu Rev Entomol. 2004;49:71–92.PubMedView ArticleGoogle Scholar
- Hansen AK, Moran NA. The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol. 2014;23:1473–96.PubMedView ArticleGoogle Scholar
- Sudakaran S, Salem H, Kost C, Kaltenpoth M. Geographical and ecological stability of the symbiotic mid-gut microbiota in European firebugs, Pyrrhocoris apterus (Hemiptera, Pyrrhocoridae). Mol Ecol. 2012;21:6134–51.PubMedView ArticleGoogle Scholar
- Bai Y, Müller DB, Srinivas G, Garrido-Oter R, Potthoff E, Rott M, et al. Functional overlap of the Arabidopsis leaf and root microbiota. Nature. 2015;528:364–9.
- Beckers B, Op De Beeck M, Weyens N, Van Acker R, Van Montagu M, Boerjan W, et al. Lignin engineering in field-grown poplar trees affects the endosphere bacterial microbiome. Proc Natl Acad Sci. 2016;113:2312–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Beckers B, Op De Beeck M, Thijs S, Truyens S, Weyens N, Boerjan W, et al. Performance of 16s rDNA primer pairs in the study of rhizosphere and endosphere bacterial microbiomes in metabarcoding studies. Front Microbiol. 2016;7:650.PubMedPubMed CentralView ArticleGoogle Scholar
- Berg G, Smalla K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol Ecol. 2009;68:1–13.PubMedView ArticleGoogle Scholar
- Bonito G, Reynolds H, Robeson MS, Nelson J, Hodkinson BP, Tuskan G, et al. Plant host and soil origin influence fungal and bacterial assemblages in the roots of woody plants. Mol Ecol. 2014;23:3356–70.PubMedView ArticleGoogle Scholar
- Bulgarelli D, Rott M, Schlaeppi K, Ver Loren van Themaat E, Ahmadinejad N, Assenza F, et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature. 2012;488:91–5.PubMedView ArticleGoogle Scholar
- Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014;8:790–803.PubMedView ArticleGoogle Scholar
- Edwards J, Johnson C, Santos-Medellín C, Lurie E, Podishetty NK, Bhatnagar S, et al. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc Natl Acad Sci. 2015;112:E911–20.PubMedPubMed CentralView ArticleGoogle Scholar
- Shakya M, Gottel N, Castro H, Yang ZK, Gunter L, Labbé J, et al. A multifactor analysis of fungal and bacterial community structure in the root microbiome of mature Populus deltoides trees. PLoS One. 2013;8:e76382.PubMedPubMed CentralView ArticleGoogle Scholar
- Lundberg DS, Lebeis SL, Paredes SH, Yourstone S, Gehring J, Malfatti S, et al. Defining the core Arabidopsis thaliana root microbiome. Nature. 2012;488:86–90.PubMedPubMed CentralView ArticleGoogle Scholar
- Schlaeppi K, Dombrowski N, Oter RG, Ver Loren van Themaat E, Schulze-Lefert P. Quantitative divergence of the bacterial root microbiota in Arabidopsis thaliana relatives. Proc Natl Acad Sci. 2014;111:585–92.PubMedView ArticleGoogle Scholar
- Gil R, Latorre A, Moya A. Evolution of prokaryote-animal symbiosis from a genomics perspective. In: Hackstein JHP, editor. (Endo)symbiotic methanogenic Archaea, microbiology monographs. Berlin: Springer; 2010. p. 207–33.View ArticleGoogle Scholar
- Hacquard S, Garrido-Oter R, González A, Spaepen S, Ackermann G, Lebeis S, et al. Microbiota and host nutrition across plant and animal kingdoms. Cell Host Microbe. 2015;17:603–16.PubMedView ArticleGoogle Scholar
- Jones JDG, Dangl JL. The plant immune system. Nature. 2006;444:323–9.PubMedView ArticleGoogle Scholar
- Kau AL, Ahern PP, Griffin NW, Goodman AL, Gordon JI. Human nutrition, the gut microbiome and the immune system. Nature. 2011;474:327–36.PubMedPubMed CentralView ArticleGoogle Scholar
- Lebeis SL, Paredes SH, Lundberg DS, Breakfield N, Gehring J, McDonald M, et al. Salicylic acid modulates colonization of the root microbiome by specific bacterial taxa. Science. 2015;349:860–4.PubMedView ArticleGoogle Scholar
- Lee YK, Mazmanian SK. Has the microbiota played a critical role in the evolution of the adaptive immune system? Science. 2010;330:1768–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Brader G, Compant S, Mitter B, Trognitz F, Sessitsch A. Metabolic potential of endophytic bacteria. Curr Opin Biotechnol. 2014;27:30–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Hardoim PR, van Overbeek LS, Van Elsas JD. Properties of bacterial endophytes and their proposed role in plant growth. Trends Microbiol. 2008;16:463–71.PubMedView ArticleGoogle Scholar
- Weyens N, van der Lelie D, Taghavi S, Newman L, Vangronsveld J. Exploiting plant–microbe partnerships to improve biomass production and remediation. Trends Biotechnol. 2009;27:591–8.PubMedView ArticleGoogle Scholar
- Berendsen RL, Pieterse CMJ, Bakker PAHM. The rhizosphere microbiome and plant health. Trends Plant Sci. 2012;17:478–86.PubMedView ArticleGoogle Scholar
- Berg G, Grube M, Schloter M, Smalla K. Unraveling the plant microbiome: looking back and future perspectives. Front Microbiol. 2014;5:148.PubMedPubMed CentralGoogle Scholar
- Mendes R, Kruijt M, de Bruijn I, Dekkers E, van der Voort M, Schneider JHM, et al. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science. 2011;332:1097–100.PubMedView ArticleGoogle Scholar
- Turner TR, James EK, Poole PS. The plant microbiome. Genom Biol. 2013;14:1–10.View ArticleGoogle Scholar
- Bulgarelli D, Schlaeppi K, Spaepen S, Ver Loren van Themaat E, Schulze-Lefert P. Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol. 2013;64:807–38.PubMedView ArticleGoogle Scholar
- Shade A, McManus PS, Handelsman J. Unexpected diversity during community succession in the apple. MBio. 2013;4:1–12.View ArticleGoogle Scholar
- Nelson EB. Microbial dynamics and interactions in the spermosphere. Annu Rev Phytopathol. 2004;42:271–309.PubMedView ArticleGoogle Scholar
- Truyens S, Weyens N, Cuypers A, Vangronsveld J. Bacterial seed endophytes: genera, vertical transmission and interaction with plants. Environ Microbiol Rep. 2014;7:40–50.View ArticleGoogle Scholar
- Compant S, Kaplan H, Sessitsch A, Nowak J, Ait Barka E, Clément C. Endophytic colonization of Vitis vinifera L. by Burkholderia phytofirmans strain PsJN: from the rhizosphere to inflorescence tissues. FEMS Microbiol Ecol. 2008;63:84–93.PubMedView ArticleGoogle Scholar
- Gottel NR, Castro HF, Kerley M, Yang Z, Pelletier DA, Podar M, et al. Distinct microbial communities within the endosphere and rhizosphere of Populus deltoides roots across contrasting soil types distinct. Appl Environ Microbiol. 2011;77:5934–44.PubMedPubMed CentralView ArticleGoogle Scholar
- Inceoğlu O, Salles JF, van Overbeek L, van Elsas JD. Effects of plant genotype and growth stage on the betaproteobacterial communities associated with different potato cultivars in two fields. Appl Environ Microbiol. 2010;76:3675–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Weinert N, Piceno Y, Ding G-C, Meincke R, Heuer H, Berg G, et al. PhyloChip hybridization uncovered an enormous bacterial diversity in the rhizosphere of different potato cultivars: many common and few cultivar-dependent taxa. FEMS Microbiol Ecol. 2011;75:497–506.PubMedView ArticleGoogle Scholar
- Coleman-Derr D, Desgarennes D, Fonseca-Garcia C, Gross S, Clingenpeel S, Woyke T, et al. Plant compartment and biogeography affect microbiome composition in cultivated and native Agave species. New Phytol. 2016;209:798–811.PubMedView ArticleGoogle Scholar
- Fonseca-García C, Coleman-Derr D, Garrido E, Visel A, Tringe SG, Partida-Martínez LP. The cacti microbiome: interplay between habitat-filtering and host-specificity. Front Microbiol. 2016;7:150.
- Tardif S, Yergeau É, Tremblay J, Legendre P, Whyte LG, Greer CW. The willow microbiome is influenced by soil petroleum-hydrocarbon concentration with plant compartment-specific effects. Front Microbiol. 2016;7:1363.
- Ottesen AR, González Peña A, White JR, Pettengill JB, Li C, Allard S, et al. Baseline survey of the anatomical microbial ecology of an important food plant: Solanum lycopersicum (tomato). BMC Microbiol. 2013;13:114.PubMedPubMed CentralView ArticleGoogle Scholar
- Wagner MR, Lundberg DS, del Rio TG, Tringe SG, Dangl JL, Mitchell-Olds T. Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. Nat Commun. 2016;7:12151.PubMedPubMed CentralView ArticleGoogle Scholar
- Bodenhausen N, Horton MW, Bergelson J. Bacterial communities associated with the leaves and the roots of Arabidopsis thaliana. PLoS One. 2013;8:e56329.PubMedPubMed CentralView ArticleGoogle Scholar
- Hur M, Kim Y, Song H-R, Kim JM, Choi YI, Yi H. Effect of genetically modified poplars on soil microbial communities during the phytoremediation of waste mine tailings. Appl Environ Microbiol. 2011;77:7611–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Taghavi S, Garafola C, Monchy S, Newman L, Hoffman A, Weyens N, et al. Genome survey and characterization of endophytic bacteria exhibiting a beneficial effect on growth and development of poplar trees. Appl Environ Microbiol. 2009;75:748–57.PubMedView ArticleGoogle Scholar
- Ulrich K, Ulrich A, Ewald D. Diversity of endophytic bacterial communities in poplar grown under field conditions. FEMS Microbiol Ecol. 2008;63:169–80.PubMedView ArticleGoogle Scholar
- Sannigrahi P, Ragauskas AJ. Poplar as a feedstock for biofuels: a review of compositional characteristics. Biofuels, Bioprod, Bioref. 2010;4:209–26.View ArticleGoogle Scholar
- Searchinger T, Heimlich R, Houghton RA, Dong F, Elobeid A, Fabiosa J, et al. Use of U.S. croplands for biofuels increases greenhouse gases through emissions from land-use change. Science. 2008;319:1238–40.PubMedView ArticleGoogle Scholar
- Naylor RL, Liska AJ, Burke MB, Falcon WP, Gaskell JC, Rozelle SD, et al. The ripple effect: biofuels, food security, and the environment. Environ Sci Policy Sustain Dev. 2007;49:30–43.View ArticleGoogle Scholar
- Cassman KG. Climate change, biofuels, and global food security. Environ Res Lett. 2007;2:011002.View ArticleGoogle Scholar
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006;313:1596–604.PubMedView ArticleGoogle Scholar
- Tsai C-J, Xue L-J. CRISPRing into the woods. GM Crops Food. 2015;6:206–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Acker R, Leplé J-C, Aerts D, Storme V, Goeminne G, Ivens B, et al. Improved saccharification and ethanol yield from field-grown transgenic poplar deficient in cinnamoyl-CoA reductase. Proc Natl Acad Sci U S A. 2014;111:845–50.PubMedView ArticleGoogle Scholar
- Vanholme B, Cesarino I, Goeminne G, Kim H, Marroni F, Van Acker R, et al. Breeding with rare defective alleles (BRDA): a natural Populus nigra HCT mutant with modified lignin as a case study. New Phytol. 2013;198:765–76.PubMedView ArticleGoogle Scholar
- Harfouche A, Meilan R, Kirst M, Morgante M, Boerjan W, Sabatti M, et al. Accelerating the domestication of forest trees in a changing world. Trends Plant Sci. 2012;17:64–72.PubMedView ArticleGoogle Scholar
- Evans LM, Slavov GT, Rodgers-Melnick E, Martin J, Ranjan P, Muchero W, et al. Population genomics of Populus trichocarpa identifies signatures of selection and adaptive trait associations. Nat Genet. 2014;46:1089–96.PubMedView ArticleGoogle Scholar
- Liu D, Hu R, Palla KJ, Tuskan GA, Yang X. Advances and perspectives on the use of CRISPR/Cas9 systems in plant genomics research. Curr Opin Plant Biol. 2016;30:70–7.PubMedView ArticleGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.PubMedPubMed CentralView ArticleGoogle Scholar
- Dyall SD, Brown MT, Johnson PJ. Ancient invasions: from endosymbionts to organelles. Science. 2004;304:253–7.PubMedView ArticleGoogle Scholar
- Raven PH. A multiple origin for plastids and mitochondria: many independent symbiotic events may have been involved in the origin of these cellular organelles. Science. 1970;169:641–6.PubMedView ArticleGoogle Scholar
- Buée M, Reich M, Murat C, Morin E, Nilsson RH, Uroz S, Martin F. 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol. 2009;184:449–56.PubMedView ArticleGoogle Scholar
- Coddington JA, Agnarsson I, Miller JA, Kuntner M, Hormiga G, et al. Undersampling bias: the null hypothesis for singleton species in tropical arthropod surveys. J Anim Ecol. 2009;78:573–84.PubMedView ArticleGoogle Scholar
- Dickie IA. Letters. Insidious effects of sequencing errors on perceived diversity in molecular surveys. New Phytol. 2010;188:916–8.PubMedView ArticleGoogle Scholar
- Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM. Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007;8:R143.PubMedPubMed CentralView ArticleGoogle Scholar
- Huse SM, Welch DM, Morrison HG, Sogin ML. Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol. 2010;12:1889–98.PubMedPubMed CentralView ArticleGoogle Scholar
- Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, et al. Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc Natl Acad Sci. 2006;103:12115–20.PubMedPubMed CentralView ArticleGoogle Scholar
- Tedersoo L, Nilsson RH, Abarenkov K, Jairus T, Sadam A, Saar I. Methods. 454 Pyrosequencing and Sanger sequencing of tropical mycorrhizal fungi provide similar results but reveal substantial methodological biases. New Phytol. 2010;188:291–301.PubMedView ArticleGoogle Scholar
- Kunin V, Engelbrektson A, Ochman H, Hugenholtz P. Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. Environ Microbiol. 2010;12:118–23.PubMedView ArticleGoogle Scholar
- Quince C, Lanzén A, Curtis TP, Davenport RJ, Hall N, Head IM, et al. Accurate determination of microbial diversity from 454 pyrosequencing data. Nat Methods. 2009;6:639–41.PubMedView ArticleGoogle Scholar
- Reeder J, Knight R. Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat Methods. 2010;7:668–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Brooks JP, Edwards DJ, Harwich MD, Rivera MC, Fettweis JM, Serrano MG, et al. The truth about metagenomics: quantifying and counteracting bias in 16S rRNA studies. BMC Microbiol. 2015;15:66.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerasimidis K, Bertz M, Quince C, Brunner K, Bruce A, Combet E, et al. The effect of DNA extraction methodology on gut microbiota research applications. BMC Res Notes. 2016;9:365.PubMedPubMed CentralView ArticleGoogle Scholar
- Gans J, Wolinsky M, Dunbar J. Computational improvements reveal great bacterial diversity and high metal toxicity in soil. Science. 2005;309:1387–90.PubMedView ArticleGoogle Scholar
- Roesch LFW, Fulthorpe RR, Riva A, Casella G, Hadwin AKM, Kent AD, et al. Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 2007;1:283–90.PubMedPubMed CentralGoogle Scholar
- Medinger R, Nolte V, Pandey RV, Jost S, Ottenwälder B, Schlötterer C, et al. Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms. Mol Ecol. 2010;19:32–40.PubMedPubMed CentralView ArticleGoogle Scholar
- Shade A, Jones SE, Caporaso JG, Handelsman J, Knight R, Fierer N, et al. Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. MBio. 2014;5:1–9.View ArticleGoogle Scholar
- Bais HP, Weir TL, Perry LG, Gilroy S, Vivanco JM. The role of root exudates in rhizosphere interactions with plants and other organisms. Annu Rev Plant Biol. 2006;57:233–66.PubMedView ArticleGoogle Scholar
- Lugtenberg BJJ, Dekkers LC. Minireview. What makes Pseudomonas bacteria rhizosphere competent ? Environ Microbiol. 1999;1:9–13.PubMedView ArticleGoogle Scholar
- Lugtenberg B, Kamilova F. Plant-growth-promoting rhizobacteria. Annu Rev Microbiol. 2009;63:541–56.PubMedView ArticleGoogle Scholar
- Walker TS, Bais HP, Grotewold E, Vivanco JM. Update on root exudation and rhizosphere biology: root exudation and rhizosphere biology. Plant Physiol. 2003;132:44–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Compant S, Clément C, Sessitsch A. Plant growth-promoting bacteria in the rhizo- and endosphere of plants: their role, colonization, mechanisms involved and prospects for utilization. Soil Biol Biochem. 2010;42:669–78.View ArticleGoogle Scholar
- Schloss PD, Handelsman J. Toward a census of bacteria in soil. PLoS Comput Biol. 2006;2:e92.PubMedPubMed CentralView ArticleGoogle Scholar
- DeAngelis KM, Brodie EL, DeSantis TZ, Andersen GL, Lindow SE, Firestone MK. Selective progressive response of soil microbial community to wild oat roots. ISME J. 2009;3:168–78.PubMedView ArticleGoogle Scholar
- Grayston SJ, Vaughan D, Jones D. Rhizosphere carbon flow in trees, in comparison with annual plants: the importance of root exudation and its impact on microbial activity and nutrient availability. Appl Soil Ecol. 1997;5:29–56.View ArticleGoogle Scholar
- Hirano SS. Bacteria in the leaf ecosystem with emphasis on Pseudomonas syringae—a pathogen, ice nucleus, and epiphyte. Microbiol Mol Biol Rev. 2000;64:624–53.PubMedPubMed CentralView ArticleGoogle Scholar
- McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.PubMedPubMed CentralView ArticleGoogle Scholar
- Lindow SE, Brandl MT. Microbiology of the phyllosphere. Appl Environ Microbiol. 2003;69:1875–83.PubMedPubMed CentralView ArticleGoogle Scholar
- McCully ME. Niches for bacterial endophytes in crop plants: a plant biologist’s view. Aust J Plant Physiol. 2001;28:983–90.Google Scholar
- Fahlgren C, Hagström A, Nilsson D, Zweifel UL. Annual variations in the diversity, viability, and origin of airborne bacteria. Appl Environ Microbiol. 2010;76:3015–25.
- Castro HF, Classen AT, Austin EE, Norby RJ, Schadt CW. Soil microbial community responses to multiple experimental climate change drivers. Appl Environ Microbiol. 2010;76:999–1007.PubMedView ArticleGoogle Scholar
- Smit E, Leeflang P, Gommans S, Mil SVAN, Wernars K. Diversity and seasonal fluctuations of the dominant members of the bacterial soil community in a wheat field as determined by cultivation and molecular methods. Appl Environ Microbiol. 2001;67:2284–91.PubMedPubMed CentralView ArticleGoogle Scholar
- Ofek-Lalzar M, Sela N, Goldman-Voronov M, Green SJ, Hadar Y, Minz D. Niche and host-associated functional signatures of the root surface microbiome. Nat Commun. 2014;5:4950.PubMedView ArticleGoogle Scholar
- Peiffer JA, Spor A, Koren O, Jin Z, Tringe SG, Dangl JL, et al. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proc Natl Acad Sci. 2013;110:6548–53.PubMedPubMed CentralView ArticleGoogle Scholar
- Zarraonaindia I, Owens SM, Weisenhorn P, West K, Hampton-Marcell J, Lax S, et al. The soil microbiome influences grapevine-associated microbiota. MBio. 2015;6:e02527–e14.PubMedPubMed CentralView ArticleGoogle Scholar
- Albertsen M, Hugenholtz P, Skarshewski A, Nielsen KL, Tyson GW, Nielsen PH. Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat Biotechnol. 2013;31:533–8.PubMedView ArticleGoogle Scholar
- Ferrari B, Winsley T, Ji M, Neilan B. Insights into the distribution and abundance of the ubiquitous candidatus Saccharibacteria phylum following tag pyrosequencing. Sci Rep. 2014;4:3957.PubMedView ArticleGoogle Scholar
- Haichar FEZ, Marol C, Berge O, Rangel-Castro JI, Prosser JI, Balesdent J, et al. Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2008;2:1221–30.PubMedView ArticleGoogle Scholar
- Innerebner G, Knief C, Vorholt JA. Protection of Arabidopsis thaliana against leaf-pathogenic Pseudomonas syringae by Sphingomonas strains in a controlled model system. Appl Environ Microbiol [Internet]. 2011;77:3202–10.View ArticleGoogle Scholar
- Delmotte N, Knief C, Chaffron S, Innerebner G, Roschitzki B, Schlapbach R, et al. Community proteogenomics reveals insights into the physiology of phyllosphere bacteria. Proc Natl Acad Sci U S A. 2009;106:16428–33.PubMedPubMed CentralView ArticleGoogle Scholar
- Mark G, Morrissey JP, Higgins P, O’gara F. Molecular-based strategies to exploit Pseudomonas biocontrol strains for environmental biotechnology applications. FEMS Microbiol Ecol. 2006;56:167–77.PubMedView ArticleGoogle Scholar
- Patten CL, Glick BR. Role of Pseudomonas putida indoleacetic acid in development of the host plant root system. Appl Environ Microbiol. 2002;68:3795–801.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu X, Monchy S, Taghavi S, Zhu W, Ramos J, van der Lelie D. Comparative genomics and functional analysis of niche-specific adaptation in Pseudomonas putida. FEMS Microbiol Rev. 2011;35:299–323.
- Vorholt JA. Microbial life in the phyllosphere. Nat Rev Microbiol. 2012;10:828–40.PubMedView ArticleGoogle Scholar
- Sy A, Timmers ACJ, Knief C, Vorholt JA. Methylotrophic metabolism is advantageous for Methylobacterium extorquens during colonization of Medicago truncatula under competitive conditions. Appl Environ Microbiol. 2005;71:7245–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Custers R. First GM trial in Belgium since 2002. Nat Biotechnol. 2009;27:506.
- Feinstein LM, Sul WJ, Blackwood CB. Assessment of bias associated with incomplete extraction of microbial DNA from soil. Appl Environ Microbiol. 2009;75:5428–33.PubMedPubMed CentralView ArticleGoogle Scholar
- Op De Beeck M, Lievens B, Busschaert P, Declerck S, Vangronsveld J, Colpaert JV. Comparison and validation of some ITS primer pairs useful for fungal metabarcoding studies. PLoS One. 2014;9:e97629.View ArticleGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.PubMedPubMed CentralView ArticleGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.PubMedPubMed CentralView ArticleGoogle Scholar
- Simpson EH. Measurement of diversity. Nature. 1949;163:688–8.View ArticleGoogle Scholar
- Pielou EC. The measurement of diversity in different types of biological collections. J Theor Biol. 1966;13:131–44.View ArticleGoogle Scholar
- R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2011. http://www.R-project.org/. ISBN 3-900051-07-0.Google Scholar
- Clarke KR. Non-parametric multivariate analyses of changes in community structure. Austral Ecol. 1993;18:117–43.View ArticleGoogle Scholar
- De Cáceres M, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90:3566–74.PubMedView ArticleGoogle Scholar
- Letunic I, Bork P. Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011;39:475–8.View ArticleGoogle Scholar