Evidence for isolated evolution of deep-sea ciliate communities through geological separation and environmental selection
© Stock et al.; licensee BioMed Central Ltd. 2013
Received: 27 November 2012
Accepted: 15 May 2013
Published: 8 July 2013
Deep hypersaline anoxic basins (DHABs) are isolated habitats at the bottom of the eastern Mediterranean Sea, which originate from the ancient dissolution of Messinian evaporites. The different basins have recruited their original biota from the same source, but their geological evolution eventually constituted sharp environmental barriers, restricting genetic exchange between the individual basins. Therefore, DHABs are unique model systems to assess the effect of geological events and environmental conditions on the evolution and diversification of protistan plankton. Here, we examine evidence for isolated evolution of unicellular eukaryote protistan plankton communities driven by geological separation and environmental selection. We specifically focused on ciliated protists as a major component of protistan DHAB plankton by pyrosequencing the hypervariable V4 fragment of the small subunit ribosomal RNA. Geospatial distributions and responses of marine ciliates to differential hydrochemistries suggest strong physical and chemical barriers to dispersal that influence the evolution of this plankton group.
Ciliate communities in the brines of four investigated DHABs are distinctively different from ciliate communities in the interfaces (haloclines) immediately above the brines. While the interface ciliate communities from different sites are relatively similar to each other, the brine ciliate communities are significantly different between sites. We found no distance-decay relationship, and canonical correspondence analyses identified oxygen and sodium as most important hydrochemical parameters explaining the partitioning of diversity between interface and brine ciliate communities. However, none of the analyzed hydrochemical parameters explained the significant differences between brine ciliate communities in different basins.
Our data indicate a frequent genetic exchange in the deep-sea water above the brines. The “isolated island character” of the different brines, that resulted from geological events and contemporary environmental conditions, create selective pressures driving evolutionary processes, and with time, lead to speciation and shape protistan community composition. We conclude that community assembly in DHABs is a mixture of isolated evolution (as evidenced by small changes in V4 primary structure in some taxa) and species sorting (as indicated by the regional absence/presence of individual taxon groups on high levels in taxonomic hierarchy).
KeywordsCiliates Hypersaline Deep-sea anoxic basins DHABs Brine Species sorting Environmental filtering Niche separation
DHABs in the Eastern Mediterranean Sea, located at depths of more than 3000 m below sea level, result from the dissolution of outcropping ancient subterranean salt deposits from the Messinian salinity crisis (late Miocene period, > 5 million years ago, ). Due to the high densities of the brines (up to 1.23 kg m-3, ), mixing of these water masses with overlying deep-sea water (average density: 1.03 kg m-3) is restricted, resulting in anoxic conditions in these brines. An interface (halocline: depending on the basin, typically 1 to 3 m thick) separates the anoxic brine from the normoxic and normsaline deep-sea water. Due to the dissolution of different strata of the evaporites from the Messinian salinity crisis, the hydrochemistries of the Eastern Mediterranean Sea DHABs differ significantly. For example, while salinity in some basins (Thetis, L’Atalante, Bannock and Tyro) ranges between 321 and 352 g l-1 (nearly 10 times higher than average seawater salinity), others exhibit a much lower salinity (such as Urania brine 240 g l-1). Potassium ions range between 19 and 300 mmol l-1, magnesium ions between 71 and 792 mmol l-1 sulfate between 52 and 323 mmol l-1, sulfide between 2.1 and 15 mmol l-1 and methane between 0.4 and 5.6 mmol l-1. Because of their unique hydrochemistries and physical separation for thousands of years, the DHABs may serve as island habitats and provide an ideal scenario to test the hypothesis that species sorting of planktonic ciliate communities results from environmental filtering through niche separation.
Molecular diversity surveys of protists, employing domain-specific PCR primers for the amplification of taxonomic marker genes (small subunit ribosomal RNA, SSU rRNA), clone library construction and Sanger sequencing revealed, that ciliates are among the most diverse and abundant plankton taxa thriving in some of the Eastern Mediterranean DHABs [2, 3]. Ciliates, through their grazing activities on bacteria, archaea and smaller eukaryotes are central players in the marine microbial loop [7–9] and species composition of ciliates can serve as an indicator of environmental health . They have been used extensively as model organisms to develop and test ideas about microbial biodiversity and biogeography (e.g. [11–17]). One major reason for this is that compared to amoeboid and flagellated organisms, they are morphologically diverse [18, 19] and there is a long history of their taxonomic and phylogenetic study (reviewed in ). The extensive foundation of knowledge on ciliate species and their inferred relationships facilitates data evaluation and hypothesis testing for studies that aim to explore ciliate biodiversity, evolution and biogeography.
None of the previous taxon samplings of SSU rRNA signatures in initial DHAB protistan diversity surveys reached saturation [2, 3], as is generally the case in cloning and Sanger sequencing-based strategies [20–24]. Therefore, it was not possible to observe many patterns of diversity in previous studies of DHABs, nor was it possible to perform comparative statistical analyses of specific taxa in the DHABs. In this study we applied a high-throughput next generation sequencing strategy (pyrosequencing) and a ciliate-specific primer set in order to recover a comprehensive dataset on this target group. The resulting data from deep sequencing enabled us to address basic ecological questions. Our first hypothesis was that the distinct chemistries of the different basins would drive species sorting in planktonic ciliate communities in the brines and interfaces of each basin. If this hypothesis is true, we would expect (i) that interface communities will differ decisively from brine communities (environmental filtering) and (ii) that ciliate communities in interfaces are more similar to each other than in the brines (isolated island character of brine basins). The brines of the different basins are isolated from one another due to the sharp density gradient that exists between these hypersaline basins and overlying Mediterranean seawater. In contrast, exchange may be possible between interface populations in different DHABs since some exchange is possible between seawater and the typically ca. 2 m-thick interfaces (haloclines). Our second hypothesis was that ciliate community composition in the brines and interfaces of these four DHABs, separated by up to 500 km, would not be significantly affected by distance between basins. If this hypothesis is true, we would expect no significant correlation between pairs of samples and geographic distance between the respective sampling sites, therefore, no isolation with distance.
Summary of ciliate V4 SSU rRNA amplicon data for each sample including the ciliate cluster numbers at a level of 95% cutoff
Number of V4-amplicons
Number of ciliate clusters
Before quality control
After quality control
Cutoff level 95%
Ciliate diversity in the DHABs
Alpha diversity indices (data normalized to 32,663 sequences in each sample) of ciliate communities in DHAB interfaces and brines
1.285 ± 0.002
58.0 ± 3.3
1.477 ± 0.004
44.6 ± 3.5
1.139 ± 0.004
42.4 ± 3.4
1.067 ± 0.003
42.9 ± 2.0
0.142 ± 0.001
16.9 ± 1.2
0.895 ± 0.004
33.9 ± 6.5
1.730 ± 0.004
47.5 ± 3.0
Putative taxonomy of ciliate amplicons
The V4-amplicons analyzed in this study were related to a total of 102 identified ciliate genera and one unclassified ciliate taxon (Additional file 3: Table S1). The unique character of the Medee brine ciliate community can be inferred from Figure 2b, which displays the taxonomy assigned to the ciliate amplicons obtained from each sampling site. Medee brine was dominated by amplicons (n = 33,961; 97% of all amplicons), which were all related to the genus Anoplophrya (Astomatida) as closest BLAST match in NCBIs GenBank nr database. The sequence similarities of these amplicons to Anoplophrya ranged between 80 and 89% (Additional file 4: Table S2). The remaining 1021 ciliate amplicons from Medee brine were related to a few other taxon groups belonging predominantly to the Peniculida (2.7%), other Astomatida (0.06%), and Pleuronematida (0.03%). Thetis brine and Tyro brine had a relatively similar ciliate community composition, both of which were dominated by amplicons that have Strombidium as the closest BLAST match in the GenBank nucleotide database (64% and 45%, of all amplicons, respectively). Other abundant taxon groups shared by these two samples were Novistrombidium (30% in Tyro brine and 9% in Thetis brine), and Pseudotontonia (4% in Tyro brine and 8% in Thetis brine). While Laboea accounted for 11% of all amplicons in Thetis brine, this taxon group was absent in Tyro brine. A tintinnid ciliate taxon related to Salpingella as closest database relative occured exclusively in Tyro brine (4% of all amplicons), but not in Thetis (Additional file 3: Table S1). The ciliate community composition in Urania brine was dissimilar to the brines in Tyro and Thetis basins. One striking quantitative difference was the high proportion of Pseudotontonia-related amplicons (40%) in Urania brine. However, while most of the relatively abundant taxon-groups were shared between these three brine samples (but in different quantities), most qualitative differences between Tyro, Thetis and Urania brines were attributed to taxon groups with lower abundances. Medee brine was distinct in its ciliate composition from other brines.
Tyro interface stood out from the other interface samples. The most significant difference was the occurrence of 14,337 amplicons (41%), with Apocoleps (Prorodontida) as the best BLAST match. The proportion of amplicons in Thetis, Urania and Medee interfaces related to this taxon was less than 0.5%. Also the proportion of Strombidium-like amplicons in Tyro interface (40%) was decisively higher compared to the other interfaces (4-21%). Thetis interface and Urania interface had a very similar taxon composition, dominated by amplicons most closely related to Pleuronema (Pleuronematida) (70% in UIF and 57% in ThIF). This taxon was also highly represented in Medee interface (49%). The second most abundant taxon group in Medee interface were clevelandellids, represented with 43%. This taxon was underrepresented in the interfaces of other basins (0.02% in UIF – 4% in ThIF).
Four taxa occured in all eight samples analyzed (closest BLAST matches: Pleuronema, Strombidium, Omegastrombidium, Apocoleps). Four taxa were exclusive to all interfaces (Palgiopyliella, Cyclidium, Schizocalpytra, Isochonida). Interestingly, not a single taxon occured exclusively in all brines simultaneously. However, 28 taxon groups were absent from interfaces but present in at least one of the brines. The same number of taxon groups was absent from all brines but occured in at least one of the interfaces. The majority of taxon groups had abundances accounting for less than 5% of all amplicons obtained within a sample.
Relating community patterns to environmental variables
The ciliate communities in the DHAB interfaces showed only small variation along the first axis, while brine samples spread across a wider range of this first axis, with Medee brine and Thetis brine defining the longest distance. None of the CCAs conducted found a meaningful correlation of this axis with any environmental variable that we have measured and tested explaining this first axis. However, it must be a factor that only separates niches for the brine communities, but not for interface communities.
Distance effect on DHAB ciliate community profiles
Fluorescent in situ hybridization and scanning electron microscopy
Deep hypersaline anoxic basins (DHABs) in the Eastern Mediterranean Sea are ideally suited for testing the effect of historical contingencies on the evolution of protist communities. The distance between individual basins is variable, and each basin is characterized by hydrochemical gradients (interfaces to brines), and slightly different origins, leading to differences in physicochemical factors of the brines and interfaces in each of the different basins. Due to the steep density gradients along the interfaces of these basins, there is little connectivity between basin brines and overlying seawater, and therefore, between basin brines. First insights into the ciliate communities in the mesopelagic realm above the brine basins came from a Sanger sequencing-based approach . Because of the relatively small amount of data (four ciliate OTUs in the mesopelagic reference and 10 in the brine) it is not a reliable dataset for comparison to the high throughput sequencing data from this study. However, the data from that preliminary study did indicate a significant community shift between the water column and the basin brines. We assessed ciliate community structures in the interfaces and brines of several basins in order to determine the degree to which these environmental barriers and basin chemistries influenced the ciliate plankton.
The proportion of rare versus abundant ciliate taxa in this study (Additional file 2: Figure S2) is comparable to previous findings reporting protistan communities with few abundant and many rare taxa [16, 24]. Several molecular diversity surveys over different spatial scales ranging from centimeters to tens of thousands of kilometers have supported distance-decay relationships (effect of distance on spatial interactions) for microbial organisms, including bacteria (e.g. [26, 27]), archaea (e.g. ), fungi (e.g. ) and also protists (e.g. [30–32]). Even organisms with large population sizes and the potential to spread globally using spores, which were assumed to be cosmopolitan [13, 33], show significant non-random spatial distribution patterns . However, in our study of ciliate communities in these DHABs, a similar distance-decay relationship was not observed (insignificant correlation between Bray-Curtis and geographic distances in Pearson correlation and Mantel test). A potential explanation could be that the small number of compared locations may have masked true patterns. Alternatively, the presence of a metacommunity  within the Mediterranean Sea could cause the absence of a significant heterogeneous distribution [36, 37]. In limnic systems geographic distance has been found to influence asymmetric latitudinal genus richness patterns between 42° S and the pole . However, this seems to be a fundamental difference between marine and “terrestrial” (land-locked) systems. Furthermore, on a global scale, historical factors were significantly more responsible for the geographic patterns in community composition of diatoms than environmental conditions . In other marine studies ciliates showed variations in taxonomic composition between closely related samples, which were explained by environmental factors rather than distance . Similarly, in our study geographic distance could not explain the variations observed between the ciliate communities. Instead, hydrochemistry explained some of the variation in observed ciliate community patterns, and there was a strong separation of halocline interface and brine communities (Figure 3). The DHAB interfaces are characterized by extremely steep physicochemical gradients on a small spatial scale typically less than a couple of meters (for example, only 70 cm in Medee, ). The concentrations of salt and oxygen are the most prominent environmental factors that change dramatically along the interfaces into the brines. In a recent metadata-analysis of environmental sequence data, these two factors were identified as strong selection factors for ciliates . Also for bacterial communities, salt concentration emerged as the strongest factor influencing global distribution . Likewise, the bacterioplankton community composition in coastal Antarctic lakes was weakly related with geographical distance, but strongly correlated with salinity . Accordingly, Logares et al. assumed that a long-term salinity change ‘promoted the diversification of bacterioplankton communities by modifying the composition of ancestral communities and by allowing the establishment of new taxa’ . Salinity shifts characterize a boundary which is one of the most difficult barriers to cross for organisms in all three domains of life . While mechanisms to cope with high salt concentrations are relatively well studied in prokaryotes, they are still largely unknown in protists (with the exception of the model algae Dunaliella salina). While there is evidence that many protists have narrow ranges of salt tolerance [45, 46], some taxa are known to occur under a wide range of salinities, from freshwater to hypersaline . One example is the ciliate Cyclidium glaucoma, which may explain the occurrence of some of the same phylotypes in haloclines and brines of specific DHABs. Other examples are likely to exist.
In contrast, adaptations to anoxia in ciliates are well known. Ciliates are one of the most successful eukaryotic taxon groups in hypoxic and anoxic habitats. In their long evolutionary history, they have acquired several strategies that allow for an anaerobic lifestyle, including hydrogenosomes [49, 50], anaerobic mitochondria , and/or symbiotic networks [52, 53]. The high taxonomic diversity of anaerobe ciliates includes taxa such as Nyctotherus, Loxodes, Pleuronema, Strombidium, Trimyema, Cyclidium and Metopus, some of which were also detected in our genetic diversity survey. Electron microscopy and fluorescence in situ hybridization assays provide unbiased evidence that the genetic signatures we detected in our rRNA-targeted gene survey can be assigned to ciliates living in the DHABs rather than reflecting ancient nucleic acids. (Figure 5, [25, 54]). Taking advantage of phylotypes that we detected exclusively in specific habitats and phylotypes that can be found in several habitats with distinct hydrochemical characteristics, we may assume that the latter have a character of more generalist taxa compared to the more locally restricted phylotypes. The total number of observed taxon groups is 102 distributed over eight different datasets (samples or habitats) (Additional file 1: Figure S1). In those eight samples there are 13 generalist taxonomic groups that appeared simultaneously in at least six of the datasets. Only four taxonomic groups appeared in all of the eight datasets. Specialists, i.e. taxa that are restricted to a single unique habitat account for 34 different taxonomic groups. This results in a specialist/generalist ratio of 8.5 to 1, indicating a high specialization of taxa in the habitats under study. However, there is a limitation to infer the autecology of specific evolutionary lineages based on sequence data and microscopy evidence . We do not make any attempt to explain the presence or absence of specific phylotypes in individual samples, and we instead focus only on community level ciliate diversity.
Hydrochemical gradients have been identified as environmental barriers in previous protistan diversity surveys including the Cariaco basin , Framvaren Fjord , Mariager Fjord , Baltic Sea  and the Black Sea  and also stratified lakes . Therefore, the clear distinction of halocline ciliate communities from brine communities is not an unexpected result. However, it is surprising that the environmental variables we measured had a minor contribution to differences among the individual brine ciliate communities. In the CCA analyses (Figure 3) the different brine communities were spread out along the y-axis. This axis, however, does not represent an environmental gradient. This is surprising, considering that different types of salts may have different physiological effects  and therefore, should require different adaptation strategies in halophiles. Basically, we can assume two scenarios: first, for isolated evolution as described in , the scenario starts with a seed taxon. After physical separation of the original habitat into two habitats neutral mutations are changing the seed taxon in these habitats independently. These neutral mutations are of minor nature considering the time scale of the basins’ geological histories. From this event we would expect similar taxon groups with only minor genetic changes in both habitats. As mentioned above, each eighth taxon recorded in our study (Additional file 3: Table S1) falls into this category. In the second scenario (environmental filtering) we have the same ‘seed bank’ community for different basins. Through environmental filtering (different hydrochemistries of the basins) some taxa may go extinct, others have the genomic potential to adapt to some specific hydrochemistries, while others are genomically equipped for adaptation to other environmental conditions. In this case we would find taxa differing on higher taxonomic (genetic) hierarchies. This is the case for 34 of 102 detected taxon groups (Additional file 3: Table S1). We cannot rule out all environmental factors from causing differences between the ciliate communities because we did not measure all possible environmental factors, but only the hydrogeochemical factors that account for the most pronounced and obvious differences. This suggests that (1) other hydrochemical variables we did not measure are leading to this separation, or (2) that biotic interactions may explain some of the differences between brine ciliate communities. Even though interactions of top-down and bottom-up factors in shaping community structures of aquatic microbes are still poorly understood  some well known biotic interactions could be considered. Such biotic interactions may be, for example, parasitic relationships between organisms like amoeboid parasitic forms that can shape the composition of cyanobacterial species in lakes (Rohrlack et al., unpublished data). Furthermore, symbiotic relationships include ciliates, associated with epibiontic prokaryotes, were found to be the dominant eukaryotic morphotypes in the haloclines of DHABs in the Mediterranean Sea . Biotic interaction between protists and viruses are also known and have been shown . Viruses specifically infect protists, e.g. the Coccolithovirus and it’s host, the calicifying haptophyte Emiliania huxleyi. Additionally, viruses can also have an an indirect influence on protists by infecting the bacteria on which the protistan grazers feed or protistan grazers can even feed directly on viruses even though the carbon transfer to the higher trophic level is of minor importance . Furthermore, different bacterioplankton communities can produce a bottom-up control on grazing protists. Namely, the growth efficiency of protists can relate strongly to the available bacterial prey [63, 67]. This is highly likely because differences in bacterial community composition in DHABs have been shown before [68, 69]. That leads to the assumption that different bacterial communities support different phagotrophic protists that show strong preferences for particular prey species [63, 67, 70, 71] or morphotypes [72, 73]. Other possible explanations are founder effects, which describe a genetic deviation of an isolated population or founder population (on an island for example) compared to the original population based on a low number of alleles within the founders individuals , random effects or genetic drift is the change in the frequency of a gene in a population due to random sampling  and random extinctions that describe when a gene causes its carriers to have a deviating fitness from unity, its frequency will be determined by selection  in different basins. For protists in particular there is no literature available on this topic to our knowledge. At last, the Monoplization Hypothesis by De Meester et al.  could be relevant to protist biogeography stating that a fast population growth and local adaptation and colonization of a new habitat result in the monopolization of resources, which yields a strong priority effect. The effect is even enhanced when a locally adapted population can provide a ‘large resting propagule bank’ as a strong buffer against new genotypes invading. This holds true especially for species that reproduce asexually and form resting stages.
Even though mass effect and dispersal  cannot be ruled out, these are unlikely alternatives to explain the observed community patterns. The habitats of the water column above the DHABs represent a potential source habitat with ‘high quality’. In comparison, the narrow interphase and the brine show ‘low quality’ conditions because these habitats harbor high gradients of change, anoxia, high salt concentration up to saturation and therefore require a high degree of physiological adaptation for microbial colonization. Chances for highly specialized organisms to cross environmental barriers outside their habitat and to disperse beyond their specific habitat are very low. Evidence for this hypothesis comes, for example, from the very rare and infrequent dispersal of protists between marine and freshwater habitats [43, 79]. One possibility may be the dispersal of spores and/or cysts (resting stages), however, our knowledge about the number of ciliates that can form such resting stages in nature is very limited . Furthermore, physical mechanisms of transport for resting stages between different basins are difficult to imagine, considering the lack of fluid flow, high density, and lack of animal vectors in the brines. In contrast, this scenario may be more plausible for cysts/spores in halocline/interphase habitats. Physical transport of resting stages between haloclines at different basin sites could explain the observed similarities in ciliate interphase communities (Figure 3).
The deep basins in the eastern Mediterranean Sea may have recruited their protistan seed communities from Atlantic Sea water during the Zenclean Flood (~5.3 mya), when the Strait of Gibraltar opened permanently and refilled the mostly dried out Mediterranean Sea . Subsequently, due to the dissolution of evaporites and the rise of anoxia in deep basins the water masses became physically separated from each other. Anoxia and hydrochemistry likely exerted an increased pressure on the original protistan communities. Species sorting may have been driven through environmental filtering [37, 42, 62, 82]. This is a predictable and fundamental process of community assembly , that allows only those taxa with the genomic and physiological potential to cope with each specific set of environmental conditions. This has been evidenced for recent ciliate communities . The normsaline and normoxic deep-sea water separating the different hypersaline anoxic basins from each other then became an environmental barrier for most protists (with the exception of cyst-forming taxa), with the consequence that genetic exchange among the different brines was no longer likely. Changes in the SSU are presumably neutral, therefore, these changes would be due to random mutations. However, it is reasonable to assume that changes in the SSU rDNA are occurring in congruency with whole genome changes and not independent of evolutionary genome processes. Evolution over geological time may have resulted in significantly different ciliate communities in the brines. Divergence of species occurring in isolation through adaptive shifts that occurs in common seed species populations has been demonstrated for a number of taxa, including several macro- and microinvertebrates using molecular as well as taxonomic studies [84–87]. Based on our data, it is not unreasonable to assume that protists are also subjected to such evolutionary processes. Our study strongly suggests that evolutionary time scales combined with physical and hydrochemical isolation can explain, in part, the observed evolutionary differences in the ciliate communities in the different DHABs studied here.
The data presented here suggests that ancient isolated habitats, like deep-sea brine lakes, which also occur in the Red Sea  and the Gulf of Mexico  or habitats like aquifers and cave systems will represent ‘hot spots’ for the discovery of as-of-yet unstudied, possibly highly divergent, and endemic microbial eukaryotes.
Coordinates, sampling depths and physico-chemical data of the brines (B) and halocline interfaces (IF) of the different DHABs under study
Coordinates (Long, Lat)
22.312124 E, 34.19468 N
26.21962 E, 33.524236 N
22.084368 E, 34.401134 N
21.283252 E, 35.13528 N
22.312124 E, 34.19468 N
26.21962 E, 33.524236 N
22.084368 E, 34.401134 N
21.283252 E, 35.13528 N
Environmental RNA Isolation, transcription and PCR amplification of ciliate SSU rRNAs
Primer sets used in this study for the specific amplification of ciliate V4-SSU rRNA fragments using a two-step (nested) PCR reaction
1. Reaction: PCR ciliate specific primers
2. Reaction: PCR hyper variable V4-region
Pyrosequencing and sequence data processing
The DNA sequencing of the V4-amplicons was conducted by Engencore (University of South Carolina, USA) using Roche’s Titanium chemistry. One half plate was sequenced with the eight different samples with individual MIDs. The number of amplicons obtained after sequencing ranged between 33,634 (Thetis brine) and 80,650 (Urania interface) sequences.
For sequence data quality control and processing, we used the program JAguc . All tags that met any of the following conditions were considered as “low quality” and removed from further analyses: sequences <200 nucleotides, sequences containing an inaccurate calibration key, incomplete or erroneous forward and reverse primer sequences, presence of an ambiguity code. Sequences were then clustered. A cluster included sequences that shared at least 95% similarity in their primary structures. This conservative cluster threshold was chosen, because it accounts for sequencing errors and for intraspecific variability in the hypervariable SSU rDNA V4 region of ciliates [91, 92]. Single singletons (unique amplicons after 95% clustering that occurred exclusively in only one of the eight samples) were removed from downstream analyses as they are most likely erroneous sequencing products [91, 93].
We assigned taxonomy to each amplicon by conducting BLASTn searches implemented in JAguc (using parameters -m 7 -r 5 -q −4 -G 8 -E 6 -b 50) of each unique tag against a local installation of NCBI’s nucleotide database (nr/nt, release 187). Only unique tags with a best BLAST hit of at least 80% sequence similarity were assigned to a taxonomic category. The remaining tags were assigned to an artificial category “others”. This information was stored in JAguc’s database. We only assigned taxonomic labels to the genus level, because taxon assignments on lower taxonomic levels become inaccurate and biased due to the relatively limited sequence information provided in short amplicons . Taxonomy of ciliates follows the compendium “The Ciliated Protozoa” by D. Lynn .
Statistical analyses of ciliate amplicon profiles
To assess the ciliate diversity within a particular sample (alpha-diversity, ), we normalized the data (to the smallest number of sequences: 32,663 sequences were picked randomly 10,000 times in each of the samples with the software R ). We used the Shannon index (combining richness and relative abundance; ) and the non-parametric richness estimator ACE  as calculated with R . The partitioning of diversity between communities (beta-diversity) was calculated with the Bray-Curtis index, using the software EstimateS v.8 , and then translated into distance matrixes (1 minus Bray-Curtis index value) for UPGMA cluster analyses. Bray-Curtis similarity index is a modified version of the Sørensen index, which considers abundance distribution (also known as the Sørensen abundance Index or the quantitative Sørensen index [99, 100].
To assess an effect of distance on community similarities, Jaccard and Chao-Sørensen indices were plotted against distance data among individual sample sites in a Pearson-rank correlation using the Statistica software package. A Student’s t-test for paired samples was used for significance testing. A Mantel test between the geographic distance and the Bray Curtis distance matrices was conducted to evaluate the significance of the correlation coefficient between geographic and genetic distance. The Mantel test was conducted using the software add-in for Microsoft Excel XLSTAT (http://www.xlstat.com) with 10000 permutations. Geographical distances were calculated via the subtraction of different depths on a single geographical position, which resulted in the altitude difference within the same basin. For the calculation of the 2-dimensional great-circle distance between two points on a sphere from their longitudes and latitudes (same depth) the haversine formula  was implemented in the script as provided by Chris Veness (2002–2011) at http://www.movable-type.co.uk/scripts/latlong.html.
A canonical correspondence analysis (CCA) of quantitative amplicon profiles was conducted to describe the relationships between ciliate community composition patterns and underlying environmental gradients, which shape these diversity patterns. Data were log-transformed  and unconstrained permutations (n = 499) were run under a reduced model. Monte Carlo significance tests of first ordination axes and of all canonical axes together were performed. Initially, all available environmental variables (see above) were included in the model. In order to develop a robust model explaining as much variance as possible while avoiding multi-colinearity, individual variables were removed in a step-wise manner. We used the Canoco software (Microcomputer Power, Ithaca, NY, USA) for the ordination analysis.
Scanning electron microscopy (SEM) preparation and enumeration of ciliates
We used SEM to visualize ciliate morphotypes and to amend the molecular diversity survey with imaging analyses. We followed the method for SEM described in [25, 103]. In short, fixed samples were filtered onto 0.4-μm polycarbonate Transwell membrane filters (Corning, USA) and washed with 1X PBS (pH 7.4) that were taken through a dehydration series and fixed with 100% hexamethyldisilizane (Electron Microscopy Sciences, Hatfield, Pennsylvania) before air-drying. Transwell filters were not exposed to air at any point during the protocol, until the final step to prevent collapse of fixed protists. Filters were attached to a carbon adhesive tab and mounted on a SEM specimen holder. Mounted specimens were then sputter coated with 10–15 nm of gold and palladium (60:40) using a Tousimis Samsputter 2A and visualized with a Hitachi S4800 scanning electron microscope. A minimum of 50 microscopic fields (0.5 × 1.0 mm) were observed to count ciliates in each of the samples, and ciliates were counted on at least three different filters.
Fluorescence in situ hybridization
In order to evaluate the relative abundance of ciliates as part of the protistan assemblages we used fluorescence in situ hybridization with a specific oligonucleotide probe. FISH followed the protocol of . In short, 100–150 ml of paraformaldehyde-fixed (2% final concentration) seawater was filtered onto 0.65 μm filters and frozen at −20°C. Filters were thawed, and cut into small triangles before the hybridization step and ~20 μl of pre-heated 0.2% metaphor agarose were pipetted onto filters. After the metaphor agarose had dried, filter pieces were transferred to 0.5 ml sterile tubes containing the hybridization mix. All hybridizations were carried out using the universal eukaryotic FISH probe Euk1209  with 40% formamide for 2 hours at 46°C. After the hybridization step, filter pieces were incubated at 48°C in preheated washing buffer for 10 minutes in sterile 50 ml tubes. Filter pieces were then washed with distilled water and placed into sterile 0.5 ml tubes containing DAPI (2 μg/ml) and incubated for 5 minutes in the dark. Filter pieces were then washed with sterile water and incubated for 2 minutes in 70% ethanol, followed by a 2-minute wash with 100% ethanol. Filters were air dried and mounted on glass slides with a Citifluor/Vectashield mix (4:1) to prevent bleaching. Cells were enumerated under epifluorescence using a Zeiss Axioplan 2 microscope and photographed with a Hamamatsu digital camera.
The authors would like to thank Dr. Maria Pachiadaki, Dr. Matthias Engel and Melanie Müller for assistance with sample collection during the R/V Urania cruise, and the captains and crew of the R/V Urania and R/V Oceanus for their tireless assistance with sample collection. We thank Dominik Forster for help with R. This work was funded by NSF grants OCE-0849578 and OCE-1061774 to VE and support from Carl Zeiss fellowship to AS and from the Deutsche Forschungsgemeinschaft (grants STO414/3-2 and STO414/7-1) to TS.
- Ovreas L, Curtis TP: Microbial diversity and ecology. Biological Diversity: frontiers in measurement and assessment. Edited by: Magurran AE, McGill BJ. 2011, Oxford: Oxford University Press, 221-236.Google Scholar
- Alexander E, Stock A, Breiner HW, Behnke A, Bunge J, Yakimov MM, Stoeck T: Microbial eukaryotes in the hypersaline anoxic L'Atalante deep-sea basin. Environ Microbiol. 2009, 11: 360-381. 10.1111/j.1462-2920.2008.01777.x.PubMedGoogle Scholar
- Edgcomb V, Orsi W, Leslin C, Epstein S, Bunge J, Jeon SO, Yakimov MM, Behnke A, Stoeck T: Protistan community patterns within the brine and halocline of deep hypersaline anoxic basins in the eastern Mediterranean Sea. Extremophiles. 2009, 13: 151-167. 10.1007/s00792-008-0206-2.PubMedGoogle Scholar
- Camerlenghi A: Anoxic basins of the eastern Mediterranean: geological framework. Mar Chem. 1990, 31: 1-19. 10.1016/0304-4203(90)90028-B.Google Scholar
- La Cono V, Smedile F, Bortoluzzi G, Arcadi E, Maimone G, Messina E, Borghini M, Oliveri E, Mazzola S, L'Haridon S, et al.: Unveiling microbial life in new deep-sea hypersaline Lake Thetis. Part I: Prokaryotes and environmental settings. Environ Microbiol. 2011, 13 (8): 2250-2268. 10.1111/j.1462-2920.2011.02478.x.PubMedGoogle Scholar
- van der Wielen PW, Bolhuis H, Borin S, Daffonchio D, Corselli C, Giuliano L, D'Auria G, de Lange GJ, Huebner A, Varnavas SP, et al.: The enigma of prokaryotic life in deep hypersaline anoxic basins. Science. 2005, 307 (5706): 121-123. 10.1126/science.1103569.PubMedGoogle Scholar
- Azam F, Fenchel T, Field J, Gray J, Meyer-Reil L, Thingstad F: The ecological role of water column microbes in the sea. Mar Ecol Prog Ser. 1983, 10: 257-263.Google Scholar
- Corliss JO: Biodiversity and biocomplexity of the protists and an overview of their significant roles in maintenance of our biosphere. Acta Protozool. 2002, 41 (3): 199-220.Google Scholar
- Finlay BJ, Corliss JO, Esteban G, Fenchel T: Biodiversity at the microbial level: the number of free-living ciliates in the biosphere. Ouart Rev Biol. 1996, 71: 221-237. 10.1086/419370.Google Scholar
- Lynn DH, Gilron GL: A brief review of approaches using ciliated protists to assess aquatic ecosystem health. J Aquatic Ecosyst Health. 1992, 1: 263-270. 10.1007/BF00044168.Google Scholar
- Doherty M, Cosatas BA, McManus GB, Katz LA: Culture independent assessment of planktonic ciliate diversity in coastal northwest Atlantic waters. Aquat Microb Ecol. 2007, 48: 141-154.Google Scholar
- Fenchel T, Finlay BJ: The diversity of microbes: resurgence of the phenotype. Phil Trans Roy Soc Lond B Biol Sci. 2006, 361 (1475): 1965-1973. 10.1098/rstb.2006.1924.Google Scholar
- Finlay BJ: Global dispersal of free-living microbial eukaryote species. Science. 2002, 296 (5570): 1061-1063. 10.1126/science.1070710.PubMedGoogle Scholar
- Foissner W, Chao A, Katz LA: Diversity and geographic distribution of ciliates (Protista: Ciliophora). Biodiv Conserv. 2008, 17: 345-363. 10.1007/s10531-007-9254-7.Google Scholar
- Katz LA, McManus GB, Snoeyenbos-West OLO, Pirog K, Griffin A, Foissner W: Reframing the Microbial 'Everything is Everywhere' Debate: evidence for high gene flow and diversity in ciliate morphospecies. Aquat Microb Ecol. 2005, 41: 55-65.Google Scholar
- Stoeck T, Bass D, Nebel M, Christen R, Jones MD, Breiner HW, Richards TA: Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol Ecol. 2010, 19: 21-31.PubMedGoogle Scholar
- Weisse T: Distribution and diversity of aquatic protists: an evolutionary and ecological perspective. Biodiv Conserv. 2008, 17: 243-259. 10.1007/s10531-007-9249-4.Google Scholar
- Dunthorn M, Foissner W, Katz LA: Molecular phylogenetic analysis of class Colpodea (phylum Ciliophora) using broad taxon sampling. Mol Phylogenet Evol. 2008, 46 (1): 316-327. 10.1016/j.ympev.2007.08.006.PubMedGoogle Scholar
- Lynn DH: The Ciliated Protozoa. 2008, New York: Springer, ThirdGoogle Scholar
- Christen R: Global sequencing: a review of current molecular data and new methods available to assess microbial diversity. Microb Environ. 2008, 23 (4): 253-268. 10.1264/jsme2.ME08525.Google Scholar
- Epstein S, Lopez-Garcia P: “Missing” protists: a molecular prospective. Biodivers Conserv. 2008, online early (17): 261-276.Google Scholar
- Jeon S, Bunge J, Leslin C, Stoeck T, Hong S, Epstein SS: Environmental rRNA inventories miss over half of protistan diversity. BMC Microbiol. 2008, 8: 222-10.1186/1471-2180-8-222.PubMedPubMed CentralGoogle Scholar
- Moreira D, Lopez-Garcia P: The molecular ecology of microbial eukaryotes unveils a hidden world. Trends Microbiol. 2002, 10 (1): 31-38. 10.1016/S0966-842X(01)02257-0.PubMedGoogle Scholar
- Pedros-Alio C: Ecology. Dipping into the rare biosphere. Science. 2007, 315 (5809): 192-193. 10.1126/science.1135933.PubMedGoogle Scholar
- Orsi WD, Charvet S, Vdacny P, Bernhard JM, Edgcomb VP: Prevalence of partnerships between bacteria and ciliates in oxygen-depleted marine water columns. Front Microbiol. 2012, 3: 341-PubMedPubMed CentralGoogle Scholar
- Yetinson T, Shilo M: Seasonal and geographic distribution of luminous bacteria in the eastern mediterranean sea and the gulf of elat. Appl Environ Microbiol. 1979, 37 (6): 1230-1238.PubMedPubMed CentralGoogle Scholar
- Inagaki F, Nunoura T, Nakagawa S, Teske A, Lever M, Lauer A, Suzuki M, Takai K, Delwiche M, Colwell FS, et al.: Biogeographical distribution and diversity of microbes in methane hydrate-bearing deep marine sediments on the Pacific Ocean Margin. Proc Natl Acad Sci U S A. 2006, 103 (8): 2815-2820. 10.1073/pnas.0511033103.PubMedPubMed CentralGoogle Scholar
- Whitaker RJ, Grogan DW, Taylor JW: Geographic barriers isolate endemic populations of hyperthermophilic archaea. Science. 2003, 301 (5635): 976-978. 10.1126/science.1086909.PubMedGoogle Scholar
- Jones EBG, Pang KL: Tropical aquatic fungi. Biodiv Conserv. 2012, 21: 2403-2423. 10.1007/s10531-011-0198-6.Google Scholar
- Dolan JR: An introduction to the biogeography of aquatic microbes. Aquat Microb Ecol. 2005, 41 (1): 39-48.Google Scholar
- Martiny JBH, Bohannan BJM, Brown JH, Colwell RK, Fuhrman JA, Green JL, Horner-Devine MC, Kane M, Krumins JA, Kuske CR, et al.: Microbial biogeography: putting microorganisms on the map. Nat Rev Microbiol. 2006, 4 (2): 102-112. 10.1038/nrmicro1341.PubMedGoogle Scholar
- Vyverman W, Verleyen E, Sabbe K, Vanhoutte K, Sterken M, Hodgson DA, Mann DG, Juggins S, Van de Vijver B, Jones V, et al.: Historical processes constrain patterns in global diatom diversity. Ecology. 2007, 88 (8): 1924-1931. 10.1890/06-1564.1.PubMedGoogle Scholar
- Fenchel T, Esteban GF, Finlay BJ: Local versus global diversity of microorganisms: cryptic diversity of ciliated protozoa. Oikos. 1997, 80 (2): 220-225. 10.2307/3546589.Google Scholar
- Stephenson SL, Schnittler M, Novozhilov YK: Myxomycete diversity and distribution from the fossil record to the present. Biodivers Conserv. 2008, 17 (2): 285-301. 10.1007/s10531-007-9252-9.Google Scholar
- Wilson DS: Complex interactions in metacommunities, with implications for biodiversity and higher levels of selection. Ecology. 1992, 73: 1984-2000. 10.2307/1941449.Google Scholar
- Leibold MA, Holyoak M, Moquet N, Amarasekare P, Chase JM, Hoopes MF, Holt RD, Shurin JB, Law R, Tilman D, Loreau M, Gonzalez A: The metacommunity concept: a framework for multi-scale community ecology. Ecol Lett. 2004, 7: 601-613. 10.1111/j.1461-0248.2004.00608.x.Google Scholar
- Holyoak M, Leibold MA, Holt RD: Metacommunities: Spatial Dynamics and Ecological Communities. 2005, Chicago, IL, USA: The University of Chicago PressGoogle Scholar
- Santangelo G, Lucchesi P: Spatial distribution pattern of ciliated protozoa in a Mediterranean interstitial environment. Aquat Microb Ecol. 1995, 9: 47-54.Google Scholar
- Albuquerque L, Taborda M, La Cono V, Yakimov M, da Costa MS: Natrinema salaciae sp. nov., a halophilic archaeon isolated from the deep, hypersaline anoxic Lake Medee in the Eastern Mediterranean Sea. Syst Appl Microbiol. 2012, 35 (6): 368-373. 10.1016/j.syapm.2012.06.005.PubMedGoogle Scholar
- Forster D, Behnke A, Stoeck T: Meta-analyses of environmental sequence data identify anoxia and salinity as parameters shaping ciliate communities. Systematics and Biodiversity. 2012, 10 (3): 277-288. 10.1080/14772000.2012.706239.Google Scholar
- Lozupone CA, Knight R: Global patterns in bacterial diversity. Proc Natl Acad Sci U S A. 2007, 104 (27): 11436-11440. 10.1073/pnas.0611525104.PubMedPubMed CentralGoogle Scholar
- Logares R, Lindstrom ES, Langenheder S, Logue JB, Paterson H, Laybourn-Parry J, Rengefors K, Tranvik L, Bertilsson S: Biogeography of bacterial communities exposed to progressive long-term environmental change. ISME J. 2013, 7 (5): 937-948. 10.1038/ismej.2012.168.PubMedPubMed CentralGoogle Scholar
- Logares R, Brate 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.PubMedGoogle Scholar
- Oren A, Larimer F, Richardson P, Lapidus A, Csonka LN: How to be moderately halophilic with broad salt tolerance: clues from the genome of Chromohalobacter salexigens. Extremophiles. 2005, 9 (4): 275-279. 10.1007/s00792-005-0442-7.PubMedGoogle Scholar
- Ramos-Cormenzana A: Halophilic organisms and their environment. General and Applied Aspects of Halophilic Microorganisms. Edited by: Rodriguez-Valera F. 1991, New York: Plenum Press, 15-24.Google Scholar
- Pedros-Alio C, Calderon-Paz JI, MacLean MH, Medina G, Marrase C, Gasol JM, Guixa-Boixereu N: The microbial food web along salinity gradients. FEMS Microbiol Ecol. 2000, 32 (2): 143-155.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Finlay BJ, Esteban GF, Brown S, Fenchel T, Hoef-Emden K: Multiple cosmopolitan ecotypes within a microbial eukaryote morphospecies. Protist. 2006, 157 (4): 377-390. 10.1016/j.protis.2006.05.012.PubMedGoogle Scholar
- Embley TM, Finlay BJ, Dyal PL, Hirt RP, Wilkinson M, Williams AG: Multiple Origins of Anaerobic Ciliates with Hydrogenosomes within the Radiation of Aerobic Ciliates. Phil Trans Roy Soc Lond B Biol Sci. 1995, 262 (1363): 87-93.Google Scholar
- Hjort K, Goldberg AV, Tsaousis AD, Hirt RP, Embley TM: Diversity and reductive evolution of mitochondria among microbial eukaryotes. Phil Trans Roy Soc Lond B Biol Sci. 2010, 365 (1541): 713-727. 10.1098/rstb.2009.0224.Google Scholar
- Boxma B, de Graaf RM, van der Staay GWM, van Alen TA, Ricard G, Gabaldon T, van Hoek AHAM, der Staay SY M-v, Koopman WJH, van Hellemond JJ, et al.: An anaerobic mitochondrion that produces hydrogen. Nature. 2005, 434 (7029): 74-79. 10.1038/nature03343.PubMedGoogle Scholar
- Fenchel T, Perry T, Thane A: Anaerobiosis and symbiosis with bacteria in free-living ciliates. J Eukaryot Microbiol. 1977, 24: 154-163. 10.1111/j.1550-7408.1977.tb05294.x.Google Scholar
- van Hoek AH, van Alen TA, Sprakel VS, Leunissen JA, Brigge T, Vogels GD, Hackstein JH: Multiple acquisition of methanogenic archaeal symbionts by anaerobic ciliates. Mol Biol Evol. 2000, 17 (2): 251-258. 10.1093/oxfordjournals.molbev.a026304.PubMedGoogle Scholar
- Edgcomb V, Orsi W, Breiner HW, Stock A, Filker S, Yakimov MM, Stoeck T: Novel active kinetoplastids associated with hypersaline anoxic basins in the Eastern Mediterranean deep-sea. Deep-Sea Res I. 2011, 58: 1040-1048. 10.1016/j.dsr.2011.07.003.Google Scholar
- Stoeck T, Taylor GT, Epstein SS: Novel eukaryotes from the permanently anoxic Cariaco Basin (Caribbean Sea). Appl Environ Microbiol. 2003, 69 (9): 5656-5663. 10.1128/AEM.69.9.5656-5663.2003.PubMedPubMed CentralGoogle Scholar
- Behnke A, Bunge J, Barger K, Breiner HW, Alla V, Stoeck T: Microeukaryote community patterns along an O2/H2S gradient in a supersulfidic anoxic Fjord (Framvaren, Norway). Appl Environ Microbiol. 2006, 72 (5): 3626-3636. 10.1128/AEM.72.5.3626-3636.2006.PubMedPubMed CentralGoogle Scholar
- Zuendorf A, Behnke A, Bunge J, Barger K, Stoeck T: Diversity estimates of microeukaryotes below the chemocline of the anoxic Mariager Fjord, Denmark. FEMS Microbiol Ecol. 2006, 58: 476-491. 10.1111/j.1574-6941.2006.00171.x.PubMedGoogle Scholar
- Stock A, Jurgens K, Bunge J, Stoeck T: Protistan diversity in suboxic and anoxic waters of the Gotland Deep (Baltic Sea) as revealed by 18S rRNA clone libraries. Aquat Microb Ecol. 2009, 55 (3): 267-284.Google Scholar
- Wylezich C, Jurgens K: Protist diversity in suboxic and sulfidic waters of the Black Sea. Environ Microbiol. 2011, 13 (11): 2939-2956. 10.1111/j.1462-2920.2011.02569.x.PubMedGoogle Scholar
- Casamayor EO, Garcia-Cantizano J, Pedros-Alio C: Carbon dioxide fixation in the dark by photosynthetic bacteria in sulfide-rich stratified lakes with oxic-anoxic interfaces. Limnol Oceanogr. 2008, 53 (4): 1193-1203. 10.4319/lo.2008.53.4.1193.Google Scholar
- Oren A: Thermodynamic limits to microbial life at high salt concentrations. Environ Microbiol. 2011, 13 (8): 1908-1923. 10.1111/j.1462-2920.2010.02365.x.PubMedGoogle Scholar
- Rengefors K, Logares R, Laybourn-Parry J: Polar lakes may act as ecological islands to aquatic protists. Mol Ecol. 2012, 21 (13): 3200-3209. 10.1111/j.1365-294X.2012.05596.x.PubMedPubMed CentralGoogle Scholar
- Pernthaler J: Predation on prokaryotes in the water column and its ecological implications. Nat Rev Microbiol. 2005, 3 (7): 537-546. 10.1038/nrmicro1180.PubMedGoogle Scholar
- Yoon HS, Price DC, Stepanauskas R, Rajah VD, Sieracki ME, Wilson WH, Yang EC, Duffy S, Bhattacharya D: Single-cell genomics reveals organismal interactions in uncultivated marine protists. Science. 2011, 332 (6030): 714-717. 10.1126/science.1203163.PubMedGoogle Scholar
- Coolen MJ: 7000 years of Emiliania huxleyi viruses in the Black Sea. Science. 2011, 333 (6041): 451-452. 10.1126/science.1200072.PubMedGoogle Scholar
- Miki T, Jacquet S: Complex interactions in the microbial world: underexplored key links between viruses, bacteria and protozoan grazers in aquatic environments. Aquat Microb Ecol. 2008, 51 (2): 195-208.Google Scholar
- Verity PG: Feeding in planktonic protozoans, evidence for non-random acquisition of prey. J Protozool. 1991, 38: 69-76.Google Scholar
- Yakimov MM, Giuliano L, Cappello S, Denaro R, Golyshin PN: Microbial community of a hydrothermal mud vent underneath the deep-sea anoxic brine lake Urania (eastern Mediterranean). Orig Life Evol Biosph. 2007, 37 (2): 177-188. 10.1007/s11084-006-9021-x.PubMedGoogle Scholar
- Yakimov MM, La Cono V, Denaro R, D'Auria G, Decembrini F, Timmis KN, Golyshin PN, Giuliano L: Primary producing prokaryotic communities of brine, interface and seawater above the halocline of deep anoxic lake L'Atalante, Eastern Mediterranean Sea. ISME J. 2007, 1 (8): 743-755. 10.1038/ismej.2007.83.PubMedGoogle Scholar
- Beardsley C, Pernthaler J, Wosniok W, Amann R: Are Readily Culturable Bacteria in Coastal North Sea Waters Suppressed by Selective Grazing Mortality?. Appl Environ Microbiol. 2003, 69 (5): 2624-2630. 10.1128/AEM.69.5.2624-2630.2003.PubMedPubMed CentralGoogle Scholar
- Matz C, Boenigk J, Arndt H, Jurgens K: Role of bacterial phenotypic traits in selective feeding of the heterotrophic nanoflagellate Spumella sp. Aquat Microb Ecol. 2002, 27 (2): 137-148.Google Scholar
- Gonzalez JM, Sherr EB, Sherr BF: Size-Selective Grazing on Bacteria by Natural Assemblages of Estuarine Flagellates and Ciliates. Appl Environ Microbiol. 1990, 56 (3): 583-589.PubMedPubMed CentralGoogle Scholar
- Simek K, Vrba J, Hartman P: Size-Selective Feeding by Cyclidium Sp on Bacterioplankton and Various Sizes of Cultured Bacteria. FEMS Microbiol Ecol. 1994, 14 (2): 157-167.Google Scholar
- James JW: The founder effect and response to artificial selection. Genet Res. 1970, 16 (3): 241-250. 10.1017/S0016672300002500.PubMedGoogle Scholar
- Masel J: Genetic drift. Curr Biol. 2011, 21 (20): R837-R838. 10.1016/j.cub.2011.08.007.PubMedGoogle Scholar
- Fisher RA: The genetical theory of natural selection. 1930, Oxford: The Clarendon PressGoogle Scholar
- De Meester L, Gómez A, Okamura B, Schwenk K: The Monopolization Hypothesis and the dispersal–gene flow paradox in aquatic organisms. Acta Oecol. 2002, 23: 121-135. 10.1016/S1146-609X(02)01145-1.Google Scholar
- Urban MC, Leibold MA, Amarasekare P, De Meester L, Gomulkiewicz R, Hochberg ME, Klausmeier CA, Loeuille N, de Mazancourt C, Norberg J, et al.: The evolutionary ecology of metacommunities. Trends Ecol Evol. 2008, 23 (6): 311-317. 10.1016/j.tree.2008.02.007.PubMedGoogle Scholar
- Brate 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 rDNA. ISME J. 2010, 4 (9): 1144-1153. 10.1038/ismej.2010.39.PubMedGoogle Scholar
- Foissner W: Biogeography and dispersal of micro-organisms: A review emphasizing protists. Acta Protozool. 2006, 45 (2): 111-136.Google Scholar
- Garcia-Castellanos D, Estrada F, Jimenez-Munt I, Gorini C, Fernandez M, Verges J, De Vicente R: Catastrophic flood of the Mediterranean after the Messinian salinity crisis. Nature. 2009, 462 (7274): 778-U796. 10.1038/nature08555.PubMedGoogle Scholar
- Whittaker RH: Classification of natural communities. Bot Rev. 1962, 28: 1-239. 10.1007/BF02860872.Google Scholar
- Lebrija-Trejos E, Perez-Garcia EA, Meave JA, Bongers F, Poorter L: Functional traits and environmental filtering drive community assembly in a species-rich tropical system. Ecology. 2010, 91 (2): 386-398. 10.1890/08-1449.1.PubMedGoogle Scholar
- Humphreys WF, Watts CHS, Cooper SJB, Leijs R: Groundwater estuaries of salt lakes: buried pools of endemic biodiversity on the western plateau, Australia (vol 626, pg 79, 2009). Hydrobiologia. 2009, 632 (1): 377-10.1007/s10750-009-9864-z.Google Scholar
- Juan C, Guzik MT, Jaume D, Cooper SJ: Evolution in caves: Darwin's 'wrecks of ancient life' in the molecular era. Mol Ecol. 2010, 19 (18): 3865-3880. 10.1111/j.1365-294X.2010.04759.x.PubMedGoogle Scholar
- Leijs R, van Nes EH, Watts CH, Cooper SJB, Humphreys WF, Hogendoorn K: Evolution of Blind Beetles in Isolated Aquifers: a Test of Alternative Modes of Speciation. PLoS One. 2012, 7 (3): e34260-10.1371/journal.pone.0034260.PubMedPubMed CentralGoogle Scholar
- Leys R, Watts CH, Cooper SJ, Humphreys WF: Evolution of subterranean diving beetles (Coleoptera: Dytiscidae: Hydroporini, Bidessini) in the arid zone of Australia. Evolution. 2003, 57 (12): 2819-2834.PubMedGoogle Scholar
- Degens ET, Ross AT: Hot Brines and Recent Heavy Metal Deposits in the Red Sea. A Geochemical and Geophysical Account. 1969, Berlin/Heidelberg/New York: SpringerGoogle Scholar
- Shokes RF, Trabant PK, Presley BJ, Reid DF: Anoxic, hypersaline basin in the northern gulf of Mexico. Science. 1977, 196 (4297): 1443-1446. 10.1126/science.196.4297.1443.PubMedGoogle Scholar
- Nebel ME, Wild S, Holzhauser M, Huttenberger L, Reitzig R, Sperber M, Stoeck T: JAGUC–a software package for environmental diversity analyses. J Bioinform Comput Biol. 2011, 9 (6): 749-773. 10.1142/S0219720011005781.PubMedGoogle Scholar
- Behnke A, Engel M, Christen R, Nebel M, Klein RR, Stoeck T: Depicting more accurate pictures of protistan community complexity using pyrosequencing of hypervariable SSU rRNA gene regions. Environ Microbiol. 2011, 13 (2): 340-349. 10.1111/j.1462-2920.2010.02332.x.PubMedGoogle Scholar
- Dunthorn M, Klier J, Bunge J, Stoeck T: Comparing the Hyper-Variable V4 and V9 Regions of the Small Subunit rDNA for Assessment of Ciliate Environmental Diversity. J Eukaryot Microbiol. 2012, 59 (2): 185-187. 10.1111/j.1550-7408.2011.00602.x.PubMedGoogle Scholar
- Kunin V, Engelbrektson A, Ochman H, Hugenholtz P: Wrinkles in the rare biosphere: pyrosequencing errors can lead to artificial inflation of diversity estimates. Environ Microbiol. 2010, 12 (1): 118-123. 10.1111/j.1462-2920.2009.02051.x.PubMedGoogle Scholar
- Whittaker RH: Evolution and measurements of species diversity. Taxon. 1972, 21: 213-251. 10.2307/1218190.Google Scholar
- Team RC: R: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical Computing,http://www.R-project.org/, 3-900051-07-0,Google Scholar
- Shannon CE: A mathematical theory of communication. The Bell System Technical Journal. 1948, 27: 379-423. 623–656Google Scholar
- Chao A, Lee SM, Jeng SL: Estimating population size for capture-recapture data when capture probabilities vary by time and individual animal. Biometrics. 1992, 48 (1): 201-216. 10.2307/2532750.PubMedGoogle Scholar
- Colwell RK: EstimateS: Statistical estimation of species richness and shared species from samples. Version 8.2. User's Guide and application. 2009,http://viceroy.eeb.uconn.edu/estimates,Google Scholar
- Bray RJ, Curtis JT: An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957, 27: 325-349. 10.2307/1942268.Google Scholar
- Magurran AE: Measuring biological diversity. 2004, Oxford: Blackwell PublishingGoogle Scholar
- Sinnott RW: Virtues of the Haversine. Sky Telescope. 1984, 68: 1-159.Google Scholar
- Grant A, Ogilvie LA: Terminal restriction fragment length polymorphism data analysis. Appl Environ Microbiol. 2003, 69 (10): 6342-10.1128/AEM.69.10.6342-6343.2003. author reply 6342–6343PubMedPubMed CentralGoogle Scholar
- Edgcomb V, Leadbeater ER, Bourland W, Beaudoin D, Bernhard JM: Structured multiple endosymbiosis of bacteria and archaea in a ciliate from marine sediments: a survival mechanism in low oxygen, sulfidic sediments?. Front Microb Physiol Metabol. 2011, 2: 55-Google Scholar
- Stoeck T, Fowle WH, Epstein SS: Methodology of protistan discovery: from rRNA detection to quality scanning electron microscope images. Appl Environ Microbiol. 2003, 69 (11): 6856-6863. 10.1128/AEM.69.11.6856-6863.2003.PubMedPubMed CentralGoogle Scholar
- Lara E, Berney C, Harms H, Chatzinotas A: Cultivation-independent analysis reveals a shift in ciliate 18S rRNA gene diversity in a polycyclic aromatic hydrocarbon-polluted soil. FEMS Microbiol Ecol. 2007, 62 (3): 365-373. 10.1111/j.1574-6941.2007.00387.x.PubMedGoogle Scholar