BRIEF RESEARCH REPORT article

Front. Mar. Sci., 11 October 2019

Sec. Marine Molecular Biology and Ecology

Volume 6 - 2019 | https://doi.org/10.3389/fmars.2019.00635

Microbiota Differences of the Comb Jelly Mnemiopsis leidyi in Native and Invasive Sub-Populations

  • 1. GEOMAR Helmholtz Centre for Ocean Research Kiel, Marine Evolutionary Ecology, Kiel, Germany

  • 2. Institute for General Microbiology, Christian-Albrechts-Universität zu Kiel, Kiel, Germany

  • 3. Max Planck Institute for Evolutionary Biology, Plön, Germany

Abstract

The translocation of non-indigenous species around the world, especially in marine systems, is a matter of concern for biodiversity conservation and ecosystem functioning. While specific physical and ecological traits are often recognized to influence the success in the establishment of non-indigenous species, the impact of the associated microbiota for the fitness and performance of basal marine metazoans remains vastly unknown. However, for the elucidation of the microbial impact on host fitness, first the diversity of the associated microbiota has to be characterized. In this study, we compared bacterial composition patterns of the invasive ctenophore Mnemiopsis leidyi in different native and invasive sub-populations along with the genetic structure of the host (polymorphic microsatellite markers). Amplicon sequencing of the 16S rRNA gene (V1-V2 hypervariable regions) revealed that M. leidyi as representative of the phylum Ctenophora, the sister group to all metazoans, harbored a distinct microbiota compared to the ambient seawater, which significantly differed across two major tissues, namely epidermis and gastrodermis. Moreover, we identified significant differences in bacterial community compositions between native and invasive sub-populations of M. leidyi. This might indicate, that the microbiota can be influenced by the genotypic background of the host. To detangle such interactions, laboratory controlled experiments in a common garden set up are needed.

Introduction

All multicellular organisms are intimately associated with microbes and thus jointly constitute metaorganisms, a fact that is attaining large public and scientific attention. For model organisms it has been shown that the microbiota contributes to host fitness and can be essential for host functioning (; ; ). Ctenophores, being a basal metazoan phylum (; ), have so far not been considered in the metaorganism context (). However, it is expected that a more detailed understanding on host-microbe interactions at the base of the metazoan tree can reveal important information about metaorganism evolution. For example, in other marine basal metazoans such as corals and jellyfish, bacteria and their exudates have a crucial impact on life history decision such as metamorphosis and settlement induction to the benthos (; ). Additionally, bacterial exudates have been documented to trigger sexual reproduction in choanoflagellates () as well as the evolution of multicellularity ().

The ctenophore Mnemiopsis leidyi is not only a basal metazoan species, but is becoming an ecological model organism for biological invasions () with large and widespread impacts on invaded ecosystems (e.g., ; ). For western Eurasia, two distinct invasion events can be differentiated with a first introduction into the Black Sea with subsequently spreading throughout adjacent waters including the Mediterranean Sea, and a more recent introduction into north-western Europe (; ). Both introductions can be traced to different sub-population origins in the native habitat with the north European invasion originating from the north-east coast of the United States, while the southern invasion can be traced back to the Gulf of Mexico and south-east coast of the United States (; ). Some earlier studies already characterized the associated microbiome of M. leidyi in invaded () and native () sub-populations, but came to inconclusive results. Moreover, the study by even claimed that M. leidyi epithelia in invaded habitats do not contain an associated microbiome (), which was, however, not confirmed by others (). Thus, the objectives of this study are to investigate (1) the effect of M. leidyi sub-population origin on the microbiota diversity, (2) the bacterial community differentiation between specific M. leidyi tissue compartments and 3) assess differences in the microbiome of native versus invasive sub-populations by using standardized sampling and sequencing methods.

Materials and Methods

Sampling of Mnemiopsis leidyi

