Skip to main content

Telonemia-specific environmental 18S rDNA PCR reveals unknown diversity and multiple marine-freshwater colonizations

Abstract

Background

Recent surveys of eukaryote 18S rDNA diversity in marine habitats have uncovered worldwide distribution of the heterotrophic eukaryote phylum Telonemia. Here we investigate the diversity and geographic distribution of Telonemia sequences by in-depth sequencing of several new 18S rDNA clone libraries from both marine and freshwater sites by using a Telonemia-specific PCR strategy.

Results

In contrast to earlier studies that have employed eukaryote-wide PCR design, we identified a large and unknown diversity of phylotypes and the first rigorous evidence for several freshwater species, altogether comprising 91 unique sequences. Phylogenies of these and publicly available sequences showed 20 statistically supported sub-clades as well as several solitary phylotypes with no clear phylogenetic affiliation. Most of these sub-clades were composed of phylotypes from different geographic regions.

Conclusions

By using specific PCR primers we reveal a much larger diversity of Telonemia from environmental samples than previously uncovered by eukaryote-wide primers. The new data substantially diminish the geographic structuring of clades identified in earlier studies. Nevertheless, since these clades comprise several distinct phylotypes we cannot exclude endemicity at species level. We identified two freshwater clades and a few solitary phylotypes, implying that Telonemia have colonized freshwater habitats and adapted to the different environmental and ecological conditions at independent occasions.

Background

Microorganisms usually exist in populations of huge sizes and are highly prone to long-distance dispersal by vectors such as wind, water, animals and humans [1–5]. Obvious barriers to dispersal are lacking, especially in the marine habitat [4–8]. The ubiquitous dispersal of microorganisms has been a prevalent view since the turn of the last century, summarized in the statement "everything is everywhere, but, the environment selects" [9, 10]. This view has been challenged however, by investigations of environmental DNA clone libraries as a large number of cryptic species and restricted biogeographies have been revealed [11–20]. High levels of genetic diversity have been found, even within the slowly evolving small ribosomal subunit gene [21, 22]. However, as more localities are being investigated and the variety of sampling strategies increase, the geographic ranges of many microorganisms have been expanded, showing that under-sampling of the diversity can cause a false impression of endemism [see [4, 5]]. Some surveys have therefore interpreted the diversity as consistent with the "Moderate Endemicity Model" (MEM), which states that some microbial lineages do in fact have a global distribution, but that there also exists species with restricted dispersal and local adaptations [4, 23–25].

The vast majority of 18S rDNA environmental surveys conducted so far have involved universal primers designed to capture the broadest diversity of eukaryotes possible. However, much diversity is most likely overlooked by applying only a single pair of universal primers [26–28]. This could be due to a number of reasons, e.g. the primers are less suitable for some groups of organisms, there are great variations in rDNA copy number, as well as bias introduced in the PCR reaction. One of the most efficient approaches to address these problems has been to apply a group-specific PCR strategy with primers targeting the particular taxonomic group of interest [29–32]. These studies have shown that the use of such primers is detecting far more diversity than the universal approach.

Telonemia is one of the groups of unicellular eukaryotes that are frequently detected in marine 18S rDNA environmental clone libraries, but usually represents only a relatively small part of the total diversity [11, 33–36]. So far, only two species have been described on the basis of morphology, Telonema subtilis Griessmann and T. antarcticum Thomsen in Klaveness et al. [20, 37], but the size range of the identified species is large (3.5 - 15 μm long and 4-20 μm wide). It was recently discovered by Shalchian-Tabrizi et al. [36] that the 18S rDNA sequences formed two major groups, Group 1 and 2, including T. subtilis and T. antarcticum respectively, and that these were further sub-divided into several statistically supported clades of sequences with restricted geographic distribution.

Species of Telonemia are heterotrophic predators, feeding on a wide range of bacteria and pico- to nano-sized phytoplankton. They are globally distributed in marine waters and are frequently encountered in environmental clone libraries e.g. [34, 38]. Telonemia are present throughout the year and are considered to play an important ecological role, as they have been found to dominate the heterotrophic protist community on certain occasions [37]. Very little is known about the life cycle and reproduction of Telonemia. Asexual reproduction occurs by cell division and the possible presence of cysts has been indicated by Vørs [39], but this is yet to be verified.

Telonemia has also been reported from fresh water habitats. Tong et al. [40] identified a freshwater T. subtilis in an Antarctic lake, Sombre Lake, but it is unclear if this specimen is truly freshwater as the lake has been classified as maritime [41]. A survey of Finnish lakes recorded Telonema sp. on a number of occasions (Liisa Lepistö, personal communication). The ability to survive under low salinity conditions have also been shown in culture experiments done on T. subtilis from Norwegian coastal waters [42]. Although Telonemia has been observed at several occasions in freshwater, only a few 18S rDNA sequences appear to be related to the group [43]. Therefore, it is still unclear how large the diversity of Telonemia might be in these habitats and what phylogenetic relationship they have to marine species. It is also unclear whether Telonemia have colonized these habitats at one or several independent occasions, and if both the two major groups related to T. subtilis and T. antarcticum have been successfully established in freshwater.

Here, we have designed Telonemia-specific 18S rDNA primers in order to investigate (i) whether group-specific environmental PCR will uncover a larger diversity of Telonemia than so far uncovered by universal primers, (ii) whether increased taxon sampling will affect the geographic structuring observed for many clades of marine Telonemia [36], and (iii) to examine whether one or several species exist in freshwater, and whether both Group 1 and 2 comprise species from freshwater. We address these questions by sequencing clone libraries from 4 marine and 3 freshwater localities, as well as including all available Telonemia sequences already published. We identify a large and unknown diversity of phylogenetically distinct sequences (phylotypes) in marine waters and provide the first rigorous evidence for the presence of several freshwater species related to both T. subtilis and T. antarcticum resulting from independent colonisations of freshwater.

