Emergence Shapes the Structure of the Seed Microbiota - Applied and ...

7 downloads 48049 Views 4MB Size Report
Nov 12, 2014 - cedure of mothur at http://www.mothur.org/wiki/MiSeq_SOP (31). 16S ... procedure of the Fungal ITS analysis tutorial (GitHub, Inc.). Briefly, se- ..... during community succession in the apple flower microbiome. mBio4(2):.
Emergence Shapes the Structure of the Seed Microbiota Matthieu Barret,a,b,c Martial Briand,a,b,c Sophie Bonneau,a,b,c Anne Préveaux,a,b,c Sophie Valière,d,e Olivier Bouchez,d,f Gilles Hunault,g Philippe Simoneau,a,b,c Marie-Agnès Jacquesa,b,c INRA, UMR1345 Institut de Recherche en Horticulture et Semences, Beaucouzé, Francea; Agrocampus Ouest, UMR1345 Institut de Recherche en Horticulture et Semences, Beaucouzé, Franceb; Université d’Angers, UMR1345 Institut de Recherche en Horticulture et Semences, SFR4207 QUASAV, Beaucouzé, Francec; GeT-PlaGe, Genotoul, INRA Auzeville, Castanet-Tolosan, Franced; INRA, UAR1209, Département de Génétique Animale, INRA Auzeville, Castanet Tolosan, Francee; UMR INRA/INPT ENSAT/INPT ENVT, Génétique, Physiologie et Systèmes d’Élevage, INRA Auzeville, Castanet Tolosan, Francef; Université d’Angers, Laboratoire d’Hémodynamique, Interaction Fibrose et Invasivité Tumorale Hépatique, UPRES 3859, IFR 132, Angers, Franceg

Seeds carry complex microbial communities, which may exert beneficial or deleterious effects on plant growth and plant health. To date, the composition of microbial communities associated with seeds has been explored mainly through culture-based diversity studies and therefore remains largely unknown. In this work, we analyzed the structures of the seed microbiotas of different plants from the family Brassicaceae and their dynamics during germination and emergence through sequencing of three molecular markers: the ITS1 region of the fungal internal transcribed spacer, the V4 region of 16S rRNA gene, and a species-specific bacterial marker based on a fragment of gyrB. Sequence analyses revealed important variations in microbial community composition between seed samples. Moreover, we found that emergence strongly influences the structure of the microbiota, with a marked reduction of bacterial and fungal diversity. This shift in the microbial community composition is mostly due to an increase in the relative abundance of some bacterial and fungal taxa possessing fast-growing abilities. Altogether, our results provide an estimation of the role of the seed as a source of inoculum for the seedling, which is crucial for practical applications in developing new strategies of inoculation for disease prevention.

S

eeds are sexually derived structures of spermatophytes, which can, under appropriate conditions, germinate to produce new plants. Seeds, like other plant organs such as roots (1), leaves (2), and flowers (3), have evolved in association with a diverse microbial community, also known as the microbiota, which can play a role in plant growth and health. Indeed, numerous plant growthpromoting bacteria and phytopathogenic microorganisms have been isolated from a wide range of seeds using classical cultivation-based methods (4, 5). The significance of seed transmission of pathogens in the emergence of diseases in new planting areas has been recognized for decades (4). In consequence, the processes involved in the transmission of microorganisms from plant to seed have been documented mainly for phytopathogenic microorganisms. Three major pathways of transmission have been described for seedborne pathogens: (i) internal transmission through the vascular system, (ii) floral transmission by the stigma, and (iii) external transmission via contact of the seed with microorganisms present on fruits, flowers, or residues (6). According to the transmission pathway, seed-borne microorganisms can be located on the seed surface or imbedded in the tissue of the seed. While the internal transmission by the host xylem is probably restricted to pathogens or endophytes, many plant-associated microorganisms are potentially transmitted to the seed by the floral pathway. Indeed, the floral pathway allows the transmission of biocontrol microorganisms (7) and phytopathogenic bacteria in nonhost seeds (8, 9). Finally, the external pathway is probably the most permissive way of microorganism transmission from plant to seed, although very few data are currently available in the literature (10). Although seeds are carriers of multiple microorganisms, this does not necessary imply that the seed-borne microbes will colonize seedlings. These seed-borne microorganisms must display great physiological adaptation capacity to the changing conditions encountered during seed germination (11). By definition, germi-

February 2015 Volume 81 Number 4

nation starts immediately after seeds imbibe water and is completed when the radicle penetrate the structures that surround the embryo (11). During this physiological process, a range of nutrients are released in the soil surrounding the seed. The availability of these carbon compounds in the zone surrounding the germinating seed creates a new habitat called the spermosphere, which is a site of intense competition between seed-borne and soil microorganisms (5). While the transmission of microorganisms from seed to seedling is the primary source of inoculum for the plant, relatively few research groups have investigated the composition of the seed microbiota (12–15) and its dynamics during germination and emergence (16, 17). Pioneer studies have highlighted (i) that the endophytic bacterial communities are relatively well conserved from one generation to another (15) and (ii) that germinating seeds are mostly colonized by soil microorganisms (13, 17). However, the vast majority of these results have been obtained with surface-disinfected seeds, thus ignoring the influence of seed epiphytes as a source of inoculum. More recently, the composition of the seed-associated epiphytic microbiota of Brassica and Triticum

Received 12 November 2014 Accepted 4 December 2014 Accepted manuscript posted online 12 December 2014 Citation Barret M, Briand M, Bonneau S, Préveaux A, Valière S, Bouchez O, Hunault G, Simoneau P, Jacques M-A. 2015. Emergence shapes the structure of the seed microbiota. Appl Environ Microbiol 81:1257–1266. doi:10.1128/AEM.03722-14. Editor: H. L. Drake Address correspondence to Matthieu Barret, [email protected]. Supplemental material for this article may be found at http://dx.doi.org/10.1128 /AEM.03722-14. Copyright © 2015, American Society for Microbiology. All Rights Reserved. doi:10.1128/AEM.03722-14

Applied and Environmental Microbiology

aem.asm.org

1257

Barret et al.

