Exploring the plant-associated bacterial communities in Medicago sativa L
© Pini et al.; licensee BioMed Central Ltd. 2012
Received: 5 August 2011
Accepted: 20 May 2012
Published: 20 May 2012
Skip to main content
© Pini et al.; licensee BioMed Central Ltd. 2012
Received: 5 August 2011
Accepted: 20 May 2012
Published: 20 May 2012
Plant-associated bacterial communities caught the attention of several investigators which study the relationships between plants and soil and the potential application of selected bacterial species in crop improvement and protection. Medicago sativa L. is a legume crop of high economic importance as forage in temperate areas and one of the most popular model plants for investigations on the symbiosis with nitrogen fixing rhizobia (mainly belonging to the alphaproteobacterial species Sinorhizobium meliloti). However, despite its importance, no studies have been carried out looking at the total bacterial community associated with the plant. In this work we explored for the first time the total bacterial community associated with M. sativa plants grown in mesocosms conditions, looking at a wide taxonomic spectrum, from the class to the single species (S. meliloti) level.
Results, obtained by using Terminal-Restriction Fragment Length Polymorphism (T-RFLP) analysis, quantitative PCR and sequencing of 16 S rRNA gene libraries, showed a high taxonomic diversity as well as a dominance by members of the class Alphaproteobacteria in plant tissues. Within Alphaproteobacteria the families Sphingomonadaceae and Methylobacteriaceae were abundant inside plant tissues, while soil Alphaproteobacteria were represented by the families of Hyphomicrobiaceae, Methylocystaceae, Bradyirhizobiaceae and Caulobacteraceae. At the single species level, we were able to detect the presence of S. meliloti populations in aerial tissues, nodules and soil. An analysis of population diversity on nodules and soil showed a relatively low sharing of haplotypes (30-40%) between the two environments and between replicate mesocosms, suggesting drift as main force shaping S. meliloti population at least in this system.
In this work we shed some light on the bacterial communities associated with M. sativa plants, showing that Alphaproteobacteria may constitute an important part of biodiversity in this system, which includes also the well known symbiont S. meliloti. Interestingly, this last species was also found in plant aerial part, by applying cultivation-independent protocols, and a genetic diversity analysis suggested that population structure could be strongly influenced by random drift.
Similar to the intensively studied animal microbioma, plants harbor a wide range of diverse bacteria forming a complex biological community, which includes pathogens, mutualists (symbionts), and commensals [1, 2]. Depending on the colonized compartment, these bacteria are rhizospheric (root colonizers), endophytic (colonizing the endosphere, the bulk of internal tissues) and phyllospheric or epiphytic (leaf or stem surface). In recent years plant-associated bacteria (endophytic, epiphytic and rhizospheric) have been widely studied, mainly as promising tools for biotechnological applications [3–7], but investigations have also been carried out on the ecology and taxonomy of plant-associated bacterial communities [8–11]. Despite a high taxonomic diversity, only few bacterial taxa have been found characteristically associated to the majority of plant species, notably members of the Alphaproteobacteria class [2, 7, 8, 12, 13]. Consequently, the generally accepted idea is that the ability to colonize a plant is not a common, widespread feature present in the soil bacterial community, but preferentially resides in specific taxa which may be considered more ecologically versatile or genetically prone to the association with plants. This last hypothesis has recently been supported by the finding that, at least in the class of Alphaproteobacteria, a common gene repertoire seems to be present in all of its plant-associated members .
Medicago sativa L. (alfalfa) is one of the most important legume crop in temperate areas throughout the world, commonly used as forage or in crop rotation practices to contribute organic nitrogen to the soil via its symbiosis with the nitrogen fixing bacteria . Moreover, it is important also for bioenergy production  and is one of the most suited plant species for land restoration . Finally, this species, and the diploid relative M. truncatula Gaertn. (barrel medic), are among the most studied model species regarding the molecular aspects of plant-bacteria symbiosis, particularly in relation with the alphaproteobacterium Sinorhizobium (syn. Ensifer) meliloti[18–20]. Concerning S. meliloti, this species is present in most temperate soils, and, when conditions are suitable, it forms specialized structures, called nodules, in the roots of alfalfa plants where it differentiates into bacteroids . It is assumed that a fraction of bacterial cells is released from dehiscent nodules to soil, giving rise to new free-living rhizobial clones . In the last years S. meliloti has been found able to also endophytically colonize the aerial part of other plant species, as rice , suggesting the presence of several ecological niches for this species (soil, nodule, other plant tissues).
While the plant-associated bacterial flora of M. sativa has never been investigated at the community level, S. meliloti population genetics have been extensively studied in the past [23–28], but only on strains isolated from nodules, with a few early studies performed on bacteria directly recovered from soil [29, 30], due to the lack of efficient selective culture media. No data have been reported on the presence in natural conditions of S. meliloti as endophytes in other plant compartments (such as leaves) and no comparison of soil vs. plant-associated populations has been done.
Based on the above mentioned considerations, there is a need to characterize the bacterial community associated with M. sativa in relation to both the potentially important role the class of Alphaproteobacteria seems to have as main component of a “core plant-associated bacterial community” in several different plant species [13, 31–33], and to the relationships of soil vs. plant-associated populations of the symbiotic alphaproteobacterial partner S. meliloti.
In this work we investigated the bacterial communities associated with the legume M. sativa, focusing on both the total bacterial community composition and on the presence and populations structure of the symbiotic partner S. meliloti in soil and plant tissues.
The analysis was conducted by cultivation-independent techniques on alfalfa (M. sativa) plants grown in mesocosm pots. The bacterial community associated with M. sativa and that of the surrounding soil were analyzed at high (class, family) and low (single species, S. meliloti) taxonomic levels by employing Terminal-Restriction Fragment Length Polymorphism (T-RFLP) profiling , 16 S rRNA library screening and S. meliloti-specific markers [34, 35]. These approaches allowed us to explore for the first time the bacterial community composition of such important plant species and the populations of S. meliloti without cultivation.
The ribotype variability of bacterial communities present in soil and associated to plant tissues (nodules, stems and leaves) was investigated by T-RFLP analysis. A total of 43 samples was analyzed: in particular one pooled soil sample for each one of the three pots, one pooled sample from all the nodules found in each pot and four plants per pot (one stem and 2–3 pools of leaves per plant). T-RFLP profiles on these samples produced 253 Terminal-Restriction Fragments (T-RFs) or ribotypes after the restriction digestion with two restriction enzymes, HinfI and TaqI. 16 S rRNA gene amplification and T-RFLP profiling was also performed on DNA extracted from surface-sterilized seeds, but no bands of 16 S rRNA gene amplification were recovered (data not shown), suggesting a very low bacterial titre in seeds.
Interestingly, stem and leaf communities showed a significant (P < 0.0001), though small (pairwise F ST = 0.05) separation (Additional file 2: Table S2). Moreover, AMOVA on stems and leaves community revealed a statistically significant differentiation between the three pots (P < 0.0001), irrespective of possible grouping (either plant genotype-related or unrelated), suggesting a pot-effect over the taxonomic shaping of the leaf-associated community and no effect of plant genotypes. These data confirmed a previous long-term experiment only addressing S. meliloti species .
T-RFLP analysis has shown that bacterial communities clustered in three groups (soil, nodules and plant aerial parts). In order to elucidate which taxa are mainly represented in the bacterial communities of these assemblies, three 16 S rRNA gene clone libraries were constructed pooling together the DNAs extracted from the samples of each environment; additionally, leaves and stems samples were also pooled together due to their high similarities as mentioned above. Pooled samples did conceivably result in an enrichment of the more shared taxa possibly preventing the detection of taxa associated only to a few individual samples. DNA was used as template to construct three 16 S rRNA libraries; a total of 276 clones (from 78 to 116 per library) were sequenced. Sequence analysis revealed, as expected, that the soil community was the most diverse (Shannon H’ = 4.63; Chao1 = 168), while the nodule-associated community was less diverse (Shannon H’ = 1.98; Chao1 = 30), (Additional file 3: Table S3). As a consequence, the library of nodules showed a coverage (85.9%) higher than those of stems + leaves (74.1%) and soil (47.1%).
Titers of S. meliloti in soil and plant tissues§
4.92 ± 2.82 x 104
2.78 ± 0.63 x 104
3.07 ± 0.67 x 109
4.25 ±1.24 x 108 **
2.73 ± 1.21 x 104
3.22 ±2.4 x 103 *
8.65 ± 4.04 x 103
4.28 ± 1.23 x 103
1.16 ± 0.33 x 104
2.88 ± 1.09 x 104
1.20 ± 0.50 x 1010
1.01 ± 0.10 x 109 **
2.37 ± 0.49 x 103
1.13 ± 0.15 x 103
9.74 ± 5.08 x 102
2.34 ±0.78 x 102
2.70 ± 0.41 x 105
7.42 ±0.93 x 104 *
6.02 ± 1.45 x 109
2.02 ± 3.22 x 107 **
4.91 ± 0.95 x 105
1.07 ± 3.74 x 105
5.54 ± 2.83 x 103
5.21 ± 3.01 x 103
In recent years there has been an increasing interest in exploring the bacterial flora associated with plants [37–41]. A recent field survey indicates  that plant aerial parts (leaves) harbor complex, but highly variable, bacterial communities, and that only a small number of bacterial taxa (mainly belonging to Alphaproteobacteria) are plant-specific. In the experiments reported here, as in the majority of the reports on endophytic microflora, we refer to endophytic and epiphytic bacteria indicating all those that are inside the plant tissue or strongly adhering to the plant surface, such as they resist washing and sterilization (or their DNA is retained by plant tissue), therefore a more correct definition could be “plant-associated bacteria”.
The present study shows that root nodules and aerial parts of Medicago sativa plants grown in mesocosm conditions, harbor distinct bacterial communities with specific signatures at the class, family and species levels and that these communities do not mirror soil bacterial communities.
Initially, T-RFLP profiles allowed us to show that bacterial communities present in the different environments (soil, nodules, stems and leaves) were strongly differentiated and in particular that a large heterogeneity was present between leaves of individual plants, though soil profiles were highly similar. Moreover a clear separation between above-ground (stem and leaves) and below-ground environments (soil and nodules) was detected. An analysis of the clone libraries, prepared from above-ground and below-ground pooled samples, revealed an uneven distribution of bacterial classes, with a marked pattern highlighting the class of Alphaproteobacteria as the more abundant in plant tissues (this class represented half of the clones in the stem + leaf library). The same uneven pattern was then observed, at lower taxonomic ranks, within the Alphaproteobacteria, with sequences of clones belonging to members of the Methylobacteriaceae and Sphingomonadaceae families being more abundant in stem than in soil and nodules. Methylobacteria and Sphingomonadaceae have been found as endophytes in a number of plants [8, 12, 31, 33, 42–45] and it is believed that this group of bacteria may take advantage from living as plant-associated, thanks to its ability to utilize the one-carbon alcohol methanol discharged by wall-associated pectin metabolism of growing plant cells.
Concerning root nodule bacterial communities, obtained data indicated that very diverse bacterial taxa are associated with nodules, the most represented being the specific rhizobial host of M. sativa, the alphaproteobacterium S. meliloti. However, additional taxa have been found, including members of Actinobacteria Flavobacteria Gammaproteobacteria and Betaproteobacteria, which may have some additional plant growth-promoting activities (see for instance [46, 47]).
In soil, Acidobacteria was one of the most important divisions (in terms of number of clones in the library) and was present exclusively in the soil clone library, in agreement with many previous observations [48, 49]. A relatively high presence of Archaea (Thermoprotei) was also found. Checking the 16 S rRNA gene sequences present in the Ribosomal Database for 799f/pHr primer annealing, we found that PCR amplification from Thermoprotei was theoretically possible with this primer pair (data not shown). The presence of Archaea in the soil is not unexpected  and could be linked also to the anoxic or nearly anoxic conditions present in the bottom of the pot. However, since the low coverage of soil clone library, the presence of many other additional taxa, as well of different proportions of those found here cannot be excluded. In addition, it should be mentioned that differences between soil and plant-tissues bacterial communities could also be ascribed to the different DNA extraction protocols we were obliged to use, since a unique protocol (bead-beading protocol for both soil DNA and plant DNA) failed in a successful extraction of DNA from both soil and plant tissues (data not shown). A similar technical need was encountered by other authors also , which renders the study of the relationships between plant-associated and soil bacterial communities still at its beginning.
At the lowest taxonomic rank here investigated, within the species S. meliloti, we detected the presence of this species in all environment analyzed (soil, nodules and plant aerial tissues). This finding is confirming earlier reports on the ability of S. meliloti to behave as an endophytic strain, colonizing all plant compartments, besides being a root symbiont of legumes , and suggest a potential higher genetic variability of S. meliloti population, and, from the other side, potential new ecological and functional roles for this species, not investigated so far[29, 51, 52]. Unfortunately, the low population size of S. meliloti in stems and leaves and the possible presence of PCR inhibitors (plant DNA or phenolic compounds, for instance) did not permit the amplification of 16 S-23 S rRNA intergenic region from plant aerial parts to obtain information about the genetic diversity and structure of S. meliloti population resident in plant aerial part. No hypothesis could then be drawn about the relationships between this population and those of soil and nodules. Concerning S. meliloti populations present in soil and nodules, similar values for diversity were detected in nodules and in soil, suggesting that both environments harbor a consistent fraction of the population’s genetic diversity. Interestingly, most of the T-RFs were detected in one sample only, and a very small fraction of T-RFs was shared among all samples, though the original soil material was homogeneous and should, in theory, contain the same S. meliloti haplotypes. Therefore, S. meliloti populations from all the three mesocosms investigated were highly differentiated between each other and, as expected from previous studies on S. meliloti and on Bradyrhizobium, no statistically significant plant genotype- related haplotypes were detected.
A possible explanation of such findings could be linked to the relatively low titers of S. meliloti in soil (104-105 cells/g), which is roughly 1/10,000 of the total bacterial community of soil (estimated at ~109 16 S rRNA gene copies/g of soil by qPCR, data not shown). Such estimated S. meliloti titers were similar to those previously observed in other soil and plant tissues  and in line with those normally found in soil with viable (Most Probable Number, MPN) estimates [26, 54]. As a consequence of this low population size, founder effect and genetic drift are likely to be among the main shaping forces of S. meliloti population in this experimental set-up, perhaps permitting the fixation of sample-specific haplotypes by simple chance . Regarding the nodule-soil relationships, though our experiments did not directly address this issue, the reported S. meliloti population analysis suggests the presence of somewhat nonoverlapping soil and nodule population fractions, even if no specific patterns of soil and nodule populations were detected. The presence of different rhizobial haplotypes in nodules and soil was previously found in chickpea  and clover , though no simple conclusion could be drawn, because of limited sampling. However, as for total bacterial community analysis, it should be mentioned that the use of two different DNA extraction protocols for soil and plant DNA may have produced some bias in the proportion of the different haplotyes detected.
In conclusion, we show on M. sativa that its associated microflora, though highly variable, is mainly related to the presence of Alphaproteobacteria. This class has an uneven presence of families in stems + leaves, nodules and soil. We then speculated that a sort of “pan-plant-associated bacterial community” may be composed of a large plethora of “accessory” taxa, which are occasionally associated with plants, and a small number of “core” taxa (e.g. Alphaproteobacteria families) which, on the contrary, are consistently found in the plants. Moreover, within Alphaproteobacteria the specific alfalfa symbiotic species S. meliloti, abundant as symbiont in root nodules, was also detected in soil and in leaves, with potentially different populations, suggesting a more complex interplay of colonization of multiple environments (soil, root nodules, other plant tissues) by this species.
A controlled experiment was set-up in mesocosms composed of three pots (numbered 1, 2, 3) containing Medicago sativa (alfalfa) plants grown at CRA-FLC Lodi, Italy, in outdoor conditions. Two of the three pots were planted with the same line of alfalfa (1x5) while the third pot was planted with a different line (5x7). The pots (cylinders of 25 cm diameter x 80 cm depth) with a drainage layer on the bottom, were filled with a sandy loam non-calcareous soil (57.8% sand, 32% silt, 10.2% clay, 1.7% organic matter and 0.09% total N; pH 6.7) in which alfalfa has never been grown. Phosphorus and potassium equivalent to 120 Kg ha-1 of P2O5 and 180 Kg ha-1 of K2O were distributed into the soil, while no mineral N was added; irrigation was not limiting. Twenty plants/pot (density equivalent to 400 plants m-2) were transplanted in March 2008 and allowed to grow until the 2nd year (the end of September 2009), when plant aerial parts of 12 plants were harvested and the pots were opened to allow sampling of the whole eye-detectable nodules present (approximately 80–100 of various sizes per pot) and of bulk soil. Roots were excluded from the analysis since the presence of small nodules or nodule primordia could not be excluded, possibly inducing a strong bias in the estimation of “non-nodule-associated root colonizers”. The plant sample size was chosen on the basis of a previous analysis of plant-by-plant variation in which the overall diversity of communities did not change from 2 to 30 plants (unpublished data and ). Stems, leaves (pools of around 10 leaves per plant) and nodules were washed with water and with 10 mM MgSO4 twice to remove most soil and dust particles and eliminate bacteria loosely adhering to the surface and then surface sterilized with 1% HClO for 1 min. Samples of soil, nodules, stem and leaves were then stored at −80°C from 1–2 weeks before DNA extraction.
A control of seed-borne bacteria was also prepared with seeds of M. sativa surface sterilized with 1% HgCl2.
S. meliloti viable titres in sterilized nodules have been estimated by serial dilution of crushed nodules as previously reported .
DNA was extracted from soil by using a commercial kit (Fast DNA Spin kit for soil, QBiogene, Cambridge, UK) following the manufacturer’s instructions. DNA extraction from plant tissues and surface sterilized control seeds was performed by a 2X CTAB protocol as previously described . The 16 S rRNA gene pool of total bacterial community was amplified from the extracted DNA with primer pairs 799f (labeled with HEX) and pHr which allow the amplification of most bacterial groups without targeting chloroplast DNA . PCR conditions and Terminal-Restriction Fragment Length Polymorphism (T-RFLP) profiling were as previously reported , by using HinfI and TaqI restriction enzymes. For sinorhizobial populations, T-RFLP was carried out on 16 S-23 S ribosomal intergenic spacer amplified from total DNA (IGS-T-RFLP) with S. meliloti specific primers and AluI and HhaII restriction enzymes, as already reported . Real-Time PCR (qPCR) for quantification of S. meliloti DNA was carried out on rpoE1 and nodC loci, as previously reported ; two different calibration curves were constructed, one for soil samples and the other one for plant samples, by using as template DNA extracted from sterile soil (without presence of S. meliloti) and from sterile plant (grown in petri dishes), both spiked with serial dilutions of known titres of S. meliloti cells, as previously reported . Controls with S. medicae WSM419 DNA were included in both IGS-T-RFLP and qPCR, for S. meliloti species-specificity check .
Amplified (with 799f and pHr primer pair) 16 S rRNA genes from DNA extracted from soil, nodules, pooled stems and leaves of a 1:1:1 mix of all pots were inserted into a pGemT vector (Promega, Fitchburg, WI, USA) and cloned in E. coli JM109 cells. Positive clones were initially screened by white/blue coloring and the inserted amplified 16SrRNA genes sequenced. Plasmid purification and sequencing reactions were performed by Macrogen Europe Inc. (Amsterdam, The Netherlands). The nucleotide sequences obtained were deposited in Gen- Bank/DDBJ/EMBL databases under accession numbers from HQ834968 to HQ835246.
For qPCR data, 1-way ANOVA with Tukey post hoc test was employed. Analyse-it 2.0 software (Analyse-It, Ldt., Leeds, UK) was used for both tests. For T-RFLP, chromatogram files from automated sequencer sizing were imported into GeneMarker ver. 1.71 software (SoftGenetics LLC, State College, PA, USA) by filtering with the default options of the module for AFLP analysis. Peaks above 100 fluorescence units and whose size ranged from 35 to 500 nt were considered for profile analysis. Only the presence/absence of peaks was considered as informative data from the chromatograms. Statistical analyses were performed on a binary matrix obtained as previously reported . Past 2.02  software package was used to compute Non-Metric Multidimensional Scaling (N-MDS). To test the distribution of the variance of T-RFLP profiles within plant tissues and among pots, Analysis of Molecular Variance (AMOVA ) was applied using Arlequin 220.127.116.11 software (http://cmpg.unibe.ch/software/arlequin3/). Although developed for population genetic analysis, the general procedure implemented by AMOVA is flexible enough to estimate the statistical significance of groups of bacterial communities as reported previously [13, 42, 59]. Pairwise F ST distances  between T-RFLP profiles of plant tissues and soils were used to infer a Neighbor-Joining dendrogram with the MEGA4 software .
Partial 16 S rRNA sequences were manually inspected for quality, then aligned and clustered with the furthest neighboring algorithm with the module present in Mothur v.1.12.3 . Diversity indices (Shannon H’ and Chao-1) were calculated with the same software. Library coverage was estimated with the formula C = 1-(n/N), where n is the number of singletons (defined at 97% sequence identity in Mothur) that are encountered only once in the library and N is the total number of sequenced clones. Taxonomic assignment was performed with the Classifier module present in Ribosomal Database Project 10 website  (http://rdp.cme.msu.edu/) at 80% confidence threshold. Sequences with 97% similarity were treated as a single Operational Taxonomic Unit (OTU). Sequences (one for each OTU) were aligned with the 16 S rRNA gene sequences of the closest match retrieved from NCBI databases, using MUSCLE  and a Neighbor-Joining dendrogram was constructed using MEGA4 . Phylogenetic inference and evolutionary distance calculations were generated using the Maximum Composite Likelihood; 1000 bootstrap replicates were used to obtain confidence estimates for the phylogenetic trees.
This work was partially supported by intramural funding of the University of Florence to AM and MB and by Soil-Sink (FISR-MIUR) project funding to MB. We are grateful to Mary Forrest, Professor of Scientific English, Department of Pharmacology, Faculty of Medicine and Surgery, University of Florence for editing English language. We acknowledge two anonymous reviewers for helpful suggestions for the improvement of the manuscript.