Next Article in Journal
Reactivity of a Recombinant Esterase from Thermus thermophilus HB27 in Aqueous and Organic Media
Previous Article in Journal
Comparative Genomics of Xylella fastidiosa Explores Candidate Host-Specificity Determinants and Expands the Known Repertoire of Mobile Genetic Elements and Immunity Systems
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 10
Issue 5
10.3390/microorganisms10050913
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity Patterns of Protists Are Highly Affected by Methods Disentangling Biological Variants: A Case Study in Oligotrich (s.l.) Ciliates

by
Jiahui Xu
1,†,
Jianlin Han
1,†,
Hua Su
1,
Changyu Zhu
1,2,
Zijing Quan
1,
Lei Wu
1 and
Zhenzhen Yi
1,*
1
Guangzhou Key Laboratory of Subtropical Biodiversity and Biomonitoring, Guangdong Provincial Key Laboratory for Healthy and Safe Aquaculture, School of Life Science, South China Normal University, Guangzhou 510631, China
2
Institute of Evolution & Marine Biodiversity, College of Fisheries, Ocean University of China, Qingdao 266003, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Microorganisms 2022, 10(5), 913; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10050913
Submission received: 28 February 2022 / Revised: 15 April 2022 / Accepted: 25 April 2022 / Published: 27 April 2022
(This article belongs to the Section Environmental Microbiology)
Browse Figures
Versions Notes

Abstract

:
Protists are a dominant group in marine microplankton communities and play important roles in energy flux and nutrient cycling in marine ecosystems. Environmental sequences produced by high-throughput sequencing (HTS) methods are increasingly used for inferring the diversity and distribution patterns of protists. However, studies testing whether methods disentangling biological variants affect the diversity and distribution patterns of protists using field samples are insufficient. Oligotrich (s.l.) ciliates are one group of the abundant and dominant planktonic protists in coastal waters and open oceans. Using oligotrich (s.l.) ciliates in field samples as an example, the present study indicates that DADA2 performs better than SWARM, UNOISE, UPARSE, and UCLUST for inferring diversity patterns of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions. UPARSE and UNOISE might underestimate species richness. SWARM might not be suitable for the resolution of alpha diversity owing to its rigorous clustering and sensitivity to sequence variations. UCLUST with 99% clustering threshold overestimates species richness, and the beta diversity pattern inferred by DADA2 is more reasonable than that of the other methods. Additionally, salinity is shown to be one of the key factors responsible for variations in the community distribution of ciliates, but infrequent marine–freshwater transitions occurred during evolutionary terms of this group.

1. Introduction

Protists, single-celled eukaryotes, are widely distributed in soil, marine environments, and freshwater worldwide and play key roles in energy flux as well as trophic interactions and nutrient cycling [1,2,3]. They are a dominant group in marine microzooplankton communities and act as primary producers, predators, decomposers, and/or parasites in marine ecosystems [2,4,5]. However, protists with low abundance in a region are easily ignored in observation in vivo due to their small size and difficulty in identification and cultivation [6,7]. Foissner [8] assumed that more than a half of the diversity of many protist groups has not been revealed. Consequently, insufficient sampling of protists in studies using isolated cells might halt revealing their accurate patterns of diversity, phylogeny, and transition. Fortunately, rare and cryptic protists could be revealed by high-throughput sequencing (HTS) with large data output, which provides us a chance to investigate more accurate diversity and distribution patterns of protists [9,10,11,12,13].
Considering that both real biological variants (inter-/intra-specific variants) and spurious sequences introduced by sequencing are included in HTS datasets, a method that disentangles biological variation is one of the main factors affecting the accuracy of downstream analyses [14,15,16,17]. Nowadays, two sequence grouping approaches are popular for disentangling biological variations [14]. One is the construction of operational taxonomic units (OTUs) with different quality filtering and clustering algorithms, e.g., UPARSE [16], UNOISE [18], UCLUST [19], and SWARM [20]; different clustering thresholds of sequences have been reported to define OTUs of protists (e.g., 95% in [21]; 97% in [22]; 98% in [23]; 99% in [24]). Another is correcting amplicon errors by generating Amplicon Sequence Variants (ASVs) at single nucleotide resolution, e.g., DADA2 [25]. To our knowledge, many references focusing on protist diversity based on HTS data are published every year, but there are only a few studies that compare the effects of different HTS data processing methods on analyses [17,26,27]. Some studies suggested that DADA2-derived ASVs could more accurately reproduce a known alpha diversity than SWARM-derived OTUs for ciliate species [27], and USEARCH and QIIME strongly affected the number of predicted OTUs but not the biogeographical patterns of protists [26]. Additionally, the numbers of OTUs varied with the change in clustering thresholds in prokaryotes and protists [28]. Consequently, protist diversity inferred from HTS data is highly affected not only by sequence grouping approaches but also by clustering thresholds producing OTUs. However, some questions are still unclear. For instance, testing whether ASVs and/or OTUs affect beta diversity patterns of protists in real case scenarios using field samples [27].
Oligotrich (s.l.) ciliates are a major group of the abundant and dominant planktonic protist communities in coastal waters and open oceans [1,29,30,31,32,33]. Salinity gradients of the Pearl River Estuary (PRE) and surrounding coasts range from 0.1‰ to 32.0‰ [34]. In the present study, we will explore whether methods disentangling biological variants will affect the diversity patterns of protists by using oligotrich (s.l.) ciliates collected from this region as an example. Our aim is to provide suggestions for data analyses in future studies. Moreover, we investigate the ecological transitions of oligotrich (s.l.) ciliates using both released sequences of identified species in GenBank and environmental sequences from this region since expanded taxa sampling has been proven to be beneficial for classifying transition patterns of ciliates [35].

2. Materials and Methods

2.1. Environmental Sequences of Oligotrich (s.l.) Ciliates