Results and discussion

Large cryptic diversity of Telonemia in marine habitats

Despite the huge amount of environmental 18S rDNA sequences from numerous diversity studies available in public databases, only 33 were found to belong to Telonemia in Shalchian-Tabrizi et al. [36], all amplified by universal eukaryotic primers. These sequences were divided into two main groups, Group 1 and Group 2, including T. subtilis and T. antarcticum respectively [36]. Within these groups, twelve distinct sub-groups or independent phylotypes were identified, each possibly representing several species or populations. The majority of these clades were composed of sequences from single localities, suggesting a considerable geographic structuring of Telonemia [36].

By using group-specific primers we generated 145 18S rDNA sequences affiliated to Telonemia. No sequences from other eukaryote groups were generated. The evolutionary origin of these sequences was inferred by phylogenetic analyses of an alignment containing a broad diversity of eukaryotic lineages (alignment 1) that included our new data and all putative Telonemia sequences downloaded from GenBank (result not shown). Hence, the group specific PCR strategy for Telonemia clearly improves our knowledge about the diversity of the group.

To better resolve the phylogeny of the Telonemia sequences we removed all other eukaryote groups (except haptophytes, cryptophytes and katablepharids used as outgroups) that allowed for inclusion of more unambiguously aligned nucleotide characters (i.e. alignment 2). This phylogeny recovered Group 1 and 2, here renamed to TEL 1 and TEL 2 respectively, with high support (1.00 posterior probability (pp) and >99% bootstrap support (%); Figure 1). Furthermore at least 20 sub-groups (1a-1d and 2a-2p in Figure 1) were supported with substantial statistical support. Several of these groups could perhaps be even further subdivided, based on the internal support values (e.g. groups 1b and 2i) but are here treated as single groups for simplicity. The naming of the groups follows that of Shalchian-Tabrizi et al. [36] and has been extended here to include the new sub-groups.

Figure 1
figure 1

Bayesian phylogeny showing the relationship of the Telonemia 18S rDNA sequences. Numbers at the nodes represent Bayesian and Maximum Likelihood support values respectively. Names in brackets indicate sub-groups recognized in [36] that are referred to in the text. Only values above 50/0.70 are shown and thick branches indicate full statistical support (100/1.00). Blue lines show freshwater sequences and dashed blue lines indicate possible freshwater origin. An asterisk (*) indicates that branch length has been cut in half.

As previously recognized, TEL 1 contained fewer clades than TEL 2 and is here divided into 4 sub-groups. Earlier, TEL 1 was found to be restricted to the English Channel and Danish waters [36], but is here extended to include phylotypes from the Antarctic and Arctic regions as well as warmer waters such as the Indian Ocean, Hawaiian waters, the Mediterranean Sea and Chinese waters. In TEL 2 we could now identify 16 sub-groups sampled from worldwide locations such as the Arctic and Antarctic oceans, the English Channel, Danish and German waters, the Indian Ocean, Sargasso Sea, Mediterranean Sea and Hawaiian waters.

Implications on the geographic structuring of Telonemia

The geographic structuring shown by Shalchian-Tabrizi et al. [36] is here diminished by the addition of more environmental sequences (Figure 1). Several of the sub-groups previously found to have restricted geographic ranges now includes sequences from new locations. For instance the sub-groups 2m and 2o (earlier 2c and 2b; Figure 1), previously found to be restricted to the Arctic Ocean, are now extended to the Indian Ocean, Hawaiian waters and the Mediterranean Sea. The sub-group 2k (earlier 2g), which was previously restricted to the English Channel, Oslo Fjord and Helgoland (i.e. southern parts of the North Sea/Skagerak), now includes sequences from the Mediterranean Sea as well. Additionally, most of the sub-groups new to this study have widespread distributions; e.g. sub-group 1b is composed of sequences from the Indian Ocean, the Mediterranean Sea and the Pacific Ocean, while sub-group 2l is composed of sequences from as distant locations as the Arctic Ocean, Western Pacific, Mediterranean Sea and the Indian Ocean. Although the majority of the subgroups show little geographic structuring, the high diversity uncovered here implies that geographical isolation has existed at some point. The combination of high diversity and low geographic structuring show that subsequent dispersal rates have been higher than speciation rates over the history of Telonemia.

The existence of endemicity cannot be completely excluded however. One important reason is that each clade may represent higher order taxonomic units, like genera or families, and each phylotype can in principle represent separate species (or even several species as 18S rDNA may be too conserved to demarcate species boundaries [21, 25, 44]). Hence, the widespread geography of the subgroups may be hiding endemicity at strain or species level; in fact we could not identify the same phylotypes from different localities. Sampling of DNA from more sites and a larger variety of marine habitats, as well as the use of faster evolving genetic markers, such as the internal transcribed spacer (ITS) of the ribosomal operon, would be necessary to resolve this question. On the other hand, any putative geographic restriction of species or groups should be interpreted with caution because endemicity in general is difficult to prove, as there will always be a possibility of undersampling and absence of species at times of sampling due to seasonal variations.

There are a few tendencies of endemicity or restricted biogeography in the inferred phylogeny (Figure 1) that deserve to be noted here and should be subject for future studies of Telonemia. In particular, a striking pattern is seen for sequences from Antarctic and Arctic regions clustering into sub-group 1a, which opens up the possibility for a bi-polar or anti-tropical distribution. If further diversity studies confirm this pattern, it would be congruent with geographic distribution of dinoflagellates and foraminiferans [45, 46]. Three other clades also appear to be endemic; the clade 2i from the Sargasso Sea and 2h and 2f, are only composed of Indian Ocean and the Norwegian Framvaren Fjord sequences respectively (Figure 1). In addition there is a large assembly of sequences from the Svalbard region that could indicate the presence of a Norwegian-Barents Sea population, but this assembly is only moderately supported (Figure 1).

