Molecular phylogenetic status of the Korean goral and Japanese serow based on partial sequences of the mitochondrial cytochrome b gene

To investigate the molecular phylogenetic status of the Korean goral, Nemorhaedus caudatus raddeanus, and Japanese serow, Capricornis crispus, we determined partial sequences of the mitochondrial cytochrome b gene of twelve Korean gorals and sixteen Japanese serows, and compared them with those of the major lineages of Rupicaprini species including two other Nemorhaedus species and two other Capricornis species. The Korean gorals examined possessed two haplotypes with only one nucleotide difference between them, while the Japanese serows showed slightly higher sequence diversity with five haplotypes. Genetic distances and molecular phylogenetic trees indicated that there is considerable genetic divergence between the Korean goral and N. caudatus (the Chinese goral) [Groves and Shields (1996)], but virtually none between Korean and Russian gorals. The Korean and Russian gorals may therefore be distinct from the Chinese goral. The data highlight the importance of conservation of the goral populations of these regions, and the need to reconsider the taxonomic status of Korean and Russian gorals. Our study also clearly demonstrated sufficient genetic distance between serows and gorals to justify their assignment to separate genera. Of the three species of Capricornis, the Formosan serow, C. swinhoei is more closely related to C. sumatraensis than to the Japanese serow, suggesting that the Formosan serow is a distinct species. Preliminary data on intraspecific genetic variation in the Japanese serow are also presented.     

Reclassification of the serows and gorals: the end of a neverending story?

相似文献   

Complete mitochondrial genome of the Himalayan serow (Capricornis thar) and its phylogenetic status within the genus Capricornis

The genus Capricornis (Caprinae, Bovidae) includes six species, i.e., Capriconis sumatraensis, Capricornis swinhoei, Capricornis crispus, Capricornis rubidus, Capricornis milneedwardsii and Capricornis thar. Based on morphological, anatomical and physiological evidence, the Himalayan serow (C. thar), which is distributed in the mountain regions of the Tibetan Plateau and adjacent areas, is generally classified as a species of the genus Capricornis. However, considering the lack of molecular phylogenetic evidence, a heated controversy still exists over the taxonomic position of the Himalayan serow. To determine the phylogenetic status of C. thar, we sequenced the complete mitochondrial genome (mtDNA) of one Himalayan serow using the polymerase chain reaction (PCR) technique. Then, we constructed a molecular phylogenetic tree of the genus Capricornis and its related species in the Caprinae subfamily based on maximum likelihood and Bayesian interference methods. The results reflect that the structural characteristics, nucleotide content and codon usage bias of the mitochondrial genome of the Himalayan serow were similar to those of other serows. A phylogenetic analysis using two rRNA genes and 12 encoded protein-coding genes indicated that C. thar is clustered within the C. milneedwardsii clade, which is a sister group of C. sumatraensis. The present study provides useful information about the evolution of the Himalayan serow, which will be essential for conservation genetic studies.     

Population stock structure of leatherback turtles (Dermochelys coriacea) in the Atlantic revealed using mtDNA and microsatellite markers

This study presents a comprehensive genetic analysis of stock structure for leatherback turtles (Dermochelys coriacea), combining 17 microsatellite loci and 763 bp of the mtDNA control region. Recently discovered eastern Atlantic nesting populations of this critically endangered species were absent in a previous survey that found little ocean-wide mtDNA variation. We added rookeries in West Africa and Brazil and generated longer sequences for previously analyzed samples. A total of 1,417 individuals were sampled from nine nesting sites in the Atlantic and SW Indian Ocean. We detected additional mtDNA variation with the longer sequences, identifying ten polymorphic sites that resolved a total of ten haplotypes, including three new variants of haplotypes previously described by shorter sequences. Population differentiation was substantial between all but two adjacent rookery pairs, and F _ST values ranged from 0.034 to 0.676 and 0.004 to 0.205 for mtDNA and microsatellite data respectively, suggesting that male-mediated gene flow is not as widespread as previously assumed. We detected weak (F _ST = 0.008 and 0.006) but significant differentiation with microsatellites between the two population pairs that were indistinguishable with mtDNA data. POWSIM analysis showed that our mtDNA marker had very low statistical power to detect weak structure (F _ST < 0.005), while our microsatellite marker array had high power. We conclude that the weak differentiation detected with microsatellites reflects a fine scale level of demographic independence that warrants recognition, and that all nine of the nesting colonies should be considered as demographically independent populations for conservation. Our findings illustrate the importance of evaluating the power of specific genetic markers to detect structure in order to correctly identify the appropriate population units to conserve.     