Oligotrich (s.l.) ciliates consist of two subclasses, Oligotrichia and Choreotrichia [36]. The environmental SSU rDNA V4 region sequences were downloaded from GenBank under the accession number PRJNA646537. The geographic locations of sampling sites and the salinity of each site in this study are shown in Figure 1. In summary, PRE1–PRE11 are sampling sites along the flow direction of the Pearl River Estuary (PRE), and sampling sites DY1–DY5, GZ1–GZ5, SZ1–SZ6, ZH1–ZH5, and ZJ1–ZJ6 are in the sampling areas Daya Bay (DY), Guangzhou (GZ), Shenzhen (SZ), Zhuhai (ZH), and Zhanjiang (ZJ), respectively. Among these sampling areas, SZ and DY are located on the eastern side of PRE, ZH and ZJ are on the western side of PRE, and GZ is in the upper reaches of PRE. The hydrological and nutritional conditions from the western and eastern sides of PRE are different, and PRE provides a complex ecosystem with various salinity gradients from the upper to lower reaches [34]. According to the Venice Salinity Classification System (Source: Limnology and Oceanography, 1958), the water salinity is classified as freshwater (≤0.5‰, low-brackish (0.5–18.0‰), high-brackish (18.0–30.0‰), and marine (30.0–40.0‰).
To test whether the methods of disentangling biological variants affect diversity patterns, all downloaded data were processed using five frequently used methods: UPARSE and UNOISE in USEARCH v11.0 [16,18], UCLUST in QIIME v1.9.1 [15,19], SWARM v3 [20], and DADA2 in QIIME2 v2019.10 [14,25]. In total, six oligotrich (s.l.) datasets, viz., UPARSE–97, UNOISE–100, UCLUST–97, UCLUST–99, SWARM–100, and DADA2–100, were produced from environmental sequences as follows. For datasets UPARSE–97, UNOISE–100, SWARM–100, UCLUST–97, and UCLUST–99, which produce OTUs, the same quality filtering protocol was performed in order to only compare differences among clustering algorithms. Briefly, the paired-end reads were merged with FLASH [37] and filtered with the following settings: sequences of length <200 or >500, average quality <20, ambiguous bases >0, or homopolymer length >6 were removed using QIIME v1.8.0 [15]. Dereplication and the discarding of singletons were performed by using USEARCH v11.0 [16]. Sequences of UPARSE–97 were clustered into OTUs at the default 97% similarity. Sequences of UNOISE–100 were denoised into zero-radius operational taxonomic units (ZOTUs), which are valid OTUs superior to conventional 97% OTUs, using unoise3 with default parameters [18]. Sequences of SWARM–100 were clustered using SWARM v3 with d = 1 [20] and subsequently subjected to chimera detection using VSEARCH v2.15.0 [38]. And 97% and 99% similarity thresholds were set for the UCLUST algorithm as an example to check whether diversity patterns varied with changes in clustering thresholds, considering that both of them have been reported to define OTUs of protists [22,24]. For DADA2–100, the demultiplexed reads with QIIME2 v2019.10 [25] were filtered and denoised using DADA2 [14] with “––p–trunc–len–f” and “––p–trunc–len–r” parameters set to 240 and 210, respectively. The sequences of each ASV with one nucleotide difference were produced. For all datasets, the Silva 138 SSURef database (https://0-doi-org.brum.beds.ac.uk/10.5281/ZENODO.3891931, accessed on 17 May 2021) was utilized to annotate taxonomic information for each representative sequence. After that, only sequences annotated as oligotrich (s.l.) ciliates (similarity > 90% as in [26]) were kept in each dataset for downstream analyses. We assumed that these sequences might reflect the abundance and community diversity of oligotrich (s.l.) ciliates even though these sequences are not biologically realistic. DY3, DY4, PRE11, SZ5, SZ6, and ZJ5 were deleted from subsequent analyses because fewer than five oligotrich (s.l.) sequences were detected in these sampling sites by all methods.
In order to explore the diversity patterns of oligotrich (s.l.) ciliates among different datasets, a Bray–Curtis distance matrix was used to quantify the community dissimilarity. Hierarchical clustering and Principal Co-ordinates Analysis (PCoA) were performed in the R platform [39] using the packages “ape” [40] and “vegan” [41]. The sequences of OTUs/ASVs from different analytical methods were log (x + 1) transformed to improve normality and homoscedasticity before downstream analysis. Hierarchical clustering was performed using the “ward.D2” algorithm of the “hclust” function. PCoA was performed based on the Bray–Curtis dissimilarity, and the “stat_ellipse” (level = 0.95) command was used to add ellipses by group except DY (too few points to calculate an ellipse) to visualization results from PCoA. The differences between predefined groups based on sampling areas were statistically tested by permutational multivariate analysis of variance (ADONIS) [42] using 1000 permutations.

2.2. Sequence Alignment and Phylogenetic Analyses

All available SSU rDNA sequences (>1000 bp, 191 in total) of identified oligotrich (s.l.) species as of August 2021, as well as six sequences of the subclass Hypotrichia used as the outgroup, were downloaded from GenBank (accession numbers in Supplementary Table S1). The dataset DADA2–100 performed better for inferring diversity patterns of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions than other five datasets due to more reliable alpha diversity and beta diversity (detailed information in discussion). Hence, the dataset containing 197 identified sequences listed in Supplementary Table S1 and 103 representative sequences from dataset DADA2–100 were aligned using the GUIDANCE2 server [43]. The resulting alignment was manually checked in SeaView v4 [44] for trimming two ends and ambiguous sites.
Maximum likelihood (ML) analyses were carried out using RAxML-HPC2 on XSEDE on CIPRES Science Gateway (http://www.phylo.org/sub_sections/portal, accessed on 30 August 2021) [45] using the GTRGAMMA model, and support for the best-scoring ML tree was assessed by 1000 bootstrap replicates. Bayesian inference (BI) analysis was also performed on CIPRES Science Gateway using MrBayes on XSEDE v3.2.7a, using the GTR+I+G model, which was selected by jModeltest v2.1.10 [46]. Markov chain Monte Carlo simulations were run for 10,000,000 generations with four chains. Trees were sampled every 100 generations, and the first 25% of trees were discarded as burn-in. The 50% majority rule consensus tree was used to calculate the posterior probabilities (PP) for each node. Trees were viewed and edited with FigTree v. 1.4.4 [47]. Finally, Mesquite v3.6 [48] was used to infer the most parsimonious pattern of marine–freshwater transitions of oligotrich (s.l.) ciliates using the ML tree above.

3. Results

3.1. Diversity Patterns of Oligotrich (s.l.) Ciliates

Since the total number of oligotrich (s.l.) OTUs/ASVs contained in each dataset varied greatly from 63 (UPARSE–97) to 248 (UCLUST–99) (Table 1), the ratio of oligotrich (s.l.) OTUs/ASVs number in each sampling area to the total number of oligotrich (s.l.) OTUs/ASVs within each dataset was also compared (Table 2). The proportions of oligotrich (s.l.) OTUs/ASVs in each sampling area/site in UPARSE–97 and UNOISE–100 were generally higher than those in the other four datasets (Table 2 and Supplementary Table S2). Among the six sampling areas, the largest proportion of oligotrich (s.l.) OTUs/ASVs was detected in PRE, with largest in UPARSE–97 (71.43%) and lowest in DADA2–100 (45.63%). The lowest proportion of oligotrich (s.l.) OTUs/ASVs was detected in DY based on DADA2–100 (5.83%), UCLUST–97 (13.46%), and UCLUST–99 (6.85%) (Table 2). However, the lowest one was detected in SZ based on SWARM–100 (10.95%), UNOISE–100 (13.48%), and UPARSE–97 (15.87%) (Table 2).
The hierarchical clustering analyses show that oligotrich (s.l.) samples grouped into two major clades in DADA2–100 (Figure 2a). Generally, Clade A consisted of most samples from sampling areas ZH and GZ, as well as samples of PRE1–PRE4. Clade B contained most samples from sampling areas ZJ, SZ, and DY, as well as samples of PRE5–PRE10. Different from the general clustering pattern in DADA2–100 (Figure 2a), samples ZJ1 and DY5 were far from Clade A and Clade B in SWARM–100, UPARSE–97, and UCLUST–97 (Figure 2b,d,e), and samples ZJ1, DY2, and DY5 were far from Clade A and Clade B in UNOISE–100 and UCLUST–99 (Figure 2c,f). Additionally, sample PRE5 fell into Clade A instead of Clade B in UPARSE–97 (Figure 2d). Notably, community structures of ZJ samples were more similar to each other than to PRE samples in DADA2–100 (Figure 2a). By contrast, some ZJ samples clustered with PRE samples first and then with other ZJ samples in all other five datasets (Figure 2b–f).
PCoA results showed that oligotrich (s.l.) samples were basically grouped depending on the six sampling areas (DY, GZ, PRE, SZ, ZH, ZJ), and PRE was divided into two groups (PRE1–PRE4, PRE5–PRE10) (Figure 3). ADONIS results supported the significant (p < 0.001) differentiations of community structures among these seven groups based on all six datasets (Figure 3). Consistent with the hierarchical clustering results (Figure 2), samples of PRE1–PRE4 were clearly separated from those of PRE5–PRE10, ZJ, SZ, and DY in all six datasets (Figure 3). However, for DADA2–100, samples of PRE1–PRE4, ZH, and GZ tended to group together and were separated from other sampling sites (Figure 3a).

3.2. Phylogeny and Transition Patterns of Oligotrich (s.l.) Ciliates

To investigate the ecological transition patterns of oligotrich (s.l.) ciliates, a phylogenetic tree including 197 sequences from identified species (Supplementary Table S1) and 103 ASVs from dataset DADA2–100 was constructed (Figure 4), as DADA2–100 yielded more reliable diversity patterns of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions than the other five datasets. Monophyly of the subclass Choreotrichia was supported by high support (98% ML, 1.00 BI), but that of subclass Oligotrichia was not. In the phylogenetic trees, crown clades usually contained identified oligotrich (s.l.) species and ASVs from various habitats. Though most ASV groups had identified species, some ASVs formed isolated clades (Clade 1, Clade 2 in Figure 4) without identified species. This indicates that oligotrich (s.l.) diversity in low-brackish habitats might have been underestimated in previous morphological investigations. Transition pattern analyses showed a high-brackish ancestor for oligotrich (s.l.) ciliates (Figure 4). The subclass Choreotrichia appeared to have evolved from the high-brackish ancestor with high support values (98% ML, 1.00 BI), with some transitions to freshwater, low-brackish, and marine areas and even back to high-brackish for some species (e.g., Tintinnopsis radix and Leprotintinnus nordqvisti) (Figure 4). By contrast, the subclass Oligotrichia seemed to be derived from an equivocal ancestor with poor support values (ML < 50%, BI < 0.50), with some transition to freshwater, low brackish, high-brackish, and marine areas and possible additional transition back to low-brackish areas for some species (e.g., Parallelostrombidium paraellipticum and P. obesum) (Figure 4).

4. Discussion

4.1. Methods Disentangling Biological Variants Highly Affect Diversity Patterns of Oligotrich (s.l.) Ciliates

The alpha and beta diversities of six oligotrich (s.l.) datasets were compared in order to check which method performed best for inferring diversity patterns of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions. The detailed discussion is as following:
Among six oligotrich (s.l.) datasets generated by different algorithms and thresholds (viz. DADA2–100, SWARM–100, UNOISE–100, UPARSE–97, UCLUST–97, UCLUST–99), the actual number of oligotrich (s.l.) OTU/ASV as well as ratio of the oligotrich (s.l.) OTUs/ASVs number in each sample area/site were different (Table 1 and Table 2, Supplementary Table S2). This is consistent with previous investigations that the alpha diversity of protists is highly variable depending on sequence grouping approaches, as well as software and clustering thresholds producing OTUs [26,27,28,49]. Recent investigations revealed that compared with OTUs, ASVs could more accurately reproduce a known alpha diversity from mock communities of various groups, e.g., bacteria [14,50,51], fungi [50], ciliated protists [27]. In this study, the number of oligotrich (s.l.) OTUs/ASVs inferred by UCLUST–99 (248) is unreasonable (Table 1), because only 288 morphological oligotrich (s.l.) ciliate species were reported in the South China Sea from 1991 to 2018 [52]. The numbers of oligotrich (s.l.) OTUs/ASVs inferred by UPARSE–97 (63) and UNOISE–100 (89) were much lower than those in the other three datasets (103–137) (Table 1). One possible explanation is that UPARSE and UNOISE could not detect fine-scale or low-abundance biological variations [14,18]. Previous studies showed that SWARM was not suitable for the resolution of genetic diversity and alpha diversity in the samples with high intraspecific sequence variations due to its rigorous clustering and sensitivity to sequence variations [20,27,53]. Only alpha diversity patterns revealed by DADA2–100 and UCLUST–97 seemed be reliable. Interestingly, numbers of ASVs/OTUs are comparable between DADA2–100 (103) and UCLUST–97 (104), but the ratio of oligotrich (s.l.) OTU/ASV numbers in each sampling area to the total number of oligotrich (s.l.) OTUs/ASVs within each dataset was rather different between these two datasets (Table 1 and Table 2). This indicates that beta diversity patterns should be compared between these two datasets in order to check which one is better.
In this study, a general beta diversity pattern of oligotrich (s.l.) ciliates was revealed in the six oligotrich (s.l.) datasets, but methods disentangling biological variants also had an impact on the beta diversity patterns (Figure 2 and Figure 3). Clearly, community variations were observed between PRE1–PRE4 and PRE5–PRE10 in DADA2–100, SWARM–100, UNOISE–100, UCLUST–97 and UCLUST–99, and PRE5 groups with PRE1–PRE4 instead of PRE6–PRE10 in UPARSE–97. This might be explained by the sharp increase in salinity between sampling sites PRE1–PRE4 (0.3–0.9‰) and PRE5–PRE10 (4.7–12.2‰). It is possible that the mixture of freshwater and seawater formed a low-salinity front between sampling sites PRE4 and PRE5 [54]. Numerous studies have proven that salinity appeared to be the factor that correlated best with distributions of phytoplankton and bacterioplankton in estuaries (e.g., [55,56,57,58,59]) and hence community compositions of oligotrich (s.l.) ciliates with phytoplankton and bacterioplankton as food also changed greatly between PRE1–PRE4 and PRE5–PRE10. Similar to PRE1–PRE4, sampling sites GZ1–GZ5 are in the upper estuary of the Pearl River Estuary and are highly influenced by large freshwater discharge from the Pearl River. Theoretically, community structures of oligotrich (s.l.) ciliates in the upper Pearl River Estuary with lower salinity and higher nutrient content should be much different from those in the lower Pearl River Estuary. In the present study, the GZ samples were grouped with PRE1–PRE4 in all six datasets (Figure 2 and Figure 3). Hierarchical clustering analyses of DADA2–100, SWARM–100, UPARSE–97, UCLUST–97, and UCLUST–99 showed that most ZH samples grouped with PRE1–PRE4 in Clade A, and most SZ samples fell into Clade B including PRE5–PRE10 (Figure 2a,b,d,f), although sampling areas ZH and SZ are located at a similar latitude of PRE5. Different community structures of oligotrich (s.l.) ciliates in ZH and SZ samples might be due to following reason. The surface flow velocity on the western side (ZH) is usually greater than that on the eastern side (SZ) of the Pearl River Estuary [60]. This indicates that community structures of oligotrich (s.l.) ciliates in ZH might be more highly influenced by river runoff than in SZ. Thus, cluster patterns of most SZ and ZH samples falling into Clade A in UNOISE–100 (Figure 2c) were less reliable than in the other five datasets (Figure 2a,b,d,f). Sampling areas ZJ and DY are located in the nearshore area of the South China Sea (surrounding regions of the Pearl River Estuary), and most samples from these sampling areas clustered with PRE5–PRE10 in all six datasets (Figure 2 and Figure 3). Possibly, salinity also plays an important role in this cluster pattern, since salinity values in sampling areas ZJ, DY and sampling sites PRE5–PRE10 were generally higher than those in other areas (Figure 1). Notably, both the geographical locations and salinity of sampling sites ZJ2–ZJ4 and ZJ6 were very similar to each other (Figure 1), but their close relationships based on hierarchical clustering analysis were only as obvious in DADA2–100 (Figure 2a). Additionally, three samples collected from DY area fell into a subclade of Clade B in DADA2–100, while DY5 always formed a separate basal clade in the other five datasets (Figure 2b–f). Hence, DADA2–100 showed the most reliable clustering patterns for ZJ and DY samples, although a general beta diversity pattern of oligotrich (s.l.) ciliates was revealed in the six oligotrich (s.l.) datasets. All these results reveal that the beta diversity pattern inferred by DADA2–100 was more reasonable in real case scenarios using field samples of oligotrich (s.l.) ciliates.
In sum, among the six datasets compared in this study, DADA2–100 performed best for inferring the diversity pattern of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions. As described in the Introduction, this is the first study to test whether ASVs and/or OTUs affect the beta diversity patterns of protists in real case scenarios. In future, more studies using filed samples are expected to test whether ASVs perform best for inferring diversity pattern of protists in various regions.

4.2. Community Distribution and Ecological Transitions of Oligotrich (s.l.) Ciliates in Environments with Various Salinity Gradients

As revealed in various groups of archaea, bacteria, and protists (e.g., [34,55,57,58,61]), the present study also revealed that salinity gradients play a vital role in structuring patterns for the community distribution of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions that encompass an entire freshwater–marine salinity gradient. First, communities of oligotrich (s.l.) ciliates in sampling sites PRE1–PRE10 along the flow direction of the Pearl River Estuary were divided into two distinct groups PRE1–PRE4 (Clade A) and PRE5–PRE10 (Clade B) (Figure 2a), which is consistent with the sharp change in salinity gradients between sampling sites PRE1–PRE4 (0.3–0.9‰) and PRE5–PRE10 (4.7–12.2‰). Additionally, all samples from high-brackish and marine habitats fell into Clade B, though these samples were collected from two sampling areas (ZJ and DY) of around 500 km. It is believed that both the physicochemical barrier of salinity gradients and the presence of locally adapted taxa limit the colonization success of microbes in different habitats with salinity gradients [62]. Within Clade B, several samples (DY2, ZJ1, PRE5–PRE10, SZ2–SZ4) were from freshwater and low-brackish habitats, which indicates that in addition to salinity, other environmental factors and geographical distance also shape the community distribution of oligotrich (s.l.) ciliates. This was also revealed in various prokaryotic and eukaryotic microbial groups in estuaries (e.g., [34,52,63]).
As mentioned above, salinity gradients have been proven to be important physicochemical factor in structuring community distribution and limiting transitions of microbes including protists. Ancestors of different ciliate groups seem to originate in various habitats. Oligotrich (s.l.) and hypotrich (present study, [35]) ciliates are ancestrally high-brackish, and peritrich and colpodean ciliates have freshwater/terrestrial ancestors [64,65]. However, species in various habitats are reported for each of these ciliate groups. That is because infrequent marine–freshwater transitions always occurred during their evolutionary terms (present study; [35,64,65,66]) due to the ability of microbes to rapidly adapt and highly colonize to new environments [62,67]. Environmental SSU rDNA V4 region sequences of the Pearl River Estuary and surrounding coasts cover oligotrich (s.l.) ciliates in various salinity gradients (0.1–32.0‰) [34], providing us a good chance to classify ecological transition patterns of this taxonomy group with expanded taxa sampling [35]. Although the short fragments of SSU rDNA have limited phylogenetic signals [64], previous studies have proven that the addition of SSU rDNA V4 region sequences produced from amplicons could improve the ecological transition patterns of ciliates by broad taxa sampling in various habits [35,64]. In our phylogenetic trees (Figure 4), some ASVs formed separate clades (Clade 1, Clade2) and might be new taxa. This is consistent with a previous report that the oligotrich (s.l.) morphospecies diversity was underestimated [68]. Additionally, these ASVs representing new taxa were from low-brackish habitats (Figure 4). By contrast, among all SSU rDNA sequences of identified oligotrich (s.l.) species deposited in GenBank, only 4% (8 out of 191) and 8% (15 out of 191) were from freshwater and low-brackish habitats (Supplementary Table S1). In future, more detailed transitions and evolutionary patterns of oligotrich (s.l.) ciliates are expected be outlined with broader sampling, especially freshwater and low-brackish habitats.

5. Conclusions

This study investigates whether methods disentangling biological variants will affect the diversity patterns of protists by using oligotrich (s.l.) ciliates in field samples as a case. Our work demonstrates that DADA2 performed better than SWARM, UNOISE, UPARSE, and UCLUST for inferring diversity patterns of oligotrich (s.l.) ciliates. In addition, salinity was shown to be one of the key factors responsible for variations in the community distribution of ciliates but infrequent marine–freshwater transitions occurred during the evolutionary terms of this group.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms10050913/s1, Table S1: Detailed information of identified oligotrich (s.l.) ciliate species from GenBank. All sequences are SSU rDNA sequences. Ecological habitat is classified as F = freshwater (≤0.5‰), L = low-brackish (0.5–18‰), H =high-brackish (18–30‰), and M = marine (30–40‰). References [66,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152] are cited in the Supplementary Materials. Table S2: Number and proportion of oligotrich (s.l.) OTUs/ASVs in each sampling site from the six datasets.

Author Contributions

Conceptualization, Z.Y.; formal analysis, J.X. and J.H.; data curation, J.X. and J.H.; writing—original draft preparation, J.X., J.H., H.S. and Z.Y.; writing—review and editing, C.Z., Z.Q. and L.W.; visualization, J.X., J.H., H.S. and C.Z.; funding acquisition, Z.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (grant number 32070406), and the Science and Technology Project of Guangzhou (No. 202102080168).

Data Availability Statement

The data generated and/or analyzed in this study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank Yusen Xu, former undergraduate of South China Normal University, and Weiwei Liu from South China Sea Institute of Oceanology, Chinese Academy of Sciences, for their help collecting samples.

Conflicts of Interest

The authors declare they had no relevant financial or non-financial interests to disclose.

References

  1. Pierce, R.W.; Turner, J.T. Ecology of planktonic ciliates in marine food webs. Rev. Aquat. Sci. 1992, 6, 139–181. [Google Scholar]
  2. Fenchel, T. The microbial loop-25 years later. J. Exp. Mar. Biol. Ecol. 2008, 366, 99–103. [Google Scholar] [CrossRef]
  3. Worden, A.Z.; Follows, M.J.; Giovannoni, S.J.; Wilken, S.; Zimmerman, A.E.; Keeling, P.J. Rethinking the marine carbon cycle: Factoring in the multifarious lifestyles of microbes. Science 2015, 347, 1257594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. De Vargas, C.; Audic, S.; Henry, N.; Decelle, J.; Mahé, F.; Logares, R.; Lara, E.; Berney, C.; Le Bescot, N.; Probert, I. Eukaryotic plankton diversity in the sunlit ocean. Science 2015, 348, 1261605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lynn, D. The Ciliated Protozoa: Characterization, Classification, and Guide to the Literature, 3rd ed.; Springer: Dordecht, The Netherlands, 2008. [Google Scholar]
  6. Amann, R.I.; Ludwig, W.; Schleifer, K.H. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol. Rev. 1995, 59, 143–169. [Google Scholar] [CrossRef]
  7. Eland, L.E.; Davenport, R.; Mota, C.R. Evaluation of DNA extraction methods for freshwater eukaryotic microalgae. Water Res. 2012, 46, 5355–5364. [Google Scholar] [CrossRef]
  8. Foissner, W. Biogeography and dispersal of micro-organisms: A review emphasizing protists. Acta Protozool. 2006, 45, 111–136. [Google Scholar] [CrossRef]
  9. Doherty, M.; Tamura, M.; Vriezen, J.A.; McManus, G.B.; Katz, L.A. Diversity of Oligotrichia and Choreotrichia ciliates in coastal marine sediments and in overlying plankton. Appl. Environ. Microbiol. 2010, 76, 3924–3935. [Google Scholar] [CrossRef] [Green Version]
  10. Li, G.; Su, L.; Zhang, Q.; Zhang, X.; Gong, J. Molecular diversity and biogeography of benthic ciliates in the Bohai Sea and Yellow Sea. Acta Oceanol. Sin. 2019, 38, 78–86. [Google Scholar] [CrossRef]
  11. Massana, R.; Gobet, A.; Audic, S.; Bass, D.; Bittner, L.; Boutte, C.; Chambouvet, A.; Christen, R.; Claverie, J.M.; Decelle, J. Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing. Environ. Microbiol. 2015, 17, 4035–4049. [Google Scholar] [CrossRef] [Green Version]
  12. Zhu, C.; Bass, D.; Wang, Y.; Shen, Z.; Song, W.; Yi, Z. Environmental parameters and substrate type drive microeukaryotic community structure during short-term experimental colonization in subtropical eutrophic freshwaters. Front. Microbiol. 2020, 11, 2252. [Google Scholar] [CrossRef] [PubMed]
  13. Zinger, L.; Boetius, A.; Ramette, A. Bacterial taxa–area and distance–decay relationships in marine environments. Mol. Ecol. 2014, 23, 954–964. [Google Scholar] [CrossRef] [PubMed]
  14. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  17. Forster, D.; Dunthorn, M.; Stoeck, T.; Mahé, F. Comparison of three clustering approaches for detecting novel environmental microbial diversity. PeerJ 2016, 4, e1692. [Google Scholar] [CrossRef] [Green Version]
  18. Edgar, R.C. UNOISE2: Improved error-correction for Illumina 16S and ITS amplicon sequencing. BioRxiv 2016, 081257. [Google Scholar] [CrossRef] [Green Version]
  19. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  20. Mahé, F.; Czech, L.; Stamatakis, A.; Quince, C.; de Vargas, C.; Dunthorn, M.; Rognes, T. Swarm v3: Towards tera-scale amplicon clustering. Bioinformatics 2022, 38, 267–269. [Google Scholar] [CrossRef]
  21. Countway, P.D.; Gast, R.J.; Dennett, M.R.; Savai, P.; Rose, J.M.; Caron, D.A. Distinct protistan assemblages characterize the euphotic zone and deep sea (2500 m) of the western North Atlantic (Sargasso Sea and Gulf Stream). Environ. Microbiol. 2007, 9, 1219–1232. [Google Scholar] [CrossRef]
  22. Wu, W.; Liu, H. Disentangling protist communities identified from DNA and RNA surveys in the Pearl River–South China Sea Continuum during the wet and dry seasons. Mol. Ecol. 2018, 27, 4627–4640. [Google Scholar] [CrossRef] [PubMed]
  23. 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. Rep. 2011, 3, 154–158. [Google Scholar] [CrossRef] [PubMed]
  24. Doherty, M.; Costas, B.A.; McManus, G.B.; Katz, L.A. Culture-independent assessment of planktonic ciliate diversity in coastal northwest Atlantic waters. Aquat. Microb. Ecol. 2007, 48, 141–154. [Google Scholar] [CrossRef] [Green Version]
  25. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef] [PubMed]
  26. Boscaro, V.; Rossi, A.; Vannini, C.; Verni, F.; Fokin, S.I.; Petroni, G. Strengths and biases of high-throughput sequencing data in the characterization of freshwater ciliate microbiomes. Microb. Ecol. 2017, 73, 865–875. [Google Scholar] [CrossRef]
  27. Forster, D.; Lentendu, G.; Filker, S.; Dubois, E.; Wilding, T.A.; Stoeck, T. Improving eDNA-based protist diversity assessments using networks of amplicon sequence variants. Environ. Microbiol. 2019, 21, 4109–4124. [Google Scholar] [CrossRef]
  28. Logares, R.; Mangot, J.F.; Massana, R. Rarity in aquatic microbes: Placing protists on the map. Res. Microbiol. 2015, 166, 831–841. [Google Scholar] [CrossRef]
  29. Bai, Y.; Wang, R.; Song, W.; Suzuki, T.; Hu, X. Redescription of five tintinnine ciliates (Alveolata: Ciliophora: Oligotrichea) from coastal waters of Qingdao, China. Mar. Life Sci. Technol. 2020, 2, 209–221. [Google Scholar] [CrossRef]
  30. Dale, T.; Dahl, E. Mass occurrence of planktonic oligotrichous ciliates in a bay in southern Norway. J. Plankton Res. 1987, 9, 871–879. [Google Scholar] [CrossRef]
  31. Edwards, E.; Burkill, P. Abundance, biomass and distribution of microzooplankton in the Irish Sea. J. Plankton Res. 1995, 17, 771–782. [Google Scholar] [CrossRef]
  32. Huang, H.; Yang, J.; Huang, S.; Gu, B.; Wang, Y.; Wang, L.; Jiao, N.; Xu, D. Spatial distribution of planktonic ciliates in the western Pacific Ocean: Along the transect from Shenzhen (China) to Pohnpei (Micronesia). Mar. Life Sci. Technol. 2021, 3, 103–115. [Google Scholar] [CrossRef]
  33. Liu, W.; Shin, M.K.; Yi, Z.; Tan, Y. Progress in studies on the diversity and distribution of planktonic ciliates (Protista, Ciliophora) in the South China Sea. Mar. Life Sci. Technol. 2021, 3, 28–43. [Google Scholar] [CrossRef]
  34. Zhu, C.; Liu, W.; Li, X.; Xu, Y.; El-Serehy, H.A.; Al-Farraj, S.A.; Ma, H.; Stoeck, T.; Yi, Z. High salinity gradients and intermediate spatial scales shaped similar biogeographical and co-occurrence patterns of microeukaryotes in a tropical freshwater-saltwater ecosystem. Environ. Microbiol. 2021, 23, 4778–4796. [Google Scholar] [CrossRef] [PubMed]
  35. Xu, J.; Wilkinson, M.; Chen, M.; Zhang, Q.; Yang, R.; Yi, Z. Concatenated data and dense taxon sampling clarify phylogeny and ecological transitions within Hypotricha. Zool. Scr. 2021, 50, 125–139. [Google Scholar] [CrossRef]
  36. Gao, F.; Warren, A.; Zhang, Q.; Gong, J.; Miao, M.; Sun, P.; Xu, D.; Huang, J.; Yi, Z.; Song, W. The all-data-based evolutionary hypothesis of ciliated protists with a revised classification of the phylum Ciliophora (Eukaryota, Alveolata). Sci. Rep. 2016, 6, 24874. [Google Scholar] [CrossRef] [Green Version]
  37. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  38. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef]
  39. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 22 August 2021).
  40. Paradis, E.; Claude, J.; Strimmer, K. APE: Analyses of phylogenetics and evolution in R language. Bioinformatics 2004, 20, 289–290. [Google Scholar] [CrossRef] [Green Version]
  41. Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 2003, 14, 927–930. [Google Scholar] [CrossRef]
  42. Clarke, K.R.; Somerfield, P.J.; Gorley, R.N. Testing of null hypotheses in exploratory community analyses: Similarity profiles and biota-environment linkage. J. Exp. Mar. Biol. Ecol. 2008, 366, 56–69. [Google Scholar] [CrossRef]
  43. Penn, O.; Privman, E.; Ashkenazy, H.; Landan, G.; Graur, D.; Pupko, T. GUIDANCE: A web server for assessing alignment confidence scores. Nucleic Acids Res. 2010, 38 (Suppl. S2), W23–W28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gouy, M.; Guindon, S.; Gascuel, O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 2010, 27, 221–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In Proceedings of the 2010 Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; pp. 1–8. [Google Scholar] [CrossRef] [Green Version]
  46. Posada, D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef] [PubMed]
  47. Rambaut, A. FigTree v1.4.2, 2012. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 27 August 2021).
  48. Maddison, W.P.; Maddison, D.R. Mesquite: A Modular System for Evolutionary Analysis. Version 3.6. 2018. Available online: http://www.mesquiteproject.org/ (accessed on 30 August 2021).
  49. Kopylova, E.; Navas-Molina, J.A.; Mercier, C.; Xu, Z.; Mahé, F.; He, Y.; Zhou, H.; Rognes, T.; Caporaso, J.G.; Knight, R. Open-source sequence clustering methods improve the state of the art. mSystems 2016, 1, e00003-15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Nearing, J.T.; Douglas, G.M.; Comeau, A.M.; Langille, M.G. Denoising the Denoisers: An independent evaluation of microbiome sequence error-correction approaches. PeerJ 2018, 6, e5364. [Google Scholar] [CrossRef] [Green Version]
  51. Xue, Z.; Kable, M.E.; Marco, M.L. Impact of DNA sequencing and analysis methods on 16S rRNA gene bacterial community analysis of dairy products. mSphere 2018, 3, e00410-18. [Google Scholar] [CrossRef] [Green Version]
  52. Liu, W.; McManus, G.B.; Lin, X.; Huang, H.; Zhang, W.; Tan, Y. Distribution patterns of ciliate diversity in the South China Sea. Front. Microbiol. 2021, 12, 689688. [Google Scholar] [CrossRef]
  53. Zhao, F.; Filker, S.; Xu, K.; Li, J.; Zhou, T.; Huang, P. Effects of intragenomic polymorphism in the SSU rRNA gene on estimating marine microeukaryotic diversity: A test for ciliates using single-cell high-throughput DNA sequencing. Limnol. Oceanogr. Methods 2019, 17, 533–543. [Google Scholar] [CrossRef]
  54. Huang, S.; Li, K.; Jiang, G.; Lu, J. Influence of salinity fronts on nutrition substance and biology in Pearl River estuary. Water Resour. Prot. 2011, 27, 18–25. [Google Scholar] [CrossRef]
  55. Zhang, X.; Shi, Z.; Huang, X.; Ye, F.; Liu, Q. Phytoplankton abundance and size-fractionated structure in three contrasting periods in the Pear River Estuary. J. Mar. Res. 2013, 71, 187–210. [Google Scholar] [CrossRef]
  56. Zhu, J.; Hong, Y.; Zada, S.; Hu, Z.; Wang, H. Spatial variability and co-acclimation of phytoplankton and bacterioplankton communities in the Pearl River Estuary, China. Front. Microbiol. 2018, 9, 2503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Liu, J.; Fu, B.; Yang, H.; Zhao, M.; He, B.; Zhang, X. Phylogenetic shifts of bacterioplankton community composition along the Pearl Estuary: The potential impact of hypoxia and nutrients. Front. Microbiol. 2015, 6, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Qiu, D.; Huang, L.; Zhang, J.; Lin, S. Phytoplankton dynamics in and near the highly eutrophic Pearl River Estuary, South China Sea. Cont. Shelf. Res. 2010, 30, 177–186. [Google Scholar] [CrossRef]
  59. Zhu, C.; Lu, K.; Yi, Z. Molecular diversities of plankotic microbial eukaryotes in the pearl river and their relationship with water environment (in Chinese). Acta Hydrobiol. Sin. 2020, 44, 187–196. [Google Scholar] [CrossRef]
  60. Hu, J.; Chivas, A.R. Molecular biomarker evidence of origins and transport of organic matter in sediments of the Pearl River estuary and adjacent South China Sea. Appl. Geochem. 2009, 24, 1666–1676. [Google Scholar] [CrossRef]
  61. Webster, G.; O’Sullivan, L.A.; Meng, Y.; Williams, A.S.; Sass, A.M.; Watkins, A.J.; Parkes, R.J.; Weightman, A.J. Archaeal community diversity and abundance changes along a natural salinity gradient in estuarine sediments. FEMS Microbiol. Ecol. 2015, 91, 1–18. [Google Scholar] [CrossRef]
  62. Logares, R.; Bråte, J.; Bertilsson, S.; Clasen, J.L.; Shalchian-Tabrizi, K.; Rengefors, K. Infrequent marine–freshwater transitions in the microbial world. Trends Microbiol. 2009, 17, 414–422. [Google Scholar] [CrossRef]
  63. Mai, Y.; Lai, Z.; Li, X.; Peng, S.; Wang, C. Structural and functional shifts of bacterioplanktonic communities associated with spatiotemporal gradients in river outlets of the subtropical Pearl River Estuary, South China. Mar. Pollut. Bull. 2018, 136, 309–321. [Google Scholar] [CrossRef]
  64. Dunthorn, M.; Otto, J.; Berger, S.A.; Stamatakis, A.; Mahé, F.; Romac, S.; de Vargas, C.; Audic, S.; Consortium, B.; Stock, A. Placing environmental next-generation sequencing amplicons from microbial eukaryotes into a phylogenetic context. Mol. Biol. Evol. 2014, 31, 993–1009. [Google Scholar] [CrossRef] [Green Version]
  65. Sun, P.; Clamp, J.; Xu, D.; Huang, B.; Shin, M.K. An integrative approach to phylogeny reveals patterns of environmental distribution and novel evolutionary relationships in a major group of ciliates. Sci. Rep. 2016, 6, 21695. [Google Scholar] [CrossRef]
  66. Bachy, C.; Gómez, F.; López-García, P.; Dolan, J.R.; Moreira, D. Molecular phylogeny of tintinnid ciliates (Tintinnida, Ciliophora). Protist 2012, 163, 873–887. [Google Scholar] [CrossRef] [PubMed]
  67. Zhao, Y.; Yi, Z.; Warren, A.; Song, W. Species delimitation for the molecular taxonomy and ecology of the widely distributed microbial eukaryote genus Euplotes (Alveolata, Ciliophora). Proc. R. Soc. B 2018, 285, 20172159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Lu, K.; Liu, W.; Warren, A.; Xu, Y.; Zhu, C.; Zhao, Y.; Yi, Z. Diversity of oligotrich ciliates (Ciliophora, Spirotrichea) in the northern coast of South China Sea as revealed in LSU rDNA sequences. J. Oceanol. Limnol. 2020, 38, 156–167. [Google Scholar] [CrossRef]
  69. Agatha, S. New observations on the tontoniid ciliate Spirotontonia grandis (Suzuki and Han, 2000) Agatha, 2004 (Ciliophora, Oligotrichida, Tontoniidae); comparison with the similar Laboea strobila. Eur. J. Protistol. 2004, 40, 295–301. [Google Scholar] [CrossRef]
  70. Agatha, S. Redescription of the tintinnid ciliate Tintinnopsis fimbriata Meunier, 1919 (Spirotricha, Choreotrichida) from coastal waters of northern Germany. Denisia 2008, 23, 261–272. [Google Scholar] [PubMed]
  71. Agatha, S.; Riedel-Lorjé, J.C. Morphology, infraciliature, and ecology of Halteriids and Strombidiids (Ciliophora, Oligotrichea) from coastal brackish water basins. Arch. Protistenkd. 1997, 148, 445–459. [Google Scholar] [CrossRef]
  72. Agatha, S.; Strüder-Kypke, M.C.; Beran, A.; Lynn, D.H. Pelagostrobilidium neptuni (Montagnes and Taylor, 1994) and Strombidium biarmatum nov. spec. (Ciliophora, Oligotrichea): Phylogenetic position inferred from morphology, ontogenesis, and gene sequence data. Eur. J. Protistol. 2005, 41, 65–83. [Google Scholar] [CrossRef] [Green Version]
  73. Agatha, S.; Tsai, S.F. Redescription of the tintinnid Stenosemella pacifica Kofoid and Campbell, 1929 (Ciliophora, Spirotricha) based on live observation, protargol impregnation, and scanning electron microscopy. J. Eukaryot. Microbiol. 2008, 55, 75–85. [Google Scholar] [CrossRef] [Green Version]
  74. Bai, Y.; Wang, R.; Liu, W.; Warren, A.; Zhao, Y.; Hu, X. Redescriptions of three tintinnine ciliates (Ciliophora: Tintinnina) from coastal waters in China based on lorica features, cell morphology, and rDNA sequence data. Eur. J. Protistol. 2020, 72, 125659. [Google Scholar] [CrossRef]
  75. Bai, Y.; Wang, R.; Song, W.; Li, L.; Santoferrara, L.F.; Hu, X. Three redescriptions in Tintinnopsis (Protista: Ciliophora: Tintinnina) from coastal waters of China, with cytology and phylogenetic analyses based on ribosomal RNA genes. BMC Microbiol. 2020, 20, 374. [Google Scholar] [CrossRef]
  76. Bardele, C.F.; Stockmann, N.; Agatha, S. Some ultrastructural features of the planktonic freshwater ciliate Limnostrombidium viride (Alveolata, Ciliophora, Oligotrichida) and improved diagnoses of oligotrich taxa. Acta Protozool. 2019, 57, 169–193. [Google Scholar] [CrossRef] [PubMed]
  77. Bernard, C.; Fenchel, T. Chemosensory behaviour of Strombidium purpureum, an anaerobic oligotrich with endosymbiotic purple non-sulphur bacteria. J. Eukaryot. Microbiol. 1994, 41, 391–396. [Google Scholar] [CrossRef]
  78. Cariou, J.B.; Dolan, J.; Dallot, S. A preliminary study of tintinnid diversity in the NW Mediterranean Sea. J. Plankton Res. 1999, 21, 1065–1075. [Google Scholar] [CrossRef] [Green Version]
  79. Cedrola, F.; Senra, M.V.X.; D’Agosto, M.; Dias, R.J.P. Phylogenetic analyses support validity of genus Eodinium (Ciliophora, Entodiniomorphida, Ophryoscolecidae). J. Eukaryot. Microbiol. 2017, 64, 242–247. [Google Scholar] [CrossRef] [PubMed]
  80. Chen, P.C.; Chiang, K.P.; Tsai, S.F. Pelagostrobilidium liui n. sp. (Ciliophora, Choreotrichida) from the coastal waters of northeastern Taiwan and an improved description of Pelagostrobilidium minutum Liu et al., 2012. J. Eukaryot. Microbiol. 2017, 64, 579–587. [Google Scholar] [CrossRef] [PubMed]
  81. Coats, D.W. Duboscquella cachoni n. sp. a parasitic dinoflagellate lethal to its tintinnine host Eutintinnus pectinis. J. Eukaryot. Microbiol. 1988, 35, 607–617. [Google Scholar] [CrossRef]
  82. Coats, D.W.; Heinbokel, J.A. Study of reproduction and other life cycle phenomena in planktonic protists using an acridine orange fluorescence technique. Mar. Biol. 1982, 67, 71–79. [Google Scholar] [CrossRef]
  83. de Cao, M.B.; Beigt, D.; Piccolo, C. Temporal variability of diversity and biomass of tintinnids (Ciliophora) in a southwestern Atlantic temperate estuary. J. Plankton Res. 2005, 27, 1103–1111. [Google Scholar] [CrossRef] [Green Version]
  84. de Cao, M.B. Abundance and species composition of Tintinnina (Ciliophora) in Bahía Blanca estuary, Argentina. Estuar. Coast. Shelf Sci. 1992, 34, 295–303. [Google Scholar] [CrossRef]
  85. Santoferrara, L.F.; Grattepanche, J.D.; Katz, L.A.; McManus, G.B. Patterns and processes in microbial biogeography: Do molecules and morphologies give the same answers? ISME J. 2016, 10, 1779–1790. [Google Scholar] [CrossRef]
  86. Feng, M.; Zhang, W.; Wang, W.; Zhang, G.; Xiao, T.; Xu, H. Can tintinnids be used for discriminating water quality status in marine ecosystems? Mar. Pollut. Bull. 2015, 101, 549–555. [Google Scholar] [CrossRef] [PubMed]
  87. Feng, M.; Zhang, W.; Yu, Y.; Xiao, T.; Sun, J. Horizontal distribution of tintinnids in the western South China Sea during summer 2007 (in Chinese). Trop. Oceanogr. 2013, 32, 86–92. [Google Scholar] [CrossRef]
  88. Gao, S.; Gong, J.; Lynn, D.; Lin, X.; Song, W. An updated phylogeny of oligotrich and choreotrich ciliates (Protozoa, Ciliophora, Spirotrichea) with representative taxa collected from Chinese coastal waters. Syst. Biodivers. 2009, 7, 235–242. [Google Scholar] [CrossRef]
  89. Graziano, C. On the ecology of tintinnids (Ciliophora: Oligotrichida) in the North Irish Sea. Estuar. Coast. Shelf Sci. 1989, 29, 233–245. [Google Scholar] [CrossRef]
  90. Gruber, M.S.; Strüder-Kypke, M.; Agatha, S. Redescription of Tintinnopsis everta Kofoid and Campbell 1929 (Alveolata, Ciliophora, Tintinnina) based on taxonomic and genetic analyses—Discovery of a new complex ciliary pattern. J. Eukaryot. Microbiol. 2018, 65, 484–504. [Google Scholar] [CrossRef] [Green Version]
  91. Urrutxurtu, I. Seasonal succession of tintinnids in the Nervión River estuary, Basque Country, Spain. J. Plankton Res. 2004, 26, 307–314. [Google Scholar] [CrossRef]
  92. Jeong, H.J.; Kim, J.S.; Kim, S.; Song, J.; Lee, I.; Lee, G.H. Strombidinopsis jeokjo n. sp. (Ciliophora: Choreotrichida) from the coastal waters off western Korea: Morphology and small subunit ribosomal DNA gene sequence. J. Eukaryot. Microbiol. 2004, 51, 451–455. [Google Scholar] [CrossRef]
  93. Jiang, J.; Xing, Y.; Miao, M.; Shao, C.; Warren, A.; Song, W. Two new marine ciliates, Caryotricha rariseta n. sp. and Discocephalus pararotatorius n. sp. (Ciliophora, Spirotrichea), with phylogenetic analyses inferred from the small subunit rRNA gene sequences. J. Eukaryot. Microbiol. 2013, 60, 388–398. [Google Scholar] [CrossRef]
  94. Jung, J.H.; Moon, J.H.; Park, K.M.; Kim, S.; Dolan, J.R.; Yang, E.J. Novel insights into the genetic diversity of Parafavella based on mitochondrial CO1 sequences. Zool. Scr. 2018, 47, 743–755. [Google Scholar] [CrossRef]
  95. Jung, J.H.; Choi, J.M.; Coats, D.W.; Kim, Y.O. Euduboscquella costata n. sp.(Dinoflagellata, Syndinea), an intracellular parasite of the ciliate Schmidingerella arcuata: Morphology, molecular phylogeny, life cycle, prevalence, and infection intensity. J. Eukaryot. Microbiol. 2016, 63, 3–15. [Google Scholar] [CrossRef]
  96. Jyothibabu, R.; Madhu, N.; Jayalakshmi, K.; Balachandran, K.; Shiyas, C.; Martin, G.; Nair, K. Impact of freshwater influx on microzooplankton mediated food web in a tropical estuary (Cochin backwaters–India). Estuar. Coast. Shelf Sci. 2006, 69, 505–518. [Google Scholar] [CrossRef]
  97. Kamiyama, T. Growth and grazing responses of tintinnid ciliates feeding on the toxic dinoflagellate Heterocapsacircularisquama. Mar. Biol. 1997, 128, 509–515. [Google Scholar] [CrossRef]
  98. Kazama, T.; Ishida, S.; Shimano, S.; Urabe, J. Discrepancy between conventional morphological systematics and nuclear phylogeny of tintinnids (Ciliophora: Choreotrichia). Plankton Benthos Res. 2012, 7, 111–125. [Google Scholar] [CrossRef] [Green Version]
  99. Kim, J.H.; Omar, A.; Jung, J.H. Morphology and phylogeny of the soil ciliate Parabistichella multilineae sp. nov. (Protozoa: Ciliophora: Hypotricha). Zool. Sci. 2019, 36, 242–249. [Google Scholar] [CrossRef] [PubMed]
  100. Kim, J.S.; Jeong, H.J.; Strüder-Kypke, M.C.; Lynn, D.H.; Kim, S.; Kim, J.H.; Lee, S.H. Parastrombidinopsis shimi n. gen., n. sp. (Ciliophora: Choreotrichia) from the coastal waters of Korea: Morphology and small subunit ribosomal DNA sequence. J. Eukaryot. Microbiol. 2005, 52, 514–522. [Google Scholar] [CrossRef] [PubMed]
  101. Kim, S.Y.; Choi, J.K.; Dolan, J.R.; Shin, H.C.; Lee, S.; Yang, E.J. Morphological and ribosomal DNA-based characterization of six antarctic ciliate morphospecies from the Amundsen Sea with phylogenetic analyses. J. Eukaryot. Microbiol. 2013, 60, 497–513. [Google Scholar] [CrossRef]
  102. Kim, S.Y.; Yang, E.J.; Gong, J.; Choi, J.K. Redescription of Favella ehrenbergii (Claparède and Lachmann, 1858) Jörgensen, 1924 (Ciliophora: Choreotrichia), with phylogenetic analyses based on small subunit rRNA gene sequences. J. Eukaryot. Microbiol. 2010, 57, 460–467. [Google Scholar] [CrossRef]
  103. Kim, Y.O.; Kim, S.Y.; Lee, W.J.; Choi, J.K. New observations on the choreotrich ciliate Strombidinopsis acuminata Fauré-Fremiet 1924, and comparison with Strombidinopsis jeokjo Jeong et al., 2004. J. Eukaryot. Microbiol. 2010, 57, 48–55. [Google Scholar] [CrossRef]
  104. Kim, Y.O.; Shin, K.; Jang, P.G.; Choi, H.W.; Noh, J.H.; Yang, E.J.; Kim, E.; Jeon, D. Tintinnid species as biological indicators for monitoring intrusion of the warm oceanic waters into Korean coastal waters. Ocean Sci. J. 2012, 47, 161–172. [Google Scholar] [CrossRef]
  105. Kršinić, F. On vertical distribution of tintinnines (Ciliata, Oligotrichida, Tintinnina) in the open waters of the South Adriatic. Deep-Sea Res. Part B Oceanogr. Lit. Rev. 1982, 29, 83–90. [Google Scholar] [CrossRef]
  106. Lee, E.S.; Kim, Y.O.; Agatha, S.; Jung, J.H.; Xu, D.; Shin, M.K. Revision of Strombidium paracalkinsi (Ciliophora: Oligotrichea: Oligotrichia), with comparison of Strombidiids bearing thigmotactic membranelles. J. Eukaryot. Microbiol. 2015, 62, 400–409. [Google Scholar] [CrossRef] [PubMed]
  107. Lee, E.S.; Shin, M.K.; Kim, Y.O. Morphological descriptions of four oligotrich ciliates (Ciliophora: Oligotrichia) from southern coast of Korea. Korean J. Syst. Zool. 2011, 27, 131–141. [Google Scholar] [CrossRef] [Green Version]
  108. Lee, K.W.; Choi, Y.U. Population growth of a tropical tintinnid, Metacylis tropica on different temperature, salinity and diet. J. Korea Acad. Industr. Coop. Soc. 2016, 17, 322–328. [Google Scholar] [CrossRef]
  109. Li, F.; Huang, Y.; Yu, Y.; Liu, W.; Lin, X. Taxonomy and phylogeny of a new marine planktonic ciliate, Strombidium pseudorapulum sp. n. (Protozoa, Ciliophora, Oligotrichia) (in Chinese). J. Ocean Univ. 2020, 19, 954–960. [Google Scholar] [CrossRef]
  110. Li, H.; Xu, Z.; Zhang, W.; Wang, S.; Zhang, G.; Xiao, T. Boreal tintinnid assemblage in the Northwest Pacific and its connection with the Japan Sea in summer 2014. PLoS ONE 2016, 11, e0153379. [Google Scholar] [CrossRef]
  111. Li, L.; Shao, C.; Song, W.; Lynn, D.H.; Chen, Z.; Shin, M.K. Does Kiitricha (Protista, Ciliophora, Spirotrichea) belong to Euplotida or represent a primordial spirotrichous taxon? With suggestion to establish a new subclass Protohypotrichia. Int. J. Syst. Evol. Microbiol. 2009, 59, 439–446. [Google Scholar] [CrossRef]
  112. Liu, W.; Xu, D.; Lin, X.; Li, J.; Gong, J.; Al-Rasheid, K.A.; Song, W. Novistrombidium sinicum n. sp. and Novistrombidium orientale n. sp.(Protozoa: Ciliophora): Two new oligotrich ciliates from a mangrove wetland, South China. J. Eukaryot. Microbiol. 2009, 56, 459–465. [Google Scholar] [CrossRef]
  113. Liu, W.; Yi, Z.; Li, J.; Warren, A.; Al-Farraj, S.A.; Lin, X. Taxonomy, morphology and phylogeny of three new oligotrich ciliates (Protozoa, Ciliophora, Oligotrichia) from southern China. Int. J. Syst. Evol. Microbiol. 2013, 63, 4805–4817. [Google Scholar] [CrossRef]
  114. Liu, W.; Yi, Z.; Lin, X.; Al-Rasheid, K.A. Morphologic and molecular data suggest that Lynnella semiglobulosa ng, n. sp. represents a new family within the subclass Choreotrichia (Ciliophora, Spirotrichea). J. Eukaryot. Microbiol. 2011, 58, 43–49. [Google Scholar] [CrossRef]
  115. Liu, W.; Yi, Z.; Lin, X.; Warren, A.; Song, W. Phylogeny of three choreotrich genera (Protozoa, Ciliophora, Spirotrichea), with morphological, morphogenetic and molecular investigations on three strobilidiid species. Zool. Scr. 2012, 41, 417–434. [Google Scholar] [CrossRef]
  116. Liu, W.; Yi, Z.; Lin, X.; Li, J.; Al-Farraj, S.A.; Al-Rasheid, K.A.; Song, W. Morphology and molecular phylogeny of three new oligotrich ciliates (Protozoa, Ciliophora) from the South China Sea. Zool. J. Linn. Soc. 2015, 174, 653–665. [Google Scholar] [CrossRef] [Green Version]
  117. Liu, W.; Yi, Z.; Warren, A.; Al-Rasheid, K.A.; Al-Farraj, S.A.; Lin, X.; Song, W. Taxonomy, morphology and molecular systematics of a new oligotrich ciliate, Williophrya maedai gen. nov., sp. nov., with redescriptions of Strombidium basimorphum and Pseudotontonia simplicidens (Protozoa, Ciliophora, Oligotrichia). Syst. Biodivers. 2011, 9, 247–258. [Google Scholar] [CrossRef]
  118. Liu, W.; Yi, Z.; Xu, D.; Clamp, J.C.; Li, J.; Lin, X.; Song, W. Two new genera of planktonic ciliates and insights into the evolution of the family Strombidiidae (Protista, Ciliophora, Oligotrichia). PLoS ONE 2015, 10, e0131726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  119. Lynn, D.; Montagnes, D.; Dale, T.; Gilron, G.; Strom, S. A reassessment of the genus Strombidinopsis (Ciliophora, Choreotrichida) with descriptions of four new planktonic species and remarks in its taxonomy and phylogeny. J. Mar. Biol. Assoc. UK 1991, 71, 597–612. [Google Scholar] [CrossRef]
  120. Miao, M.; Shao, C.; Jiang, J.; Li, L.; Stoeck, T.; Song, W. Caryotricha minuta (Xu et al., 2008) nov. comb., a unique marine ciliate (Protista, Ciliophora, Spirotrichea), with phylogenetic analysis of the ambiguous genus Caryotricha inferred from the small-subunit rRNA gene sequence. Int. J. Syst. Evol. Microbiol. 2009, 59, 430–438. [Google Scholar] [CrossRef]
  121. Modeo, L.; Petroni, G.; Rosati, G.; Montagnes, D.J. A multidisciplinary approach to describe protists: Redescriptions of Novistrombidium testaceum Anigstein 1914 and Strombidium inclinatum Montagnes, Taylor, and Lynn 1990 (Ciliophora, Oligotrichia). J. Eukaryot. Microbiol. 2003, 50, 175–189. [Google Scholar] [CrossRef]
  122. Modigh, M.; Castaldo, S.; Saggiomo, M.; Santarpia, I. Distribution of tintinnid species from 42° N to 43° S through the Indian Ocean. Hydrobiologia 2003, 503, 251–262. [Google Scholar] [CrossRef]
  123. Park, T.; Yu, Z. Do ruminal ciliates select their preys and prokaryotic symbionts? Front. Microbiol. 2018, 9, 1710. [Google Scholar] [CrossRef] [Green Version]
  124. Pierce, R.W.; Turner, J.T. Plankton studies in buzzards Bay, Massachusetts, USA. IV. Tintinnids, 1987 to 1988. Mar. Ecol-Prog. Ser. 1994, 112, 235–240. [Google Scholar] [CrossRef]
  125. Putt, M. Metabolism of photosynthate in the chloroplast-retaining ciliate Laboea strobila. Mar. Ecol.-Prog. Ser. 1990, 60, 271–282. [Google Scholar] [CrossRef]
  126. Sitran, R.; Bergamasco, A.; Decembrini, F.; Guglielmo, L. Temporal succession of tintinnids in the northern Ionian Sea, Central Mediterranean. J. Plankton Res. 2007, 29, 495–508. [Google Scholar] [CrossRef]
  127. Santoferrara, L.F.; Alder, V.V.; McManus, G.B. Phylogeny, classification and diversity of Choreotrichia and Oligotrichia (Ciliophora, Spirotrichea). Mol. Phylogenet. Evol. 2017, 112, 12–22. [Google Scholar] [CrossRef] [PubMed]
  128. Santoferrara, L.F.; McManus, G.B.; Alder, V.A. Phylogeny of the order Tintinnida (Ciliophora, Spirotrichea) inferred from small-and large-subunit rRNA genes. J. Eukaryot. Microbiol. 2012, 59, 423–426. [Google Scholar] [CrossRef] [PubMed]
  129. Santoferrara, L.F.; McManus, G.B.; Alder, V.A. Utility of genetic markers and morphology for species discrimination within the order Tintinnida (Ciliophora, Spirotrichea). Protist 2013, 164, 24–36. [Google Scholar] [CrossRef] [PubMed]
  130. Santoferrara, L.F.; Tian, M.; Alder, V.A.; McManus, G.B. Discrimination of closely related species in tintinnid ciliates: New insights on crypticity and polymorphism in the genus Helicostomella. Protist 2015, 166, 78–92. [Google Scholar] [CrossRef] [PubMed]
  131. Smith, S.A.; Song, W.; Gavrilova, N.A.; Kurilov, A.V.; Liu, W.; McManus, G.B.; Santoferrara, L.F. Dartintinnus alderae n. g., n. sp., a brackish water tintinnid (Ciliophora, Spirotrichea) with dual-ended lorica collapsibility. J. Eukaryot. Microbiol. 2018, 65, 400–411. [Google Scholar] [CrossRef] [PubMed]
  132. Song, W.; Li, J.; Liu, W.; Al-Rasheid, K.A.; Hu, X.; Lin, X. Taxonomy and molecular phylogeny of four Strombidium species, including description of S. pseudostylifer sp. nov. (Ciliophora, Oligotrichia). Syst. Biodivers. 2015, 13, 76–92. [Google Scholar] [CrossRef]
  133. Song, W.; Li, J.; Liu, W.; Jiang, J.; Al-Rasheid, K.A.; Hu, X. Taxonomy, morphology and molecular systematics of three oligotrich ciliates, including a description of Apostrombidium parakielum spec. nov.(Ciliophora, Oligotrichia). Int. J. Syst. Evol. Microbiol. 2013, 63, 1179–1191. [Google Scholar] [CrossRef]
  134. Song, W.; Pan, B.; El-Serehy, H.A.; Al-Farraj, S.A.; Liu, W.; Li, L. Morphology and molecular phylogeny of two freshwater oligotrich ciliates (Protozoa, Ciliophora, Oligotrichia), Pelagostrombidium fallax (Zacharias, 1895) Krainer, 1991 and Limnostrombidium viride (Stein, 1867) Krainer, 1995, with brief notes on stomatogenesis. J. Eukaryot. Microbiol. 2020, 67, 232–244. [Google Scholar] [CrossRef]
  135. Song, W.; Xu, D.; Zhang, Q.; Liu, W.; Warren, A.; Song, W. Taxonomy and phylogeny of two poorly studied genera of marine oligotrich ciliates including descriptions of two new species: Cyrtostrombidium paraboreale sp. n. and Apostrombidium orientale sp. n.(Ciliophora: Spirotrichea). Eur. J. Protistol. 2019, 70, 1–16. [Google Scholar] [CrossRef]
  136. Song, W.; Wang, L.; Li, L.; Al-Farraj, S.A.; Aleidan, A.; Smith, S.; Hu, X. Morphological characterizations of four species of Parallelostrombidium (Ciliophora, Oligotrichia), with a note on the phylogeny of the genus. J. Eukaryot. Microbiol. 2018, 65, 679–693. [Google Scholar] [CrossRef] [PubMed]
  137. Song, W.; Wang, M.; Warren, A. Redescriptions of three marine ciliates, Strombidium elegans Florentin, 1901, Strombidium sulcatum Claparede & Lachmann, 1859 and Heterostrombidium paracalkinsi Lei, Xu & Song, 1999 (Ciliophora, Oligotrichida). Eur. J. Protistol. 2000, 36, 327–342. [Google Scholar] [CrossRef]
  138. Song, W.; Zhao, X.; Liu, W.; Hu, X.; Al-Farraj, S.A.; Al-Rasheid, K.A.; Song, W.; Warren, A. Biodiversity of oligotrich ciliates in the South China Sea: Description of three new Strombidium species (Protozoa, Ciliophora, Oligotrichia) with phylogenetic analyses. Syst. Biodivers. 2015, 13, 608–623. [Google Scholar] [CrossRef]
  139. Song, Y.; Liu, Y.; Pan, B.; Luo, X.; Song, W.; Warren, A. Morphological studies on four brackish water ciliates of the class Spirotrichea (Protista, Ciliophora). J. Ocean Univ. 2019, 18, 663–674. [Google Scholar] [CrossRef]
  140. Strüder-Kypke, M.C.; Lynn, D.H. Sequence analyses of the small subunit rRNA gene confirm the paraphyly of oligotrich ciliates sensu lato and support the monophyly of the subclasses Oligotrichia and Choreotrichia (Ciliophora, Spirotrichea). J. Zool. 2003, 260, 87–97. [Google Scholar] [CrossRef]
  141. Tsai, S.F.; Chen, J.Y.; Chiang, K.P. Spirotontonia taiwanica n. sp.(Ciliophora: Oligotrichida) from the coastal waters of Northeastern Taiwan: Morphology and nuclear small subunit rDNA sequence. J. Eukaryot. Microbiol. 2010, 57, 429–434. [Google Scholar] [CrossRef]
  142. Tsai, S.F.; Chen, W.T.; Chiang, K.P. Phylogenetic position of the genus Cyrtostrombidium, with a description of Cyrtostrombidium paralongisomum nov. spec. and a redescription of Cyrtostrombidium longisomum Lynn & Gilron, 1993 (Protozoa, Ciliophora) based on live observation, protargol impregnation, and 18S rDNA sequences. J. Eukaryot. Microbiol. 2015, 62, 239–248. [Google Scholar] [CrossRef]
  143. Tsai, S.F.; Xu, D.; Chung, C.C.; Chiang, K.P. Parastrombidinopsis minima n. sp. (Ciliophora: Oligotrichia) from the coastal waters of northeastern Taiwan: Morphology and small subunit ribosomal DNA sequence. J. Eukaryot. Microbiol. 2008, 55, 567–573. [Google Scholar] [CrossRef]
  144. Wang, R.; Bai, Y.; Hu, T.; Xu, D.; Suzuki, T.; Hu, X. Integrative taxonomy and molecular phylogeny of three poorly known tintinnine ciliates, with the establishment of a new genus (Protista; Ciliophora; Oligotrichea). BMC Ecol. Evol. 2021, 21, 115. [Google Scholar] [CrossRef]
  145. Weisse, T.; Rammer, S. Pronounced ecophysiological clonal differences of two common freshwater ciliates, Coleps spetai (Prostomatida) and Rimostrombidium lacustris (Oligotrichida), challenge the morphospecies concept. J. Plankton Res. 2006, 28, 55–63. [Google Scholar] [CrossRef] [Green Version]
  146. Xu, D.; Song, W.; Lin, X.; Warren, A. On two marine oligotrich ciliates, Spirostrombidium agathae n. sp. and S. schizostomum (Kahl, 1932) n. comb. from China, with a key to the identification of seven well-characterized Spirostrombidium spp. (Ciliophora: Oligotrichida). Acta Protozool. 2006, 45, 433–442. [Google Scholar]
  147. Xu, D.; Song, W.; Warren, A. Morphology and infraciliature of two new species of marine oligotrich ciliates (Ciliophora: Oligotrichida) from China. J. Nat. Hist. 2006, 40, 1287–1299. [Google Scholar] [CrossRef]
  148. Xu, D.; Sun, P.; Clamp, J.C.; Ma, H.; Song, W. The establishment of a new oligotrich genus Varistrombidium gen. nov. and the morphology and phylogeny of a marine ciliate, Varistrombidium kielum (Maeda and Carey, 1985) nov. comb. (Protista, Ciliophora). Acta. Zootaxa. Sin. 2011, 36, 502–510. [Google Scholar]
  149. Xu, D.; Sun, P.; Warren, A.; Hoon Noh, J.; Lim Choi, D.; Kyoon Shin, M.; Ok Kim, Y. Phylogenetic investigations on ten genera of tintinnid ciliates (Ciliophora: Spirotrichea: Tintinnida), based on small subunit ribosomal RNA gene sequences. J. Eukaryot. Microbiol. 2013, 60, 192–202. [Google Scholar] [CrossRef] [PubMed]
  150. Zhang, Q.; Agatha, S.; Zhang, W.; Dong, J.; Yu, Y.; Jiao, N.; Gong, J. Three rDNA loci-based phylogenies of tintinnid ciliates (Ciliophora, Spirotrichea, Choreotrichida). J. Eukaryot. Microbiol. 2017, 64, 226–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  151. Zhang, W.; Sun, J.; Sun, S. Spatial distribution of tintinnidous ciliate CodoneHopsis mobilis (Protozoa, Ciliophora) in April, 1999 in Bohai Sea, China. Mar. Sci. 2004, 28, 67–69. [Google Scholar]
  152. Zhang, W.; Wang, R. Summertime ciliate and copepod nauplii distributions and micro-zooplankton herbivorous activity in the Laizhou Bay, Bohai Sea, China. Estuar. Coast. Shelf Sci. 2000, 51, 103–114. [Google Scholar] [CrossRef]