has been investigated by sequencing a portion of cpn60 (18). Although this molecular marker has provided interesting information on bacterial and fungal taxa present at the seed surface, the lack of 16S rRNA gene and internal transcribed spacer (ITS) sequence data impaired comparison of these results with those of studies performed on other plant compartments. The production of seeds with optimal sanitary quality often depends on application of pesticides, including fungicides and insecticides. Since reducing pesticide use is a key objective for sustainable agriculture, alternative seed treatments have been designed. One of these alternatives is to perform seed coating with microorganisms having plant growth-promoting or biocontrol activities (19). Application of biocontrol microorganisms may have the potential to improve seedlings’ health by protecting them against seed-borne or soilborne pathogens (20). However, the performance of these biological treatments is generally inconsistent. Variations of efficiency can be partly explained by the activity of the seed-borne microbial community, which can sometimes limit the installation of exogenous microorganisms (21). Therefore, it is necessary to improve our knowledge of the nature, succession, and activities of seed-borne microorganisms during germination and emergence. The objective of this study was to investigate the dynamics of the seed microbiota during germination and emergence. To get insights into this dynamics, the compositions of the bacterial and fungal communities associated with 28 plants genotypes affiliated mostly to the Brassicaceae were evaluated at three physiological stages: seed, germinating seed, and seedlings. Sequencing of two molecular markers classically employed in microbial ecology, namely, the V4 region of 16S rRNA (22) and the ITS1 region of the fungal internal transcribed spacer (23), was performed on the MiSeq platform. Since 16S rRNA genes have relatively low sequence divergence among related bacterial taxa and are discriminant mostly at the genus level (24), an alternative bacterial marker was developed and employed in this work. This molecular marker is based on a portion of gyrB, a gene encoding the ␤ subunit of the DNA gyrase, which is frequently employed as a phylogenetic marker for many bacterial genera (25, 26). The different molecular markers used in this study provide valuable insights into the taxonomic composition of the seed microbiota. In addition, we identified key microbial taxa enriched during emergence, which could be promising candidates as seed inoculants. MATERIALS AND METHODS Experimental design. Twenty-eight seed samples (S01 to S28) were obtained from various plants belonging to different varieties, species, genera, and families (see Table S1 in the supplemental material). These seed samples were chosen to represent a large range of members of the family Brassicaceae. The structure of the seed microbiota on approximately 1,000 seeds per seed sample, as assessed by 1,000-seed weights, was studied. In parallel, 250 seeds of each sample were incubated in duplicate in sterile plastic boxes containing either crepe paper or blotter paper according to the standard germination methods of the International Seed Testing Association (ISTA). Plastic boxes were incubated at 20°C in obscurity for 24 and 96 h after imbibition in sterile distilled water. Seeds collected after 24 and 96 h of imbibition were defined as germinating seeds and seedlings, respectively. Although this is clearly an oversimplification of the physiological state of the seed, the 24 h time point (H24) was chosen as a proxy for germination since the radicles of most seed samples (including all the Brassicaceae) had emerged by that time. In addition, 96 h (H96)

1258

aem.asm.org

after imbibition corresponded to the apparition of the cotyledon for most seedlings and was thus chosen as a proxy for emergence. Sample preparation. A total of 1,000 seeds (H0), 200 germinating seeds (H24), and 100 seedlings (H96) of each sample were transferred in sterile tubes containing phosphate-buffered saline supplemented with 0.05% (vol/vol) of Tween 20. Samples were incubated either overnight at 4°C or for 2 h 30 min at room temperature according to standard protocols of ISTA (see Table S1 in the supplemental material). Suspensions were centrifuged (6,000 ⫻ g, 10 min, 4°C), and pellets were resuspended in approximately 2 ml of supernatant and transferred to Eppendorf tubes. Total genomic DNA was extracted from 84 different samples using the PowerSoil DNA isolation kit according to the manufacturer’s protocol. In addition, an artificial community sample was prepared by mixing equal amounts of genomic DNA from 15 bacterial strains (see Table S3 in the supplemental material). Bacterial strains were provided by the Collection for Plant-Associated Bacteria (CIRM-CFBP, IRHS, Beaucouzé, France). Genomic DNA of each bacterial strain was extracted with the Wizard genomic DNA purification kit according to the manufacturer’s protocol. All DNAs were pooled at an equimolar concentration with a final concentration of 20 ng · ␮l⫺1. Amplicon library construction and sequencing. Amplicon libraries were constructed following two rounds of PCR amplification. The first step was performed with the PCR primers 515f/806r (22) and ITS1F/ITS2 (27), which target the V4 region of 16S rRNA gene and ITS1, respectively. In addition, primers gyrB_aF64 (5=-MGNCCNGSNATGTAYATHGG3=) and gyrB_aR353 (5=-ACNCCRTGNARDCCDCCNGA-3=) were designed to amplify a portion of gyrB, which encodes subunit B of the bacterial gyrase. The primers’ binding sites correspond to Escherichia coli E22 (IMG taxon ID, 638341087) nucleotide positions 64 to 353 (see “gyrB sequence collection and analysis” below for further information). Forward and reverse primers carry the 5=-CTTTCCCTACACGACGCTCTT CCGATCT-3= and 5=-GGAGTTCAGACGTGTGCTCTTCCGATCT-3= tails, respectively. All PCRs were performed with a high-fidelity polymerase (AccuPrime Taq DNA polymerase system; Invitrogen) using the manufacturer’s protocol and 2 ␮l of environmental DNA (approximately 10 ng). The cycling conditions for 515f/806r and ITS1F/ITS2 were adapted from those described in references 22 and 27. Briefly, reaction mixtures were held at 94°C for 2 min, followed by 30 cycles of amplification at 94°C (30 s), 50°C (60 s), and 68°C (90 s), with a final extension step of 10 min at 68°C. Amplification of gyrB was performed as follows: 94°C (2 min) followed by 35 cycles of amplification at 94°C (30 s), 55°C (60 s), and 68°C (90 s), with a final extension step of 10 min at 68°C. All amplicons were purified with the Agencourt AMPure XP system and quantified with QuantIT PicoGreen. A second round of amplification was performed with 5 ␮l of purified amplicons and primers containing the Illumina adapters and indexes. PCR cycling conditions were 94°C (2 min), followed by 12 cycles of amplification (94°C for 1 min, 55°C for 1 min, 68°C for 1 min) and a final extension step at 68°C (10 min). All amplicons were purified and quantified as previously described. The purified amplicons were then pooled in equimolar concentrations, and the final concentration of the library was determined using a quantitative PCR (qPCR) next-generation sequencing (NGS) library quantification kit. Amplicon libraries were mixed with 10% PhiX control according to Illumina’s protocols. Three sequencing runs were performed for this study (see Table S3 in the supplemental material) with MiSeq reagent kit v2 (500 cycles). gyrB sequence collection and analysis. The prevalence of gyrB was investigated in 19,345 genomic sequences publicly available in the IMG database v4.3 (28) at the time of analysis. Coding sequences (CDSs) that exclusively belong to the protein family TIGR01059 were defined as GyrB homologs and retrieved for further analysis (18,572 hits found in 18,185 genomic sequences). These amino acid sequences were aligned with EINS-i of MAFFT v7.012 (29). The corresponding nucleotide sequences were aligned with TranslatorX (30) using the protein alignment. Following this step, conserved nucleotide blocks were visualized with Bioedit v7.2.5, and primers gyrB_aF64 and gyrB_aR353 were designed. The tax-

Applied and Environmental Microbiology

February 2015 Volume 81 Number 4

Dynamics of the Seed Microbiota