Specimen of M. leidyi (8 to 17 individuals) were sampled during summer 2012 (approximately during temperature peak) in northern and southern areas of its native as well as invasive distribution range (see ). In more detail, samples from native habitats originated from Mosquito Lagoon, Florida (17 ctenophores, 11 water samples), Atlantic Coast of south-east United States; Woods Hole (Cape Cod), Massachusetts, Atlantic Coast of north-east United States (10 ctenophores, nine water samples); whereas invasive M. leidyi were sampled in Marseille, France, north-west Mediterranean Sea (10 ctenophores, six water samples); and Helgoland, Germany, German Bight/North Sea (eight ctenophores, five water samples). Animals were caught by use of a dip net and animals were individually kept in sterile filtered seawater for up to 6 h before taxidermy of epi- and gastrodermis was conducted. All animals were ensured to have empty guts upon dissection and were washed thoroughly with sterile filtered seawater to remove non-associated microbes. Epi- and gastrodermis samples were attained by cutting approximately 1 × 1 cm of the respective tissue out of the animals on sterile petri dishes using a sterile scalpel. Excised samples comprising epithelial tissue as well as mucus were transferred into reaction tubes and homogenized by use of sterile pipet tips for later DNA extraction resulting in isolation of 1 mL v/v samples from the outer epithelia and 200 μL v/v samples from the gastric cavity. Ambient seawater was collected from the sampling sites (5 to 11 replicates per site) by bucket sampling and 1L was subsequently and successively filtered through 10 μm isopore (TCTP 04700) membranes and 0.22 μm polyethersulfone (GPWP 04700) membrane filters (Millipore Corporation, Bedford, MA, United States). All samples were immediately frozen at −80°C until processing.

Nucleic Acid Isolation and 16S rRNA Gene Amplicon Sequencing

DNA was isolated from two different tissue types of M. leidyi, namely mucus and epidermis tissue samples from the outer epithelia (1 mL v/v) as well as tissue samples from the gastric cavity (pharynx: 200 μL v/v), and ambient seawater (from filters - 0.22 μm). DNA was isolated using the Wizard Genomic DNA Purification Kit (Promega, Madison, United States) according the manufacturers protocol. PCR amplification of the V1–V2 hypervariable region of the 16S rRNA gene was performed with primers V2_A_Pyro_27F (5′-CGTATCGCCTCCCTCGCGCCATCAGTCAGAGTTTGATCCTGGCTCAG-3′) and V2_B_Pyro_27F (5′-CTATGCGCCTTGCCAGCCCGCTCAGTCAGAGTTTGATCCTGGCTCAG-3′). Amplicons were size checked and purified using MinElute Gel Extraction Kit (Qiagen, Hilden, Germany). Purified amplicon DNAs were quantified using Quant-iT PicoGreen kit (Thermo Fisher Scientific, Waltham, United States) and equal amounts of purified PCR products were pooled for subsequent sequencing using the Illumina MiSeq platform with library preparation according to manufacturer’s instructions at the Max-Planck Institute for Evolutionary Biology.

Data Processing and Bioinformatics

All steps of sequence processing were conducted with the program mothur v1.27.0 (). The Greengenes reference taxonomy (), available on http://www.mothur.org/wiki/Greengenes-formatted_databases, was used with reference sequences trimmed to the V1-V2 hypervariable region of the 16S rRNA gene to improve accuracy of classification (). The sequence depth of all 117 analyzed samples ranged from 23,622 to 675,908 single sequences. A random subset of 23,622 sequences per sample (corresponding to the smallest number of reads across samples) was generated to eliminate bias due to unequal sampling effort. The aligned dataset was used to compute a distance matrix (dist.seqs) for binning sequences into operational taxonomic units (OTU) by average neighbor clustering (cluster.split). 2,829 OTUs at a 97% similarity threshold (roughly corresponding to species level) distinction were considered. All downstream computations were performed in R v2.15.1 (), following published protocols (). Briefly, the extent of change in OTU composition was explored by beta diversity analysis in terms of redundancy analysis (RDA) with Hellinger-transformed OTU abundances (; ; ). In general, the RDA model was formulated as “OTU abundances F,” with factor F being either “provenance,” or “compartment.” Models were evaluated with the vegan function rda, followed by anova.cca, using 10,000 random permutations. Homoscedasticity was checked by using the vegan betadisper function, followed by anova.betadisper to assess significance. Upon determination of the significance (P < 0.05) of the overall model, pairwise factor-level combinations were tested for significance. OTUs significantly correlated with any axis in a significant RDA model were determined by using the envfit function with 105 permutations, followed by Benjamini-Hochberg correction.

