Eutrophication has no short-term effect on the Cymbastela stipitata holobiont

Levels of nitrogen in coastal areas have been rapidly increasing due to accumulative inputs of sewage and terrigenous sediments carrying fertilizers. Sponges have an immense filtering capacity and may be directly impacted (positively or negatively) by elevated concentrations of nitrogen. Sponges also host a wide diversity of microbes involved in nitrogen metabolism, yet little is known about the effects of nitrogen loading on these symbiotic partnerships. Manipulative experiments were undertaken to examine the potential effects of excess nitrogen (up to 240 μM) on microbial symbiosis in the abundant sponge species Cymbastela stipitata. Microbial composition and activity were examined using 454-pyrotag sequencing of DNA- and RNA-derived samples. Despite the high levels of nitrogen exposure (up to 124-fold above ambient), sponges appeared visibly unaffected at all treatment concentrations. At the phylum level, the microbial community was consistent between all sponge samples regardless of nitrogen treatment, with Cyanobacteria and Thaumarchaeota being the dominant taxa. Higher microbial diversity was observed at the operational taxonomic units (OTU) level (97% sequence similarity), with only 40% of OTUs shared between samples from all treatments. However, a single cyanobacterial OTU dominated the community of all individuals (average 73.5%) and this OTU did not vary with nitrogen treatment. The conserved microbial community in all sponges irrespective of nitrogen treatment highlights the stability of the sponge-microbe relationship and indicates that the holobiont is resistant to short pulses of nitrogen at levels mimicking sewage effluent.

Introduction

Nitrogen is generally considered limiting in marine systems although levels in coastal areas have been increasing due to inputs of terrigenous sediments, carrying fertilizers and sewage (Brodie and Mitchell, 2005; Burford et al., 2008). Increased eutrophication has led to the degradation of coastal reefs worldwide via: (1) increased nutrient levels causing trophic shifts, (2) light attenuation from increased turbidity which reduces photosynthesis, and (3) increased sedimentation which reduces larval settlement and increases mortality (Fabricius, 2011). In particular, sewage effluent has been linked to phytoplankton blooms (Caperon et al., 1971) and increased benthic microalgae biomass (Smith et al., 1981) both of which contribute to reef degradation. Darwin Harbor, located in the Northern Territory of Australia, has generally been considered a pristine Harbor with low nutrient inputs (McKinnon et al., 2006). However, recent studies have linked high nitrogen levels in tidal creeks to sewage effluent (Smith et al., 2011; Burford et al., 2012), making point-source sewage discharges and increased urbanization the largest water quality management issues in the Harbor (Fortune and Maly, 2009).

Sponges in Darwin Harbor comprise a diverse and ecologically significant component of benthic communities (Alvarez and Hooper, 2011). Due to their immense filtering capacity [an individual sponge can filter 1000s of liters of seawater per day (Reiswig, 1971)], sponges will be highly exposed to these elevated levels of nitrogenous compounds. Other benthic invertebrates, such as corals, display adverse effects from increased nutrients including dissolved inorganic nitrogen levels above 4 μM which has been shown to cause reduced reproduction and growth (Koop et al., 2001) and increases in disease severity and progression (Bruno et al., 2003; Voss and Richardson, 2006). However, it is still unclear whether sponges share this sensitivity to elevated levels of nitrogen.

Sponges are amongst the oldest metazoans and their evolutionary success has been linked to their intimate association with microorganisms (Webster and Taylor, 2012). Sponges host a diverse array of microorganisms, with 32 different bacterial phyla and both major lineages of archaea identified to date (Taylor et al., 2007; Schmitt et al., 2012; Webster and Taylor, 2012). While it is known that sponges host microorganisms involved in nitrogen metabolism (Wilkinson and Fay, 1979; Southwell et al., 2008; Weisz et al., 2008; Hoffmann et al., 2009; Schläppy et al., 2009; Mohamed et al., 2010; Off et al., 2010), little is known about the effects of nitrogen loading on the sponge-microbe symbiosis, particularly at the high levels found in sewage effluent. However, in a recent study of sponges in Brazil, a higher diversity of sponge-associated Crenarchaeota was observed from the more polluted areas of Rio de Janeiro, suggesting that sponges may have the ability to alter their community in response to eutrophic environments (Turque et al., 2010).

Here we exposed the sponge Cymbastela stipitata (Bergquist and Tizard, 1967), one of the most abundant sponge species found in Darwin Harbor (Alvarez and Hooper, 2009), to concentrations of nitrogen ranging from 2 μM (ambient) to 240 μM and performed 454-pyrotag sequencing of the 16S rRNA gene in DNA- and RNA-derived samples to determine how anthropogenic inputs of nitrogen affect sponge microbial symbiosis.