onomic resolution of the gyrB region targeted by the primers gyrB_aF64 and gyrB_aR353 was assessed as follows: nucleotide regions located between nucleotide positions 64 to 353 were selected in the gyrB alignment, primer sequences were trimmed, pairwise distances between all gyrB regions were computed, and sequences were grouped according to different genetic distances (from 0.01 to 0.10). These gyrB groups were then compared to cliques obtained with whole-genome-based average nucleotide identity (gANI) values, available at http://ani.jgi-psf.org. These gANI cliques can be used as a proxy for species delineation (N. Varghese, S. Mukherjee, N. Ivanova, K. Konstantinidis, K. Mavrommatis, N. Kyrpides, and A. Pati, unpublished data). At each genetic distance, the sensitivity, precision, and F1 score were calculated. Using these criterions, we defined the distance of 0.02 as the most representative of the species level, with an F1 score of 0.959 (see Fig. S2 in the supplemental material). Clustering MiSeq reads into OTUs. Raw reads were assembled in quality sequences using the steps described in the standard operating procedure of mothur at http://www.mothur.org/wiki/MiSeq_SOP (31). 16S rRNA gene and gyrB sequences were aligned against the 16S rRNA gene Silva alignment and a gyrB reference alignment, respectively. All sequences that did not align correctly were removed from the data sets. Chimeric sequences were detected with Uchime (32) and subsequently removed from the data set. Moreover, errors in coding sequences were assessed by translation of gyrB in amino acid sequences. Any sequence possessing a stop codon was discarded. Taxonomic affiliation of 16S rRNA gene and gyrB sequences was performed with a Bayesian classifier (33) (80% bootstrap confidence score) against the 16S rRNA gene training set (v9) of the Ribosomal Database Project (34) or against the gyrB database created with sequences retrieved from the IMG database (see “gyrB sequence collection and analysis” above). Unclassified sequences or sequences belonging to Eukaryota or Archaea, chloroplasts, or mitochondria were discarded. Sequences were divided into groups according to their taxonomic rank (level of order) and then assigned to operational taxonomic units (OTUs) at a 97% identity cutoff for 16S rRNA gene and 98% identity for gyrB sequences. ITS read pairs were initially combined in contigs with the command “make.contigs” of mothur v1.33, using the same parameters as those described for 16S rRNA gene and gyrB sequences. The variable ITS1 regions of ITS sequences were extracted with the Perl-based software ITSx (35). Then, sequences were processed using the Quantitative Insight Into Microbial Ecology (QIIME v1.7.0) software (36) according to the procedure of the Fungal ITS analysis tutorial (GitHub, Inc.). Briefly, sequences were clustered at a 97% identity cutoff using Uclust (37) and taxonomic affiliation was performed with a Bayesian classifier (33) (80% bootstrap confidence score) against the Unite database (38). Definition of a minimum threshold of 1‰ relative abundance for reproducible detection of OTUs. Incorporation of an artificial community sample in each sequencing run is useful to assess sequencing error rates (39). The artificial microbial community sequenced in this study is composed of genomic DNA from 15 bacterial isolates belonging to 13 distinct families (see Table S2 in the supplemental material). Therefore, in theory, the number of bacterial OTUs observed in this sample should be 13 after filtering of noisy sequences, removal of chimeras, and OTU picking. However, 183 and 281 OTUs were detected in the artificial community samples after analyses of sequencing runs 1 and 2, respectively (see Fig. S1 in the supplemental material). While the majority of these supernumerary OTUs were represented by few quality sequences in the artificial community samples, they were composed of numerous quality sequences in other sequenced samples. Thus, the additional OTUs observed in our artificial community sample are probably due to false assignation of index reads (40). Consequently, we used the artificial community sequence data sets obtained in sequencing runs 1 and 2 to define a threshold at which OTUs were considered of low abundance and then removed from the sample studied. Due to differences in the number of quality sequences per sample, the number of quality sequences per OTU was divided by the total size of the library (41). At a threshold of ⱖ1 ‰ of the library size, OTUs were defined as abundant OTUs (aOTUs) and conserved for further anal-

February 2015 Volume 81 Number 4

yses. Following this procedure, 14 and 13 aOTUs were observed in the artificial community samples sequenced in runs 1 and 2, respectively (see Fig. S1A and B in the supplemental material). Therefore, subsequent analyses of 16S rRNA gene, gyrB, and ITS were systematically performed on aOTUs. Microbial community analyses. A recent study has highlighted that the rarefying procedure classically performed in microbial ecology to normalize library size could undermine the performance of sample clustering (41). Since “omitted read counts added noise from the random sampling step” (41), we decided to present in this paper the data obtained with the normalization procedure related to the proportion. Both ␣ and ␤ diversity indexes were calculated with mothur v1.33 (42). Richness was assessed with the number of observed aOTUs, and diversity was assessed with Simpson’s inverse index. Beta diversity was assessed using the Bray-Curtis dissimilarity matrix (43). Nonmetric multidimensional scaling (NMDS) plots were generated for ␤-diversity analyses. Analysis of molecular variance (AMOVA) was performed to assess the effects of the different factors on the microbial community structure (P ⬍ 0.001). Differences in taxonomic abundance were highlighted with Krona radial space-filling displays (44). Differences in the relative abundance of aOTUs between the different factors were assessed with the R package edgeR (45). This type of analysis is usually employed to detect differential expressions of genes in wholetranscriptome shotgun sequencing (RNA-seq) data sets but has been recently proven to be effective for detection of enriched or depleted OTUs between treatments (41). Sequence counts were first normalized with the relative log expression (RLE) method (46), which is implemented in edgeR. Exact binomial tests corrected for multiples inferences with the Benjamini-Hochberg method (47) were then performed to detect differences in relative aOTU abundance between factors. aOTUs were defined as significantly enriched or depleted in one treatment with a corrected P value of ⬍0.001 and a log2-fold change in magnitude of ⱖ2. Correlations between aOTUs were calculated with the Sparse Correlations for Compositional data algorithm (SparCC) (48) implemented in mothur. The effect of uneven sampling was corrected by dividing sequence counts by total library size (proportion). Only correlations with values less than ⫺0.30 or larger than 0.30 were represented in the network using the R package qgraph (49). Nucleotide sequence accession number. All sequences have been deposited in the ENA database under the accession number ERP006367.

RESULTS

The structure of the seed microbiota was assessed in a total of 28 seed samples (S01 to S28) harvested from plants belonging to different varieties, species, genera, and families (see Table S1 in the supplemental material). Samples of 1,000 seeds (H0), 200 germinating seeds (after 24 h of water imbibition [H24]), and 100 seedlings (after 96 h of water imbibition [H96]) were used to study the dynamics of the composition of the seed microbiota during germination and emergence. Amplicon libraries of (i) the V4 region of 16S rRNA gene (22), (ii) the ITS1 region of the fungal ITS (23), and (iii) a portion of gyrB were sequenced in three independent runs using the Illumina MiSeq platform. Overall, 25,371,440 pairs of reads (see Table S2 in the supplemental material) from 85 samples corresponding to H0, H24, H96, and an artificial microbial community control (see Table S3 in the supplemental material) were obtained. Sequences corresponding to 16S and ITS amplicons were clustered into abundant operational taxonomic units (aOTUs) at ⱖ97% sequence identity, while gyrB amplicons were grouped at ⱖ98% sequence identity, which corresponded approximately to the bacterial species level (see Materials and Methods). Microbial ␣ diversity decreased during the transition from germinating seed to seedling. Bacterial and fungal richness levels were first assessed in seed samples. Median bacterial richness levels of 55 and 72 aOTUs were observed for 16S rRNA gene and gyrB

Applied and Environmental Microbiology

aem.asm.org

1259

Barret et al.

gyrB

H96

5 1

1

H0

H24

H96

H0

Time

H96

50 2 1

1 H24

H96

20

invsimpson

20 2

5

invsimpson

50 20 5 2 1

H0

H24 Time

50

Time

5

H24

50

aOTUs

500 5

50

aOTUs

500 50

aOTUs

