Impact Factor 5.640 | CiteScore 7.3
More on impact ›


Front. Microbiol., 13 July 2017 |

Eelgrass Leaf Surface Microbiomes Are Locally Variable and Highly Correlated with Epibiotic Eukaryotes

Mia M. Bengtsson1*, Anton Bühler2, Anne Brauer1, Sven Dahlke3, Hendrik Schubert2 and Irmgard Blindow3
  • 1Institute of Microbiology, University of Greifswald, Greifswald, Germany
  • 2Institut für Biowissenschaften, University of Rostock, Rostock, Germany
  • 3Biological Station of Hiddensee, University of Greifswald, Kloster, Germany

Eelgrass (Zostera marina) is a marine foundation species essential for coastal ecosystem services around the northern hemisphere. Like all macroscopic organisms, it possesses a microbiome (here defined as an associated prokaryotic community) which may play critical roles in modulating the interaction of eelgrass with its environment. For example, its leaf surface microbiome could inhibit or attract eukaryotic epibionts which may overgrow the eelgrass leading to reduced primary productivity and subsequent eelgrass meadow decline. We used amplicon sequencing of the 16S and 18S rRNA genes of prokaryotes and eukaryotes to assess the leaf surface microbiome (prokaryotes) as well as eukaryotic epibionts in- and outside lagoons on the German Baltic Sea coast. Prokaryote microbiomes varied substantially both between sites inside lagoons and between open coastal and lagoon sites. Water depth, leaf area and biofilm chlorophyll a concentration explained a large amount of variation in both prokaryotic and eukaryotic community composition. The prokaryotic microbiome and eukaryotic epibiont communities were highly correlated, and network analysis revealed disproportionate co-occurrence between a limited number of eukaryotic taxa and several bacterial taxa. This suggests that eelgrass leaf surfaces are home to a mosaic of microbiomes of several epibiotic eukaryotes, in addition to the microbiome of the eelgrass itself. Our findings thereby underline that eukaryotic diversity should be taken into account in order to explain prokaryotic microbiome assembly and dynamics in aquatic environments.


Seagrasses are aquatic flowering plants that form underwater meadows critically important for coastal ecosystems around the world. Seagrass meadows are nursing grounds for juvenile fish which can hide and forage between the seagrass leaves, sediments are stabilized by the seagrass roots and the biomass of seagrass and associated organisms sequester carbon with implications for climate change mitigation (Fourqurean et al., 2012). The seagrass eelgrass (Zostera marina) is an important foundation species along soft-bottom coasts in the northern hemisphere. Eelgrass meadows are biodiversity hotspots, providing a home to a myriad of invertebrates, algae, fish and microorganisms. Due to its ability to tolerate low and fluctuating salinity levels, eelgrass is common in estuaries as well as in the biggest brackish water habitat in the world, the Baltic Sea. However, significant declines in the depth limit and areal cover of eelgrass meadows in the Baltic Sea have been observed since several decades (Boström et al., 2014). This echoes the distressing situation for seagrass ecosystems around the world, which experience threats by human activities and global change (Orth et al., 2006).

The mechanisms that are responsible for eelgrass meadow decline appear to depend on many different factors. Historically, outbreaks of the pathogenic protist Labyrinthula zosterae (“wasting disease”) has caused catastrophic die-off of eelgrass (Muehlstein et al., 1991), although its importance as a pathogen under current conditions in Europe seems limited (Brakel et al., 2014). Instead, factors such as water clarity, eutrophication, grazing pressure, and interspecific competition have been identified as culprits for eelgrass productivity decline (Baden et al., 2010; Duffy et al., 2015). Especially important is the competition between eelgrass and algae that grow within meadows, often as epiphytes on the eelgrass itself. In some cases, eelgrass becomes extensively covered with a mixture of algae, bacteria and sessile animals, which shade it and inhibit transport of solutes (Brodersen et al., 2015) and may over time cause eelgrass meadows to degrade. Recent research suggests that relative success of eelgrass and associated algae is determined by complex interactions between biotic processes such as grazing on both algae and eelgrass by invertebrates, and environmental factors such as nutrient concentrations and temperature (Eklöf et al., 2012; Alsterberg et al., 2013). However, these mechanisms are not well understood and their complexity necessitates a holistic ecosystem approach taking into account several organism groups that inhabit seagrass meadows and their interactions to be resolved (Boström et al., 2014; Maxwell et al., 2016).

An overlooked group of organisms in eelgrass meadows is the prokaryotes (bacteria and archaea), which cover eelgrass leaves and roots forming biofilms, are responsible for degradation of eelgrass detritus and are also associated to all other organisms in eelgrass beds including algae and other epiphytes. Prokaryotes are the first colonizers on new seagrass leaves and thus initiate a successional process that may end with severe epiphytic overgrowth. They are also an important food source for microbial eukaryotic and invertebrate grazers and may therefore support a substantial part of the foodweb found within eelgrass meadows. Research on eelgrass leaf prokaryotic communities is very limited, but early work has determined that prokaryotic abundance and productivity on eelgrass leaves vary during the year, with a peak in early autumn (Törnblom and Søndergaard, 1999), and that eelgrass leaves have an associated prokaryotic community (the “microbiome”), with only some compositional overlap with other aquatic macrophytes (Crump and Koch, 2008).