Materials and Methods

Sponge Collection

Individuals of C. stipitata were collected from Channel Island, Northern Territory (12° 33′ 02. 4′′ S, 130° 52′ 31. 3′′ E) during the lowest tide of June 2013. Sponges were collected on foot using a hammer and chisel, taking special care to ensure some substrate was included to avoid injuring the sponge tissue. Sponges were immediately transported to the Darwin Aquaculture Centre (DAC) on Channel Island where they were placed into 12 × 30 l flow-through aquaria (flow rate of 400 ml min^-1) and maintained under natural lighting conditions. DAC seawater supply is pumped in from a pipe 40 m off Channel Island and filtered to 5 μm to remove large particulates, yet leaving sponges with a sufficient food supply (Reiswig, 1971). Sponges were kept under these conditions for one week to allow them to acclimatize prior to commencing the experiment.

Experimental Design

To assess the effect of elevated nitrogen levels on the sponge-microbe symbiosis, C. stipitata individuals were exposed to four different nutrient treatments using the water-soluble plant fertilizer Thrive^® (Yates, NPK; 27:5:5:9 and trace elements). For the first two days, three different stock concentrations (low, medium and high) of Thrive were continuously pumped into treatment tanks at 4 ml min^-1 at final total inorganic nitrogen concentrations of 2 μM (Ambient), 120 μM (Low), 160 μM (Medium) and 240 μM (High), mimicking a sewage discharge event (Power and Water Corporation, 2006). After 48 h, Thrive input was decreased for an additional 5 days (2 μM (Ambient), 3, 5, and 15 μM), reflecting a flood plume event (Kroon et al., 2012). The experimental design comprised three replicate tanks per treatment, each holding three sponges. One individual sponge per tank/treatment was sacrificed at each sampling day: 0, 2, and 7. Samples were snap frozen in liquid nitrogen and stored at -80°C for further analysis. Histological specimens were prepared as previously described in Luter et al. (2011) to compare the internal tissue structure between samples in different treatments.

After the first day of the experiment, one of the sponges from the high nitrogen treatment displayed substantial tissue discoloration. It was immediately removed from the experiment to avoid compromising the health of the other sponges in the tank.

Nutrient Analysis

Dissolved nutrient levels [ammonium, nitrite, nitrate, phosphate, dissolved organic carbon (DOC)] were monitored throughout the experiment to ensure treatment levels were maintained. Seawater samples were hand filtered through 0.45 μm cellulose acetate filter cartridges (Sartorius MiniStart) into acid-washed 15 ml tubes and stored at -20° C prior to dissolved nutrient analysis. Filtered samples for DOC were acidified with 100 μl HCl prior to freezing. Nutrient samples were analyzed by the analytical services laboratory at the Australian Institute of Marine Science (AIMS, Townsville).

DNA/RNA Extractions and 454-Pyrosequencing

