Next Article in Journal
Zymobacter palmae Pyruvate Decarboxylase is Less Effective Than That of Zymomonas mobilis for Ethanol Production in Metabolically Engineered Synechocystis sp. PCC6803
Previous Article in Journal
Genetic and Phenotypic Characterisation of a Saccharomyces cerevisiae Population of ‘Merwah’ White Wine
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 11
10.3390/microorganisms7110493
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

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

by
Zifeng Zhan
1,2,3,
Ju Li
1,3 and
Kuidong Xu
1,2,3,4,*
1
Laboratory of Marine Organism Taxonomy and Phylogeny, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266071, China
3
Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
4
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Microorganisms 2019, 7(11), 493; https://doi.org/10.3390/microorganisms7110493
Submission received: 22 September 2019 / Revised: 22 October 2019 / Accepted: 25 October 2019 / Published: 26 October 2019
(This article belongs to the Section Environmental Microbiology)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

DNA barcoding is an exciting tool to assess biodiversity, and the mitochondrial cytochrome c oxidase subunit 1 (COI) gene as a barcode has been successfully used in many animal taxa [1,2,3]. In the case of eukaryotic microorganisms, the Protist Working Group proposed a two-step barcoding pipeline for eukaryotic microorganisms: first, the hypervariable V4 region of small subunit ribosomal DNA (SSU-V4 rDNA) is applied to assign organisms to a high-ranking taxonomic group; second, a group-specific barcode is used for species-level assignments [4]. Many microbial eukaryotic groups, however, still lack a suitable barcode [5]. At the current stage, small subunit ribosomal DNA (SSU rDNA) is still the primary marker used to identify species and to examine their phylogenetic positions [6,7]. On the basis of the large sequence database, the SSU-V4 rDNA has been widely chosen as a metabarcode for high-throughput sequencing analyses to assess the environmental diversity of the eukaryotic microorganisms [8,9,10]. However, both the partial and complete SSU rDNA have been questioned as being too conserved to distinguish closely related species and to uncover cryptic species [9,11,12,13,14,15]. Thus, there is a pressing need to test other gene markers as suitable barcodes and metabarcodes.
Pseudokeronopsis Borror and Wicklow, 1983 is a conspicuous ciliate genus, in which most species have a colourful body and occur at high frequencies in marine or limnetic habitats worldwide [7,12,16,17,18,19,20]. Nowadays, there is growing interest in using Pseudokeronopsis in various research fields including chemical defense mechanisms, actin evolution, and morphogenesis [18,21,22,23]. Pseudokeronopsids are distinguished mainly by the morphological combination of ciliary pattern, cortical granules, and contractile vacuole [12,16,19]. However, Pseudokeronopsis is one of the most difficult genera of ciliates for species discrimination due to its high interspecific similarities [12,16,19,20]. Many populations were often misidentified as certain “well-known” species, for example, the type species Pseudokeronopsis rubra and Pseudokeronopsis carnea [12,16,19]. Consequently, many Pseudokeronopsis sequences in GenBank are dubious, and even the correct sequences of the type species are awaiting being submitted [9]. On the other hand, the inner phylogenetic relationships are still relatively unknown due to the incomplete or inaccurate morphological and molecular data [12,16,19,20,24]. Pseudokeronopsis frequently occurs in marine habitats, and 7 of 12 known species could be sampled in the same location, that being the coastal area of Qingdao, China (Table 1). Hence, the ciliate genus Pseudokeronopsis is an ideal taxon that can be utilized to test the potential of different markers for species delimitation, environmental diversity estimation, and phylogenetic reconstruction of the eukaryotic microorganisms.
No universal barcode has been accepted for the species delimitation of ciliates [14]. The mitochondrial COI and cytochrome b and the nuclear 5.8S rDNA combining the internal transcribed spacer regions (ITS1-5.8S-ITS2) and the hypervariable D1-D2 region of large subunit ribosomal DNA (LSU rDNA) have showed varying degrees of barcoding potential [5,13,14,15,25,26,27,28,29]. Among them, COI is most frequently tested as a candidate barcode, providing much better resolution at the intraspecific level than that of the nuclear genes [13,14,25]. Recently, the D1-D2 region of LSU (LSU-D1-D2) rDNA was proposed as a good candidate barcode for ciliates due to the fact that it shows clear barcoding gaps within Euplotes, Paramecium, and Tintinnida [5,9,15]. Due to the ITS1 and ITS2 regions having shown much higher rates of evolutionary changes than the code regions of the ribosomal subunits, the ITS1-5.8S-ITS2 and ITS2 regions are also popular markers for the ciliate identification and investigation of intraspecific genetic diversity [30,31,32,33]. However, the barcoding efficiency of the mitochondrial and nuclear markers was only assessed in several oligohymenophorean genera (e.g., Paramecium, Frontonia, and Tetrahymena) with congeneric and conspecific sequences [15,25,34]. No such assessment has been made for other groups, for example, the most diverse ciliate class Spirotrichea.
In the present study, we assessed the potential of the genetic distance of the ITS1-5.8S-ITS2, ITS2, LSU-D2, SSU-V4 rDNAs, and the mitochondrial COI, as well as the secondary structures of ITS2 and LSU-D2 for species discrimination of the spirotrichean genus Pseudokeronopsis. Subsequently, we revised all available Pseudokeronopsis sequences in GenBank on the basis of the barcoding analyses combined with morphological data. On the basis of the comparative analysis of the present and previously published data, we examined which markers were appropriate barcodes and metabarcodes for ciliates. Simultaneously, we compared the resolving power of the selected candidate barcodes and SSU rDNA as well as the concatenated nuclear and mitochondrial regions for the phylogenetic reconstruction. By combining the morphological and molecular data, we inferred the putative phylogenetic relationships for all known species of Pseudokeronopsis. The objective of this study was to shed more light on the screening of barcodes and metabarcodes for the eukaryotic microorganisms.

2. Materials and Methods

2.1. Specimen Collection and Identification

Pseudokeronopsis rubra, the type species of the genus, was collected in April 2017 from the Qingdao Port (36°05′17″ N, 120°19′21″ E). Pseudokeronopsis cf. songi was collected in April 2018 from the treated sewage outfall near Qingdao No. 1 Bathing Beach (36°03′25″ N, 120°20′10″ E) on the China coast of the Yellow Sea. Pseudokeronopsis erythrina was collected in November 2014 from the Sheyang Estuary (33°49′00″ N, 120°28′05″ E) on the China coast of the Yellow Sea. The genomic DNAs of P. carnea, P. pararubra, and P. songi were obtained from [12].
Cells from field samples were isolated and observed in vivo at 100–1000× magnifications using bright field and differential interference contrast microscopy. The ciliates were impregnated with protargol following [38] to reveal the ciliary and nuclear pattern for morphological identification. Terminology and systematics mainly followed [7,16].

2.2. DNA Extraction, PCR Amplification, and Sequencing

Monoclonal cultures for Pseudokeronopsis cf. songi and P. erythrina, and rough culture for P. rubra were established. Live cells were maintained in filtered seawater overnight, and then rinsed in filtered seawater four times. One cell of P. cf. songi and P. rubra, and three cells of P. erythrina were transferred into three microfuge tubes for DNA extraction using DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. PCR amplification for small subunit ribosomal DNA (SSU rDNA) was conducted using primers 16S-F (5′-AACCTGGTTGATCCTGCCAGT-3′) and 16S-R (5′-TGATCCTTCTGCAGGTTCACCTAC-3′) as follows: denaturation at 98 °C for 2 min, followed by 35 cycles of denaturation at 98 °C for 10 s, annealing at 58 °C for 12 s, extension at 72 °C for 20 s, and a final extension step at 72 °C for 5 min [39]. PCR amplification for internal transcribed spacer region ITS1-5.8S-ITS2 and a 5′-end region of large subunit ribosomal DNA (ITS1-5.8S-ITS2-5′LSU rDNA) sequences was conducted using the ITS-F (5′-GTAGGTGAACCTGCGGAAGGATCATTA-3′) [40] and CILI28S-1000 (5′-CATTCGGCAGGTGAGTTGTTACACTCC-3′) [41] primers that were complementary to conserved regions and encompassed the 3′-end of the SSU rDNA (27 bp), the entire ITS1-5.8S-ITS2 region, and the complete D1-D2 and partial D3 regions of the large subunit rRNA. The amplification cycle conditions were the same with that of SSU rDNA. Mitochondrial COI gene was amplified with the primer set (CiCOIFv2: 5′-GWTGRGCKATGATYACACC-3′; CiCOIRv2: 5′-ACCATRTACATATGATGWCC-3′), with the following conditions: 2 min at 94°C followed by 40 cycles (95 °C for 15 s, 53 °C for 30 s, and 72 °C for 60 s) and a final extension at 72 °C for 5 min [27]. PCR reactions were performed using I-5 2x High-Fidelity Master Mix DNA polymerase (TP001, Beijing Tsingke Biotech Co., Ltd.). PCR purification and sequencing were performed by Tsingke (Beijing, China).

2.3. Genetic Distances, Operated Taxonomic Unite Delimitation, and Phylogenetic Analyses

