- Short Communication
A comparison of microeukaryote communities inhabiting sponges and seawater in a Taiwanese coral reef system
Annals of Microbiology volume 69, pages 861–866 (2019)
Assess microeukaryote community composition in seawater and sponge samples from Taiwanese coral reefs.
In the present study, we used Illumina sequencing to explore the microeukaryote communities of seven biotopes (six sponge species and seawater) sampled in the Penghu archipelago of Taiwan.
Microeukaryote communities were dominated by Dinoflagellates with Dinophyceae and Syndiniales well represented in all biotopes. Other abundant taxa included metazoa, red and green algae and Radiolaria. The only significant differences were a significantly higher relative abundance of Picobiliphyta and Stramenopiles_X in seawater and Metamonada in the sponge Acanthostylotella cornuta. There was also a significant difference in composition among biotopes with samples from sponges and seawater forming distinct clusters. There was, however, no congruence between prokaryote and microeukaryote community composition. After removing all OTUs < 100 sequences, more than 90% of remaining OTUs representing > 99.5% of sequences were shared between sponge and seawater samples.
This data in the present study would appear to suggest that marine microeukaryote communities in sponges are largely derived from the surrounding seawater. Abundant OTUs were also related to organisms previously retrieved from seawater. A number of these OTUs though had relatively low sequence similarity to organisms in GenBank suggesting that more research of the microeukaryote communities in the Penghu archipelago may yield novel organisms in this relatively unexplored area.
Sponges are sessile, benthic and mainly filter-feeding organisms (Diaz and Rützler 2001) that can be found in several aquatic habitats (freshwater, marine, warm or cold water). They are involved in the bioerosion and consolidation of the coral reef carbonate framework, modification of the water column due to filtration and secondary metabolite emanation and benthic-pelagic coupling (Bell 2008; Maldonado et al. 2016). Sponges are also noteworthy for the range and abundance of natural products (Hochmuth et al. 2010) and their prokaryote diversity and abundance (Taylor et al. 2007), which are widely accepted to be inter-connected (Piel 2009). In addition to being involved in the biosynthesis of secondary metabolites (Fan et al. 2012; Hentschel et al. 2012), sponge microbial communities also play a critical role in the sponge metabolism (Hentschel et al. 2006; Taylor et al. 2007). Prokaryotes are not the only inhabitants of the sponge holobiont. Another important, but largely overlooked group is the sponge microeukaryote community. Relatively little research has focused on this group compared to the relatively well-known prokaryote community (Rodríguez-Marconi et al. 2015).
In the present study, we used Illumina sequencing to explore and compare the microeukaryote communities of 6 sponge species and seawater collected in the Penghu archipelago, Taiwan from 25 to 29 of July, 2014 (Supplementary Table 1; see Huang et al. (2016) and Coelho et al. (2018) for detailed descriptions of the sampling area). The sponge species sampled were Acanthostylotella cornuta, Agelas cavernosa, Echinodictyum asperum, Hyrtios erectus, Stylissa carteri and Suberites diversicolor. Two to three replicates were sampled per species. Sponges were photographed in situ, collected using scuba diving, brought back to the laboratory and preserved in 95% ethanol for further identification and molecular work. All specimens have been deposited at the Naturalis Biodiversity Center, the Netherlands. Our main objectives were to (1) test for compositional differences among samples from different biotopes (sponge species and seawater) and (2) test for significant compositional congruence between prokaryote and microeukaryote communities.
Prokaryotes and microeukaryotes were sequenced, but the present study will mainly focus on microeukaryotes and only use the prokaryote data to test for compositional congruence between both groups. Briefly, for prokaryotes, the 16S rRNA gene V3V4 variable region was amplified with PCR primers 341F 5′-CCTACGGGNGGCWGCAG-3′ and 785R 3′-GACTACHVGGGTATCTAATCC-5′ and with barcode on the forward primer. For microeukaryotes, 18S rDNA gene fragments were amplified using primers TAReuk454FWD1 (5′-CCAGCA(G/C)C(C/T)GCGGTAATTCC-3′ and TAReukREV3 (5′-ACTTTCGTTCTTGAT(C/T)(A/G)A-3′ (Stoeck et al. 2010). DNA extraction and sequencing using an Illumina MiSeq device followed Coelho et al. (2018) and Swierts et al. (2018). Resultant files were analysed using QIIME (Quantitative Insights Into Microbial Ecology; (Caporaso et al. 2010) (http://www.qiime.org/; last checked 2017-01-20) and UPARSE (Edgar 2013) with OTU selection at 97% similarity cut-off. Taxonomy for microeukaryotes was assigned using the PR2 database (http://ssu-rrna.org/pr2) (Guillou et al. 2013) and for prokaryotes using the SILVA_128_QIIME_release database (Quast et al. 2013). All statistical analyses were performed in the R environment (R Core Team 2013). For a more detailed description of the sequencing and statistical analyses, see Coelho et al. (2018), Cleary et al. (2018) and Cleary and Polónia (2018). The sequences generated in this study can be downloaded from the NCBI SRA: SRP109605 (SUB2790757; SUB4184084).
Sequencing yielded 1093775 sequences binned into 2882 microeukaryote OTUs after quality control. OTU richness varied widely among and within biotopes with the most prominent feature, a higher richness for samples collected in the southern Penghu islands including a much higher richness for the water sample collected in the southern islands versus both water samples collected in the northern islands (Supplementary Table 1). The most abundant taxa overall was the phylum Dinophyta (dinoflagellates, 545564 sequences; 1593 OTUs) followed by the kingdom Metazoa (302301 sequences; 286 OTUs), phylum/division Chlorophyta (green algae, 89355 sequences; 96 OTUs), phylum Rhodophyta (red algae, 88745 sequences; 244 OTUs), the subphylum Radiolaria (phylum Retaria, 24724 sequences; 109 OTUs) and the phylum Apicomplexa (16205 sequences; 71 OTUs) (Supplementary Fig. 1).
Relative abundance varied considerably among and within biotopes, but this variation was only significant for Picobiliphyta (Fig. 1h), Stramenopiles_X (Fig. 1j) and Metamonada (Fig. 1u) after applying a Bonferroni correction for multiple comparisons. The relative abundance of OTUs assigned to the Picobiliphyta and Stramenopiles_X was significantly higher in seawater than other biotopes, whereas the relative abundance of OTUs assigned to the Metamonada was significantly higher in A. cornuta than all other biotopes (Supplementary Table 2). Microeukaryote OTU richness was highest in E. asperum and lowest in S. diversicolor (Fig. 1x). Evenness, in turn, was highest in seawater and lowest in A. cornuta (Fig. 1w).
There was a highly significant difference in microeukaryote composition among biotopes (Adonis, F6,11 = 1.73; P < 0.001; R2 = 0.486; Fig. 2a). The first axis separated samples of sponges and water collected from the northern Penghu islands at positive axis 1 values from samples of sponges and water collected in the southern Penghu islands at negative axis 1 values (Supplementary Table 1). The only exception was a single specimen of H. erectus collected in Si Ji Yu (southern Penghu). The second axis separated samples of water from samples of sponges. There was no significant congruence between prokaryote (Fig. 2b) and microeukaryote (Fig. 2a) composition (Fig. 2c; Procrustes correlation = 0.386, P = 0.130).
A number of OTUs significantly discriminated (Simper analysis; P < 0.01) between pairs of biotopes (Supplementary Fig. 2). By far, the most abundant OTU was OTU-3, which significantly discriminated between A. cornuta and Suberites sp. (Simper analysis, P < 0.001; Supplementary Table 3). OTU-3 was assigned to the coral genus Montipora (Supplementary Table 4). OTUs 36 and 306 significantly discriminated between seawater and all sponge species (Simper analysis, P < 0.05). OTUs 36 and 306 were assigned to the Dinophyceae and closely related to organisms obtained from seawater (Supplementary Table 4). OTU-59, assigned to the Dinophyceae, significantly discriminated between the sponge species A. cornuta and all other biotopes (Simper analysis; P < 0.05; Supplementary Table 3). OTU-8, assigned to the Syndiniales, significantly discriminated between H. erectus and the sponge species A. cavernosa and S. diversicolor in addition to seawater. The highly variable nature of OTUs within biotopes can be seen in Supplementary Figs. 2 and 3.
The microeukaryote OTU subset of all OTUs < 100 sequences (OTUs100) included 390 OTUs100 and 1059147 sequences (96.8% of all sequences). Of the 390 OTUs100, 363 (93.1%) representing 1053869 sequences (99.5%) were found in seawater samples and 27 (278 sequences; 0.5% of total sequences) were only found in sponge samples. After controlling for multiple tests, there was no significant difference in the relative abundance of seawater OTUs among sponge species (Fig. 1v).
Dinoflagellates were by far the most abundant eukaryote taxon in terms of relative abundance and number of OTUs. Both of the major taxa within the dinoflagellates, the Dinophyceae and Syndiniales were well represented in all biotopes. Dinoflagellates are dominant unicellular members of marine and freshwater plankton with high morphological and trophic diversity that play critical roles as predators, primary producers, parasites and symbionts. They are also responsible for toxic algal blooms with deleterious effects to aquatic and human life (Gómez 2012; Price and Bhattacharya 2017). Only one abundant OTU, OTU-28, was assigned to the genus Symbiodinium, and was highly similar to an organism identified as Symbiodinium goreaui strain CCATM-210 obtained from Iranian seawater (Supplementary Table 4). Although well known as symbionts of scleractinian hosts, species of Symbiodinium inhabit a wide range of metazoan hosts (including foraminifera, black corals and sponges) where they perform important functional tasks including transmembrane transport and translocation of photosynthesis products to their host; they have also been shown to be enriched for functions including reactive oxygen species response and UV radiation protection (González-Pech et al. 2017; Ramsby et al. 2017). Although stable associations with species are known from just a few sponge species, Strehlow et al. (2016) showed that normally non-symbiont-bearing sponge species were able to acquire and maintain Symbiodinium for a number of days after experimental inoculation. Wecker et al. (2015) also recorded dinoflagellates including Symbiodinium from a number of nudibranchs and sponges in French Polynesia. Symbiodinium Clade C was found in the sponge species Lamellodysidea herbacea while other Dinophyceae members were recorded in the sponge species Haliclona sp., Phycopsis sp. and Flabellina sp. Dinoflagellates have also been recorded in freshwater sponges inhabiting lake Baikal including sequences assigned to the order Suessiales, known for establishing symbiotic relationships with other invertebrates (Annenkova et al. 2011). Other numerically important taxa included metazoa, red and green algae and Radiolaria.
Although there was compositional variation in eukaryote composition among biotopes, the degree of differentiation was less pronounced than that observed for prokaryote communities. In contrast to prokaryote communities, there was a pronounced effect of sampling location (northern versus southern islands). This would appear to suggest that microeukaryote communities in sponges are more susceptible to local environmental conditions than prokaryote communities. A more detailed and well-designed sampling strategy is, however, needed to explicitly test this hypothesis.
Our results are also in line with Chaib De Mares et al. (2017) who also found no significant congruence between bacterial and microeukaryote composition and noted that the compositional structure of dinophytes, ciliophores (alveolates) and stramenopiles could not be explained by host sponge or geographic location at large geographic scales (Caribbean versus Mediterranean). In the present study, more than 90% of microeukaryote OTUs100 representing more than 99.5% of all sequences were also recorded in the seawater samples. This would appear to suggest that microeukaryote communities in sponges are largely derived from seawater and may just represent environmental DNA. That being said, microeukaryote communities are far less known than prokaryote communities and the relatively low sequence similarities of abundant OTUs sampled in this study to organisms in GenBank suggest that novel taxa remain to be discovered.
Annenkova NV, Lavrov DV, Belikov SI (2011) Dinoflagellates associated with freshwater sponges from the ancient Lake Baikal. Protist 162:222–236. https://doi.org/10.1016/j.protis.2010.07.002
Bell J (2008) The functional roles of marine sponges. Estuar Coast Shelf Sci 79:341–353. https://doi.org/10.1016/j.ecss.2008.05.002
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK et al (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336
Chaib De Mares M, Sipkema D, Huang S, Bunk B, Overmann J and van Elsas JD (2017) Host specificity for bacterial, archaeal and fungal communities determined for high- and low-microbial abundance sponge species in two genera. Front. Microbiol 8:2560. https://doi.org/10.3389/fmicb.2017.02560
Cleary DFR, Polónia ARM (2018) Bacterial and archaeal communities inhabiting mussels, sediment and water in Indonesian anchialine lakes. Antonie Van Leeuwenhoek 111:237–257. https://doi.org/10.1007/s10482-017-0944-1
Cleary DFR, Polónia ARM, de Voogd NJ (2018) Bacterial communities inhabiting the sponge Biemna fortis, sediment and water in marine lakes and the open sea. Microb Ecol 76:610–624. https://doi.org/10.1007/s00248-018-1156-6
Coelho FJRC, Cleary DFR, Gomes NCM, Pólonia ARM, Huang YM, Liu LL, de Voogd NJ (2018) Sponge prokaryote communities in Taiwanese coral reef and shallow hydrothermal vent ecosystems. Microb Ecol 75:239–254. https://doi.org/10.1007/s00248-017-1023-x
Diaz MC, Rützler K (2001) Sponges: an essential component of Caribbean coral reefs. Bull Mar Sci 69:535–546
Edgar RC (2013) UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods 10:996–998
Fan L, Reynolds D, Liu M, Stark M, Kjelleberg S, Webster NS, Thomas T (2012) Functional equivalence and evolutionary convergence in complex communities of microbial sponge symbionts. Proc Natl Acad Sci U S A 109:E1878–E1887. https://doi.org/10.1073/pnas.1203287109
Gómez F (2012) A quantitative review of the lifestyle, habitat and trophic diversity of dinoflagellates (Dinoflagellata, Alveolata). Syst Biodivers 10:267–275
González-Pech RA, Ragan MA, Chan CX (2017) Signatures of adaptation and symbiosis in genomes and transcriptomes of Symbiodinium. Sci Rep 7:15021. https://doi.org/10.1038/s41598-017-15029-w
Guillou L, Bachar D, Audic S et al (2013) The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy. Nucleic Acids Res 41(Database issue):D597–D604. https://doi.org/10.1093/nar/gks1160
Hentschel U, Usher KM, Taylor MW (2006) Marine sponges as microbial fermenters. FEMS Microbiol Ecol 55:167–177
Hentschel U, Piel J, Degnan SM, Taylor MW (2012) Genomic insights into the marine sponge microbiome. Nat Rev Microbiol 10:641–654
Hochmuth T, Niederkrüger H, Gernert C, Siegl A, Taudien S, Platzer M, Crews P, Hentschel U, Piel J (2010) Linking chemical and microbial diversity in marine sponges: possible role for poribacteria as producers of methyl-branched fatty acids. Chembiochem 11:2572–2578. https://doi.org/10.1002/cbic.201000510
Huang YM, de Voogd NJ, Cleary DFR, Li T-H, Mok HK, Ueng JP (2016) Biodiversity pattern of subtidal sponges (Porifera: Demospongiae) in the Penghu Archipelago (Pescadores), Taiwan. J Mar Biol Assoc UK 96:417–427
Maldonado M, Aguilar R, Bannister RJ, Bell D, Conway KW, Dayton PK, Díaz C, Gutt J, Kelly M et al (2016) Sponge grounds as key marine habitats: a synthetic review of types, structure, functional roles, and conservation concerns. Marine Animal Forests. Springer, Berlin. https://doi.org/10.1007/978-3-319-17001-5_24-1
Piel J (2009) Metabolites from symbiotic bacteria. Nat Prod Rep 26:338–362
Price DC, Bhattacharya D (2017) Robust Dinoflagellata phylogeny inferred from public transcriptome databases. J Phycol 53:725–729
Quast C et al (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41(D1):D590–D596
R Core Team (2013) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna 3-900051-07-0. http://www.R-project.org
Ramsby BD, Hill MS, Thornhill DJ, Steenhuizen SF, Achlatis M, Lewis AM, LaJeunesse TC (2017) Sibling species of mutualistic Symbiodinium clade G from bioeroding sponges in the western Pacific and western Atlantic oceans. J Phycol 53:951–960. https://doi.org/10.1111/jpy.12576
Rodríguez-Marconi S, De la Iglesia R, Díez B, Fonseca CA, Hajdu E, Trefault N (2015) Characterization of bacterial, archaeal and eukaryote symbionts from antarctic sponges reveals a high diversity at a three-domain level and a particular signature for this ecosystem. PLoS One 10:e0138837. https://doi.org/10.1371/journal.pone.0138837
Stoeck T, Bass D, Nebel M, Christen R, Jones MD, Breiner HW, Richards TA (2010) Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol Ecol 19:21–31. https://doi.org/10.1111/j.1365-294X.2009.04480.x
Strehlow B, Friday S, McCauley M, Hill M (2016) The potential of azooxanthellate poriferan hosts to assess the fundamental and realized Symbiodinium niche: evaluating a novel method to initiate Symbiodinium associations. Coral Reefs 35:1201–1212. https://doi.org/10.1007/s00338-016-1465-5
Swierts T, Cleary DFR, de Voogd NJ (2018) Biogeography of prokaryote communities in closely related giant barrel sponges across the Indo-Pacific. FEMS Microbiol Ecol 94(12):fiy194. https://doi.org/10.1093/femsec/fiy194
Taylor MW, Radax R, Steger D, Wagner M (2007) Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol Mol Biol Rev 71:295–347
Wecker P, Fournier A, Bosserelle P, Debitus C, Lecellier G, Berteaux-Lecellier V (2015) Dinoflagellate diversity among nudibranchs and sponges from French Polynesia: insights into associations and transfer. C R Biol 338:278–283. https://doi.org/10.1016/j.crvi.2015.01.005
Support in the field and lab was provided by Julian Cleary, Floris Cleary, Yusheng Huang, Kate, Ana R.M. Polónia and Nicole J de Voogd.
Financial support was provided to CESAM (UID/AMB/50017 – POCI-01-0145-FEDER-007638) and for the project LESS CORAL (PTDC/AAC-AMB/115304/2009) by FCT/MEC through national funds, and co-funding by FEDER, within the PT2020 Partnership Agreement and Compete 2020.
Conflict of interest
The author declares that he has no conflict of interest.
Research involving human participants and/or animals
This article does not contain any studies with human participants or laboratory animals.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Relative abundance of the most abundant microeukaryote higher taxa in sponges: Ac - Agelas cavernosa, Ar - Acanthostylotella cornuta, He - Hyrtios erectus, Xt - Xestospongia testudinaria, Echinodictyum asperum, Sc - Stylissa carteri, Su - Suberites diversicolor and Wt - water. (PDF 6.94 kb)
Relative abundance of significantly discriminating OTUs between pairs of biotopes identified using Simper analysis (P < 0.001) and colour-coded according to microeukaryote taxon for Ac - Agelas cavernosa, Ar - Acanthostylotella cornuta, He - Hyrtios erectus, Xt - Xestospongia testudinaria, Ea - Echinodictyum asperum, Sc - Stylissa carteri, Su - Suberites diversicolor and Wt – water. The circle size of the OTU is proportional to the mean percentage of sequences per sample as indicated by the symbol legend in the bottom right corner. (PDF 66.6 kb)
Heatmap showing the abundance of abundant microeukaryote OTUs (≥ 4000 sequences). The heatmap was generated using the function heatmap2() in the R package gplots (http://www.cran.r-project.org/). The OTUs were log-transformed and clustered according to their occurrence by UPGMA hierarchical clustering. Sponge species: Ac - Agelas cavernosa, Ar - Acanthostylotella cornuta, He - Hyrtios erectus, Xt - Xestospongia testudinaria, Ea - Echinodictyum asperum, Sc - Stylissa carteri, Su - Suberites diversicolor and Wt - water. (PDF 11.9 kb)
A summary of the samples collected is presented study including the name of the site, the location of the site in the northern or southern Penghu islands, the GPS coordinates, the biotope (host species, sediment or seawater), order and family of the host species. (XLS 20.5 kb)
Results of emmeans analysis showing pairwise comparisons of differences in the relative abundances of selected eukaryote higher taxa and the percentage of OTUs100 recorded in seawater as a percentage of total OTUs100 (WtToPr) between biotopes based on the ‘fdr’ test. Estimate (Estimated marginal means), SE (standard error), P (probability), Sig. (Significance): * 0.01 < Pr < 0.05 ** 0.001 < Pr < 0.01; *** Pr < 0.001. Variable: the dependent variable, Contrast: Contrasts between pairs of biotopes, SE: standard error. (XLS 78.5 kb)
Results of Simper analysis showing the contribution of microeukaryote OTUs to differences in similarity between pairs of biotopes. Contrast: contrasts between pairs of biotopes. Average: average contribution to overall dissimilarity. Sd: Standard deviation of contribution. Ratio: Average to sd ratio. Ava, Avb: average abundances per biotope. CumSum: ordered cumulative contribution. P: permutation p value. OTUs that contributed significantly to differences are indicated by significance (Sig.): * 0.01 < P < 0.05 ** 0.001 < P < 0.01; *** P < 0.001. (XLS 340 kb)
List of abundant (≥ 4000 sequence reads) OTUs and closely related organisms identified using BLAST search. OTU: OTU number; Sum: number of sequence reads; Group: biotope or biotopes where the OTUs were mainly found; Acc: Genbank accession numbers of closely related organisms identified using BLAST; Seq: sequence similarity of these organisms with our representative OTU sequences; Source: isolation source of organisms identified using BLAST. In the Group column the abbreviations stand for sponge species: Ac - Agelas cavernosa, Ar - Acanthostylotella cornuta, He - Hyrtios erectus, Xt - Xestospongia testudinaria, Ea - Echinodictyum asperum, Sc - Stylissa carteri, Su - Suberites diversicolor and Wt – seawater. (XLS 13.5 kb)
About this article
Cite this article
Cleary, D.F.R. A comparison of microeukaryote communities inhabiting sponges and seawater in a Taiwanese coral reef system. Ann Microbiol 69, 861–866 (2019). https://doi.org/10.1007/s13213-019-01476-5
- Coral reefs
- Penghu islands