# Diversity of Bacterial Photosymbionts inLubomirskiidae Sponges from Lake Baikal

Authors: Nina V. Kulakova; Natalia N. Denikina; Sergei I. Belikov
Journal: International Journal of Biodiversity
 (2014)

---

## Abstract

Sponges are permanent benthos residents which establish complex associations with a variety of microorganisms that raise interest in the nature of sponge-symbionts interactions. A molecular approach, based on the identification of the 16S rRNA and ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit genes, was applied to investigate diversity and phylogeny of bacterial phototrophs associated with four species ofLubomirskiidae in Lake Baikal. The phylogeny inferred from both genes showed three main clusters of Synechococcus associated with Baikalian sponges. One of the clusters belonged to the cosmopolitan Synechococcus rubescens group and the two other were not related to any of the assigned phylogenetic groups but placed as sister clusters to S. rubescens. These results expanded the understanding of freshwater sponge-associated photoautotroph diversity and suggested that the three phylogenetic groups of Synechococcus are common photosynthetic symbionts in Lubomirskiidae sponges.


---

## Body

## 1. Introduction

Sponges are an important component of the marine and freshwater benthos ecosystems that establish associations with a great diversity of unicellular and multicellular organisms [1]. At the photosynthetic zone, sponges can benefit from phototrophic symbionts which fix carbon using the Calvin-Benson cycle and provide products of photosynthesis to the host [2–4].Photosynthetic symbionts are prevalent in marine sponges of coastal regions worldwide where they contribute significantly to net primary production [5, 6]. From one-third to more than half of the sponges of tropical and temperate regions harbor a high level of photosynthetic symbionts [7, 8]. In Lake Baikal, sponges are necessary components of the benthos and ubiquitous on rocky grounds in the littoral zone. Sponges from the endemic familyLubomirskiidae are widely distributed in Lake Baikal and often harbor photosynthetic symbionts. From 14 described species ofLubomirskiidae [9], there are three common species among which photosyntheticLubomirskia baicalensis (L. baicalensis) andBaikalospongia bacillifera (B. bacillifera) are widely distributed in the photic zone of Lake Baikal. In contrast to marine sponges, there is not a lot of data on photosynthetic symbionts of freshwater sponges, although associations with unicellular green algae, includingChlorella spp.,Choricystis minor, yellow-green algae, andChloroflexi have been shown in cosmopolitan sponges [4, 10–13] and cyanobacterial sequences detected inL. baicalensis [14]. Nevertheless, the identification and diversity of sponge-associated phototrophs in Lake Baikal are undetermined to date.We analysed 16S rRNA and ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL) genes to identify and assess the diversity of photosynthetic symbionts in Baikalian sponges. Both molecular markers are widely used for analysis of genetic diversity [15–21].


## 2. Materials and Methods

Four species of sponges from the familyLubomirskiidae:L. baicalensis,Lubomirskia abietina (L. abietina),B. bacillifera, andBaikalospongia martinsoni (B. martinsoni) and water samples were collected from the Southern Basin of Lake Baikal at depths of 10–20 m by scuba diving (Table S1, see Supplementary Material available online at http://dx.doi.org/10.1155/2014/152097). Sponge tissue samples (3–5 cm3) were rinsed twice in 96% ethanol and stored at 4°C until DNA extraction. Species identification was based on external morphological characteristics and morphology of spicules according to the guides of Rezvoi (1936) and Efremova (2001) [22, 23].Chlorophyll concentration was measured after 96% ethanol extraction as described by Bergmann and Peters [24] and Webb et al. [25] from fresh frozen (−20°C) samples (two replicates). Spectral absorbance scans were performed from 300–800 nm using a UV-visible spectrophotometer (Cintra-10E, GBC, Australia).Total DNA was extracted using the RiboSorb kit (AmpliSens, Russia) according to the manufacturer’s protocol. TheRbcL gene was amplified in 15 μL of PCR reaction mix (Screen Mix, Evrogen, Russia) with 10 pmol of each primer, cbbL 595f (5′-GACTTCACCAAAGACGACGA-3′) and 1387r (5′-TCGAACTTGATTTCTTTCCA-3′), as described by Elsaied and Naganuma [16]. Touchdown PCR was done in four repeats and mixed for each tissue sample. The amplification profile consisted of an initial denaturation at 95°C for 3 min followed by 5 touchdown cycles of 94°C for 20 s, 55°C for 20 s, and 72°C for 1 min and then 20 cycles of 94°C for 20 s, 52°C for 20 s, and 72°C for 1 min and 10 cycles of 94°C for 20 s, 50°C for 20 s, and 72°C for 1 min. 35 PCR cycles were followed by a final extension at 72°C for 15 min. The 16S rRNA gene was amplified with primers CYA106Fc (5′ CGGACGGGTGAGTAACGCGTGA 3′) and CYA781R (5′ GACTACWGGGGTATCTAATCCCWTT 3′) [26]. The amplification profile consisted of an initial denaturation at 95°C for 3 min followed by 30 cycles of 94°C for 20 s, 60°C for 20 s, 72°C for 1 min and a final extension at 72°C for 15 min.The PCR products were detected by electrophoresis in 0.8% agarose gel and purified using a DNA purification kit (Cytokine, Russia) and cloned into the pAL-TA vector (Evrogen, Russia). Inserts were detected by amplification with M13 primers and digested with endonucleases HhaI and HaeIII. Clone inserts with different restriction profiles were sequenced using ABI 3130xl Genetic Sequencer (Applied Biosystem, USA) (Table S1). Sequences were checked for chimeras with the Decipher tool [27].The basic local alignment search tool (BLAST,http://blast.ncbi.nlm.nih.gov/) was used to compare sample sequences to closely related sequences from the NCBI database. Pairwise alignment was accomplished using BioEdit alignment editor version 7.0.9.0. [28] and phylogenetic reconstructions were performed in Mr. Bayes 3.2.1 and Mega 5 programs [29, 30].For phylogenetic tree reconstruction maximum parsimony (MP), maximum likelihood (ML), and Bayesian analysis were performed. Kimura 2-parameter model was used for estimating of genetic distances [31]. The first two codon positions were used for analysis ofrbcL dataset (MP, ML). Bootstrap analysis was performed with 1000 replicates. In the Bayesian analysis, the MCMC chain was run for 1,000,000, sampled every 100th step.Sequences were deposited in the GenBank database under the following accession numbers: JX570937–JX570966, KF856235–KF856242, and JX570967–JX571009 and KF856243–KF856253 for the 16S rRNA andrbcL genes, respectively.


## 3. Results

Clone library screening was performed to identify phototrophs associated with Baikalian sponges. The 16S rRNA andrbcL gene fragments were analysed inL. baicalensis,L. abietina,B. bacillifera andB. martinsoni. The chlorophyll level was estimated for three of these species and chlorophyll A was detected inL. baicalensis (122
±
32 μg/g),B. bacillifera (37
±
12 μg/g), andL. abietina(20
±
0.5 μg/g).In total, 71 and 126 inserts were sequenced from the 16S rRNA andrbcL clone libraries, respectively. Identical sequences were removed from further analysis and 38 16S rRNA and 54rbcL gene sequences were deposited in the GenBank database. Sequences obtained from 16S rRNA gene were 95–100% identical, while sequences obtained from therbcL gene were more heterogeneous (83–100%).RbcL sequences with nucleotide identities greater than 90% were combined into operational taxonomic units (OTUs). There were three α-cyanobacteria OTUs resolved from analysis of the 733 bprbcL gene alignment. The first OTU was 98% identical toSynechococcus rubescens (AM701775) while the maximum nucleotide sequence identities from two other OTUs were 87% and 89% similar toCyanobium gracile 6307 (CP003495), and closer relatives could not be found in the GenBank database.The phylogenetic analysis inferred from 16S rRNA andrbcL genes showed three main clusters forSynechococcus which included sequences from all species of sponges and lake water analysed (Figures 1 and 2). Better resolved trees were based onrbcL in comparison with 16S rRNA based phylogeny. The phylogenetic reconstruction inferred from both genes showed that phylogenetic positions for two of three clusters were similar. There was a fully resolved cluster ofS. rubescens that combined highly identical (99%) sequences from Lake Baikal, as well asSynechococcus strains from European lakes. The two other groups were placed as sister clusters toS. rubescens (Figure 1). Clusters BL1 and BS1 included sequences from Lake Baikal only (Figures 1 and 2).Figure 1
Phylogenetic tree of cyanobacteria associated with photosynthetic sponges from Lake Baikal. Consensus tree inferred from 733 nucleotide positions of therbcL gene. Clusters fromLubomirskiidae sponges are colored in grey. Bootstrap support values and posterior probability (%) are shown at the nodes (Maximum Parsimony/Maximum Likelihood/Bayesian analysis). Values below 50% are not shown. Scale bar represents 0.02 substitutions per site.Figure 2
Phylogenetic tree of cyanobacteria associated with photosynthetic sponges from Lake Baikal inferred from 522 nucleotide positions of the 16S rRNA gene. Sequences from this study and reference sequences ofS. rubescens andC. gracile are shown in bold. Bayesian posterior probability (%) is shown at the nodes. Scale bar represents 0.2 substitutions per site. GenBank accession numbers are given in the parentheses. Accession numbers for 6 sequences from sponges and water, Baikal: JX570937, JX570945, JX570948,  KF856235, KF856236, and KF856237; 23 sequences from sponges and water, Baikal: JX570939, JX570957–JX570966, KF856238–KF856242, and uncultured bacterium (Pyrenean, Tibetan, Mongolian lakes): HE857287, HE857263, HM129960, DQ297463, DQ422951, DQ297460, and DQ297461. In the Bayesian analysis, the MCMC chain was run for 2,000,000, sampled every 100th step.The phylogenetic positions of identified clusters were slightly different from the 16S rRNA based phylogeny. Using 16S rRNA data, there was a cluster related to theS. rubescens phylogenetic lineage that included the clade of 23Synechococcus sequences (intragroup identity 98-99%) from Lake Baikal, sequences from lakes in Mongolia, and high altitude oligotrophic lakes of the Tibetan plateau and the Pyrenees, whilerbcL-based phylogeny showed a specific cluster of 25 sequences that was sister toS. rubescens cluster (Figures 1 and 2). The rest of the sequences (16S rRNA tree) formed single branches or were clustered within Baikalian clades of picocyanobacteria from lake water (Figure 2).The onlyrbcL sequence fromB. bacillifera (JX570973) that showed a low-level nucleotide identity (<83%) with other sequences from the GenBank database formed a distinct branch within the β-cyanobacteria (Figure 1).Thus, the majority of sequences from sponges and water belonged toα-cyanobacteria with form IA RuBisCO. Sequences from Baikalian sponges did not form distinct phylogenetic clusters and no major sponge-specific clusters have been found.


## 4. Discussion

The phylogeny inferred from 16S rRNA andrbcL genes generally agrees that three main phylogenetic clusters for chlorophyll IA-containing picocyanobacteria are present in Baikalian sponges. The only one sequence belonging toOscillatoria was detected inB. bacillifera. The finding ofSynechococcus andOscillatoria groups in freshwater sponges is similar to the results obtained from marine sponges [7, 8, 17–19]. In contrast with marine sponges, no specific phylogenetic groups for sponge species or sponges with various levels of chlorophyll A were found, and sequences derived from sponges (as well as from lake water) were dispersed across the clusters described.The different levels of chlorophyll A in sponges are also demonstrated in marine environments, where the major role for sponge-specificSynechococcus in ecology and primary productivity of sponges with high level of chlorophyll A has been described [7]. The highest level of chlorophyll A we found inL. baicalensis may be dependent on better light condition available to branched structures that grow up to one meter high in comparison to fully attached cushion- and crust-like sponges (L. abietina, B. bacillifera, andB. martinsoni), though the presence of eukaryotic photosymbionts could also be significant.TheS. rubescens cluster detected in Baikalian sponges was affiliated with widely distributed Subalpine cluster I (group B) [20] proving high ecological plasticity of this group. The presence of phylogenetic clusters that included only Baikalian sequences, as well as clusters with no attribution to theSynechococcus phylogenetic group described, confirms the results of Pommier et al. [21] suggesting that endemic clusters are formed in the local bacterioplankton community in response to their adaptation to unique ecological conditions.These results give a new insight into biodiversity of phototrophs associated with Baikalian sponges that could be further used for purpose of ecosystem monitoring. This is especially important with regard to increasing anthropogenic pressure as a result of tourism development in the Baikalian region.


## 5. Conclusions

Molecular detection of 16S rRNA andrbcL genes from clone libraries were performed for the identification of phototrophs associated with four species ofLubomirskiidae from Lake Baikal. The phylogeny inferred from both genes showed thatS. rubescens and two specific clusters ofSynechococcus were associated with sponges and also found in lake water. Three major clusters ofSynechococcus that did not form sponge-specific groups were detected in all studied species of Baikalian sponges. These results added to the understanding of symbiotic associations in freshwater photosynthetic sponges and could be further applied to the assessment of ecological impacts on the ecosystems in Lake Baikal.



---

*Source: 152097-2014-11-20.xml*