Cryptic diversity of Telonemia in freshwater

In order to investigate the putative existence of Telonemia in freshwater we had to use a nested PCR amplification strategy. This could explain why so little sequence data from Telonemia in freshwater has been generated previously and confirm visual observations that freshwater Telonemia exists only in minute quantities (L. Lepistö unpublished). The sequences obtained from the three different Norwegian freshwater lakes, Lake Lutvann, Lake Sværsvann and Lake Pollen, together with a few publicly available freshwater environmental sequences, formed three clades (1d, 2e and 2p) and two single phylotypes with representatives in both TEL 1 and TEL 2 (Figure 1). In Lake Lutvann we sampled both the sediment and the water column. Strikingly, these sequences formed two distantly related habitat-specific clades, in which all the benthic sequences clustered into one group (1d) and the pelagic sequences into another (2e), highlighting a vertical stratification of phylotypes or populations within this lake at the time of sampling (Figure 1). Sub-group 2e was in addition composed of sequences from the pelagic zone of the two other Norwegian lakes as well as three other freshwater sequences from Svalbard and France.

A few other phylotypes in TEL 1 may represent additional successful transitions from marine to freshwater lakes. One sequence (DGGE band 20) is sampled from a hyperhaline lake in Chile, Lake Tebenquiche that is situated in the Andes at 2500 m.a.s.l. The lake is classified as hyperhaline but has extreme variations in salinity, ranging from 1% to 30% [47]; hence the potential Telonemia species from this lake could be adapted to any of these salinity conditions or could simply be a marine species that have dispersed into the lake. Another sequence (B-2-8), is sampled from the Bayelva River in Svalbard, which is composed of glacial melt water as well as water from nearby freshwater lakes [48], and discharges into the Kings Bay delta in Spitsbergen. Interestingly the two sequences (SS1E0120 and SS1E0155) that have been obtained from the Kings Bay is branching off close to the sequence from the Bayela River in Figure 1, and may therefore actually have freshwater origin.

A third cluster of freshwater sequences (2p), entirely composed of sequences sampled from a glacier in Svalbard, belonged to TEL 2. This cluster was distantly related to the other freshwater group (2e) and was embedded in a large assembly of Arctic and Antarctic sequences, although this relationship was weakly supported (Figure 1). T. subtilis is commonly observed inhabiting the sea-ice in the Baltic Sea [49] and it is therefore possible that these sequences originate from a marine species transported onto the glacier from marine waters by aerosols or other vectors. On the other hand, if these represent an actual freshwater species this would be a second freshwater species within TEL 2, distantly related to the Bayelva River sequences. It remains to be verified that these are actually living cells and whether these have been transported from freshwater sources or dispersed on to the glacier from marine habitats via aerosols or other vectors. So far, we have not detected sequences from the marine samples that are identical to these glacier phylotypes, which could indicate such freshwater dispersal, but as only few samples have been made in these areas we cannot exclude this possibility.

Few marine-freshwater cross-colonizations

In Figure 1 the freshwater sequences form distinct clusters and phylotypes, suggesting the existence of several different freshwater species. These are placed within both TEL 1 and TEL 2, demonstrating that relatively distantly related species of Telonemia exists in freshwater. This diversity is detected even with a very limited number of samples; we therefore expect future surveys of other types of freshwaters at other continents to uncover an even larger diversity. The clustering pattern of the Telonemia sequences is in accordance with recent studies of other protist groups showing that freshwater species form distinct clades in phylogenetic trees, i.e. they are more closely related to each other than to marine species [reviewed in [50]]. Such clustering pattern of freshwater phylotypes has in these studies been interpreted as successful marine-freshwater transitions. These transitions have often been ancient and rare events, resulting in most of the extant species being restricted to either of the two habitats: e.g. in bodonids [51], goniomonas [52], cryptomonads [53], dinoflagellates [54] and Perkinsea [55]. If further examinations of freshwater with the use of Telonemia-specific PCR approaches confirms the clustering pattern shown here (see Figure 1), it would imply that the biogeophysical differences between marine and fresh waters constitutes a significant ecological barrier for dispersal of Telonemia that affects diversification of the lineage. In addition to more samples from various freshwater localities, a better identification of the species is necessary in order to further reveal the species diversity. For instance, analyzing RNA to confirm that the species are alive and metabolize in the habitat, and fluorescence in situ hybridization (FISH) would be helpful in relating the sequences to the actual cells in the habitats, and to better understand whether the dispersal of species between different and similar habitats take place in the form of spores or as active cells.

Each of the freshwater clades in our tree are habitat-specific in that they only contain phylotypes from either sediment (clade 1d), pelagic (clade 2e), and potentially also glacier (clade 2p) and hyperhaline habitats, implying that each of these habitats have possibly been colonized independently by marine species and adapted to different environmental and ecological conditions. Interestingly, this clustering pattern indicate the existence of ecological barriers also between freshwaters habitats, but as this study primarily has focused on revealing the existence of Telonemia in freshwater, the geographic distribution of the various strains and species should be addressed by much more extensive sampling and adequate molecular methods.

Conclusions

