Hostname: page-component-586b7cd67f-t7czq Total loading time: 0 Render date: 2024-11-24T10:28:05.268Z Has data issue: false hasContentIssue false

Morphological and molecular phylogenetic analyses of Zepedanulus ishikawai (Arachnida: Opiliones: Laniatores: Epedanidae) in the southern part of the Ryukyu Archipelago

Published online by Cambridge University Press:  16 September 2021

Yoshimasa Kumekawa*
Affiliation:
Ehime University, The United Graduate School of Agricultural Sciences, Environmental Sciences, 200 Monobe, Nankoku, Kochi, 783-8502, Japan
Haruka Fujimoto
Affiliation:
Kochi University, Graduate School of Integrated Arts and Sciences, 200 Monobe, Nankoku, Kochi, 783-8502, Japan
Osamu Miura
Affiliation:
Kochi University, Faculty of Agriculture and Marine Science, 200 Monobe, Nankoku, Kochi, 783-8502, Japan
Ryo Arakawa
Affiliation:
Kochi University, Graduate School of Integrated Arts and Sciences, 200 Monobe, Nankoku, Kochi, 783-8502, Japan
Jun Yokoyama
Affiliation:
Yamagata University, Faculty of Science, 1-4-12 Kojirakawa-machi, Yamagata, 990-8560, Japan
Tatsuya Fukuda
Affiliation:
Tokyo City University, Faculty of Knowledge Engineering, 1-28-1, Tamazutsumi, Setagaya, Tokyo, 158-8557Japan Tokyo City University, Faculty of Science and Engineering, 1-28-1, Tamazutsumi, Setagaya, Tokyo, 158-8557Japan
*
*Corresponding author. Email: [email protected]

Abstract

Harvestmen (Arachnida: Opiliones) are soil animals with extremely low dispersal abilities that experienced allopatric differentiation. To clarify the morphological and phylogenetic differentiation of the endemic harvestman Zepedanulus ishikawai (Suzuki, 1971) (Laniatores: Epedanidae) in the southern part of the Ryukyu Archipelago, we conducted molecular phylogenetic analyses and divergence time estimates based on CO1 and 16S rRNA sequences of mtDNA, the 28S rRNA sequence of nrDNA, and the external morphology. A phylogenetic tree based on mtDNA sequences indicated that individuals of Z. ishikawai were monophyletic and were divided into clade I and clade II. This was supported by the nrDNA phylogenetic tree. Although clades I and II were distributed sympatrically on all three islands examined (Ishigaki, Iriomote, and Yonaguni), heterogeneity could not be detected by polymerase chain reaction–restriction fragment length polymorphism of nrDNA, indicating that clades I and II do not have a history of hybridisation. Also, several morphological characters differed significantly between individuals of clade I and clade II. The longstanding isolation of the southern Ryukyus from the surrounding islands enabled estimation of the original morphological characters of both clades of Z. ishikawai.

Type
Research Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of the Entomological Society of Canada

Introduction

Continental archipelagos are suitable for examining causal connections between geographical barriers, disjunct distributions, and genetic differentiation (Hedges et al. Reference Hedges, Hass and Maxon1992; Hisheh et al. Reference Hisheh, Westerma and Schmitt1998; Atkins et al. Reference Atkins, Presto and Cronk2001; Ohdachi et al. Reference Ohdachi, Dokuchaev, Hasegawa and Masuda2001; Poulakakis et al. Reference Poulakakis, Lymberakis, Antonioub, Chalkia, Zourosb, Mylonasa and Valakos2003; Bittkau and Comes Reference Bittkau and Comes2005; Velo-Antón et al. Reference Velo-Antón, Zamudio and Cordero-Rivera2012). The Ryukyu Archipelago was the eastern margin of continental East Asia in the middle to late Miocene (15.97–5.333 million years ago). In the late Miocene to the early Pleistocene (5.333–0.774 million years ago), it fragmented into large islands and now includes approximately 140 islands scattered in an arc between Kyushu Island of Japan and Taiwan (Kimura Reference Kimura and Kimura2002; Osozawa et al. Reference Osozawa, Pavlis, Flower, Wakabayashi and Dilek2011). The islands are the product of the repeated formation and division of a land bridge to continental Asia that extended through Taiwan and Kyushu. Expansion and regression of the land bridge have been attributed to climatic oscillations in the Pliocene and Quaternary periods (Kimura Reference Kimura1996, Reference Kimura2000; Kizaki and Ohshiro Reference Kizaki and Oshiro1977, Reference Kizaki, Oshiro and Kizaki1980; Ujiie Reference Ujiie and Ujiie1990; Ohshiro Reference Ohshiro, Nishijima, Nishida, Shikatani and Shokita2003; Keally Reference Keally2005; Shinjo Reference Shinjo, Fujita, Arakaki, Denda, Hidaka, Hirose and Reimer2015). Plants and animals on the islands consist largely of relictual taxa that presumably differentiated from their relatives on the adjacent continent (Ota Reference Ota1998; Reference Ota2012). A number of biogeographic and evolutionary studies on the organisms of the Ryukyu Archipelago have been conducted and have enhanced our understanding of the biogeographic and evolutionary histories of organisms in continental archipelagos (e.g., Ota Reference Ota1998; Maekawa et al. Reference Maekawa, Lo, Kitade, Miura and Matsumoto1999; Seo et al. Reference Seo, Watabanabe, Hotta and Murakami2004; Chiang and Schaal Reference Chiang and Schaal2006; Seki et al. Reference Seki, Sakanashi, Kawaji and Kotaka2007; Nakamura et al. Reference Nakamura, Denda, Kokubugata, Huang, Peng and Yokota2015; Kaito and Toda Reference Kaito and Toda2016).

In general, the low dispersal ability of flightless organisms leads to a low rate of gene flow, suggesting accelerated differentiation among populations (Zera Reference Zera1981; Hansen Reference Hansen1983; Jablonski Reference Jablonski1986; Smith and Farrell Reference Smith and Farrell2006; Miura et al. Reference Miura, Torchin, Bermingham, Jacobs and Hechinger2012). For example, Ikeda et al. (Reference Ikeda, Nishikawa and Sota2011) reported that the loss of flight accelerates allopatric speciation among carrion beetles (Coleoptera: Silphidae). Flightless soil animals have also been found to exhibit extensive differentiation among allopatric populations (Tsurusaki Reference Tsurusaki2006a). In the Ryukyu Archipelago, for example, Kilungius insulanus (Hirst, Reference Hirst1911) (Opiliones: Epedanidae) is a species of harvestman that is genetically differentiated between Amami-Oshima and Okinawa Islands in the middle of the archipelago (Kumekawa et al. Reference Kumekawa, Ito, Miura, Yokoyama, Tebayashi, Arakawa and Fukuda2015). Zepedanulus ishikawai (Suzuki, Reference Suzuki1971) (Opiliones: Epedanidae) is another species of harvestman, belonging to the same family as K. insulanus. It is approximately 5 mm in body length, lives in humid places under forest litter, and is distributed in the southern Ryukyu Archipelago, namely Miyako, Ishigaki, Iriomote, Yonaguni, and Senkaku Islands (Fig. 1; Suzuki Reference Suzuki1973; Tsurusaki and Suzuki, Reference Tsurusaki, Suzuki and Aoki2015). This scattered distribution of Z. ishikawai might promote the differentiation of populations on different islands. In fact, allopatric speciation, where geographic isolation results in speciation, is widely accepted (Mayr Reference Mayr1963; Coyne and Orr Reference Coyne and Orr2004; Bolnick and Fitzpatrick Reference Bolnick and Fitzpatrick2007), and the accidental differentiation such as a genetic drift found between allopatric populations has led to substantial postzygotic isolation it is also known to bring about (Dobzhansky Reference Dobzhansky1937). The presence of physical barriers, such as oceans, rivers, or high mountains, will cause populations isolated by these boundaries to undergo different accumulations of variation over time from isolation. If reproductive isolation is acquired through this, it is likely that even if the isolation is subsequently resolved, the populations will not be able to interbreed with each other, and each population will have its own history.

Fig. 1. Distribution of Z. ishikawai in the Ryukyu Archipelago, Japan.

Recently developed molecular techniques allow phylogenies to be reconstructed and the systematic relationships among species or populations within a given species to be evaluated. Application of mitochondrial DNA as a marker for phylogenetic and phylogeographic surveys enables evaluation of intraspecific population structures and differentiation (e.g., Avise Reference Avise2000, Reference Avise2009). Several phylogenetic and phylogeographic studies of harvestman have been based on mitochondrial DNA sequences (Giribet et al. Reference Giribet, Rambla, Carranza, Baguna, Riutort and Ribera1999; Thomas and Hedin Reference Thomas and Hedin2008; Sharma and Giribet Reference Sharma and Giribet2009, Reference Sharma and Giribet2011; Derkarabetian et al. Reference Derkarabetian, Steinmann and Hedin2010, Reference Derkarabetian, Ledford and Hedin2011; Giribet et al. Reference Giribet, Vogt, Perez, Sharma and Kury2010; Hedin and Thomas Reference Hedin and Thomas2010; Burns et al. Reference Burns, Hedin and Shultz2012; Sharma Reference Sharma2010; Sharma et al. Reference Sharma, Buenavente, Clouse, Diesmos and Giribet2012; Kumekawa et al. Reference Kumekawa, Ito, Tsurusaki, Hayakawa, Ohga and Yokoyama2014, Reference Kumekawa, Ito, Miura, Yokoyama, Tebayashi, Arakawa and Fukuda2015; Didodemico and Hedin Reference Didodemico and Hedin2016; Emata and Hedin Reference Emata and Hedin2016; Cruz-López et al. Reference Cruz-López, Monjaraz-Ruedas and Francke2019). It is also widely accepted that selecting suitable genetic markers is important in phylogeny and phylogeography studies, and a number of studies using mitochondrial DNA data have demonstrated the utility of the mitochondrial genes cytochrome c oxidase subunit I (CO1) and 16S ribosomal RNA (rRNA). Particularly for shallow-level animal molecular phylogeny and phylogeography, these genes enable accurate identification at the species (or lower) level through comparative analyses of sequence variations in short, standardised fragments of the genome (Hebert et al. Reference Hebert, Cywinska, Ball and deWaard2003; Jinbo et al. Reference Jinbo, Kato and Ito2011). We thus assessed the differentiation among Z. ishikawai populations on different islands using the sequences of CO1 and 16S rRNA in mitochondrial DNA.