There is currently little known about what factors influence eelgrass microbiomes, or what functions the prokaryotes have on eelgrass leaves. Bacteria and archaea are likely to play a fundamental role in the competition of eelgrass and epibionts such as algae. For example, certain bacteria may inhibit the attachment of algal spores while others may promote further colonization (Celdrán et al., 2012; Mieszkin et al., 2013). Conversely, epibiotic eukaryotes including algae may in turn shape the prokaryotic communities on eelgrass leaves by release of chemical attractants or deterrents (Steinberg and de Nys, 2002) or through selective grazing (Huws et al., 2005), for example. Biotic and abiotic environmental factors as well as host-related factors such as eelgrass productivity and genotype are also likely to shape the eelgrass microbiome, and thereby modulate its function.

In this study, we aimed to obtain a first view onto eelgrass leaf microbiomes by investigating the community composition of prokaryotic and eukaryotic epibiotic organisms in relation to abiotic and biotic environmental variables. We sampled eelgrass in semi-sheltered lagoons and along exposed open shorelines around the Island of Hiddensee on the eastern German Baltic Sea coast and used high-throughput Illumina amplicon sequencing of the 16S and 18S rRNA genes. We hypothesized that (1) both prokaryotic and eukaryotic communities would vary substantially between lagoon and open coast sites due to different abiotic conditions and that (2) prokaryotic microbiome composition would be influenced by the composition of epibiotic eukaryotes.

Materials and Methods

Study Area

We chose to perform our survey in the area around the island of Hiddensee on the German Baltic Sea coast as it offers variable yet locally representative environments colonized by eelgrass meadows. The area is characterized by semi-enclosed, shallow lagoons (German: “Bodden”) which differ from the open coast with respect to depth distribution, salinity and nutrient concentrations. The surveyed lagoons (the Vitter Bodden and Schaproder Bodden) feature somewhat lower average salinity (8.8 ± 1.0‰, range 6.5–13.4‰) compared to the open coast (Libben bay, 9.4 ± 1.6‰, range 6.8–15.9‰). Nutrient levels are elevated in lagoons (total N: 38.2 ± 12.3, total P: 1.2 ± 0.62 μmol l-1) compared to open coast waters (total N 19.9 ± 4.3, total P: 0.91 ± 0.4 μmol l-1) mainly due to agricultural runoff and other human activities (Schiewer, 2008). Salinity and nutrient values are yearly averages from monitoring data 2005 – 2014 (Landesamt für Umwelt, Naturschutz und Geologie Mecklenburg-Vorpommern, unpublished data).


Zostera marina shoots were sampled via scuba diving at seven sites around the island of Hiddensee (Figure 1). At every site, five replicate shoots (above-ground parts) were collected within an area of approximately 1 m2. Shoots were kept cool and in clean plastic bags until sample processing which was carried out within a few hours. Environmental variables such as water depth and parameters relating to the surrounding macrophyte vegetation (eelgrass and co-occurring macroalgae and aquatic plants) were determined as part of a simultaneous macrophyte inventory (Bühler, 2016).


FIGURE 1. Location of sampling sites around the island of Hiddensee on the German Baltic Sea coast. The inset shows the geographic placement of the area.

Leaf Biofilm Sample Processing

The sampled eelgrass shoots were rinsed with sterile filtered seawater and handled under sterile conditions. Leaf areas with heavy macroscopic epibiosis and visibly degrading tissue were removed before biofilm sampling. Epibiotic biofilm from the entire leaf surface of each shoot was scraped off using sterile cotton swabs. The scraped off material was suspended in sterile-filtered seawater (10 ml) and kept cool. The suspensions were aliquoted into separate tubes and centrifuged (10000 × g, 10 min, 4°C) to pellet biofilm material. Pellets were frozen at -20°C until DNA and chlorophyll a extraction. Chlorophyll a was extracted overnight in acetone at +4°C and absorbance at 665 and 750 nm (A665 and A750) was measured spectrophotometrically (Genesys 20, Thermo Fisher Scientific). Biofilm chlorophyll a content (Cchla) was calculated as Cchla = 11.4 (A665A750) and was expressed in μg chlorophyll a cm-2 of biofilm.

DNA Extraction and Illumina Amplicon Sequencing

DNA was extracted from biofilm pellets using the MoBio PowerSoil kit for DNA (MOBIO Laboratories). Mechanical lysis was achieved by bead beating in a FastPrep 24 5G (MP Biomedicals). Extracted DNA was amplified with primer pairs targeting the V4 region of the 16S rRNA gene [515f: 5′-GTGYCAGCMGCCGCGGTAA-3′, 806r: 5′-GGACTACNVGGGTWTCTAAT-3′ (Walters et al., 2016)] and the V7 region of the 18S rRNA gene [F-1183mod: 5′-AATTTGACTCAACRCGGG-3′, R-1443mod: 5′-GRGCATCACAGACCTG-3′ (Ray et al., 2016)] coupled to custom adaptor-barcode constructs. PCR amplification and Illumina MiSeq library preparation and sequencing (V3 chemistry) was carried out by LGC Genomics in Berlin. Sequences have been submitted to the NCBI short read archive under bioproject number PRJNA389390, and accession numbers SAMN07197184 – SAMN07197218 (16S rRNA gene sequences) and SAMN07197255 – SAMN07197289 (18S rRNA gene sequences).

Sequence Processing

Sequences clipped from adaptor and primer sequence remains were processed using the DADA2 package in R (version 1.2.0) (Callahan et al., 2016; R Development Core Team, 2017). Briefly, forward and reverse Illumina reads were truncated to 200 bp, filtered (maxEE = 2, truncQ = 2), dereplicated and error rates were estimated using the maximum possible error estimate from the data as an initial guess. Sample sequences were inferred, and paired forward and reverse reads were merged. Chimeric sequences were removed using the removeBimeraDenovo function. The resulting unique sequence variants (analogous to operational taxonomic units) were used to construct a table containing relative abundances of sequence variants across all samples. Sequence variants were taxonomically classified using a lowest common ancestor approach (Lanzén et al., 2012) based on the Silva database (Pruesse et al., 2007).

Statistical Analyses

Multivariate statistical analyses were carried out using functions in the vegan package [version 2.4-1, (Oksanen et al., 2016)] in R [version 3.3.1, (R Development Core Team, 2017)]. To visualize similarities in sequence variant composition (community composition) between sampling sites, non-metric multidimensional scaling (vegan function metaMDS) was performed on Hellinger-transformed sequence variant counts using Bray–Curtis distance. To explain the variation in community composition in response to the environment, abiotic and biotic variables were selected based on their explanatory power in PERMANOVA tests (vegan function adonis, variables were transformed (loge) to achieve normal-distribution if needed). A final model including the two most relevant environmental variables was selected and distance-based redundancy analysis ordination (vegan function capscale, Bray–Curtis distance) constrained to these environmental variables was used to visualize variation.

Sequence variant richness was calculated by rarefying read counts to the lowest number of reads in a sample of each dataset (16S: 5861, 18S: 25780). Evenness was calculated as J = H/logS, where H is the Shannon diversity index and S is rarefied richness (Pielou, 1977). To identify sequence variants that were differently abundant in lagoon sites and open coast sites, differential abundance analysis as implemented in the DEseq2 R package (Love et al., 2014).

A partial mantel test was used to assess correlation between prokaryotic and eukaryotic community composition while correcting for environmental variables [vegan function mantel.partial, Bray–Curtis distance for community data, euclidean distance for environmental variables, (Legendre and Legendre, 2012)]. To analyze pairwise correlations between 16S and 18S sequence variants, co-occurrence networks were calculated in R using Pearson correlation. P-values were adjusted for multiple testing using the Benjamini–Hochberg method (Benjamini and Hochberg, 1995). Only correlations with adjusted p-values < 0.01 were treated as significant and included in the network. The network was plotted and edited for clarity (symbols, colors etc.) in Cytoscape version 3.3.0 (Shannon et al., 2003).


Due to the different depth distributions of eelgrass on open coasts and in lagoons (Bühler, 2016), all open coast sites where eelgrass was found were deeper than lagoon sites. Eelgrass biomass (leaf surface area and dry weight) varied widely among sampling sites (Table 1). Leaf surface biofilms were visually different in color and density (results not shown) and biofilm chlorophyll a content varied substantially among sampling sites (Table 1).


TABLE 1. An overview of the sampling sites and key abiotic and biotic parameters.

Illumina amplicon sequencing of 16S and 18S rRNA fragments resulted in 1.0 million 16S reads and 1.9 million 18S reads with an average of 29564 16S reads (min = 5861, max = 67456) and an average of 56886 18S reads per sample (min = 25780 max = 126938). One sample from site SBL was excluded from all analyses because of low 16S read count (861 reads). In total, 4409 prokaryotic and 1946 eukaryotic sequence variants (analogous to operational taxonomic units) were identified. Among prokaryotes, Bacteria were overwhelmingly abundant over Archaea (99.2% of sequence variants and 99.9% of reads). Due to this, 16S sequence variants are referred to as “bacteria” from here on.

Bacterial community composition varied substantially among sites (PERMANOVA R2 = 0.82, p < 0.001), and the main variation was between lagoon sites and open coast sites (PERMANOVA R2 = 0.42, p < 0.001, Figure 2A). Eukaryotic community composition also varied among individual sites (PERMANOVA R2 = 0.83, p < 0.001) and varied strongly between lagoon and open coast sites (PERMANOVA R2 = 0.48, p < 0.001, Figure 2B). Water depth, eelgrass shoot leaf area and biofilm chlorophyll a content were the strongest environmental predictors of community variation in both bacteria and eukaryotes, with depth explaining about 30%, leaf area 11% and chlorophyll a 10% of total variation (PERMANOVA p < 0.001). There was also a significant interaction term between depth and chlorophyll a content for both bacterial and eukaryote communities (PERMANOVA, p < 0.01), indicating that chlorophyll a affects communities differently depending on depth, or vice-versa. The variation explained by depth and chlorophyll a content is illustrated in Figures 2C,D.


FIGURE 2. Ordinations based on 16S rRNA (A,C) and 18S rRNA (B,D) gene sequence variant relative abundances illustrate the variability in bacterial and eukaryotic community composition between sampling sites. In (A,B), non-metric multidimensional scaling (nMDS) display the distinct clustering of lagoon (blue ellipse) and open coast sites (green ellipse). Ellipses represent 95% confidence intervals of site position. In (C,D), distance-based RDA (Bray–Curtis distance metric used) constrained to water depth and leaf surface biofilm chlorophyll a content show the variability that can be explained by these two environmental variables.

