Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Among-Population Variation in Microbial Community Structure in the Floral Nectar of the Bee-Pollinated Forest Herb Pulmonaria officinalis L

  • Hans Jacquemyn ,

    hans.jacquemyn@bio.kuleuven.be

    Affiliation Division of Plant Ecology and Systematics, Biology Department, KU Leuven, Heverlee, Belgium

  • Marijke Lenaerts,

    Affiliations Laboratory for Process Microbial Ecology and Bioinspirational Management, Thomas Moore University College, Campus De Nayer, Department of Microbial and Molecular Systems (M2S), KU Leuven Association, Sint-Katelijne-Waver, Belgium, Scientia Terrae Research Institute, Sint-Katelijne-Waver, Belgium

  • Rein Brys,

    Affiliation Division of Plant Ecology and Systematics, Biology Department, KU Leuven, Heverlee, Belgium

  • Kris Willems,

    Affiliations Laboratory for Process Microbial Ecology and Bioinspirational Management, Thomas Moore University College, Campus De Nayer, Department of Microbial and Molecular Systems (M2S), KU Leuven Association, Sint-Katelijne-Waver, Belgium, Scientia Terrae Research Institute, Sint-Katelijne-Waver, Belgium

  • Olivier Honnay,

    Affiliation Division of Plant Ecology and Systematics, Biology Department, KU Leuven, Heverlee, Belgium

  • Bart Lievens

    Affiliations Laboratory for Process Microbial Ecology and Bioinspirational Management, Thomas Moore University College, Campus De Nayer, Department of Microbial and Molecular Systems (M2S), KU Leuven Association, Sint-Katelijne-Waver, Belgium, Scientia Terrae Research Institute, Sint-Katelijne-Waver, Belgium

Abstract

Background

Microbial communities in floral nectar have been shown to be characterized by low levels of species diversity, yet little is known about among-plant population variation in microbial community composition.

Methodology/Principal Findings

We investigated the microbial community structure (yeasts and bacteria) in floral nectar of ten fragmented populations of the bee-pollinated forest herb Pulmonaria officinalis. We also explored possible relationships between plant population size and microbial diversity in nectar, and related microbial community composition to the distance separating plant populations. Culturable bacteria and yeasts occurring in the floral nectar of a total of 100 plant individuals were isolated and identified by partially sequencing the 16S rRNA gene and D1/D2 domains of the 26S rRNA gene, respectively. A total of 9 and 11 yeast and 28 and 39 bacterial OTUs was found, taking into account a 3% (OTU0.03) and 1% sequence dissimilarity cut-off (OTU0.01). OTU richness at the plant population level (i.e. the number of OTUs per population) was low for yeasts (mean: 1.7, range: 0–4 OTUs0.01/0.03 per population), whereas on average 6.9 (range: 2–13) OTUs0.03 and 7.9 (range 2–16) OTUs0.01 per population were found for bacteria. Both for yeasts and bacteria, OTU richness was not significantly related to plant population size. Similarity in community composition among populations was low (average Jaccard index: 0.14), and did not decline with increasing distance between populations.

Conclusions/Significance

We found low similarity in microbial community structure among populations, suggesting that the assembly of nectar microbiota is to a large extent context-dependent. Although the precise factors that affect variation in microbial community structure in floral nectar require further study, our results indicate that both local and regional processes may contribute to among-population variation in microbial community structure in nectar.

Introduction

Floral nectar is a sweet, aqueous secretion containing sugars and amino acids that is offered by flowering plants to attract pollinators, mostly insects [1], [2]. Traditionally it has been assumed that nectar properties represent intrinsic plant features that are stable in time. However, recent studies have indicated that nectar is often contaminated with microorganisms, most often yeasts and bacteria, which may change the chemistry and attractiveness of nectar, potentially affecting pollination success and plant fitness [3], [4]. Although it has already been known since the early 1920’s that yeasts are common inhabitants of floral nectars [5], [6], only recently the microbial community structure in nectar and its ecological impact have been explored in more detail [7][12]. These studies have highlighted that the floral nectar of animal-pollinated plants often harbors highly specialized yeast communities. In most cases, species richness tends to be low and often only a few species are found within the nectar of a single nectary [8], [12], [13], suggesting that important filtering mechanisms (e.g. priority effects and nectar chemistry) determine community composition of nectar-inhabiting microorganisms in floral nectar [11], [14], [15].

Bacteria, on the other hand, have been less frequently studied, and there are only a few studies that have characterized bacterial communities in floral nectar [16], [17]. A recent study, investigating phylogenetic diversity of culturable bacteria in 27 South African plant species, revealed that bacteria are common in floral nectar, but that their phylogenetic diversity is rather restricted, with most isolates belonging to three major bacterial phyla, including Actinobacteria, Firmicutes, and Proteobacteria (Alphaproteobacteria, Betaproteobacteria and Gammaproteobacteria) [18]. Similar to yeast communities, species richness was also found to be low (18 operational taxonomic units (OTUs) at a 16S rRNA gene sequence dissimilarity cut-off of 3%) (but see [19]).

