- Research article
- Open Access
Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean
BMC Microbiology volume 17, Article number: 192 (2017)
Benthic ciliates and the environmental factors shaping their distribution are far from being completely understood. Likewise, deep-sea systems are amongst the least understood ecosystems on Earth. In this study, using high-throughput DNA sequencing, we investigated the diversity and community composition of benthic ciliates in different sediment layers of a seamount and an adjacent abyssal plain in the tropical Western Pacific Ocean with water depths ranging between 813 m and 4566 m. Statistical analyses were used to assess shifts in ciliate communities across vertical sediment gradients and water depth.
Nine out of 12 ciliate classes were detected in the different sediment samples, with Litostomatea accounting for the most diverse group, followed by Plagiopylea and Oligohymenophorea. The novelty of ciliate genetic diversity was extremely high, with a mean similarity of 93.25% to previously described sequences. On a sediment depth gradient, ciliate community structure was more similar within the upper sediment layers (0-1 and 9-10 cm) compared to the lower sediment layers (19-20 and 29-30 cm) at each site. Some unknown ciliate taxa which were absent from the surface sediments were found in deeper sediments layers. On a water depth gradient, the proportion of unique OTUs was between 42.2% and 54.3%, and that of OTUs shared by all sites around 14%. However, alpha diversity of the different ciliate communities was relatively stable in the surface layers along the water depth gradient, and about 78% of the ciliate OTUs retrieved from the surface layer of the shallowest site were shared with the surface layers of sites deeper than 3800 m. Correlation analyses did not reveal any significant effects of measured environmental factors on ciliate community composition and structure.
We revealed an obvious variation in ciliate community along a sediment depth gradient in the seamount and the adjacent abyssal plain and showed that water depth is a less important factor shaping ciliate distribution in deep-sea sediments unlike observed for benthic ciliates in shallow seafloors. Additionally, an extremely high genetic novelty of ciliate diversity was found in these habitats, which points to a hot spot for the discovery of new ciliate species.
Deep-sea systems (> 1000 m) cover more than 65% of the Earth’s surface and fulfill a range of key ecosystem functions . Yet, they are amongst the least understood ecosystems on Earth . Among various deep-sea habitats, seamounts are widespread and prominent features of the world’s underwater landscape [3, 4]. Their specific topography creates distinct habitats, characterized by particular hydrography and substrate types, which influence the diversity of the benthos . So far, most biological research related to seamounts focused on patterns of mega- and macrobenthic biodiversity and their biogeography . Information about the diversity and distribution of benthic protists from seamount sediments are scarce, with the exception of benthic foraminifera, whose seamount communities have been described by several authors using morphological criteria . In the past decade, molecular techniques and sequencing the 18S rDNA as a taxonomic marker were used to investigate the diversity of protists in deep-sea sediments. These studies did not only reveal high protistan diversities in various deep-sea habitats [7, 8], but also numerous undescribed taxa and even several early branching eukaryotic lineages . In some of these habitats, ciliates were the predominant and often the most diverse microeukaryotic group [10,11,12,13]. However, knowledge on their distribution patterns and the major factors enforcing their dispersal in deep-sea sediments are far from complete.
In a microscopy study, water depth emerged as an important factor structuring the distribution of benthic protists, including ciliates, in the deep sea . However, because protists are not highly abundant in the deep sea , it is not surprising that microscopy studies miss a large proportion of the deep-sea protist diversity. Accordingly, using molecular techniques, more than 125 ciliate OTUs were detected in less than 1 g of abyssal sediments . This exceeded the number of ciliate OTUs detected in the coastal sediments . Thus, it is reasonable to assume that a large proportion of ciliate diversity in the deep-sea is still unknown to science.
Moreover, studies of benthic protists mainly focused on the top 1-2 cm surface layer of deep-sea sediments . However, it is well known that different ciliate species are not equally distributed along a vertical sediment gradient, mostly owed to oxygen, organic carbon and grain-size gradients [17,18,19,20]. Thus, some unknown ciliate taxa which are absent from the surface sediments are likely to be found in deeper sediment layers.
In this study, we collected sediment samples from four sites in the tropical Western Pacific Ocean with water depths ranging between 813 m and 4566 m as well as from different sediment layers. This sampling strategy was to maximize the number of different deep-sea habitats and thus, the proportion of novel ciliate diversity to be found. These habitats included a seamount and the adjacent abyssal plain. Statistical analyses were then used to assess shifts in ciliate communities across vertical sediment gradients and water depth.
Environmental parameters of sampling sites
For each sediment sample, the median grain size, the proportion of each sediment component and the total organic carbon (TOC) content were analyzed. The median grain size (mgs) of sediments differed notably among the different layers and sites and ranged between 4.5 * 10−13 μm (DS2.2) and 799.4 μm (S-B.1; Table 1). The mgs of surface sediments was highest at all sampling sites, but site DS1. Here, the highest mgs was detected at 20 cm depth. Sediment composition was different for each sampling site (Table 1). The TOC content of the sediments at sites S-M, S-B and DS1 varied between 0.1% (S-M.1) and 0.82% (DS1.3). The TOC content at site DS2 was, in contrast, more than 8 times higher than at the other sampling sites (e.g. DS2.2: 7.95%; Table 1).
Overview of the sequencing data
For the 12 sediment samples analyzed in this study, we obtained a total of 640,860 high-quality ciliate V4 sequences, which clustered into 104 distinct operational taxonomic units (OTUs) based on a 97% sequence similarity. The number of reads per sample varied between 25,169 (DS1.2) and 137,990 (S-B.2), with an average of 56,820 reads. The total number of OTUs obtained from each sample and used in downstream analyses ranged between 12 (S-B.2) and 55 (S-M.3), with an average of 33 OTUs (Table 2). Rarefaction analyses indicated near-saturated sampling for all samples except of S-B.1 (Additional file 1).
Distribution of ciliate OTUs along water depth gradients and sediment layers
Observed benthic ciliate alpha diversity varied remarkably among the different sites and sediment depths. The effective number of species was on average highest at sampling site S-M and lowest at site S-B (S-M: 13 ± 1.9; DS2: 9.9 ± 2.4; DS1: 3.1 ± 0.8; S-B: 2.8 ± 1.1), whereas OTU richness was on average highest at site S-M and lowest at DS1 (S-M: 46.3 ± 7.1; DS2: 41 ± 0.0; S-B: 21.7 ± 12.9, DS1: 21.3 ± 7.5; Table 2). The number of OTUs in surface sediments along water depth was relatively stable and ranged from 33 to 41. Except of sampling site S-M, the effective number of species and the OTU richness were higher in the surface sediment samples than in the deeper layers (Table 2).
About 24% of the ciliate OTUs detected in this study were unique to one sample. No single OTU was present in all the 12 samples, whereas 20 OTUs (19%) were detected at more than six sampling sites. Abundant OTUs in the dataset accounted for 20% (DS2.2) to 45% (S-M.2), whereas rare OTUs contributed 26% (S-M.2) to 69% (DS1.2; Table 2).
At a local scale (i.e. sediment depth gradient), the major proportion of OTUs was unique to one sample and accounted for 30.1% (S-M) to 75% (S-B), whereas OTUs shared by 4 (S-M), 3 (S-B, DS1) or 2 (DS2) samples only contributed 1.9% (S-B) to 46.4% (DS2) to the ciliate OTU distribution (Fig. 1 a-d). Also on a water depth gradient, the proportion of unique OTUs was between 42.2% (Fig. 1 f) and 54.3% (Fig. 1 e) and that of shared OTUs (by 4 sites) only around 14%. 29 OTUs (78%) of the surface layer at the shallow site S-M were shared with the surface layers at sites deeper than 3800 m (Fig. 1 e).
Partitioning of diversity revealed that the ciliate community structure was more similar within the upper sediment layers (0-1 and 9-10 cm) compared to the lower sediment layers (19-20 and 29-30 cm) at each site. Ciliate communities in the upper sediment layers (0-1 and 9-10 cm) of sites DS1, DS2 and S-M grouped together, respectively, with a high mean Jaccard similarity indicating a high stability of these clusters (Fig. 2 a). Ciliate communities from site S-B are distinct to the other ciliate communities and do not group into the clustering pattern (Fig. 2 a). Taking sequence abundances into account, this pattern is confirmed, although the overall grouping of the samples is slightly different (Additional file 2 a).
Ciliate diversity, i.e. OTU richness and effective number of species, and ciliate community composition in the surface sediments were not significantly correlated with any of the measured environmental parameters (Additional file 3). Additionally, changes of ciliate community composition in the deep-sea sediments were not significantly related to any combination of the environmental factors. The Mantel test showed that the pairwise community dissimilarity had no significant correlation with the geographic distance at a 350 km distance scale (r = −0.432, p = 0.83).
Taxonomic classification of ciliate OTUs
Ciliate OTUs were assigned to 9 (out of 12) ciliate classes (Fig. 2 b; Additional file 2 b), originating from 40 families and 51 genera (Additional file 4). The classes Armophorea, Cariacotrichea and Colpoda were not detected. Most OTUs belonged to the class Litostomatea (on average 23.3% of total OTUs at each site). The Litostomatea accounted for up to 33.3% (S-M.4: 17 OTUs) of the total ciliate OTU diversity in the investigated sediment samples. However, litostomatean sequences contributed only 4.7% to the total ciliate sequence abundance (up to 15.6% in sample S-M.2; Additional file 2 b). All of these OTUs were affiliated to 8 known genera. Within the Litostomatea, sequences related to Phialina, Loxophyllum and Litonotus were most abundant (Additional file 4). No litostomatean sequences could be retrieved from sample DS1.3.
Plagiopylea were the second most diverse ciliate group, contributing on average 18.6% to the OTU diversity at each site, but 47.8% to the total sequence abundance. Plagiopylean ciliates were detected at 11 of 12 sites, with the exception of site S-B.3. The proportion of plagiopylean OTUs varied between 0 (S-B.3) and 26.8% (DS2.1: 11 OTUs) and the proportion of sequences between 0 (S-B.3) and 78% (S-B.2) (Fig. 2 b, Additional file 2 b). All OTUs/sequences were related to the anaerobic ciliate Epalxella (Additional file 4). The sequence similarities of these OTUs to Epalxella ranged between 85% and 94%.
Oligohymenophorean and spirotrichean OTUs accounted for 16% on average of all detected OTUs at each site, respectively. The proportion of OTUs affiliated to Oligohymenophorea ranged between 4.9% (S-B.3: 2 OTUs) and 20% (S-M.3: 11 OTUs) (Fig. 2 b). Oligohymenophorean sequences accounted for 35.4% of all ciliate sequences, with the lowest proportion of sequences detected in sample DS1.3 (3%) and the highest in sample S-B.3 (90%) (Additional file 2 b). All oligohymenophorean OTUs/sequences were assigned to 10 known genera and one unclassified taxon. Within this class, sequences related to Pseudocyclidium, Trichodina and Pleuronema were most abundant. The most cosmopolitan oligohymenophorean OTU, which could be detected in 11 of 12 samples, was closely related to Pleuronema setigerum with the similarity of 99% (Additional file 4).
The proportion of spirotrichean OTUs ranged between 4.9% (DS2.2) and 37.5% (S-M.3) in the different benthic ciliate communities, all of these OTUs were affiliated to 15 known genera (Additional file 4). Spirotrichean sequences contributed 8.9% to the total number of ciliate sequences, the sequence proportion accounted for up to 35.6% in sample DS1.3 (Additional file 2 b). Within this class, sequences related to the genera Amphisiella, Apokeronopsis and Aspidisca were most abundant. Ten of 21 OTUs, which were observed in six or more samples, belonged to Spirotrichea (Additional file 4), of which five were affiliated to the order Urostylida. The OTUs affiliated to oligotrich ciliates of Novistrombidium orientale and Parallelostrombidium obesum were detected in the down-core sediments only (Additional file 4).
Classes Karyorelictea (on average 8.8% of total number of OTUs in each sample), Prostomatea (5.3%), Phyllopharyngea (4.3%), Nassophorea (4%) and Heterotrichea (3.8%) were relatively rare (Fig. 2 b, Additional file 2 b) at all sampling sites.
Degree of novel diversity
The observation of numerous low-identity (< 97%) OTUs in our dataset points to a high genetic novelty within the deep-sea habitats. The mean similarity of all detected ciliate OTUs to previously described sequences was only 93.25%. 73% of OTUs had a sequence similarity of less than 97% to reference sequences and 25% of the OTUs had an identity match of less than 90%. The mean similarities for each sample were relative stable and ranged between 92.86% (S-M.1) and 94.55% (DS1.1). The lowest overall mean sequence similarities, i.e. the highest degree of novel diversity occurred within the class Plagiopylea (89%). All of the 18 plagiopylean OTUs had a sequence similarity < 95%. The highest overall mean sequence similarities were observed within the classes Spirotrichea and Karyorelictea (95.7%, respectively; Fig. 3). Eleven of 23 spirotrichean OTUs and 4 of 8 karyorelictean OTUs had a sequence similarity of equal or more than 97%.
Comparison of benthic ciliate communities from different marine habitats
OTU richness peaked in the continental shelf samples (mean OTU richness: 105.8 ± 44), followed by the intertidal samples (mean OTU richness: 93.3 ± 25) and the hydrothermal vent samples (mean OTU richness: 40 ± 21). Alpha diversity was lowest in the seamount and abyssal plain samples (mean OTU richness: 34 ± 4). Changes in alpha diversity were significant between ciliate communities from the intertidal zone (p = 0.01) or continental shelf (p = 0.01), respectively, and the communities of the seamount and abyssal plain. No significant difference was detected between the OTU richness of the seamount and abyssal plain communities and the hydrothermal vent communities.
The proportion of spirotrichean OTUs was highest in the communities from the intertidal zone and continental shelf (Fig. 4). Almost no Plagiopylea ciliate was detected in the coastal samples, whereas Plagiopylea was one of the most diverse ciliate group in the seamount, abyssal plain, and hydrothermal vent samples (Fig. 4). The proportion of oligohymenophorean OTUs was the highest in the communities of the hydrothermal vent samples. UPGMA clustering based on the Sorensen-Dice similarity coefficient between the ciliate communities revealed two large clusters, separating intertidal and offshore communities from the deep-sea seamount, abyssal plain and hydrothermal vent communities (Fig 4). ANOSIM analyses confirmed this observation showing a significant disparity among different habitat types (r = 0.92; p = 0.001).
Ciliate diversity and their distribution along the water depth gradient
Our study revealed the diversity of ciliates along a wide water depth gradient (813-4566 m) in the sediments from a seamount and an adjacent abyssal plain in the Western Pacific Ocean. Although sample saturation was reached for all but one sample, we found a notably lower number of distinct ciliate OTUs (33-41 OTUs at each site) in the surface sediment samples analyzed in this study compared to a previous molecular diversity survey of protists in deep-sea sediments (68-172 ciliate OTUs at each site) . This could be partially explained by methodological differences, such as the use of another barcode marker (V9 vs. V4 region of SSU rDNA) and/or the deletion of singletons (OTUs containing only one sequence) in our study, which aimed to improve specificity at the cost of a possible loss in sensitivity. Instead, our results corroborate well with a microscopy study , suggesting a lower diversity of ciliates in deep-sea sediments compared to that in sediments from an intertidal zone and a continental shelf (e.g. 58-115 OTUs in the Yellow Sea and the East China Sea [21, 22]; 77-365 OTUs in the European coastal regions ). However, if water depth is deep enough, e.g. 813 m or more than 3800 m, it seems that water depth has little to no effect on the benthic ciliate OTU richness. This is different from the conclusion of the investigation in coastal sediments, where water depth overrides both geographic distance and environmental heterogeneity in determining the community composition of protists . Likewise, ciliate biomass and the number of species showed a significantly negative correlation with water depth in the offshore sediments of the Yellow Sea . However, little variation was observed in microeukaryotic communities from abyssal depths ranging between 5033 m and 5655 m . Likewise, using the data provided by the study of Pawlowski et al. for water depths ranging between 2292 and 6326 m, we could not detect a significant correlation between the benthic ciliate OTU richness and water depth . This indicates that water depth is likely more important on influencing the distribution of microbial eukaryotes in coastal sediments than in deep-sea sediments.
In our study, other factors than water depth might be more relevant for the structuring of ciliates in deep-sea sediments, e.g. oxygen saturation, which is a decisive variable for shaping protist community in shallow water sediments [17,18,19,20]. Unfortunately, this parameter could not be measured in our study.
Sediment patchiness at the deep-sea bottom is extremely high , possibly explaining the high dissimilarity of ciliate communities between all samples and the low number of shared OTUs between samples. Since rarefaction analyses confirm near-saturated sampling profiles for all but one sample, undersampling is no reason for the observed differences in ciliate communities.
Previous studies mainly focused on high taxonomic levels, such as the phylum level, when investigating the proportion of ciliates in a protist community (e.g. ). However, there is no detailed information on the assemblage of benthic ciliates in deep sea environments, except of the deep-sea hydrothermal vents  and cold seeps . In our study, we, for the first time, revealed that the pattern in the sediments of a seamount and the adjacent abyssal plain is quite different from that observed in intertidal zones and offshore areas, where spirotrichs are predominant in terms of relative abundance and species richness [20,21,22, 29]. By contrast, in the sediments of the seamount and the abyssal plain, Plagiopylea ciliates were obviously abundant and diverse. The high sequence contribution implies that they might have contributed a high biomass in the ciliate community in these habitats . Plagiopylea ciliates have generally been found in high sulfide, anoxic sediments . Previous studies indicate that these ciliates are abundant in sediments from the hydrothermal vents [22, 27]. Our study, therefore, extends the knowledge on the distribution of Plagiopylea ciliates, and reveals their potential importance in the benthic ciliate communities of seamounts and abyssal plains. However, some ciliate genera, like Pleuronema, which contains more than 20 morphotypes reported from various geographical locations of intertidal and coastal zones , were found to be widely distributed along the pronounced range of water depths, revealing their cosmopolitan character.
Regarding the proportion of rare OTUs in the sediments of a seamount and adjacent abyssal plain, there was no obvious shift along the water depth gradient. The proportion of rare OTUs detected in our study was obviously lower than what is known from intertidal and coastal sediments . Rare taxa can be considered as a seed bank of genetic resources, and are hypothesized to include ecologically redundant taxa that could increase in abundance following environmental perturbation to maintain continuous ecosystem functioning [33, 34]. Compared to coastal zones, deep-sea benthic habitats are relatively stable environments, explaining the observed lower proportion of rare taxa.
Vertical distribution of ciliates in the deep-sea sediments
In intertidal sediments, previous studies of benthic protists have shown that the abundance and species richness of protozoa were highest at the top 1 to 2 cm [17, 18]. Likewise, ciliate abundance and biomass in the 0-2 cm layers of coastal sediments accounted for 77% and 81%, respectively . Several species, however, were also found to live below 5 cm, even in 8 cm sediment depth [17, 18]. Previous studies showed that the respective proportions of ciliate abundance were about 23% in the 3–8-cm layers . Our study detected a high diversity of ciliate sequences in 20 cm sediment depth. When working with DNA, this is certainly no proof for the presence of living and active ciliate species, as DNA can be preserved in marine sediments over time. However, ciliates are very well adapted to a variety of different environment types, being able to not only tolerate anoxia, high salt or hydrogen sulfide concentrations or high pressure but also combinations of those [22, 27, 28]. It, therefore, would not be surprising if ciliates evolved to live in this water and sediment depths, especially in the presence of large bacterial populations in these deep sediment spheres [35, 36], which could serve as food source for heterotrophic ciliates.
The observation of OTUs affiliated to oligotrich ciliates in the deep layer of the deep-sea sediments confirmed a fraction of ciliate DNA sequences which likely do not belong to living organisms, but represent either extracellular DNA or cysts of planktonic species. The deep-sea floor, therefore, appears as a global DNA repository, which preserves molecular information of organisms living in the sediments, as well as in the overlying water column . In contrast, OTUs, which are found in RNA surveys, were more likely active species [38, 39]. Thus, RNA surveys should provide a better representation of in situ protist biomass and diversity. Clearly, in the future, the analyses of deep-sea RNA will be helpful to identify metabolically active organisms.
High novelty of ciliate genetic diversity
A high degree of genetic novelty was uncovered within the ciliates from the sediments of a seamount and the adjacent abyssal plain, although ciliates are considered as the best-known group within the protists. The application of high throughput DNA sequencing contributed substantially to the detection of a broader protist diversity [7, 40]. The extremely high molecular diversity largely exceeded the one detected by morphological methods. Thus, there might be some artifacts resulting from limitations of the techniques. Previous studies have indicated that protists are well known for having extremely high rDNA copy numbers [41, 42], which might result in a high polymorphism of SSU rDNA in a single cell and lead to the overestimation of molecular diversity in environmental samples. Gong et al.  showed that the minimal similarity between two copies in a single ciliate cell was 99.1%. The commonly used cutoff of 97% for creating OTUs based on the ciliate V4 fragment of their 18S rDNA is very well analyzed and established and excludes the effect of intragenomic variations [23, 43,44,45,46].
Moreover, there is still an obvious gap between the obtained sequences and the reference databases, which are used for taxonomic assignment of the sequences. The ciliate class Plagiopylea is such an example. Only one species, i.e. Epalxella antiquorum, within the order Odontostomatida and a total of eight plagiopylean species had their 18S rDNA sequences deposited in the NCBI database (searched on July 4th, 2017). The lack of reference sequences in the database resulted in all plagiopylean OTUs having a sequence similarity of less than 95% to the closest reference sequence. Additionally, even the most abundant plagiopylean OTU had only a sequence similarity of less than 90% to the reference sequences. More efforts in the isolation, cultivation and description of protists are necessary to link the environmental sequences to the real protist inventory [23, 47]. The design of novel species-specific primers and probes based on the retrieved sequences will also help to identify the target species by molecular techniques [48, 49]. In addition, different ciliate species were not equally distributed along a vertical sediment gradient and a large proportion of unknown ciliate diversity was found in deeper sediment layers. In the future, a more adequate sampling strategy is needed, which could maximize the number of different deep-sea habitats, and thus, detect more new ciliate species.
Our data point to obvious variations in ciliate communities along a sediment depth gradient in a seamount and its adjacent abyssal plain in the tropical Western Pacific Ocean and reflect the heterogeneity and diversity of benthic habitats. Water depth occurred as a less important factor impacting ciliate distribution in deep-sea sediments unlike observed for benthic ciliates in shallow seafloor habitats. Furthermore, our results indicate that ciliate diversity is not as well known as previously assumed and further efforts have to be made towards the identification of ciliate species. In this respect, deep-sea sediments appear as a hotspot for capturing and investigating novel ciliate species.
Study sites and sampling
Sediment samples were collected from four sites in the tropical Western Pacific Ocean (Fig. 5). Among these, site S-M is located at the middle of a seamount adjacent to the Yap Trench at a water depth of 813 m, and site S-B at the foot of the seamount at a water depth of 3812 m (Table 1). Sediment samples of S-M and S-B were collected using a push-corer operated by the remotely operated vehicle (ROV) “Discovery” onboard the R/V KEXUE in December 2014. Sampling sites DS1 and DS2 are located in the adjacent abyssal plain with water depths of 4042 m and 4566 m, respectively (Fig. 5). Sediment samples of DS1 and DS2 were taken using a 0.25-m2 modified Gray-O’Hara box corer in December 2014. Sites S-M and S-B are about 20 km distant from each other, and about 200 km apart from site DS2 and about 350 km apart from site DS1. Site DS1 is about 170 km distant from site DS2.
Four layers of the sediment cores (0-1 cm, 9-10 cm, 19-20 cm and 29-30 cm; labelled as 1, 2, 3 and 4, respectively) were sampled at site S-M, three layers (0-1 cm, 9-10 cm and 19-20 cm; labelled as 1, 2 and 3, respectively) at sites S-B and DS1, and two layers (0-1 cm, 9-10 cm; labelled as 1 and 2, respectively) at site DS2.
For each site and each layer, three replicate samples, each containing about 15 g sediments, were taken. Samples were stored at −80 °C for DNA extraction, at −20 °C for examination of total organic carbon (TOC), and at 4 °C for grain size measurements. TOC was determined in a Vario TOC cube (Elementar, Germany). Grain size analyses were performed using a Laser Diffraction Particle Size Analyzer (Cilas 940 L).
DNA extraction, PCR amplification and high-throughput sequencing
Environmental DNA was extracted from 0.3 g sediment of each replicate sample (a total of 3 DNA samples for each sediment sample) using the Power Soil DNA isolation kit (MoBio Laboratories, USA) according to the manufacturer’s protocol. Ciliate sequences were amplified by a nested PCR approach , using ciliate-specific 18S rRNA gene primers in the first reaction. The second PCR reaction employed a primer set specific for the hypervariable V4 region, which was included in the first PCR product. To minimize PCR errors, we used the PrimeSTAR GXL DNA High Fidelity Polymerase (TAKARA BIO INC., Japan). Three products from each sample were pooled.
Sequencing libraries were constructed using the NEB Next® Ultra™ DNA Library Prep Kit for Illumina (NEB, USA). Quality of the libraries was assessed with an Agilent Bioanalyzer 2100 system. Finally, libraries were sequenced on an Illumina MiSeq platform, generating 300-bp paired-end reads. The ciliate sequence reads have been deposited at the National Center for Biotechnology Information (NCBI) Sequence Read Archive under the accession number SRP101585.
Sequence data processing
Paired-end reads were merged and then filtered. Quality filtering of the raw sequences was performed according to the QIIME quality control process . Afterwards, a set of unique sequences was identified and the number of occurrences for each sequence was recorded by UPARSE software, followed by discarding all singleton sequences . A sequence similarity of 97% was used to delineate ciliate OTUs [43,44,45,46] by the UPARSE-OTU algorithm. Representative sequences from each OTU were extracted and subjected to the basic local alignment search tool (BLAST) analyses against the Silva database (v. 123) as implemented in the QIIME pipeline (v. 1.9.0) .
In each sample, OTUs were defined as rare when their number was equal to or less than 0.1% of all the sequences in the sample, and abundant, when their number was equal to or more than 1% of all the sequences in the sample [34, 53].
Statistical analyses were conducted in R using the vegan and fossil packages, unless stated otherwise. Rarefaction analyses were conducted in order to investigate the degree of sample saturation. Alpha diversity for each sample was calculated by the exponential of the Shannon index H′ (exp(H′), effective number of species) . After conversion, the effective numbers of species lets us avoid the serious misinterpretations spawned by the nonlinearity of the Shannon index H′ . Additionally, the ciliate OTU richness was determined as the number of OTUs in each sample. Spearman’s coefficient was used to relate ciliate alpha diversity and the environmental parameters, including water depth, TOC, the median grain size and the proportions of sand, silt and clay.
Partitioning of diversity (beta diversity) was investigated by an unweighted pair-group method with arithmetic means (UPGMA) cluster analysis based on the incidence Jaccard index. Stability in cluster analysis was evaluated by bootstrap resampling (100 times) using the clusterboot function in the fpc package of R . Here, the Jaccard similarities of the original clusters to the most similar clusters in the resampled data were computed. The mean over these similarities was then used as an index of the stability of a cluster. Generally, “highly stable” clusters should yield average Jaccard similarities of 0.85 and above. Clusters with a Jaccard similarity value smaller or equal to 0.5 were considered as a “dissolved cluster” . Prior to alpha and beta diversity analyses, the OTU table was randomly subsampled to the lowest number of sequences present in a sample (n = 25,169).
RALATE function in PRIMER v6 (Plymouth Marine Laboratory, UK) was applied to calculate Spearman’s correlation coefficients between the ciliate community composition and the environmental parameters. The BIOENV function in PRIMER v6 was used to calculate the correlation between the ciliate community and the environmental similarity matrix to determine the variables that best explain variation in the ciliate community.
The significance of the relationship between community Bray-Curtis dissimilarity and geographical distance (distance in kilometers between pairs of sites) within the detected distance-decay extent was assessed by a Mantel test for each data set. The Mantel tests were performed using the vegan package in R .
Identification of novel diversity
Identification of novel diversity followed the description of Filker et al. 2015 . Briefly, representative sequences of all detected OTUs were aligned with Seaview v.4.6.1 using “clusto” . Based on these alignments, pairwise similarities were calculated (custom script) and used for network construction in R (“igraph” package) . The resulting network was visualized and modified with GEPHI v.0.9.1. according to the OTU taxonomic affiliation and the Blast Hit values . In the network, two nodes were connected by an edge if they shared a sequence similarity of at least 90%.
Comparison of benthic ciliate communities from different marine habitats
We collected a set of publicly available data that were related to high throughput sequencing of the ciliate V4 fragment of 18S rDNA. These sequences were obtained from the surface sediments of the intertidal zone (SPR068269, ), the continental shelf of the Yellow Sea and East China Sea, and the hydrothermal vents in Okinawa Trough (SPR064020, ). Sequence data were processed as the protocol mentioned above. The significance of differences in the OTU richness among different habitats were estimated with Tukey’s HSD test. To investigate the partitioning of ciliate diversity, the abundance-based OTU table was transformed into a presence/absence table prior to calculating Sorensen-Dice indices. Sorensen-Dice values were then used for UPGMA cluster analyses. The ANOSIM function in PRIMER v6 was used to test the differences among different habitats. ANOSIM provides an R-statistic to evaluate the dissimilarity of groups, thus, groups are dissimilar if R-statistic is close to 1. Prior to alpha and beta diversity analyses, the OTU table was randomly subsampled to the lowest number of sequences present in a sample (n = 22,000).
Basic local alignment search tool
Median grain size
Operational taxonomic unit
Remotely operated vehicle
Total organic carbon
Unweighted pair-group method with arithmetic means
Danovaro R. Extending the approaches of biodiversity and ecosystem functioning to the deep ocean. In: Solan M, Aspden RJ, Paterson DM, editors. Marine biodiversity and ecosystem functioning: frameworks, methodologies, and integration. Oxford: Oxford University Press; 2012. p. 115–26.
Creer S, Sinniger F. Cosmopolitanism of microbial eukaryotes in the global deep seas. Mol Ecol. 2012;21:1033–5.
Pitcher TJ, Morato T, Hart PJB, Clark MR, Haggan N, Santos RS. Seamounts: ecology, fisheries, and conservation. Oxford: Blackwell; 2007.
Clark MR, Rowden AA, Schlacher T, Williams A, Consalvey M, Stocks KI, Rogers AD, O’Hara TD, White M, Shank TM, Hall-Spencer JM. The ecology of seamounts: structure, function, and human impacts. Annu Rev Mar Sci. 2010;2:253–78.
Ramirez-Llodra E, Brandt A, Danovaro R, De Mol B, Escobar E, German CR, Levin LA, Martinez Arbizu P, Menot L, Buhl-Mortensen P, Narayanaswamy BE, Smith CR, Tittensor DP, Tyler PA, Vanreusel A, Vecchione M. Deep, diverse and definitely different: unique attributes of the world’s largest ecosystem. Biogeosciences. 2010;7:2851–99.
Stefanoudis PV, Bett BJ, Gooday AJ. Abyssal hills: influence of topography on benthic foraminiferal assemblages. Prog Oceanogr. 2016;148:44–55.
Pawlowski J, Christen R, Lecroq B, Bachar D, Shahbazkia HR, Amaral-Zettler L, Guillou L. Eukaryotic richness in the abyss: insights from pyrotag sequencing. PLoS One. 2011;6:4.
Bik HM, Sung W, De Ley P, Baldwin JG, Sharma J, Rocha-Olivares A, Thomas WK. Metagenetic community analysis of microbial eukaryotes illuminates biogeographic patterns in deep-sea and shallow water sediments. Mol Ecol. 2012;21:1048–59.
Edgcomb VP, Kysela DT, Teske A, Gomez AD, Sogin ML. Benthic eukaryotic diversity in the Guaymas Basin hydrothermal vent environment. Proc Natl Acad Sci U S A. 2002;99:7658–62.
López-García P, Philippe H, Gail F, Moreira D. Autochthonous eukaryotic diversity in hydrothermal sediment and experimental microcolonizers at the mid-Atlantic ridge. Proc Natl Acad Sci U S A. 2003;100:697–702.
López-García P, Vereshchaka A, Moreira D. Eukaryotic diversity associated with carbonates and fluid-seawater interface in lost City hydrothermal field. Environ Microbiol. 2007;9:546–54.
Sauvadet AL, Gobet A, Guillou L. Comparative analysis between protist communities from the deep-sea pelagic ecosystem and specific deep hydrothermal habitats. Environ Microbiol. 2010;12:2946–64.
Urich T, Lanzen A, Stokke R, Pedersen RB, Bayer C, Thorseth IH, Schleper C, Steen IH, Ovreas L. Microbial community structure and functioning in marine sediments associated with diffuse hydrothermal venting assessed by integrated meta-omics. Environ Microbiol. 2014;16:2699–710.
Hausmann K, Hulsmann N, Polianski I, Schade S, Weitere M. Composition of benthic protozoan communities along a depth transect in the eastern Mediterranean Sea. Deep-Sea Res I. 2002;49:1959–70.
Pachiadaki MG, Taylor C, Oikonomou A, Yakimov MM, Stoeck T, Edgcomb V. In situ grazing experiments apply new technology to gain insights into deep-sea microbial food webs. Deep Sea Res II. 2016;129:223–31.
Gong J, Shi F, Ma B, Dong J, Pachiadaki M, Zhang X, Edgcomb VP. Depth shapes alpha- and beta-diversities of microbial eukaryotes in surficial sediments of coastal ecosystems. Environ Microbiol. 2015;17:3722–37.
Fenchel T. The ecology of marine microbenthos. IV structure and function of the benthic ecosystem, its chemical and physical factors and the microfauna communities with special reference to the ciliated protozoa. Ophelia. 1969;6:1–182.
Berninger UG, Epstein SS. Vertical distribution of benthic ciliates in response to the oxygen concentration in an intertidal North Sea sediment. Aquat Microb Ecol. 1995;9:229–36.
Carey PC. Marine interstitial ciliates: an illustrated key. New York: Chapman and Hall; 1992.
Hamels I, Muylaert K, Sabbe K, Vyverman W. Contrasting dynamics of ciliate communities in sandy and silty sediments of an estuarine intertidal flat. Eur J Protistol. 2005;41:241–50.
Zhao F, Xu K. Distribution of ciliates in intertidal sediments across geographic distances: a molecular view. Protist. 2017;168:172–81.
Zhao F, Xu K. Molecular diversity and distribution pattern of ciliates in sediments from deep-sea hydrothermal vents in the Okinawa trough and adjacent sea areas. Deep-Sea Res I. 2016;116:22–32.
Forster D, Dunthorn M, Mahe F, Dolan JR, Audic S, Bass D, Bittner L, Boutte C, Christen R, Claverie JM, Decelle J, Edvardsen B, Egge E, Eikrem W, Gobet A, Kooistra W, Logares R, Massana R, Montresor M, Not F, Ogata H, Pawlowski J, Pernice MC, Romac S, Shalchian-Tabrizi K, Simon N, Richards TA, Santini S, Sarno D, Siano R, Vaulot D, Wincker P, Zingone A, de Vargas C, Stoeck T. Benthic protists: the under-charted majority. FEMS Microbiol Ecol. 2016;92:8.
Meng Z, Xu K, Dai R, Lei Y. Ciliate community structure, diversity and trophic role in offshore sediments from the Yellow Sea. Eur J Protistol. 2012;48:73–84.
Scheckenbach F, Hausmann K, Wylezich C, Weitere M, Arndt H. Large-scale patterns in biodiversity of microbial eukaryotes from the abyssal sea floor. Proc Natl Acad Sci U S A. 2010;107:115–20.
Lejzerowicz F, Esling P, Pawlowski J. Patchiness of deep-sea benthic foraminifera across the Southern Ocean: insights from high-throughput DNA sequencing. Deep-Sea Res II. 2014;108:17–26.
Coyne KJ, Countway PD, Pilditch CA, Lee CK, Caron DA, Cary SC. Diversity and distributional patterns of ciliates in Guaymas Basin hydrothermal vent sediments. J Eukaryot Microbiol. 2013;60:433–47.
Takishita K, Kakizoe N, Yoshida T, Maruyama T. Molecular evidence that phylogenetically diverged ciliates are active in microbial mats of deep-sea cold-seep sediment. J Eukaryot Microbiol. 2010;57:76–86.
Du Y, Xu K, Warren A, Lei Y, Dai R. Benthic ciliate and meiofaunal communities in two contrasting habitats of an intertidal estuarine wetland. J Sea Res. 2012;70:50–63.
Fu R, Gong J. Single cell analysis linking ribosomal (r)DNA and rRNA copy numbers to cell size and growth rate provides insights into molecular protistan ecology. J Eukaryot Microbiol. 2017; 10.1111/jeu.12425.
Esteban G, Finlay BJ, Embley TM. New species double the diversity of anaerobic ciliates in a Spanish lake. FEMS Microbiol Lett. 1993;109:93–100.
Wang Y, Song W, Warren A, Al-Rasheid KAS, Al-Quraishy SA, Al-Farraj SA, Hu X, Pan H. Descriptions of two new marine scuticociliates, Pleuronema sinica n. Sp and P. wilberti n. Sp (Ciliophora: Scuticociliatida), from the Yellow Sea, China. Eur J Protistol. 2009;45:29–37.
Caron DA, Countway PD. Hypotheses on the role of the protistan rare biosphere in a changing world. Aquat Microb Ecol. 2009;57:227–38.
Pedrós-Alió C. The rare bacterial biosphere. Annu Rev Mar Sci. 2012;4:449–66.
Kallmeyer J, Pockalny R, Adhikari RR, Smith DC, D'Hondt S. Global distribution of microbial abundance and biomass in subseafloor sediment. Proc Natl Acad Sci U S A. 2012;109:16213–6.
Parkes RJ, Cragg BA, Wellsbury P. Recent studies on bacterial populations and processes in subseafloor sediments: a review. Hydrogeol J. 2000;8:11–28.
Coolen MJL, Saenz JP, Giosan L, Trowbridge NY, Dimitrov P, Dimitrov D, Eglinton TI. DNA and lipid molecular stratigraphic records of haptophyte succession in the Black Sea during the Holocene. Earth Planet Sci Lett. 2009;284:610–21.
Stoeck T, Zuendorf A, Breiner HW, Behnke A. A molecular approach to identify active microbes in environmental eukaryote clone libraries. Microb Ecol. 2007;53:328–39.
Jones SE, Lennon JT. Dormancy contributes to the maintenance of microbial diversity. Proc Natl Acad Sci U S A. 2010;107:5881–6.
Stoeck T, Bass D, Nebel M, Christen R, Jones MDM, 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.
Gong J, Dong J, Liu XH, Massana R. Extremely high copy numbers and polymorphisms of the rDNA operon estimated from single cell analysis of Oligotrich and Peritrich ciliates. Protist. 2013;164:369–79.
Weber AAT, Pawlowski J. Wide occurrence of SSU rDNA intragenomic polymorphism in foraminifera and its implications for molecular species identification. Protist. 2014;165:645–61.
Balzano S, Abs E, Leterme SC. Protist diversity along a salinity gradient in a coastal lagoon. Aquat Microb Ecol. 2015;74:263–77.
Massana R, Gobet A, Audic S, Bass D, Bittner L, Boutte C, Chambouvet A, Christen R, Claverie JM, Decelle J, Dolan JR, Dunthorn M, Edvardsen B, Forn I, Forster D, Guillou L, Jaillon O, Kooistra W, Logares R, Mahe F, Not F, Ogata H, Pawlowski J, Pernice MC, Probert I, Romac S, Richards T, Santini S, Shalchian-Tabrizi K, Siano R, Simon N, Stoeck T, Vaulot D, Zingone A, de Vargas C. Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing. Environ Microbiol. 2015;17:4035–49.
Nebel M, Pfabel C, Stock A, Dunthorn M, Stoeck T. Delimiting operational taxonomic units for assessing ciliate environmental diversity using small-subunit rRNA gene sequences. Environ Microbiol Reports. 2011;3:154–8.
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:185–7.
Filker S, Gimmler A, Dunthorn M, Mahe F, Stoeck T. Deep sequencing uncovers protistan plankton diversity in the Portuguese Ria Formosa solar saltern ponds. Extremophiles. 2015;19:283–95.
Orsi W, Edgcomb V, Faria J, Foissner W, Fowle WH, Hohmann T, Suarez P, Taylor C, Taylor GT, Vd'acny P, Epstein SS. Class Cariacotrichea, a novel ciliate taxon from the anoxic Cariaco Basin. Venezuela Int J Syst Evol Microbiol. 2012;62:1425–33.
Gimmler A, Stoeck T. Mining environmental high-throughput sequence data sets to identify divergent amplicon clusters for phylogenetic reconstruction and morphotype visualization. Environ Microbiol Reports. 2015;7:679–86.
Stock A, Edgcomb V, Orsi W, Filker S, Breiner HW, Yakimov MM, Stoeck T. Evidence for isolated evolution of deep-sea ciliate communities through geological separation and environmental selection. BMC Microbiol. 2013;13:150.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Tumbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.
Fuhrman JA. Microbial community structure and its functional implications. Nature. 2009;459:193–9.
Jost L. Partitioning diversity into independent alpha and beta components. Ecology. 2007;88:2427–39.
Hennig C. Cluster-wise assessment of cluster stability. Comput Stat Data An. 2007;52:258–71.
Hennig C. Dissolution point and isolation robustness: robustness criteria for general cluster analysis methods. J Multivar Anal. 2008;99:1154–76.
Oksanen J, Guillaume BF, Kindt R, Legendre P, O'Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H. Vegan: community ecology package. R package version 1.17-3. 2010. http://CRAN.R-project.org/package=vegan
Galtier N, Gouy M, Gautier C. SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996;12:543.
Csárdi G, Nepusz T. The igraph software package for complex network research. Inter J Complex Syst. 2006;1695:1–9.
Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. ICWSM. 2009;8:361–2.
We thank the WPOS sample center of Institute of Oceanology, Chinese Academy of Sciences and RV KEXUE for providing the sediments samples. Thanks are extended to Dr. Junlong Zhang for his help in making the figure of sampling sites.
This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (No. XDA11030201), the National Basic Research Program of China (973 Program) (No. 2015CB755902), the Scientific and Technological Innovation Project Financially Supported by Qingdao National Laboratory for Marine Science and Technology (No. 2016ASKJ14) and the China Scholarship Council (CSC, No. 201604910395).
Availability of data and materials
All relevant data are within the paper and its additional files. The raw sequence reads are deposited at the National Center for Biotechnology Information (NCBI) Sequence Read Archive under the accession number SRP101585.
Ethics approval and consent to participate
Our study did not involve the use of any animal or human data or tissue. No specific permission was required for this study according to the Regulations for the Administration of Affairs Concerning Experimental Animals approved by the State Council of the People’s Republic of China (CLI.2.4104).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rarefaction curves for all samples under study. (TIFF 3182 kb)
UPGMA clustering analysis based on the Bray-Curtis similarity coefficient (a) and the relative proportion of sequences (b) related to the major taxonomic groups of ciliates detected in the 12 sediment samples. (TIFF 713 kb)
Spearman’s rank correlation coefficients between ciliate alpha diversity in the surface layer sediments and the environmental parameters. P-value <0.05 are considered as significant. (DOCX 14 kb)
OTU-table displaying, for each sample, the number of sequences per OTU, the BLAST similarity of the representative sequence of each OTU to the reference sequence and the taxonomic assignment of the representative sequence. (XLSX 18 kb)