Sequence data were submitted to the NCBI Sequence Read Archive under Accession No. PRJNA548849 and are additionally stored in a sequence data management system at the Institute of Clinical Molecular Biology (IKMB) at Kiel University, Germany.

Host Genotyping

Genomic DNA from M. leidyi epithelia tissue samples, isolated with the Wizard Genomic DNA Purification Kit (Promega, Madison, United States), were used for genotype analyses. Sampled M. leidyi representing different sub-populations were examined using 7 highly polymorphic microsatellite markers (). For genotyping, tissue was lysated in 100 μL low TE buffer (10 mM Tris, 0.1 mM EDTA) for 1 h at 56°C, while shaking and 1 μl extract used as template for the PCR reaction. Primer concentrations and PCR conditions followed published protocols (; ; ). Amplified fragments were analyzed on an ABI 3130 genetic analyzer using Rox-350 (Applied Biosystems) as internal size standard. Allele sizes were scored using the software GENEMARKER v1.91 (SoftGenetics, LLC). Pairwise FST values were computed in ARLEQUIN v3.5.1.2 () using non-parametric permutation procedures. All tests were performed with an initial alpha level of 0.05 and corrected following the Benjamini-Hochberg False Discovery Rate (B-H FDR) procedure ().

Results

Taxonomic Composition of Ctenophore Associated Bacteria Shows Difference to Ambient Seawater

The microbiota associated with different epithelia (epi- and gastrodermis) and sub-populations of M. leidyi as well as the bacterial communities of ambient water were analyzed by Illumina sequencing the hypervariable V1-V2 regions of the bacterial 16S rRNA gene. Sequence depths of 117 analyzed samples ranged between 23,622 to 675,908 single sequences. In total, 2,829 OTUs at a 97% similarity threshold (roughly corresponding to species level) were considered. Abundances of OTUs (97% sequence similarity) between individual samples are depicted at genus level in Figure 1. Consequently, main representative bacterial classes in the samples were Alphaproteobacteria (e.g., including genera Nautella, uncl. Pelagibacteraceae, uncl. Rhodobacteraceae, Octadecabacter), Gammaproteobacteria (e.g., including genera Marinomonas, Halomonas, Shewanella), Mollicutes (e.g., including genera uncl. Entoplasmatales, Mycoplasma), Cyanobacteria (e.g., including genus Synechococcus) and Bacteroidia (e.g., including genus Bacteroides). While OTU abundances between the different sample types were significant, in the following the different comparisons are considered in more detail.

FIGURE 1

Significant Differences Between Associated Bacterial Communities of M. leidyi and Ambient Seawater

Beta diversity analysis demonstrated that the associated bacterial communities of the ctenophore sub-populations were significantly different from the community compositions present in the corresponding ambient water (Figure 2 and Supplementary Table S1). Water samples were mainly dominated by uncultured Pelagibacteraceae and uncultured Rhodobacteraceae of class Alphaproteobacteria, and Synechococcus of class Cyanobacteria, whereas the ctenophore tissues were mostly colonized by Halomonas, Shewanella, and Marinomonas of Gammaproteobacteria, and Mollicutes.

FIGURE 2

Tissue Variation of Ctenophore-Associated Bacteria