Given the mounting number of studies that have investigated microbial community composition in floral nectar [7], [9], [10], [12], [18][22], surprisingly little is known about among-plant population variation in community composition of these nectar-inhabiting microorganisms. Although precipitation and microorganisms in the air can be considered as constant sources of microorganisms in flowers, yeasts and bacteria are most likely transported to flowers by pollinating insects or small birds [9], [14], [15]. It is therefore reasonable to assume that limitations in their dispersal capacity may lead to significant spatial turnover of microbial community composition in floral nectar, especially when plant populations occur in highly fragmented habitats, surrounded by an inhospitable urban or agricultural landscape matrix. Recent research investigating community organization of nectar-inhabiting microorganisms in the hummingbird-pollinated shrub Mimulus aurantiacus has indeed shown significant turnover of microbial community composition, even at a very small scale [14]. However, there is currently no information available regarding differences in microbial community structure among plant populations that occur in discrete habitat fragments within a hostile matrix.

In this paper, we investigated the community structure of nectar-inhabiting microorganisms in ten fragmented populations of the bee-pollinated understory forest herb Pulmonaria officinalis (Common lungwort) in northern Belgium. Previous genetic marker based research on this plant species in the same area has shown strong genetic differentiation and significant isolation-by-distance [23]. Because gene flow in this species occurs mainly through pollen [24], these results indicate that pollen dispersal is mainly restricted to neighboring populations. Assuming that pollinators are the main dispersal agents of nectar-inhabiting microorganisms [7] and given that local populations occur in forest fragments with pronounced differences in local environmental conditions, it can be expected that community composition of nectar-inhabiting microorganisms is more dependent on population characteristics, such as nectar quality, population size or local plant community composition, than on geographic isolation. As a result, similarity in community composition between populations is expected to be low. To test these general predictions, the presence of culturable yeasts and bacteria was assessed for each population and microbial species richness at the plant population level was related to the size of the plant population. Finally, similarity in microbial community composition between populations was assessed and related to the distance between plant populations.

Materials and Methods

Ethics Statement

All necessary permits were obtained for the described field studies.

Study Species

Pulmonaria officinalis L. is a perennial forest herb that grows in species-rich mixed and open forests, characterized by relatively humid, wet and loamy soils. Its distribution range is located in Mid-East Europe, but fragmented populations reach till Britain and Denmark. The species is wintergreen and flowers early in the growth season, from March until the end of April. Flowers exhibit reciprocal herkogamy and several ancillary polymorphisms [25]. During anthesis, the colour of the corolla gradually changes from red through purple to violet and finally blue, offering a visual sign to pollinators which flowers are most rich in nectar [26]. Nectar is secreted at the bottom of the corolla tube, where it accumulates. Within the study area (northern Belgium, Flanders, Fig. 1), flowers are visited by generalist insect species, including Bombus terrestris, B. pascuorum, B. pratorum and Bombylius major, but only the long-tongued Anthophora plumipes was shown to serve as an efficient pollinator [25].

thumbnail
Figure 1. Geographic distribution of ten populations of the bee-pollinated forest herb Pulmonaria officinalis used in this study (indicated in red) among several other P. pulmonaria populations in the same study area (northern Belgium).

https://doi.org/10.1371/journal.pone.0056917.g001

Study Region and Nectar Sampling

This study was conducted in the region around Brakel, where several P. officinalis populations are known to occur in small, highly fragmented forest patches [23], [27]. Within this region, ten populations of P. officinalis were randomly selected (Fig. 1). The distances separating the studied populations varied between 0.74 and 16.73 km (mean: 7.30 km). All populations were located in small forest fragments embedded within a hostile agricultural and urban landscape matrix. At the beginning of April 2012, the number of flowering individuals was counted for each population. Simultaneously, in each population ten individuals were randomly selected for nectar sampling, and from each individual five flowers were harvested, yielding a total of 50 flowers per population.

Within 24 h after harvesting, nectar was extracted using sterile 5-µl microcapillaries. Nectar of flowers from the same plant was pooled and diluted in 150 µl of sterile distilled H2O [18], yielding a total of 100 nectar samples. Since floral nectar usually contains high concentrations of sucrose and other sugars and can also contain high levels of inorganic ions, nectar dilutions (even in distilled H2O) are not hypotonic and both bacteria and yeasts have been shown to remain viable in nectar dilutions in distilled H2O for several months [18]. Subsequently, fifty microliters was plated on both trypticase soy agar (TSA; Oxoid) and yeast extract peptone dextrose agar (YPDA; Difco), representing a general growth medium for bacteria and yeasts, respectively. These media have been used successfully for isolating microorganisms from nectar previously [10][12], [18]. Plates were incubated at 25°C for 10 days. For each plate one colony was picked for each morphologically distinct colony type, and further subcultivated to obtain pure cultures [18]. In addition, a preliminary screen of several morphologically identical colonies from the same plate had revealed that they all belonged to the same species, illustrating the suitability of the used approach. The obtained bacterial and yeast isolates were stored at −80°C in trypticase soy broth (Oxoid) and yeast extract peptone dextrose broth (Difco) containing 37.5% glycerol, respectively.

DNA Extraction, PCR Amplification and Sequencing

For each isolate, genomic DNA was extracted from five-day old cultures grown on either TSA (bacteria) or YPDA (yeasts) using the phenol–chloroform extraction method [28]. Samples were amplified in a reaction volume of 20 µl, containing 312.5 µM of each dNTP, 1.0 µM of each primer, 1.25 units TaKaRa ExTaq polymerase, 1× Ex Taq Buffer (Clontech Laboratories, Palo Alto, CA, USA), and 5 ng genomic DNA (as measured by a Nanodrop spectrophotometer). DNA amplification of the D1/D2 domain of the large subunit rRNA and 16S rRNA gene was performed using the primer sets NL1–NL4 [29]) and 27F-1492R [30] for yeasts and bacteria, respectively. When amplification failed using the latter pair, primers 1387R [31] or 1541R [32] were used as reverse primer. Before amplification, DNA samples were denatured at 94°C for 2 min. Next, 30 cycles were run consisting of 45 s at 94°C, 45 s at 55°C (for NL1–NL4) or 59°C (for 27F-1492R/1387R/1541R), and 45 s at 72°C, with a final extension at 72°C for 10 min. Sequencing was performed using the reverse primer used for DNA amplification.

Data Analysis

The obtained sequences were compared with reference sequences using BLAST software [33] and the Ribosomal Database Project (RDP) website [34] (http://rdp.cme.msu.edu/). Isolates were assigned to the highest taxonomic rank possible (generally the species level) using the RDP classifier, BLAST analysis (uncultured/environmental sample sequences excluded), and based on the nearest neighbors in a phylogenetic tree containing GenBank sequences from several type strains showing the highest sequence homology to our sequences. To this end, both our sequences and the reference sequences were aligned with Clustal W implemented in MEGA5 [35], followed by trimming to consensus start and end motifs. Phylogenetic trees were calculated by the neighbour-joining method [36] implemented in Clustal X and displayed by TreeView version 1.6.6 [37]. Support of internal nodes was assessed using bootstrap analysis performed with 1000 replications. For ease of visualization, highly similar sequences (>99% sequence identity) were restricted to one representative sequence per OTU. In all cases, presumptive identifications based on top BLAST hits were confirmed by the nearest neighbor in the phylogenetic tree containing type strain sequences.

For subsequent analyses, bacterial and yeast OTUs were assigned using the Mothur v.1.23.1 software program [38]. DNA dissimilarity cut-offs of 1% and 3% were used in these analyses. For each OTU, the capability to grow in nectar was verified for a few isolates according [7], [18], [39]. All isolates tested were found to tolerate sucrose concentrations of at least 50% (w/v). In addition, all examined bacterial isolates showed catalase activity, suggesting that the detected OTUs are physiologically capable to overcome the presence of toxic hydrogen peroxide in nectar. Representative sequences for each OTU were deposited in GenBank (accession numbers KC433478–KC433527).

In order to assess the overall richness of microbial OTUs in the whole study region, sample-based rarefaction methods were applied to species presence-absence [40], [41]. Since the nectar of multiple flowers from a single plant was combined, individual plants rather than nectar drops were considered as sample units [12]. In this analysis, OTU occurrence data from all individuals were analyzed together, irrespective of the population of origin, yielding a rarefaction curve that assesses overall species richness of nectar yeasts and bacteria at the landscape scale. Rarefaction curves were computed using EstimateS version 8.2 [40], with 50 randomizations and sampling without replacement. Analyses were performed for bacteria and yeasts separately. Additionally, as our taxa richness data are based on incidence, the expected yeast and bacterial OTU richness in nectar was also determined using the nonparametric estimator Chao2 [42]. Rarefaction generates the expected number of species (OTUs) in a small collection of n samples drawn at random from the large pool of N samples [43]. In contrast, richness estimators predict the total richness of a community from samples [40].

For each population, OTU richness of bacteria and yeasts was determined by counting the total number of different bacterial and yeast OTUs. The observed richness was related to the size of the population using the Pearson product-moment correlation coefficient. Population size was log transformed prior to analysis. To visualize differences in microbial community structure among populations, we applied non-metric multidimensional scaling (NMDS) ordination techniques using the program PC-ORD version 6 [44]. As distance measure, we used the Bray-Curtis coefficient. This coefficient is also known as the Sørensen or Czekanowski coefficient and is considered as one of the most robust measures for this purpose [45]. In addition, pairwise similarity matrices were created to determine microbial community similarity between populations. The Jaccard index was used to describe the similarity in composition of the bacterial and yeast nectar communities [46]. For each population, the nearest distance (bird’s eye view) to any other population was also determined. A Mantel test was used to test the hypothesis that community similarity was related to the distance separating populations. Statistical significance was determined using 9999 randomizations in PopTools [47].

Results

Bacterial and yeast isolates were obtained from both TSA and YPDA. Following isolation and purification, a total of 37 yeast and 152 bacterial isolates was obtained from nectar samples from 24 and 59 P. officinalis plants respectively, with 18 plants containing both yeasts and bacteria in their nectar (Fig. 2). For 35 out of the 100 nectar samples, no microbial growth was observed. Yeasts were recovered from nine out of ten populations, whereas bacteria were found in all sampled populations (Fig. 2).

thumbnail
Figure 2. Distribution of culturable nectar yeasts and bacteria in floral nectar of ten plant individuals from ten Pulmonaria officinalis populations.

https://doi.org/10.1371/journal.pone.0056917.g002

Using a 3% sequence dissimilarity cut-off value, nine yeast OTUs (OTUs0.03) were identified, comprising both ascomycetous and basidiomycetous yeasts (Table 1). Two additional OTUs were identified when the dissimilarity cut-off was lowered to 1% (OTUs0.01), corresponding to 11 different yeast species (Table 1). These included for example Metschnikowia reukaufii, Candida bombi, Sporobolomyces roseus and several Cryptococcus species (Fig. 3). Rarefaction curves showed that the number of OTUs was relatively close to saturation (Fig. 4a). However, the nonparametric richness estimator Chao2 gradually shifted from the observed species richness, indicating that our sampling only detected a part of the total estimated yeast species richness. Most likely, the erratic behavior of the Chao2 estimator was caused by the overall low yeast abundance. M. reukaufii and C. bombi were recorded as the most common yeast species, occurring in five and three populations, respectively. All other yeast OTUs were only observed in a single population.

thumbnail
Figure 3. Neighbour-joining phylogram showing phylogenetic relationships between different large subunit rRNA gene sequences from Pulmonaria officinalis nectar-inhabiting yeasts and reference sequences of the most related type strains found in GenBank.

Bootstrap percentages based on 1000 replications are shown at the major nodes.

https://doi.org/10.1371/journal.pone.0056917.g003

thumbnail
Figure 4. Rarefaction curves (bold, solid line) and the nonparametric estimator Chao2 (thin solid line) of microbial OTU richness found in the floral nectar of 100 sampled individuals of Pulmonaria officinalis.

Dotted lines represent 95% confidence intervals. Rarefaction curves are given for a) yeast and b) bacterial OTUs (based on a DNA dissimilarity cut-off value of 1%).