Richness of bacterial sequence variants was positively correlated with eelgrass dry weight at the sampling site (R2 = 0.69, p < 0.001, Figure 3A) as well as with eelgrass leaf area, although somewhat more weakly (R2 = 0.29, p < 0.001). Eukaryotic sequence variant richness did not correlate with eelgrass dry weight or other environmental gradients, yet was significantly higher in lagoon sites compared to open coast sites (ANOVA F = 208.6, p < 0.001, Figure 3B). Richness and evenness were correlated for both bacteria (R2 = 0.34, p < 0.001) and eukaryotes (R2 = 0.84, p < 0.001).


FIGURE 3. Richness of 16S (A) and 18S (B) rRNA gene sequence variants in relation to the biomass (dry weight) of Zostera marina at all sampling sites. Symbols and colors are the same as in Figure 1.

A differential abundance analysis revealed sequence variants that differed significantly (adjusted p < 0.01) in relative abundance between lagoon and outer coast sites (Figure 4). On a broad taxonomical level, several bacterial sequence variants from Cyanobacteria, Planctomycetes, Gemmatimonadetes, and Acidobacteria were overrepresented in lagoon sites (Figure 4A). Among eukaryotic taxa, one highly abundant individual sequence variant, classified as Hydrozoa (Metazoa: Cnidaria) was significantly overrepresented in open coast sites, while several other metazoan (e.g., classified as Dorvillea and Rotifera) and diatom sequence variants (Cocconeidaceae) were overrepresented in lagoon sites. Fourteen bacterial and 3 eukaryotic sequence variants were detected in all samples, and were defined as the “core community” of eelgrass leaf biofilms for the purpose of this study (Table 2).


FIGURE 4. Differential abundance analysis of 16S (A) and 18S (B) rRNA gene sequence variants comparing open coast and lagoon sites. Positive log2 fold-change values indicate a significantly (adjusted p < 0.01) higher abundance in open coast sites, while negative values indicate significantly higher abundances in lagoon sites. The area of each circle representing an individual sequence variant is proportional to the relative abundance of that sequence variant across all samples. Prominent eukaryotic sequence variants are labeled with more detailed taxonomic classification.


TABLE 2. Mean relative abundance (across all samples) and taxonomic classifications of bacterial and eukaryotic rRNA gene sequence variants detected in all samples, i.e., the “core community” of eelgrass leaf biofilms in the present study.

A partial Mantel test, correcting for all measured environmental variables, confirmed a strong correlation between bacterial and eukaryote community composition (partial mantel r = 0.92, p = 0.001), which was suggested by the similar clustering in the nMDS ordinations Figures 2A,B. To further explore correlations between individual bacterial and eukaryotic sequence variants, a co-occurrence network was constructed (Figure 5). The network features several prominent sequence variants that correlate with many other taxa and in some cases form distinct network modules. In particular, a sequence variant classified as Hydrozoa (Metazoa: Cnidaria) displayed both positive (co-occurrence) and negative (mutual exclusion) correlations with several bacterial taxa. Other prominent eukaryotic taxa were classified as Mytilus sp. (mussel), Dorvillea sp. (polychaete) and Bacillariophyta (diatoms). In some cases, network modules centered around bacteria, such as Bacteroidetes, Proteobacteria, and Cyanobacteria.


FIGURE 5. Co-occurrence of bacterial and eukaryotic rRNA gene sequence variants presented as a network. Solid lines (network edges) indicate significant positive Pearson correlations (r > 0.8, p < 0.01) between bacterial (circles) and eukaryotic (triangles) rRNA gene sequence variants (i.e., taxa, network nodes) while dashed lines indicate significant negative correlations (r > –0.8, p < 0.01). Nodes are sized according to the number of correlating neighbors (degree) and color-coded according to taxonomic affiliation. Selected prominent nodes (same as in Figure 4) are further labeled with more detailed taxonomic classifications.


Eelgrass Bacterial Microbiomes Are Locally Variable

This study represents the first effort to characterize both bacterial and eukaryotic epibiotic communities on eelgrass leaves and to relate their communities to potentially important environmental drivers and to each other. We found that the local abiotic environment, in the form of water depth, explained a great deal of variation for both communities. Water depth was one of the main differences between lagoon and open coast sites which likely in part explains the clear separation observed between these sites. However, water depth encompasses separate physical factors such as light penetration and wave shear forces which are difficult to further separate in this study. In addition, several other factors such as mean salinity, nutrient levels and water temperature differ on seasonal time scales in lagoons compared to the open coast (Schiewer, 2008; Blindow et al., 2016) and could also contribute to this variation.

Of the biotic environmental variables measured, eelgrass shoot leaf area (a measure of the size of each shoot) and biofilm chlorophyll a content explained the most variation in community composition of both bacteria and eukaryotes. The explanatory power of eelgrass leaf area, which is related to eelgrass productivity (Echavarría-Heras et al., 2010), suggests an influence of host physiology on epibiotic biofilm composition. The amount of variation explained by chlorophyll a content, a proxy for the contribution of photosynthetic organisms (algae and cyanobacteria) in the biofilms, further indicates that interactions within the biofilms, for example between eukaryotic algae and bacteria play an important role in shaping community composition.