5 1

H0

invsimpson

ITS1 500

16S rRNA gene

H0

Time

H24

H96

Time

H0

H24

H96

Time

FIG 1 Estimation of bacterial and fungal diversity. Richness (aOTUs) and diversity (Simpson’s inverse index [invsimpson]) were estimated in seeds (H0), germinating seeds (H24), and seedlings with 16S rRNA gene, gyrB, and ITS sequences. Each sample is represented by a green line, while the gray area represents the estimation of the distribution (created via Beanplot [71]).

sequences in the different seed samples (H0). However, important variations of bacterial richness were observed between seed samples (4 to 196 aOTUs for 16S rRNA gene sequences and 9 to 240 aOTUs for gyrB sequences) (Fig. 1). In comparison, fungal richness was less dispersed, ranging from 14 to 74 aOTUs with a median richness of 40 (Fig. 1). To assess if the observed variation was due to the plant genotype, the location of the seed production region, or the harvesting year, one-way analysis of variance (ANOVA) tests were performed. According to the analysis, none of these factors explained the variability of bacterial and fungal richness (P ⬎ 0.01) (see Fig. S1 in the supplemental material). Therefore, it is tempting to conclude that the variability of aOTU richness probably reflects heterogeneity between seed lots. The influence of germination and emergence on microbial diversity on germinating seed (H24) and seedling (H96), respectively, was then investigated. According to 16S rRNA gene sequences, bacterial richness significantly decreases during germination and emergence (one-way ANOVA, P ⬍ 0.001) (Fig. 1). However, gyrB sequences highlighted a significant decrease in aOTU richness only during the transition from germinating seed to seedling. The discrepancies observed with the two bacterial molecular markers employed in this work are probably due to differences in the level of taxonomic resolution obtained with OTU clustering. Indeed, at a cutoff of 98% identity, gyrB is species specific for more than 95% of the genomic sequences examined (see Materials and Methods and Fig. S2 in the supplemental material), while at a 97% identity cutoff, numerous 16S rRNA gene sequences are grouped at the genus level (24). A significant decline of bacterial diversity was also observed on seedlings with both bacterial molecular markers, although this decrease was observed earlier (H24) with 16S rRNA gene sequences (P ⬍ 0.001) (Fig. 1). Regarding fungal diversity, results obtained with ITS1 indicated a significant reduction in fungal richness and evenness during emergence (P ⬍ 0.001) (Fig. 1). Altogether, these

1260

aem.asm.org

results suggest that bacterial and fungal diversity decreases mainly during emergence. Emergence influences microbial ␤ diversity. Microbial community diversity between samples was estimated with Bray-Curtis dissimilarity. Nonmetric multidimensional scaling ordinations were used to plot a Bray-Curtis dissimilarity matrix obtained with 16S rRNA gene, gyrB, and ITS aOTUs (Fig. 2). NMDS plots indicate a spatial separation between the microbial communities associated with seedlings and the microbial communities associated with seeds and germinating seeds (Fig. 2A, B, and F). To test whether this observed clustering was statistically significant, AMOVA tests (50) were performed on each factor present in our sample collection (see Table S4 in the supplemental material). According to the magnitude of the F test statistic (Fs) values, the spatial separation observed between the fungal communities could be primarily explained by emergence (H96) and plant genotypes (at the family and genus level) (see Table S4 in the supplemental material). Regarding bacterial diversity, conclusions are less straightforward. Indeed, when bacterial community similarity was examined at the genus level with 16S rRNA gene sequences, samples were significantly clustered (P ⬍ 0.001) by physiological stage (H96) and location of the seed production region. However, samples were significantly clustered neither by emergence nor by production region at the bacterial species level with gyrB sequences (P ⬎ 0.001) (see Table S4 in the supplemental material). The seed microbiota is composed mainly of Gammaproteobacteria and Dothideomycetes. The composition of the microbial community was initially studied on seeds (H0). According to 16S rRNA gene sequences, bacterial aOTUs belonged mainly to Proteobacteria (13.1%, 5.8%, and 56.1% in the Alpha-, Beta- and Gammaproteobacteria classes, respectively), Firmicutes (11.3%), and Actinobacteria (9.1%) (Fig. 3). Despite differences in copy

Applied and Environmental Microbiology

February 2015 Volume 81 Number 4

Dynamics of the Seed Microbiota

FIG 2 Emergence influences the microbial ␤-diversity NMDS ordination of Bray-Curtis dissimilarity matrix obtained with 16S (A, B), gyrB (C, D), and ITS (E, F) aOTUs. Each dot represents a microbial community observed in samples derived from H0 (red), H24 (blue), and H96 (green). Stress values and total variance are indicated for each NMDS ordination.

number per genome, the bacterial community composition obtained with gyrB sequences is overall very similar to results obtained with 16S rRNA gene sequences, with the exception of Firmicutes (21.5%) (Fig. 3). The fungal community associated with seeds contained four main classes related to the Ascomycota phylum, i.e., Dothideomycetes (60.7%), Eurotiomycetes (5.5%), Leotiomycetes (6.3%), and Sordariomycetes (4.9%), and the Tremellomycetes class (15.5%), belonging to the Basidiomycota (Fig. 3). Few bacterial and fungal aOTUs were systematically detected in all seed samples (see Table S5 in the supplemental material). For instance, 3 aOTUs belonging to Pantoea (Otu00002), Pseudomonas (Otu00001), and Xanthomonas (Otu00007) were obtained with 16S rRNA gene sequences. This value dropped to one bacterial aOTU (Pantoea agglomerans) and one fungal aOTU (unclassified member of Mycosphaerellaceae) with species-specific markers (gyrB and ITS). Even the core community associated with 9 seed samples harvested from the same plant variety, namely, Brassica oleracea var. Capitata, was composed of only 7 fungal aOTUs and 6 bacterial aOTUs or 1 bacterial aOTU depending on the molecular marker employed (see Table S5 in the supplemental material). Enrichment of specific aOTUs during germination and emergence. Changes in the microbial community composition were observed mainly during emergence (H96). For instance, the relative abundance of Gammaproteobacteria increased to 85% and 83.2% with 16S rRNA gene and gyrB sequences, respectively (Fig. 3). This increase in relative abundance of Gammaproteobacteria was mainly due to the dominance of aOTUs related to the Pseudomonas genus, which represented approximately 57% and 49% of the total number of 16S rRNA gene and gyrB sequences, respectively (Fig. 3). Variations of the fungal community composition were also observed during emergence, with a simultaneous decrease of the relative abundance of Dothideomycetes (46.8%) and an increase in the relative abundance of Eurotiomycetes (17.7%).

February 2015 Volume 81 Number 4