https://doi.org/10.1371/journal.pone.0056917.g004

thumbnail
Table 1. Yeast operational taxonomic units (OTUs) identified in this study.

https://doi.org/10.1371/journal.pone.0056917.t001

Using a 3% cut-off value, a total of 28 bacterial OTUs (OTUs0.03) was detected (Table 2). By lowering the dissimilarity cut-off to 1% 11 additional OTUs were found, resulting in a total of 39 OTUs0.01 (Table 2). As for the yeasts, the rarefaction curves were relatively close to reach a plateau. However, in contrast with the yeasts, the Chao 2 estimator gave a predicted OTU richness which was close to the number of observed OTUs, resulting in an estimated richness of the nectar bacterial community of 39 OTU0.03 and 54 OTU0.01 (Fig. 4b). The recovered bacteria belonged to three major phyla, including Actinobacteria (18 OTUs0.03; 25 OTUs0.01), Firmicutes (4 OTUs0.03; 7 OTUs0.01) and Proteobacteria (Alpha and Gamma subdivisions; 6 OTUs0.03; 7 OTUs0.01) (Table 2). The most common bacteria were species from the genera Rhodococcus, Microbacterium and Methylobacterium which were retrieved in five or more populations (Table 2). Other OTUs that were identified (>97.5% sequence homology with GenBank sequence) included members from the genera Arthrobacter, Bacillus, Brachybacterium, Brevibacterium, Devosia, Erwinia, Enhydrobacter, Flexivirga, Gordonia, Janibacter, Luteipulveratus, Micrococcus, Moraxella, Nocardioides, Okibacterium, Plantibacter, Ponticoccus, Pseudomonas, Rhodanobacter, Saxeibacter, Staphylococcus and Streptomyces (Table 2; Fig. 5). Although the presence of Micrococcus and Staphylococcus may suggest possible contamination as these bacteria may also occur on the skin of humans and animals, we clearly showed that the detected species were able to resist high sugar concentrations typically experienced in nectar. In addition, members of Staphylococcus and the Micrococcaceae family have been isolated from other nectar sources as well [18], [19]. Therefore, we can reasonably assume that all bacteria obtained in our study can be considered as true nectar-inhabiting microbes.

thumbnail
Figure 5. Neighbour-joining phylogram showing phylogenetic relationships between 16S rRNA gene sequences from Pulmonaria officinalis nectar-inhabiting bacteria and reference sequences of the most related type strains found in GenBank.

For ease of visualization, the dataset was limited to one representative sequence (see Table 1 and 2) for each operational taxonomic unit (OTU) found in this study at a DNA dissimilarity cut-off value of 1%. Bootstrap percentages based on 1000 replications are shown at the major nodes.

https://doi.org/10.1371/journal.pone.0056917.g005

thumbnail
Table 2. Bacterial operational taxonomic units (OTUs) identified in this study.

https://doi.org/10.1371/journal.pone.0056917.t002

The sampled P. officinalis populations differed in size between 98 and >5000 flowering individuals. The number of yeast OTUs observed per population varied between 0 and 4 (mean: 1.7) (irrespective of the cut-off value used), and was not significantly related to population size (r = −0.18, P>0.05) (Fig. 5). The number of bacterial OTUs observed per population varied between 2 and 13 OTUs0.03 (mean: 6.9) and 2 and 16 OTUs0.01 (mean: 7.9), and was also not significantly related to the size of the plant population (r = −0.03, P>0.05) (Fig. 6). Although the NMDS analysis showed that some geographic clustering in microbial community (taking into account both bacteria and yeasts) was present (Fig. 7), the overall similarity in community composition was low (average Jaccard index: 0.14). Populations 7, 8 and 9 and populations 2, 4 and 5 formed distinct clusters on the NMDS graph which to some extent coincided with their geographic location in the landscape (Fig. 1). However, there was no significant relationship between community similarity and geographic distance (rM = −0.22, P = 0.09).

thumbnail
Figure 6. Relationship between population size (number of flowering ramets) and the number of bacterial and yeast OTUs (based on a DNA dissimilarity cut-off value of 1%) in the floral nectar of Pulmonaria officinalis.

https://doi.org/10.1371/journal.pone.0056917.g006

thumbnail
Figure 7. NMDS ordination of the total microbial community composition (bacteria and yeasts) in the floral nectar of Pulmonaria officinalis obtained from sampling flowers from ten individuals in ten populations.

Numbers refer to the populations depicted in Fig. 1.

https://doi.org/10.1371/journal.pone.0056917.g007

Discussion

Microorganisms

In this study a wide variety of nectar-inhabiting microorganisms was found in the floral nectar of the early-flowering forest herb P. officinalis. Using a 1% dissimilarity cut-off, a total of eleven different yeast species was identified, many of which have been recorded in nectar before, including for example M. reukaufii, C. bombi, C. victoriae, C. macerans, and S. roseus [8], [12]. Although the number of species retrieved seems not to be unusual compared to results reported in similar studies [8], [12], it is surprising to see that within a single plant species almost as much yeast species were found as were reported in previously published datasets, which mostly covered a large number of plant species. For example, studying 128 nectar drops from 24 plant species in Spain yielded 216 yeast isolates and a total of 12 different yeast species [12]. Similarly, 11 different yeast species in 42 isolates were found in a total of 11 tropical plant species from the Yucatan Peninsula, Mexico [8].

There is less information about the occurrence and distribution of bacterial species within the floral nectar of individual plants and species. Recently, 53 bacterial isolates were recovered from 38 nectar samples (53.5% of all investigated nectar drops) from 27 plant species belonging to 13 plant families occurring in South Africa [18]. This yielded a total of 18 and 24 bacterial OTUs at the 3% and 1% 16S rRNA gene dissimilarity cut-off, respectively. In P. officinalis, in 59% of the sampled plants bacterial OTUs were found. However, the number of bacterial OTUs was about twice as large as that found in [18], i.e. 28 and 39 at 3% and 1% dissimilarity cut-off, respectively. Similar to the results of [18], all OTUs recovered belonged to only three bacterial phyla (Actinobacteria, and to a lesser extent to Firmicutes and Proteobacteria), confirming previous findings that microbial communities in nectar are characterized by low phylogenetic diversity. Using both catalase activity and sucrose tolerance tests we also showed that all recovered species can hydrolyze hydrogen peroxide and display high osmotolerance, and thus are able to overcome some of the main stressors found in floral nectar. Interestingly, whereas Proteobacteria were the dominant bacteria in South-African plants, about 50% of the retrieved OTUs belonged to Actinobacteria in P. officinalis. The most prevalent genera of this phylum found in this study included Microbacterium, Rhodococcus and Streptomyces. Species belonging to these genera are known to thrive in a broad range of environments, including soils and many plant-associated environments such as roots and leaves. However, as far as we know, these bacteria have not been associated with nectar so far.

Species Richness