Here we have applied a group-specific PCR approach to better understand the diversity of Telonemia and to investigate whether the geographic structuring observed in earlier studies has been affected by undersampling. Our results show that the use of group-specific primers will uncover a much larger diversity from environmental samples compared to eukaryote-wide primers. The Telonemia-specific primers and the PCR protocol presented here were highly specific for the Telonemia group as no sequences from other eukaryote groups were identified in our sequence libraries. Further, the geographic structuring of marine groups found in earlier analyses is clearly diminished by the addition of the newly generated sequences, showing that undersampling of the diversity may lead to a false impression of endemicity. However, as only two species of Telonemia are defined on basis of morphology, it is not clear what taxonomic units the identified clades represent. Most likely each of these sub-groups are composed of many distinct species, as they comprise phylotypes with different 18S rDNA sequences. If each phylotype is representing separate species, it will be a tremendous task to understand the geographic distribution of each. Nevertheless, congruent with other recent studies [29–32] we have clearly shown the importance of using a group-specific PCR approach to better understand the cryptic diversity of protist groups. Studies of endemicity could be further undertaken by designing procedures that target each of the subgroups detected here and complemented with FISH and RNA sequencing strategies to verify that the species actually inhabit the location. For species or population demarcation, other faster evolving markers, such as ITS, may be needed. There is however a challenge when studying uncultured species to directly compare the phylogenetic clades based on different genes.

We provide here the first rigorous evidence for the existence of freshwater Telonemia. Two groups of freshwater sequences are identified showing that multiple and independent transitions from marine to freshwater have taken place during the evolutionary history of the group. It is obvious that the diversity of freshwater Telonemia is highly underestimated, and the ecological roles of Telonemia in these habitats are so far very much unclear. The possible stratification of species in freshwater is a first glimpse of potential differences in ecological adaptations - more studies combining molecular and microscopy approaches are clearly necessary to assess the diversity and dispersal patterns of Telonemia.

Methods

Environmental samples

Freshwater samples were collected from three different Norwegian lakes in May 2007; Lake Lutvann (59°54'N and 10°52'E) a small and deep (Zmax = 52 m) clearwater oligotrophic lake with long retention time, Lake Pollen (59°44'N and 10°45'E) a small and meromictic lake of intermediate depth (Zmax = 18 m) with only 7 m of freshwater and seawater in the monimolimnion, and Lake Sværsvann (59°48'N and 10°53'E) a small and shallow (Zmax = 11 m) meso- to polyhumic lake of complex morphology. Two litres of surface water (down to 50 cm) was collected from each lake and filtered through a Whatman GF/C glass-fiber filter with pore sizes of approximately 1 μm. Filters were dried and stored at -20°C. Sediment samples from Lake Lutvann were collected with a simple gravity corer at three depths, 50 m, 20 m and 5 m. The sediment samples from Lake Lutvann, including up to 500 ml of lake water were kept at 17°C with a 14/10 h light/dark cycle. 100 ml of culture of the cryptomonad species Plagioselmis nannoplanctica was added on average every three days for the Telonemia species to feed upon for seven days. P. nannoplanktica was grown in the freshwater media of Guillard & Lorenzen [56] without organic buffer.

Marine DNA was sampled from the following locations; Antarctica (59°22'S, 55°46W, December 1998), The Arctic Ocean (NOR26 and PD6 samples: 76°19'N, 23°45'E and NOR46 and AD6 samples: 76°20'N, 03°59'E, August 2002), The Mediterranean Sea (41°40'N, 2°48'E, January 2004) and the Indian Ocean (31°45'S, 52°37'E, May/June 1999). For sampling and DNA isolation methods see [11, 57–59].

DNA isolation and sequencing

DNA was isolated from the different freshwater samples by using the Power Max Soil DNA Isolation kit (MoBio, USA) following the manufacturers instructions. For DNA isolation from the sediments, 15 ml of sediment from the top layer were collected and centrifuged at 4000 rpm for 10 minutes. The isolated DNA was stored at -20°C.

Nested PCR was used to amplify the 18S rDNA gene from the freshwater samples with universal eukaryotic primers (based on PrimerA and PrimerB by Medlin et al. [60]; here called 1F and 1528R) used in the first round of amplification and Telonemia-specific primers used in the second round (for primer sequences, see Figure 2). The amplifications were done on an Eppendorf Mastercycler ep (Eppendorf, Germany) and a Biometra Thermocycler (Biometra, Germany) with a sample volume of 25 μl containing 10 - 200 ng of template DNA, 1 Ă— HotMaster Taq Buffer with 2.5 mM Mg2+ (5 Prime, USA), 200 μM dNTPs, 0.2 μM of each primer and 1.5 U HotMaster Taq DNA Polymerase (5 Prime, USA). The reaction mixture was incubated at 94°C for 2 min, followed by 30 - 34 cycles of 45 s at 94°C, 45 s at 60°C, 135 s at 72°C with a final extension at 72°C for 10 min. The PCR products were gel-extracted and purified using Wizard SV Gel and PCR Clean-Up System (Promega, USA), and cloned using TOPO TA Cloning Kit (Invitrogen, USA) following the manufacturers instructions. Colonies were checked for positive inserts by PCR amplification with the primers TopoF (5'-GGCTCGTATGTTGTGTGGAATTGT-3') and TopoR (5'-CCGTCGTTTTACAACGTCGTGACT-3') and identical reaction mixtures as described above, except that DynaZymeII (Finnzymes, Finland) DNA polymerase (1.5 U) and 1 Ă— DynaZyme buffer (F-511) were used. The PCR program was as follows: Initial denaturation at 95°C for 5 min, 34 cycles of 15 s at 95°C, 30 s at 60°C, 120 s at 72°C with a final extension at 72°C for 7 min. The positive inserts were sequenced on an ABI 3730 DNA Analyzer (Applied Biosystems, USA) with the primers M13F and M13R (Invitrogen, USA) using the ABI BigDye terminator v3.1 kit (Applied Biosystems, USA). 183 clones were randomly picked from the generated libraries and sequenced with the M13F primer (Invitrogen, USA). Identical, or nearly identical, sequences were not sequenced further. 82 of the inserts were full-length sequenced (approximately 1500 bp) with the M13R primer (Invitrogen, USA). Accession numbers for sequences generated in this study [GenBank: GQ365764-GQ365903 and GU117661-GU117693].

Figure 2
figure 2

Primers used in this study and their relative position in 18S rDNA gene. * indicates that primer is based on PrimerA and ** indicates that primer is based on PrimerB designed by Medlin et al. [55]. The 18S rDNA gene in the figure is based on the Telonema antarcticum sequence AJ564773 (1787 bp) in GenBank [62].

Phylogenetic analyses

Available sequences of possible Telonemia origin were identified by BLAST searches against the Entrez Nucleotide database [61, 62] using sequences of known Telonemia origin as query. The sequences identified from the BLAST searches were downloaded and pooled into a local database together with the sequences generated in this study. These sequences were added to an 18S rDNA alignment of all the major eukaryotic groups (hereafter called alignment 1) to confirm relationship to Telonemia. After removal of ambiguously aligned characters using the program MacClade version 4.07 [63], alignment 1 consisted of 374 taxa and 1465 characters. Alignment 1 was subjected to maximum likelihood (ML) analyses by using the program RAxML v.6 [64]. The topology with the highest likelihood score out of 100 heuristic searches, each from a random starting tree, was selected, and bootstrapping was done with 100 pseudoreplicates and one heuristic search per replicate. In the ML analyses, the General Time Reversible (GTR) model, with a gamma-distributed rate of variation across sites (G), was employed.

The ML analyses of alignment 1 showed that 198 sequences grouped together within Telonemia (results not shown). To be able to include more unambiguously aligned characters, a second alignment (alignment 2) was created with MacClade version 4.07 [63], consisting of the Telonemia sequences identified in the analysis of alignment 1. Identical sequences were excluded and the putative closest sister groups of Telonemia, the cryptomonads, haptophytes and katablepharids, were used as an outgroup [20]. Chimeric sequences were identified as described in [65]. The sequence NW614.39 is chimeric with the last 100 bp from a diatom. This part of the sequence was not included in the analyses. Accession numbers and clone names of sequences in alignment 2 are given in Additional file 1. Alignment 2 consisted of 159 taxa and 1758 characters. This alignment was analysed by ML (as for alignment 1) and Bayesian inferences. The Bayesian inferences were done with the program MrBayes [66] as follows: two independent runs, each with three cold and one heated MCMC (Markov Chain Monte Carlo) chains were started from a random starting tree. The two runs lasted for 4,000,000 generations. The covarion (COV) model was used together with the GTR+G+I to accommodate for different substitution rates across sites (G + proportion of invariable sites (I)) and across sequences (COV). The covarion model included two parameters, sites being on > off and off > on. All phylogenetic analyses were done on the freely available Bioportal at University of Oslo http://www.bioportal.uio.no.

References

  1. Lynch M: The Origins of Eukaryotic Gene Structure. Mol Biol Evol. 2006, 23 (2): 450-468. 10.1093/molbev/msj050.

    Article  CAS  PubMed  Google Scholar 

  2. Wilson AE, Sarnelle O, Neilan BA, Salmon TP, Gehringer MM, Hay ME: Genetic variation of the bloom-forming cyanobacterium Microcystis aeruginosa within and among lakes: Implications for harmful algal blooms. Appl Environ Microbiol. 2005, 71 (10): 6126-6133. 10.1128/AEM.71.10.6126-6133.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Snoke MS, Berendonk TU, Barth D, Lynch M: Large global effective population sizes in Paramecium. Mol Biol Evol. 2006, 23 (12): 2474-2479. 10.1093/molbev/msl128.

    Article  CAS  PubMed  Google Scholar 

  4. Foissner W: Biogeography and dispersal of micro-organisms: a review emphasizing protists. Acta Protozool. 2006, 45: 111-136.

    Google Scholar 

  5. Weisse T: Distribution and diversity of aquatic protists: an evolutionary and ecological perspective. Biodivers Conserv. 2008, 17 (2): 243-259. 10.1007/s10531-007-9249-4.

    Article  Google Scholar 

  6. Kristiansen J: 16. Dispersal of freshwater algae -- a review. Hydrobiologia. 1996, 336 (1): 151-157.

    Article  Google Scholar 

  7. Finlay BJ: Global Dispersal of Free-Living Microbial Eukaryote Species. Science. 2002, 296 (5570): 1061-1063. 10.1126/science.1070710.

    Article  CAS  PubMed  Google Scholar 

  8. Fenchel T, Finlay BJ: The Ubiquity of Small Species: Patterns of Local and Global Diversity. Bioscience. 2004, 54: 777-784. 10.1641/0006-3568(2004)054[0777:TUOSSP]2.0.CO;2.

    Article  Google Scholar 

  9. Baas-Becking LGM: Geobiologie of Inleiding Tot de Milieukunde. 1934, The Hague: Van Stockkum & Zoon

    Google Scholar 

  10. de Wit R, Bouvier T: 'Everything is everywhere, but, the environment selects'; what did Baas Becking and Beijerinck really say?. Environ Microbiol. 2006, 8 (4): 755-758. 10.1111/j.1462-2920.2006.01017.x.

    Article  PubMed  Google Scholar 

  11. Massana R, Balague V, Guillou L, Pedros-Alio C: Picoeukaryotic diversity in an oligotrophic coastal site studied by molecular and culturing approaches. FEMS Microbiol Ecol. 2004, 50 (3): 231-243. 10.1016/j.femsec.2004.07.001.

    Article  CAS  PubMed  Google Scholar 

  12. Casamatta DA, Vis ML, Sheath RG: Cryptic species in cyanobacterial systematics: a case study of Phormidium retzii (Oscillatoriales) using RAPD molecular markers and 16S rDNA sequence data. Aquat Bot. 2003, 77 (4): 295-309. 10.1016/j.aquabot.2003.08.005.

    Article  CAS  Google Scholar 

  13. Pawlowski J, Holzmann M: Molecular phylogeny of Foraminifera a review. Eur J Protistol. 2002, 38 (1): 1-10. 10.1078/0932-4739-00857.

    Article  Google Scholar 

  14. Moon-van der Staay SY, De Wachter R, Vaulot D: Oceanic 18S rDNA sequences from picoplankton reveal unsuspected eukaryotic diversity. Nature. 2001, 409 (6820): 607-610. 10.1038/35054541.

    Article  CAS  PubMed  Google Scholar 

  15. Not F, Valentin K, Romari K, Lovejoy C, Massana R, Tobe K, Vaulot D, Medlin LK: Picobiliphytes: A Marine Picoplanktonic Algal Group with Unknown Affinities to Other Eukaryotes. Science. 2007, 315 (5809): 253-255. 10.1126/science.1136264.

    Article  CAS  PubMed  Google Scholar 

  16. Dawson SC, Pace NR: Novel kingdom-level eukaryotic diversity in anoxic environments. Proc Natl Acad Sci USA. 2002, 99 (12): 8324-8329. 10.1073/pnas.062169599.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Habura A, Pawlowski JAN, Hanes SD, Bowser SS: Unexpected Foraminiferal Diversity Revealed by Small-subunit rDNA Analysis of Antarctic Sediment. J Eukaryot Microbiol. 2004, 51 (2): 173-179. 10.1111/j.1550-7408.2004.tb00542.x.

    Article  CAS  PubMed  Google Scholar 

  18. Holzmann M, Habura A, Giles H, Bowser SS, Pawlowski JAN: Freshwater Foraminiferans Revealed by Analysis of Environmental DNA Samples. J Eukaryot Microbiol. 2003, 50 (2): 135-139. 10.1111/j.1550-7408.2003.tb00248.x.

    Article  PubMed  Google Scholar 

  19. LĂ³pez-GarcĂ­a P, RodrĂ­guez-Valera F, PedrĂ³s-AliĂ³ C, Moreira D: Unexpected diversity of small eukaryotes in deep-sea Antarctic plankton. Nature. 2001, 409 (6820): 603-607. 10.1038/35054537.

    Article  PubMed  Google Scholar 

  20. Shalchian-Tabrizi K, Eikrem W, Klaveness D, Vaulot D, Minge MA, Le Gall F, Romari K, Throndsen J, Botnen A, Massana R: Telonemia, a new protist phylum with affinity to chromist lineages. Proc Biol Sci. 2006, 273 (1595): 1833-1842. 10.1098/rspb.2006.3515.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Boenigk J, Pfandl K, Stadler P, Chatzinotas A: High diversity of the 'Spumella-like' flagellates: an investigation based on the SSU rRNA gene sequences of isolates from habitats located in six different geographic regions. Environ Microbiol. 2005, 7 (5): 685-697. 10.1111/j.1462-2920.2005.00743.x.

    Article  CAS  PubMed  Google Scholar 

  22. Koch TA, Ekelund F: Strains of the heterotrophic flagellate Bodo designis from different environments vary considerably with respect to salinity preference and SSU rRNA gene composition. Protist. 2005, 156 (1): 97-112. 10.1016/j.protis.2004.12.001.

    Article  CAS  PubMed  Google Scholar 

  23. Foissner W: Protist diversity: estimates of the near-imponderable. Protist. 1999, 150 (4): 363-368. 10.1016/S1434-4610(99)70037-4.

    Article  CAS  PubMed  Google Scholar 

  24. Foissner W: Protist diversity and distribution: some basic considerations. Biodivers Conserv. 2008, 17 (2): 235-242. 10.1007/s10531-007-9248-5.

    Article  Google Scholar 

  25. Bass D, Richards TA, Matthai L, Marsh V, Cavalier-Smith T: DNA evidence for global dispersal and probable endemicity of protozoa. BMC Evol Biol. 2007, 7: 162-10.1186/1471-2148-7-162.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Jeon S, Bunge J, Leslin C, Stoeck T, Hong S, Epstein S: Environmental rRNA inventories miss over half of protistan diversity. BMC Microbiol. 2008, 8 (1): 222-10.1186/1471-2180-8-222.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Stoeck T, Hayward B, Taylor GT, Varela R, Epstein SS: A Multiple PCR-primer Approach to Access the Microeukaryotic Diversity in Environmental Samples. Protist. 2006, 157 (1): 31-43. 10.1016/j.protis.2005.10.004.

    Article  CAS  PubMed  Google Scholar 

  28. Potvin M, Lovejoy C: PCR-Based Diversity Estimates of Artificial and Environmental 18S rRNA Gene Libraries. J Eukaryot Microbiol. 2009, 56 (2): 174-181. 10.1111/j.1550-7408.2008.00386.x.

    Article  CAS  PubMed  Google Scholar 

  29. Lin S, Zhang H, Hou Y, Miranda L, Bhattacharya D: Development of a Dinoflagellate-Oriented PCR Primer Set Leads to Detection of Picoplanktonic Dinoflagellates from Long Island Sound. Appl Environ Microbiol. 2006, 72 (8): 5626-5630. 10.1128/AEM.00586-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Bass D, Cavalier-Smith T: Phylum-specific environmental DNA analysis reveals remarkably high global biodiversity of Cercozoa (Protozoa). Int J Syst Evol Microbiol. 2004, 54 (6): 2393-2404. 10.1099/ijs.0.63229-0.

    Article  CAS  PubMed  Google Scholar 

  31. Viprey M, Guillou L, Ferréol M, Vaulot D: Wide genetic diversity of picoplanktonic green algae (Chloroplastida) in the Mediterranean Sea uncovered by a phylum-biased PCR approach. Environ Microbiol. 2008, 10 (7): 1804-1822. 10.1111/j.1462-2920.2008.01602.x.

    Article  CAS  PubMed  Google Scholar 

  32. Lara E, Moreira D, Vereshchaka A, LĂ³pez-GarcĂ­a P: Pan-oceanic distribution of new highly diverse clades of deep-sea diplonemids. Environ Microbiol. 2009, 11 (1): 47-55. 10.1111/j.1462-2920.2008.01737.x.

    Article  CAS  PubMed  Google Scholar 

  33. Zuendorf A, Bunge J, Behnke A, Barger KJA, Stoeck T: Diversity estimates of microeukaryotes below the chemocline of the anoxic Mariager Fjord, Denmark. FEMS Microbiol Ecol. 2006, 58 (3): 476-491. 10.1111/j.1574-6941.2006.00171.x.

    Article  CAS  PubMed  Google Scholar 

  34. Lovejoy C, Massana R, Pedros-Alio C: Diversity and Distribution of Marine Microbial Eukaryotes in the Arctic Ocean and Adjacent Seas. Appl Environ Microbiol. 2006, 72 (5): 3085-3095. 10.1128/AEM.72.5.3085-3095.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Not F, Latasa M, Scharek R, Viprey M, Karleskind P, BalaguĂˆ V, Ontoria-Oviedo I, Cumino A, Goetze E, Vaulot D: Protistan assemblages across the Indian Ocean, with a specific emphasis on the picoeukaryotes. Deep Sea Research Part I: Oceanographic Research Papers. 2008, 55 (11): 1456-1473. 10.1016/j.dsr.2008.06.007.

    Article  Google Scholar 

  36. Shalchian-Tabrizi K, Kauserud H, Massana R, Klaveness D, Jakobsen KS: Analysis of environmental 18S ribosomal RNA sequences reveals unknown diversity of the cosmopolitan phylum Telonemia. Protist. 2007, 158 (2): 173-180. 10.1016/j.protis.2006.10.003.

    Article  PubMed  Google Scholar 

  37. Klaveness D, Shalchian-Tabrizi K, Thomsen HA, Eikrem W, Jakobsen KS: Telonema antarcticum sp. nov., a common marine phagotrophic flagellate. Int J Syst Evol Microbiol. 2005, 55 (Pt 6): 2595-2604. 10.1099/ijs.0.63652-0.

    Article  CAS  PubMed  Google Scholar 

  38. Countway PD, Gast RJ, Dennett MR, Savai P, Rose JM, Caron DA: Distinct protistan assemblages characterize the euphotic zone and deep sea (2500 m) of the western North Atlantic (Sargasso Sea and Gulf Stream). Environ Microbiol. 2007, 9 (17472636): 1219-1232. 10.1111/j.1462-2920.2007.01243.x.

    Article  CAS  PubMed  Google Scholar 

  39. Vørs N: Heterotrofe protister (ekskl. dinoflagellater, loricabærende choanoflagellater og ciliater). 1992, Copenhagen: Havforskning fra Miljøstyrelsen

    Google Scholar 

  40. Tong S, Vørs N, Patterson DJ: Heterotrophic flagellates, centrohelid heliozoa and filose amoebae from marine and freshwater sites in the Antarctic. Polar Biol. 1997, 18 (2): 91-106. 10.1007/s003000050163.

    Article  Google Scholar 

  41. Laybourn-Parry J, Ellis-Evans JC, Butler H: Microbial dynamics during the summer ice-loss phase in maritime Antarctic lakes. J Plankton Res. 1996, 18 (4): 495-511. 10.1093/plankt/18.4.495.

    Article  Google Scholar 

  42. Throndsen J: Flagellates of Norwegian coastal waters. Nytt Magasin Botanikk. 1969, 161-216. 16

  43. Lefèvre E, Roussel B, Amblard C, Sime-Ngando T: The Molecular Diversity of Freshwater Picoeukaryotes Reveals High Occurrence of Putative Parasitoids in the Plankton. PLoS ONE. 2008, 3 (6): e2324-10.1371/journal.pone.0002324.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Logares R, Rengefors K, Kremp A, Shalchian-Tabrizi K, Boltovskoy A, Tengs T, Shurtleff A, Klaveness D: Phenotypically Different Microalgal Morphospecies with Identical Ribosomal DNA: A Case of Rapid Adaptive Evolution?. Microb Ecol. 2007, 53 (4): 549-561. 10.1007/s00248-006-9088-y.

    Article  CAS  PubMed  Google Scholar 

  45. Montresor M, Lovejoy C, Orsini L, Procaccini G, Roy S: Bipolar distribution of the cyst-forming dinoflagellate Polarella glacialis. Polar Biol. 2003, 26 (3): 186-194.

    Google Scholar 

  46. Darling KF, Wade CM, Stewart IA, Kroon D, Dingle R, Brown AJL: Molecular evidence for genetic mixing of Arctic and Antarctic subpolar populations of planktonic foraminifers. Nature. 2000, 405 (6782): 43-47. 10.1038/35011002.

    Article  CAS  PubMed  Google Scholar 

  47. ZĂºĂ±iga L, Campos V, Pinochet H, Prado B: A limnological reconnaissance of Lake Tebenquiche, Salar de Atacama, Chile. Hydrobiologia. 1991, 210 (1): 19-24.

    Article  Google Scholar 

  48. Krawczyk WE, Lefauconnier B, Pettersson LE: Chemical denudation rates in the Bayelva Catchment, Svalbard, in the Fall of 2000. Physics and Chemistry of the Earth, Parts A/B/C. 2003, 28 (28-32): 1257-1271. 10.1016/j.pce.2003.08.054.

    Article  Google Scholar 

  49. Ikävalko J, Thomsen HA: The Baltic Sea ice biota (March 1994): study of the protistan community. Eur J Protistol. 1998, 33: 229-243.

    Article  Google Scholar 

  50. Logares R, BrĂ¥te J, Bertilsson S, Clasen JL, Shalchian-Tabrizi K, Rengefors K: Infrequent marine-freshwater transitions in the microbial world. Trends Microbiol. 2009, 17 (9): 414-422. 10.1016/j.tim.2009.05.010.

    Article  CAS  PubMed  Google Scholar 

  51. von der Heyden S, Chao EE, Vickerman K, Cavalier-Smith T: Ribosomal RNA phylogeny of bodonid and diplonemid flagellates and the evolution of euglenozoa. J Eukaryot Microbiol. 2004, 51 (4): 402-416. 10.1111/j.1550-7408.2004.tb00387.x.

    Article  CAS  PubMed  Google Scholar 

  52. von der Heyden S, Chao EE, Cavalier-Smith T: Genetic diversity of goniomonads: an ancient divergence between marine and freshwater species. Eur J Phycol. 2004, 39 (4): 343-350. 10.1080/09670260400005567.

    Article  CAS  Google Scholar 

  53. Shalchian-Tabrizi K, BrĂ¥te J, Logares R, Klaveness D, Berney C, Jakobsen KS: Diversification of unicellular eukaryotes: Cryptomonad colonisations of marine and fresh waters inferred from revised 18S rRNA phylogeny. Environ Microbiol. 2008, 10 (10): 2635-2644. 10.1111/j.1462-2920.2008.01685.x.

    Article  CAS  PubMed  Google Scholar 

  54. Logares R, Shalchian-Tabrizi K, Boltovskoy A, Rengefors K: Extensive dinoflagellate phylogenies indicate infrequent marine-freshwater transitions. Mol Phylogenet Evol. 2007, 45 (3): 887-903.

    Article  CAS  PubMed  Google Scholar 

  55. BrĂ¥te J, Logares R, Berney C, Ree DK, Klaveness D, Jakobsen KS, Shalchian-Tabrizi K: Freshwater Perkinsea and marine-freshwater colonizations revealed by pyrosequencing and phylogeny of environmental DNA. ISME Journal. 2010

    Google Scholar 

  56. Guillard RRL, Lorenzen CJ: Yellow-green algae with chlorophyllide c. J Phycol. 1972, 8 (1): 10-14.

    CAS  Google Scholar 

  57. Diez B, Pedros-Alio C, Massana R: Study of genetic diversity of eukaryotic picoplankton in different oceanic regions by small-subunit rRNA gene cloning and sequencing. Appl Environ Microbiol. 2001, 67 (7): 2932-2941. 10.1128/AEM.67.7.2932-2941.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Not F, Massana R, Latasa M, Marie D, Colson C, Eikrem W, Pedros-Alio C, Vaulot D, Simon N: Late summer community composition and abundance of photosynthetic picoeukaryotes in Norwegian and Barents Seas. Limnol Oceanogr. 2005, 50 (5): 1677-1686. 10.4319/lo.2005.50.5.1677.

    Article  CAS  Google Scholar 

  59. Massana R, Guillou L, Terrado R, Forn I, Pedros-Alio C: Growth of uncultured heterotrophic flagellates in unamended seawater incubations. Aquat Microb Ecol. 2006, 45 (2): 171-180. 10.3354/ame045171.

    Article  Google Scholar 

  60. Medlin L, Elwood HJ, Stickel S, Sogin ML: The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene. 1988, 71 (2): 491-499. 10.1016/0378-1119(88)90066-2.

    Article  CAS  PubMed  Google Scholar 

  61. Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  62. Entrez Nucleotide database.http://www.ncbi.nlm.nih.gov/sites/entrez?db=nuccore

  63. Maddison D, Maddison W: MacClade 4: Analysis of Phylogeny and Character Evolution. 2000, Sinauer Associates, Sunderland, MA, 4

    Google Scholar 

  64. Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.

    Article  CAS  PubMed  Google Scholar 

  65. Berney C, Fahrni J, Pawlowski J: How many novel eukaryotic 'kingdoms'? Pitfalls and limitations of environmental DNA surveys. BMC Biol. 2004, 2: 13-10.1186/1741-7007-2-13.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Ramon Massana for marine DNA samples, Liisa Lepistö for providing unpublished data and Cédric Berney for identification of chimeric sequences. We would also like to thank the Bioportal http://www.bioportal.uio.no for computer resources. This work was supported by grants from the Norwegian Research Council to KSJ and UiO grants to KST and JB.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kamran Shalchian-Tabrizi.

Additional information

Authors' contributions

JB collected the freshwater samples, generated the sequence data, performed the phylogenetic analyses, and wrote the manuscript. TR collected environmental samples and designed and tested the Telonemia specific primers. DK supervised and participated in the sample collection and manuscript writing. KSJ funded and coordinated the study and contributed to writing the manuscript. All authors have read and approved the final manuscript. KST designed the project, supervised the analyses and interpretation of the molecular phylogenies and participated in writing the manuscript.

Electronic supplementary material

12866_2009_1124_MOESM1_ESM.DOCX

Additional file 1: Supplementary table Description of sequences used in the phylogenetic analyses in Figure 1. Sequences in bold are generated in this study. (DOCX 105 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

BrĂ¥te, J., Klaveness, D., Rygh, T. et al. Telonemia-specific environmental 18S rDNA PCR reveals unknown diversity and multiple marine-freshwater colonizations. BMC Microbiol 10, 168 (2010). https://doi.org/10.1186/1471-2180-10-168

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2180-10-168

Keywords