The congruence of multiple genetic trees enables assessment of the reliability of a hypothesis regarding the evolutionary relationships of organisms. Indeed, even phylogenetic incongruence may provide a window into evolutionary processes that cannot be accessed using a single gene data set. Although a mitochondrial DNA phylogenetic tree allows clarification of the evolutionary relationships of a given species, it cannot clarify certain evolutionary processes because of its maternal inheritance. Theoretically, incongruence between a genetic tree and species phylogeny, or among different genetic trees, can result from a variety of evolutionary processes, such as conflation of orthology and paralogy, incomplete lineage sorting of ancestral polymorphism, and introgressive hybridisation (Rieseberg and Soltis Reference Rieseberg and Soltis1991; Wendel and Doyle Reference Wendel, Doyle, Soltis, Soltis and Doyle1998; Maddison and Knowles Reference Maddison and Knowles2006; Toews and Brelsford Reference Toews and Brelsford2012). Of these, introgressive hybridisation and incomplete lineage sorting frequently lead to phylogenetic incongruence among different data sets at lower taxonomic levels (Avise et al. Reference Avise, Shapiro, Daniel, Aquadro and Lansman1983; Cronn and Wendel Reference Cronn and Wendel2004). Incongruence among genetic trees should be confirmed using data from nuclear DNA in addition to mitochondrial DNA. In fact, Schlick-Steiner et al. (Reference Schlick-Steiner, Steiner, Seifert, Stauffer, Christian and Crozier2010) claimed that it is necessary to analyse various types of information to clarify the relationships among and differentiation within, populations within of a species. Among nuclear DNA sequences, 28S rRNA is particularly useful for differentiating taxa among groups of organisms (Nunn et al. Reference Nunn, Theisen, Christensen and Arctander1996; Friedrich and Tautz Reference Friedrich and Tautz1997; Babcock et al. Reference Babcock, Heraty, De Barro, Driver and Schmidt2001; Schmidt et al. Reference Schmidt, Driver and De Barro2006; Taylor and Rogers Reference Taylor and Rogers2015; Che et al. Reference Che, Wang, Shi, Du, Zhao, Lo and Wang2016), suggesting that this gene could be used to resolve the phylogenetic relationships of Z. ishikawai. In addition, the analysis of polymerase chain reaction–restriction fragment length polymorphisms has been used to show whether hybridisation and introgression occur in a variety of species (Hollingsworth et al. Reference Hollingsworth, Bailey, Hollingsworth and Ferris1999; Nijman et al. Reference Nijman, Bradley, Hanotte, Otsen and Lenstra1999; Mráz et al. Reference Mráz, Chrtek, Fehrer and Plačková2005; Oleksa et al. Reference Oleksa, Chybicki, Tofilski and Burczyk2011; Vaini et al. Reference Vaini, Grisolia, Prado and Porto-Foresti2014; Kumekawa et al. Reference Kumekawa, Miura, Fujimoto, Ito, Arakawa, Yokoyama and Fukuda2019).

Therefore, we evaluated the differentiation of Z. ishikawai populations among islands using 28S rRNA in nuclear DNA and compared the results with gene fragments of mitochondrial DNA. We also analysed the morphological characters of Z. ishikawai to determine whether the populations could be distinguished morphologically.

Materials and methods

Sample collection

We collected Z. ishikawai by hand from three major islands of the Yaeyama Islands, in the southernmost part of the Ryukyu Archipelago (Ishigaki, Iriomote, and Yonaguni Islands; Table 1; Fig. 2).

Table 1. Summary of Zepedanulus ishikawai analysed in this study. Specimen number corresponds to our phylogenetic results (Figs. 3, 4, 5, and 6).

Fig. 2. (A) Spined type and (B) nonspined type, and characters of Zepedanulus ishikawai (collected from Yonaguni Island): A, BL – body length; B, CL – carapace length; C, CW – carapace width; D, Fe1L – femur length, first leg; E, Fe2L – femur length, second leg; F, Fe3L – femur length, third leg; G, Fe4L – femur length, fourth leg; H, SL – spine length; I, Chel.L – chelicerae length; and J, PTL – palpal–tarsus length, and the morphology of the second scutal area. Arrowhead indicates the spine. Scale bar = 1 mm.

DNA extraction, amplification, sequencing, and polymerase chain reaction–restriction fragment length polymorphism analysis

We extracted DNA from a leg using DNeasy kits (Qiagen, Valencia, California, United States of America), according to the manufacturer’s protocol for animal tissue samples. The isolated DNA was resuspended in an appropriate volume of Tris-ethylenediaminetetraacetic acid buffer (pH 8.0) and stored at −20 °C until use. For all specimens, we amplified the CO1 and 16S rRNA genes in mitochondrial DNA and the 28S rRNA gene in nuclear DNA. The CO1 gene fragments were amplified and sequenced using primers published previously (Kumekawa et al. Reference Kumekawa, Ito, Tsurusaki, Hayakawa, Ohga and Yokoyama2014). The 16S rRNA and 28S rRNA gene fragments were amplified and sequenced using the primers of Xiong and Kocher (Reference Xiong and Kocher1991) and Mallatt and Sullivan (Reference Mallatt and Sullivan1998), respectively. Amplification of DNA followed the method of Kumekawa et al. (Reference Kumekawa, Murjoko, Hayakawa, Ohga, Mori and Miyazaki2013). The reaction solution (50 µL) contained 50 ng of total DNA, 10 mM Tris (hydroxymethyl) aminomethane hydrochloride buffer (pH 8.3) with 50 mM potassium chloride and 1.5 mM magnesium dichloride, 0.2 mM of each deoxyribonucleotide triphosphate, 1.25 U Taq DNA polymerase (TaKaRa, Tokyo, Japan), and 0.5 µM of each primer. The following thermal cycle was applied: incubation at 94 °C for 10 seconds, followed by 45 cycles of incubation at 94 °C for 1.5 minutes, 48 °C for 2 minutes, and 72 °C for 3 minutes, with a final extension at 72 °C for 15 minutes. After amplification, the reaction mixtures were subjected to electrophoresis in 1% low melting–temperature agarose gels, and the products were purified using QuickSpin kits (Qiagen) according to the manufacturer’s specifications. We sequenced the purified polymerase chain reaction products using a BigDye Terminator Cycle Sequencing Kit (Applied Biosystems, Tokyo, Japan) and ABI PRISM 3100-Avant Genetic Analyzer (Applied Biosystems) according to the manufacturers’ instructions.

We also carried out polymerase chain reaction–restriction fragment length polymorphism analysis, because autapomorphic characters (clade II, see Results) of 28S rRNA are involved in the restriction site HaeIII (GGCC). We designed the following primers to assess the length of the restriction fragments and conducted polymerase chain reaction as above, using the primers Zepe28S-RFLP-F: 5'-GCC CTC GAA GTT CGG AGG CG-3' and Zepe28S-RFLP-R: 5'-CAC CCG GCT ACC TAC CGG TT-3'. The amplified products were digested by HaeIII at 37 °C for more than 1 hour. The digested DNAs were separated in a 1.5% agarose gel and the size of each band was determined.

Data analysis

To construct phylogenetic trees for Z. ishikawai and related species (Table 1), 58 nucleotide sequences were aligned using the ClustalW program (Thompson et al. Reference Thompson, Higgins and Gibson1994) implemented in MEGA, version 6, software (Tamura et al. Reference Tamura, Stecher, Peterson, Filipski and Kumar2013). The positions containing gaps and missing data were eliminated manually. A Bayesian analysis was conducted using MrBayes, version 3.1.2, software (Ronquist and Huelsenbeck Reference Ronquist and Huelsenbeck2003). We considered the model with the lowest log-likelihood score to be the most accurate description of the substitution pattern. The Akaike information criterion value of each model was calculated. The best-fit model (mitochondrial DNA, GTR + I + G model; nuclear DNA, F81 model) selected by hLRT in MrModeltest 2.3 (Nylander Reference Nylander2004) and PAUP * 4.0b10 (Swofford Reference Swofford2002) software was used. We ran four Markov chains for 50 million generations and sampled Markov chains every 100 generations; 20% of the runs were discarded as burn-in. The effective sample size (ESS) of each parameter indicated by Tracer version 1.6, software (Rambaut et al. Reference Rambaut, Suchard, Xie and Drummond2014) confirmed that the number of Markov chain generations was sufficient (ESS > 200). The maximum-likelihood method was conducted using MEGA, version 6.06, software (Tamura et al. Reference Tamura, Stecher, Peterson, Filipski and Kumar2013). The model test function in MEGA was used to select the appropriate model for maximum-likelihood analysis by considering the model with the lowest Bayesian information criterion scores as the most accurate description of the substitution pattern. For each model, the Akaike information criterion–corrected value, maximum-likelihood value (lnL), and the number of parameters (including branch lengths) are presented (Nei and Kumar Reference Nei and Kumar2000). The nonuniformity of evolutionary rates among sites may be modelled using a discrete gamma distribution (+ G) with five rate categories and a parameter for invariable sites (+ I). The relative values of instantaneous r should be considered when evaluating them. For simplicity, the sum of r values is made equal to 1 for each model. To estimate maximum-likelihood values, a tree topology was automatically computed. Nodal support in maximum-likelihood trees was measured by bootstrapping (1000 samples).

