Phylogenetic analyses

For phylogenetic relationships, analyses were based on 18S and 28S rDNA. The newly obtained sequences were manually edited using BioEdit 7.2.6 (Hall, 1999) and aligned with another 18S or 28S rRNA gene sequences available in GenBank using Muscle alignment tool implemented in the MEGA7 (Kumar et al., 2016). The ambiguously aligned parts and divergent regions were known using the online version of Gblocks 0.91b (Castresana, 2000) (http://molevol.cmima.csic.es/castresana/Gblocks_server.html) and were removed from the alignments using MEGA7. The best-fit model of nucleotide substitution used for the phylogenetic analysis was statistically selected using jModelTest 2.1.10 (Darriba et al., 2012). Phylogenetic tree was generated with Bayesian inference method using MrBayes 3.2.6 (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003). Myolaimus byersi ({"type":"entrez-nucleotide","attrs":{"text":"KU180665","term_id":"1003466857","term_text":"KU180665"}}KU180665 for 18S and {"type":"entrez-nucleotide","attrs":{"text":"KU180676","term_id":"1003466868","term_text":"KU180676"}}KU180676 for 28S) was chosen as outgroup according to previous results by Kanzaki et al. (2009). The analysis under GTR+I+G model was initiated with a random starting tree and run with the Markov Chain Monte Carlo (MCMC) for 1 × 10⁶ generations. The tree was visualized and saved with FigTree 1.4.3 (Rambaut, 2014).