All the available SSU rDNA, ITS1-5.8S-ITS2, LSU rDNA, and COI sequences of Pseudokeronopsis spp. except for those considered as non-Pseudokeronopsis by [12] and the out-group species from related hypotrich genera were downloaded from GenBank (Table 1). On the basis of the morphological classification of [9], the previously sampled populations of P. rubra and P. carnea were misidentified and marked with asterisks in the tables and figures of the following analyses. These sequences and those noted as Pseudokeronopsis sp. or submitted without associated morphological publications were excluded from the initial statistics of intra- and interspecific distances.
Sequences of each region were aligned independently using MAFFT v.7 with the G-INS-i and Q-INS-i (secondary structure of RNA is considered) algorithms for the COI and ribosomal regions [42]. The SSU-V4, ITS1-5.8S-ITS2, ITS2, and LSU-D2 sequence datasets were constructed according to the previous secondary structure predictions by [43,44,45]. For the barcoding analysis, genetic distances between species/populations were calculated with MEGA 6.0 [46] using the Kimura two-parameter model. Considering that 7 of 12 known Pseudokeronopsis species can be sampled in the same location, that is, the coastal area of Qingdao, China (Table 1), we also compared the species diversities detected by the metabarcoding, barcoding and morphological methods. Using UPARSE for the metabarcoding analysis, sequence similarities of 97% and 99% of SSU-V4 were applied to delineate operated taxonomic unites (OTUs) of Pseudokeronopsis [47].
Phylogenetic analyses were performed on the datasets of SSU, SSU-V4, and LSU-D2 rDNA; ITS1-5.8S-ITS2, ITS2, COI, and the concatenated SSU-ITS1-5.8S-ITS2-5′LSU rDNA; and rDNA-COI. The best-fit models, GTR + I + G for SSU-V4, LSU-D2 rDNA, ITS1-5.8S-ITS2, ITS2, COI and the concatenated regions, and TIM3 + I + G for SSU rDNA were selected with the Akaike information criterion in jModelTest2 v2.1.6 [48]. For the amino acid alignment of COI, the MtZoa + G evolutionary model was the best-fitted model selected by SMS in PhyML [49,50]. Maximum likelihood (ML) analyses were carried out using PhyML-3.1 [49] with node support from a majority-rule consensus tree of 1000 bootstrap replicates. For the ML bootstraps, we considered values <70% as low, 70–94% as moderate, and ≥95% as high, following [51]. Bayesian inference (BI) was carried out using MrBayes v3.2.6 [52]. Posterior probability was estimated by using four chains running 10,000,000 generations sampling every 1000 generations. The first 25% of sampled trees were considered burn-in trees. For the Bayesian posterior probabilities, we considered values <0.95 as low and ≥0.95 as high, following [53].

2.4. Secondary Structure Prediction