Richness of bacterial sequence variants was correlated with leaf area and with eelgrass dry weight at the sampling sites. Larger leaf area corresponds directly to a larger area of biofilm sampled, since biofilm from entire eelgrass shoots was removed. This observed richness-area relationship agrees with similar biodiversity-area relationships in other microbial habitats (reviewed in Martiny et al., 2006) and among macroscopic organisms. However, an alternative explanation for this pattern relates to the classical notion that energy supply limits biodiversity (Currie, 1991). Since leaf area is a proxy for eelgrass productivity (Echavarría-Heras et al., 2010), its correlation with bacterial diversity could indicate that higher resource supply in the form of for example organic exudates promotes leaf surface bacterial diversity. In contrast, eukaryotic biodiversity (measured as sequence variant richness) seemed to be controlled by factors specific to lagoon vs. open coast habitats, rather than area.

Despite the small geographic area surveyed in this study, variability in community composition of both bacteria and eukaryotes was substantial among sites, which shows that eelgrass leaf biofilms are dynamic habitats on a regional scale. This observation is consistent with, and complementary to, recent results that revealed that eelgrass leaf surface microbiomes were more variable than root-associated microbiomes on a global scale (Fahimipour et al., 2017), as well as within an individual eelgrass patch (Ettinger et al., 2017). Nevertheless, we detected several sequence variants in all samples analyzed, supporting the possible existence of a core microbiome of eelgrass leaves. Several core sequence variants were classified to the family Rhodobacteraceae (Alphaproteobacteria), which are often surface-attached, chemoorganotrophic bacteria common in the marine environment (Dang and Lovell, 2002). Interestingly, this family was also consistently detected in the rhizospheres of eelgrass and its close relative Z. noltii (Cúcio et al., 2016; Ettinger et al., 2017) as well as on leaves and below-ground parts of the seagrass Halophila stipulacea in the red sea (Mejia et al., 2016) indicating that their members indeed play important roles in association with various seagrasses. The observation of a Methylophilaceae (Betaproteobacteria) sequence variant as an abundant core community member agrees with detection of this taxon on eelgrass and other co-occurring macrophytes in the brackish-water Chesapeake Bay (Crump and Koch, 2008) and suggests importance of one-carbon metabolism on eelgrass leaf surfaces (Lapidus et al., 2011).

Eelgrass Bacterial Microbiomes Correlate with Eukaryotic Epibiont Communities

The strong correlation between the bacterial and eukaryotic communities observed in this study has at least two possible causes. First, the abiotic and biotic environment experienced by both communities could influence them in similar ways. This explanation is partially supported by the similar amounts of variation explained by water depth, leaf area and chlorophyll a content for both communities. Second, bacterial and eukaryotic taxa may interact directly, resulting in a bacterial community that is driven by the composition of the eukaryotic community, or vice versa. The topology of the co-occurrence network (Figure 5) suggests some importance of the latter scenario, where certain eukaryotic taxa (e.g., classified as Hydrozoa, Mytilus sp. and Cocconeidaceae) feature a disproportionate number of significant correlations with bacterial taxa. This indicates that bacterial communities in mature biofilms on eelgrass leaves are to some extent a reflection of the combined microbiomes of several epibiotic eukaryotes such as hydrozoans, mussels and diatoms in addition to any eelgrass-specific microbiome that may exist. Conversely, some bacterial taxa correlate with several eukaryotic taxa. For example, a cyanobacterial sequence variant classified as Rivularia sp. correlates with several ciliate taxa which may for example be attracted to this cyanobacterium as grazers. In a previous study, bacterial and eukaryotic communities of the tropical seagrass Enhalus acoroides were found to respond similarly to environmental CO2 exposure, suggesting that seagrass bacterial and eukaryotic epibiotic communities are highly correlated also in other seagrass species and environments. However, potential direct interactions between bacterial and eukaryotic taxa were not addressed in this study (Hassenrück et al., 2015).


To date, many studies addressing the microbiomes of marine organisms have focused exclusively on prokaryote communities (see e.g., Hentschel et al., 2012; Wahl et al., 2012; Bourne et al., 2016 and references therein). In contrast, free-living microbial eukaryotes are receiving increasing attention and molecular community analysis has led to several recent discoveries that emphasize their underexplored diversity and crucial roles in the environment (Liu et al., 2009; Lima-Mendez et al., 2015). It is likely that insights of a similar caliber can be gained by applying molecular tools to eukaryotes associated to macroscopic hosts (Andersen et al., 2013), beyond relatively well-described eukaryotic symbionts such as zooxanthella in corals (Muller-Parker et al., 2015). We have shown that the bacterial microbiome is strongly correlated with epibiotic eukaryote communities, including microbial and macroscopic taxa, indicating that interactions between eukaryotes and bacteria drive community assembly on eelgrass surfaces to some extent. By ignoring eukaryotes in bacterial microbiome studies, a large part of the variation in bacterial communities could remain unexplained or misinterpreted. Indeed, the high diversity of both eukaryotic primary producers such as diatoms and heterotrophic taxa such as ciliates, rhizarians and many macroscopic and microscopic animals paint a picture of complex microbial ecosystems encompassing several trophic levels on eelgrass leaves. In these “microbial jungles,” bacteria are bound to play important roles which may be especially relevant in interaction with their eukaryote neighbors.

In addition to considering the understudied eukaryotes in microbiome studies, gaining a mechanistic understanding into microbiome assembly requires means to infer causative relationships between microbiome dynamics and abiotic and biotic factors (Widder et al., 2016). The present study included a limited number of samples from a small geographic area, providing a valuable first assessment of environmental drivers of community assembly and potential microbial interactions. Further insight could be gained by manipulative experiments, whereby for example, seagrass surfaces are exposed to colonization by microorganisms under controlled conditions. In addition, large observational microbiome datasets from a range of eelgrass meadows encompassing wider environmental gradients would allow more statistical power to tease apart different drivers of community assembly. Modeling approaches such as structural equation modeling (Pugesek et al., 2003) could allow assessment of causative relationships between the environment, the seagrass host and its microbiome. A major challenge lies in linking microbiome dynamics to processes that are relevant on the ecosystem level. Manipulative experiments that operate on an ecosystem scale, i.e., mesocosm experiments mimicking seagrass meadows, offer a promising approach to manipulate the environment while tracking microbiome dynamics and their links to critical ecosystem fluxes such as primary production and respiration. Ultimately, understanding assembly and function of seagrass microbiomes could contribute to tackling some of the challenges faced by seagrass meadow ecosystems under global climate change and other anthropogenic stressors.

Author Contributions

MB conceived the study, participated in the field work, processed the eelgrass samples, led the laboratory work, analyzed the data and wrote the paper. ABü planned and executed the scuba diving, gathered the macrophyte data and prepared Figure 1. ABr performed molecular laboratory work and contributed to data analysis. SD participated in planning and execution of the field work and drove the research vessel. HS planned the fieldwork. IB planned and supervised the fieldwork and macrophyte data collection. All authors edited the manuscript.


The field sampling campaign was funded through a grant from the German Federal Ministry of Education and Research (Project BACOSA, 03F0665C) awarded to HS and IB. Microbial community analysis was funded through the University of Greifswald.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The National Park administration (Nationalparkamt Vorpommern) kindly permitted us to dive and sample within the lagoon area. Tim Urich provided essential infrastructure and support for the microbial community analysis.


Alsterberg, C., Eklöf, J. S., Gamfeldt, L., Havenhand, J. N., and Sundbäck, K. (2013). Consumers mediate the effects of experimental ocean acidification and warming on primary producers. Proc. Natl. Acad. Sci. U.S.A. 110, 8603–8608. doi: 10.1073/pnas.1303797110/-/DCSupplemental

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersen, L. O., Nielsen, H. V., and Stensvold, C. R. (2013). Waiting for the human intestinal Eukaryotome. ISME J. 7, 1253–1255. doi: 10.1038/ismej.2013.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Baden, S., Boström, C., Tobiasson, S., Arponen, H., and Moksnes, P.-O. (2010). Relative importance of trophic interactions and nutrient enrichment in seagrass ecosystems: a broad-scale field experiment in the Baltic–Skagerrak area. Limnol. Oceanogr. 55, 1435–1448. doi: 10.4319/lo.2010.55.3.1435

CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300. doi: 10.2307/2346101

CrossRef Full Text | Google Scholar

Blindow, I., Dahlke, S., Dewart, A., Flügge, S., Hendreschke, M., Kerkow, A., et al. (2016). Long-term and interannual changes of submerged macrophytes and their associated diaspore reservoir in a shallow southern Baltic Sea bay: influence of eutrophication and climate. Hydrobiologia 778, 121–136. doi: 10.1007/s10750-016-2655-4

CrossRef Full Text | Google Scholar

Boström, C., Baden, S., Bockelmann, A.-C., Dromph, K., Fredriksen, S., Gustafsson, C., et al. (2014). Distribution, structure and function of Nordic eelgrass (Zostera marina) ecosystems: implications for coastal management and conservation. Aquat. Conserv. Mar. Freshw. Ecosyst. 24, 410–434. doi: 10.1002/aqc.2424

PubMed Abstract | CrossRef Full Text | Google Scholar

Bourne, D. G., Morrow, K. M., and Webster, N. S. (2016). Insights into the coral microbiome: underpinning the health and resilience of reef ecosystems. Annu. Rev. Microbiol. 70, 317–340. doi: 10.1146/annurev-micro-102215-095440

PubMed Abstract | CrossRef Full Text | Google Scholar

Brakel, J., Werner, F. J., Tams, V., Reusch, T. B. H., and Bockelmann, A.-C. (2014). Current European Labyrinthula zosterae are not virulent and modulate seagrass (Zostera marina) defense gene expression. PLoS ONE 9:e92448. doi: 10.1371/journal.pone.0092448

PubMed Abstract | CrossRef Full Text | Google Scholar

Brodersen, K. E., Lichtenberg, M., Paz, L.-C., and Kühl, M. (2015). Epiphyte-cover on seagrass (Zostera marina L.) leaves impedes plant performance and radial O2 loss from the below-ground tissue. Front. Mar. Sci. 2:58. doi: 10.3389/fmars.2015.00058

CrossRef Full Text | Google Scholar

Bühler, A. (2016). Distribution of Zostera marina in the Coastal Waters of the Baltic Sea. Master thesis, University of Rostock, Rostock, 1–70.

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Celdrán, D., Espinosa, E., Sánchez-Amat, A., and Marín, A. (2012). Effects of epibiotic bacteria on leaf growth and epiphytes of the seagrass Posidonia oceanica. Mar. Ecol. Prog. Ser. 456, 21–27. doi: 10.3354/meps09672

CrossRef Full Text | Google Scholar

Crump, B. C., and Koch, E. W. (2008). Attached bacterial populations shared by four species of aquatic angiosperms. Appl. Environ. Microbiol. 74, 5948–5957. doi: 10.1128/AEM.00952-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Cúcio, C., Engelen, A. H., Costa, R., and Muyzer, G. (2016). Rhizosphere microbiomes of European seagrasses are selected by the plant, but are not species specific. Front. Microbiol. 7:440. doi: 10.3389/fmicb.2016.00440

PubMed Abstract | CrossRef Full Text | Google Scholar

Currie, D. J. (1991). Energy and large-scale patterns of animal-species and plant-species richness. Am. Nat. 137, 27–49. doi: 10.1086/285144

CrossRef Full Text | Google Scholar

Dang, H., and Lovell, C. (2002). Numerical dominance and phylotype diversity of marine Rhodobacter species during early colonization of submerged surfaces in coastal marine waters as determined by 16S ribosomal DNA sequence analysis and fluorescence in situ hybridization. Appl. Environ. Microbiol. 68, 496–504. doi: 10.1128/AEM.68.2.496-504.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Duffy, J. E., Reynolds, P. L., Boström, C., Coyer, J. A., Cusson, M., Donadi, S., et al. (2015). Biodiversity mediates top-down control in eelgrass ecosystems: a global comparative-experimental approach. Ecol. Lett. 18, 696–705. doi: 10.1111/ele.12448

PubMed Abstract | CrossRef Full Text | Google Scholar

Echavarría-Heras, H., Solana-Arellano, E., and Franco-Vizcaíno, E. (2010). An allometric method for the projection of eelgrass leaf biomass production rates. Math. Biosci. 223, 58–65. doi: 10.1016/j.mbs.2009.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Eklöf, J. S., Alsterberg, C., Havenhand, J. N., Sundbäck, K., Wood, H. L., and Gamfeldt, L. (2012). Experimental climate change weakens the insurance effect of biodiversity. Ecol. Lett. 15, 864–872. doi: 10.1111/j.1461-0248.2012.01810.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ettinger, C. L., Voerman, S. E., Lang, J. M., Stachowicz, J. J., and Eisen, J. A. (2017). Microbial communities in sediment from Zostera marina patches, but not the Z. marina leaf or root microbiomes, vary in relation to distance from patch edge. PeerJ 5:e3246. doi: 10.7717/peerj.3246

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahimipour, A. K., Kardish, M. R., Lang, J. M., Green, J. L., Eisen, J. A., and Stachowicz, J. J. (2017). Global-scale structure of the eelgrass microbiome. Appl. Environ. Microbiol. 83:e03391-16. doi: 10.1128/AEM.03391-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Fourqurean, J. W., Duarte, C. M., Kennedy, H., Marbà, N., Holmer, M., Mateo, M. A., et al. (2012). Seagrass ecosystems as a globally significant carbon stock. Nat. Geosci. 5, 505–509. doi: 10.1038/ngeo1477

CrossRef Full Text | Google Scholar

Hassenrück, C., Hofmann, L. C., Bischof, K., and Ramette, A. (2015). Seagrass biofilm communities at a naturally CO2-rich vent. Environ. Microbiol. Rep. 7, 516–525. doi: 10.1111/1758-2229.12282

PubMed Abstract | CrossRef Full Text | Google Scholar

Hentschel, U., Piel, J., Degnan, S. M., and Taylor, M. W. (2012). Genomic insights into the marine sponge microbiome. Nat. Rev. Microbiol. 10, 641–654. doi: 10.1038/nrmicro2839

PubMed Abstract | CrossRef Full Text | Google Scholar

Huws, S. A., McBain, A. J., and Gilbert, P. (2005). Protozoan grazing and its impact upon population dynamics in biofilm communities. J. Appl. Microbiol. 98, 238–244. doi: 10.1111/j.1365-2672.2004.02449.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lanzén, A., Jorgensen, S. L., Huson, D. H., Gorfer, M., Grindhaug, S. H., Jonassen, I., et al. (2012). CREST – classification resources for environmental sequence tags. PLoS ONE 7:e49334. doi: 10.1371/journal.pone.0049334

PubMed Abstract | CrossRef Full Text | Google Scholar

Lapidus, A., Clum, A., Labutti, K., Kaluzhnaya, M. G., Lim, S., Beck, D. A. C., et al. (2011). Genomes of three methylotrophs from a single niche reveal the genetic and metabolic divergence of the Methylophilaceae. J. Bacteriol. 193, 3757–3764. doi: 10.1128/JB.00404-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Legendre, P., and Legendre, L. (2012). Numerical Ecology. Amsterdam: Elsevier.

Google Scholar

Lima-Mendez, G., Faust, K., Henry, N., Decelle, J., Colin, S., Carcillo, F., et al. (2015). Ocean plankton. Determinants of community structure in the global plankton interactome. Science 348:1262073. doi: 10.1126/science.1262073

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Probert, I., Uitz, J., Claustre, H., Aris-Brosou, S., Frada, M., et al. (2009). Extreme diversity in noncalcifying haptophytes explains a major pigment paradox in open oceans. Proc. Natl. Acad. Sci. U.S.A. 106, 12803–12808. doi: 10.1073/pnas.0905841106

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Martiny, J. B. H., Bohannan, B. J. M., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341

PubMed Abstract | CrossRef Full Text | Google Scholar

Maxwell, P. S., Eklöf, J. S., van Katwijk, M. M., O’Brien, K. R., de la Torre-Castro, M., Boström, C., et al. (2016). The fundamental role of ecological feedback mechanisms for the adaptive management of seagrass ecosystems – a review. Biol. Rev. 92, 1521–1538. doi: 10.1111/brv.12294

PubMed Abstract | CrossRef Full Text | Google Scholar

Mejia, A. Y., Rotini, A., Lacasella, F., Bookman, R., Thaller, M. C., Shem-Tov, R., et al. (2016). Assessing the ecological status of seagrasses using morphology, biochemical descriptors and microbial community analyses. A study in Halophila stipulacea (Forsk.) Aschers meadows in the northern Red Sea. Ecol. Indic. 60, 1150–1163. doi: 10.1016/j.ecolind.2015.09.014

CrossRef Full Text | Google Scholar

Mieszkin, S., Callow, M. E., and Callow, J. A. (2013). Interactions between microbial biofilms and marine fouling algae: a mini review. Biofouling 29, 1097–1113. doi: 10.1080/08927014.2013.828712

PubMed Abstract | CrossRef Full Text | Google Scholar

Muehlstein, L. K., Porter, D., and Short, F. T. (1991). Labyrinthula zosterae sp. nov., the causative agent of wasting disease of eelgrass. Zostera marina. Mycologia 83, 180–191. doi: 10.2307/3759933

CrossRef Full Text | Google Scholar

Muller-Parker, G., D’Elia, C. F., and Cook, C. B. (2015). “Interactions between corals and their symbiotic algae,” in Coral Reefs in the Anthropocene, ed. C. Birkeland (Dordrecht: Springer), 99–116. doi: 10.1007/978-94-017-7249-5_5

CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2016). vegan: Community Ecology Package. R Package Version 2.2-1.

Google Scholar

Orth, R. J., Carruthers, T. J. B., Dennison, W. C., Duarte, C. M., Fourqurean, J. W., Heck, K. L. J., et al. (2006). A global crisis for seagrass ecosystems. Bioscience 56, 987–996. doi: 10.1641/0006-3568 (2006)56[987:AGCFSE]2.0.CO;2

CrossRef Full Text | Google Scholar

Pielou, E. C. (1977). Mathematical Ecology. New York, NY: Wiley.

Google Scholar

Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864

PubMed Abstract | CrossRef Full Text | Google Scholar

Pugesek, B. H., Tomer, A., and von Eye, A. (2003). Structural Equation Modeling. Cambridge: Cambridge University Press. doi: 10.1017/CBO9780511542138

CrossRef Full Text | Google Scholar

R Development Core Team (2017). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Ray, J. L., Althammer, J., Skaar, K. S., Simonelli, P., Larsen, A., Stoecker, D., et al. (2016). Metabarcoding and metabolome analyses of copepod grazing reveal feeding preference and linkage to metabolite classes in dynamic microbial plankton communities. Mol. Ecol. 25, 5585–5602. doi: 10.1111/mec.13844

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiewer, U. (ed.). (2008). Ecology of Baltic Coastal Waters. Berlin: Springer Science & Business Media, doi: 10.1007/978-3-540-73524-3

CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Steinberg, P., and de Nys, R. (2002). Chemical mediation of colonization of seaweed surfaces. J. Phycol. 38, 621–629. doi: 10.1046/j.1529-8817.2002.02042.x

CrossRef Full Text | Google Scholar

Törnblom, E., and Søndergaard, M. (1999). Seasonal dynamics of bacterial biomass and production on eelgrass Zostera marina leaves. Mar. Ecol. Prog. Ser. 179, 231–240. doi: 10.3354/meps179231

CrossRef Full Text | Google Scholar

Wahl, M., Goecke, F., Labes, A., Dobretsov, S., and Weinberger, F. (2012). The second skin: ecological role of epibiotic biofilms on marine organisms. Front. Microbiol. 3:292. doi: 10.3389/fmicb.2012.00292

PubMed Abstract | CrossRef Full Text | Google Scholar

Walters, W., Hyde, E. R., Berg-Lyons, D., Ackermann, G., Humphrey, G., Parada, A., et al. (2016). Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems 1:e00009-15. doi: 10.1128/mSystems.00009-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Widder, S., Allen, R. J., Pfeiffer, T., Curtis, T. P., Wiuf, C., Sloan, W. T., et al. (2016). Challenges in microbial ecology: building predictive understanding of community function and dynamics. ISME J. 10, 2557–2568. doi: 10.1038/ismej.2016.45

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: seagrass, microbial interactions, microbial eukaryotes, community assembly, biofilms, Baltic Sea, microalgae, epibiosis

Citation: Bengtsson MM, Bühler A, Brauer A, Dahlke S, Schubert H and Blindow I (2017) Eelgrass Leaf Surface Microbiomes Are Locally Variable and Highly Correlated with Epibiotic Eukaryotes. Front. Microbiol. 8:1312. doi: 10.3389/fmicb.2017.01312

Received: 24 February 2017; Accepted: 28 June 2017;
Published: 13 July 2017.

Edited by:

Télesphore Sime-Ngando, Centre National de la Recherche Scientifique (CNRS), France

Reviewed by:

Susannah Green Tringe, Joint Genome Institute (DOE), United States
Lucas Stal, Royal Netherlands Institute for Sea Research (NWO), Netherlands

Copyright © 2017 Bengtsson, Bühler, Brauer, Dahlke, Schubert and Blindow. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Mia M. Bengtsson,