Sequencing results based on the 16S rRNA gene amplicon sequencing revealed differences in the bacterial composition among two different tissues, the outer mucus layer and epidermis, and the epithelium in the gastric cavity (Figures 1, 2 and Supplementary Table S1). For all sub-populations, the tissues harbored a significantly different and specific microbiota (P = 0.001), except for the North Sea, where the microbiota composition differed only marginally (P = 0.059). Epidermis from native north-east animals (Atlantic, Massachusetts, United States) was dominated by Marinomonas and unclassified Pelagibacteraceae, whereas their gastric cavity was mainly colonized by Shewanella, Mycoplasma, and Bacteroides, unclassified Erythrobacteraceae and a multitude of low abundant OTUs. Both, epidermis and gastric cavity from native south-east animals (Atlantic, Florida, United States) were dominated by Nautella; epidermis additionally comprised unclassified Entoplasmatales, whereas gastric cavity was additionally colonized by unclassified Enterobacteriaceae, Bacteroides and Arcobacter. In comparison, epidermis and gastric samples of the invasive Mediterranean animals mostly comprised Synechococcus, while epidermis samples additionally contained unclassified Entoplasmatales and Candidatus Aquiluna, whereas gastric samples comprised Halomonas, Shewanella, and Mycoplasma, unclassified Erythrobacteraceae and Bacteroides. Main representatives in epidermis as well as gastric samples of invasive North Sea animals were Halomonas, Shewanella, and Mycoplasma, unclassified Enterobacteriaceae and Bacteroides. Epidermis and gastric cavity community patterns of those animals were almost similar, however, slight differences were observed in abundances of main representatives and presence of different low abundant taxa.

Sub-Population Specific Community Patterns of M. leidyi

Population structure analyses using the polymorphism of seven microsatellite loci revealed that all tested sub-populations were significantly differentiated (FST range 0.03 to 0.33; P < 0.05 after correction for multiple testing, Table 1). The largest observed population differentiation was found between the southern native sub-population (Atlantic, Florida, United States) and the northern invasive sub-population (Helgoland, North Sea) with FST = 0.333 (P < 0.0001, see Table 1), while the smallest population difference was observed between native northern sub-populations (Atlantic, Massachusetts, United States) in comparison with their recently introduced counterparts in Northern Europe (Helgoland, North Sea) with FST = 0.032 (P < 0.05, Table 1). Based on the revealed genotyping of M. leidyi, bacterial community analysis of native northern (Atlantic, Massachusetts, United States) and southern (Atlantic, Florida, United States) as well as invasive (Mediterranean Sea and North Sea) M. leidyi sub-populations, bacterial community patterns where analyzed concerning ctenophore provenance (Figure 3 and Supplementary Table S1). RDA plots demonstrated that every single sub-population harbored a specific microbiota. Native sub-populations are significantly different and share just few OTUs. Invasive sub-populations shared a more similar microbiota and additionally were more similar to native northern sub-population animals (Atlantic, Massachusetts, United States), however, significantly different to native southern sub-population animals (Atlantic, Florida, United States) (Figure 3).

TABLE 1

MSFLWHHE
MS
FL0.078∗∗∗
WH0.277∗∗∗0.157∗∗∗
HE0.333∗∗∗0.205∗∗∗0.032

Pairwise comparison of genetic differentiation (FST) between native and invasive Mnemiopsis leidyi sub-populations using highly polymorphic microsatellite markers.

The four tested sub-populations are abbreviated with southern native, FL (Florida, United States, south-east Atlantic coast); northern native, WH (Woods Hole, Massachusetts, United States, north-east Atlantic coast); southern invasive, MS (Marseille, France, north-east Mediterranean Sea); and northern invasive, HE (Helgoland, German Bight, North Sea). Significant comparisons are highlighted in bold with P-levels abbreviated by P < 0.05 and ∗∗∗P < 0.0001 (P-values corrected for multiple testing).