The divergence times of the clades were estimated for all individuals used for phylogenetic analyses. Divergence times were inferred using an uncorrelated relaxed lognormal clock, implemented in BEAST, version 2.6.2, software (Bouckaert et al. Reference Bouckaert, Vaughan, Barido-Sottani, Duchêne, Fourment and Gavryushkina2019). Because there are no mitochondrial DNA molecular clocks available for Zepedanulus or related harvestman species, we conservatively adopted a wide range of estimated divergence rates for arthropods from 0.21% (Clouse and Wheeler Reference Clouse and Wheeler2014) to 3.34% per million years (Pons and Vogler Reference Pons and Vogler2005) for the CO1 gene, and 0.45% (Gomez-Zurita et al. Reference Gómez-Zurita, Juan and Petitpierre2000) to 1.9% per million years (Percy et al. Reference Percy, Page and Cronk2004) for the 16S gene. Prior distributions for mean substitution rates were sampled from a log normal distribution (clock.rate = 0.006, sigma = 0.8 for the CO1 gene; clock.rate = 0.005, sigma = 0.5 for the 16S gene). We used the same substitution model as the Bayesian analysis and used the Yule process as a tree prior. We constrained the monophyly of two clades. The analysis was run for 500 million generations, sampled every 50 000 steps, and the first 1000 samples were discarded as burn-in. We used Tracer, version 1.6 (Rambaut et al. Reference Rambaut, Suchard, Xie and Drummond2014), and FIGTREE, version 1.4.2 (Rambaut Reference Rambaut2014), software to check for convergence and to visualise the results.

Morphological analysis

For morphological analysis, the body length, carapace length, carapace width, femur lengths of the first to fourth legs, spine length, chelicerae length, and palpal–tarsus length were measured for each Z. ishikawai specimen (Fig. 2). We measured the morphological characters of 139 Z. ishikawai individuals (39 from Ishigaki Island, 34 from Iriomote Island, and 66 from Yonaguni Island). Measurements were performed using cellSens standard imaging software (Olympus Co., Tokyo). Statistical analysis was conducted by Tukey–Kramer test (P < 0.05).

Results

Molecular analysis of Zepedanulus ishikawai and related species

We reconstructed phylogenetic relationships using 43 samples (27 from Ishigaki Island, eight from Yonaguni Island, and six from Iriomote Island) of Z. ishikawai and related species (Table 1), including previously published sequences of the CO1 and 16S rRNA genes in mitochondrial DNA. Of these, even though only one sample of Epedanellus tuberculatus was used in the analysis of the 16S rRNA and CO1 sequences, four samples of Penaeus japonicus and two of K. insulanus were analysed because previous studies indicated that these species are divided into four and two groups, respectively (Kumekawa et al. Reference Kumekawa, Ito, Tsurusaki, Hayakawa, Ohga and Yokoyama2014, Reference Kumekawa, Ito, Miura, Yokoyama, Tebayashi, Arakawa and Fukuda2015). The combined length of the CO1 and 16S rRNA sequences in Z. ishikawai was 765 bp (16S rRNA, 287 bp; CO1, 478 bp). We reconstructed phylogenetic trees based on the maximum-likelihood (Fig. 3) and Bayesian inference (Fig. 4) methods. The branching patterns of the Z. ishikawai maximum-likelihood phylogenetic tree were identical to those of the tree that was based on Bayesian inference (Figs. 3 and 4). The phylogenetic trees based on the CO1 and 16S rRNA genes indicated that individuals of Z. ishikawai composed a monophyletic group with a high bootstrap value (96%) and high posterior probability (1.00). The species Z. ishikawai was divided into two clades, henceforth referred to as clade I and clade II (Fig. 3). Clade I was not strongly supported in the maximum-likelihood tree (68%), but Bayesian analysis indicated a high posterior probability for the clade (1.00). Clade I involved three monophyletic subgroups clades IA, IB, and IC. Clade IA consisted of individuals from Ishigaki Island and was the sister group of clade IB, which comprised from Iriomote Island. Clade IC was the sister group to the combined clade of IA + IB, and comprised samples from Yonaguni Island. These three groups had high bootstrap values (70–100%) and posterior probabilities (0.99–1.00). Clade II, which had high support values (100% maximum-likelihood bootstrap and 1.00 Bayesian posterior probability values), also comprised three monophyletic subgroups clades IIA, IIB, and IIC. The branching patterns of clade II were identical to those of clade I, and individuals of Yonaguni Island, Ishigaki Island, and Iriomote Island constituted monophyletic groups of clades IIA, IIB, and IIC with high support values (100% maximum-likelihood bootstrap and 1.00 Bayesian posterior probability values), respectively.

Fig. 3. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by the maximum-likelihood method based on the sequences of CO1 and 16S rRNA in mtDNA. Numbers above or below branches indicate bootstrap values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Fig. 4. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by Bayesian inference based on the sequences of CO1 and 16S rRNA in mtDNA. Numbers above or below branches indicate posterior probabilities. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

The phylogenetic tree based on the 28S rRNA gene in nuclear DNA was constructed using Z. ishikawai and Tokunosia tenuipes, which was revealed to be a sister species of Z. ishikawai based on the mitochondrial DNA phylogenetic tree (Figs. 5 and 6). The length of the Z. ishikawai 28S rRNA sequence was 978 base pairs. As for mitochondrial DNA data sets, we reconstructed phylogenetic trees based on the maximum-likelihood (Fig. 5) and Bayesian inference (Fig. 6) methods. The major branching patterns of the phylogenetic trees based on 28S rRNA were similar to those of the mitochondrial DNA trees. The phylogenetic trees indicated that individuals of Z. ishikawai comprised a monophyletic group and were divided into two groups. These groups corresponded to clades I and II in the mitochondrial DNA phylogenetic tree.

Fig. 5. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by the maximum-likelihood method based on the sequence of 28S rRNA in nrDNA. Numbers above or below branches indicate the bootstrap values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Fig. 6. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by Bayesian inference based on the sequence of 28S rRNA in nrDNA. Numbers above or below branches indicate the posterior probabilities. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Our phylogenetic analyses of mitochondrial DNA and 28S rRNA sequences both indicated the presence of two identical clades. Clades I and II were sympatrically distributed on all three islands examined, and no incongruence of mitochondrial DNA and nuclear DNA was observed. Therefore, whether hybrids of individuals of clades I and II on each of the three islands could form is unclear because assessing heterogeneity at polymorphic sites in nuclear DNA using the results of Sanger sequencing can be problematic (Simsek et al. Reference Simsek, Tanira, Al-Baloushi, Al-Barwani, Lawatia and Bayoumi2001). Based on the above analyses, we found that the 28S rRNA sequences of clade II had HaeIII restriction sites, whereas those of clade I did not (Fig. 7). We conducted polymerase chain reaction–restriction fragment length polymorphism to evaluate hybridisation and introgression between clade I and clade II of Z. ishikawai. No hybridisation or introgression between clade I and clade II of Z. ishikawai was detected on the three islands (Fig. 7).

Fig. 7. A, Restriction sites in the 28S rRNA sequence determined by polymerase chain reaction–restriction fragment length polymorphism (M, HaeIII site) and B, polymerase chain reaction–restriction fragment length polymorphism profile of clades I and II of Zepedanulus ishikawai (M, size markers).

The dated tree of Z. ishikawai is shown in Figure 8. The tree exhibited that the divergence of clade I and clade II occurred about 28 million years ago (47–13 million years ago: 95% 4-hydroxyphenylpyruvate dioxygenase). In addition, the divergence time of the three groups of clade I occurred about 25 million years ago, and the divergence time of the three groups of clade II occurred about 10 million years ago.

Fig. 8. Estimated divergence time (over millions of years ago; MYA) for the two clades of Zepedanulus ishikawai. Horizontal bars represent the upper and lower 95% of highest posterior density values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Morphological analysis of Zepedanulus ishikawai

All examined morphological characters of clade I were significantly larger than those of clade II (Table 2; Figs. 911). No significant differences within the same clade were found among the islands. The spine length could not be measured for individuals of clade I because they lacked an eye-spine and spine on the second suctal area. Therefore, the spine length was measured only for individuals of clade II.

Table 2. The measurement (average ± standard deviation) of Zepedanulus ishikawai.

Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; PTL, palpal–tarsus length; Chel.L, chelicerae length.

Fig. 9. Histograms of the morphological characteristics of Zepedanulus ishikawai in Iriomote Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Fig. 10. Histograms of the morphological characteristics of Zepedanulus ishikawai in Ishigaki Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Fig. 11. Histograms of the morphological characteristics of Zepedanulus ishikawai in Yonaguni Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Discussion

