Symbiotic Archaea in Marine Sponges Show Stability and Host Specificity In
Total Page:16
File Type:pdf, Size:1020Kb
RESEARCH ARTICLE Symbiotic archaea in marine sponges show stability and host specificity in community structure and ammonia oxidation functionality Fan Zhang1, Lucıa Pita2, Patrick M. Erwin3, Summara Abaid1, Susanna Lopez-Legentil 2,3 & Russell T. Hill1 1Institute of Marine and Environmental Technology (IMET), University of Maryland Center for Environmental Science (UMCES), Baltimore, MD, USA; 2Department of Animal Biology, Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona (UB), Barcelona, Spain; and 3Department of Biology & Marine Biology, Center for Marine Science, University of North Carolina Wilmington, Wilmington, NC, USA Correspondence: Russell T. Hill, 701 E Pratt Abstract St, Baltimore, MD 21202, USA. Tel.: +1 410 234 8802; Archaea associated with marine sponges are active and influence the nitrogen fax: +1 410 234 8818; metabolism of sponges. However, we know little about their occurrence, speci- e-mail: [email protected] ficity, and persistence. We aimed to elucidate the relative importance of host specificity and biogeographic background in shaping the symbiotic archaeal Received 4 April 2014; revised 2 September communities. We investigated these communities in sympatric sponges from 2014; accepted 10 September 2014. Final the Mediterranean (Ircinia fasciculata and Ircinia oros, sampled in summer and version published online 20 October 2014. winter) and from the Caribbean (Ircinia strobilina and Mycale laxissima). PCR DOI: 10.1111/1574-6941.12427 cloning and sequencing of archaeal 16S rRNA and amoA genes showed that the archaeal community composition and structure were different from that in Editor: Gary King seawater and varied among sponge species. We found that the communities were dominated by ammonia-oxidizing archaea closely related to Nitrosopumi- Keywords lus. The community in M. laxissima differed from that in Ircinia spp., includ- Porifera; ammonia-oxidizing archaea; 16S ing the sympatric sponge I. strobilina; yet, geographical clusters within Ircinia rRNA gene; microbial biogeography; spp. were observed. Whereas archaeal phylotypes in Ircinia spp. were persistent symbiosis; archaeal communities. and belong to ‘sponge-enriched’ clusters, archaea in M. laxissima were closely related with those from diverse habitats (i.e. seawater and sediments). For all four sponge species, the expression of the archaeal amoA gene was confirmed. Our results indicate that host-specific processes, such as host ecological strategy and evolutionary history, control the sponge–archaeal communities. previously known as Marine Group 1 Crenarchaeota (Bro- Introduction chier-Armanet et al., 2008; Pester et al., 2011). Based on Marine sponges are one of the most significant groups in genomic information and physiology experiments, mem- benthic ecosystems in terms of diversity, abundance, and bers of this phylum of archaea are chemolithoautotrophic function (Diaz & Rutzler, 2001; Ribes et al., 2012; De ammonia oxidizers; that is, they can convert ammonium Goeij et al., 2013); notably, they are key players in the to nitrite and fix CO2 using the energy obtained from recycling of nutrients such as nitrogen (Jimenez & Ribes, ammonia oxidation. Archaeal members of sponge symbi- 2007; Southwell et al., 2008). Interestingly, their contribu- otic microbial communities may play an important role tion to nutrient cycles in the ecosystem relies on the in the nitrogen cycle of the ocean. complex microbiota associated with them, including bac- Despite the interesting metabolism of these organisms, teria and archaea (Mohamed et al. 2010; Lopez-Legentil the symbiotic archaeal communities in sponges have MICROBIOLOGY ECOLOGY MICROBIOLOGY et al., 2010; Radax et al., 2012). Archaea have been received less attention than bacterial communities and we detected in sponges from different oceans (Webster et al., know little about the specificity, persistence, or resilience 2001; Margot et al., 2002; Lee et al., 2003) and belong of archaea in sponges. With one exception (Lee et al., mostly to the newly defined kingdom Thaumarchaeota, 2011), sponge symbiotic archaeal communities have been FEMS Microbiol Ecol 90 (2014) 699–707 ª 2014 Federation of European Microbiological Societies. Published by John Wiley & Sons Ltd. All rights reserved 700 F. Zhang et al. found in low diversity, are metabolically active, and may collected, maintained at the ambient temperature, and reach high densities in the mesohyl in some sponge spe- transported rapidly (within 1 h) back to a shore labora- cies (Webster et al., 2001; Bayer et al., 2008; Radax et al., tory for processing. Seawater was drained and the sponges 2012). The first study of archaea in sponges found a sole were rinsed three times with sterile artificial seawater. Tis- archaeal phylotype persistently associated with Axinella sue samples (1-cm cubes) were sterilely excised from mexicana that was named Cenarchaeum symbiosum (Pres- sponges. Samples for DNA and RNA extraction were pre- ton et al., 1996); and, since then, closely related phylo- served in RNAlater stabilization solution (Qiagen, Valen- types have been found in other sponge species (Hentschel cia, CA) on board, transferred to a À20 °C freezer on site et al., 2006). The presence of archaea in sponge larvae then kept in a À80 °C freezer for long-term storage. suggests a tight link between archaeal symbionts and the sponge host phylogeny (Schmitt et al., 2008; Steger et al., Genomic DNA/RNA extraction 2008). Conversely, environmental factors (e.g. seasonal changes in seawater conditions) or host biogeography Total DNA and RNA from sponges and the seawater may also affect the structure of archaeal symbionts in samples from Key Largo were extracted using a TissueLy- marine sponges (Turque et al., 2010). ser System (Qiagen) and AllPrep DNA/RNA Mini Kit Our study aims to elucidate the relative importance of (Qiagen) with RNAase-free DNase treatment steps (Qia- host species-specific factors and biogeographic back- gen) for RNA samples following the manufacturer’s pro- ground in shaping sponge-associated archaeal communi- tocol. Reverse transcription of RNA was performed using ties. To achieve our objective, we investigated the RevertAid Reverse Transcriptase (Thermo Scientific, Wal- archaeal communities in congeneric sponges and distantly tham, MA) with random primers and following manufac- related sympatric sponges. Two Mediterranean species, turer’s protocol. RNA samples without the RT step were Ircinia fasciculata and Ircinia oros, and two Caribbean included as PCR template to check for residual DNA in species, Ircinia strobilina and Mycale laxissima, were tar- the RNA samples. geted. Mediterranean species were sampled in summer and winter, as this region shows a marked seasonality in Archaeal 16S rRNA gene clone library seawater conditions (Erwin et al., 2012b). We used clon- construction ing and Sanger sequencing techniques to study archaeal 16S rRNA and amoA genes as phylogenetic markers as We constructed archaeal 16S rRNA gene clone libraries well as amoA transcripts to assess the active ammonia- from total DNA for each sponge specimen and seawater oxidizing archaea communities. The archaeal amoA gene samples from Key Largo. Partial archaeal 16S rRNA gene encodes the a-subunit of ammonia mono-oxygenase, sequences (c. 950 bp) were amplified using the archaea- which is essential in the nitrification process and is well specific primer set 21F (DeLong, 1992): (50-TTC CGG conserved. The archaeal communities from seawater were TTG ATC CYG CCG GA-30) and 915R (Siboni et al., also included in the analysis. 2008): (50-GTG CTC CCC CGC CAA TCC-30) at 0.2 lM each. 50 lL of PCR mix included: 5 lLof109 high- fidelity PCR buffer (Invitrogen Life Technologies, Carls- Material and methods bad, CA), 2 lL MgCl2,1lL of a mix of deoxynucleoside â triphosphates (dNTP) at 0.2 mM, 0.2 lL of Platinum Sampling strategy Taq DNA Polymerase (Invitrogen), and 1 lL of each pri- All sponge and seawater samples were collected by mer. Thermocycler parameters were set to 95 °C for SCUBA diving. Tissue samples from three I. fasciculata 5 min; 30 cycles of 94 °C for 45 s, 56 °C for 60 s, and individuals and three I. oros individuals were collected at 72 °C for 60 s; then, a final extension at 72 °C for Tossa de Mar (Girona), NE Spain, NW Mediterranean 15 min. Amplification products were analyzed by electro- Sea (41°43.230N, 2°56.450E) during the summer (Septem- phoresis in 1.0% (w/v) agarose gels in 19 TAE buffer. ber 2012) and winter (March 2013) seasons. Three PCR products were ligated into PCR-XL-TOPO vectors M. laxissima individuals and three I. strobilina individuals and transformed into ONEShot TOP10 chemically com- were collected from Conch Reef, Key Largo, FL, NE petent Escherichia coli cells using the TOPO XL PCR Caribbean (24°57.110N, 80°27.570W), in July 2011. Three Cloning Kit (Invitrogen). Plasmid DNA was isolated from seawater samples from Key Largo (3 L for each sample) individual clones, purified using Mini prep spin kit were collected simultaneously in close proximity (1 m) to (Qiagen), and sequenced using an ABI PRISM 3130XL sampled sponges and filtered through 0.22 lm Sterivex genetic analyzer (Applied Biosystems, Foster City, CA). filter units (Millipore, Billerica, MA). In all cases, intact Sequences were screened for chimeras in MOTHUR software sponges were held in the seawater in which they were (Schloss et al., 2009). ª 2014 Federation of European Microbiological Societies. FEMS Microbiol Ecol 90 (2014) 699–707 Published by John Wiley & Sons Ltd. All rights reserved Stability and specificity of sponge-associated archaea 701 improved manually using the ARB editor. The aligned ar- Diversity and structure of 16S rRNA gene clone chaeal 16S rRNA gene (699 bp) or amoA gene (489 bp) libraries was imported into PHYML 3.1 software package to construct All archaeal 16S rRNA gene sequences were aligned in a tree based on maximum likelihood method (Jukes–Can- ARB (Ludwig et al., 2004) and ascribed to 97% opera- tor correction). The robustness of the resulting tree topolo- tional taxonomic units (OTUs) using the MOTHUR soft- gies was evaluated by 100 bootstrap replicates for archaeal ware package (Schloss et al., 2009).