This shift in fungal community composition was mostly due to an increase in the relative abundance of aOTUs related to the Penicillium genus (Fig. 3). To gain more insights into the dynamics of the microbial community during emergence, we investigated changes in the relative abundance of aOTUs between samples. Differences in the relative abundance of aOTUs were assessed with the R package edgeR (45) using exact binomial tests corrected for multiple inferences with the Benjamini-Hochberg method (47). Pairwise comparisons were performed for each physiological stage (H0, H24, and H96). aOTUs were defined as significantly enriched or depleted in one treatment at a corrected P value of ⬍0.001 and a log2-fold change magnitude of ⱖ2. Changes in the relative abundance of aOTUs were detected mostly in the comparison of H96 versus H0 results. Indeed, the relative abundances of 199 (16S rRNA gene), 106 (gyrB), and 167 (ITS) aOTUs were decreased during the transition from seeds to seedlings, confirming the decrease in bacterial and fungal richness previously observed (Fig. 1). In contrast, 13 (16S rRNA gene), 47 (gyrB), and 50 (ITS) aOTUs were significantly enriched in seedlings (see Table S6 in the supplemental material). Interestingly, bacterial aOTUs significantly enriched in the comparison of H96 with H0 were related to bacterial taxa (e.g., Massilia, Pantoea, or Pseudomonas) frequently found in plantassociated environments. In the same vein, several fungal aOTUs significantly enriched in seedlings compared to seeds correspond to ubiquitous cosmopolitan taxa (e.g., Penicillium, Chaetomium globosum, Rhizopus oryzae) often found in soil and on seed surfaces. These bacterial and fungal taxa have in common a fastgrowing ability (51, 52, 53, 54, 55, 56) that may explain their rapid colonization of seedlings. Cooperation between microbial taxa is frequent within the seed microbiota. Cooperation and competition between aOTUs were monitored by generating correlations networks with SparCC (48) on 16S rRNA gene, gyrB, and ITS sequences. Considering

Applied and Environmental Microbiology

aem.asm.org

1261

Barret et al.

FIG 3 Dynamics of microbiota composition during germination-emergence. Krona radial space-filling (44) charts show the mean relative abundances of bacterial and fungal taxa in seeds (H0), germinating seeds (H24), and seedlings (H96).

only correlations with values less than ⫺0.30 or greater than 0.30, we identified 236, 652, and 98 associations between 41 (16S rRNA gene), 146 (gyrB), and 33 (ITS) aOTUs. Positive associations between aOTUs were mostly observed (Fig. 4A, D, and G), suggesting weak competition and frequent cooperation between microorganisms. Overall, the cooccurrence of aOTUs across networks was not driven by the physiological stage of the sample, as aOTUs enriched during emergence were not systematically grouped together (Fig. 4B, E, and H). In contrast, clustering was observed for aOTUs affiliated to the same taxonomic class (Fig. 4C, F, and I).

1262

aem.asm.org

For instance, bacterial aOTUs affiliated to the class Bacilli were clustered in the networks drawn with 16S rRNA gene and gyrB sequences (Fig. 4C and F). Moreover, fungal aOTUs belonging to the Dothideomycetes were also grouped in a module of network (Fig. 4I). Taken together, these results suggest that cooperation occurs between aOTUs that belong to phylogenetically related taxa. DISCUSSION

Seeds harbor a complex microbial community, which may exert beneficial or deleterious effects on plant growth and health. To

Applied and Environmental Microbiology

February 2015 Volume 81 Number 4

Dynamics of the Seed Microbiota

FIG 4 Correlation networks observed between aOTUs. Correlation networks obtained with 16S rRNA gene (A, B, C), gyrB (D, E, F), and ITS (G, H, I) sequences. Nodes correspond to aOTUs, and connecting edges indicate correlations between them. Correlations between aOTUs were calculated with the Sparse Correlations for Compositional data algorithm (SparCC) (48) implemented in mothur. The effect of uneven sampling was corrected by dividing sequence counts by total library size. Only correlations with values less than ⫺0.30 (blue) or larger than 0.30 (orange) were represented in the network using the R package qgraph (49). While graphics in panels A, D, and G represent all the aOTUs, graphics in panels B, C, E, F, H, and I are restricted to aOTUs with negative or positive correlations. Blue and orange nodes (B, E, and F) represent aOTUs with decrease and increase in relative abundance during transition from seeds to seedlings. Node colors (C, F, and I) represent the different bacterial and fungal classes.

date, the structure of the seed microbiota has been explored mainly through culture-based diversity studies and therefore remains largely unknown. In this work, a comprehensive analysis of the seed microbiota and its dynamics during germination and emergence was performed through an amplicon sequencing approach. The amplicon sequencing approach relies on sequencing of a genomic region (also known as molecular markers) to identify the different species present in environmental samples. In microbial ecology, the most frequently employed markers are located in ribosomal genes (16S rRNA gene for Archaea and Bacteria) or in flanking regions (ITS for Fungi) since these regions are amplified with universal primers in a wide range of microorganisms (23, 57). However, variability of rRNA copy number per genome along with intragenomic polymorphisms within these regions may lead

February 2015 Volume 81 Number 4

to an overestimation of species diversity (23, 24). In addition, 16S rRNA gene sequences often fail to resolve bacterial species and are frequently limited to taxonomic affiliation at the genus level (24). Therefore, methods using alternative molecular markers (i.e., cpn60 and rpoB) present in all bacterial genomes as a single-copy gene have been successfully developed (58, 59). In this work, we have designed another universal bacterial marker based on a portion of gyrB, a gene encoding the ␤ subunit of the DNA gyrase that is frequently employed as a phylogenetic marker for numerous bacterial genera (25, 26). In comparison to cpn60 and rpoB markers, the region of gyrB can be successfully amplified with a single primer set in multiple bacterial phyla. Although some PCR bias probably affected the representation of the bacterial community composition, the results obtained at the phylum level with gyrB sequences are similar to those obtained with 16S rRNA gene se-

Applied and Environmental Microbiology

aem.asm.org

1263

Barret et al.

quences (Fig. 1). Moreover, the gyrB marker developed allows the affiliation of 18,185 genomic sequences to the species level with a sensitivity and precision of ⬎0.95 (see Materials and Methods; see also Fig. S2 in the supplemental material), which outperformed the taxonomic classification achieved with 16S rRNA markers (24, 57). In consequence, gyrB is useful to monitor changes in the relative abundance of bacterial species and affords an alternative tool to 16S rRNA gene for description of bacterial community composition. The microbial communities associated with seeds examined in this study are composed of approximately 70 bacterial and 50 fungal species. Although fungal richness is relatively constant across seed samples, strong variations in bacterial richness are observed with both molecular markers. Direct comparison of diversity indexes between studies is not straightforward since different molecular markers, sequencing technologies, and analysis workflows are employed. For instance, in order to decrease the false-assignation rate of index reads in MiSeq experiments, we have considered only aOTUs in our analysis. Despite discrepancies between different studies, the range of bacterial richness observed in our study seems to be in accordance with previous results obtained on Brassica, Spinacia oleracea, and Triticum seeds (16, 18). In addition, it is tempting to conclude that the seed microbiota contains on average fewer bacterial and fungal taxa than the microbial community associated with the rhizosphere (1, 60, 61) and a comparable level of diversity in comparison to the phyllosphere (2, 61). Differences in seed microbiota composition could be explained by the plant genotype itself (e.g., seed size, seed anatomy) but also by abiotic factors, such as field management practices, harvesting methods, seed processing, and storage (62). Based on our samples, the structure of the seed microbiota seems to be indeed driven by abiotic factors, such as the geographic location of the production region and the harvesting year. Similarly, the composition of the seed endophyte community of Zea and Oryza sativa is influenced mostly by soil types and water regime treatments (13, 15). In contrast, the effect of the host species is significant only on fungal community composition but does not seem to impact the structure of the bacterial community associated with seeds. Additional studies with experimental design dedicated to the evaluation of the effect of the host genotype, soil types, and anthropogenic activities on the structure of the seed microbiota will be necessary to further support these observations. According to our DNA profiling results, germination does not seem to affect the bacterial and fungal diversity. This is quite surprising since chemical properties drastically changed within germinating seeds as a result of intense exudation (11) and consequently one might expect a selection of copiotrophic microorganisms on germinating seeds. These results have now to be confirmed by RNA profiling approaches to specifically monitor active microbial populations and avoid the amplification of extracellular DNA. On the other hand, a major decrease of bacterial and fungal diversity is observed during transition from germinating seeds to seedlings, which is probably indicative of a strong selective force exerted by the young plant on seed-borne microorganisms. The decrease of diversity observed in our study is undoubtedly counterbalanced by the acquisition of additional microbial members when seeds are sown in soils (62). Indeed, the spermosphere of various germinating seeds is frequently colonized by soil bacteria (5, 12, 63). This soil inoculum along with airborne mi-

1264

aem.asm.org

croorganisms (64) will then compete with seed-borne microorganisms for colonization of different plant compartments. Overall, the seed microbiota is composed mainly of three bacterial phyla, namely, Actinobacteria, Firmicutes, and Proteobacteria, as well as two fungal classes, the Dothideomycetes and Tremellomycetes. Interestingly, these taxa are also frequently associated with the rhizosphere and the phyllosphere of various plant species (60, 61, 65), which suggests that the seed might provide an important source of microbial inoculum for other plant compartments, as already highlighted by Van Overbeek et al. (66). Based on gyrB and ITS sequences, one fungal aOTU affiliated to the Mycosphaerellaceae and one bacterial aOTU affiliated to Pantoea agglomerans (Otu00001) were ubiquitous on seeds examined in this study. Some strains of P. agglomerans possess plant growth-promoting activities and are therefore of interest for seed treatment (18). In order to have further insights on seed-borne bacterial and fungal taxa selected during germination and emergence, we analyzed changes in relative abundance of aOTUs through edgeR, a statistical method employed in RNA-seq analysis (45). We detected some aOTUs significantly enriched on seedlings that belong to Bacillus, Massilia, Pantoea, and Pseudomonas. These bacterial taxa are frequently encountered in several plant compartments, including seeds, leaves, and roots (18, 60, 65, 67). In addition, fungal aOTUs enriched during emergence, such as Trichoderma viride (68) and Chaetomium globosum (69), corresponded to biocontrol fungal species, which have the ability to colonize plant tissues both as epiphytes and as endophytes. All these microbial taxa enriched in seedlings are capable of rapid growth in response to increase in nutrient availability (52, 54, 55, 56), which may indicate that r-strategists (70) are selected during emergence. Enhanced colonization of seedlings by copiotrophic taxa suggests that competitive exclusion probably occurs between functional equivalent species. However, this hypothesis is not supported by analysis of correlation between aOTUs. Indeed, negative correlations between aOTUs were sparse, while frequent cooccurrences of aOTUs between microbial taxa were highlighted (Fig. 4). In summary, our detailed analysis of epiphytic and endophytic microbial communities associated with seeds through amplicon sequencing approaches revealed that the plant genotype had a strong effect on the dynamics of the seed microbiota during germination and emergence and selected key microbial taxa frequently associated with other plant compartments. Through the design of an alternative bacterial taxonomic marker based on a portion of gyrB, we were able to monitor changes in the relative abundance of bacterial species. This marker is a valuable alternative tool to 16S rRNA gene for the description of bacterial community composition. ACKNOWLEDGMENTS This research was supported in parts by grants awarded by the Region des Pays de la Loire (Qualisem, 2009 05369 and metaSEED, 2013 10080) and the European Commission (TESTA, FP7-KBBE-2012-6, 311875). We thank Steven Lindow for valuable comments on the manuscript, Thomas Baldwin (Vilmorin), Julia Buitink (IRHS), Régine Delourme (IGEPP), and Hubert Lybeert (HM-Clause) for providing seed samples, Amrita Pati for providing access to whole-genome-based average nucleotide identity values, Geraldine Taghouti, Perrine Portier, and CIRMCFBP (IRHS, UMR 1345 INRA-ACO-UA) for providing bacterial strains, and Muriel Bahut and Laurence Hibrand-Saint Oyant from the platform ANAN of SFR Quasav and the Genotoul Get-PlaGe sequencing facilities for their help on the MiSeq experiments.

Applied and Environmental Microbiology

February 2015 Volume 81 Number 4

Dynamics of the Seed Microbiota

REFERENCES 1. Philippot L, Raaijmakers JM, Lemanceau P, van der Putten WH. 2013. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol 11:789 –799. http://dx.doi.org/10.1038/nrmicro3109. 2. Rastogi G, Coaker GL, Leveau JH. 2013. New insights into the structure and function of phyllosphere microbiota through high-throughput molecular approaches. FEMS Microbiol Lett 348:1–10. http://dx.doi.org/10 .1111/1574-6968.12225. 3. Shade A, McManus PS, Handelsman J. 2013. Unexpected diversity during community succession in the apple flower microbiome. mBio 4(2): e00602-12. http://dx.doi.org/10.1128/mBio.00602-12. 4. Baker KF, Smith SH. 1966. Dynamics of seed transmission of plant pathogens. Annu Rev Phytopathol 4:311–332. http://dx.doi.org/10.1146 /annurev.py.04.090166.001523. 5. Nelson EB. 2004. Microbial dynamics and interactions in the spermosphere. Annu Rev Phytopathol 42:271–309. http://dx.doi.org/10.1146 /annurev.phyto.42.121603.131041. 6. Maude RB. 1996. Seedborne diseases and their control: principles and practice. CAB International, Wallingford, United Kingdom. 7. Spinelli F, Ciampolini F, Cresti M, Geider K, Costa G. 2005. Influence of stigmatic morphology on flower colonization by Erwinia amylovora and Pantoea agglomerans. Eur J Plant Pathol 113:395– 405. http://dx.doi.org /10.1007/s10658-005-4511-7. 8. Darrasse A, Darsonval A, Boureau T, Brisset M-N, Durand K, Jacques M-A. 2010. Transmission of plant-pathogenic bacteria by nonhost seeds without induction of an associated defense reaction at emergence. Appl Environ Microbiol 76:6787– 6796. http://dx.doi.org /10.1128/AEM.01098-10. 9. Darsonval A, Darrasse A, Meyer D, Demarty M, Durand K, Bureau C, Manceau C, Jacques MA. 2008. The type III secretion system of Xanthomonas fuscans subsp. fuscans is involved in the phyllosphere colonization process and in transmission to seeds of susceptible beans. Appl Environ Microbiol 74:2669 –2678. http://dx.doi.org/10.1128/AEM.02906-07. 10. Ngugi HK, Scherm H. 2006. Biology of flower-infecting fungi. Annu Rev Phytopathol 44:261–282. http://dx.doi.org/10.1146/annurev.phyto.44 .070505.143405. 11. Bewley JD. 1997. Seed germination and dormancy. Plant Cell 9:1055– 1066. http://dx.doi.org/10.1105/tpc.9.7.1055. 12. Buyer JS, Roberts DP, Russek-Cohen E. 1999. Microbial community structure and function in the spermosphere as affected by soil and seed type. Can J Microbiol 45:138 –144. http://dx.doi.org/10.1139/w98-227. 13. Hardoim PR, Hardoim CCP, van Overbeek LS, van Elsas JD. 2012. Dynamics of seed-borne rice endophytes on early plant growth stages. PLoS One 7:e30438. http://dx.doi.org/10.1371/journal.pone.0030438. 14. Ikeda S, Ytow N, Ezura H, Minamisawa K, Fujimura T. 2006. Discrimination of the commercial seeds of forage crops using ribosomal intergenic spacer analysis. Breed Sci 56:185–188. http://dx.doi.org/10.1270/jsbbs.56 .185. 15. Johnston-Monje D, Raizada MN. 2011. Conservation and diversity of seed associated endophytes in Zea across boundaries of evolution, ethnography and ecology. PLoS One 6:e20396. http://dx.doi.org/10.1371/journal .pone.0020396. 16. Lopez-Velasco G, Carder PA, Welbaum GE, Ponder MA. 2013. Diversity of the spinach (Spinacia oleracea) spermosphere and phyllosphere bacterial communities. FEMS Microbiol Lett 346:146 –154. http://dx.doi .org/10.1111/1574-6968.12216. 17. Ofek M, Hadar Y, Minz D. 2011. Colonization of cucumber seeds by bacteria during germination. Environ Microbiol 13:2794 –2807. http://dx .doi.org/10.1111/j.1462-2920.2011.02551.x. 18. Links MG, Demeke T, Grafenhan T, Hill JE, Hemmingsen SM, Dumonceaux TJ. 2014. Simultaneous profiling of seed-associated bacteria and fungi reveals antagonistic interactions between microorganisms within a shared epiphytic microbiome on Triticum and Brassica seeds. New Phytol 202:542–553. http://dx.doi.org/10.1111/nph.12693. 19. Deaker R, Roughley RJ, Kennedy IR. 2004. Legume seed inoculation technology—a review. Soil Biol Biochem 36:1275–1288. http://dx.doi.org /10.1016/j.soilbio.2004.04.009. 20. Berg G. 2009. Plant-microbe interactions promoting plant growth and health: perspectives for controlled use of microorganisms in agriculture. Appl Microbiol Biotechnol 84:11–18. http://dx.doi.org/10.1007/s00253 -009-2092-7. 21. Bacilio-Jimenez M, Aguilar-Flores S, del Valle MV, Perez A, Zepeda A,