A phylogenetic tree of the mitochondrial DNA sequences of Z. ishikawai indicated that the species consisted of two phylogenetic groups clades I and II which was supported by the phylogenetic tree based on the 28S rRNA gene. Our morphological analyses indicated that individuals of clade I could be distinguished from those of clade II not only quantitatively (significant differences in body length, carapace length, carapace width, chelicerae length, palpal–tarsus length, spine length, and the femur lengths of the first to fourth legs) but also qualitatively (presence of an eye-spine and a spine in the second suctal area). Both groups contain individuals on three islands (Ishigaki, Iriomote, and Yonaguni), suggesting that both groups are distributed sympatrically throughout the southern Ryukyu Archipelago. This has consequences for the diversification of Z. ishikawai, because two morphologically and genetically differentiated groups were apparent at the species level. We did not detect incongruence of the branching patterns of the phylogenetic trees based on mitochondrial DNA and nuclear DNA sequences (Figs. 3, 4, 5, and 6) or on nuclear DNA heterogeneity by polymerase chain reaction–restriction fragment length polymorphism, despite collecting samples from sympatric populations on each island (Fig. 7). In addition, the body colour differed between clades on all the islands (clade I, yellow; clade II, dark brown on Ishigaki and Iriomote Islands, and clade I, brown; clade II, black on Yonaguni Island; Fig. 12). We therefore infer that these two groups represent reproductively isolated species. Suzuki (Reference Suzuki1971) reported that the total body length of Z. ishikawai was 3.65 mm in males and 2.81–3.24 mm in females and that the legs of the type specimen of Z. ishikawai were comparatively long. Although a new taxonomic treatment is needed, based on Suzuki’s description of legs and total body length, clade I is Z. ishikawai and a new name is needed for clade II.

Fig. 12. Distribution of morphological characteristics in a simplified phylogenetic tree of Zepedanulus ishikawai and related species.

The phylogeographic patterns of organisms on the islands have been influenced by the formation of the continental island systems, which offer relatively simple arenas for the evolutionary dynamics of community assembly, because harvestmen are generally small and spatially isolated beyond their vagility. Prior phylogeographic studies reported that the current distribution of biological diversity cannot be understood without information about how organisms responded to historical changes in geological and climatic conditions (e.g., Arbogast and Kenagy Reference Arbogast and Kenagy2001; Fukuda et al. Reference Fukuda, Yokoyama and Ohashi2001; Griffin and Barrett Reference Griffin and Barrett2004; Hewitt Reference Hewitt2004; Vink et al. Reference Vink, Thomas, Paquin, Hayashi and Hedin2005; Matsumura et al. Reference Matsumura, Yokoyama, Fukuda and Maki2009; Minamiya et al. Reference Minamiya, Yokoyama and Fukuda2009; Fukuda et al. Reference Fukuda, Song, Ito, Hayakawa, Minamiya, Kanno, Hayakawa, Minamiya and Yokoyama2011; Maura et al. Reference Maura, Salvi, Bologna, Nascetti and Canestrelli2014; Mezzasalma et al. Reference Mezzasalma, Dall’Asta, Loy, Cheylan, Lymberakis and Zuffi2015). Kumekawa et al. (Reference Kumekawa, Ito, Tsurusaki, Hayakawa, Ohga and Yokoyama2014) suggested that Pseudobiantes japonicus comprised four allopatric groups with morphological and genetic differences, and based on the current distribution and historical changes of suitable habitats, these groups had undergone independent northward spread by various routes from allopatric refugia distributed along the coasts of the Pacific Ocean in the western part of the main island of Japan. It is intriguing to consider how Z. ishikawai colonised the southern part of the Ryukyu Archipelago despite the presence of a barrier to expansion. The southern part of the Ryukyu Archipelago formed the eastern margin of continental East Asia in the mid to late Miocene (15.97–5.333 million years ago) and then fragmented into large islands, breaking the direct land connection to China, in the late Miocene to the early Pleistocene (5.333–0.774 million years ago; Kimura Reference Kimura and Kimura2002; Osozawa et al. Reference Osozawa, Shinjo, Armid, Watanabe, Horiguchi and Wakabayashi2012). This suggests that Z. ishikawai colonised the archipelago region from continental East Asia during this period. Subsequently, sea levels ranged from 200 m higher to 200 m lower than present-day sea level, likely due to climatic oscillations during the Quaternary period. Lower islands were then either largely or completely submerged during the interglacial age (Kimura Reference Kimura1996). During this period, subsidence created two deep-water passages through the island arc, namely the Tokara Gap and the Kerama Gap, which divide the Ryukyu area into three groups: the northern, central, and southern Ryukyus (Osozawa et al. Reference Osozawa, Shinjo, Armid, Watanabe, Horiguchi and Wakabayashi2012). Although the islands were repeatedly fragmented into smaller islands and reconnected in the middle to late Pleistocene (0.774–0.0017 million years ago), the existence of land connections across the two gaps after their opening is unlikely, even during the Quaternary glacial sea level minima, because the ocean depth is > 1000 m between the constituent island groups (Ota Reference Ota1998; Kawana Reference Kawana and Kimura2002; Osozawa et al. Reference Osozawa, Shinjo, Armid, Watanabe, Horiguchi and Wakabayashi2012). These gaps acted as strong barriers, even for flying insect species, and resulted in the generation of endemic species (e.g., Osozawa et al. Reference Osozawa, Su, Oba, Yagi, Watanabe and Wakabayashi2013, Reference Osozawa, Takáhashi and Wakabayashi2015). Tsurusaki (Reference Tsurusaki2006b) reported that Opiliones fauna is divided into three categories along the above island groups, based on cluster analysis using the similarity index. Therefore, Z. ishikawai could not cross the Kerama Gap and is included in the third category of Tsurusaki (Reference Tsurusaki2006b). Alternatively, the land connection between Yonaguni Island, which is about 130 km from Taiwan Island and is the westernmost area in the distribution of Z. ishikawai, and Taiwan Island ceased in the early Pleistocene (2.58 million years ago), roughly contemporaneous with disconnection of the southern Ryukyus from continental East Asia (Koba Reference Koba1992; Osozawa et al. Reference Osozawa, Shinjo, Armid, Watanabe, Horiguchi and Wakabayashi2012). The distance between the two islands could not prevent the movement of flying organisms, such as winged insects and birds, nor the spread of plants with long-distance dispersal abilities (Hatusima Reference Hatusima1975; Ota Reference Ota1998). However, the flightlessness and low vagility of Z. ishikawai could restrict its movement between Yonaguni and Taiwan Islands, and this species likely experienced the long-standing isolation of the southern Ryukyus from the surrounding islands. The deep genetic divergences we observed within the Z. ishikawai species complex suggest that the geological history of the Ryukyus played key roles in driving geographic isolation of populations since their divergence from their common ancestor in the archipelago.

The endemism of Z. ishikawai in the southern Ryukyus is supported by the historical isolation of the area from surrounding regions. However, our phylogenetic results indicated that Z. ishikawai comprises two groups and three subgroups; the former are considered different species, and the latter are considered differentiation among islands. Based on the geographic history of the area and the divergence time, we considered the following hypotheses for the formation of the current distribution of Z. ishikawai: (1) single colonisation by a common ancestor and speciation into clades I and II and (2) independent colonisation of two differentiated ancestors. The former hypothesis seems unlikely, based on the estimated divergence time of clade I and clade II of > 20 million years ago (Fig. 8), when the southern part of the Ryukyu Archipelago was connected to Taiwan and the eastern margin of continental East Asia (Koba Reference Koba1992; Osozawa et al. Reference Osozawa, Shinjo, Armid, Watanabe, Horiguchi and Wakabayashi2012). In this geographical setting, speciation within the limited prospective area of the southern part of Ryukyu archipelago was unlikely even with the limited dispersal ability of harvestmen, thereby supporting the latter hypothesis. The ancient divergence between clade I and clade II suggests that each clade of Z. ishikawai had accumulated variations independently, leading to considerable morphological differentiation. We estimate that the colonising ancestors of both Z. ishikawai clades originated from continental Southeast Asia, because the genus Zepedanulus is also recorded in, for instance, northern Thailand and Malacca (Roewer Reference Roewer1927; Suzuki Reference Suzuki1981). However, it is unclear whether the ancestors of Z. ishikawai colonised the Ryukyu Archipelago region via Taiwan or invaded the archipelago directly from continental Asia. Systematic studies and phylogenetic analyses using a range of Epedanidae species are required to assess the history of diversification of Z. ishikawai and related species.

Our phylogenetic results suggest that Z. ishikawai did not frequently colonise islands in the southern part of the Ryukyu Archipelago after divergence of the two clades, because all monophyletic subgroups within each clade of this species comprised samples from an island, and individuals from the same island did not show paraphyletic or polyphyletic assemblies. Therefore, each subgroup of Z. ishikawai has an independent history of differentiation on each island, without additional introductions from other islands. Indeed, geographical variations among islands are also seen in other organisms in the Ryukyu Archipelago (Maki et al. Reference Maki, Yamashiro and Matsumura2003; Kiyoshi Reference Kiyoshi2008; Osozawa et al. Reference Osozawa, Su, Oba, Yagi, Watanabe and Wakabayashi2013, Reference Osozawa, Takáhashi and Wakabayashi2015; Hosokawa et al. Reference Hosokawa, Nikoh and Fukatsu2014; Osozawa et al. Reference Osozawa, Takáhashi and Wakabayashi2015). In the case of harvestmen, Kumekawa et al. (Reference Kumekawa, Ito, Miura, Yokoyama, Tebayashi, Arakawa and Fukuda2015) reported that K. insulanus in the central part of the Ryukyu Archipelago exhibited genetic differentiation between Amami-Oshima and Okinawa Islands based on the sequence of the CO1 in mitochondrial DNA. Although we detected morphological differences between clades of Z. ishikawai, whether differentiation occurred among subgroups of each clade is unclear because few individuals were collected from each island. Further analyses using larger numbers of samples from the islands are needed to clarify the morphological differences within each clade of Z. ishikawai.

Acknowledgements

The authors thank Drs K. Ito, Y. Minamiya, and H. Hayakawa for discussion in this study, and C. Uemoto, Y. Ozaki, K. Yoshioka, Y. Kamakura, Y. Mori, Y. Nibuno, T. Miyashita, and F. Maekawa for their field and laboratory assistance. The authors also thank Drs. N. Tsurusaki and T. Suguro for providing samples. This study was partly supported by a Grant-in-Aid for Scientific Research from the Ministry of Education, Science and Culture of Japan.

Footnotes

Subject editor: Marla Schwarzfeld

References

Arbogast, B.S. and Kenagy, G.J. 2001. Comparative phylogeography as an integrative approach to historical biogeography. Journal of Biogeography, 28: 819825. https://doi.org/10.1046/j.1365-2699.2001.00594.x.CrossRefGoogle Scholar
Atkins, H., Presto, J., and Cronk, Q.C.B. 2001. A molecular text of Huxley’s line: Cyrtandra (Gesneriaceae) in Borneo and the Philippines. Biological Journal of the Linnean Society, 72: 143159. https://doi.org/10.1006/bijl.2000.0500.Google Scholar
Avise, J.C., Shapiro, J.F., Daniel, S.W., Aquadro, C.F., and Lansman, R.A. 1983. Mitochondrial DNA differentiation during the speciation process in Peromyscus. Molecular Biology and Evolution, 1: 3856. doi: 10.1093/oxfordjournals.molbev.a040301 https://doi.org/10.1093/oxfordjournals.molbev.a040301.Google ScholarPubMed
Avise, J.C. 2000. Phylogeography: the history and formation of species. Harvard University Press, Cambridge, Massachusetts, United States of America.10.2307/j.ctv1nzfgj7CrossRefGoogle Scholar
Avise, J.C. 2009. Phylogeography: retrospect and prospect. Journal of Biogeography, 36: 315. https://doi.org/10.1111/j.1365-2699.2008.02032.x.CrossRefGoogle Scholar
Babcock, C.S., Heraty, J.M., De Barro, P.J., Driver, F., and Schmidt, S. 2001. Preliminary phylogeny of Encarsia Förster (Hymenoptera: Aphelinidae) based on morphology and 28S rDNA. Molecular Phylogenetics and Evolution, 18: 306323. https://doi.org/10.1006/mpev.2000.0875.CrossRefGoogle ScholarPubMed
Bittkau, C. and Comes, H.P. 2005. Evolutionary processes in a continental island system: molecular phylogeography of the Aegean Nigella arvensis alliance (Ranunculaceae) inferred from chloroplast DNA. Molecular Ecology, 14: 40654083. https://doi.org/10.1111/j.1365-294X.2005.02725.x.CrossRefGoogle Scholar
Bolnick, D.I. and Fitzpatrick, B.M. 2007. Sympatric speciation: models and empirical evidence. Annual Review of Ecology, Evolution and Systematics, 38: 459487.10.1146/annurev.ecolsys.38.091206.095804CrossRefGoogle Scholar
Bouckaert, R., Vaughan, T.G., Barido-Sottani, J., Duchêne, S., Fourment, M., and Gavryushkina, A. 2019. BEAST 2.5: an advanced software platform for Bayesian evolutionary analysis. PLOS Computational Biology, 15: e1006650.10.1371/journal.pcbi.1006650CrossRefGoogle ScholarPubMed
Burns, M., Hedin, M., and Shultz, J.W. 2012. Molecular phylogeny of the leiobunine harvestmen of eastern North America (Opiliones: Sclerosomatidae: Leiobuninae). Molecular Phylogenetics and Evolution, 63: 291298. https://doi.org/10.1016/j.ympev.2011.12.025.CrossRefGoogle Scholar
Che, Y., Wang, D., Shi, Y., Du, X., Zhao, Y., Lo, N., and Wang, Z. 2016. A global molecular phylogeny and timescale of evolution for Cryptocercus woodroaches . Molecular Phylogenetics and Evolution, 98: 201209. https://doi.org/10.1016/j.ympev.2016.02.005.CrossRefGoogle ScholarPubMed
Chiang, T.Y. and Schaal, B.A. 2006. Phylogeography of plants in Taiwan and the Ryukyu Archipelago. Taxon, 55: 3141. https://doi.org/10.2307/25065526.CrossRefGoogle Scholar
Clouse, R.M. and Wheeler, W.C. 2014. Descriptions of two new, cryptic species of Metasiro (Arachnida: Opiliones: Cyphophthalmi: Neogoveidae) from South Carolina, USA, including a discussion of mitochondrial mutation rates. Zootaxa, 3814: 177201. http://dx.doi.org/10.11646/zootaxa.3814.2.2.CrossRefGoogle Scholar
Coyne, J.A. and Orr, H.A. 2004. Speciation. Sinauer Associates Inc., Sunderland, Massachusetts, United States of America.Google Scholar
Cronn, R. and Wendel, J.F. 2004. Cryptic trysts, genomic mergers, and plant speciation. New Phytologist, 161: 133142. https://doi.org/10.1111/j.1469-8137.2004.00947.x.CrossRefGoogle Scholar
Cruz-López, J.A., Monjaraz-Ruedas, R., and Francke, O.F. 2019. Turning to the dark side: evolutionary history and molecular species delimitation of a troglomorphic lineage of armoured harvestman (Opiliones: Stygnopsidae). Arthropod Systematics and Phylogeny, 77: 285302. https://doi.org/10.26049/ASP77-2-2019-06.Google Scholar
Derkarabetian, S., Ledford, J., and Hedin, M. 2011. Genetic diversification without obvious genitalic morphological divergence in harvestmen (Opiliones, Laniatores, Sclerobunus robustus) from montane sky islands of western North America. Molecular Phylogenetics and Evolution, 63: 844853. https://doi.org/10.1016/j.ympev.2011.08.004.CrossRefGoogle Scholar
Derkarabetian, S., Steinmann, D.B., and Hedin, M. 2010. Repeated and time-correlated morphological convergence in cave-dwelling harvestmen (Opiliones, Laniatores) from Montane Western North America. PLOS One, 5: e10388. https://doi.org/10.1371/journal.pone.0010388.CrossRefGoogle ScholarPubMed
Didodemico, A. and Hedin, M. 2016. New species in the Sitalcina sura species group (Opiliones, Laniatores, Phalangodidae), with evidence for a biogeographic link between California desert canyons and Arizona sky islands. Zookeys, 586: 136. https://doi.org/10.3897/zookeys.586.7832.Google Scholar
Dobzhansky, T. 1937. Genetic and the origin of species. Columbia University Press, New York, New York, United States of America.Google Scholar
Emata, K.N. and Hedin, M. 2016. From the mountains to the coast and back again: ancient biogeography in a radiation of short-range endemic harvestmen from California. Molecular Phylogenetics and Evolution, 98: 233243. https://doi.org/10.1016/j.ympev.2016.02.002.CrossRefGoogle Scholar
Friedrich, M. and Tautz, D. 1997. Evolution and phylogeny of the Diptera: a molecular phylogenetic analysis using 28S rDNA sequences. Systematic Biology, 46: 674698. https://doi.org/10.2307/2413500.CrossRefGoogle ScholarPubMed
Fukuda, T., Song, I.J., Ito, T., Hayakawa, H., Minamiya, Y., Kanno, A., Hayakawa, H., Minamiya, Y., and Yokoyama, J. 2011. Phylogeography of Asparagus schoberioides Kunth (Asparagaceae) in Japan. American Journal of Plant Science, 2: 781789. https://doi.org/10.4236/ajps.2011.26093.CrossRefGoogle Scholar
Fukuda, T., Yokoyama, J., and Ohashi, H. 2001. Phylogeny and biogeography in Lycium (Solanaceae): inferences from chloroplast DNA sequences. Molecular Phylogenetics and Evolution, 19: 246258. https://doi.org/10.1006/mpev.2001.0921.CrossRefGoogle ScholarPubMed
Giribet, G., Rambla, M., Carranza, S., Baguna, J., Riutort, M., and Ribera, C. 1999. Phylogeny of the arachnid order Opiliones (Arthropoda) inferred from a combined approach of complete 18S and partial 28S ribosomal DNA sequences and morphology. Molecular Phylogenetics and Evolution, 11: 296307. https://doi.org/10.1006/mpev.1998.0583.CrossRefGoogle ScholarPubMed
Giribet, G., Vogt, L., Perez, G.A., Sharma, P., and Kury, A. 2010. A multilocus approach to harvestman (Arachnida: Opiliones) phylogeny with emphasis on biogeography and the systematics of Laniatores. Cladistics, 26: 408437. https://doi.org/10.1111/j.1096-0031.2009.00296.x.CrossRefGoogle Scholar
Gómez-Zurita, J., Juan, C., and Petitpierre, E. 2000. The evolutionary history of the genus Timarcha (Coleoptera, Chrysomelidae) inferred from mitochondrial CO1 gene and partial 16S rDNA sequences. Molecular Phylogenetics and Evolution, 14: 304317. https://doi.org/10.1006/mpev.1999.0712.CrossRefGoogle Scholar
Griffin, S.R. and Barrett, S.C.H. 2004. Post-glacial history of Trillium grandiflorum (Melanthiaceae) in eastern North America: inferences from phylogeography. American Journal of Botany, 91: 465473. https://doi.org/10.3732/ajb.91.3.465.CrossRefGoogle ScholarPubMed
Hansen, T.A. 1983. Modes of larval development and rates of speciation in early Tertiary neogastropods. Science, 220: 501502. https://doi.org/10.1126/science.220.4596.501.CrossRefGoogle ScholarPubMed
Hatusima, S. 1975. Flora of the Ryukyus, added and corrected. Okinawan Association of Biology Education, Naha, Japan. [Japanese].Google Scholar
Hebert, P.D., Cywinska, A., Ball, S.L., and deWaard, J.R. 2003. Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Series B, Biological Sciences, 270: 313321. https://doi.org/10.1098/rspb.2002.2218.CrossRefGoogle ScholarPubMed
Hedges, S.B., Hass, C.A., and Maxon, L.R. 1992. Caribbean biogeography: molecular evidence for dispersal in West Indian terrestrial vertebrates. Proceedings of the National Academy of Sciences, 89: 19091913. https://doi.org/10.1073/pnas.89.5.1909.CrossRefGoogle ScholarPubMed
Hedin, M. and Thomas, S.M. 2010. Molecular systematics of eastern North American Phalangodidae (Arachnida: Opiliones: Laniatores), demonstrating convergent morphological evolution in caves. Molecular Phylogenetics and Evolution, 54: 107121. https://doi.org/10.1016/j.ympev.2009.08.020.CrossRefGoogle Scholar
Hewitt, G.M. 2004. Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 359: 183195. https://doi.org/10.1098/rstb.2003.1388.CrossRefGoogle ScholarPubMed
Hirst, A.S. 1911. On some new Opiliones from Japan and the Loo-Choo Islands. Annals and Magazine of Natural History, 8: 625636.10.1080/00222931108693070CrossRefGoogle Scholar
Hisheh, S., Westerma, M., and Schmitt, L.H. 1998. Biogeography of the Indonesian archipelago: mitochondrial DNA variation in the fruit bat, Eonycteris spelaea . Biological Journal of the Linnean Society, 65: 329345. https://doi.org/10.1006/bijl.1998.0244.Google Scholar
Hollingsworth, M.L., Bailey, J.P., Hollingsworth, P.M., and Ferris, C. 1999. Chloroplast DNA variation and hybridisation between invasive populations of Japanese knotweed and giant knotweed (Fallopia, Polygonaceae). Botanical Journal of the Linnean Society, 129: 139154.10.1111/j.1095-8339.1999.tb00496.xCrossRefGoogle Scholar
Hosokawa, T., Nikoh, N., and Fukatsu, T. 2014. Fine-scale geographical origin of an insect pest invading North America. PLOS One, 9: e89107. https://doi.org/10.1371/journal.pone.0089107.CrossRefGoogle ScholarPubMed
Ikeda, H., Nishikawa, M., and Sota, T. 2011. Loss of flight promote beetle diversification. Nature Communications, 3: 648. https://doi.org/10.1038/ncomms1659.CrossRefGoogle Scholar
Jablonski, D. 1986. Larval ecology and macroevolution in marine invertebrates. Bulletin of Marine Science, 39: 565587.Google Scholar
Jinbo, U., Kato, T., and Ito, M. 2011. Current progress in DNA barcoding and future implications for entomology. Entomological Science, 14: 107124. https://doi.org/10.1111/j.1479-8298.2011.00449.x.CrossRefGoogle Scholar
Kaito, T. and Toda, M. 2016. The biogeographical history of Asian keelback snakes of the genus Hebius (Squamata: Colubridae: Natricinae) in the Ryukyu Archipelago, Japan. Biological Journal of the Linnean Society, 118: 187199. https://doi.org/10.1111/bij.12726.CrossRefGoogle Scholar
Kawana, T. 2002. Neotectonics in the Ryukyu Arc. In Establishment of the Ryukyu Arc and colonization of biota. Edited by Kimura, M.. Okinawa Times, Naha, Okinawa, Japan. Pp. 5983. [Japanese].Google Scholar
Keally, C.T. 2005. Japanese Pleistocene landbridges and the earliest watercraft [online]. Available from http://www.t-net.ne.jp/˜keally/MiddlePalaeol/landbridges.html [accessed 15 June 2019].Google Scholar
Kimura, M. 1996. Quaternary paleogeography of the Ryukyu Arc. Journal of Geography, 105: 259285. [Japanese]. https://doi.org/10.5026/jgeography.105.3_259.CrossRefGoogle Scholar
Kimura, M. 2000. Paleogeography of the Ryukyu Islands. Tropics, 10: 524. https://doi.org/10.3759/tropics.10.5.CrossRefGoogle Scholar
Kimura, M. 2002. Formation and paleogeography of the Ryukyu Arc. In Establishment of the Ryukyu Arc and colonization of biota. Edited by Kimura, M.. Okinawa Times, Naha, Japan. Pp. 1954. [Japanese].Google Scholar
Kiyoshi, T. 2008. Differentiation of golden-ringed dragonfly Anotogaster sieboldii (Selys, 1854) (Cordulegastridae: Odonata) in the insular east Asia revealed by the mitochondrial gene genealogy with taxonomic implications. Journal of Zoological Systematics and Evolutionary Research, 46: 105109. https://doi.org/10.1111/j.1439-0469.2007.00445.x.CrossRefGoogle Scholar
Kizaki, K. and Oshiro, I. 1977. Palaeogeography of the Ryukyu Islands. Marine Science Monthly, 9: 542549. [In Japanese with English summary].Google Scholar
Kizaki, K. and Oshiro, I. 1980. The origin of the Ryukyu Islands. In Natural history of the Ryukyus. Edited by Kizaki, K.. Tsukiji-shokan, Tokyo, Japan [Japanese].Google Scholar
Koba, M. 1992. Influx of Kuroshio current into the Okinawa trough and inauguration of Quaternary coral-reef building in the Ryukyu Island Arc, Japan. Quaternary Research, 31: 359373. https://doi.org/10.4116/jaqua.31.359.CrossRefGoogle Scholar
Kumekawa, Y., Ito, K., Miura, O., Yokoyama, J., Tebayashi, S., Arakawa, R., and Fukuda, T. 2015. Molecular phylogeny of Kilungius insulanus (Arachnida: Opiliones: Epedanidae) in Amami-Oshima Island and Okinawa Island. Edaphologia, 96: 17. https://doi.org/10.20695/edaphologia.96.0_1.Google Scholar
Kumekawa, Y., Ito, K., Tsurusaki, N., Hayakawa, H., Ohga, K., Yokoyama, J., et al. 2014. Morphological and molecular analyses of Pseudobiantes japonicus Hirst, 1911 and its allied species. Annals of the Entomological Society of America, 107: 756772.10.1603/AN13107CrossRefGoogle Scholar
Kumekawa, Y., Miura, O., Fujimoto, H., Ito, K., Arakawa, R., Yokoyama, J., and Fukuda, T. 2019. Possible occurrence of reproductive isolation between two geographical clades of a laniatorid harvestman Pseudobiantes japonicus (Arachnida: Opiliones: Epedanidae) in Shikoku. Edaphologia, 104: 1924.Google Scholar
Kumekawa, Y., Murjoko, A., Hayakawa, H., Ohga, K., Mori, M., Miyazaki, A., et al. 2013. Molecular analyses of folk varieties of the sago palm (Metroxylon sagu Rottb) using the internal transcribed spacer (ITS) region and nuclear microsatellite DNA. Sago Palm, 21: 1419.Google Scholar
Maddison, W.P. and Knowles, L.L. 2006. Inferring phylogeny despite incomplete lineage sorting. Systematic Biology, 55: 2130. https://doi.org/10.1080/10635150500354928.CrossRefGoogle ScholarPubMed
Maekawa, K., Lo, N., Kitade, O., Miura, T., and Matsumoto, T. 1999. Molecular phylogeny and geographic distribution of wood-feeding cockroaches in East Asian islands. Molecular Phylogenetics and Evolution, 13: 360376. https://doi.org/10.1006/mpev.1999.0647.CrossRefGoogle ScholarPubMed
Maki, M., Yamashiro, T., and Matsumura, S. 2003. High levels of genetic diversity in island populations of the island endemic Suzukia luchuensis (Labiatae). Heredity, 91: 300306. https://doi.org/10.1038/sj.hdy.6800311.CrossRefGoogle Scholar
Mallatt, J. and Sullivan, J. 1998. 28S and 18S rDNA sequences support the monophyly of lampreys and hagfishes. Molecular Biology and Evolution, 15: 17061718.10.1093/oxfordjournals.molbev.a025897CrossRefGoogle Scholar
Matsumura, S., Yokoyama, J., Fukuda, T., and Maki, M. 2009. Intraspecific differentiation of Limonium wrightii (Plumbaginaceae) on northwestern Pacific Islands: rate heterogeneity in nuclear rDNA and its distance-independent geographic structure. Molecular Phylogenetics and Evolution, 53: 10321036. https://doi.org/10.1016/j.ympev.2009.06.011.CrossRefGoogle ScholarPubMed
Maura, M., Salvi, D., Bologna, M.A., Nascetti, G., and Canestrelli, D. 2014. Northern richness and cryptic refugia: phylogeography of the Italian smooth newt Lissotriton vulgaris meridionalis . Biological Journal of the Linnean Society, 113: 590603. https://doi.org/10.1111/bij.12360.CrossRefGoogle Scholar
Mayr, E. 1963. Animal Species and Evolution. Harvard University Press, Cambridge, Massachusetts, United States of America.10.4159/harvard.9780674865327CrossRefGoogle Scholar
Mezzasalma, M., Dall’Asta, A., Loy, A., Cheylan, M., Lymberakis, P., Zuffi, M.A.L., et al. 2015. A sister’s story: comparative phylogeography and taxonomy of Hierophis viridiflavus and H. gemonensis (Serpentes, Colubridae). Zoologica Scripta, 44: 495508. https://doi.org/10.1111/zsc.12115.CrossRefGoogle Scholar
Minamiya, Y., Yokoyama, J., and Fukuda, T. 2009. A phylogeographic study of the Japanese earthworm, Metaphire sieboldi (Horst, 1883) (Oligochaeta: Megascolecidae): inferences from mitochondrial DNA sequences. European Journal of Soil Biology, 45: 423430. https://doi.org/10.1016/j.ejsobi.2009.06.004.CrossRefGoogle Scholar
Miura, O., Torchin, M.E., Bermingham, E., Jacobs, D.K., and Hechinger, R.F. 2012. Flying shells: historical dispersal of marine snails across Central America. Proceedings of the Royal Society B: Biological Sciences, 279: 10611067. https://doi.org/10.1098/rspb.2011.1599.CrossRefGoogle ScholarPubMed
Mráz, P., Chrtek, J., Fehrer, J., and Plačková, I. 2005. Rare recent natural hybridisation in Hieracium s. str.: evidence from morphology, allozymes, and chloroplast DNA. Plant Systematics and Evolution, 255: 177192.10.1007/s00606-005-0329-1CrossRefGoogle Scholar
Nakamura, K., Denda, T., Kokubugata, G., Huang, C., Peng, C.I., and Yokota, M. 2015. Phylogeny and biogeography of the Viola iwagawae-tashiroi species complex (Violaceae, section Plagiostigma) endemic to the Ryukyu archipelago. Plant Systematics and Evolution, 301: 337351. https://doi.org/10.1007/s00606-014-1076-y.CrossRefGoogle Scholar
Nei, M. and Kumar, S. 2000. Molecular Evolution and Phylogenetics. Oxford University Press, New York, United States of America.Google Scholar
Nijman, I.J., Bradley, D.G., Hanotte, O., Otsen, M., and Lenstra, J.A. 1999. Satellite DNA polymorphisms and AFLP correlate with Bos indicus-taurus hybridization. Animal Genetics, 30: 265273.10.1046/j.1365-2052.1999.00475.xCrossRefGoogle ScholarPubMed
Nunn, G.B., Theisen, B.F., Christensen, B., and Arctander, P. 1996. Simplicity-correlated size growth of the nuclear 28S ribosomal RNA D3 expansion segment in the crustacean order Isopoda. Journal of Molecular Evolution, 42: 211223. https://doi.org/10.1007/BF02198847.CrossRefGoogle ScholarPubMed
Nylander, J.A.A. 2004. MrModeltest 2.3: programme distributed by the author. Evolutionary Biology Centre, Uppsala University, Uppsala, Sweden.Google Scholar
Ohdachi, S., Dokuchaev, N.E., Hasegawa, M., and Masuda, R., 2001. Intraspecific phylogeny and geographical variation of six species of northeastern Asiatic Sorex shrews based on the mitochondrial cytochrome b sequences. Molecular Ecology, 10: 21992213. https://doi.org/10.1046/j.1365-294X.2001.01359.x.CrossRefGoogle ScholarPubMed
Ohshiro, I. 2003. Geology and paleontology of Ryukyu Islands. In The flora and fauna of inland waters in the Ryukyu islands. Edited by Nishijima, S., Nishida, M., Shikatani, N., and Shokita, S.. Tokai University Press, Hiratsuka, Kanagawa, Japan. Pp. 1116. [Japanese].Google Scholar
Oleksa, A., Chybicki, I., Tofilski, A., and Burczyk, J. 2011. Nuclear and mitochondrial patterns of introgression into native dark bees (Apis mellifera mellifera) in Poland. Journal of Apicultural Research, 50: 116129.10.3896/IBRA.1.50.2.03CrossRefGoogle Scholar
Osozawa, S., Pavlis, T., and Flower, M. 2011. Sedimentary block-in-matrix fabric affected by tectonic shear, Miocene Nabae complex, Japan. In Melanges: processes of formation and societal significance. Edited by Wakabayashi, J. and Dilek, Y.. Special Paper of Geological Society of America, 480:189–206. https://doi.org/10.1130/2011.2480(08).Google Scholar
Osozawa, S., Shinjo, R., Armid, A., Watanabe, Y., Horiguchi, T., and Wakabayashi, J. 2012. Palaeogeographic reconstruction of the 1.55 Ma synchronous isolation of the Ryukyu Islands, Japan and Taiwan and inflow of the Kuroshio warm current. International Geology Review, 54: 13691388. https://doi.org/10.1080/00206814.2011.639954.CrossRefGoogle Scholar
Osozawa, S., Su, Z.H., Oba, Y., Yagi, T., Watanabe, Y., and Wakabayashi, J. 2013. Vicariant speciation due to 1.55 Ma isolation of the Ryukyu islands, Japan, based on geological and GenBank data. Entomological Science, 16:267277. https://doi.org/10.1111/ens.12037.CrossRefGoogle Scholar
Osozawa, S., Takáhashi, M., and Wakabayashi, J. 2015. Ryukyu endemic Mycalesis butterflies, speciated vicariantly due to isolation of the islands since 1.55 Ma. Lepidoptera Science, 66: 814. https://doi.org/10.18984/lepid.66.1_8.Google Scholar
Ota, H. 1998. Geographic patterns of endemism and speciation in amphibians and reptiles of the Ryukyu Archipelago, Japan, with special reference to their paleogeographical implications. Researches on Population Ecology, 40: 189204. https://doi.org/10.1007/BF02763404.CrossRefGoogle Scholar
Ota, H. 2012. Phylogeography of terrestrial biota and geohistory of the Ryukyu Archipelago, Japan: hypothetical scenarios hitherto proposed and directions of future studies. Chikyu Monthly, 34: 427436.Google Scholar
Percy, D.M., Page, R.D., and Cronk, Q.C. 2004. Plant–insect interactions: double-dating associated insect and plant lineages reveals asynchronous radiations. Systematic Biology, 53: 120127. https://doi.org/10.1080/10635150490264996.CrossRefGoogle ScholarPubMed
Pons, J. and Vogler, A.P. 2005. Complex pattern of coalescence and fast evolution of a mitochondrial rRNA pseudogene in a recent radiation of tiger beetles. Molecular Biology and Evolution, 22: 9911000. https://doi.org/10.1093/molbev/msi085.CrossRefGoogle Scholar
Poulakakis, N., Lymberakis, P., Antonioub, A., Chalkia, D., Zourosb, E., Mylonasa, M., and Valakos, E. 2003. Molecular phylogeny and biogeography of the wall-lizard Podarcis erhardii (Squamata: Lacertidae). Molecular Phylogenetics and Evolution, 28: 3846.10.1016/S1055-7903(03)00037-XCrossRefGoogle Scholar
Rambaut, A. 2014. FIGTREE. Version 1.4.2 [online]. Available from http://tree.bio.ed.ac.uk/software/figtree/ [accessed 25 June 2018].Google Scholar
Rambaut, A., Suchard, M.A., Xie, D., and Drummond, A.J. 2014. Tracer v1.6 [online]. Available from http://beast.bio.ed.ac.uk/Tracer [accessed 25 June 2018].Google Scholar
Rieseberg, L.H. and Soltis, D.E. 1991. Phylogenetic consequences of cytoplasmic gene flow in plants. Evolutionary Trends in Plants, 5: 6584.Google Scholar
Roewer, C.F. 1927. Weitere Weberknechte I. (1. Erganzung der : “Weberknechte der Erde”, 1923) Abhandlungen der Naturwissenschaftlichen Verein zu Bremen, Bremen 26: 261–402.Google Scholar
Ronquist, F. and Huelsenbeck, J.P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. BioInformatics, 19: 15721574. https://doi.org/10.1093/bioinformatics/btg180.CrossRefGoogle ScholarPubMed
Schlick-Steiner, B.C., Steiner, F.M., Seifert, B., Stauffer, C., Christian, E., and Crozier, T.H. 2010. Integrative taxonomy: a multisource approach to exploring biodiversity. Annual Review of Entomology, 55: 421438. https://doi.org/10.1146/annurev-ento-112408-085432.CrossRefGoogle ScholarPubMed
Schmidt, S., Driver, F., and De Barro, P. 2006. The phylogenetic characteristics of three different 28S rRNA gene regions in Encarsia (Insecta, Hymenoptera, Aphelinidae). Organisms Diversity and Evolution, 6: 127139. https://doi.org/10.1016/j.ode.2005.07.002.CrossRefGoogle Scholar
Seki, S., Sakanashi, M., Kawaji, N., and Kotaka, N. 2007. Phylogeography of the Ryukyu robin (Erithacus komadori): population subdivision in land-bridge islands in relation to the shift in migratory habit. Molecular Ecology, 16: 101113. https://doi.org/10.1111/j.1365-294X.2006.03117.x.CrossRefGoogle ScholarPubMed
Seo, A., Watabanabe, M., Hotta, M., and Murakami, N. 2004. Geographical patterns of allozyme variation in Angelica japonica (Umbelliferae) and Farfugium japonicum (Compositae) on the Ryukyu Islands, Japan. Acta Phytotaxonomica et Geobotanica, 55: 2944. https://doi.org/10.18942/apg.KJ00004622802.Google Scholar
Sharma, P. and Giribet, G. 2009. Sandokanid phylogeny based on eight molecular markers: the evolution of a Southeast Asian endemic family of Laniatores (Arachnida, Opiliones). Molecular Phylogenetics and Evolution, 52: 432447. https://doi.org/10.1016/j.ympev.2009.03.013.CrossRefGoogle Scholar
Sharma, P.P. 2010. New Australasian Zalmoxidae (Opiliones: Laniatores) and a new case of male polymorphism in Opiliones. Zootaxa, 3236: 135. http://dx.doi.org/10.11646/zootaxa.3236.1.1.CrossRefGoogle Scholar
Sharma, P.P. and Giribet, G. 2011. The evolutionary and biogeographic history of the armoured harvestmen: Laniatores phylogeny based on 10 molecular markers, with the description of two new families of Opiliones (Arachnida). Invertebrate Systematics, 25: 106142. https://doi.org/10.1071/IS11002.CrossRefGoogle Scholar
Sharma, P.P., Buenavente, P.A.C., Clouse, R.M., Diesmos, A.C., and Giribet, G. 2012. Forgotten gods: Zalmoxidae of the Philippines and Borneo (Opiliones: Laniatores). Zootaxa, 3280: 2955. http://dx.doi.org/10.11646/zootaxa.3280.1.2.CrossRefGoogle Scholar
Shinjo, R. 2015. The geology of Okinawa-Jima Island. In Nature in the Ryukyu Archipelago: coral reefs, biodiversity, and the natural environment. Edited by Fujita, K., Arakaki, T., Denda, T., Hidaka, M., Hirose, E., and Reimer, J.D.. Faculty of Science, University of the Ryukyus, Nishimura, Okinawa, Japan. Pp. 12119.Google Scholar
Simsek, M., Tanira, M.O.M., Al-Baloushi, K.A., Al-Barwani, H.S., Lawatia, K.M., and Bayoumi, R.A.L. 2001. A precaution in the detection of heterozygotes by sequencing: comparison of automated DNA sequencing and PCR–restriction fragment length polymorphism methods. Clinical Chemistry, 47: 134137. https://doi.org/10.1093/clinchem/47.1.134.CrossRefGoogle ScholarPubMed
Smith, C.I. and Farrell, B.D. 2006. Evolutionary consequences of dispersal ability in cactus-feeding insects. Genetica, 126: 323334. https://doi.org/10.1007/s10709-005-0714-2.CrossRefGoogle ScholarPubMed
Suzuki, S. 1971. Opiliones of the Ryukyus. Journal of Science of the Hiroshima University, Series B, Division 1 (Zoology), 23: 187213.Google Scholar
Suzuki, S. 1973. Opiliones from the south-west islands. Journal of Science of the Hiroshima University, Series B, Division 1 (Zoology), 24: 205279.Google Scholar
Suzuki, S. 1981. Three opilionids from Thailand. Annotationes Zoologicae Japonenses, 54: 267272.Google Scholar
Swofford, D.L. 2002. PAUP 4.0b10: phylogenetic analysis using parsimony Sinauer Associates, Sunderland, Massachusetts, United States of America.Google Scholar
Tamura, K., Stecher, G., Peterson, N., Filipski, A., and Kumar, S. 2013. MEGA6: molecular evolutionary genetics analysis, version 6.0. Molecular Biology and Evolution, 30: 27252729. https://doi.org/10.1093/molbev/mst197.CrossRefGoogle Scholar
Taylor, M.L. and Rogers, A.D. 2015. Evolutionary dynamics of a common sub-Antarctic octocoral family. Molecular Phylogenetics and Evolution, 84: 185204. https://doi.org/10.1016/j.ympev.2014.11.008.CrossRefGoogle ScholarPubMed
Thomas, S.M. and Hedin, M. 2008. Multigenic phylogeographic divergence in the paleoendemic southern Appalachian opilionid Fumontana deprehendor Shear (Opiliones, Laniatores, Triaenonychidae). Molecular Biology and Evolution, 46: 645658. https://doi.org/10.1016/j.ympev.2007.10.013.Google Scholar
Thompson, J.D., Higgins, D.G., and Gibson, T.J. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Research, 22: 46734680. https://doi.org/10.1093/nar/22.22.4673.CrossRefGoogle ScholarPubMed
Toews, D.P.L. and Brelsford, A. 2012. The biogeography of mitochondrial and nuclear discordance in animals. Molecular Ecology, 21: 39073930. https://doi.org/10.1111/j.1365-294X.2012.05664.x.CrossRefGoogle ScholarPubMed
Tsurusaki, N. 2006a. Geographic variation of chromosome and somatic morphology in the Japanese polymorphic species Leiobunum hiraiwai (Arachnida: Opiliones: Sclerosomatidae). Zootaxa, 1325: 157190. http://dx.doi.org/10.11646/zootaxa.1325.1.10.CrossRefGoogle Scholar
Tsurusaki, N. 2006b. Biogeography of Opiliones in the Ryukyu Islands. Nature and Insects, 41: 3033.Google Scholar
Tsurusaki, N. and Suzuki, S. 2015. Pictorial keys to soil animals of Japan. Second edition. Edited by Aoki, J.. Tokai University Press, Hadano, Japan. Pp. 121145. [Japanese].Google Scholar
Ujiie, H. 1990. Geological history of the Ryukyu Island Arc. In Nature of Okinawa geomorphology and geology. Edited by Ujiie, H.. Hirugisha, Naha, Okinawa, Japan. Pp. 251255. [Japanese].Google Scholar
Vaini, J.O., Grisolia, A.B., Prado, F.D., and Porto-Foresti, F. 2014. Genetic identification of interspecific hybrid of Neotropical catfish species (Pseudoplatystoma corruscans vs. Pseudoplatystoma reticulatum) in rivers of Mato Grosso do Sul State Brazil. Neotropical Ichthyology, 12: 635641.10.1590/1982-0224-20130169CrossRefGoogle Scholar
Velo-Antón, G., Zamudio, K.R., and Cordero-Rivera, A. 2012. Genetic drift and rapid evolution of viviparity in insular fire salamanders (Salamandra salamandra). Heredity, 108: 410418. https://doi.org/10.1038/hdy.2011.91.CrossRefGoogle Scholar
Vink, C.J., Thomas, S.M., Paquin, P., Hayashi, C.Y., and Hedin, M. 2005. The effects of preservatives and temperatures on arachnid DNA. Invertebrate Systematics, 19: 99104. https://doi.org/10.1071/IS04039.CrossRefGoogle Scholar
Wendel, J.F. and Doyle, J.J. 1998. Phylogenetic incongruence: windows into genome history and molecular evolution. In Molecular systematics of plants II DNA sequencing. Edited by Soltis, D.E., Soltis, P.S., and Doyle, J.J.. Kluwer, Dordrecht, The Netherlands. Pp. 265296. https://doi.org/10.1007/978-1-4615-5419-6_10.CrossRefGoogle Scholar
Xiong, B. and Kocher, T.D. 1991. Comparison of mitochondrial DNA sequences of seven morphospecies of black flies (Diptera: Simuliidae). Genome, 34: 306311. https://doi.org/10.1139/g91-050.CrossRefGoogle Scholar
Zera, A.J. 1981. Genetic structure of two species of waterstriders (Gerridae: Hemiptera) with differing degrees of winglessness. Evolution, 35: 218225. https://doi.org/10.2307/2407832.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Distribution of Z. ishikawai in the Ryukyu Archipelago, Japan.