DNA and RNA were simultaneously extracted from individual sponge tissue samples (≤30 mg) using the AllPrep DNA/RNA Mini kit (Qiagen), following the manufacture’s protocols. DNA contamination was eliminated from RNA samples using the TURBO DNA-free kit (Ambion). The quantity and purity of DNA and RNA was determined using gel electrophoresis [1.1% agarose gels containing GelRed (Biotium)] and a NanoDrop 2000 spectrophotometer (Thermo Scientific). DNA and RNA were sent to MR DNA (http://www.mrdnalab.com) where RNA was reverse transcribed into cDNA and the 16S rRNA gene was amplified using both bacterial (27Fmod and 530 R; Dowd et al., 2008) and archaeal [arch344F (Giovannoni et al., 1988) and arch915R (Stahl and Amann, 1991)] primers and sequenced using the Roche GS-FLX platform. Raw sequence data have been submitted to the NCBI Sequence Read Archive (SRA) under the accession number: SRP038768.

Processing of Raw Sequence Data and Taxonomic Assignment

Pyrosequencing flowgram files (SFF) were processed using Mothur (Schloss et al., 2009). Flowgrams were filtered and denoised using the AmpliconNoise (Quince et al., 2011) function in Mothur. If sequences were less than 200 bp, contained ambiguous characters, or had homopolymers longer than 8 bp, more than one MID mismatch, or more than two mismatches to the reverse primer sequence, they were removed from the analysis. Sequences deemed unique by Mothur were aligned against a SILVA alignment (http://www.mothur.org/wiki/Silva_reference_alignment). Chimeric sequences were removed using UCHIME (Edgar et al., 2011), and samples were sub sampled down to the lowest read number and grouped into 97% OTUs based on pairwise distance matrices created in Mothur. OTUs were taxonomically classified at 97% sequence similarity with the RDP Classifier (Wang et al., 2007) using the SILVA database as a reference.

Data Analyses

The unconstrained principal coordinate analysis (PCO) was used to visually compare communities, while PERMNOVA was used to test differences in community structure. All analyses were performed using PRIMER/PERMANOVA+ (Plymouth, UK).

Clone Library Construction and Phylogenetic Analysis

A clone library was constructed to obtain a full length sequence of the dominant cyanobacterial OTU identified in 454-pyrosequencing analysis. The 16S rRNA gene was amplified using the universal bacterial primers 63f (Marchesi et al., 1998) and 1492r (Lane, 1991). PCR reactions contained 5 μl Bio-X-Act buffer, 0.10 μl of each primer (100 pmol μl^-1), 0.4 μl BSA (10 mg ml^-1), 0.5 μl Bio-X-act short DNA polymerase (Bioline, London, UK) and 1 μl DNA. Reactions were made up to 50 μl using high pure water. The conditions consisted of 1 cycle at 95°C for 5 min, followed by 30 cycles at 95°C for 30 s, 56°C for 30 s and 72°C for 2 min, and a final elongation at 72°C for 10 min. PCR products from the three control sponges samples were pooled and gel purified using a Wizard SV Gel and PCR Clean-Up Kit following the manufacture’s protocol (Promega, Madison, WI, USA). The purified PCR product was cloned using a TOPO-TA cloning kit (Invitrogen, Carlsbad, CA, USA) according to the manufacture’s protocol. Plasmids were checked for inserts by PCR amplification using M13 forward and reverse primers. Restriction digests using HAEIII and HhaI (New England Biolabs Inc., Ipswich, MA, USA) were used to determine OTUs. Twenty clones were screened and triplicates of the dominant OTU pattern was sent for sequencing using the forward and reverse primers at the Australian Genomic Research Facilities (Brisbane, Australia). The sequence was submitted to Genbank under the accession number: KJ719259.

The clone sequence was compared to available databases by using BLAST to determine the nearest relatives and the percent similarity. Sequences from the nearest relatives, other sponge-specific cyanobacterial sequences and sequences representing the diversity of Cyanobacteria were compiled and aligned using the SINA web aligner (Pruesse et al., 2012). Maximum likelihood analysis and tree construction were performed using MEGA version 6.06 (Tamura et al., 2013), following the procedures of Hall (2013). The Kimura 2-parameter model with a gamma distribution (four categories) and invariant sites was found to be the most appropriate model. The reliability of the tree was tested by computing 500 bootstrap replicates starting with a neighbor-joining tree, and using the nearest-neighbor interchange (NNI) tree search option.

Chlorophyll (CHL) a Extraction and Analysis

As an estimate of potential photosymbiont abundance, Chl a was extracted from pre-weighed tissue samples (stored at -80°C) using 90% acetone:water overnight at 4°C. Samples were centrifuged and the supernatant was transferred to a spectrophotometer cuvette. Absorbance readings were taken at 750, 664, 647, and 630 nm using the Helios Gamma2 spectrometer. Chl a levels were calculated following equations from Parsons et al. (1984) and standardized using tissue weight.

Results

Nutrient levels from the ambient treatment were used to determine fold differences in total inorganic nutrient levels achieved in the treatment exposures over the course of the experiment. During the first 48 h of the experiment, C. stipitata individuals in the highest treatment were exposed to more than a 100-fold enrichment of total inorganic nitrogen (Table 1). For the remainder of the experiment levels were reduced, with sponges in the highest treatment exposed to a 13.5-fold enrichment of total inorganic nitrogen (Table 1). Despite the extremely high nitrogen exposures, C. stipitata individuals showed no visible effects on the external (Figures 1A–C) or internal (Figures 1D–F) tissues.

FIGURE 1

FIGURE 1. Representative images, both external and internal, of C. stipitata from the day 0 control (A,D), T = 2, High (B,E) and T = 7, High (C,F). A representative choanocyte chamber (Cc), sponge fiber (Sf), and spicule (S) are depicted.

TABLE 1

TABLE 1. Nutrient parameters measured over the course of the 7 day experiment (three replicate water samples averaged per treatment).

After filtering and quality control, 316,494 bacterial sequences and 252,640 archaeal sequences were obtained from 52 samples (26 DNA- and 26 RNA-derived). Bacterial sequences averaged 6,086 (±3,322 1SD) reads per sample, while archaeal sequences averaged 3,770 (±2,912 1SD) reads per sample. We identified 1,471 bacterial OTUs and 141 archaeal OTUs (97% sequence similarity) spanning nine bacterial phyla and two major lineages of the archaea (Thaumarchaeota and Euryarchaeota).

Taxonomic assignment of the OTUs revealed that the bacterial community from both the DNA- and RNA-derived samples was conserved between treatments and time points at the phyla level (Figures 2A,B). Both communities were dominated by Cyanobacteria, with 70% of 16S rRNA-genes, and 86% of transcripts (averaged across treatments) being assigned to Cyanobacteria. The next most abundant phyla observed across all DNA- and RNA-derived samples was the Proteobacteria, particularly members of the Alphaproteobacteria (14 and 4%) and Gammaproteobacteria (14 and 5%; Figures 2A,B). All other Phyla observed (Bacteroidetes, Deltaproteobacteria, and Verrucomicrobia) occurred in lower abundances and varied between treatments, but not consistently across replicates.

FIGURE 2

FIGURE 2. Relative abundance of each bacterial phyla, plus class for Proteobacteria, within each treatment for (A) DNA-derived and (B) RNA-derived sequences. The three replicate individuals per treatment were averaged. For clarity, only the 50 most abundant operational taxonomic units (OTUs) were used to create the graphs.

Interestingly, of the 169 cyanobacterial-assigned OTUs, a single OTU (OTU0001) that was 94% similar to the Cyanobacteria symbiont of Mycale (Carmia) hentscheli (AJ292195) comprised 99% of the total Cyanobacteria reads. A similar trend was also observed for the other dominant phyla. OTU0002, which was 98% similar to a Phyllobacterium of Muricea elongata (DQ917852) comprised 71% of the Alphaproteobacteria reads, while OTU003 comprised 47% of the Gammaproteobacteria reads and shared 99% similarity with a symbiont of C. concentrica (AY942761).

Phylogenetic analysis of the full length sequence of the dominant cyanobacterial OTU (OTU0001) revealed its position within group 5 Cyanobacteria (Figure 3), comprising the Chroococcales, Oscillatoriales, Pleurocapsales, and Prochlorales (Honda et al., 1999). While most sponge-derived cyanobacteria fall within clusters of sponge-specific sequences (Simister et al., 2012a) and many affiliate with group 6 Cyanobacteria (Steindler et al., 2005), symbionts within group 5 Cyanobacteria have previously been reported from Cymbastela sp., M. hentsheli, and Lendenfeldia dendyi (Figure 3).

FIGURE 3

FIGURE 3. Maximum-likelihood phylogenetic tree of the C. stipitata cyanobacterial clone, indicated by boldfacing, based on 16S rRNA sequences. Other sponge-associated Cyanobacteria are indicated by asterisks (*). The numbers at the nodes are percentages, indicating the levels of bootstrap support based on analysis of 500 resampled data sets. Only values >50% are shown. The scale bar represents 0.1 substitutions per nucleotide position.

The DNA-derived archaeal community of C. stipitata from every treatment was almost exclusively comprised of Thaumarchaeota (Figure 4). The other archaeal OTUs represented members of the Euryachaeota; however, combined they comprised <1% of the total community. Interestingly, OTUs classified as Euryachaeota were not observed in samples from the highest nitrogen treatment. The dominant Thaumarchaeota OTU shared 99% similarity with a Montipora hispida symbiont (JX021862).

FIGURE 4

FIGURE 4. Relative abundance of archaeal OTUs from DNA-derived sequences. The three replicate individuals per treatment were averaged.

Principal coordinate analysis of the OTUs from DNA-derived samples revealed a large degree of variation between samples, with 27.7% of the total variation in community composition explained within the first two factors (Figure 5). All samples shared 40% similarity in community composition and no particular groupings according to nitrogen treatment were evident. The DNA-derived community was not significantly different between nutrient treatments (PERMNOVA, p = 0.075). PCO of the OTUs from RNA-derived samples also showed no distinctive clustering of samples according to treatment or time point (Figure 6) and explained a lower amount of the total variation in the first two factors (24.3%). The RNA-derived community shared only 20% similarity between all samples and was also not significantly different between nitrogen treatments (PERMNOVA, p = 0.216).

FIGURE 5

FIGURE 5. Principal coordinate analysis (PCO) of the OTUs from the DNA-derived samples. Control samples are from day 0, solid symbols represent day 2 samples and open symbols represent day 7 samples.

FIGURE 6

FIGURE 6. PCO of the OTUs from the RNA-derived samples. Control samples are from day 0, solid symbols represent day 2 samples and open symbols represent day 7 samples.

Chlorophyll a concentrations of individuals among treatments were not significantly different (PERMNOVA, p = 0.286). Over the duration of the experiment concentrations ranged from 9 ± 0.2 to 27 ± 0.6 μg Chl a g^-1 sponge tissue (Figure 7). Both the lowest and highest levels of Chl a were observed at day 7, with sponges from the Low treatment displaying the lowest concentrations, while sponges from the Ambient treatment displayed the highest concentrations.

FIGURE 7

FIGURE 7. Mean (±SE) concentration of Chl a in C. stipitata individuals from each treatment (n = 3 for each).

Discussion

This is the first study to comprehensively describe the microbial community of C. stipitata at the gene and transcript level using 454-pyrotag sequencing. Both the DNA- and RNA-derived samples revealed that the microbial communities (both total and active) were conserved in all individuals at the phyla level and dominated by Cyanobacteria and Thaumarchaeota, although Alphaproteobacteria and Gammaproteobacteria were also abundant. Previous studies of sponges from the Cymbastela genus described similar microbial compositions, with representatives of Cyanobacteria, Alphaproteobacteria, and Gammaprotobacteria being reported in both C. coralliophila (Webster et al., 2013) and C. concentrica (Taylor et al., 2004, 2005). While the representation of phyla was similar across the genus, Cyanobacteria appeared to dominate a larger percentage of the total community of C. stipitata (70% total, 86% active) than C. coralliophila (20%) or C. concentria (<5%; Fan et al., 2012), although differences in methodology preclude a direct comparison. While it is clear that the Cymbastela genus forms an association with cyanobaterial symbionts, phylogenetic analysis revealed the dominant symbiont of C. stipitata to be more closely related to M. hentscheli than the other Cymbastela species; however, both cyanobacterial symbionts fall into group 5 Cyanobacteria (Honda et al., 1999). The archaeal community of the Cymbastela genus appears more variable, with metagenomic studies revealing the conditional association of Thaumarchaeota in C. concentrica, yet, Thaumarchaeota were absent in C. coralliophilia (Fan et al., 2012). Interestingly, Cyanobacteria and Thaumarchaeota are well known for their roles in nitrogen metabolism including nitrogen fixation (Wilkinson and Fay, 1979) and ammonia oxidation (Pester et al., 2011), respectively. The dominance of these two groups in the active community of C. stipitata suggests they may be playing an important functional role in this sponge.

Examination of the DNA- and RNA-derived microbial communities at the OTU level revealed a low degree of similarity between individual samples. However, given that the majority (73.5%) of the total community (averaged across DNA- and RNA-derived samples) was comprised of a single cyanobacterial OTU (OTU0001), the variation between individuals is likely due to rare OTUs found in lower abundances. The dominance of a single symbiont is not uncommon among sponge species. For instance, a single gammaproteobacterial OTU made up 49% of all 454 tag sequences (Webster et al., 2010) and >90% of all clone library sequences (Luter et al., 2010) in the sponge Ianthella basta. Similarly, a specialist Gammaproteobacteria symbiont dominated more than 50% of the Haliclona tubifera microbial community (Erwin et al., 2011). The presence of a few dominant symbionts in the sponge hosts studied may be a key factor in their ability to withstand a range of environmental stressors. For instance, the three dominant symbionts of I. basta remain stable regardless of their disease status (Luter et al., 2010) or exposure to elevated temperature and sedimentation (Luter et al., 2012).

Cymbastela stipitata exposed to levels of nitrogen consistent with sewage effluent discharge and flood plume events remained visibly unaffected after 7 days. In addition, phyla level comparisons revealed that both the total and active microbial community was conserved between individuals, regardless of the nitrogen treatment. It is possible that this high tolerance is linked to adaptation mechanisms of the host and/or microbial symbionts that could not be detected with amplicon sequencing. The microbial stability observed in C. stipitata during elevated nutrient exposure is consistent with findings for Rhopaloeides odorabile (Simister et al., 2012b) and Aplysina cauliformis (Gochfeld et al., 2012). However, in both of those studies sponges were exposed to inorganic nitrogen levels 11.4-fold (Simister et al., 2012b) and 4.8-fold (nitrate; Gochfeld et al., 2012) above ambient, whereas sponges from the highest treatment in this experiment experienced enrichments ten times higher (124-fold above ambient), demonstrating an incredibly high resistance to elevated nitrogen concentrations. Given sponges can be abundant on inshore reefs (Wilkinson and Cheshire, 1989) more commonly associated with eutrophic conditions (Brodie and Mitchell, 2005), it is possible that they and their symbionts have developed a tolerance to higher levels of nitrogenous inputs.

This study provides a comprehensive microbial community analysis of an abundant sponge species, revealing the dominance of two symbionts known to be involved in nitrogen metabolism and highlighting the nitrogen tolerance of the C. stipitata holobiont.

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.

Acknowledgments

This study was supported by a NAMRA Postdoctoral Fellowship awarded to Heidi M. Luter. We acknowledge the staff of DAC for their assistance with aquarium facilities and A. Rose, S. Every, A. Skillingtion, D. Blondeau-Patissier, A. McInnes, and C. Cameron for their assistance with sponge collections. W. Whitburn, S. Harries, and K. Makarous assisted with sample processing. Nicole S. Webster was funded through an Australian Research Council Future Fellowship (FT120100480).

References

Alvarez, B., and Hooper, J. N. A. (2009). Taxonomic revision of the order Halichondrida (Porifera: Demospongiae) from northern Australia. Family Axinellidae. Beagle 25,17–42.

Alvarez, B., and Hooper, J. N. A. (2011). Taxonomic revision of the order Halichondrida (Porifera: Demospongiae) from northern Australia. Family Halichondriidae. Beagle 27, 55–84.

Bergquist, P. R., and Tizard, C. A. (1967). Australian intertidal sponges from the Darwin area. Micronesica 3, 175–202.

Brodie, J. E., and Mitchell, A. W. (2005). Nutrients in Australian tropical rivers: changes with agricultural development and implications for receiving environments. Mar. Freshwater Res. 56, 279–302. doi: 10.1071/MF04081

CrossRef Full Text

Bruno, J. F., Petes, L. E., Harvell, D. C., and Hettinger, A. (2003). Nutrient enrichment can increase the severity of coral diseases. Ecol. Lett. 6, 1056–1061. doi: 10.1046/j.1461-0248.2003.00544.x

CrossRef Full Text

Burford, M. A., Alongi, D. M., McKinnon, A. D., and Trott, L. A. (2008). Primary production and nutrients in a tropical macrotidal estuary, Darwin Harbour, Australia. Estuar. Coast. Shelf Sci. 79, 440–448. doi: 10.1016/j.ecss.2008.04.018

CrossRef Full Text

Burford, M. A., Revill, A. T., Smith, J., and Clementson, L. (2012). Effect of sewage nutrients on algal production, biomass and pigments in tropical tidal creeks. Mar. Pollut. Bull. 64, 2671–2680. doi: 10.1016/j.marpolbul.2012.10.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Caperon, J., Cattell, S. A., and Krasnick, G. (1971). Phytoplankton kinetics in a subtropical estuary: eutrophication. Limnol. Oceanogr. 16, 599–607. doi: 10.4319/lo.1971.16.4.0599

CrossRef Full Text

Dowd, S. E., Callaway, T. R., Wolcott, R. D., Sun, Y., McKeehan, T., Hagevoort, R. G., et al. (2008). Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol. 8:125. doi: 10.1186/1471-2180-8-125

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Erwin, P. M., Olson, J. B., and Thacker, R. W. (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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fabricius, K. (2011). “Factors determining the resilience of coral reefs to eutrophication: a review and conceptual model,” in Coral Reefs: An Ecosystem in Transition, eds Z. Dubinsky and N. Stambler (Houten: Springer), 493–505.

Fan, L., Reynolds, D., Liu, M., Stark, M., Kjelleberg, S., Webster, N. S., et al. (2012). Functional equivalence and evolutionary convergence in complex communities of microbial sponge symbionts. Proc. Natl. Acad. Sci. U.S.A. 109, E1878–E1887. doi: 10.1073/pnas.1203287109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fortune, J., and Maly, G. (2009). Phase One Report – Towards the Development of a Water Quality Protection Plan for the Darwin Harbour Region. Darwin, NT: Aquatic Health Unit, Department of Natural Resources, Environment, the Arts and Sport.

Giovannoni, S. J., DeLong, E. F., Olsen, G. J., and Pace, N. R. (1988). Phylogenetic group-specific oligodeoxynucleotide probes for identification of single microbial cells. J. Bacteriol. 170, 720–726.

Pubmed Abstract | Pubmed Full Text

Gochfeld, D. J., Easson, C. G., Freeman, C. J., Thacker, R. W., and Olson, J. B. (2012). Disease and nutrient enrichment as potential stressors on the Caribbean sponge Aplysina cauliformis and its bacterial symbionts. Mar. Ecol. Prog. Ser. 456, 101–111. doi: 10.3354/meps09716

CrossRef Full Text

Hall, B. G. (2013). Building phylogenetic trees from molecular data with MEGA. Mol. Biol. Evol. 30, 1229–1235. doi: 10.1093/molbev/mst012

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hoffmann, F., Radax, R., Woebken, D., Holtappels, M., Lavik, G., Rapp, H. T., et al. (2009). Complex nitrogen cycling in the sponge Geodia barretti. Environ. Microbiol. 11, 2228–2243. doi: 10.1111/j.1462-2920.2009.01944.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Honda, D., Yokota, A., and Sugiyama, J. (1999). Detection of seven major evolutionary lineages in cyanobacteria based on the 16S rRNA gene sequence analysis with new sequences of five marine Synechococcus strains. J. Mol. Evol. 48, 723–739. doi: 10.1007/PL00006517

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Koop, K., Booth, D., Broadbent, A., Brodie, J., Bucher, D., Capone, D., et al. (2001). ENCORE: the effect of nutrient enrichment on coral reefs. Synthesis of results and conclusions. Mar. Pollut. Bull. 42, 91–120. doi: 10.1016/S0025-326X(00)00181-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kroon, F. J., Kuhnert, P. M., Henderson, B. L., Wilkinson, S. N., Kinsey-Henderson, A., Abbott, B., et al. (2012). River loads of suspended solids, nitrogen, phosphorus and herbicides delivered to the Great Barrier Reef lagoon. Mar. Pollut. Bull. 65, 167–181. doi: 10.1016/j.marpolbul.2011.10.018

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lane, D. J. (1991). “16S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackerbrandt and M. Goodfellow (New York, NY: John Wiley and Sons, Inc.), 115–148.

Luter, H. M., Whalan, S., and Webster, N. S. (2010). Exploring the role of microorganisms in the disease-like syndrome affecting the sponge Ianthella basta. Appl. Environ. Microbiol. 76, 5736–5744. doi: 10.1128/AEM.00653-10

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Luter, H. M., Whalan, S., and Webster, N. S. (2011). The marine sponge Ianthella basta can recover from stress-induced tissue regression. Hydrobiologia 687, 227–235. doi: 10.1007/s10750-011-0887-x

CrossRef Full Text

Luter, H. M., Whalan, S., and Webster, N. S. (2012). Thermal and sedimentation stress are unlikely causes of brown spot syndrome in the coral reef sponge, Ianthella basta. PLoS ONE 7:e39779. doi: 10.1371/journal.pone.0039779

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Marchesi, J. R., Sata, T., Weightman, A. J., Martine, T. A., Fry, J. C., Hiom, S. J., et al. (1998). Design and evolution of useful bacterium-specific PCR primers that amplify genes coding for bacterial 16S rRNA. Appl. Environ. Microbiol. 64, 795–799.

Pubmed Abstract | Pubmed Full Text

McKinnon, A. D., Smit, N., Townsend, S., and Duggan, S. (2006). “Darwin Harbour: water quality and ecosystem structure on a tropical Harbour in the early stages of urban development,” in The Environment in Asia Pacific Harbours, ed. E. Wolanski ( Netherlands: Springer), 433–459.

Mohamed, N. M., Saito, K., Tal, Y., and Hill, R. T. (2010). Diversity of aerobic and anaerobic ammonia-oxidizing bacteria in marine sponges. ISME J. 4, 38–48. doi: 10.1038/ismej.2009.84

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Off, S., Alawi, M., and Spieck, E. (2010). Enrichment and physiological characterization of a novel Nitrospira-like bacterium obtained from a marine sponge. Appl. Environ. Microbiol. 76, 4640–4646. doi: 10.1128/AEM.00320-10

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Parsons, T. R., Maita, Y., and Lalli, C. (1984). A Manual of Chemical and Biological Methods for Seawater Analysis. New York: Pergamon Press.

Pester, M., Schleper, C., and Wagner, M. (2011). The Thaumarchaeota: an emerging view of their phylogeny and ecophysiology. Curr. Opin. Microbiol. 14, 300–306. doi: 10.1016/j.mib.2011.04.007

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Power and Water Corporation. (2006). Wastewater Treatment, Reuse and Discharge Report. Darwin, NT: Power and Water Corporation.

Pruesse, E., Peplies, J., and Glöckner, F. O. (2012). SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 28, 1823–1829. doi: 10.1093/bioinformatics/bts252

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Quince, C., Lanzen, A., and Davenport, R. J. (2011). Removing noise from pyrosequenced amplicons. BMC Bioinformatics 12:38. doi: 10.1186/1471-2105-12-38

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Reiswig, H. M. (1971). Particle feeding in natural populations of three marine Demosponges. Biol. Bull. 141, 568–591. doi: 10.2307/1540270

CrossRef Full Text

Schläppy, M. L., Schöttner, S. I., Lavik, G., Kuypers, M. M. M., Beer, D., and Hoffmann, F. (2009). Evidence of nitrification and denitrification in high and low microbial abundance sponges. Mar. Biol. 157, 593–602. doi: 10.1007/s00227-009-1344-5

CrossRef Full Text

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister. E. B., et al. (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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Schmitt, S., Hentschel, U., and Taylor, M. W. (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

CrossRef Full Text

Simister, R. L., Deines, P., Botté, E. S., Webster, N. S., and Taylor, M. W. (2012a). 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

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Simister, R., Taylor, M. W., Tsai, P., and Webster, N. (2012b). Sponge-microbe associations survive high nutrients and temperatures. PLoS ONE 7:e52220. doi: 10.1371/journal.pone.0052220

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Smith, J., Burford, M. A., Revill, A. T., Haese, R. R., and Fortune, J. (2011). Effect of nutrient loading on biogeochemical processes in tropical tidal creeks. Biogeochemistry 108, 359–380. doi: 10.1007/s10533-011-9605-z

CrossRef Full Text

Smith, V. S., Kimmerer, W. J., Laws, E. A., Brock, R. E., and Walsh, T. W. (1981). Kaneohe Bay sewage diversion experiment: perspectives on ecosystem responses to nutritional perturbation. Pac. Sci. 35, 279–340.

Southwell, M. W., Weisz, J. B., Martens, C. S., and Lindquist, N. (2008). In situ fluxes of dissolved inorganic nitrogen from the sponge community on Conch Reef, Key Largo, Florida. Limnol. Oceanogr. 53, 986–996. doi: 10.4319/lo.2008.53.3.0986

CrossRef Full Text

Stahl, D. A., and Amann, R. I. (1991). “Development and application of nucleic acid probes in bacterial systematics,” in Nucleic Acid Techniques in Bacterial Systematics, eds E. Stackebrandt and M. Goodfellow (Chichester: John Wiley and Sons), 205–248.

Steindler, L., Huchon, D., Avni, A., and Ilan, M. (2005). 16S rRNA phylogeny of sponge-associated Cyanobacteria. Appl. Environ. Microbiol. 71, 4127–4131. doi: 10.1128/AEM.71.7.4127-4131.2005

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tamura, K., Strecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Taylor, M. W., Radax, R., Steger, D., and Wagner, M. (2007). Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol. Mol. Biol. Rev. 71, 295–347. doi: 10.1128/MMBR.00040-06

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Taylor, M. W., Schupp, P. J., Dahllöf, I., Kjelleberg, S., and Steinberg, P. D. (2004). Host specificity in marine sponge-associated bacteria, and potential implications for marine microbial diversity. Environ. Microbiol. 6, 121–130. doi: 10.1046/j.1462-2920.2003.00545.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Taylor, M. W., Schupp, P. J., de Nys, R., Kjelleberg, S., and Steinberg, P. D. (2005). Biogeography of bacteria associated with the marine sponge Cymbastela concentrica. Environ. Microbiol. 7, 419–433. doi: 10.1111/j.1462-2920.2004.00711.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Turque, A. S., Batista, D., Silveira, C. B., Cardoso, A. M., Vieira, R. P., Moraes, F. C., et al. (2010). Environmental shaping of sponge associated archaeal communities. PLoS ONE 5:e15774. doi: 10.1371/journal.pone.0015774

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Voss, J. D., and Richardson, L. L. (2006). Nutrient enrichment enhances black band disease progression in corals. Coral Reefs 25, 569–576. doi: 10.1007/s00338-006-0131-8

CrossRef Full Text

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Webster, N. S., Luter, H. M., Soo, R. M., Botté, E. S., Simister, R. L., Abdo, D., et al. (2013). Same, same but different: symbiotic bacterial associations in GBR sponges. Front. Microbiol. 3:444. doi: 10.3389/fmicb.2012.00444

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Webster, N. S., and Taylor, M. W. (2012). Marine sponges and their microbial symbionts: love and other relationships. Environ. Microbiol. 14, 335–346. doi: 10.1111/j.1462-2920.2011.02460.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Webster, N. S., Taylor, M. W., Behnam, F., Lücker, S., Rattei, T., Whalan, S., et al. (2010). Deep sequencing reveals exceptional diversity and modes of transmission for bacterial sponge symbionts. Environ. Microbiol. 12, 2070–2082. doi: 10.1111/j.1462-2920.2009.02065.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Weisz, J. B., Lindquist, N., and Martens, C. S. (2008). Do associated microbial abundances impact marine demosponge pumping rates and tissue densities? Oecologia 155, 367–376. doi: 10.1007/s00442-007-0910-0

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wilkinson, C. R., and Cheshire, A. C. (1989). Patterns in the distribution of sponge populations across the central Great Barrier Reef. Coral Reefs 8, 127–134. doi: 10.1007/BF00338268

CrossRef Full Text

Wilkinson, C. R., and Fay, P. (1979). Nitrogen fixation in coral reef sponges with symbiotic cyanobacteria. Nature 279, 527–529. doi: 10.1038/279527a0

CrossRef Full Text