Microorganisms 10 00913 g001 550
Figure 1. The location of sampling regions. Sampling areas and their abbreviations are as follows: Daya Bay (DY1, DY2, DY5), Pearl River Estuary (PRE1–PRE10), Guangzhou (GZ1–GZ5), Shenzhen (SZ1–SZ4), Zhuhai (ZH1–ZH5), and Zhanjiang (ZJ1–ZJ4, ZJ6). The numbers shown in parentheses represent the salinity of the sampling sites (unit: ‰), and salinity values of freshwater, low-brackish, high-brackish, and marine areas are colored in blue, green, orange, and black, respectively. Adapted with permission from [34] published by John Wiley and Sons, 2022.
Figure 1. The location of sampling regions. Sampling areas and their abbreviations are as follows: Daya Bay (DY1, DY2, DY5), Pearl River Estuary (PRE1–PRE10), Guangzhou (GZ1–GZ5), Shenzhen (SZ1–SZ4), Zhuhai (ZH1–ZH5), and Zhanjiang (ZJ1–ZJ4, ZJ6). The numbers shown in parentheses represent the salinity of the sampling sites (unit: ‰), and salinity values of freshwater, low-brackish, high-brackish, and marine areas are colored in blue, green, orange, and black, respectively. Adapted with permission from [34] published by John Wiley and Sons, 2022.
Microorganisms 10 00913 g001
Microorganisms 10 00913 g002 550
Figure 2. Hierarchical clustering analysis of oligotrich (s.l.) community of 32 samples from six datasets. (a) DADA2–100; (b) SWARM–100; (c) UNOISE–100; (d) UPARSE–97; (e) UCLUST–97; (f) UCLUST–99. Samples from the same sampling area are indicated by the same color except for samples from PRE. Samples named PRE1–4 and PRE5–10 are indicated by different colors.
Figure 2. Hierarchical clustering analysis of oligotrich (s.l.) community of 32 samples from six datasets. (a) DADA2–100; (b) SWARM–100; (c) UNOISE–100; (d) UPARSE–97; (e) UCLUST–97; (f) UCLUST–99. Samples from the same sampling area are indicated by the same color except for samples from PRE. Samples named PRE1–4 and PRE5–10 are indicated by different colors.
Microorganisms 10 00913 g002
Microorganisms 10 00913 g003 550
Figure 3. The distributions of oligotrich (s.l.) communities from six datasets ((a). DADA2–100; (b). SWARM–100; (c). UNOISE–100; (d). UPARSE–97; (e). UCLUST–97; (f). UCLUST–99) based on Principal Co-ordinates Analysis (PCoA). Samples from the same sampling area are indicated by the same color except for samples from PRE. Samples named PRE1–4 and PRE5–10 are indicated by different colors. The ellipses represent 95% confidence intervals except for DY (too few samples to calculate an ellipse). P represents global significance among oligotrich (s.l.) communities from different sampling areas based on ADONIS. p values < 0.05 indicate significant differences among community structures.
Figure 3. The distributions of oligotrich (s.l.) communities from six datasets ((a). DADA2–100; (b). SWARM–100; (c). UNOISE–100; (d). UPARSE–97; (e). UCLUST–97; (f). UCLUST–99) based on Principal Co-ordinates Analysis (PCoA). Samples from the same sampling area are indicated by the same color except for samples from PRE. Samples named PRE1–4 and PRE5–10 are indicated by different colors. The ellipses represent 95% confidence intervals except for DY (too few samples to calculate an ellipse). P represents global significance among oligotrich (s.l.) communities from different sampling areas based on ADONIS. p values < 0.05 indicate significant differences among community structures.
Microorganisms 10 00913 g003
Microorganisms 10 00913 g004 550
Figure 4. Maximum likelihood (ML) tree based on SSU rDNA of 103 environmental sequences (ASVs in DADA2–100) and 197 sequences of identified species from GenBank (accession numbers in Supplementary Table S1). The tree shows the topology and transition pattern; support values (ML/BI) > 50%/0.50 are labeled. The number of collapsed sequences and ecological characteristics are shown in parentheses. F, freshwater; L, low-brackish; H, high-brackish; M, marine.
Figure 4. Maximum likelihood (ML) tree based on SSU rDNA of 103 environmental sequences (ASVs in DADA2–100) and 197 sequences of identified species from GenBank (accession numbers in Supplementary Table S1). The tree shows the topology and transition pattern; support values (ML/BI) > 50%/0.50 are labeled. The number of collapsed sequences and ecological characteristics are shown in parentheses. F, freshwater; L, low-brackish; H, high-brackish; M, marine.
Microorganisms 10 00913 g004
Table
Table 1. The number of total and oligotrich (s.l.) OTUs/ASVs in six datasets.
Table 1. The number of total and oligotrich (s.l.) OTUs/ASVs in six datasets.
DADA2–100SWARM–100UNOISE–100UPARSE–97UCLUST–97UCLUST–99
Number of total OTUs/ASVs3890665629153095761319,993
Number of oligotrich (s.l.) OTUs/ASVs1031378963104248
Table
Table 2. Number and proportion of oligotrich (s.l.) OTUs/ASVs in each sampling area from the six datasets. The proportions in parentheses indicate the ratio of the oligotrich (s.l.) OTU/ASV number in each sampling area to the total number of oligotrich (s.l.) OTUs/ASVs within each dataset.
Table 2. Number and proportion of oligotrich (s.l.) OTUs/ASVs in each sampling area from the six datasets. The proportions in parentheses indicate the ratio of the oligotrich (s.l.) OTU/ASV number in each sampling area to the total number of oligotrich (s.l.) OTUs/ASVs within each dataset.
Sampling AreasDADA2–100SWARM–100UNOISE–100UPARSE–97UCLUST–97UCLUST–99
DY6 (5.83%)17 (12.41%)14 (15.73%)15 (23.81%)14 (13.46%)17 (6.85%)
GZ31 (30.10%)54 (39.42%)49 (55.06%)31 (49.21%)39 (37.50%)88 (35.48%)
PRE47 (45.63%)76 (55.47%)58 (65.17%)45 (71.43%)67 (64.42%)137 (55.24%)
SZ8 (7.77%)15 (10.95%)12 (13.48%)10 (15.87%)16 (15.38%)30 (12.10%)
ZH23 (22.33%)40 (29.20%)40 (44.94%)28 (44.44%)27 (25.96%)51 (20.56%)
ZJ32 (31.07%)63 (45.99%)46 (51.69%)32 (50.79%)45 (43.27%)110 (44.35%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, J.; Han, J.; Su, H.; Zhu, C.; Quan, Z.; Wu, L.; Yi, Z. Diversity Patterns of Protists Are Highly Affected by Methods Disentangling Biological Variants: A Case Study in Oligotrich (s.l.) Ciliates. Microorganisms 2022, 10, 913. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10050913

AMA Style

Xu J, Han J, Su H, Zhu C, Quan Z, Wu L, Yi Z. Diversity Patterns of Protists Are Highly Affected by Methods Disentangling Biological Variants: A Case Study in Oligotrich (s.l.) Ciliates. Microorganisms. 2022; 10(5):913. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10050913

Chicago/Turabian Style

Xu, Jiahui, Jianlin Han, Hua Su, Changyu Zhu, Zijing Quan, Lei Wu, and Zhenzhen Yi. 2022. "Diversity Patterns of Protists Are Highly Affected by Methods Disentangling Biological Variants: A Case Study in Oligotrich (s.l.) Ciliates" Microorganisms 10, no. 5: 913. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10050913

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop