The genetic variation in Eumeces latiscutatus from the main islands of Japan and the northern Ryukyus was investigated and compared with that of its close relatives (E. okadae and E. marginatus) using allozyme data. All three samples from the northern Tokara Island Group, currently identified as E. latiscutatus, were shown to belong to E. marginatus and not to E. latiscutatus. The non-monophyly of the northern Tokara samples and the great genetic differences within these samples may have resulted from colonization from more than one source population by northeastward overseas dispersal. The samples from the Izu Peninsula were genetically much closer to E. okadae than the other conspecific ones. This indicates that the samples from the Izu Peninsula and the other samples of E. latiscutatus should be treated as distinct species, and that E. latiscutatus from the Izu Peninsula and E. okadae from the Izu Island Group may be treated as conspecific. Samples from western Japan were genetically well differentiated from those of eastern Japan. Within the western group, the samples from the Osumi Island Group was genetically distinct from those from the other regions, by possessing unique alleles. Our phenograms also reveal a distant affinity between samples from the Danjo Island Group and the main islands of Japan. This may be the result of long geographic isolation of the Osumi and Danjo Island Groups from Kyushu. By contrast, samples from Sapporo and Aomori were poorly differentiated genetically in spite of the long separation of these two localities by the Tsugaru Strait. This suggests that overseas dispersal of E. latiscutatus occurred across this strait after its formation.
The Japanese five-lined skink, Eumeces latiscutatus (Hallowell, 1861), is a moderate-sized scincine lizard that occurs on the main islands of Japan (Hokkaido, Honshu, Shikoku, Kyushu), the adjacent islands, the northern Ryukyus, and the coastal region of Primorsky Kray, eastern Russia. Within East Asian Eumeces, this species is considered most closely related to E. okadae, which is distributed in the Izu Island Group (Hikida, 1993; Kato et al., 1994).
Hikida et al. (1992) reported that specimens of E. latiscutatus collected from the northern Tokara Island Group in the northern Ryukyus differed from those from the main islands of Japan in scale row number. They also noted the presence of remarkable differentiation in a few external characters among insular populations within this region. Based on allozyme data, Kato et al. (1994) demonstrated that the Kuchinoshima sample in the northern Tokara Island Group is genetically closer to E. marginatus than to E. latiscutatus from the main islands of Japan. However, no genetic data are available for the remaining samples from this area.
On the main islands of Japan, E. latiscutatus shows variation in scutellation, such as in the postnasal condition and scale row number, even within a population (Stejneger, 1907; Taylor, 1936), but there have been no attempts to study the geographic variation of the morphological characters in this species. Kato et al. (1994) found only slight genetic differentiation between two samples, Kyoto and Kagoshima, from western Japan. However, recent genetic studies revealed substantial east-west divergences in other non-volant terrestrial species that show little morphological variation (Saitoh et al., 1989; Matsui et al., 2000; Harada et al., 2001). Considering the wide distribution of E. latiscutatus, more comprehensive analyses of the genetic variation in this species are required.
This study analyzed genetic variation and differentiation in E. latiscutatus based on allozyme data for a number of samples from throughout its range, and discusses its evolution and taxonomic implications.
MATERIALS AND METHODS
A total of 303 skinks from 20 samples of E. latiscutatus (Sample 1–20), a sample of E. okadae (Sample 21) and two samples of E. marginatus (Sample 22 and 23) were examined (Fig. 1 and Appendix).
The liver, skeletal muscles and heart were removed from each specimen, and their extracts were subjected to a horizontal starch-gel electrophoretic assay. Sample preparation and electrophoretic procedures followed Kato et al. (1994). Voucher specimens were deposited in the collection of the Department of Zoology, Kyoto University (KUZ: see Appendix).
The enzymes examined and their Enzyme Commission (E.C.) numbers, presumptive loci, tissues used, and the buffer systems employed in the electrophoretic assay are given in Table 1. Enzyme nomenclature, E.C. numbers, and abbreviations of enzymes follow Murphy et al. (1996). Electromorphs were designated alphabetically in order of their anodal mobility.
Enzymes, presumptive loci, tissues, and buffer systems used in the analysis of the variation among populations of E. latiscutatus, E. okadae, and E. marginatus.
All loci were evaluated for genetic polymorphism (mean number of alleles per locus, percentage of polymorphic loci, and mean heterozygosity by direct count) and for genetic structuring using Wright's (1965) F-statistics (Fst). In order to estimate the overall genetic differentiation among samples, Nei's (1978) unbiased genetic distance and the modified Rogers' distance (Wright, 1978) were calculated using allele frequencies. Of these, Nei's distance is considered to be well correlated with the divergence time of two given samples, although it is non-metric (Nei, 1987). Estimations of the genetic relationships among samples were obtained using the UPGMA algorithm (Sneath and Sokal, 1973) based on Nei's distance, and the neighbor-joining (NJ) procedure (Saitou and Nei, 1987) based on the modified Rogers' distance. The UPGMA method assumes equal rates of molecular evolution along all branches (Nei, 1987), whereas the NJ analysis does not need such a presumption. The NJ phenogram was rooted at the mid-point of the longest pathway. NJ analysis was applied to the data set using the program included with PHYLIP (Felsenstein, 1993). The remaining computations were made with BIOSYS-1 (Swofford and Selander, 1981). In addition, we conducted a principal component analysis (PCA) using the PRINCOMP procedure (SAS Inst. Inc., 1990) to detect regional patterns of geographical variation. PCA was applied to a correlation matrix based on alleles occurring at frequencies greater than 0.05 in one or more samples. The number of alleles (k) analyzed at a locus was k - 1; therefore, we arbitrarily dropped the most common alleles, for a total of 31 alleles from 12 loci. To render the scale of variation more linear, we used an arcsine-square-root transformation for the allele frequencies (Sokal and Rohlf, 1981). A two-dimensional projection of the samples on the first two PC axes was made to illustrate relationships among the samples.
Of the 20 presumptive loci examined, four (Aat-2, Fumh, Gtdh, and Pgm-1) were fixed for the same allele in all samples, and the remaining 16 loci were polymorphic in at least one sample (Table 2). At four of the 16 variable loci, a single allele was predominant (frequencies of 50% or more) in all samples (Idh-1, Mdh-2, Mpi, and Sod-1); the predominant allele for the remaining loci differed among the 23 samples.
Allele frequencies at 16 polymorphic loci, mean number of alleles per locus (A), percentage of polymorphic loci (P), and mean hetero zygosity by direct count (H) in E. latiscutatus (Sample 1-20), E. okadae (Sample 21), and E. marginatus (Sample 22 and 23). Sample numbers correspond to those in Fig. 1. Sample size is given in parentheses.
The samples from the northern Tokara Island Group (Sample 18–20), while showing complete allelic displacement from all other samples of E. latiscutatus (Sample 1–17) at three loci (Ck-2, Est-4, and Gpi), exhibited no such displacements at all with the E. marginatus samples (Sample 22 and 23). Among the samples from the northern Tokara Island Group, remarkable allele frequency differences were observed at five (Aat-1, Acoh-1, Ada, Pgdh, and Pgm-2) of the nine polymorphic loci. Based on the pattern of shared alleles at the Idh-2 locus, all samples except for the northern Tokara samples of E. latiscutatus were clearly divided into three groups: an Izu Peninsula group consisting of two samples from the Izu Peninsula (Sample 4 and 5: Fig. 1) possessing Idh-2c, an eastern group consisting of six samples from the eastern parts of the main islands (Sample 1–3, 6–8) possessing Idh-2b, a western group consisting of nine samples from the western parts of the main islands, including the Danjo and Osumi Island Groups (Sample 9–17) possessing Idh-2d. Samples of the Izu Peninsula group, while showing complete allelic displacement from those of the eastern and western groups at two loci (Ada and Idh-2), showed no such displacement from the sample of E. okadae (Sample 21). At Est-4, obvious allele frequency differences were observed between samples from the Osumi Island Group (Sample 15–17) and the other region (Sample 9–14) within the western group. Moreover, the former had a unique allele (Idh-1c) and two alleles (Pgdhc and Pgm-2d) at noticeable frequencies at highly polymorphic loci, which were absent or rare in the latter.
Among the samples of E. latiscutatus, the mean number of alleles per locus (A) varied from 1.0 to 1.5, the percentage of polymorphic loci (P) from 0.0 to 35.0, and the mean heterozygosity (H) values from 0.000 to 0.116 (Table 2). The greatest heterozygosity was found for the Yakushima sample. The mean Fst calculated for all samples of E. latiscutatus, for 17 samples excluding the northern Tokara samples, and 15 samples from eastern and western Japan was 0.806, 0.740 and 0.700, respectively.
Pairwise comparisons of Nei's (1978) and the modified Rogers' (Wright, 1978) genetic distance coefficients, calculated over all 20 loci, are presented in Table 3. Within E. latiscutatus, high values of Nei's genetic distance were found, ranging from 0.001 to 0.799. This is because the northern Tokara samples were genetically highly differentiated from all other samples (D=0.301–0.799). Compared to these values, the values obtained between the northern Tokara samples and E. marginatus were relatively small (0.062–0.295). When the three samples from the northern Tokara Island Group were excluded, the D values among 17 samples of E. latiscutatus ranged from 0.001 to 0.381, with the highest D value between Sample 1 and 17. The D values of the three groups (Izu Peninsula, eastern, and western) varied considerably, showing 0.178–0.227 between the Izu Peninsula and the eastern groups, 0.163–0.361 between Izu Penisula and western groups, and 0.114–0.381 between eastern and western groups. Two samples from the Izu Peninsula were genetically very similar to each other (D=0.001) and to E. okadae (D=0.014–0.015), whereas relatively large D values were observed within the eastern and western groups (0.001–0.172 and 0.001–0.198, respectively). Furthermore, among three samples from the Osumi Island Group in the western group, the D values were very low (0.005–0.014), when compared with values between these samples and the samples in this group from other regions (0.037–0.198). Within the eastern group, an extremely small D value was obtained between the samples from Sapporo and Aomori (0.001).
Matrix of Nei's (1978) unbiased genetic distance (below diagonal) and the modified Rogers' distance (Wright, 1978: above diagonal) coefficients for all pairwise comparisons between samples of E. latiscutatus (Sample 1–20), E. okadae (Sample 21), and E. marginatus (Sample 22 and 23).
Fig. 2A shows the UPGMA phenogram based on Nei's genetic distance. In this phenogram, the three samples from the northern Tokara Island Group (Sample 18–20) cluster with the E. marginatus samples, and are clearly separated from another cluster that includes all the remaining samples. In the latter cluster, three subclusters corresponding to the above-mentioned groups (i.e., Izu Peninsula, eastern and western groups) are recognized. The samples from Izu Penisula and E. okadae consititute in a subcluster. Within the western group, three samples from the Osumi Island Group (Sample 15–17) and two samples from the Danjo Island Group (Sample 13 and14) form a compact subgroup.
The topology of the neighbor-joining phenogram (Fig. 2B) is similar to that of the UPGMA one in that two distinct clusters, one consisting of three subclusters, are recognized, although these two phenograms differ slightly in the position of the subcluster consisting of the Izu Peninsula group, and in the internal topology of the two subclusters corresponding to the eastern and western groups.
Principal component analysis
The relationships among samples were also visualized by projecting the sample positions onto the plane defined by the first two principal components of the PCA, which explain 38.7% of the total variation (Fig. 3 and Table 4). PC1 accounted for 22.5% and PC2 for 16.2%. The relationships in the PCA are remarkably similar to those obtained by clustering (Fig. 2). On the first component axis, the northern Tokara samples (Sample 18–20) are completely distinguishable from the remaining samples of E. latiscutatus (Sample 1–17) with a high positive score. The second component completely distinguishes the three groups (Izu Peninsula, eastern, and western) from each other. The western group has a positive PC2 score, whereas the Izu Peninsula and eastern groups have negative scores. The PC2 scores in the eastern group are intermediate between those in the western and Izu Peninsula groups. Moreover, within the western group, the samples from the Osumi Island Group (Sample 15–17) were separated from the remainder by fairly large PC2 scores.
Factor scores of the principal component analysis.
Populations of the northern Tokara Island Group
The distributions of Eumeces latiscutatus and E. marginatus are separated by the Tokara Tectonic Strait, which is located between the northern and southern Tokara Island Group. Morphologically, this is supported by a few external characters, such as the presence or absence of the postnasals (Nagai, 1928; Toyama, 1989; Hikida et al., 1992). However, Kato et al. (1994) demonstrated that the sample of E. latiscutatus from Kuchinoshima within the northern Tokara Island Group is genetically poorly diverged from the samples of E. marginatus, despite their presumed long geographical isolation by the strait. Their result was interpreted as indicative of the northeastward overseas dispersal of E. marginatus across the strait after its formation (Kato et al., 1994). In our study, the UPGMA and NJ phenograms (Fig. 2) show two distinct genetic groups, samples of E. latiscutatus from the northern Tokara Island Group and E. marginatus, and samples of E. latiscutatus and E. okadae. These results strongly suggest that all the northern Tokara populations, formerly assigned to E. latiscutatus, actually belong to E. marginatus, and support the hypothesized northeastward overseas dispersal of E. marginatus (Kato et al., 1994; Hikida and Motokawa, 1999).
Very high levels of genetic differentiation were found among samples within the northern Tokara Island Group (D=0.161–0.317). This result coincides with the high morphological differentiation among insular populations of the northern Tokara Island Group (Hikida et al., 1992), indicating a large amount of genetic divergence among them. Moreover, the Kuchinoshima sample (Sample 18) is most closely related to E. marginatus from Okinawajima (D=0.081), whereas the Suwanosejima sample (Sample 20) is closest to E. marginatus from Amamioshima (D=0.062). Compared with these close genetic affinities, the genetic divergence of the Nakanoshima sample (Sample 19) from both samples of E. marginatus is relatively large (D=0.243 and 0.277). Kato et al. (1994) surmised the founder's effect on a small insular population for explanation of the close genetic affinty between the samples of Kuchinoshima and Okinawajima. However, this is not supported by our data, which indicate that the levels of genetic variability in the samples from Kuchinoshima and Suwanosejima (H=0.060 and 0.084) are higher than those in the samples of E. marginatus (H=0.013 and 0.051). To explain our results, it is necessary to postulate that the northern Tokara populations were colonized from more than one source population by rafting or artificial transportation. Further extensive electro-phoretic surveys of E. marginatus from throughout its range are needed to verify this hypothesis.
Populations of the Izu Peninsula
The high Fst values found in all but the northern Tokara samples of E. latiscutatus (0.740) indicate that E. latiscutatus is highly diverged genetically. Within these samples, the results of our UPGMA and NJ analyses (Fig. 2) and the PC2 values (Fig. 3) reveal three distinct groups (Izu Peninsula, eastern, and western), with a fixed allelic difference at the Idh-2 locus (Table 2). Of these, the samples of the Izu Peninsula group are genetically similar to E. okadae. The samples from the Izu Peninsula possessed alleles unique to E. okadae. Nei's D values between them were much lower (0.014–0.015) than any intraspecific values in this genus estimated previously (0.09–0.68: Murphy et al., 1983; Kato et al., 1994). These results clearly indicate that the samples in the Izu Peninsula belong to the E. okadae lineage and not to any lineage of the other groups in E. latiscutatus, and that E. okadae, now known only from the Izu Isalnd Group, occurs on the Izu Peninsula, at least in part.
Therefore, the samples of E. latiscutatus from the Izu Peninsula and from the other areas, currently treated as a single species, should be treated as two different species. The Izu Peninsula population was assigned to E. latiscutatus because Hallowell (1861) described this species from Shimoda on the Izu Peninsula (sample 5 in this study). The name E. japonicus, which was originally described from Nagasaki in Kyushu (close to sample 11 and 12) by Peters (1864), may be applied to the remaining populations now called E. latiscutatus. Further morphological examination of these two species (E. latiscutatus and E. japonicus) is strongly desired to produce a key for their identification.
By contrast, from the small genetic distance obtained between E. latiscutatus from the Izu Peninsula and E. okadae, it does not follow that the former population is con-specific with E. okadae, because these two species differ in the number of midbody scale rows and the line pattern on the head in juveniles (Taylor, 1936; Nakamura and Uéno, 1963). However, E. okadae, which is distributed on most islands throughout the Izu Island Group, is known to be geographically highly variable, even in the key characters mentioned above (Taylor, 1936; Hikida, 1993). Since Stejneger (1907) first described E. okadae as a subspecies of E. latiscutatus (type locality = Miyakejima), our results seem to imply a conspecific relationship between E. latiscutatus from the Izu Peninsula and E. okadae. In addition, our results alone do not seem to justify the taxonomic treatment of Stejneger (1907), because we examined only one sample from the Izu Island Group. In conclusion, we treat them as a single species, and consider E. okadae a junior synonym of E. latiscutatus. In order to revise the subspecific division of E. latiscutatus recognized by Stejneger (1907), the genetic relationships among insular samples throughout the Izu Island Group and the Izu Peninsula samples of this species need to be clarified, with more detailed study of the morphological variation in this species.
Geographic variation in Eumeces latiscutatus
As shown above, the samples of E. latiscutatus from eastern and western Japan formed two distinct groups in our phenograms (Fig. 2). Phenetic clustering methods such as the UPGMA and NJ methods yield misleading results, such as a large geographic gap in samples exhibiting isolation-by-distance within a single genetically continuous unit, if the distribution of samples was uneven (de Queiroz and Good, 1997). Considering this phenomenon, the apparent existence of two distinct groups (eastern and western) in E. latiscutatus might result solely from our sampling, which is not absolutely regular geographically. However, the high Fst value (0.700) calculated for 15 samples of this species from eastern and western Japan indicates that the samples show a fairly large interpopulational differentiation, which is largely attributed to the intraspecific subdivision of the eastern and western groups. Moreover, the PC2 values completely discriminated the eastern and western groups (Fig. 3), and complete allelic displacement between them was observed (Table 2). These results suggest that little gene flow occurs between the eastern and western groups. Therefore, we believe that there is fairly large genetic divergence between E. latiscutatus from eastern and western Japan.
While there was poor genetic differentiation among the Osumi Island Group samples, they were genetically well diverged from the remaining samples within the western group (Figs. 2 and 3). This genetic divergence is also supported by obvious allele frequency differences and the possession of characteristic alleles. Ota et al. (1993) postulated that the Osumi Island Group has been connected to Kyushu by a narrow land bridge during the most recent glaciation, when the sea level dropped by no less than 120 m. However, the mean D value between samples from the Osumi Island Group and Kagoshima in Kyushu (0.069) seems to be too large to reflect this paleogeographical hypothesis. Based on the evidence of submarine topography, Ohshima (1990) estimated that the separation of the Osumi Island Group from Kyushu occurred during the Riss-Würm interglacial period in the late Pleistocene (120,000 years ago). He also estimated that the separation of the islands in the Osumi Island Group followed the most recent glaciation (i.e., the Würm: 15,000 years ago). If these estimates are correct, the divergence time between samples from the Osumi Island Group and Kyushu might be roughly eight times as long as that within the Osumi Island Group. The mean D value between samples from the Osumi Island Group and Kagoshima (0.069) is eight times as large as that among the Osumi Island Group samples (0.009). Assuming a gross correlation between the genetic distance value and the time of population divergence (Nei, 1987), it is likely that the genetic relationship among samples in the western group reflects the paleogeography in this area proposed by Ohshima (1990).
The populations of several species, including Amphiesma vibakari from the Danjo Island Group, are reported to be very divergent morphologically from conspecific populations from the main islands of Japan (Toriba, 1986; Shibata et al., 1988; Nakanishi and Ejima, 1989). Therefore, we expected the samples from the Danjo Island Group to be genetically distinct from those from the main islands within the western group. However, the PCA did not detect marked differences between samples from the Danjo Island Group and the main islands (Fig. 3). This is attributed to the absence of alleles unique to the samples from the Danjo Island Group. Conversely, our phenograms reveal the relatively distant affinity between samples from the Danjo Island Group and the main islands, supporting their relatively long isolation (Fig. 2). The mean D value observed between samples from the Danjo Island Group and Usuki, Kyushu (0.042) is nearly as large as that between samples from the Osumi Island Group. This suggests that, like the Osumi Island Group, the Danjo Island Group was isolated from Kyushu during the Riss-Würm interglacial period.
Hokkaido and Honshu are assumed to have been separated by the formation of the Tsugaru Strait during the Riss-Würm interglacial period (Ohshima, 1990). However, the sample from Sapporo in Hokkaido proved to be almost identical genetically to that from Aomori in northernmost Honshu (D=0.001), which is incompatible with the paleogeo-graphical scenario (Ohshima, 1990). Therefore, it is plausible that there has been of E. latiscutatus across the Tsugaru Strait overseas dispersal by rafting. This explanation requires substantial corroboration based on more detailed surveys.
We wish to express our sincere thanks to H. Ota and M. Motokawa for pertinent discussion. We also thank M. Toda for the critical reading of the manuscript, and M. Matsui, H. Ota, A. Mori, K. Araya, Y. Misawa, K. Nishikawa, K. Matsuo, S. Narumi and M. Honda for their help in collecting materials examined. This study was partially supported by a grant from the Fujiwara Natural History Foundation (to J. Motokawa) and a Grant-in-Aid from the Japan Ministry of Education, Science, Sports and Culture (C-10836010 to T. Hikida).
Sampling locality for each population and catalogue number of specimens examined.
Sample 1: Sapporo, Hokkaido Prefecture, 43°10′N, 141°15′E, KUZ 45151-4, 50355; sample 2: Aomori, Aomori Pref., 40°38′N, 140°50′E, KUZ 28003, 30349-53, 30384, 30674-7, 35199-200, 35222-3; sample 3: Kamakura, Kanagawa Pref., 35°20′N, 139°33′E, KUZ 45679, 45681-5, 45691; sample 4: Odawara, Kanagawa Pref., 35°15′N, 139°10′E, KUZ 36566, 39000-1, 39004, 39009, 39012-3; sample 5: Shimoda, Shizuoka Pref., 34°40′N, 138°58′E, KUZ 35178-83, 36303-7, 36309-12, 36314, 36318-9, 36321-4; sample 6: Yui, Shizuoka Pref., 35°05′N, 138°33′E, KUZ 45650-4; sample 7: Azumi, Nagano Pref., 36°10′N, 137°46′E, KUZ 39023-8; sample 8: Kasahara, Gifu Pref., 35°18′N, 137°08′E, KUZ 36274-5, 36277, 36279-81, 39388-91, 45929-32; sample 9: Kyoto, Kyoto Pref., 35°03′N, 135°50′E, KUZ 28000, 28185, 28573, 28577, 28583, 28586-7, 28887, 28892, 29280, 30185-6, 30316, 30354-6, 30370, 30383, 30667-8, 35923, 35925-6, 36244, 36250, 36257, 36259, 39051-2, 39316, 39376; sample 10: Miyama, Kyoto Pref., 35°20′N, 135°37′E, KUZ 35121, 35201-2, 35257-60, 35922; sample 11: Usuki, Oita Pref., 33°04′N, 131°44′E, KUZ 50377-85; sample 12: Kagoshima, Kagoshima Pref., 31°36′N, 130°34′E, KUZ 8360-1, 8363, 29301, 30346-8, 30360-2, 36295, 36297-301, 36414-5, 36452-3, 36486-91, 46919; sample 13: Oshima Is., Danjo Island Group, Nagasaki Pref., 32°02′N, 128°25′E, KUZ 50002-10; sample 14: Meshima Is., Danjo Island Group, Nagasaki Pref., 32°00′N, 128°21′E, KUZ 21799-801, 46996-9, 50000-1; sample 15: Tanegashima Is., Osumi Island Group, Kagoshima Pref., 30°24′N, 130°26′E, KUZ 27989-99; sample 16: Yakushima Is., Osumi Island Group, Kagoshima Pref., 30°32′N, 130°58′E, KUZ 35870-80, 35882-4, 36523-8; sample 17: Kuchierabujima Is., Osumi Island Group, Kagoshima Pref., 30°28′N, 130°12′E, KUZ 27978-82, 50403-4; sample 18: Kuchinoshima Is., Tokara Island Group, Kagoshima Pref., 29°39′N, 129°58′E, KUZ 8369-73, 28563-4, 28570, 28916-7; sample 19: Nakanoshima Is., Tokara Island Group, Kagoshima Pref., 29°30′N, 129°54′E, KUZ 28582, 28889; sample 20: Suwanosejima Is., Tokara Island Group, Kagoshima Pref., 29°35′N, 129°41′E, KUZ 8374-7.
Sample 21: Miyakejima Is., Izu Island Group, Tokyo Pref., 34°05′N, 139°30′E, KUZ 35123-35, 35295-8, 35747-8, 35751.
Sample 22: Amamioshima Is., Amami Island Group, Kagoshima Pref., 28°13′N, 129°27′E, KUZ 8379-81, 28165, 28566, 28914, 29295-6, 30333-4, 36421-2, 36454-7, 36459, 36462-3; sample 23: Okinawajima Is., Okinawa Island Group, Okinawa Pref., 26°40′N, 127°54′E, KUZ 27950-1, 27955-60, 28177-84, 28462, 28464, 28466-7, 28478, 28499, 28571, 28908, 30309, 30329, 30357, 30375, 30389, 30679-80.