Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean

Evaluating sediment and water sampling methods for the estimation of deep-sea biodiversity using environmental DNA

Despite representing one of the largest biomes on earth, biodiversity of the deep seafloor is still poorly known. Environmental DNA metabarcoding offers prospects for fast inventories and surveys, yet requires standardized sampling approaches and careful choice of environmental substrate. Here, we aimed to optimize the genetic assessment of prokaryote (16S), protistan (18S V4), and metazoan (18S V1-V2, COI) communities, by evaluating sampling strategies for sediment and aboveground water, deployed simultaneously at one deep-sea site. For sediment, while size-class sorting through sieving had no significant effect on total detected alpha diversity and resolved similar taxonomic compositions at the phylum level for all markers studied, it effectively increased the detection of meiofauna phyla. For water, large volumes obtained from an in situ pump (~ 6000 L) detected significantly more metazoan diversity than 7.5 L collected in sampling boxes. However, the pump being limited by larger mesh sizes (> 20 µm), only captured a fraction of microbial diversity, while sampling boxes allowed access to the pico- and nanoplankton. More importantly, communities characterized by aboveground water samples significantly differed from those characterized by sediment, whatever volume used, and both sample types only shared between 3 and 8% of molecular units. Together, these results underline that sediment sieving may be recommended when targeting metazoans, and aboveground water does not represent an alternative to sediment sampling for inventories of benthic diversity.

Ciliate Diversity From Aquatic Environments in the Brazilian Atlantic Forest as Revealed by High-Throughput DNA Sequencing

Rainforest aquatic ecosystems include complex habitats with scarce information on their unicellular eukaryote diversity and community structure. We have investigated the diversity of ciliates in freshwater and brackish environments along the Brazilian Atlantic Forest, based on the hypervariable V4 region of the 18S-rDNA obtained by high-throughput DNA sequencing. Our analyses detected 409 ciliate taxonomic units (OTUs), mostly attributed to the classes Oligohymenophorea and Spirotrichea. A total of 11 classes, 12 subclasses, 112 genera, and 144 species were reported. We found the following: (a) the ciliate communities are more diverse in freshwater- than in Atlantic Forest-associated brackish environments; (b) the ciliate communities are composed by a small amount of highly abundant OTUs, but a high number of low-abundant or rare OTUs; (c) nearly one-third of the ciliate OTUs share less than 97% sequence identity to reference sequences and (d) phylogenetic inference supports the hypothesis that the V4 region of the Ciliophora 18S-rDNA is a suitable marker for accurate evolutionary inferences at class level. Our results showed that a considerable fraction of the HTS-detected diversity of ciliates from Brazilian Atlantic Forest is not represented in the currently available molecular databases.

Bathymetric gradient shapes the community composition rather than the species richness of deep-sea benthic ciliates

The bathymetric gradient is one of the most important factors that regulate the distribution of life. However, community variations of benthic ciliates along bathymetric gradients in the deep sea remain rather unexplored. In this study, we hypothesize that in the deep sea, the bathymetric gradient shapes the benthic ciliate community composition rather than the species richness. Here, we evaluated the distribution patterns and drivers of benthic ciliate communities of an abyssal plain, a seamount, and a trench with water depths ranging from 800 m down to 6600 m by high throughput eDNA sequencing and statistical analyses. We observed no significant correlation between ciliate operated taxonomic unit (OTU) richness and water depth. A meta-analysis, which combined our previously published data from the neritic habitats, supports the notion that water depth exceeding 800 m has little effect on the richness of benthic ciliate species. In contrast, the composition of deep-sea ciliate communities was significantly distinct in different habitats along the bathymetric gradients. A SourceTracker analysis revealed extremely low connectivity among ciliate communities along the bathymetric gradients. More than 95% of the community dissimilarity in the deep-sea floor was attributed to species replacement, which might be caused by environmental sorting or historical constraints. Furthermore, the observed community variations could be ascribed more to water depth than to geographic distance. The findings imply that the strong force of environmental sorting along the bathymetric gradients and the low connectivity among the ciliate communities might lead to an isolated evolution. This could shape the community composition rather than the species richness, which is mainly determined by the limited nutrient availability and the extreme environmental conditions in the deep sea.

Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora)

Metabarcoding and high-throughput sequencing methods have greatly improved our understanding of protist diversity. Although the V4 region of small subunit ribosomal DNA (SSU-V4 rDNA) is the most widely used marker in DNA metabarcoding of eukaryotic microorganisms, doubts have recently been raised about its suitability. Here, using the widely distributed ciliate genus Pseudokeronopsis as an example, we assessed the potential of SSU-V4 rDNA and four other nuclear and mitochondrial markers for species delimitation and phylogenetic reconstruction. Our studies revealed that SSU-V4 rDNA is too conservative to distinguish species, and a threshold of 97% and 99% sequence similarity detected only one and three OTUs, respectively, from seven species. On the basis of the comparative analysis of the present and previously published data, we proposed the multilocus marker including the nuclear 5.8S rDNA combining the internal transcribed spacer regions (ITS1-5.8S-ITS2) and the hypervariable D2 region of large subunit rDNA (LSU-D2) as an ideal barcode rather than the mitochondrial cytochrome c oxidase subunit 1 gene, and the ITS1-5.8S-ITS2 as a candidate metabarcoding marker for ciliates. Furthermore, the compensating base change and tree-based criteria of ITS2 and LSU-D2 were useful in complementing the DNA barcoding and metabarcoding methods by giving second structure and phylogenetic evidence.

Further analyses of variation of ribosome DNA copy number and polymorphism in ciliates provide insights relevant to studies of both molecular ecology and phylogeny

Sequence-based approaches, such as analyses of ribosome DNA (rDNA) clone libraries and high-throughput amplicon sequencing, have been used extensively to infer evolutionary relationships and elucidate the biodiversity in microbial communities. However, recent studies demonstrate both rDNA copy number variation and intra-individual (intra-genomic) sequence variation in many organisms, which challenges the application of the rDNA-based surveys. In ciliates, an ecologically important clade of microbial eukaryotes, rDNA copy number and sequence variation are rarely studied. In the present study, we estimate the intraindividual small subunit rDNA (SSU rDNA) copy number and sequence variation in a wide range of taxa covering nine classes and 18 orders of the phylum Ciliophora. Our studies reveal that: (i) intra-individual sequence variation of SSU rDNA is ubiquitous in all groups of ciliates detected and the polymorphic level varies among taxa; (ii) there is a most common version of SSU rDNA sequence in each cell that is highly predominant and may represent the germline micronuclear template; (iii) compared with the most common version, other variant sequences differ in only 1-3 nucleotides, likely generated during macronuclear (somatic) amplification; (iv) the intra-cell sequence variation is unlikely to impact phylogenetic analyses; (v) the rDNA copy number in ciliates is highly variable, ranging from 10³ to 10⁶, with the highest record in Stentor roeselii. Overall, these analyses indicate the need for careful consideration of SSU rDNA variation in analyses of the role of ciliates in ecosystems.