Sponge-Specific Unknown Bacterial Groups Detected in Marine Sponges Collected from Korea Through Barcoded Pyrosequencing
Sponge-Specific Unknown Bacterial Groups Detected in Marine Sponges Collected from Korea Through Barcoded Pyrosequencing
Journal of Microbiology and Biotechnology. 2015. Jan, 25(1): 1-10
Copyright © 2015, The Korean Society For Microbiology And Biotechnology
  • Received : June 13, 2014
  • Accepted : October 13, 2014
  • Published : January 28, 2015
Export by style
Cited by
About the Authors
Jong-Bin, Jeong
Department of Biological Science and Biotechnology, Hannam University, Daejeon 305-811, Republic of Korea
Kyoung-Ho, Kim
Department of Microbiology, Pukyong National University, Busan 608-737, Republic of Korea
Jin-Sook, Park
Department of Biological Science and Biotechnology, Hannam University, Daejeon 305-811, Republic of Korea

The bacterial diversity of 10 marine sponges belonging to the species Cliona celata , an unidentified Cliona species, Haliclona cinerea , Halichondria okadai , Hymeniacidon sinapium , Lissodendoryx isodictyalis , Penares incrustans , Spirastrella abata , and Spirastrella panis collected from Jeju Island and Chuja Island was investigated using amplicon pyrosequencing of the 16S rRNA genes. The microbial diversity of these sponges has as of yet rarely or never been investigated. All sponges, except Cliona celata , Lissodendoryx isodictyalis , and Penares incrustans , showed simple bacterial diversity, in which one or two bacterial OTUs occupied more than 50% of the pyrosequencing reads and their OTU rank abundance curves saturated quickly. Most of the predominant OTUs belonged to Alpha -, Beta -, or Gammaproteobacteria . Some of the OTUs from the sponges with low diversity were distantly (88%~89%) or moderately (93%~97%) related to known sequences in the GenBank nucleotide database. Phylogenetic analysis showed that many of the representative sequences of the OTUs were related to the sequences originating from sponges and corals, and formed sponge-specific or -related clades. The marine sponges investigated herein harbored unexplored bacterial diversity, and further studies should be done to understand the microbes present in sponges.
Sponges are well known as reservoirs of microbes; they are inhabited by high numbers of bacteria, which were estimated to be up to 10 9 cells per milliliter of sponge tissue [14] , accounting for as much as 50% of the biomass of the sponges [31] . The microbes in sponges have various roles, including as food sources [29] , pathogenic agents [2] , symbiotic residents [49] , sponge skeleton stabilizers, metabolic waste consumers, secondary metabolite producers [11] , and nutrient cycle mediators [40] .
Numerous studies have been performed to reveal the diversity in sponges. Culture-dependent studies based on the isolation of pure cultures revealed various kinds of bacteria, [44 , 48] and novel strains were isolated and proposed as new taxa [28 , 51] . Many studies based on culture-independent methods have also been performed to reveal so-called unidentified majorities [1] and unknown bacterial diversity in sponges [13 , 36] . The massive parallel 454 pyrosequencing, a next-generation sequencing method, has been applied to investigate the bacterial diversity in various environments [7 , 30 , 34] as well as in marine sponges [17 , 18 , 24 , 33 , 34 , 45 , 46] .
Marine sponges harbor various kinds of bacteria, comprising a total of 25 bacterial phyla [34] . The most diverse phyla contained in sponges are Proteobacteria , Chloroflexi , and Poribacteria , but many less diverse phyla such as Acidobacteria , Actinobacteria , Cyanobacteria , Gemmatimonadetes , Bacteroidetes , Spirochaetes , Firmicutes , Nitrospirae , TM7, SBR1093, and OS-K have also been observed [34] .
However, some sponges, known as low microbial abundance (LMA) sponges, have low microbial content of only about 10 5 -10 6 cells per gram or milliliter [12] . LMA sponges also have lower taxonomic diversity (less number of taxa) than high microbial abundance (HMA) sponges, and were reported to contain only up to five phyla in a study of five specimens of three different LMA sponges, Callyspongia vaginalis , Niphates digitalis , and Raspailia topsenti [10] .
In the present study, the bacterial diversity in 10 sponges collected from Jeju Island and Chuja Island, Korea, was investigated using pyrosequencing to reveal the unidentified bacterial diversity of sponges, including Cliona celata , Haliclona cinerea , Lissodendoryx isodictyalis , Spirastrella abata , Spirastrella panis , and Penares incrustans , which have rarely been explored to date. Most of the sponges, except Cliona celata , Lissodendoryx isodictyalis , and Penares incrustans , were found to contain very simple diversity comparable to that of LMA sponges, and some of the predominant OTUs were found to be unrelated to known sequences.
Materials and Methods
- Sample Collection and DNA Extraction
Sponge specimens were collected from Mueung-ri, Daejeongeup, Seogwipo City, Jeju Province (Jeju Island), Korea ( n = 6) in Feb. 2011 by scuba diving, and from around Chuja Island in the South Sea of Korea ( n = 4) in Nov. 2009 by scuba diving or from the intertidal region. Samples were collected aseptically and delivered to the laboratory. Small pieces (about 1 cm 3 in volume) of sponge tissues were washed with sterilized seawater, frozen at -70℃ for 24 h, and then lysophilized at -50℃, 0.033 Mbar for 24 h. The freeze-dried tissues were aseptically homogenized in a mortar. The G-spin genomic DNA extraction kit (Intron, South Korea) was used for DNA extraction.
- Barcoded Pyrosequencing
The V1 to V3 region of the 16S rRNA genes were amplified using primer sets (V1-9F: 5’-X-AC-GAGTTTGATCMTGGCTCAG-3’ and V3-541R: 5’-X-AC-WTTACCGCGGCTGCTGG-3’; X indicates the barcode sequences that were composed of various combinations of six nucleotides to tag different samples). The Genome Sequencer FLX titanium (Roche, Germany) system was used for pyrosequencing, according to the method in the manufacturer’s manuals (Macrogen, Korea).
- Bioinformatic Analysis of Reads
After pyrosequencing, reads were sorted according to tags, and then the primer and tag sequences were cut out. Sequences less than 300 bp in size and with any ambiguous bases, “N”, were discarded from further analyses. Sequences with reverse complement orientation were adjusted into the same orientation. Sequences were aligned and the front and rear parts of sequences that did not match with the other sequences were cut. Putative chimera sequences were detected with the chimer.uchime [6] command of the Mothur package [32] . The QIIME package ver. 1.7.0 was used for to determine the OTUs and for their taxonomic assignment [22] . Sequence names were marked in accordance with the samples and merged into one file. The uclust approach in the QIIME package was used for clustering the OTUs and the determination of representative sequences. A sequence similarity of 97% was used as the criterion for taxonomic assignments. The representative sequences of OTUs were compared with the 97_otus and the 97_otu_taxonomy files of the gg_13_05 version of the Greengenes database [27] and taxonomy was assigned using the RDP classifier [23] . Sequences related to chloroplasts and mitochondria were discarded from further analyses. Diversity estimators, including Chao1, Shannon, and Simpson indices, were determined for subsamples made with sequences of the same numbers. Unifrac analysis [26] was performed to compare the samples. For the Unifrac analysis, similarities of 97%, 94%, 91%, 88%, and 85% were used as criteria for OTU clustering. The PyNAST program [3] was used for the alignment of representative sequences obtained from the clustering, after which the FastTree program [25] was used for construction of a phylogenetic tree using the aligned sequences. UPGMA trees were made from the Unifrac analysis. Some representative sequences were compared against the nucleotide database of GenBank and related sequences were used for the construction of a phylogenetic tree. A neighborjoining tree was constructed with the MEGA5 program [39] , in which the Kimura 2-parameter method [21] was used to calculate a distance matrix.
Results and Discussion
- Sponges and Their Bacterial Communities
All 10 specimens were so-called siliceous sponges (the class Demospongiae ), which have spicules made out of silicon dioxide [41] . Cliona celata is an ecologically/biotechnologically important species that is well known as an excavating and cosmopolitan sponge [50] . The sterol composition and novel fatty acids were discovered from Haliclona cinerea sponges from the Black Sea [8 , 19] . Lissodendoryx isodictyalis is an encrusting sponge with blue-gray color and a strong odor [35] . Lipids with cytotoxicity and other bioactivities were isolated from Spirastrella abata sponges from Korea [15] . The diversity of cultured and uncultured bacteria of Spirastrella abata and Spirastrella panis have been investigated with restriction fragment length polymorphism and denaturing gradient gel electrophoresis (DGGE) [4 , 16] . Bioactive triterpenoids and penaramides were discovered from Penares incrustans sponges [37 , 42] . However, the latter five species have not yet been studied extensively.
As far as we know, bacterial diversity (especially uncultured) of sponges such as Cliona celata , Haliclona cinerea , Lissodendoryx isodictyalis , Spirastrella abata , Spirastrella panis , and Penares incrustans have not yet been or have rarely been published as of yet. Nucleotide records of the GenBank database were surveyed using three or four keywords: (16S) AND (rRNA) AND (“genus name”), or (16S) AND (rRNA) AND (“genus name”) AND (“species name”) of each sponge. The numbers of records that came up for each sponge species or genus name were as follows: Cliona celata , 0; Cliona , 6; Haliclona cinerea , 0; Haliclona , 1839; Halichondria okadai , 122; Halichondria , 608; Hymeniacidon sinapium , 131; Hymeniacidon , 1,245; Lissodendoryx isodictyalis , 0; Lissodendoryx , 1,598; Penares incrustans , 3; Penares , 6; Spirastrella abata , 0; Spirastrella panis , 1; Spirastrella , 5. When species names were used, six sponge species among eight showed no or no more than three records. Even when genus names were used, three sponge genera including Cliona , Penares , and Spirastrella showed no more than six records. This shows that the bacterial diversity of most of the sponges in this study was rarely investigated. Therefore, this study may provide the first observation of the comprehensive bacterial diversity of the sponges examined.
- Unifrac Analysis of Sponges According to Bacterial Diversity
A total of 28,913 sequences were assigned taxonomically after non-bacterial reads including chloroplasts and mitochondria-related sequences were discarded. UPGMA trees from Unifrac analyses showed somewhat different tree topologies in accordance with 16S rRNA gene similarities (91%~85%). One group composed of Sppa2S, Spab1S, and Spab1M showed the same topology in all trees, and one specimen, Pein1M, was always distinguished from the others ( Fig. 1 ). The relationships among other specimens was not robust in the tree topologies that were made based on different criteria of 16S rRNA gene similarities, because the bacterial communities from each of the sponges were highly different from each other. A previous study with 32 sponges showed the bacterial patterns to have no correlation with sponge phylogeny. Sponges from different genera were grouped closer than those from the same genus based on their bacterial compositions [34] . However, Sppa2S, Spab1S, and Spab1M belonged to the same genus and showed a closely related relationship in a clade. Although Spab1S and Spab1M were collected from different sample sites, their bacterial profiles were very similar, which implies that taxon-specific bacterial communities exist in some sponges.
Description of sponge specimens.
PPT Slide
Lager Image
aC, Chuja Island; M, Mureung-ri.
PPT Slide
Lager Image
Weighted UPGMA tree from Unifrac analysis showing relationships between samples according to bacterial profile. OTUs were determined based on 97% similarity. Numbers at the branched points indicate percentages of frequency of each cluster detected in the five UPGMA trees, which were constructed based on different 16S rRNA gene similarities (97%, 94%, 91%, 88%, and 85%). The scale bar represents Unifrac distance.
- Diversity Estimation of Sponges
The sponges used in this study generally showed low and moderate diversity, except for the specimen Pein1M. Several sponges showed extremely low diversity, which could be expressed by Shannon indexes below 3.00 and even down to 0.81 ( Table 2 ). The simplest bacterial community was that of the specimen Haok1M, identified as Halichondria okadai , which contained only seven OTUs. In a previous study, specimens of the sponge Halichondria panacea in the same genus also showed simpler diversity than seawater, with some specimens showing less than five bands in the DGGE patterns [47] . Another sponge showing extremely low diversity was Haci1M, a specimen of the sponge Haliclona cinerea . Two species in the same genus, H. heliophila and H. tubifera , were also reported to have low diversity and species-specific bacterial communities [9] . Spab1M and Spab1S also had low levels of diversity with a single highly abundant OTU. Low diversity in the sponge genus Spirastrella was reported in a previous study [16] , where only five bands were observed in the DDGE patterns. The specimen Clsp1S showed much lower diversity than Clce1S, another specimen of the same genus, Cliona . Rank-accumulated relative abundance curves showed distinguished patterns among the sponge specimens ( Fig. 2 ). The sponges with low diversity had single predominant OTUs with high proportions from 50.8% to 92.6%. Some of them had only a few or dozens of OTUs, but few OTUs showed up in high proportions. The sponges with relatively high diversity such as Liis2M, Clce1S, and Pein1M showed different curve shapes. Pein1M showed evenly distributed abundances among the OTUs. Clce1S also showed a similar moderate slope, except for the OTU ranking first. Liis2M showed evenly distributed abundances among OTUs with several highly abundant OTUs, but the curve saturated quickly. The bacterial diversities of marine sponges examined in previous studies were much higher than those of the present study [17 , 18] , and many more bacterial phyla were detected. In the present study, however, several sponges had very low diversity and were composed of very few OTUs, with one or two OTUs sometimes occupying almost all reads. This might be derived from several factors, such as different sampling sites, experimental procedure, number of reads, and sponge condition. There are reports of low diversity in sponges in which very few OTUs were found in many sponges [34] , so it may not seem to be a rare phenomenon.
Diversity estimation results from pyrosequencing reads.
PPT Slide
Lager Image
aReads were subsampled from each sample and used to calculate diversity estimators. Data from Haok1M were not used for diversity index calculation because the number of reads was too small. OTUs were determined based on 16S rRNA gene similarity of 97%.
PPT Slide
Lager Image
Rank-cumulative relative abundance curve in which relative abundances of OTUs were accumulated sequentially in order of abundance from the largest. The first 19 OTUs are shown. OTUs were determined based on 97% similarity.
- Bacterial Diversity at Phylum and Class Levels
Pein1M was clearly distinguished from the other sponges at the phylum level. It had a bacterial community spanning to 14 phyla, in which the bacteria were evenly distributed (9 phyla were more than 1% in proportion). Clce1S also contained many phyla (16 phyla), but only five of the phyla were more than 1%, and Proteobacteria occupied up to 83.0%. A much higher proportion of Proteobacteria was observed in all specimens except Pein1M (12.5%), ranging from 70.7% to 100.0%. These high proportions were unusual, but not impossible according to previous studies. A high proportion of Proteobacteria (up to 72.6%) was also observed in a previous study [18] . In the low microbial abundance sponges Crella cyathophora , Stylissa carteri , and Niphates digitalis , 78%, 87%, and 93% abundances of Proteobacteria , respectively, were roughly calculated from data on the investigation of bacterial communities [10] . Diversity at the class level of bacteria detected herein is shown in Fig. 3 . Reads related to Alphaproteobacteria were predominant in Spab1S, Spab1M, Sppa2S, Clsp1S, and Haok1M. In contrast, Gammaproteobacteria was the major class in Hysi1M and Liis2M, and they also contained a significant amount of Alphaproteobacteria . Clce1S and Haci1M contained Betaproteobacteria as the major class, and Clce1S also contained Gammaproteobacteria as a significant member.
PPT Slide
Lager Image
Taxonomic assignment results at the class level. Only phyla with an abundance of more than 1% in at least one sample are shown. “Unclassified” means that the sequence was related to reads of which the position has not been verified taxonomically in that level. Square brackets indicate that the name was not verified taxonomically.
- BLASTN and Phylogenetic Analysis of Predominant OTUs
Some representative sequences of highly abundant OTUs were compared against the nucleotide database in GenBank ( Table 3 ). The criterion of 94% similarity was used for OTU determination to reduce the sequence number for comparison. Above all, representative sequences were related to sequences mainly from sponges. Some sequences were related to other marine invertebrates such as coral and octacoral, and to marine environments such as seawater, marine sediments, and microbialites. Most of the sequences were taxonomically unidentified from genus even to class level based on QIIME analysis.
Characteristics of predominant OTUs.
PPT Slide
Lager Image
Ratios in each specimen, taxonomic assignment, and results of the BLASTN comparison with the nucleotide database of GenBank are summarized. aClassification was shown at the genus level. “Unclassified” means that the sequence was related to reads of which the position has not been verified taxonomically in that level. Square brackets indicate that the name was not verified taxonomically. OTU was determined based on a 16S rRNA gene similarity of 94%.
For example, the sequence of denovo103, which was highly predominant in the clade Spirastrella , was related to a clone (Accession No. JQ515708) reported from the coral Montastraea faveolata [20] . Denovo103 was distantly related (93% similarity) to the known sequence JQ515708 and was taxonomically identified only at the class level as Alphaproteobacteria . The representative sequence of denovo694, the most abundance OTU from Haok1M, showed 100% similarity with a sequence (AB054177) obtained from the same sponge species Halichondria okadai from an unpublished study. Interestingly, the second most abundant OTU from Haok1M (denovo197) also showed similarity with a sequence (AB054180, 98%) from the same unpublished study. The two sequences, denovo694 and denovo197, showed high similarity with clones in public databases, but they were very different from strains that have been cultivated and identified taxonomically (in Alphaproteobacteria and Gammaproteobacteria with 90% and 86% similarity, respectively).
Denovo624, the representative sequence of the most abundant OTU of the sponge HaciM, showed the best hit with a sequence (JF824778, 98% similarity) obtained from the same sponge genus, Halichondria tubifera , and belonged to the Betaproteobacteria . Another sequence (KC492704) belonging to Betaproteobacteria was reported to be important in a different sponge, Crambe crambe , although the similarity was low (88%). The betaproteobacterial cells (KC492704) were proven to predominantly inhabit the sponge through cloning-sequencing, pyrosequencing, catalyzed reporter deposition fluorescent in situ hybridization, and transmission electron microscopy, implying that they have potentially important roles in sponges [5] . Denovo205 was the most abundant OTU in Clce1S and had the closest similarity with a sequence (DQ889898, 88%) that had been detected from the octocoral Erythropodium caribaeorum in an unpublished study. The low similarity of denovo205 with other sequences (<88%) implies that new sponge specimens can be reservoirs of novel bacteria that have been unexplored so far.
Many monophyletic and sponge-specific phylogenetic clusters in Bacteria and Archaea have been described in previous papers [38 , 40] . Among them, the cluster within Betaproteobacteria is one of the largest sponge-specific clusters that branched deeply from cultivated strains and is composed only of sequences obtained directly from environments [38] . Denovo205 and denovo624 were also included in this sponge-specific cluster, as shown in Fig. 4 .
PPT Slide
Lager Image
Neighbor-joining tree of predominant OTUs. A 16S rRNA gene similarity of 94% was used for OTU determination to reduce the number of OTUs. Bootstrap values above 70% are shown at the branches.
Through this study, it was revealed that the bacterial communities of some sponges were very simple but differed from one another, and that some predominant OTUs of those sponges showed high dissimilarity with known sequences, phylogenetically forming a sponge-specific clade. Considering that about 8,500 species of sponges have been confirmed taxonomically and the total number of species was estimated to be double that many [43] , while only dozens of studies have been performed to reveal their bacterial diversities, much unexplored microbial diversity of sponges remains to be revealed. Many more studies should be carried out to further understand the microbes in sponges.
This research was a part of the project titled “Depository institution of marine bio-resources”, funded by the Ministry of Oceans and Fisheries, Korea (2009B413-GV-R-0101-S000100) and by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2010-0007045).
Amann RI , Ludwig W , Schleifer KH 1995 Phylogenetic identification andin situdetection of individual microbial cells without cultivation. Microbiol. Rev. 59 143 - 169
Bavestrello G , Arillo A , Calcinai B , Cattaneo-Vietti R , Cerrano C , Gaino E 2000 Parasitic diatoms inside antarctic sponges. Biol. Bull. 198 29 - 33    DOI : 10.2307/1542801
Caporaso JG , Bittinger K , Bushman FD , DeSantis TZ , Andersen GL , Knight R 2010 PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26 266 - 267    DOI : 10.1093/bioinformatics/btp636
Cho H , Park J 2009 Comparative analysis of the community of culturable bacteria associated with sponges,Spirastrella abataandSpirastrella panisby 16S rDNA-RFLP. Kor. J. Microbiol. 45 155 - 162
Croue J , West NJ , Escande ML , Intertaglia L , Lebaron P , Suzuki MT 2013 A single betaproteobacterium dominates the microbial community of the crambescidine-containing spongeCrambe crambe. Sci. Rep. 3 2583 -    DOI : 10.1038/srep02583
Edgar RC , Haas BJ , Clemente JC , Quince C , Knight R 2011 UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27 2194 - 2200    DOI : 10.1093/bioinformatics/btr381
Edwards RA , Rodriguez-Brito B , Wegley L , Haynes M , Breitbart M , Peterson DM 2006 Using pyrosequencing to shed light on deep mine microbial ecology. BMC Genomics 7 57 -    DOI : 10.1186/1471-2164-7-57
Elenkov I , Popov S , Andreev S 1999 Sterols from two Black Sea sponges (Haliclona sp.). Comp. Biochem. Physiol. B 123 357 - 360    DOI : 10.1016/S0305-0491(99)00079-6
Erwin PM , Olson JB , Thacker RW 2011 Phylogenetic diversity, host-specificity and community profiling of sponge-associated bacteria in the northern Gulf of Mexico. PLoS One 6 e26806 -    DOI : 10.1371/journal.pone.0026806
Giles EC , Kamke J , Moitinho-Silva L , Taylor MW , Hentschel U , Ravasi T , Schmitt S 2013 Bacterial community profiles in low microbial abundance sponges. FEMS Microbiol. Ecol. 83 232 - 241    DOI : 10.1111/j.1574-6941.2012.01467.x
Hentschel U , Hopke J , Horn M , Friedrich AB , Wagner M , Hacker J , Moore BS 2002 Molecular evidence for a uniform microbial community in sponges from different oceans. Appl. Environ. Microbiol. 68 4431 - 4440    DOI : 10.1128/AEM.68.9.4431-4440.2002
Hentschel U , Usher KM , Taylor MW 2006 Marine sponges as microbial fermenters. FEMS Microbiol. Ecol. 55 167 - 177    DOI : 10.1111/j.1574-6941.2005.00046.x
Hill M , Hill A , Lopez N , Harriott O 2006 Sponge-specific bacterial symbionts in the Caribbean sponge,Chondrilla nucula(Demospongiae,Chondrosida). Mar. Biol. 148 1221 - 1230    DOI : 10.1007/s00227-005-0164-5
Hoffmann F , Larsen O , Thiel V , Rapp HT , Pape T , Michaelis W , Reitner J 2005 An anaerobic world in sponges. Geomicrobiol. J. 22 1 - 10    DOI : 10.1080/01490450590922505
Jang KH , Lee Y , Sim CJ , Oh KB , Shin J 2012 Bioactive lipids from the spongeSpirastrella abata. Bioorg. Med. Chem. Lett. 22 1078 - 1081    DOI : 10.1016/j.bmcl.2011.11.105
Jeong EJ , Im CS , Park JS 2010 A comparison of bacterial diversity associated with the spongeSpirastrella abatadepending on RFLP and DGGE. Kor. J. Microbiol. 46 366 - 374
Jeong IH , Kim KH , Lee HS , Park JS 2014 Analysis of bacterial diversity in sponges collected from Chuuk and Kosrae Islands in Micronesia. J. Microbiol. 52 20 - 26    DOI : 10.1007/s12275-014-3619-x
Jeong IH , Kim KH , Park JS 2013 Analysis of bacterial diversity in sponges collected off Chujado, an island in Korea, using barcoded 454 pyrosequencing: analysis of a distinctive sponge group containingChloroflexi. J. Microbiol. 51 570 - 577    DOI : 10.1007/s12275-013-3426-9
Joh YG , Elenkov IJ , Stefanov KL , Popov SS , Dobson G , Christie WW 1997 Novel di-, tri-, and tetraenoic fatty acids withbis-methylene-interrupted double-bond systems from the spongeHaliclona cinerea. Lipids 32 13 - 17
Kimes NE , Johnson WR , Torralba M , Nelson KE , Weil E , Morris PJ 2013 TheMontastraea faveolatamicrobiome: ecological and temporal influences on a Caribbean reefbuilding coral in decline. Environ. Microbiol. 15 2082 - 2094    DOI : 10.1111/1462-2920.12130
Kimura M 1980 A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16 111 - 120    DOI : 10.1007/BF01731581
Kuczynski J , Stombaugh J , Walters WA , Gonzalez A , Caporaso JG , Knight R 2012 Using QIIME to analyze 16S rRNA gene sequences from microbial communities. Curr. Protoc. Microbiol. Chapter 1: Unit 1E 5
Lan Y , Wang Q , Cole JR , Rosen GL 2012 Using the RDP classifier to predict taxonomic novelty and reduce the search space for finding novel organisms. PLoS One 7 e32491 -    DOI : 10.1371/journal.pone.0032491
Lee OO , Wang Y , Yang J , Lafi FF , Al-Suwailem A , Qian PY 2011 Pyrosequencing reveals highly diverse and speciesspecific microbial communities in sponges from the Red Sea. ISME J. 5 650 - 664    DOI : 10.1038/ismej.2010.165
Liu K , Linder CR , Warnow T 2011 RAxML and FastTree: comparing two methods for large-scale maximum likelihood phylogeny estimation. PLoS One 6 e27731 -    DOI : 10.1371/journal.pone.0027731
Lozupone C , Lladser ME , Knights D , Stombaugh J , Knight R 2011 UniFrac: an effective distance metric for microbial community comparison. ISME J. 5 169 - 172    DOI : 10.1038/ismej.2010.133
McDonald D , Price MN , Goodrich J , Nawrocki EP , DeSantis TZ , Probst A 2012 An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6 610 - 618    DOI : 10.1038/ismej.2011.139
Park S , Yokoa A , Ioh T , Park JS 2011 Sphingomonas jejuensissp. nov., isolated from marine spongeHymeniacidon flavia. J. Microbiol. 49 238 - 242    DOI : 10.1007/s12275-011-0500-z
Pile A , Patterson M , Witman J 1996 In situgrazing on plankton <10 µm by the boreal spongeMycale lingua. Mar. Ecol. Prog. Ser. 141 95 - 102    DOI : 10.3354/meps141095
Roh SW , Kim KH , Nam YD , Chang HW , Park EJ , Bae JW 2010 Investigation of archaeal and bacterial diversity in fermented seafood using barcoded pyrosequencing. ISME J. 4 1 - 16    DOI : 10.1038/ismej.2009.83
Santavy DL , Willenz P , Colwell RR 1990 Phenotypic study of bacteria associated with the Caribbean sclerosponge,Ceratoporella nicholsoni. Appl. Environ. Microbiol. 56 1750 - 1762
Schloss PD , Westcott SL , Ryabin T , Hall JR , Hartmann M , Hollister EB 2009 Introducing Mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75 7537 - 7541    DOI : 10.1128/AEM.01541-09
Schmitt S , Hentschel U , Taylor M 2012 Deep sequencing reveals diversity and community structure of complex microbiota in five Mediterranean sponges. Hydrobiologia 687 341 - 351    DOI : 10.1007/s10750-011-0799-9
Schmitt S , Tsai P , Bell J , Fromont J , Ilan M , Lindquist N 2012 Assessing the complex sponge microbiota: core, variable and species-specific bacterial communities in marine sponges. ISME J. 6 564 - 576    DOI : 10.1038/ismej.2011.116
Sears MA , Gerhart DJ , Rittschof D 1990 Antifouling agents from marine spongeLissodendoryx isodictyaliscarter. J. Chem. Ecol. 16 791 - 799    DOI : 10.1007/BF01016489
Sfanos K , Harmody D , Dang P , Ledger A , Pomponi S , McCarthy P , Lopez J 2005 A molecular systematic survey of cultured microbial associates of deep-water marine invertebrates. Syst. Appl. Microbiol. 28 242 - 264    DOI : 10.1016/j.syapm.2004.12.002
Shoji N , Umeyama A , Mooki S , Arihara S , Ishida T , Nomoto K 1992 Potent inhibitors of histamine release, two novel triterpenoids from the Okinawan marine spongePenares incrustans. J. Nat. Prod. 55 1682 - 1685    DOI : 10.1021/np50089a021
Simister RL , Deines P , Botte ES , Webster NS , Taylor MW 2012 Sponge-specific clusters revisited: a comprehensive phylogeny of sponge-associated microorganisms. Environ. Microbiol. 14 517 - 524    DOI : 10.1111/j.1462-2920.2011.02664.x
Tamura K , Peterson D , Peterson N , Stecher G , Nei M , Kumar S 2011 MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28 2731 - 2739    DOI : 10.1093/molbev/msr121
Taylor MW , Radax R , Steger D , Wagner M 2007 Spongeassociated microorganisms: evolution, ecology, and biotechnological potential. Microbiol. Mol. Biol. Rev. 71 295 - 347    DOI : 10.1128/MMBR.00040-06
Uriz MJ , Turon X , Becerro MA , Agell G 2003 Siliceous spicules and skeleton frameworks in sponges: origin, diversity, ultrastructural patterns, and biological functions. Microsc. Res. Tech. 62 279 - 299    DOI : 10.1002/jemt.10395
Ushio-Sata N , Matsunaga S , Fusetani N , Honda K , Yasumuro K 1996 Penaramides, which inhibit binding of ω-conotoxin GVIA to N-type Ca2+channels, from the marine spongePenaresaff.incrustans. Tetrahedron Lett. 37 225 - 228    DOI : 10.1016/0040-4039(95)02134-5
Van Soest RW , Boury-Esnault N , Vacelet J , Dohrmann M , Erpenbeck D , De Voogd NJ 2012 Global diversity of sponges (Porifera). PLoS One 7 e35105 -    DOI : 10.1371/journal.pone.0035105
Webster NS , Cobb RE , Soo R , Anthony SL , Battershill CN , Whalan S , Evans-Illidge E 2011 Bacterial community dynamics in the marine spongeRhopaloeides odorabileunderin situandex situcultivation. Mar. Biotechnol. (NY) 13 296 - 304    DOI : 10.1007/s10126-010-9300-4
Webster NS , Taylor MW , Behnam F , Lucker S , Rattei T , Whalan S 2010 Deep sequencing reveals exceptional diversity and modes of transmission for bacterial sponge symbionts. Environ. Microbiol. 12 2070 - 2082
White JR , Patel J , Ottesen A , Arce G , Blackwelder P , Lopez JV 2012 Pyrosequencing of bacterial symbionts withinAxinella corrugatasponges: diversity and seasonal variability. PLoS One 7 e38204 -    DOI : 10.1371/journal.pone.0038204
Wichels A , Wurtz S , Dopke H , Schutt C , Gerdts G 2006 Bacterial diversity in the breadcrumb spongeHalichondria panicea(Pallas). FEMS Microbiol. Ecol. 56 102 - 118    DOI : 10.1111/j.1574-6941.2006.00067.x
Wilkinson C , Nowak M , Austin B , Colwell R 1981 Specificity of bacterial symbionts in Mediterranean and Great Barrier Reef sponges. Microbial Ecol. 7 13 - 21    DOI : 10.1007/BF02010474
Wilkinson CR 1983 Net primary productivity in coral reef sponges. Science 219 410 - 412    DOI : 10.1126/science.219.4583.410
Xavier JR , Rachello-Dolmen PG , Parra-Velandia F , Schonberg CH , Breeuwer JA , van Soest RW 2010 Molecular evidence of cryptic speciation in the “cosmopolitan” excavating spongeCliona celata(Porifera,Clionaidae). Mol. Phylogenet. Evol. 56 13 - 20    DOI : 10.1016/j.ympev.2010.03.030
Yoon BJ , You HS , Lee DH , Oh DC 2011 Aquimarina spongiaesp. nov., isolated from marine spongeHalichondria oshoro. Int. J. Syst. Evol. Microbiol. 61 417 - 421    DOI : 10.1099/ijs.0.022046-0