Supplementary_materials.docx

<b>Table S1. </b>List of species analyzed in this study.<b>Table S2.</b> Genetic divergence among total 69 populations and within genera. Raw data (upper line) and changed data including cryptic species (lower line) are shown.<br><b>Figure S1. </b>Sampling locations. Detailed information is shown in Table S1. <b>Figure S2.</b> Schematic illustration of the sequencing regions of <i>CO1</i>. The primer set designed in this study results in amplification of a 478 base pair (bp) fragment of the <i>CO1</i> 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 <i>CO1</i> sequences: <i>Paramecium caudatum </i>(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<br>