FIGURE 3

Discussion

This is the first study comparing the microbiota of native and invasive sub-populations of the comb jelly M. leidyi using standardized high throughput amplicon sequencing methods based on the V1-V2 hypervariable region of the 16S rRNA gene. We demonstrated that the microbiota was significantly different from the surrounding seawater and specific to each analyzed sub-population. Contrary to the previous report of describing that the outer surface of M. leidyi is free from associated microbes (), we observed a highly diverse microbiota on epithelial tissues. Additionally, and in line with , we detected different community patterns on the outer surface of the animals compared to the gastric cavity. Two published studies (; ) dealing with microbiota patterns of invasive M. leidyi observed that the core microbiome of M. leidyi is host-specific and main representatives were identified belonging to Gamma- and Alphaproteobacteria as well as Bacteroidetes with Marinomonas spp. (Gammaproteobacteria) as one abundant and most likely important colonizer. However, further work is needed to assess the functional roles of such ctenophore-associated microbes.

Our multivariate analyses indicated that native and invasive northern sub-populations are more similar compared with southern native and invasive sub-populations. For northern native sub-populations, Marinomonas spp. dominated the bacterial community. However, for invasive northern sub-populations, this species was only observed with a limited contribution to the microbiome of M. leidyi sampled during summer 2012. This is in contrast to previous findings, where during the years 2009/2010 Marinomonas spp. had been documented to be the most abundant bacterial species on M. leidyi from the invasive northern sub-population, i.e., Helgoland Road () and the Gulmar Fjord (). A change in the bacterial community composition with reduced contribution of Marinomonas spp. could be related to a significant change in winter conditions as observed during the early 2010s, where a series of cold winters had led to a retraction of the range expansion of M. leidyi in northern Europe with the extinction of M. leidyi in, e.g., the Gulmar Fjord (). Hence, drastic physical conditions and extinction of local sub-populations could be reflected in a change of the microbiome, but deserve further investigation. It has been shown that the microbiome changes with season in southern native and northern invasive M. leidyi sub-populations (; ). In this study, we do not expect that seasonal differences impact the general conclusions of our study as all regions were sampled during the same seasonal time point.

In addition to the comparison of the microbiota on the comb jelly with the ambient water as well as between outer surface and gastric cavity, to the best of our knowledge, this is the first study focusing on the comparison between native and invasive M. leidyi. Therefore, our results are an essential first step to assess the impact of microbiota on the invasion process by describing basic diversity pattern.

A recent study investigating the microbiota of M. leidyi in comparison with other ctenophore species from the same environment documented significant difference between the microbiome of all studied ctenophore species (). All investigated M. leidyi microbiomes were dominated by Proteobacteria, though variable through time (). Our deep sequencing approach revealed that the native southern sub-population showed a large difference of the microbiota when compared to the invasive southern sub-population. Specifically, the native southern microbiome of M. leidyi was dominated by Nautella (Alphaproteobacteria), while this genus was absent from the corresponding invasive M. leidyi sub-population. In line with this, southern invasive and native M. leidyi sub-populations have also been differentiated for the longest time (), which might explain the difference in their respective microbial community. Similar to our results, in a previous study, Alphaproteobacteria had been shown to be key bacterial community members of M. leidyi from the native southern range in the Gulf of Mexico during late spring and summer (). For the Mediterranean Sea, no previous microbiome characterization has been conducted. In our sequence data, we observed a large proportion of the Cyanobacteria species Synechococcus in epidermis as well as gastric samples. Even though this is a free living bacterium with a widespread distribution in marine and freshwater systems, Cyanobacteria have also been shown to contribute to the core microbiome of native southern M. leidyi during certain periods () and might be more abundant during summer when they bloom in ocean surface waters ().