In 35 individuals no culturable microorganisms were found, which may suggest that these plants were either not visited by insects [7], that no transfer of microorganisms that can survive and develop in nectar had occurred during flower visits [15], or that the present microbiota represented only non-culturable microorganisms. The latter, for example, could be assessed using culture-independent methods such as 454 pyrosequencing. These methods generally have a higher resolution compared to culture-based approaches, which are restricted to the isolation of culturable microorganisms. So far, however, no comparison has been made between conventional methods and culture-independent methods for studying nectar microbiota. A major drawback of culture-independent methods, on the other hand, is the fact that no isolates are available to investigate or confirm specific features that allow these microorganisms to survive and grow in nectar [7], [18], [38]. We further found that most plants contained only nectar bacteria (41%), while fewer plants were found containing both bacteria and yeast (18%) and only a minority containing only yeasts (6%). These results indicate that at least in this plant species bacteria may be much more widespread in nectar than yeast species and contrast with findings of Álvarez-Pérez and Herrera [48], who recently showed that bacteria and yeasts generally coexisted in floral nectar of a selection of wild Mediterranean plants.

The number of yeast species per plant population was low (on average 1.7 OTUs per population), indicating that only a few yeast species dominated in a population. Low diversity of yeast species is in line with previous studies investigating yeast diversity in individual floral nectar samples [8], [10], [12]. For example, in two populations of the winter-blooming herb Helleborus foetidus the nectar was dominated by a single yeast species (M. reukauffi), although several different yeast species were observed on the bodies of visiting insects [15]. The dominance of a particular yeast species in the floral nectar of plant species has been explained by filtering mechanisms, such as priority effects, which predict that early-arriving species have a competitive advantage toward late-arriving species. Using laboratory experiments, priority effects appeared to be important in structuring microbial communities in floral nectar [11]. However, results depended strongly on the phylogenetic relationships of the yeast strains involved. Priority effects were particularly strong between closely related species, whereas effects were less pronounced for phylogenetically distantly related species. Overall, these results suggest that yeast species can outcompete other species, and that the first species to colonize and spread within a plant population can become the dominant yeast in the population. Nevertheless, in contrast to Herrera et al. [15], who found yeasts in 72.5% of the investigated H. foetidus nectar samples, in our study yeasts were found in less than 25% of the P. officinalis individuals tested, suggesting that other factors such as chemical nectar composition may also play an important role in the distribution of nectar yeasts in P. officinalis. Clearly, more research is needed to investigate whether priority effects really are the dominant factor determining yeast community organization in this species, or whether the occurrence of yeast species is affected by the presence of bacteria, and vice versa [48].

Community Turnover

We found low community similarity and no significant turnover in microbial community composition among populations. These results may either suggest that little exchange of microorganisms between populations occurred or that nectar conditions differed between populations thereby selecting for different microbial communities. These results are in line with genetic work that has shown a strong genetic differentiation between P. officinalis populations [23]. Alternatively, because the local community of co-flowering plant species also differed substantially between populations, this may additionally have affected the species composition of microorganisms potentially colonizing P. officinalis flowers. Since nectar is considered to be initially sterile [7], pollinators are believed to be the main vectors transferring micro-organisms from one plant to the other, and between populations [14], [15].

Given that the most common pollinators of P. officinalis in the study region are generalist pollinators [27] that visit several co-flowering species, it is reasonable to expect that microbial community composition in the floral nectar of P. officinalis populations reflects to some extent local plant species composition. Although it is likely that pollinators (most often bumblebees and bees) can cross smaller distances across agricultural landscapes, it is unlikely they fly across very large distances further contributing to the low similarity in species composition. Recent findings of Belisle et al. [14] have shown that non-random small-scale foraging of pollinators resulted in non-random distributions of nectar-inhabiting yeasts in the hummingbird-pollinated Mimulus aurantiacus, but it is unlikely that this foraging behavior contributes substantially to large-scale patterns of community turnover. The observed strong genetic differentiation between P. officinalis populations and significant isolation-by-distance [23] support this hypothesis.

Conclusion

Nectar of the bee-pollinated forest herb P. officinalis was commonly colonized by microorganisms, both bacteria and to a lesser extent yeasts. However, large variation in community structure was observed between populations. The inability of pollinators to cross larger distances across hostile agricultural and urban landscapes has probably contributed to the observed low similarity in community composition. However, the importance of variation in nectar properties between populations or differences in local species composition of co-flowering plants as drivers of microbial community composition cannot be ruled out. Especially in plant species that are pollinated by generalist pollinators that visit several co-flowering species at the same time, local microbial community structure cannot be studied independently from the local plant community. These findings thus suggest that the assembly of the nectar microbiota is context-dependent. More research, both experimental and observational studies, is therefore needed to elucidate the ecological mechanisms explaining variation in microbial community structure within and among populations and to disentangle the importance of local and regional factors.

Acknowledgments

We are grateful to Ingeborg Frans for her help with the nectar isolations and Pieter Busschaert and Sam Crauwels for their help with the sequence analyses. Sergio Álvarez-Pérez provided very useful comments that improved the quality of this manuscript.

Author Contributions

Conceived and designed the experiments: HJ RB BL. Performed the experiments: ML. Analyzed the data: HJ ML BL. Contributed reagents/materials/analysis tools: HJ BL KW. Wrote the paper: HJ RB OH BL.

References

  1. 1. Heil M (2011) Nectar: generation, regulation and ecological functions. Trends in Plant Science 16: 191–200.
  2. 2. Nicolson SW, Thornburg RW (2007) Nectar chemistry. In: Nicolson SW, Nepi M, Pacini E, editors. Nectaries and Nectar. Dordrecht: Springer-Verlag. 215–263.
  3. 3. Herrera CM, Pozo MI, Medrano M (2013) Yeasts in nectar of an early-blooming herb: sought by bumble bees, detrimental to plant fecundity. Ecology, in press. doi:10.1890/12-0595.1.
  4. 4. Vannette RL, Gauthier MPL, Fukami T (2013) Nectar bacteria, but not yeast, weaken a plant-pollinator mutualism. Proceedings of the Royal Society of London Series B - Biological Sciences 280: 20122601.
  5. 5. Guilliermond A (1920) The yeasts. London: Chapman and Hall.
  6. 6. Jimbo T (1926) Yeasts isolated from flower nectar. Science Reports of the Tohoku Imperial University 2: 161–187.
  7. 7. Brysch-Herzberg M (2004) Ecology of yeasts in plant-bumblebee mutualism in Central Europe. FEMS Microbiology Ecology 50: 87–100.
  8. 8. Canto A, Herrera CM (2012) Micro-organisms behind the pollination scenes: microbial imprint on floral sugar variation in a tropical plant community. Annals of Botany 110: 1173–1183.
  9. 9. de Vega C, Herrera CM, Johnson SD (2009) Yeasts in floral nectar of some South African plants: quantification and associations with pollinator type. South African Journal of Botany 75: 798–806.
  10. 10. Herrera CM, de Vega C, Canto A, Pozo MI (2009) Yeasts in floral nectar: a quantitative survey. Annals of Botany 103: 1415–1423.
  11. 11. Peay KG, Belisle M, Fukami T (2012) Phylogenetic relatedness predicts priority effects in nectar yeast communities. Proceedings of the Royal Society of London Series B - Biological Sciences 279: 749–758.
  12. 12. Pozo MI, Herrera CM, Bazaga P (2011) Species richness of yeast communities in floral nectar of southern Spanish plants. Microbial Ecology 61: 82–91.
  13. 13. Pozo MI, Lachance MA, Herrera CM (2012) Nectar yeasts of two southern Spanish plants: the roles of immigration and physiological traits in community assembly. FEMS Microbiology Ecology 80: 281–293.
  14. 14. Belisle M, Peay KG, Fukami T (2012) Flowers as islands: spatial distribution of nectar-inhabiting microfungi among plants of Mimulus aurantiacus: a hummingbird-pollinated shrub. Microbial Ecology 63: 711–718.
  15. 15. Herrera CM, Canto A, Pozo MI, Bazaga P (2010) Inhospitable sweetness: nectar filtering of pollinator-borne inocula leads to impoverished, phylogenetically clustered yeast communities. Proceedings of the Royal Society of London Series B - Biological Sciences 277: 747–754.
  16. 16. Ehlers BK, Olesen JM (1997) The fruit-wasp route to toxic nectar in Epipactis orchids? Flora 192: 223–229.
  17. 17. Gilliam M, Moffett JO, Kauffeld NK (1983) Examination of floral nectar of Citrus, cotton, and Arizona dessert plants for microbes. Apidologie 14: 299–302.
  18. 18. Álvarez-Pérez S, Herrera CM, de Vega C (2012) Zooming-in on floral nectar: a first exploration of nectar-associated bacteria in wild plant communities. FEMS Microbiology Ecology 80: 591–602.
  19. 19. Fridman S, Izhaki I, Gerchman Y, Halpern M (2012) Bacterial communities in floral nectar. Environmental Microbiology Reports 4: 97–104.
  20. 20. Lachance MA, Starmer WT, Rosa CA, Bowles JM, Barker JSF, et al. (2001) Biogeography of the yeasts of ephemeral flowers and their insects. FEMS Yeast Research 1: 1–8.
  21. 21. Mushtaq M, Ayesha J, Sharfun N (2006) Biodiversity of yeast mycoflora in nectar of Bombax ceiba and Canna indica flower. Pakistan Journal of Botany 38: 1279–1288.
  22. 22. Mushtaq M, Jamal A, Nahar S (2007) Biodiversity of yeast mycoflora in nectar of Hibiscus rosa-sinensis and Ixora coccinea flowers. Pakistan Journal of Botany 39: 1367–1376.
  23. 23. Meeus S, Honnay O, Brys R, Jacquemyn H (2012) Biased morph ratios and skewed mating success contribute to loss of genetic diversity in the distylous Pulmonaria officinalis. Annals of Botany 109: 227–235.
  24. 24. Brys R, Jacquemyn H (2010) Impact of floral display and spatial distribution of potential mates affect pollen deposition and female reproductive success in distylous Pulmonaria officinalis (Boraginaceae). Plant Biology 12: 597–603.
  25. 25. Brys R, Jacquemyn H, Hermy M, Beeckman T (2008) Pollen deposition rates and the functioning of distyly in the perennial Pulmonaria officinalis (Boraginaceae). Plant Systematics and Evolution 273: 1–12.
  26. 26. Oberrath R, Zanke C, Böhning-Gaese K (1995) Triggering and ecological significance of floral color change in Lungwort (Pulmonaria spec.). Flora 190: 155–159.
  27. 27. Brys R, Jacquemyn H, Beeckman T (2008) Morph ratio variation, population size and female reproductive success in the distylous Pulmonaria officinalis (Boraginaceae). Journal of Evolutionary Biology 21: 1281–1289.
  28. 28. Lievens B, Brouwer M, Vanachter ACRC, Lévesque CA, Cammue BPA, et al. (2003) Design and development of a DNA array for rapid detection and identification of multiple tomato vascular wilt pathogens. FEMS Microbiology Letters 223: 113–122.
  29. 29. O’Donnell K (1993) Fusarium and its near relatives. In: Reynolds DR, Taylor JW, editors. The Fungal Holomorph: Mitotic, Meiotic and Pleomorphic Speciation in Fungal Systematics. Wallingford: CAB International. pp.225–233.
  30. 30. Lane DJ (1991) 16S/23S rRNA sequencing. In: Stackebrandt E, Goodfellow M, editors. Nucleic Acid Techniques in Bacterial Systematics. New York: John Wiley and Sons. pp.115–175.
  31. 31. Marchesi JR, Sato T, Weightman AJ, Martin TA, Fry JC, et al. (1998) Design and evaluation of useful bacterium-specific PCR primers that amplify genes coding for bacterial 16S rRNA. Applied and Environmental Microbiology 64: 795–799.
  32. 32. Liu H, Xu Y, Ma Y, Zhou P (2000) Characterization of Micrococcus antarcticus sp. nov., a psychrophilic bacterium from Antarctica. International Journal of Systematic and Evolutionary Microbiology 50: 715–719.
  33. 33. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. Journal of Molecular Biology 215: 403–410.
  34. 34. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, et al. (2009) The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Research 37: D141–D145.
  35. 35. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution 28: 2731–2739.
  36. 36. Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution 4: 406–425.
  37. 37. Page RDM (1996) TreeView: an application to display phylogenetic trees on personal computers. Computer Applications in the Biosciences 12: 357–358.
  38. 38. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, et al. (2009) Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology 75: 37–41.
  39. 39. Álvarez-Pérez S, Lievens B, Jacquemyn H, Herrera CM (2012) Acinetobacter nectaris sp. nov. and Acinetobacter boissieri sp. nov., two novel bacterial species isolated from floral nectar of wild Mediterranean insect-pollinated plants. International Journal of Systematic and Evolutionary Microbiology, doi: 10.1099/ijs.0.043489-0.
  40. 40. Colwell RK (2009) EstimateS: statistical estimation of species richness and shared species from samples. Version 8.2. User’s guide and application published at: http://purl.oclc.org/estimates.
  41. 41. Gotelli NJ, Colwell RK (2001) Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecology Letters 4: 379–391.
  42. 42. Chao A, Chazdon RL, Colwell RK, Shen TJ (2005) A new statistical approach for assessing similarity of species composition with incidence and abundance data. Ecology Letters 8: 148–159.
  43. 43. Simberloff D (1978) Use of rarefaction and related methods in ecology. In: Dickson KL, Cairns J, Linvingston RJ, editors) Biological data in water pollution assessment: quantitative and statistical analyses. Philadelphia: American Society for Testing and Materials. 150–165.
  44. 44. McCune B, Mefford MJ (2011) PC-ORD. Multivariate Analysis of Ecological Data. Version 6. MjM Software Design, Gleneden Beach, Oregon.
  45. 45. Faith DP, Minchin PR, Belbin L (1987) Compositional dissimilarity as a robust measure of ecological distance. Vegetatio 69: 7–68.
  46. 46. Anderson MJ, Crist TO, Chase JM, Vellend M, Inouye BD, et al. (2011) Navigating the multiple meanings of beta diversity: a roadmap for the practicing ecologist. Ecology Letters 14: 19–28.
  47. 47. Hood GM (2010) PopTools version 3.2.5. Available on the internet. URL http://www.poptools.org.
  48. 48. Álvarez-Pérez S, Herrera CM (2013) Composition, richness and non-random assembly of culturable bacterial–microfungal communities in floral nectar of Mediterranean plants. FEMS Microbiology Ecology, in press.
  49. 49. Wuczkowski M, Passoth V, Turchetti B, Andersson AC, Olstorpe M, et al. (2011) Description of Holtermanniella gen. nov., including Holtermanniella takashimae sp. nov. and four new combinations, and proposal of the order Holtermanniales to accommodate tremellomycetous yeasts of the Holtermannia clade. International Journal of Systematic and Evolutionary Microbiology 61: 680–689.