Analysis of Human Mitochondrial DNA Polymorphisms in the Japanese Population

The highly polymorphic nature and high amplification efficiency of mitochondrial DNA (mtDNA) is valuable for the analysis of biological evidence in forensic casework, such as the identification of individuals and assignment of race/ethnicity. To be useful, a mtDNA polymorphism database for the Japanese population requires an understanding of the range of haplotype variation and phylogenies of mtDNA sequences. To extend current knowledge on the haplotypes in the Japanese population, this study defines new lineages and provides more detail about some of those previously described. We compared the hypervariable regions (HVRs) of 270 healthy, unrelated Japanese individuals and demonstrated 192 haplotypes. Combining HVR1 and HVR2, the genetic diversity was 0.9935, thus providing a high level of identification capability. Haplogroup status was defined for 160 individuals using HVR1, HVR2, and particular coding region polymorphisms; these individuals belonged to 94 haplotypes, four of which were new lineages. The complete mtDNA sequence was also determined from seven individuals.     

Population structure of North Atlantic and North Pacific sei whales (<Emphasis Type="Italic">Balaenoptera borealis</Emphasis>) inferred from mitochondrial control region DNA sequences and microsatellite genotypes

Currently, three stocks of sei whales (Balaenoptera borealis) are defined in the North Atlantic; the Nova Scotian, Iceland-Denmark Strait and Eastern North Atlantic stocks, which are mainly based upon historical catch and sighting data. We analyzed mitochondrial control region DNA (mtDNA) sequences and genotypes from 7 to 11 microsatellite loci in 87 samples from three sites in the North Atlantic; Iceland, the Gulf of Maine and the Azores, and compared against the North Pacific using 489 previously published samples. No statistically significant deviations from homogeneity were detected among the North Atlantic samples at mtDNA or microsatellite loci. The genealogy estimated from the mtDNA sequences revealed a clear division of the haplotypes into a North Atlantic and a North Pacific clade, with the exception of one haplotype detected in a single sample from the Azores, which was included in the North Pacific clade. Significant genetic divergence between the North Atlantic and North Pacific Oceans was detected (mtDNA Φ_ST?=?0.72, microsatellite Weir and Cockerham’s ? = 0.20; p?<?0.001). The coalescent-based estimate of the population divergence time between the North Atlantic and North Pacific populations from the sequence variation among the mtDNA sequences was at 163,000 years ago. However, the inference was limited by an absence of samples from the Southern Hemisphere and uncertainty regarding mutation rates and generation times. The estimates of inter-oceanic migration rates were low (Nm at 0.007 into the North Pacific and at 0.248 in the opposite direction). Although estimates of genetic divergence among the current North Atlantic stocks were low and consistent with the extensive range of movement observed in satellite tagged sei whales, the high uncertainty of the genetic divergence estimates precludes rejection of multiple stocks in the North Atlantic.     

Organization and variation of the mitochondrial DNA control region in five Caprinae species

The complete mitochondrial DNA (mtDNA) control region was analyzed from five species of the subfamily Caprinae; Naemorhedus caudatus, N. goral, Capra hircus, Capricornis swinhoei, and Capricornis crispus. Among these species, the control region ranged from 1,096 to 1,212 bp in length. Our results were compatible with the scheme of three domains (ETAS, Central, and CSB) within the control region. A + T < G + C was observed in all the domains. In the Korean gorals, of the 31 variable sites in the whole control region resulting in 15 haplotypes, 27 variable sites were in the ETAS domain. We found two to three tandem repeat in all five species examined in this study, three in N. caudatus and N. goral, two in Capra hircus and C. crispus, and one in C. swinhoei, respectively. All of these repeat units include two short sections of mirror symmetry (TACAT and ATGTA). Short mirror symmetries were well-resolved among five different species, although left domain has high substitution rates. By Kimura’s two parameter method, the genetic distances between the genera Naemorhedus and Capricornis were calculated and divergence time between Naemorhedus and Capricornis may be nearly 2Myr.     

Establishment of vascular endothelial cell lines from the aortas of wild Japanese serows (Capricornis crispus)

Endothelial cell lines were established from the aortas of wild Japanese serows (Capricornis crispus) by transfection of a simian virus 40 (SV40) large T antigen gene. The cloned cell lines, designed SeET (Japanese serow endothelial-SV40T) cells, express SV40T antigen and retain cobblestone-like morphology. Although von Willbrand Factor (vWF) is expressed in the cells, the expression rate and the quantity are lower than in serow primary endothelial cells. The SeET cells exhibit positive uptake of acetylated low density lipoprotein and dose-dependent cell proliferation upon exposure to vascular endothelial growth factor. These results suggest that these SeET cells have preserved endothelial phenotypes and able to function with decreased expression of vWF. The SeET cell line will be a valuable tool for in vitro studies on the physiological properties of endothelial cells and for the propagation of viruses and parasites of Japanese serows.     

Molecular phylogeographic analysis of East Asian cryptoperlan stoneflies (Insecta: Plecoptera, Peltoperlidae)

Cryptoperlan stoneflies inhabit the headwaters or upper stream areas of rivers. A molecular phylogeographic study of cryptoperlans in the Japanese archipelago and on Taiwan Island has been conducted. Altogether the mtDNA 16S rRNA region of 71 individuals from 61 populations, the mtDNA COI region of 76 individuals from 41 populations, and the nDNA Histone 3 region of 56 individuals from 52 populations were sequenced and analyzed. The respective ML, NJ, MP and Bayesian dendrograms were proposed from the sequencing data for the 16S rRNA region (362-bp), the COI region (540-bp), and the Histone 3 region (322-bp), estimated using Yoraperla uenoi as an outgroup. Based upon those data and the resulting dendrograms, it has become clear that the cryptoperlan stoneflies of the Japanese archipelago and those of Taiwan Island comprise two major clades. The first of these two major clades consists of a number of OTUs [operational taxonomic units: Cryptoperla japonica (Honshu, Shikoku and Kyushu Islands) + C. ishigakiensis (Ishigaki-jima Island) + Cryptoperla spp. (Okinawa-jima and Taiwan Islands)]. The other clade consists of the species Cryptoperla kawasawai inhabiting only Shikoku Island. Of particular note, C. kawasawai was observed to be significantly genetically differentiated from all other cryptoperlans examined. Yet, despite the fact that the specimens of C. japonica were taken from a very broad range of populations, their genetic diversity was relatively low, similar to that of C. kawasawai, which inhabits only a limited region within Shikoku Island. Furthermore, even the species C. kawasawai was revealed to be composed of two significantly genetically differentiated subclades. It is considered that this genetic structure among cryptoperlans largely reflects the geological history from the middle to upper Miocene Epoch (i.e., Tortonian stage) of the Japanese archipelago and Taiwan Island.     

Complete mitochondrial DNA sequences of two haplotypes of the smallmouth bass, <Emphasis Type="Italic">Micropterus dolomieu</Emphasis>, collected from nonindigenous populations in Japan

We determined the complete mitochondrial genome sequences of two haplotypes of the smallmouth bass, Micropterus dolomieu; individuals used in the analysis were collected from nonindigenous populations in Japan. Both genomes comprised 16,488 bp, with genome contents and gene orders being identical to those of other teleost fishes. A previous study revealed that the Japanese smallmouth bass had only two haplotypes, and the present study revealed that the complete mitochondrial DNA (mtDNA) sequences of the haplotypes differed in only one nucleotide difference. The low genetic diversity in the mtDNA of the smallmouth bass individuals in our study and the results of the mtDNA sequence comparison between the Japanese and the North American individuals suggested that the fish had been transplanted from a fish farm with a low-diversity stock.     

Limited genetic diversity in the critically endangered Mexican howler monkey (Alouatta palliata mexicana) in the Selva Zoque,Mexico

The Mexican howler monkey (Alouatta palliata mexicana) is a critically endangered primate, which is paleoendemic to Mexico. However, despite the potential significance of genetic data for its management and conservation, there have been no population genetic studies of this subspecies. To examine genetic diversity in the key remaining forest refuge for A. p. mexicana, the Selva Zoque, we amplified full-length mitochondrial control region sequences (1,100 bp) from 45 individuals and found 7 very similar haplotypes. Haplotype diversity (h = 0.486) and nucleotide diversity (π = 0.0007) were extremely low compared to other Neotropical primates. Neutrality tests, used to evaluate demographic effects (Tajima’s D = ?1.48, p = 0.05; Fu’s F _s = ?3.33, p = 0.02), and mismatch distribution (sum of squares deviation = 0.006, p = 0.38; raggedness index = 0.12, p = 0.33) were consistent with a recent and mild population expansion and genetic diversity appears to be historically low in this taxon. Future studies should use a combination of mitochondrial and nuclear markers to fully evaluate genetic diversity and to better understand demographic history in A. p. mexicana. These studies should be undertaken throughout its geographic range in order to evaluate population structure and identify management units for conservation. Due to the limited distribution and population size of A. p. mexicana, future conservation strategies may need to consider genetic management. However, a more detailed knowledge of the population genetics of the subspecies is urgently recommended to maximise the conservation impact of these strategies.     

Genetic variation of complete mitochondrial genome sequences of the Sumatran rhinoceros (<Emphasis Type="Italic">Dicerorhinus sumatrensis</Emphasis>)

The Sumatran rhinoceros (Dicerorhinus sumatrensis) is the smallest and one of the most endangered rhinoceros species, with less than 100 individuals estimated to live in the wild. It was originally divided into three subspecies but only two have survived, D. sumatrensis sumatrensis (Sumatran subspecies), and D. s. harrissoni (Bornean). Questions regarding whether populations of the Sumatran rhinoceros should be treated as different management units to preserve genetic diversity have been raised, particularly in light of its severe decline in the wild and low breeding success in captivity. This work aims to characterize genetic differentiation between Sumatran rhinoceros subspecies using complete mitochondrial genomes, in order to unravel their maternal evolutionary history and evaluate their status as separate management units. We identified three major phylogenetic groups with moderate genetic differentiation: two distinct haplogroups comprising individuals from both the Malay Peninsula and Sumatra, and a third group from Borneo. Estimates of divergence time indicate that the most recent common ancestor of the Sumatran rhinoceros occurred approximately 360,000 years ago. The three mitochondrial haplogroups showed a common divergence time about 80,000 years ago corresponding with a major biogeographic event in the Sundaland region. Patterns of mitochondrial genetic differentiation may suggest considering Sumatran rhinoceros subspecies as different conservation units. However, the management of subspecies as part of a metapopulation may appear as the last resource to save this species from extinction, imposing a conservation dilemma.     

Demographic History and Population Structure of Blackfin Flounder (Glyptocephalus stelleri) in Japan Revealed by Mitochondrial Control Region Sequences

The demographic history and population genetic structure of the blackfin flounder (Glyptocephalus stelleri) along coastal regions of Japan were investigated. Genetic variation in DNA sequences was examined from the first hypervariable region of the mitochondrial DNA control region. A high level of haplotypic diversity (h = 0.99 ± 0.004) was detected, indicating a high level of intrapopulation genetic diversity. The starburst structure of the minimum spanning tree suggested a very recent origin for most haplotypes. The demographic history of G. stelleri was examined using neutrality tests and mismatch distribution analysis, which also indicated a Pleistocene population expansion at about 124,100–413,400 years ago. Hierarchical molecular variance analysis and conventional population Fst comparisons revealed no significant genetic differentiation throughout the range examined.     

Genetic Population Structure of Chum Salmon in the Pacific Rim Inferred from Mitochondrial DNA Sequence Variation

We examined the genetic population structure of chum salmon, Oncorhynchus keta, in the Pacific Rim using mitochondrial (mt) DNA analysis. Nucleotide sequence analysis of about 500 bp in the variable portion of the 5′ end of the mtDNA control region revealed 20 variable nucleotide sites, which defined 30 haplotypes of three genealogical clades (A, B, and C), in more than 2,100 individuals of 48 populations from Japan (16), Korea (1), Russia (10), and North America (21 from Alaska, British Columbia, and Washington). The observed haplotypes were mostly associated with geographic regions, in that clade A and C haplotypes characterized Asian populations and clade B haplotypes distinguished North American populations. The haplotype diversity was highest in the Japanese populations, suggesting a greater genetic variation in the populations of Japan than those of Russia and North America. The analysis of molecular variance and contingency χ² tests demonstrated strong structuring among the three geographic groups of populations and weak to moderate structuring within Japanese and North American populations. These results suggest that the observed geographic pattern might be influenced primarily by historic expansions or colonizations and secondarily by low or restricted gene flow between local groups within regions. In addition to the analysis of population structure, mtDNA data may be useful for constructing a baseline for stock identification of mixed populations of high seas chum salmon.     

Population structure and conservation implications for the loggerhead sea turtle of the Cape Verde Islands

The Cape Verde Islands harbour the second largest nesting aggregation of the globally endangered loggerhead sea turtle in the Atlantic. To characterize the unknown genetic structure, connectivity, and demographic history of this population, we sequenced a segment of the mitochondrial (mt) DNA control region (380 bp, n = 186) and genotyped 12 microsatellite loci (n = 128) in females nesting at three islands of Cape Verde. No genetic differentiation in either haplotype or allele frequencies was found among the islands (mtDNA F _ST = 0.001, P > 0.02; nDNA F _ST = 0.001, P > 0.126). However, population pairwise comparisons of the mtDNA data revealed significant differences between Cape Verde and all previously sequenced Atlantic and Mediterranean rookeries (F _ST = 0.745; P < 0.000). Results of a mixed stock analysis of mtDNA data from 10 published oceanic feeding grounds showed that feeding grounds of the Madeira, Azores, and the Canary Islands, in the Atlantic Ocean, and Gimnesies, Pitiüses, and Andalusia, in the Mediterranean sea, are feeding grounds used by turtles born in Cape Verde, but that about 43% (±19%) of Cape Verde juveniles disperse to unknown areas. In a subset of samples (n = 145) we evaluated the utility of a longer segment (~760 bp) amplified by recently designed mtDNA control region primers for assessing the genetic structure of Atlantic loggerhead turtles. The analysis of the longer fragment revealed more variants overall than in the shorter segments. The genetic data presented here are likely to improve assignment and population genetic analyses, with significant conservation and research applications.     

Extensive Female-Mediated Gene Flow and Low Phylogeography Among Seventeen Goat Breeds in Southwest China

Indigenous Chinese goat mtDNA is highly diverse but lacks geographic specificity; however, whether gene flow or gene exchange contributed to this remains unknown. We reanalyzed a consensus fragment of 481 bp in the D-loop region from 339 individuals. The network and neighbor-joining tree revealed three divergent maternal haplogroups (A, B₁, and B₂) in 17 local breeds. Although high polymorphism resulting in 198 different haplotypes was observed (h = 0.984 ± 0.002; π = 0.0336 ± 0.0008), neither the distribution of haplotypes nor PCA analysis revealed any obvious geographic structure in the local breeds. Extensive gene flow was widely detected among breeds from southwest China. High levels of gene exchange were detected between Qianbei Brown goats and the other breeds, indicating either more contribution or introgression to their gene pools. This study will be helpful in understanding the phylogeography and gene flow among the goat breeds of southwest China.     

Population genetic structure of the endangered Eastern Bristlebird, Dasyornis brachypterus; implications for conservation

For species that are habitat specialists or sedentary, population fragmentation may lead to genetic divergence between populations and reduced genetic diversity within populations, with frequent inbreeding. Hundreds of kilometres separate three geographical regions in which small populations of the endangered Eastern Bristlebird, Dasyornis brachypterus, a small, ground-dwelling passerine that occurs in fire-prone bushland in eastern Australia, are currently found. Here, we use mitochondrial and microsatellite DNA markers to: (i) assess the sub-specific taxonomy designated to northern range-edge, and central and southern range-edge D. brachypterus, respectively, and (ii) assess levels of standing genetic variation and the degree of genetic subdivision of remnant populations. The phylogenetic relationship among mtDNA haplotypes and their spatial distribution did not support the recognised subspecies boundaries. Populations in different regions were highly genetically differentiated, but in addition, the two largest, neighboring populations (located within the central region and separated by ~50 km) were moderately differentiated, and thus are likely closed to migration (microsatellites, F _ST = 0.06; mtDNA, F _ST = 0.12, ?? _ST = 0.08). Birds within these two populations were genotypically diverse and apparently randomly mating. A long-term plan for the conservation of D. brachypterus??s genetic diversity should consider individual populations as separate management units. Moreover, managers should avoid actively mixing birds from different populations or regions, to conserve the genetic integrity of local populations and avoid outbreeding depression, should further translocations be used as a recovery tool for this species.     

Genetic variation of the East Balkan Swine (Sus scrofa) in Bulgaria,revealed by mitochondrial DNA and Y chromosomal DNA

East Balkan Swine (EBS) Sus scrofa is the only aboriginal domesticated pig breed in Bulgaria and is distributed on the western coast of the Black Sea in Bulgaria. To reveal the breed's genetic characteristics, we analysed mitochondrial DNA (mtDNA) and Y chromosomal DNA sequences of EBS in Bulgaria. Nucleotide diversity (π_n) of the mtDNA control region, including two newly found haplotypes, in 54 EBS was higher (0.014 ± 0.007) compared with that of European (0.005 ± 0.003) and Asian (0.006 ± 0.003) domestic pigs and wild boar. The median‐joining network based on the mtDNA control region showed that the EBS and wild boar in Bulgaria comprised mainly two major mtDNA clades, European clade E1 (61.3%) and Asian clade A (38.7%). The coexistence of two mtDNA clades in EBS in Bulgaria may be the relict of historical pig translocation. Among the Bulgarian EBS colonies, the geographical differences in distribution of two mtDNA clades (E1 and A) could be attributed to the source pig populations and/or historical crossbreeding with imported pigs. In addition, analysis of the Y chromosomal DNA sequences for the EBS revealed that all of the EBS had haplotype HY1, which is dominant in European domestic pigs.     

Complex patterns of population genetic structure of moose, Alces alces, after recent spatial expansion in Poland revealed by sex-linked markers

In recent years, human activity directly and indirectly influenced the demography of moose in Poland. The species was close to extinction, and only a few isolated populations survived after the Second World War; then, unprecedented demographic and spatial expansions had occurred, possibly generating a very complex pattern of population genetic structure at the present-day margins of the species range in Poland. Over 370 moose from seven populations were collected from Poland, and partial sequences of the mitochondrial control region (mtDNA-cr; 607 bp) were obtained. In addition, the entire mtDNA cytochrome b gene (1,140 bp) and Y-chromosome markers (1,982 bp in total) were studied in a chosen set of individuals. Twelve mtDNA haplotypes that all belonged to the European moose phylogroup were recorded. They could be divided into two distinct clades: Central Europe and the Ural Mountains. The first clade consists of three distinct groups/branches: Biebrza, Polesie, and Fennoscandia. The Biebrza group has experienced spatial and demographic expansion in the recent past. Average genetic differentiation among moose populations in Poland at mtDNA-cr was great and significant (Φ _ST?=?0.407, p?<?0.001). Using mtDNA-cr data, four separate groups of population were recognized using spatial analysis of molecular variance and principal coordinate analysis, including a relict population in Biebrza National Park, a reintroduced Kampinos National Park population, as well as populations that were descendants of moose that colonized Poland from the east (Lithuania, Belarus, and Ukraine) and the north (former East Prussia). Among all the sequenced Y-chromosome markers, polymorphisms were found in the DBY14 marker in three populations only; four haplotypes were recorded in total. No significant differentiation was detected for this Y-linked marker among moose populations in Poland. Our mtDNA study revealed that a variety of different factors—bottleneck, the presence of relict, autochthonous populations, translocations, limited female dispersal, and the colonization from the east and north—are responsible for the observed complex pattern of population genetic structure after demographic and spatial expansion of moose in Poland.     

Phylogeography of Cycas revoluta Thunb. (Cycadaceae) on the Ryukyu Islands: very low genetic diversity and geographical structure

Variations in the chloroplast and mitochondrial DNA of Cycas revoluta Thunb. (Cycadaceae) were examined in 22 populations distributed across the Ryukyu Islands and southern Kyushu. Among the 14,130 bp of sequence examined, only one site mutation and one indel were polymorphic. The identified polymorphisms were located in the spacers between trnS (UGA) and trnfM (CAU) of the chloroplast DNA and between nad1 exon B and exon C of the mitochondrial DNA, respectively. Three haplotypes were identified from the Ryukyu Islands and southern Kyushu. The areas of distribution of the three haplotypes were highly geographically structured. The boundaries of two of the three haplotypes were demarcated by Okinoerabujima Island in the middle Ryukyus. The northern type and southern types lay north and south of the island, respectively. The third haplotype was almost sympatrically distributed with the southern type. The genetic variation within C. revoluta was estimated to be very low (h = 0.641, π = 0.00071) in comparison to its relative in Taiwan, C. taitungensis, which possesses 97 cpDNA haplotypes and 55 mtDNA haplotypes from two relic populations. A reasonable explanation for the low genetic diversity of the cycad on the Ryukyu Islands could be severe bottleneck effects, resulting from the submersion of low islands and the diminished landmass of islands in the interglacial age in the Quaternary period. The geographically restricted nature of the haplotypes could be attributed to vicariance resulting from the land configuration of the Ryukyu Islands, including changes in geography during the interglacial age in the Quaternary.