Supplementary_materials.docx

Table S1.
List
of species analyzed in this study.
Table S2. Genetic divergence among total 69 populations and
within genera. Raw data (upper line) and changed data including cryptic species
(lower line) are shown.

Figure
S1. Sampling locations. Detailed information is
shown in Table S1.

Figure S2. Schematic illustration of the sequencing
regions of CO1. The primer set
designed in this study results in amplification of a 478 base pair (bp)
fragment of the CO1 gene; which
overlaps with 260–348 bp of a previously reported sequencing region
(Strüder-Kypke and Lynn 2010). *The range is based on complete CO1 sequences: Paramecium caudatum (FN424190). **Strüder-Kypke and Lynn (2010, p.
133) used different primer pairs that resulted in different amplicon length. Of
the primer sets F388dT and R1184dT, which were used to amplify sequences from five
classes (Oligohymenophorea, Nassophorea, Colpodea, Spirotrichea, and
Heterotrichea), were selected as the sequencing region in this study, as shown
in the illustration. The indel region is marked approximately because we could
not determine this region by means of its various alignment patterns