February 2015 Volume 81 Number 4

22.

23.

24. 25. 26.

27.

28.

29. 30. 31.

32. 33.

34.

35.

36.

Zenteno E. 2001. Endophytic bacteria in rice seeds inhibit early colonization of roots by Azospirillum brasilense. Soil Biol Biochem 33:167–172. http://dx.doi.org/10.1016/S0038-0717(00)00126-7. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. 2011. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A 108(Suppl 1):S4516 –S4522. Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, Chen W, Bolchacova E, Voigt K, Crous PW, Miller AN, Wingfield MJ, Aime MC, An KD, Bai FY, Barreto RW, Begerow D, Bergeron MJ, Blackwell M, Boekhout T, Bogale M, Boonyuen N, Burgaz AR, Buyck B, Cai L, Cai Q, Cardinali G, Chaverri P, Coppins BJ, Crespo A, Cubas P, Cummings C, Damm U, de Beer ZW, de Hoog GS, Del-Prado R, Dentinger B, Dieguez-Uribeondo J, Divakar PK, Douglas B, Duenas M, Duong TA, Eberhardt U, Edwards JE, Elshahed MS, Fliegerova K, Furtado M, Garcia MA, Ge ZW, Griffith GW, Griffiths K, et al. 2012. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc Natl Acad Sci U S A 109:6241– 6246. http://dx.doi.org/10.1073/pnas.1117018109. Veˇtrovský T, Baldrian P. 2013. The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses. PLoS One 8:e57923. http://dx.doi.org/10.1371/journal.pone.0057923. Watanabe K, Nelson JS, Harayama S, Kasai H. 2001. ICB database: the gyrB database for identification and classification of bacteria. Nucleic Acids Res 29:344 –345. http://dx.doi.org/10.1093/nar/29.1.344. Yamamoto S, Harayama S. 1995. PCR amplification and direct sequencing of gyrB genes with universal primers and their application to the detection and taxonomic analyses of Pseudomonas putida strains. Appl Environ Microbiol 61:1104 –1109. Buee M, Reich M, Murat C, Morin E, Nilsson RH, Uroz S, Martin F. 2009. 454 Pyrosequencing analyses of forest soils reveal an unexpectedly high fungal diversity. New Phytol 184:449 – 456. http://dx.doi.org/10.1111 /j.1469-8137.2009.03003.x. Markowitz VM, Chen I-MA, Palaniappan K, Chu K, Szeto E, Grechkin Y, Ratner A, Jacob B, Huang J, Williams P, Huntemann M, Anderson I, Mavromatis K, Ivanova NN, Kyrpides NC. 2012. IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res 40:D115–D122. http://dx.doi.org/10.1093/nar/gkr1044. Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. http://dx.doi.org/10.1093/molbev/mst010. Abascal F, Zardoya R, Telford MJ. 2010. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res 291:1–7. http://dx.doi.org/10.1093/nar/gkq291. Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. 2013. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol 79:5112–5120. http://dx.doi.org/10 .1128/AEM.01043-13. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. 2011. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194 –2200. http://dx.doi.org/10.1093/bioinformatics/btr381. Wang Q, Garrity GM, Tiedje JM, Cole JR. 2007. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73:5261–5267. http://dx.doi.org/10.1128/AEM .00062-07. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-SyedMohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM. 2009. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res 37:D141–D145. http://dx.doi.org/10 .1093/nar/gkn879. Bengtsson-Palme J, Ryberg M, Hartmann M, Branco S, Wang Z, Godhe A, De Wit P, Sánchez-García M, Ebersberger I, de Sousa F, Amend A, Jumpponen A, Unterseher M, Kristiansson E, Abarenkov K, Bertrand YJK, Sanli K, Eriksson KM, Vik U, Veldre V, Nilsson RH. 2013. Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol Evol 4:914 –919. http://dx.doi.org /10.1111/2041-210X.12073. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters

Applied and Environmental Microbiology

aem.asm.org

1265

Barret et al.

37. 38.

39. 40. 41. 42.

43. 44. 45.

46. 47. 48. 49. 50. 51.

52. 53. 54.

WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. http://dx.doi.org/10.1038/nmeth.f.303. Edgar RC. 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460 –2461. http://dx.doi.org/10.1093 /bioinformatics/btq461. Abarenkov K, Henrik Nilsson R, Larsson KH, Alexander IJ, Eberhardt U, Erland S, Hoiland K, Kjoller R, Larsson E, Pennanen T, Sen R, Taylor AF, Tedersoo L, Ursing BM, Vralstad T, Liimatainen K, Peintner U, Koljalg U. 2010. The UNITE database for molecular identification of fungi–recent updates and future perspectives. New Phytol 186:281–285. http://dx.doi.org/10.1111/j.1469-8137.2009.03160.x. Schloss PD, Gevers D, Westcott SL. 2011. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS One 6:e27310. http://dx.doi.org/10.1371/journal.pone.0027310. Kircher M, Sawyer S, Meyer M. 2012. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res 40:e3. http://dx.doi.org/10.1093/nar/gkr771. McMurdie PJ, Holmes S. 2014. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol 10:e1003531. http: //dx.doi.org/10.1371/journal.pcbi.1003531. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537–7541. http://dx.doi.org/10.1128/AEM.01541-09. Bray JR, Curtis JT. 1957. An ordination of the upland forest communities of southern Wisconsin. Ecol Monographs 27:325–349. http://dx.doi.org /10.2307/1942268. Ondov B, Bergman N, Phillippy A. 2011. Interactive metagenomic visualization in a Web browser. BMC Bioinformatics 12:385. http://dx.doi .org/10.1186/1471-2105-12-385. Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26:139 –140. http://dx.doi.org/10.1093 /bioinformatics/btp616. Anders S, Huber W. 2010. Differential expression analysis for sequence count data. Genome Biol 11:R106. http://dx.doi.org/10.1186/gb-2010-11 -10-r106. Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B 57:289 –300. Friedman J, Alm EJ. 2012. Inferring correlation networks from genomic survey data. PLoS Comput Biol 8:e1002687. http://dx.doi.org/10.1371 /journal.pcbi.1002687. Epskamp S, Cramer AOJ, Waldorp LJ, Schmittmann VD, Borsboom D. 2012. qgraph: network visualizations of relationships in psychometric data. J Stat Software 48:1–18. Excoffier L, Smouse PE, Quattro JM. 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131:479 – 491. Battaglia E, Benoit I, van den Brink J, Wiebenga A, Coutinho PM, Henrissat B, de Vries RP. 2011. Carbohydrate-active enzymes from the zygomycete fungus Rhizopus oryzae: a highly specialized approach to carbohydrate degradation depicted at genome level. BMC Genomics 12:38. http://dx.doi.org/10.1186/1471-2164-12-38. Ofek M, Hadar Y, Minz D. 2012. Ecology of root colonizing Massilia (Oxalobacteraceae). PLoS One 7:e40117. http://dx.doi.org/10.1371 /journal.pone.0040117. Suryanarayanan TS. 2012. The diversity and importance of fungi associated with marine sponges. Botanica Marina 55:553–564. http://dx.doi.org /10.1515/bot-2011-0086. Thirup L, Ekelund F, Johnsen K, Jacobsen CS. 2000. Population dynamics

1266

aem.asm.org

55. 56. 57.

58.

59.

60.

61.

62. 63.

64.

65. 66.

67.

68.

69.

70.

71.

of the fast-growing sub-populations of Pseudomonas and total bacteria, and their protozoan grazers, revealed by fenpropimorph treatment. Soil Biol Biochem 32:1615–1623. http://dx.doi.org/10.1016/S0038-0717(00)00075-4. Zhang L, Birch RG. 1997. Mechanisms of biocontrol by Pantoea dispersa of sugar cane leaf scald disease caused by Xanthomonas albilineans. J Appl Microbiol 82:448 – 454. http://dx.doi.org/10.1046/j.1365-2672.1997.00135.x. Kim DM, Lee MH, Suh MK, Ha GY, Kim H, Choi JS. 2013. Onychomycosis caused by Chaetomium globosum. Ann Dermatol 25:232–236. http://dx.doi.org/10.5021/ad.2013.25.2.232. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO. 2013. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 41:e1. http://dx.doi.org/10.1093/nar/gks808. Links MG, Dumonceaux TJ, Hemmingsen SM, Hill JE. 2012. The chaperonin-60 universal target is a barcode for bacteria that enables de novo assembly of metagenomic sequence data. PLoS One 7:e49755. http: //dx.doi.org/10.1371/journal.pone.0049755. Vos M, Quince C, Pijl AS, de Hollander M, Kowalchuk GA. 2012. A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity. PLoS One 7:e30600. http://dx.doi.org/10.1371 /journal.pone.0030600. Bulgarelli D, Schlaeppi K, Spaepen S, Ver Loren van Themaat E, Schulze-Lefert P. 2013. Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol 64:807– 838. http://dx.doi.org/10.1146 /annurev-arplant-050312-120106. Meiser A, Balint M, Schmitt I. 2013. Meta-analysis of deep-sequenced fungal communities indicates limited taxon sharing between studies and the presence of biogeographic patterns. New Phytol 201:623– 635. http: //dx.doi.org/10.1111/nph.12532. Aleklett K, Hart M. 2013. The root microbiota—a fingerprint in the soil? Plant Soil 370:671– 686. http://dx.doi.org/10.1007/s11104-013-1647-7. Normander B, Prosser JI. 2000. Bacterial origin and community composition in the barley phytosphere as a function of habitat and presowing conditions. Appl Environ Microbiol 66:4372– 4377. http://dx.doi.org/10 .1128/AEM.66.10.4372-4377.2000. Maignien L, DeForce EA, Chafee ME, Eren AM, Simmons SL. 2014. Ecological succession and stochastic variation in the assembly of Arabidopsis thaliana phyllosphere communities. mBio 5:e00682-13. http://dx .doi.org/10.1128/mBio.00682-13. Vorholt JA. 2012. Microbial life in the phyllosphere. Nat Rev Microbiol 10:828 – 840. http://dx.doi.org/10.1038/nrmicro2910. van Overbeek LS, Franke AC, Nijhuis EH, Groeneveld RM, da Rocha UN, Lotz LA. 2011. Bacterial communities associated with Chenopodium album and Stellaria media seeds from arable soils. Microb Ecol 62:257– 264. http://dx.doi.org/10.1007/s00248-011-9845-4. Green SJ, Michel FC, Jr, Hadar Y, Minz D. 2007. Contrasting patterns of seed and root colonization by bacteria from the genus Chryseobacterium and from the family Oxalobacteraceae. ISME J 1:291–299. http://dx.doi .org/10.1038/ismej.2007.33. Bruez E, Vallance J, Gerbore J, Lecomte P, Da Costa J-P, GuerinDubrana L, Rey P. 2014. Analyses of the temporal dynamics of fungal communities colonizing the healthy wood tissues of esca leafsymptomatic and asymptomatic vines. PLoS One 9:e95928. http://dx.doi .org/10.1371/journal.pone.0095928. Longoni P, Rodolfi M, Pantaleoni L, Doria E, Concia L, Picco AM, Cella R. 2012. Functional analysis of the degradation of cellulosic substrates by a Chaetomium globosum endophytic isolate. Appl Environ Microbiol 78: 3693–3705. http://dx.doi.org/10.1128/AEM.00124-12. Reznick D, Bryant MJ, Bashey F. 2002. r- and K-selection revisited: the role of population regulation in life-history evolution. Ecology 83:1509 – 1520. http://dx.doi.org/10.1890/0012-9658(2002)083[1509:RAKSRT]2.0 .CO;2. Kampstra P. 2008. Beanplot: a boxplot alternative for visual comparison of distributions. J Stat Software 28:1–9.

Applied and Environmental Microbiology

February 2015 Volume 81 Number 4