Figure 1

Table 1. Summary of Zepedanulus ishikawai analysed in this study. Specimen number corresponds to our phylogenetic results (Figs. 3, 4, 5, and 6).

Figure 2

Fig. 2. (A) Spined type and (B) nonspined type, and characters of Zepedanulus ishikawai (collected from Yonaguni Island): A, BL – body length; B, CL – carapace length; C, CW – carapace width; D, Fe1L – femur length, first leg; E, Fe2L – femur length, second leg; F, Fe3L – femur length, third leg; G, Fe4L – femur length, fourth leg; H, SL – spine length; I, Chel.L – chelicerae length; and J, PTL – palpal–tarsus length, and the morphology of the second scutal area. Arrowhead indicates the spine. Scale bar = 1 mm.

Figure 3

Fig. 3. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by the maximum-likelihood method based on the sequences of CO1 and 16S rRNA in mtDNA. Numbers above or below branches indicate bootstrap values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Figure 4

Fig. 4. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by Bayesian inference based on the sequences of CO1 and 16S rRNA in mtDNA. Numbers above or below branches indicate posterior probabilities. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Figure 5

Fig. 5. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by the maximum-likelihood method based on the sequence of 28S rRNA in nrDNA. Numbers above or below branches indicate the bootstrap values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Figure 6

Fig. 6. Phylogenetic tree of Zepedanulus ishikawai and related species with outgroups by Bayesian inference based on the sequence of 28S rRNA in nrDNA. Numbers above or below branches indicate the posterior probabilities. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Figure 7

Fig. 7. A, Restriction sites in the 28S rRNA sequence determined by polymerase chain reaction–restriction fragment length polymorphism (M, HaeIII site) and B, polymerase chain reaction–restriction fragment length polymorphism profile of clades I and II of Zepedanulus ishikawai (M, size markers).

Figure 8

Fig. 8. Estimated divergence time (over millions of years ago; MYA) for the two clades of Zepedanulus ishikawai. Horizontal bars represent the upper and lower 95% of highest posterior density values. Symbols indicate the origin of samples: ○, Ishigaki Island; △, Iriomote Island; □, Yonaguni Island.

Figure 9

Table 2. The measurement (average ± standard deviation) of Zepedanulus ishikawai.

Figure 10

Fig. 9. Histograms of the morphological characteristics of Zepedanulus ishikawai in Iriomote Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Figure 11

Fig. 10. Histograms of the morphological characteristics of Zepedanulus ishikawai in Ishigaki Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Figure 12

Fig. 11. Histograms of the morphological characteristics of Zepedanulus ishikawai in Yonaguni Island. White and grey colours indicate clade I and clade II, respectively. Abbreviations: Fe1L, femur length, first leg; Fe2L, femur length, second leg; Fe3L, femur length, third leg; Fe4L, femur length, fourth leg; BL, body length; CL, carapace length; CW, carapace width; SL, spine length; Chel.L, chelicerae length; PTL, palpal–tarsus length.

Figure 13

Fig. 12. Distribution of morphological characteristics in a simplified phylogenetic tree of Zepedanulus ishikawai and related species.