Consensus structures of Pseudokeronopsis ITS2 and LSU-D2 regions (Figures 3A and 5A) were predicted using the LocARNA Sever with default options (available from http://rna.informatik.uni-freiburg.de/LocARNA/Input.jsp) [54]. The consensus secondary structure of ITS2 was compared with the structural template constructed by other spirotrichean ciliates [40], manually corrected and used as a model for Pseudokeronopsis species. With the guidance of the model, putative secondary structures of the Pseudokeronopsis ITS2 molecules were constructed on RNAfold server of the ViennaRNA Web Services (available from http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi) [55] with constraint folding strategy. The folding results for minimum free energy (MFE) structure drawing encoding base-pair probabilities were chosen as outputs.
The model of LSU-D2 secondary structure was not available for any ciliate group. Consequently, the structures predicted using the default values were compared to each other as well as the consensus one to establish the model. The structures of all Pseudokeronopsis populations except P. carnea JQ424836 are nearly identical with the consensus one, which was selected as the LSU-D2 model in Pseudokeronopsis. On the basis of the model, the secondary structure of P. carnea JQ424836 was reconstructed with constraint folding strategy. Once the conserved structural models of ITS2 and LSU-D2 in Pseudokeronopsis ciliates were established, compensating base changes (CBCs, e.g., G = C ↔ C = G, A = U or U = A) were also detected with the CBCAnalyzer [56] as implemented in 4SALE, and hemi-CBCs (compensatory change on lonely one side of a helix pairing, i.e., G = C ↔ G = U) were found manually.

3. Results

3.1. Molecular Sequences, Genetic Distance, and Species Delimitation Analyses

A total of 28 new sequences from six species were deposited in GenBank and their accession numbers are provided in Table 1. The alignments of SSU-V4, ITS2, ITS1-5.8S-ITS2, LSU-D2, and COI had lengths of 226 bp, 199 bp, 494 bp, 247 bp, and 480 bp, respectively. The inter- and intraspecific genetic divergences of these alignments were calculated to investigate the genetic distances in Pseudokeronopsis. In total, seven species: P. rubra, P. cf. songi, P. songi, P. erythrina, P. carnea, P. pararubra, and P. flava were included the initial statistics of genetic distances among conspecifics and congeners. Simultaneously, sequences noted as Pseudokeronopsis sp. and those without associated accurate morphological publications (marked with gray shadows in Table S1–S5) were excluded from the initial statistics. Maximum genetic distances among conspecifics and minimum genetic distances among congeners were used to identify those uncertain species.
For the COI alignment, the genetic distances among the former six species ranged from 14.27% to 22.27%, whereas there was no genetic variation based on the initial statistics of intraspecific distances (Table S1). On the basis of the distance data and the lack of morphological information, P. rubra MG594873 and P. carnea MG594874 + MG594875 (marked as C4 and C5 in the COI trees) were regarded as two uncertain species of Pseudokeronopsis (Table 2). When the data of all known and uncertain species was included, the intraspecific distances of COI were in the range of 0–2.58%, whereas the interspecific ones were in the range of 13.25–22.28% (Table S1; Figure 1).
Like the case of COI, the genetic distances among the seven species were also distinctly higher than those among conspecifics for the ITS1-5.8S-ITS2 (3.68–13.01% vs. 0–0.42%), ITS2 (6.55–16.93% vs. 0–0.52%), and LSU-D2 (2.53–17.51% vs. zero; Table S2–S4). When these cryptic and uncertain species (marked as C1, C2, and C3 in the ITS1-5.8S-ITS2 trees) were included, no change occurred in the range of the intraspecific distances, whereas the interspecific distances were slightly extended to 3.23–13.28% for ITS1-5.8S-ITS2, 3.72–16.93% for ITS2, and 2.1–17.51% for LSU-D2 (Table 2, Table S2–S4; Figure 1).
The SSU-V4 genetic distances among the seven species ranged from zero to 2.75%, corresponding to zero to six variable sites and the sequence similarity of 97.3–100%. No intraspecific variation existed (Table S5). Pseudokeronopsis cf. songi, P. songi, the uncertain species C4, and P. rubra HQ228548 (C2) were even identical in the SSU-V4 sequence.
All available sequences of Pseudokeronopsis in GenBank were assigned to seven species (Table 1; Figure 2). The species delimitation results based on the ITS1-5.8S-ITS2, ITS2, and LSU-D2 were identical, which not only further confirmed the status of the seven species, but also uncovered the two cryptic species C1 and C3 and the uncertain species C2 (Table 2; Figure 2). The COI sequences were only available for six nominal species, and the species delimitation analysis conformed the status of the six species and distinguished the two uncertain species C4 and C5 as well (Table 2; Figure 2). In contrast, only one and three OTUs could be detected from the seven species when using the SSU-V4 sequence similarity of 97% and 99%, respectively, as a threshold (Figure 2). Because C4 and C5 had no common marker sequences comparable with C1–C3, we could not determine whether C4 and C5 belonged to part of C1–C3. Nonetheless, the species delimitation analyses indicated that the available sequences represented at least 10 Pseudokeronopsis species including 7 known species, 2 cryptic, and at least 1 uncertain species.

3.2. ITS2 and LSU-D2 Secondary Structures

Putative secondary structures of ITS2 transcripts are proposed in Figure 3 and Figure 4, possessing the following features: (1) an internal loop bearing two helices, (2) a conserved Helix A with a motif 5′-GAGA versus UCUC-3′ and a U-U mismatch in the terminal loop, and (3) Helix B being the longest one and containing three sub-helices (B-a,b,c). The lengths of the all the parts except for helix B were identical, and Helix A was much shorter than Helix B (32 vs. 114–117 nt; Figure S1). Pseudokeronopsis rubra MH513652, P. pararubra, P. flava, and P. carnea KU663903, DQ503580, DQ503582, and EF174292–7 shared the same structure pattern of Helix A without any CBC or hemi-CBC (Figure 3C). The CBCs, hemi-CBCs, and/or different patterns of the bugles were present within Helix A in other species of Pseudokeronopsis (marked with filled and hollow arrows in Figure 3C). The structures of the sub-helices B-a and c were highly similar among Pseudokeronopsis species, whereas the sub-helix B-b showed much variation in lengths, the positions of the bulge, and the numbers of CBC and hemi-CBA (Figure 4). The CBC analysis revealed a maximum of two CBCs between Pseudokeronopsis species and no CBC between the conspecifics (Table 3).
The consensus secondary structure of LSU-D2 showed a closed central loop with three helices (I, II, and III; Figure 5A). Helix I was the shortest and most conserved region, which showed no variation among Pseudokeronopsis populations (Figure 5A). Helix II contained at least two U-U mismatches and a terminal loop with 5′-AAA versus GGA-3′. Helix III included an A-A/C mismatch and A terminal loop with 5′-AAA versus G-3′, and 5′-UU versus UC (or U) U-3′ (Figure 6). The CBC analysis revealed a maximum of nine CBCs between Pseudokeronopsis species and no CBC between the conspecifics (Table 4).

3.3. Molecular Phylogenetic Analyses

The phylogenetic reconstruction of SSU-V4 rDNA showed similar results as that of the nearly full-length SSU rDNA, and the same was true for ITS2 and ITS1-5.8S-ITS2, but with lower resolution for both of the part regions. Hence, the SSU-V4 and ITS2 topologies are not provided here. The COI trees constructed using the nucleotide and the amino acid sequences were identical in topology. Thus, only the nucleotide tree, which had relatively higher support values, was selected for the analyses. Both the Bayesian inference (BI) and maximum likelihood analyses (ML) yielded a nearly identical topology by using either the concatenated ITS1-5.8S-ITS2-5′LSU or the SSU rDNA datasets. Thus, a single topology was presented for each of these datasets with support values indicated near branches (Figure 7; Figure S2).
In all the phylogenetic trees, the branch of Pseudokeronopsis consisted of two major clades/groups (I and II; Figure 7 and Figure 8; Figures S2–S4). In sum, Clade I contained P. songi, P. cf. songi, P. flava, P. erythrina, Pseudokeronopsis sp. JQ424859, three uncertain/cryptic Pseudokeronopsis species (marked with C1, C2, C4), and Uroleptopsis citrina. Clade/Group II included P. rubra, P. pararubra, and two uncertain/cryptic species (marked with C3, C5). In all the trees, the genus Pseudokeronopsis was not monophyletic due to Uroleptopsis citrina nesting into the Clade I with low to high support.

4. Discussion

4.1. Suitable Barcodes for Species Delimitation of Pseudokeronopsis

On the basis of the revision of [12,16] and the present study, the genus Pseudokeronopsis now contains 12 species, and most of them occur at high frequencies in marine or limnetic habitats worldwide [7,12,16,19,20,58]. However, Pseudokeronopsis is still one of the most difficult genera for species discrimination due to their high morphological similarity [12,16,19]. Misidentifications are frequent and thus there is a pressing need to assess DNA barcodes for the species delimitation [12].
In the present study, both the ITS1-5.8S-ITS2 and ITS2 regions showed a similar degree of genetic variations and distinct gaps between the intra- and interspecific distances (0.52–3.72% for ITS2; 0.42–3.84% for ITS-5.8S-ITS2), suggesting they are promising barcodes for Pseudokeronopsis. Considering that the almost all intraspecific genetic distances of the previously investigated ciliates are no more than 1.5% [14,15,25,26,30,31,32,43,59,60], we identified a sequence divergence of 1.5% as an ideal threshold of ITS1-5.8S-ITS2 and ITS2 to discriminate Pseudokeronopsis species. Like in the case of the ITS region markers, the LSU-D2 showed a distinct genetic distance gap (0–2.1%) as well, indicating it is also a promising barcode. Furthermore, the species delimitation results based on the ITS1-5.8S-ITS2, ITS2, and LSU-D2 regions were identical (Figure 2). Thus, to promote confidence in uncovering novel species, we proposed the multilocus barcoding marker including ITS-5.8S-ITS2 and LSU-D2 as a better choice.
COI has the highest rate of evolution among these markers, and has been frequently tested for DNA barcoding of ciliates [13,14,26,34,61,62,63]. In the present analyses, COI showed a distinct gap (2.58–13.25%) in discriminating all available species, showing the potential for species delimitation of Pseudokeronopsis. However, the exact range of intraspecific genetic distance may not have been fully uncovered because the intraspecific data were obtained from only two populations each of P. cf. songi and the uncertain species C4. It is worth noting that previous studies indicated that the intraspecific divergence range of COI was highly variable, even within the same genus. For instance, the intraspecific genetic distance range was ≤16.4% in Frontonia canadensis and ≤4.8% in F. sinica [14]. Thus, it is risky to infer a threshold of COI to uncover novel species before comprehensive studies with rich sampling are conducted.
To summarize, the ITS1-5.8S-ITS2, ITS2, and LSU-D2 regions are promising barcodes for Pseudokeronopsis, and the multilocus marker including ITS-5.8S-ITS2 and LSU-D2 is more suitable to uncover novel species. The mitochondrial COI is also a candidate barcode, but more intraspecific data are needed to infer the threshold for species discrimination. In contrast, SSU-V4 is useless in discriminating Pseudokeronopsis species due to the absence of a barcoding gap.

4.2. Appropriate Metabarcoding Markers for Ciliate Diversity Estimation

On the basis of second-generation sequencing (SGS) methods, the SSU-V4 marker has been widely used for estimating ciliate environmental diversity, with the sequence similarity of 97% used as a threshold to identify species [8,10,64]. In our test, the interspecific genetic distances of Pseudokeronopsis were in the range of 0–2.75%, corresponding to 0–6 variable sites and 97.3–100% sequence similarities. Thus, only one and three OTUs/species could be detected from the seven known species when using the SSU-V4 sequence similarity of 97% and 99%, respectively, as a threshold (Table 1; Figure 2). Similar results were obtained in Tintinnida and Frontonia [9,14]. As such, using the SSU-V4 sequence similarity of 97% or even 99% as a threshold to delineate species may underestimate the environmental diversity of ciliates. An appropriate marker is thus in pressing need for the estimation of environmental diversity.
The mitochondrial COI showed potential in discriminating Pseudokeronopsis species and several oligohymenophorean taxa, for example, with Tetrahymena, Paramecium, and Frontonia [14,25,34]. However, COI is not a good metabarcode for estimating the ciliate environmental diversity due to several issues. First, the barcoding gaps of COI are highly variable (e.g., about 1–2% in Tetrahymena vs. 16.4–25.2% in Frontonia) [14,34], making it potentially unable to generalize for different taxa. Second, the lack of universal primer sets hampers its utility as a metabarcode. Additionally, there are no mitochondria in some ecologically important taxa [7].
The D1-D2 region (>600 bp) of LSU rDNA has also been proposed as an excellent candidate barcode for ciliates due to the presence of barcoding gaps in Paramecium and Tintinnida [5,9]. Considering the second-generation sequencing (SGS) read limitation (<600 bp), Zhao et al. [14] tested the barcoding potentials of the D1 and D2 regions and indicated that only LSU-D2 showed a distinct barcoding gap (2–4.5%) in discriminating Frontonia species. In the present study, LSU-D2 also showed a barcoding gap (0–2.1%) to separate congeners. However, the overlapping part (2–2.1%) between the two gaps were too tiny to set a uniform threshold to discriminate all ciliates. On the other hand, the universal primers for the LSU-D2 amplification (<600 bp) were not available, hampering its utility as a SGS metabarcode.
The ITS1-5.8S-ITS2 region (<600 bp), by contrast, showed relatively low intraspecific distances (<1.5%) in most ciliate groups investigated, and over half of conspecific sequences were identical even in populations from distant localities (Figure 9) ([14,15,25,26,28,29,30,31,32,33,43,59,60,65,66] and the present study). For instance, no variation in the ITS1-5.8S-ITS2 region was observed for nine populations of Paramecium caudatum from China, Germany, Italy, and Australia. The same results were found in Vorticella convallaria populations from China and the USA, and in 21 strains of Philasterides dicentrarchi isolated from 3 different hosts and 11 localities in Korea and Japan [26,32,59]. When the probable misidentifications mentioned in the above studies were excluded, the interspecific distances in almost all these investigated taxa except for Paramecium and Helicostomella were >3%, showing a general barcoding gap (Figure 9).
So far, no distance-based marker could discriminate all closely related species in Paramecium, and one explanation is that some species complexes, for example, the P. bursaria complex as well as P. multimicronucelatum complex, have been officially described as a single species [5,26,66]. The tintinnid Helicostomella is a difficult genus possessing polymorphism, making it difficult in unambiguously identifying it down to a species on the basis of lorica morphology, and no single marker (<600 bp) alone can discriminate all the closely related species [33,60]. Therefore, considering the barcoding efficacy and suitable amplication length, we proposed the ITS1-5.8S-ITS2, instead of SSU-V4, as a candidate metabarcode for assessing ciliate environmental diversity. The ITS1-5.8S-ITS2 as a DNA metabarcoding marker can also be complemented by the CBC-based and tree-based criterions (see the following Section 4.3 and Section 4.4). Nevertheless, more genus data from different classes including rich congeneric and conspecific sequences are needed to confirm or reject this hypothesis.
It is worth noting that the third-generation sequencing has been used to study prokaryotic and eukaryotic microbial diversity [67,68]. Longer sequencing reads of several thousand base pairs (bp) are now possible. Considering that the combination of ITS1-5.8S-ITS2 and LSU-D1-D2 performed better than each marker in discriminating closely related species of Paramecium and tintinnids [60], the multilocus barcode including ITS1-5.8S-ITS2 and LSU-D2 may provide better taxonomic resolution. The present marker ITS1-5.8S-ITS2-5′LSU (ca. 1300 bp) including the complete ITS-5.8S-ITS2 and LSU-D2 regions is easily amplified in ciliates using the conserved primers ITSF and CILI28S-1000 [40,41], making multilocus barcoding feasible in identifying environmental ciliates with the third-generation sequencing platform. Furthermore, the multilocus metabarcoding marker can also complement the CBC-based and tree-based criteria (see the following Section 4.3 and Section 4.4).

4.3. The Species Delimitation Efficiency of CBCs in Putative Secondary Structures

Compensatory base changes (CBCs) in the ITS2 secondary structure have been found to correlate strongly with distinct biological species, and have already been used as a marker for species delimitation in many eukaryotes [43,56,69,70,71,72,73]. For plants and fungi, at least one CBC is a classifier with a 93% reliability that indicates two organisms belonging to two species [70]. This criterion has been tested for the species identification of three ciliate genera, whereas the results are discordant—it worked well for Paramecium and Vorticella [31,32,43] but was regarded as useless for Spirostomum [30]. Whether the CBC-based marker is efficient for ciliate identification in the same way as the corresponding distance-based barcode requires further evaluation.
In the case of Pseudokeronopsis, no CBC was observed between the conspecifics and between a few closely related species, and CBCs only occur between different species (Table 3; Figure 3C and Figure 4). Consequently, the presence of CBCs can be a classifier with a 100% reliability for Pseudokeronopsis to assign different populations to separate species. However, the suggested threshold is unidirectional, that is, the lack of CBCs does not necessarily indicate two organisms belonging to the same species, for instance, no CBC occurs between P. rubra and P. pararubra (Table 3; Figure 3C and Figure 4). This was also shown in Spirostomum, Paramecium, and Vorticella, in which CBCs only exist between distant species, whereas no CBC was observed among some closely related species, for example, Spirostomum minus, S. ambiguum, and S. subtilis [30]. All these test cases indicate that the CBC presence in the ITS2 secondary structure can serve as a good indicator of species separation in ciliates, but not vice versa.
For those closely related species without CBCs, there are 1–3 hemi-CBCs of ITS2 existing among congeners, whereas 0–1 hemi-CBC occurs among conspecifics (Figure 5C and Figure 6, hollow arrows). Since the taxonomic significance of the presence of hemi-CBCs is still not resolved, more data are needed to test whether a number of ITS2 hemi-CBCs can be served as an indicator of species discrimination within Pseudokeronopsis.
This study is the first attempt to construct a LSU-D2 secondary structure model for a ciliate genus. Among the helices of the model, Helixes II and III contained variant positions to provide cases of CBC and hemi-CBC, the key to “proof” of secondary structure pattern (marked as filled and hollow arrows in Figure 5 and Figure 6). Like the case of ITS2, the presence of one CBC of LSU-D2 can be used separate two organisms into different species with a 100% reliability, whereas the absence of CBCs does not indicate two organisms belonging to the same species. However, due to the conspecific secondary structures being only available from three species of Pseudokeronopsis, more conspecific data are needed to confirm this criterion.
Due to the unidirection of the thresholds, the ITS2 and LSU-D2 CBC-based markers showed lower efficiencies than those of the distance-based barcodes in species delimitation of ciliates. Nevertheless, the CBC-based markers can provide further evidence to separate some closely related species and to uncover cryptic species, for example, one ITS2 CBC occurs between the songi-like cryptic species C1 and P. songi, and at least one LSU-D2 CBC exits between the pararubra-like cryptic species C3 and congeners (Table 3 and Table 4; Figure 3C, Figure 5B and Figure 6). Therefore, the CBC-based marker is useful to complement the traditional DNA barcoding by giving second structure evidence.

4.4. Phylogenetic Relationships within Pseudokeronopsis

The non-monophyly of Pseudokeronopsis was revealed by the rDNA data because Uroleptopsis citrina clustered into this group in the SSU, LSU, and ITS-5.8S rDNA trees [24,27,74]. By contrast, in the COI trees (including U. citrina) of [27], the Pseudokeronopsis species formed a monophyletic clade but with low support. In the present study, Uroleptopsis citrina nested into the branch of Pseudokeronopsis in both the rDNA and the COI trees with moderate to high support (Figure 7 and Figure 8; Figures S2–S4). However, Uroleptopsis and Pseudokeronopsis can be easily separated via morphological features, for example, the transverse cirri (absent vs. present) and the position of the buccal cirri (in a pattern of bicorona vs. on the right of paroral) [16,24]. The question whether the morphological differences justify the separation of U. citrina and Pseudokeronopsis species at genus level has been discussed by [16]. No new sequence of Uroleptopsis is available, thus we refrain from this molecular biologically-indicated issue and focus only on the phylogeny of Pseudokeronopsis.
The topologies constructed with the selected markers showed more or less differences (Figure 7 and Figure 8; Figures S2–S5). To test which molecule is suitable for elucidating the relationships within Pseudokeronopsis properly, we compared the molecular topologies with the morphological data. The SSU rDNA trees formed many parallel clades, which are useless in clarifying the relationships among P. erythrina, P. songi, P. cf. songi, and P. flava (Figure S2). Likewise, LSU-D2 provides few solutions for this group due to the low support (BI posterior probabilities ≤0.77) for two-thirds of the Pseudokeronopsis branches (Figure S3). The ML topology of ITS1-5.8S-ITS2 was also discordant with the morphological data, as P. pararubra formed a sister clade with P. carnea instead of P. rubra (Figure S5). The present COI topologies conflicted with the morphological hypothesis. For example, among the known congeners, P. songi is morphologically most similar to P. cf. songi, but they did not form a sister clade in the COI trees; P. rubra, which morphologically resembles P. pararubra, fell outside of the clade including P. pararubra and P. carnea (Figure 8).
The BI tree of ITS1-5.8S-ITS2, the concatenated ITS1-5.8S-ITS2-5′LSU, and ITS1-5.8S-ITS2-5′LSU-COI trees showed nearly identical topology (Figure 7 and Figure S3), which corresponded well with the morphological data. These markers can provide tree-based evidence to complement the traditional DNA barcoding for ciliate identification. Among these trees, the ITS1-5.8S-ITS2-5′LSU region provided relatively higher support values (nearly all BI posterior probabilities >0.95 and ML bootstraps >70%). Thus, the ITS1-5.8S-ITS2-5′LSU trees were used for elucidating the relationships among Pseudokeronopsis species. Within the Pseudokeronopsis branch of the ITS1-5.8S-ITS2-5′LSU trees, P. carnea branched early, and P. rubra and the group P. pararubra + the C3 cryptic species formed an intermediate clade (Clade II) between P. carnea and the rest of the Pseudokeronopsis species (Clade I). Within Clade I, P. songi and the C1 cryptic species formed a terminal clade, followed by P. cf. songi, the uncertain species C2, P. erythrina, and P. flava, successively. Although the described species P. multinucleata, P. decolor, P. elongate, and P. similis have still not been barcoded, their relationships with sequenced taxa can be inferred on the basis of their morphology. These four species are highly similar in ciliary pattern, and close to the Clade II species according to the characters of the midventral complex connected to six or more transverse cirri. On the basis of the position of contractile vacuole, these four species may form a terminal subclade within Clade II (anterior 1/3 or mid-body vs. posterior 1/3–2/5 for the rest species). Within this subclade, P. decolor resembles P. multinucleata in the position of contractile vacuole (at mid-body), and P. elongata is close to P. similis due to the distinct moniliform macronucleus and freshwater habitat. Additionally, compared to congeners, the midventral cirral rows of P. multinucleata are distinctly separated, also suggesting it nests deeply within Clade II. As such, we combined all available morphological and molecular data and drew a putative evolutionary tree for the genus Pseudokeronopsis (Figure 10).

4.5. Confirmation and Revision of Pseudokeronopsis Sequences in GenBank

In the nucleotide database of NCBI (https://www.ncbi.nlm.nih.gov/nuccore/?term=txid311385[Organism:exp] on 22 July 2019), there are 140 sequences which were assigned to six known species and four uncertain species of Pseudokeronopsis. On the basis of the morphological information and the ITS1-5.8S-ITS2 genetic gap, Li et al. [9] considered that over half of these sequences were misidentified or dubious. However, with the new availability of the ITS1-5.8S-ITS2 sequences of the type species P. rubra and P. cf. songi, we found that the intraspecific genetic distance range sensu [9] is overestimated (0–4.41% vs. the present 0–0.52%). Thus, we confirmed the reliable sequences and revised the misidentified sequences again on the basis of the present genetic distance, secondary structure, and phylogenetic tree data, as well as the morphological information.
The present genetic distances and the phylogenetic analyses confirmed the validity of P. carnea, P. songi, P. pararubra, P. erythrina, and P. flava sequences, which were submitted with morphological data [9,25,27,31]. The newly sequenced P. cf. songi might represent a new species, which is indicated by the genetic distances of the ITS region (>1.5% between it and congeners). On the basis of the sequence characters of P. rubra combined with morphological analysis (not shown here), we confirmed that all previously submitted sequences of P. rubra were derived from misidentified populations (Table 2).
Without morphological and barcoding data, the identification of the Pseudokeronopsis species sensu [30,54] are in progress (Table 2). Additionally, P. erythrina was originally described as an estuarine species, whereas the SSU rDNA sequences P. erythrina KX459375 and MG994990 were obtained from a freshwater environment. Thus, other barcodes are needed to confirm their validity (Table 2) [22], Kaur et al., unpublished. In sum, a total of 25 misidentified or uncertain Pseudokeronopsis sequences were clarified and proposed to be amended, and 65 sequences are still awaiting to be confirmed (Table 2).

5. Conclusions

In the present analyses, all the test markers except SSU-V4 were useful in discriminating Pseudokeronopsis species. SSU-V4 rDNA was too conservative for the species discrimination, and only one and three OTUs could be obtained from seven Pseudokeronopsis species by using a sequence similarity of 97% and 99%, respectively, as a threshold. The nuclear markers ITS1-5.8S-ITS2, ITS2, and LSU-D2 are promising barcodes, and the multilocus marker including ITS-5.8S-ITS2 and LSU-D2 was more suitable in uncovering novel species. The mitochondrial COI was also a candidate barcode, but more intraspecific data are needed to confirm it. The presence of CBCs was an indicator for Pseudokeronopsis to assign different populations to separate species, but not vice versa. The ITS1-5.8S-ITS2-5′LSU region provided better phylogenetic resolution of Pseudokeronopsis at species level. At least 25 Pseudokeronopsis sequences are misidentified in GenBank and are proposed to be amended. Considering the overlapping barcoding gap between different ciliate taxa, universal primers, suitable amplification lengths, and regularly updated databases available, we propose the multilocus marker including the ITS1-5.8S-ITS2 and LSU-D2 as an ideal candidate barcode, and the ITS1-5.8S-ITS2 as a candidate SGS metabarcode for assessing ciliate environmental diversity. Furthermore, the CBC-based and tree-based criteria are useful to complement the DNA metabarcoding method by giving second structure and phylogenetic evidences. With the application of third-generation sequencing, the multilocus metabarcode may provide greater taxonomic resolution.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2607/7/11/493/s1: Figure S1: Alignment of ITS2 showing the lengths of the helices. Figure S2: Bayesian inference (BI) and maximum likelihood (ML) trees based on 18S rDNA sequences showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent BI and ML support values. Newly sequenced species are in bold. *, misidentified. Figure S3: Maximum likelihood (ML) and Bayesian inference (BI) trees based on the concatenated ITS1-5.8S-ITS2-5′LSU-COI (ITS-COI) sequences showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent ML and BI support values. Newly sequenced species are in bold. *, misidentified. Figure S4: Bayesian inference (BI) tree based on LSU-D2 sequences showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent BI support values. Newly sequenced species are in bold. *, misidentified. **, uncertain species. Figure S5: Maximum likelihood (ML) tree based on ITS1-5.8S-ITS2 sequences showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent BI support values. Newly sequenced species are in bold. *, misidentified. **, uncertain species. Table S1: Genetic distances of COI between Pseudokeronopsis species/populations. Table S2: Genetic distances of ITS1-5.8S-ITS2 between Pseudokeronopsis species/populations. Table S3: Genetic distances of ITS2 between Pseudokeronopsis species/populations. Table S4 Genetic distances of LSU-D2 between Pseudokeronopsis species/populations. Table S5 Genetic distances of SSU-V4 between Pseudokeronopsis species/populations.

Author Contributions

Conceptualization, Z.Z. and K.X.; methodology, Z.Z. and J.L.; validation, Z.Z., J.L. and K.X.; formal analysis, Z.Z. and J.L.; writing—original draft preparation, Z.Z. and J.L.; writing—review and editing, Z.Z and K.X.; visualization, Z.Z.; supervision, K.X.; funding acquisition, K.X.

Funding

This research was funded by National Natural Science Foundation of China, grant numbers 41876171, 41476144 and 41406171. The APC was funded by 41876171.

Acknowledgments

Many thanks to the anonymous reviewers for comments and discussion that improved the quality of this manuscript.

Conflicts of Interest

We declare that we have no competing interests.

References

  1. April, J.; Mayden, R.L.; Hanner, R.H.; Bernatchez, L. Genetic calibration of species diversity among North America’s freshwater fishes. Proc. Natl. Acad. Sci. USA 2011, 108, 10602–10607. [Google Scholar] [CrossRef]
  2. Hebert, P.D.N.; Ratnasingham, S.; de Waard, J.R. Barcoding animal life: Cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lond. B Biol. 2003, 270, S96–S99. [Google Scholar] [CrossRef]
  3. Hebert, P.D.N.; Stoeckle, M.Y.; Zemlak, T.S.; Francis, C.M. Identification of birds through DNA barcodes. PLoS Biol. 2004, 2, e312. [Google Scholar] [CrossRef] [PubMed]
  4. Pawlowski, J.; Audic, S.; Adl, S.; Bass, D.; Belbahri, L.; Berney, C.; Bowser, S.S.; Cepicka, I.; Decelle, J.; Dunthorn, M.; et al. CBOL Protist Working Group: Barcoding eukaryotic richness beyond the animal, plant, and fungal kingdoms. PLoS Biol. 2012, 10, e1001419. [Google Scholar] [CrossRef]
  5. Stoeck, T.; Przybos, E.; Dunthorn, M. The D1-D2 region of the large subunit ribosomal DNA as barcode for ciliates. Mol. Ecol. Resour. 2014, 14, 458–468. [Google Scholar] [CrossRef]
  6. Elwood, H.J.; Olsen, G.J.; Sogin, M.L. The small–subunit ribosomal RNA gene sequences from the hypotrichous ciliates Oxytricha nova and Stylonychia pustulata. Mol. Biol. Evol. 1985, 2, 399–410. [Google Scholar] [PubMed]
  7. Lynn, H.D. The Ciliated Protozoa: Characterization, Classification and Guide to the Literature; Springer: Dordrecht, The Netherlands, 2008. [Google Scholar]
  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–187. [Google Scholar] [CrossRef]
  9. 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]
  10. Zhao, F.; Filker, S.; Stoeck, T.; Xu, K. Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean. BMC Microbiol. 2017, 17, 192. [Google Scholar] [CrossRef] [PubMed]
  11. Lara, E.; Acosta-Mercado, D. A molecular perspective on ciliates as soil bioindicators. Eur. J. Soil Bol. 2012, 49, 107–111. [Google Scholar] [CrossRef]
  12. Li, J.; Zhan, Z.; Xu, K. Systematics and molecular phylogeny of the ciliate genus Pseudokeronopsis (Ciliophora, Hypotrichia). J. Eukaryot. Microbiol. 2017, 64, 850–872. [Google Scholar] [CrossRef] [PubMed]
  13. Strüder-Kypke, M.C.; Lynn, D.H. Comparative analysis of the mitochondrial cytochrome c oxidase subunit I (COI) gene in ciliates (Alveolata, Ciliophora) and evaluation of its suitability as a biodiversity marker. Syst. Biodivers. 2010, 8, 131–148. [Google Scholar] [CrossRef]
  14. Zhao, Y.; Yi, Z.; Gentekaki, E.; Zhan, A.; Al-Farraj, S.A.; Song, W. Utility of combining morphological characters, nuclear and mitochondrial genes: An attempt to resolve the conflicts of species identification for ciliated protists. Mol. Phylogenet. Evol. 2016, 94, 718–729. [Google Scholar] [CrossRef]
  15. Zhao, Y.; Yi, Z.; Warren, A.; Song, W. Species delimitation for the molecular taxonomy and ecology of a widely distributed microbial eukaryotes genus Euplotes (Alveolata, Ciliophora). Proc. Biol. Sci. 2018, 285, 20172159. [Google Scholar] [CrossRef] [PubMed]
  16. Berger, H. Monograph of the Urostyloidea (Ciliophora, Hypotricha). Monogr. Biol. 2006, 85, 731–980. [Google Scholar]
  17. Borror, A.C.; Wicklow, B.J. The suborder Urostylina, Jankowski (Ciliophora, Hypotrichida): Morphology, systematics and identification of species. Acta Protozool. 1983, 22, 97–126. [Google Scholar]
  18. Chen, X.; Clamp, J.C.; Song, W. Phylogeny and systematic revision of the family Pseudokeronopsidae (Protista, Ciliophora, Hypotricha), with description of a new estuarine species of Pseudokeronopsis. Zool. Scr. 2011, 40, 659–671. [Google Scholar] [CrossRef]
  19. Song, W.; Warren, A.; Roberts, D.; Wilbert, N.; Li, L.; Sun, P.; Hu, X.; Ma, H. Comparison and redefinition of four marine, colored Pseudokeronopsis spp. (Ciliophora, Hypotrichida), with emphasis on their living morphology. Acta Protozool. 2006, 45, 271–287. [Google Scholar]
  20. Yi, Z.; Chen, Z.; Warren, A.; Roberts, D.M.; Al-Rasheid, K.A.S.; Miao, M.; Gao, S.; Shao, C.; Song, W. Molecular phylogeny of Pseudokeronopsis (Protozoa, Ciliophora, Urostylida), with reconsideration of three closely related species at inter- and intra-specific levels inferred from the small subunit ribosomal RNA gene and the ITS1-5.8S-ITS2 region sequences. J. Zool. 2008, 275, 268–275. [Google Scholar] [CrossRef]
  21. Anesi, A.; Buonanno, F.; Di Giuseppe, G.; Ortenzi, C.; Guella, G. Metabolites from the euryhaline ciliate Pseudokeronopsis erythrina. Eur. J. Org. Chem. 2016, 7, 1330–1336. [Google Scholar] [CrossRef]
  22. Buonanno, F.; Anesi, A.; Giuseppe, G.D.; Guella, G.; Ortenzi, C. Chemical Defense by Erythrolactones in the Euryhaline Ciliated Protist, Pseudokeronopsis erythrina. Zool. Sci. 2017, 34, 42–51. [Google Scholar] [CrossRef] [PubMed]
  23. Yi, Z.; Huang, L.; Yang, R.; Lin, X.; Song, W. Actin evolution in ciliates (Protist, Alveolata) is characterized by high diversity and three duplication events. Mol. Biol. Evol. 2016, 96, 45–54. [Google Scholar] [CrossRef] [PubMed]
  24. Huang, J.; Chen, Z.; Song, W.; Berger, H. Three-gene based phylogeny of the Urostyloidea (Protista, Ciliophora, Hypotricha), with notes on classification of some core taxa. Mol. Phylogenet. Evol. 2014, 70, 337–347. [Google Scholar] [CrossRef] [PubMed]
  25. Barth, D.; Krenek, S.; Fokin, S.I.; Berendonk, T.U. Intraspecific genetic variation in Paramecium revealed by mitochondrial cytochrome c oxidase I sequences. J. Eukaryot. Microbiol. 2006, 53, 20–25. [Google Scholar] [CrossRef]
  26. Greczek-Stachura, M.; Potekhin, A.; Przyboś, E.; Rautian, M.; Skoblo, I.; Tarcz, S. Identification of Paramecium bursaria syngens through molecular markers–comparative analysis of three loci in the nuclear and mitochondrial DNA. Protist 2012, 163, 671–685. [Google Scholar] [CrossRef]
  27. Park, M.H.; Jung, J.H.; Jo, E.; Park, K.M.; Baek, Y.S.; Kim, S.J.; Min, G.S. Utility of mitochondrial COI sequences for species discrimination of Spirotrichea ciliates (Protozoa, Ciliophora). Mitochondrial DNA A DNA Mapp. Seq. Anal. 2019, 30, 148–155. [Google Scholar]
  28. Tarcz, S.; Przyboś, E.; Surmacz, M. An assessment of haplotype variation in ribosomal and mitochondrial DNA fragments suggests incomplete lineage sorting in some species of the Paramecium aurelia complex (Ciliophora, Protozoa). Mol. Phylogenet. Evol. 2013, 67, 255–265. [Google Scholar] [CrossRef]
  29. Tarcz, S.; Rautian, M.; Potekhin, A.; Sawka, N.; Beliavskava, A.; Kiselev, I.; Nekrasova, I.; Przyboś, E. Paramecium putrinum (Ciliophora, Protozoa): The first insight into the variation of two DNA fragments–molecular support for the existence of cryptic species. Mol. Phylogenet. Evol. 2014, 73, 140–145. [Google Scholar] [CrossRef]
  30. Shazib, S.U.A.; Vďačný, P.; Kim, J.H.; Jang, S.W.; Shin, M.K. Molecular phylogeny and species delimitation within the ciliate genus Spirostomum (Ciliophora, Postciliodesmatophora, Heterotrichea), using the internal transcribed spacer region. Mol. Biol. Evol. 2016, 102, 128–144. [Google Scholar] [CrossRef]
  31. Sun, P.; Clamp, J.C.; Xu, D. Analysis of the secondary structure of ITS transcripts in peritrich ciliates (Ciliophora, Oligohymenophorea): Implications for structural evolution and phylogenetic reconstruction. Mol. Phylogenet. Evol. 2010, 56, 242–251. [Google Scholar] [CrossRef]
  32. Sun, P.; Clamp, J.C.; Xu, D.; Huang, B.; Shin, M.K.; Turner, F. An ITS–based phylogenetic framework for the genus Vorticella: Finding the molecular and morphological gaps in a taxonomically difficult group. Proc. Biol. Sci. 2013, 280, 20131177. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, D.; Sun, P.; Shin, M.K.; Kim, Y.O. Species boundaries in tintinnid ciliates: A case study—Morphometric variability, molecular characterization, and temporal distribution of Helicostomella species (Ciliophora, Tintinnina). J. Eukaryot. Microbiol. 2012, 59, 351–358. [Google Scholar] [CrossRef] [PubMed]
  34. Chantangsi, C.; Lynn, D.H.; Brandl, M.T.; Cole, J.C.; Hetrick, N.; Ikonomi, P. Barcoding ciliates: A comprehensive study of 75 isolates of the genus Tetrahymena. Int. J. Syst. Evol. Microbiol. 2007, 57, 2412–2425. [Google Scholar] [CrossRef] [PubMed]
  35. Baek, Y.; Jung, J.; Min, G. Redescription of two marine ciliates (Ciliophora: Urostylida: Pseudstokesokeronopsidae), Pseudokeronopsis carnea and Uroleptopsis citrina, from Korea. Korean J. Syst. Zool. 2011, 27, 220–227. [Google Scholar] [CrossRef]
  36. Li, Y.; Niu, Y.; Liu, L. Phylogenetic studies of four species of ciliate inferred from 16S-like small subunit rRNA gene sequences. J. For. Res. 2008, 19, 119–124. [Google Scholar] [CrossRef]
  37. Smith, D.; Leary, P.; Craggs, J.; Bythell, J.; Sweet, M. Healthy and White Syndrome-Affected Echinopora lamellosa in Aquaria and Experimental Treatment with the Antibiotic Ampicillin. PLoS ONE 2015, 10, e0121780. [Google Scholar] [CrossRef]
  38. Foissner, W. An update of ‘basic light and scanning electron microscopic methods for taxonomic studies of ciliated protozoa’. Int. J. Syst. Evol. Microbiol. 2014, 64, 271–292. [Google Scholar] [CrossRef]
  39. Medlin, L.; Elwood, H.J.; Stickel, S.; Sogin, M.L. The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene 1988, 71, 491–499. [Google Scholar] [CrossRef] [Green Version]
  40. Goggin, C.L.; Murphy, N.E. Conservation of sequence in the internal transcribed spacers and 5.8S ribosomal RNA among geographically separated isolates of parasitic scuticociliates (Ciliophora, Orchitophryidae). Dis. Aquat. Org. 2000, 40, 79–83. [Google Scholar] [CrossRef]
  41. Gong, J.; Kim, S.J.; Kim, S.Y.; Min, G.S.; Roberts, D.M.; Warren, A.; Choi, J.K. Taxonomic redescriptions of two ciliates, Protogastrostyla pulchra ng, n. comb. and Hemigastrostyla enigmatica (Ciliophora: Spirotrichea, Stichotrichia), with phylogenetic analyses based on 18S and 28S rRNA gene sequences. J. Eukaryot. Microbiol. 2007, 54, 468–478. [Google Scholar]
  42. Katoh, K.; Standley, D. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  43. Coleman, A.W. Paramecium aurelia revisited. J. Eukaryot. Microbiol. 2005, 52, 68–77. [Google Scholar] [CrossRef] [PubMed]
  44. Wuyts, J.; De Rijk, P.; Van de Peer, Y.; Winkelmans, T.; De Wachter, R. The European large subunit ribosomal RNA database. Nucleic Acids Res. 2001, 29, 175–177. [Google Scholar] [CrossRef] [PubMed]
  45. Wuyts, J.; Van de Peer, Y.; Winkelmans, T.; De Wachter, R. The European database on small subunit ribosomal RNA. Nucleic Acids Res. 2002, 30, 183–185. [Google Scholar] [CrossRef]
  46. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef]
  47. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  48. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef]
  49. Guindon, S.; Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef]
  50. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [Green Version]
  51. Hillis, D.M.; Bull, J.J. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst. Biol. 1993, 42, 182–192. [Google Scholar] [CrossRef]
  52. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed]
  53. Alfaro, M.E.; Zoller, S.; Lutzoni, F. Bayes or Bootstrap? A simulation study comparing the performance of Bayesian Markov chain Monte Carlo sampling and bootstrapping in assessing phylogenetic confidence. Mol. Biol. Evol. 2003, 20, 255–266. [Google Scholar] [CrossRef] [PubMed]
  54. Will, S.; Joshi, T.; Hofacker, I.L.; Stadler, P.F.; Backofen, R. LocARNA-P: Accurate boundary prediction and improved detection of structural RNAs. RNA 2012, 5, 900–914. [Google Scholar] [CrossRef] [PubMed]
  55. Gruber, A.R.; Lorenz, R.; Bernhart, S.H.; Neuböck, R.; Hofacker, I.L. The Vienna RNA websuite. Nucleic Acids Res. 2008, 36, W70–W74. [Google Scholar] [CrossRef]
  56. Wolf, M.; Friedrich, J.; Dandekar, T.; Müller, T. CBCAnalyzer: Inferring phylogenies based on compensatory base changes in RNA secondary structures. Silico Biol. 2005, 5, 291–294. [Google Scholar]
  57. Yi, Z.; Song, W. Evolution of the order Urostylida (Protozoa, Ciliophora): New hypotheses based on multi-gene information and identification of localized incongruence. PLoS ONE 2011, 6, e17471. [Google Scholar] [CrossRef]
  58. Song, W.; Sun, P.; Ji, D. Redefinition of the yellow hypotrichous ciliate, Pseudokeronopsis flava (Hypotrichida: Ciliophora). J. Mar. Biol. Assoc. U. K. 2004, 84, 1137–1142. [Google Scholar] [CrossRef]
  59. Jung, S.J.; Bae, M.J.; Oh, M.J.; Lee, J. Sequence conservation in the internal transcribed spacers and 5.8S ribosomal RNA of parasitic scuticociliates Miamiensis avidus (Ciliophora, Scuticociliatia). Parasitol. Int. 2011, 60, 216–219. [Google Scholar] [CrossRef]
  60. 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]
  61. Gentekaki, E.; Lynn, D.H. High-level genetic diversity but no population structure inferred from nuclear and mitochondrial markers of the peritrichous ciliate Carchesium polypinum in the Grand River Basin (North America). Appl. Environ. Microbiol. 2009, 75, 3187–3195. [Google Scholar] [CrossRef]
  62. Kher, C.P.; Doerder, F.P.; Cooper, J.; Ikonomi, P.; Achilles-Day, U.; Küpper, F.C.; Lynn, D.H. Barcoding Tetrahymena: Discriminating species and identifying unknowns using the cytochrome c oxidase subunit I (cox-1) barcode. Protist 2011, 162, 2–13. [Google Scholar] [CrossRef] [PubMed]
  63. Lynn, D.H.; Strüder-Kypke, M. Species of Tetrahymena identical by small subunit rRNA gene sequences are discriminated by mitochondrial cytochrome c oxidase I gene sequences. J. Eukaryot. Microbiol. 2006, 53, 385–387. [Google Scholar] [CrossRef] [PubMed]
  64. Stoeck, T.; Breiner, H.W.; Filker, S.; Ostermaier, V.; Kammerlander, B.; Sonntag, B. A morphogenetic survey on ciliate plankton from a mountain lake pinpoints the necessity of lineage-specific barcode markers in microbial ecology. Environ. Microbiol. 2014, 16, 430–444. [Google Scholar] [CrossRef] [PubMed]
  65. Tarcz, S. Intraspecific differentiation of Paramecium novaurelia strains (Ciliophora, Protozoa) inferred from phylogenetic analysis of ribosomal and mitochondrial DNA variation. Eur. J. Protistol. 2013, 49, 50–61. [Google Scholar] [CrossRef] [PubMed]
  66. Tarcz, S.; Potekhin, A.; Rautian, M.; Przyboś, E. Variation in ribosomal and mitochondrial DNA sequences demonstrates the existence of intraspecific groups in Paramecium multimicronucleatum (Ciliophora, Oligohymenophorea). Mol. Phylogenet. Evol. 2012, 63, 500–509. [Google Scholar] [CrossRef] [PubMed]
  67. Heeger, F.; Bourne, E.C.; Baschien, C.; Yurkov, A.; Bunk, B.; Spröer, C.; Overmann, J.; Mazzoni, C.J.; Monaghan, M.T. Long-read DNA metabarcoding of ribosomal RNA in the analysis of fungi from aquatic environments. Mol. Ecol. Resour. 2018, 18, 1500–1514. [Google Scholar] [CrossRef]
  68. Wagner, J.; Coupland, P.; Browne, H.P.; Lawley, T.D.; Francis, S.C.; Parkhill, J. Evaluation of PacBio sequencing for full-length bacterial 16S rRNA gene classification. BMC Microbiol. 2016, 16, 274. [Google Scholar] [CrossRef]
  69. Coleman, A.W. Is there a molecular key to the level of “biological species” in eukaryotes? A DNA guide. Mol. Phylogenet. Evol. 2009, 50, 197–203. [Google Scholar] [CrossRef]
  70. Müller, T.; Philippi, N.; Dandekar, T.; Schultz, J.; Wolf, M. Distinguishing species. RNA 2007, 13, 1469–1472. [Google Scholar] [CrossRef] [Green Version]
  71. Ruhl, M.W.; Wolf, M.; Jenkins, T.M. Compensatory base changes illuminate morphologically difficult taxonomy. Mol. Phylogenet. Evol. 2010, 54, 664–669. [Google Scholar] [CrossRef]
  72. Sorhannus, U.; Ortiz, J.D.; Wolf, M.; Fox, M.G. Microevolution and speciation in Thalassiosira weissflogii (Bacillariophyta). Protist 2010, 161, 237–249. [Google Scholar] [CrossRef] [PubMed]
  73. Wolf, M.; Chen, S.; Song, J.; Ankenbrand, M.; Müller, T. Compensatory base changes in ITS2 secondary structures correlate with the biological species concept despite intragenomic variability in ITS2 sequences—A proof of concept. PLoS ONE 2013, 8, e6672. [Google Scholar] [CrossRef] [PubMed]
  74. Huang, J.; Yi, Z.; Al-Farraj, S.A.; Song, W. Phylogenetic positions and taxonomic assignments of the systematically controversial genera, Spirotrachelostyla, Uroleptopsis and Tunicothrix (Protozoa, Ciliophora, Stichotrichia) based on small subunit rRNA gene sequences. Syst. Biol. 2010, 8, 409–416. [Google Scholar]
Microorganisms 07 00493 g001 550
Figure 1. Genetic distance ranges of ITS1-5.8S-ITS2, ITS2, LSU-D2, SSU-V4, and COI within Pseudokeronopsis.
Figure 1. Genetic distance ranges of ITS1-5.8S-ITS2, ITS2, LSU-D2, SSU-V4, and COI within Pseudokeronopsis.
Microorganisms 07 00493 g001
Microorganisms 07 00493 g002 550
Figure 2. Horizontal bars show the species delimitation results inferred from the morphology and molecular markers (listed on the upper side) for Pseudokeronopsis species. Each color represents a species (online version in color). C1–C5 correspond to the C1–5 taxa marked in the following phylogenetic trees. Dashed line indicates no available data. *, misidentified species.
Figure 2. Horizontal bars show the species delimitation results inferred from the morphology and molecular markers (listed on the upper side) for Pseudokeronopsis species. Each color represents a species (online version in color). C1–C5 correspond to the C1–5 taxa marked in the following phylogenetic trees. Dashed line indicates no available data. *, misidentified species.
Microorganisms 07 00493 g002
Microorganisms 07 00493 g003 550
Figure 3. Secondary structure features of the Pseudokeronopsis ITS2 molecule. (A) Consensus secondary structure predicted using the LocARNA Sever. Color bars (0–1) indicate base-pairing probability. (B) Structural model for Pseudokeronopsis. (C) Helix A structures of congeners. Using the structure of P. carnea as a reference, compensating base changes (CBCs) are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable region is marked by a grey box. Newly sequenced species are in bold. *, misidentified.
Figure 3. Secondary structure features of the Pseudokeronopsis ITS2 molecule. (A) Consensus secondary structure predicted using the LocARNA Sever. Color bars (0–1) indicate base-pairing probability. (B) Structural model for Pseudokeronopsis. (C) Helix A structures of congeners. Using the structure of P. carnea as a reference, compensating base changes (CBCs) are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable region is marked by a grey box. Newly sequenced species are in bold. *, misidentified.
Microorganisms 07 00493 g003
Microorganisms 07 00493 g004 550
Figure 4. Helix B secondary structures of the Pseudokeronopsis ITS2 molecule. Using the structure of P. carnea as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. Newly sequenced species are in bold. *, misidentified.
Figure 4. Helix B secondary structures of the Pseudokeronopsis ITS2 molecule. Using the structure of P. carnea as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. Newly sequenced species are in bold. *, misidentified.
Microorganisms 07 00493 g004
Microorganisms 07 00493 g005 550
Figure 5. Secondary structure features of the Pseudokeronopsis LSU-D2 molecule. (A) Consensus secondary structure predicted using the LocARNA Sever. Color bars (0–1) indicate base-pairing probability. (B) Helix II structures of congeners. Using the structure of P. songi as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable regions are marked by grey boxes. Newly sequenced species are in bold. *, misidentified species.
Figure 5. Secondary structure features of the Pseudokeronopsis LSU-D2 molecule. (A) Consensus secondary structure predicted using the LocARNA Sever. Color bars (0–1) indicate base-pairing probability. (B) Helix II structures of congeners. Using the structure of P. songi as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable regions are marked by grey boxes. Newly sequenced species are in bold. *, misidentified species.
Microorganisms 07 00493 g005
Microorganisms 07 00493 g006 550
Figure 6. Helix III secondary structures of the Pseudokeronopsis LSU-D2 molecule. Using the structure of P. songi as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable regions are marked by grey boxes. Newly sequenced species are in bold. *, misidentified species.
Figure 6. Helix III secondary structures of the Pseudokeronopsis LSU-D2 molecule. Using the structure of P. songi as a reference, CBCs are marked with filled arrows and hemi-CBCs are indicated with hollow arrows. The variable regions are marked by grey boxes. Newly sequenced species are in bold. *, misidentified species.
Microorganisms 07 00493 g006
Microorganisms 07 00493 g007 550
Figure 7. Bayesian inference (BI) trees constructed by ITS1-5.8S-ITS2 (ITS-5.8S) and ITS1-5.8S-ITS2-5′LSU sequences showing phylogenetic relationships among the available pseudokeronopsids. The maximum likelihood (ML) tree and the BI tree of ITS1-5.8S-ITS2-5′LSU are identical in topology, and numbers at the nodes represent BI and ML support values. Newly sequenced species are in bold. *, misidentified species.
Figure 7. Bayesian inference (BI) trees constructed by ITS1-5.8S-ITS2 (ITS-5.8S) and ITS1-5.8S-ITS2-5′LSU sequences showing phylogenetic relationships among the available pseudokeronopsids. The maximum likelihood (ML) tree and the BI tree of ITS1-5.8S-ITS2-5′LSU are identical in topology, and numbers at the nodes represent BI and ML support values. Newly sequenced species are in bold. *, misidentified species.
Microorganisms 07 00493 g007
Microorganisms 07 00493 g008 550
Figure 8. Bayesian inference (BI) and maximum likelihood (ML) trees based on the COI gene showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent BI and ML support values. Newly sequenced species are in bold. *, misidentified species.
Figure 8. Bayesian inference (BI) and maximum likelihood (ML) trees based on the COI gene showing phylogenetic relationships among the available pseudokeronopsids. Numbers at the nodes represent BI and ML support values. Newly sequenced species are in bold. *, misidentified species.
Microorganisms 07 00493 g008
Microorganisms 07 00493 g009 550
Figure 9. Genetic distances of ITS1-5.8S-ITS2 from the previous and present studies. Data were obtained from or calculated on the basis of the database of (a) [25], (b) [14], (c) [31], (d) [30], (e) present study, (f) [15], (g) [60] and (h) [33].
Figure 9. Genetic distances of ITS1-5.8S-ITS2 from the previous and present studies. Data were obtained from or calculated on the basis of the database of (a) [25], (b) [14], (c) [31], (d) [30], (e) present study, (f) [15], (g) [60] and (h) [33].
Microorganisms 07 00493 g009
Microorganisms 07 00493 g010 550
Figure 10. Putative phylogenetic relationships of all available Pseudokeronopsis species. *, previously misidentified species. Synapomorphies of major clades were numbered. The number 10 indicates the feature that contractile vacuole is located at mid-body or the anterior portion of body. CV, contractile vacuole; Ma, macronuclear nodules; MC, midventral complex; PT, pretransverse cirri; TC, transverse cirri.
Figure 10. Putative phylogenetic relationships of all available Pseudokeronopsis species. *, previously misidentified species. Synapomorphies of major clades were numbered. The number 10 indicates the feature that contractile vacuole is located at mid-body or the anterior portion of body. CV, contractile vacuole; Ma, macronuclear nodules; MC, midventral complex; PT, pretransverse cirri; TC, transverse cirri.
Microorganisms 07 00493 g010
Table
Table 1. Specifications of the Pseudokeronopsis spp. in the present research. Newly obtained sequences are in bold.
Table 1. Specifications of the Pseudokeronopsis spp. in the present research. Newly obtained sequences are in bold.
OrganismsSampling LocationReferencesAccession Numbers
SSU rDNAITS1-5.8S-ITS2LSU rDNACOI
Pseudokeronopsis rubraYellow Sea, Qingdao, ChinaPresent studyMG457184MG457183MG457181MH513652, MH513653
Pseudokeronopsis pararubraYellow Sea, Qingdao, ChinaPresent studyKU663902KU663904MG457180MH513647, MH513648
Pseudokeronopsis carneaYellow Sea, Qingdao, ChinaPresent studyKU663901KU663903 MG457179MH513645, MH513646
P. carnea * (C5)Yellow Sea, Incheon harbor, Korea[27,35]JN714476--MG594874
P. carnea * (C5)Yellow Sea, Gwangyang, Korea[27]MG603616--MG594875
P. carnea *Yellow Sea, Qingdao, China[23]KT984168---
P. carnea pop-I * (C3)Yellow Sea, Qingdao, China[20]AY881633DQ503580, EF174292-4JQ424836-
P. carnea pop-II *Yellow Sea, Qingdao, China[20]-DQ503581--
P. carnea pop-III *Yellow Sea, Qingdao, China[20]-DQ503582, EF174295-7--
Pseudokeronopsis erythrinaEstuary, Yancheng, ChinaPresent studyMH478194MH481652, MH517355MH478198, MH478200MH513643, MH513644
P. erythrinaEstuary of Pearl River, Guangzhou, China[18,24]FJ775723 JQ424868JQ424852-
P. erythrinaEstuary of Pearl River, Guangzhou, China[23]KT984173 ---
P. erythrinaFreshwater River Yamuna, Delhi, IndiaKaur et al., unpublishedMG994990---
P. erythrinaLake Trasimeno, Perugia, Italy[22]KX459375---
Pseudokeronopsis cf. songiYellow Sea, Qingdao, ChinaPresent studyMH513618, MH513619MH513613, MH513614MH513615, MH513616MH513649, MH513650
Pseudokeronopsis songiYellow Sea, Qingdao, ChinaPresent studyKY313623KY313624MG457182MH513651
P. rubra * (C4)Sea of Japan, Hwajinpo, Korea[27]MG603620 -MG594873
P. rubra * (C1)Yellow Sea, Qingdao, China[20]DQ640314DQ640313--
P. rubra *-Khan and Shin, unpublishedHM140387---
P. rubra *-[36]EF535729---
P. rubra *Yellow Sea, Qingdao, China[27]KT984169---
P. rubra *South China Sea, Guangdong, China[27]KT984170---
P. rubra *Daya Bay, Guangdong, China[27]KT984171---
P. rubra *Mangrove wetland, Shenzhen, China[27]KT984172---
P. rubra * (C2)Yellow Sea, Incheon, KoreaJung et al., unpublishedHQ228548HQ228548HQ228548-
Pseudokeronopsis flavaSea of Japan, Pohang-is, Korea[27]MG603616--MG594872
P. flavaSouth China Sea, Zhanjiang, China[20]AY881634DQ503579JQ424835-
P. flavaYellow Sea, Qingdao, China[20]DQ227798---
P. flavaMangrove wetland, Shenzhen, China[27]KT984174---
P. flavaYellow Sea, Qingdao, China[27]KT984175 ---
P. flava-Khan and Shin, unpublishedHM140386 -HM122031-
Pseudokeronopsis sp.Yellow Sea, Qingdao, China[24]JQ424830 JQ424859JQ424847-
Pseudokeronopsis sp.Horniman Museum and Gardens Aquarium, London, England[37]KP793002---
* Misidentified or dubious species. C1–5 correspond to the cryptic species (C1) related to P. songi, the cryptic species (C3) related to P. pararubra, and the uncertain species (C2, C4, and C5) in Table 2 and the following phylogenetic trees. SSU rDNA: small subunit ribosomal DNA, LSU rDNA: large subunit ribosomal DNA, COI: cytochrome c oxidase subunit 1.
Table
Table 2. Emendation for the misidentified or uncertain Pseudokeronopsis sequences in GenBank.
Table 2. Emendation for the misidentified or uncertain Pseudokeronopsis sequences in GenBank.
SequencesReferenceMorphological Data *Molecular AnalysesEmendation
P. carnea DQ503581[20]Similar to P. pararubra and P. rubraNo difference from P. pararubra in ITS1-5.8S-ITS2.P. pararubra
P. carnea DQ503582, EF174295-7[20]Similar to P. pararubra and P. rubraNo difference from P. rubra in ITS2 sequence.P. rubra
P. carnea DQ503580, EF174292-4, JQ424836, AY881633 (C3)[20]Similar to P. pararubra and P. rubraThe ITS1-5.8S-ITS2 genetic distances between C3 and congeners were in the range of 3.23–10.74%, similar to those among congeners (3.68–13.01%), indicating a cryptic species. Cryptic species related to P. pararubra
P. rubra DQ640313, DQ640314 (C1)[20]Similar to P. songiThe ITS1-5.8S-ITS2 genetic distances between C1 and congeners were in the range of 3.46–11.93%, indicating a cryptic species. Cryptic species related to P. songi
P. rubra MG603620, MG594873 (C4)[27]Not availableThe genetic distances between C4 and congeners were in the range of 13.55–21.54%, indicating an uncertain or cryptic species.Misidentified, uncertain species
P. flava MG594872, MG603616[27] Not availableThe COI genetic distance between P. flava MG594872 and P. cf. songi was 2.58%, whereas those among congeners ranged from 13.25% to 21.28%, indicating it conspecific with P. cf. songi. P. cf. songi
Pseudokeronopsis sp. JQ424859, JQ424830[24]Not availableThe ITS1-5.8S-ITS2 genetic distance between JQ424859 and P. songi was 0.21% (vs. conspecific distances 0–0.42%), indicating that the former belonged to P. songi.P. songi
P. carnea JN714476, MG603616, MG594874, MG594875 (C5) [27] Not availableThe COI genetic distance between P. carnea MG594874 and MG594875 was zero, whereas those between them and congeners were much higher than conspecific ones (13.55–20.85% vs. 0–2.58%), suggesting they represented an uncertain species.Misidentified, uncertain species
P. rubra HQ228548 (C2), P. flava HM122031Jung et al. (unpublished), Khan & Shin (unpublished)Not availableThe ITS1-5.8S-ITS2 genetic distance between C2 and congeners were in the range of 3.67–9.73%, suggesting it represented an uncertain species. P. flava HM122031 and C2 were conspecific, possessing the identical LSU-D2 sequence.Misidentified, uncertain species
P. carnea KT002202-07, P. flava KT002234-44, KR263877-79, KR263883-89, P. rubra HM140387, EF535729, KT984169-72, KT002196-201, KT002222-33, KT002245-48, KR263874-76, KR263880-82, KR263888-91[23,36]Not availableOnly SSU rDNA and Actin I sequences were available for these populations, showing six, one, and three nucleotide differences of SSU-V4 region from the valid P. carnea, P. flava, and P. rubra, respectively. Dubious
P. carnea GQ258110[57]Not availableOnly alpha-tubulin gene sequence was available.Dubious
P. erythrina KX459375, MG994990 [22], Kaur et al. unpublishedNot availableOnly SSU rDNA sequences were available, showing 1 and 23 nucleotide differences of SSU rDNA from the valid P. erythrina, respectively.Dubious
The morphological data is shown in Figure 10.
Table
Table 3. Compensating base changes of ITS2 between Pseudokeronopsis species/populations.
Table 3. Compensating base changes of ITS2 between Pseudokeronopsis species/populations.
Species/Populations12345678910
1P. pararubra KU663904, P. carnea DQ503581 *0
2P. carnea DQ503580, EF174292-7 * (C3)00
3P. carnea KU663903 11-
4P. songi KY313624, Pseudokeronopsis sp. JQ4248592220
5P. rubra DQ640313 * (C1)1111-
6P. cf. songi MH513613, MH513614111000
7P. rubra HQ228548 * (C2)111000-
8P. flava DQ5035792112110-
9P. erythrinaMH481652, MH517355, JQ424868211100010
10P. rubraMG457183, P. carnea DQ503582 *0011101110
* Misidentified species. C1–3 correspond to the C1–3 taxa in the ITS1-5.8S-ITS2 tree. Newly obtained sequences are in bold.
Table
Table 4. Compensating base changes of LSU-D2 between Pseudokeronopsis species/populations.
Table 4. Compensating base changes of LSU-D2 between Pseudokeronopsis species/populations.
Species/Populations123456789
1P. songi, Pseudokeronopsis sp. JQ4248470
2P. cf. songiMH513615, MH51361600
3P. flava JQ42483511-
4P. erythrinaMH478198, MH478200, JQ4248522300
5P. carneaMG4571797879-
6P. pararubraMG45718056581-
7P. rubraMG457181786821-
8P. carnea JQ424836 * (C3)5667111-
9P. rubra HQ228548 *, P. flava HM122031 (C2)000175750
* Misidentified species. C2–3 correspond to the C2–3 taxa in the ITS1-5.8S-ITS2-5′LSU tree. Newly obtained sequences are in bold.

Share and Cite

MDPI and ACS Style

Zhan, Z.; Li, J.; Xu, K. Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora). Microorganisms 2019, 7, 493. https://doi.org/10.3390/microorganisms7110493

AMA Style

Zhan Z, Li J, Xu K. Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora). Microorganisms. 2019; 7(11):493. https://doi.org/10.3390/microorganisms7110493

Chicago/Turabian Style

Zhan, Zifeng, Ju Li, and Kuidong Xu. 2019. "Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora)" Microorganisms 7, no. 11: 493. https://doi.org/10.3390/microorganisms7110493

APA Style

Zhan, Z., Li, J., & Xu, K. (2019). Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora). Microorganisms, 7(11), 493. https://doi.org/10.3390/microorganisms7110493

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