showed significant differences in the microbiota associated to the polyp stage of the moon jellyfish Aurelia aurita among three source populations that were cultured for more than 10 years under similar artificial lab conditions. In the present study, we found the lowest host population differentiation between native north-east United States sub-population (Atlantic, Massachusetts, United States) and the invasive northern sub-population (Helgoland, North Sea). These host population structure results are in agreement with molecular results (). In detail, a recent investigation summarizing the range expansion of M. leidyi in western Eurasia has shown that the northern introduction of M. leidyi is very recent and dates back to 2005 (). Hence, genetic background along with time since first introduction has an impact on microbiome evolution. We further showed that analyzed invasive sub-populations, in general, share a more similar microbiota, most likely due to a more similar environment they share, disclosing the impact of the environment. Moreover, we observed that the microbiota of those invasive M. leidyi sub-populations is more similar to the northern native sub-population from Massachusetts than to the southern native sub-population from Florida. We hypothesize that the observed microbiota reflects the invasion process of M. leidyi, since the southern invasive sub-population has been introduced from the southern native sub-population in the 1980’s, while the northern sub-population present in European waters has been introduced from the north-east United States coast during the early 2000’s as reviewed in . Further common garden experiments have to be performed in the future to confirm this hypothesis.

In conclusion, significant variance in microbial community patterns and significant population differences between sub-populations suggest an impact of environment along with the genetic host background (host population affiliation) on the microbiota. Long-term experiments over several generations under controlled conditions are necessary to partition the influence of the environment from population affiliation and genotype of the animals. We found that M. leidyi harbors a diverse microbiota on its epithelia. The significant difference between water and the microbiota suggests that ctenophores might select the microbes in a host-specific way, thereby shaping the entire metaorganism. The direct interactions between microbes and host need to be studied in controlled environments to gain insights into the impact of the associated microbiota on the fitness and invasion success of the ctenophore M. leidyi.

Statements

Data availability statement

Sequence data were submitted to the NCBI Sequence Read Archive under Accession No. PRJNA548849 and are additionally stored in a sequence data management system at the Institute of Clinical Molecular Biology (IKMB) at Kiel University, Germany.

Author contributions

CJ, NW-B, RS, and TR: conceptualization. CJ, NW-B, and SK: sample analyses. MF with support from NW-B: bioinformatic analyses. CJ with support from TR: genotyping and analyses. CJ and NW-B with support from MF, SK, RS, and TR: manuscript drafting. All authors read and approved the final version of the manuscript draft.

Funding

This study received financial support from the Collaborative Research Center (CRC) 1182 “Origin and Function of Metaorganisms” to TR and RS, funded through the German Research Foundation (DFG, Re1108/17-1).

Acknowledgments

We thank Sören Bolte for support and data collection, Katharina Bading for collection of animals in the native southern range, and Nicole Pinnow for experimental support. This manuscript has been released as a Pre-Print at https://www.biorxiv.org/content/10.1101/601419v1.

Conflict of interest

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2019.00635/full#supplementary-material

References

Summary

Keywords

ctenophore, metaorganism, bacteria, species translocations, invasive species

Citation

Jaspers C, Weiland-Bräuer N, Fischer MA, Künzel S, Schmitz RA and Reusch TBH (2019) Microbiota Differences of the Comb Jelly Mnemiopsis leidyi in Native and Invasive Sub-Populations. Front. Mar. Sci. 6:635. doi: 10.3389/fmars.2019.00635

Received

03 June 2019

Accepted

27 September 2019

Published

11 October 2019

Volume

6 - 2019

Edited by

Sandie M. Degnan, The University of Queensland, Australia

Reviewed by

Adam Michael Reitzel, The University of North Carolina at Charlotte, United States; Cynthia B. Silveira, San Diego State University, United States

Updates

Copyright

*Correspondence: Cornelia Jaspers, ; Nancy Weiland-Bräuer,

These authors have contributed equally to this work

This article was submitted to Marine Molecular Biology and Ecology, a section of the journal Frontiers in Marine Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics