Jatropha curcas L. (Euphorbiaceae) has acquired a great importance as a renewable source of energy with a number of environmental benefits. Very few attempts were made to understand the extent of genetic diversity and its distribution. This study was aimed to study the diversity and deduce the phylogeography of Jatropha curcas L. which is said to be the most primitive species of the genus Jatropha. Here we studied the intraspecific genetic diversity of the species distributed in different parts of the globe. The study also focused to understand the molecular diversity at reported probable center of origin (Mexico), and to reveal the dispersal route to other regions based on random amplified polymorphic DNA, amplified fragment length polymorphism and nrDNA-ITS sequences data. The overall genetic diversity of J. curcas found in the present study was narrow. The highest genetic diversity was observed in the germplasm collected from Mexico and supports the earlier hypothesis based on morphological data and natural distribution, it is the center for origin of the species. Least genetic diversity found in the Indian germplasm and clustering results revealed that the species was introduced simultaneously by two distinct germplasm and subsequently distributed in different parts of India. The present molecular data further revealed that J. curcas might have spread from the center of the origin to Cape Verde, than to Spain, Portuguese to other neighboring countries and simultaneously to Africa. The molecular evidence supports the Burkill et al. (A dictionary of the economic products of the Malay Peninsula, Governments of Malaysia and Singapore by the Ministry of Agriculture and Co-operatives. Kuala Lumpur, Malaysia, 1966) view of Portuguese might have introduced the species to India. The clustering pattern suggests that the distribution was interfered by human activity.
Mycobacterium abscessus is a rapidly growing mycobacterium that is often associated with human infections. The taxonomy of this species has undergone several revisions and is still being debated. In this study, we sequenced the genomes of 12 M. abscessus strains and used phylogenomic analysis to perform subspecies classification.
Bacterial heart rot of pineapple reported in Hawaii in 2003 and reoccurring in 2006 was caused by an undetermined species of Dickeya. Classification of the bacterial strains isolated from infected pineapple to one of the recognized Dickeya species and their phylogenetic relationships with Dickeya were determined by a multilocus sequence analysis (MLSA), based on the partial gene sequences of dnaA, dnaJ, dnaX, gyrB and recN. Individual and concatenated gene phylogenies revealed that the strains form a clade with reference Dickeya sp. isolated from pineapple in Malaysia and are closely related to D. zeae; however, previous DNA-DNA reassociation values suggest that these strains do not meet the genomic threshold for consideration in D. zeae, and require further taxonomic analysis. An analysis of the markers used in this MLSA determined that recN was the best overall marker for resolution of species within Dickeya. Differential intraspecies resolution was observed with the other markers, suggesting that marker selection is important for defining relationships within a clade. Phylogenies produced with gene sequences from the sequenced genomes of strains D. dadantii Ech586, D. dadantii Ech703 and D. zeae Ech1591 did not place the sequenced strains with members of other well-characterized members of their respective species. The average nucleotide identity (ANI) and tetranucleotide frequencies determined for the sequenced strains corroborated the results of the MLSA that D. dadantii Ech586 and D. dadantii Ech703 should be reclassified as Dickeya zeae Ech586 and Dickeya paradisiaca Ech703, respectively, whereas D. zeae Ech1591 should be reclassified as Dickeya chrysanthemi Ech1591.
Phylogenetic comparisons of the different mammalian genetic transmission elements (mtDNA, X-, Y-, and autosomal DNA) is a powerful approach for understanding the process of speciation in nature. Through such comparisons the unique inheritance pathways of each genetic element and gender-biased processes can link genomic structure to the evolutionary process, especially among lineages which have recently diversified, in which genetic isolation may be incomplete. Bulldog bats of the genus Noctilio are an exemplar lineage, being a young clade, widely distributed, and exhibiting unique feeding ecologies. In addition, currently recognized species are paraphyletic with respect to the mtDNA gene tree and contain morphologically identifiable clades that exhibit mtDNA divergences as great as among many species. To test taxonomic hypotheses and understand the contribution of hybridization to the extant distribution of genetic diversity in Noctilio, we used phylogenetic, coalescent stochastic modeling, and divergence time estimates using sequence data from cytochrome-b, cytochrome c oxidase-I, zinc finger Y, and zinc finger X, as well as evolutionary reconstructions based on amplified fragment length polymorphisms (AFLPs) data. No evidence of ongoing hybridization between the two currently recognized species was identified. However, signatures of an ancient mtDNA capture were recovered in which an mtDNA lineage of one species was captured early in the noctilionid radiation. Among subspecific mtDNA clades, which were generally coincident with morphology and statistically definable as species, signatures of ongoing hybridization were observed in sex chromosome sequences and AFLP. Divergence dating of genetic elements corroborates the diversification of extant Noctilio beginning about 3 Ma, with ongoing hybridization between mitochondrial lineages separated by 2.5 myr. The timeframe of species' divergence within Noctilio supports the hypothesis that shifts in the dietary strategies of gleaning insects (N. albiventris) or fish (N. leporinus) are among the most rapid instances of dietary evolution observed in mammals. This study illustrates the complex evolutionary dynamics shaping gene pools in nature, how comparisons of genetic elements can serve for understanding species boundaries, and the complex considerations for accurate taxonomic assignment.
The authors studied the myxosporean infection of wild gobiid fishes (Perciformes: Gobioidei) in the Merang Estuary of Terengganu, Malaysia, and described Myxobolus ophiocarae sp. n. in Ophiocara porocephala. Several myxosporean plasmodia were found intralamellarly within the gill filaments. The spores differed from those of other Myxobolus species previously recorded on gobiid fishes. They were round in valvular view and lens-shaped in sutural view, and had two equal-sized, pyriform polar capsules with polar filaments having six to seven turns. The spores measured 10.34 × 8.79 × 4.53 μm. The 18S rDNA sequence of M. ophiocarae sp. n., based on a contiguous sequence of 1,789 base pairs, differed from any other Myxobolus spp. in GenBank. Phylogenetic analysis of the 18S rDNA gene revealed that this species showed the closest similarity to Myxobolus nagaraensis, Myxobolus lentisuturalis, and Myxobolus cultus.
Oil palm is the most productive oil-bearing crop. Although it is planted on only 5% of the total world vegetable oil acreage, palm oil accounts for 33% of vegetable oil and 45% of edible oil worldwide, but increased cultivation competes with dwindling rainforest reserves. We report the 1.8-gigabase (Gb) genome sequence of the African oil palm Elaeis guineensis, the predominant source of worldwide oil production. A total of 1.535 Gb of assembled sequence and transcriptome data from 30 tissue types were used to predict at least 34,802 genes, including oil biosynthesis genes and homologues of WRINKLED1 (WRI1), and other transcriptional regulators, which are highly expressed in the kernel. We also report the draft sequence of the South American oil palm Elaeis oleifera, which has the same number of chromosomes (2n = 32) and produces fertile interspecific hybrids with E. guineensis but seems to have diverged in the New World. Segmental duplications of chromosome arms define the palaeotetraploid origin of palm trees. The oil palm sequence enables the discovery of genes for important traits as well as somaclonal epigenetic alterations that restrict the use of clones in commercial plantings, and should therefore help to achieve sustainability for biofuels and edible oils, reducing the rainforest footprint of this tropical plantation crop.
To determine the geographic origin of the black fly Simulium suzukii on Okinawa Island, Japan, macrogenomic profiles derived from its polytene chromosomes were compared with those of mainland and other insular populations of S. suzukii and of the isomorphic Simulium tani species complex. The Okinawan population is a chromosomally unique cytoform, designated 'D,' which is essentially monomorphic and differs by about 27 fixed rearrangements from the chromosomal standard sequence for the subgenus Simulium and by two fixed differences from its nearest known relative, representing the type of S. suzukii, on the main islands of Japan. Chromosomal band sequences revealed two additional, sympatric cytoforms of S. suzukii, designated 'A' and 'B,' each with species status, in Korea, and a third cytoform, designated 'C,' on Hokkaido, Japan. A new cytoform, 'K,' of S. tani from Malaysia, representing the type of S. tani, is more closely related to cytoforms in Thailand, as are populations from Taiwan previously treated as S. suzukii but more closely aligned with S. tani and newly recognized as cytoform 'L' of the latter nominal species. Rooting of chromosomal band sequences by outgroup comparisons allowed directionality of chromosomal rearrangements to be established, enabling phylogenetic inference of cytoforms. Of 41 macrogenomic rearrangements discovered in the five new cytoforms, four provide evidence for a stepwise origin of the Okinawan population from populations characteristic of the main islands of Japan. The macrogenomic approach applied to black flies on Okinawa Island illustrates its potential utility in defining source areas for other species of flies including those that might pose medical and veterinary risks.
Brugia malayi is one of the parasitic worms which causes lymphatic filariasis in humans. Its geographical distribution includes a large part of Asia. Despite its wide distribution, very little is known about the genetic variation and molecular epidemiology of this species. In this study, the internal transcribed spacer 1 (ITS1) nucleotide sequences of B. malayi from microfilaria-positive human blood samples in Northeast Borneo Island were determined, and compared with published ITS1 sequences of B. malayi isolated from cats and humans in Thailand. Multiple alignment analysis revealed that B. malayi ITS1 sequences from Northeast Borneo were more similar to each other than to those from Thailand. Phylogenetic trees inferred using Neighbour-Joining and Maximum Parsimony methods showed similar topology, with 2 distinct B. malayi clusters. The first cluster consisted of Northeast Borneo B. malayi isolates, whereas the second consisted of the Thailand isolates. The findings of this study suggest that B. malayi in Borneo Island has diverged significantly from those of mainland Asia, and this has implications for the diagnosis of B. malayi infection across the region using ITS1-based molecular techniques.
This study addressed the taxonomic position and group classification of a phytoplasma responsible for virescence and phyllody symptoms in naturally diseased Madagascar periwinkle plants in western Malaysia. Unique regions in the 16S rRNA gene from the Malaysian periwinkle virescence (MaPV) phytoplasma distinguished the phytoplasma from all previously described 'Candidatus Phytoplasma' species. Pairwise sequence similarity scores, calculated through alignment of full-length 16S rRNA gene sequences, revealed that the MaPV phytoplasma 16S rRNA gene shared 96.5 % or less sequence similarity with that of previously described 'Ca. Phytoplasma' species, justifying the recognition of the MaPV phytoplasma as a reference strain of a novel taxon, 'Candidatus Phytoplasma malaysianum'. The 16S rRNA gene F2nR2 fragment from the MaPV phytoplasma exhibited a distinct restriction fragment length polymorphism (RFLP) profile and the pattern similarity coefficient values were lower than 0.85 with representative phytoplasmas classified in any of the 31 previously delineated 16Sr groups; therefore, the MaPV phytoplasma was designated a member of a new 16Sr group, 16SrXXXII. Phytoplasmas affiliated with this novel taxon and the new group included diverse strains infecting periwinkle, coconut palm and oil palm in Malaysia. Three phytoplasmas were characterized as representatives of three distinct subgroups, 16SrXXXII-A, 16SrXXXII-B and 16SrXXXII-C, respectively.
The taxonomy of the causal pathogen of basal stem rot of oil palms, Ganoderma is somewhat problematic at present. In order to determine the genetic distance relationship between G. boninense isolates and non-boninense isolates, a random amplified microsatellites DNA (RAMS) technique was carried out. The result was then compared with interfertility data of G. boninense that had been determined in previous mating studies to confirm the species of G. boninense. Dendrogram from cluster analysis based on UPGMA of RAMS data showed that two major clusters, I and II which separated at a genetic distance of 0.7935 were generated. Cluster I consisted of all the biological species G. boninense isolates namely CNLB, GSDK 3, PER 71, WD 814, GBL 3, GBL 6, OC, GH 02, 170 SL and 348781 while all non-boninense isolates namely G. ASAM, WRR, TFRI 129, G. RES, GJ, and CNLM were grouped together in cluster II. Although the RAMS markers showed polymorphisms in all the isolates tested, the results obtained were in agreement with the interfertility data. Therefore, the RAMS data could support the interfertility data for the identification of Ganoderma isolates.
Long columns of migrating larval sciarid armyworms were discovered in central and northern Japan, specifically Kanagawa, Gunma, Miyagi and Akita prefectures, as well as Hokkaido. This is the first examination of armyworms in East Asia. In Europe, armyworms have been identified as Sciara militaris, belonging to the family Sciaridae (sciarid flies or black fungus gnats), by rearing them to adulthood. In Japan, we were unable to obtain live samples for rearing; therefore, DNA barcodes were obtained from the samples of armyworms collected in the Gunma and Miyagi prefectures. The DNA barcodes were compared with those obtained from the following samples: pupae of S. militaris from UK, adults of Sciara kitakamiensis, Sciara humeralis, Sciara hemerobioides, Sciara thoracica, Sciara helvola and Sciara melanostyla from Japan, and adults of one undescribed Sciara species from Malaysia. Neighbour-joining, maximum parsimony, and maximum likelihood analyses revealed that the armyworms discovered in Japan are S. kitakamiensis. Although adults of this species have been recorded in several locations in Japan, this is the first report of migrating larval armyworms. DNA barcodes were effectively used to link different life stages of this species. The average intraspecific and interspecific pairwise genetic distances of the genus Sciara were 0.3% and 12.6%, respectively. The present study illustrates that DNA barcodes are an effective means of identifying sciarid flies in Japan.
Malaysia remains as a crossroad of different cultures and peoples, and it has long been recognized that studying its population history can provide crucial insight into the prehistory of Southeast Asia as a whole. The earliest inhabitants were the Orang Asli in Peninsular Malaysia and the indigenous groups in Sabah and Sarawak. Although they were the earliest migrants in this region, these tribes are divided geographically by the South China Sea. We analyzed DNA sequences of 18 Orang Asli using mitochondrial DNA extracted from blood samples, each representing one sub-tribe, and from five Sarawakian Iban. Mitochondrial DNA was extracted from hair samples in order to examine relationships with the main ethnic groups in Malaysia. The D-loop region and cytochrome b genes were used as the candidate loci. Phylogenetic relationships were investigated using maximum parsimony and neighbor joining algorithms, and each tree was subjected to bootstrap analysis with 1000 replicates. Analyses of the HVS I region showed that the Iban are not a distinct group from the Orang Asli; they form a sub-clade within the Orang Asli. Based on the cytochrome b gene, the Iban clustered with the Orang Asli in the same clade. We found evidence for considerable gene flow between Orang Asli and Iban. We concluded that the Orang Asli, Iban and the main ethnic groups of Malaysia are probably derived from a common ancestor. This is in agreement with a single-route migration theory, but it does not dismiss a two-route migration theory.
Termites of the genus Odontotermes are important decomposers in the Old World tropics and are sometimes important pests of crops, timber and trees. The species within the genus often have overlapping size ranges and are difficult to differentiate based on morphology. As a result, the taxonomy of Odontotermes in Peninsular Malaysia has not been adequately worked out. In this study, we examined the phylogeny of 40 samples of Odontotermes from Peninsular Malaysia using two mitochondrial DNA regions, that is, the 16S ribosomal RNA and cytochrome oxidase subunit I genes, to aid in elucidating the number of species in the peninsula. Phylogenies were reconstructed from the individual gene and combined gene data sets using parsimony and likelihood criteria. The phylogenies supported the presence of up to eleven species in Peninsular Malaysia, which were identified as O. escherichi, O. hainanensis, O. javanicus, O. longignathus, O. malaccensis, O. oblongatus, O. paraoblongatus, O. sarawakensis, and three possibly new species. Additionally, some of our taxa are thought to comprise a complex of two or more species. The number of species found in this study using DNA methods was more than the initial nine species thought to occur in Peninsular Malaysia. The support values for the clades and morphology of the soldiers provided further evidence for the existence of eleven or more species. Higher resolution genetic markers such as microsatellites would be required to confirm the presence of cryptic species in some taxa.
Plastid trnL-trnF and nuclear ribosomal ITS sequences were obtained from selected wild-type individuals of Polygonum minus Huds. in Peninsular Malaysia. The 380 bp trnL-trnF sequences of the Polygonum minus accessions were identical. Therefore, the trnL-trnF failed to distinguish between the Polygonum minus accessions. However, the divergence of ITS sequences (650 bp) among the Polygonum minus accessions was 1%, indicating that these accessions could be distinguished by the ITS sequences. A phylogenetic relationship based on the ITS sequences was inferred using neighbor-joining, maximum parsimony and Bayesian inference. All of the tree topologies indicated that Polygonum minus from Peninsular Malaysia is unique and different from the synonymous Persicaria minor (Huds.) Opiz and Polygonum kawagoeanum Makino.
Recent morphological and molecular studies led to the recognition of two extant species of clouded leopards; Neofelis nebulosa from mainland southeast Asia and Neofelis diardi from the Sunda Islands of Borneo and Sumatra, including the Batu Islands. In addition to these new species-level distinctions, preliminary molecular data suggested a genetic substructure that separates Bornean and Sumatran clouded leopards, indicating the possibility of two subspecies of N. diardi. This suggestion was based on an analysis of only three Sumatran and seven Bornean individuals. Accordingly, in this study we re-evaluated this proposed subspecies differentiation using additional molecular (mainly historical) samples of eight Bornean and 13 Sumatran clouded leopards; a craniometric analysis of 28 specimens; and examination of pelage morphology of 20 museum specimens and of photographs of 12 wild camera-trapped animals. Molecular (mtDNA and microsatellite loci), craniomandibular and dental analyses strongly support the differentiation of Bornean and Sumatran clouded leopards, but pelage characteristics fail to separate them completely, most probably owing to small sample sizes, but it may also reflect habitat similarities between the two islands and their recent divergence. However, some provisional discriminating pelage characters are presented that need further testing. According to our estimates both populations diverged from each other during the Middle to Late Pleistocene (between 400 and 120 kyr). We present a discussion on the evolutionary history of Neofelis diardi sspp. on the Sunda Shelf, a revised taxonomy for the Sunda clouded leopard, N. diardi, and formally describe the Bornean subspecies, Neofelis diardi borneensis, including the designation of a holotype (BM.3.4.9.2 from Baram, Sarawak) in accordance with the rules of the International Code of Zoological Nomenclature.
Infection with Chloromyxum careni Mutschmann, 1999 was found in the Asian horned frog Megophrys nasuta from Malaysia and Indonesia. Kidney was the only organ infected. Coelozoic plasmodia up to 300 μm were localized in Bowman's space, embracing the glomerulus from all sides, or rarely in lumina of renal tubules. Plasmodia are polysporic, containing disporic pansporoblasts. Myxospores observed by light microscopy are colorless, variable in shape and size, measuring 6.0-8.5 × 5.0-6.5 μm, composed of two symmetrical valves joined by a meridian suture, containing four pyriform polar capsules 3.0-4.0 × 2.5-3.0 μm and a single sporoplasm. Each valve possesses 14-24 (median 21) fine longitudinal ridges clearly visible only in scanning electron microscopy. Rarely, atypical spores with a markedly pointed posterior pole and only 6-10 surface ridges are present in plasmodia together with typical spores. Both small subunit (SSU) and large subunit (LSU) rRNA gene sequences possess extremely long GU-rich inserts. In all SSU and LSU rDNA-based phylogenetic analyses, C. careni clustered as a distinct basal branch to the Myxobolus+Myxidium lieberkuehni clade, out of the marine Chloromyxum clade containing Chloromyxum leydigi, the type species of the genus. These morphological and phylogenetic data suggest erection of a new genus for the C. careni lineage, but we conservatively treat it as a Chloromyxum sensu lato until more information is available.
Japanese encephalitis virus (JEV) is the most important cause of epidemic encephalitis worldwide but its origin is unknown. Epidemics of encephalitis suggestive of Japanese encephalitis (JE) were described in Japan from the 1870s onwards. Four genotypes of JEV have been characterised and representatives of each genotype have been fully sequenced. Based on limited information, a single isolate from Malaysia is thought to represent a putative fifth genotype. We have determined the complete nucleotide and amino acid sequence of Muar strain and compared it with other fully sequenced JEV genomes. Muar was the least similar, with nucleotide divergence ranging from 20.2 to 21.2% and amino acid divergence ranging from 8.5 to 9.9%. Phylogenetic analysis of Muar strain revealed that it does represent a distinct fifth genotype of JEV. We elucidated Muar signature amino acids in the envelope (E) protein, including E327 Glu on the exposed lateral surface of the putative receptor binding domain which distinguishes Muar strain from the other four genotypes. Evolutionary analysis of full-length JEV genomes revealed that the mean evolutionary rate is 4.35 × 10(-4) (3.4906 × 10(-4) to 5.303 × 10(-4)) nucleotides substitutions per site per year and suggests JEV originated from its ancestral virus in the mid 1500s in the Indonesia-Malaysia region and evolved there into different genotypes, which then spread across Asia. No strong evidence for positive selection was found between JEV strains of the five genotypes and the E gene has generally been subjected to strong purifying selection.
By investigating genealogical relationships, we estimated the phylogenetic history and biogeography in the megophryid genus Leptobrachium (sensu lato, including Vibrissaphora) from southern China, Indochina, Thailand and the Sundaland. The genealogical relationships among the 30 named and unnamed taxa were estimated using 2009 bp of sequences from the mitochondrial DNA genes 12S rRNA, tRNA(val), and 16S rRNA using maximum parsimony, maximum likelihood, and Bayesian inference methods. The genus Leptobrachium was a well-supported monophyletic group that contained two major clades. One clade had three subclades primarily from disjunct regions including Borneo, Peninsular Malaysia and Java, and Thailand. The Bornean subclade included one species each from the Philippines and Sumatra. The other major clade consisted of two subclades, one from Indochina and the other from southern China (Vibrissaphora). Divergence times estimated an old evolutionary history of each subclade, one that could not be explained by the geohistory of Southeast Asian major landmasses.
Human enterovirus 71 (EV-71) is genotyped for molecular epidemiological investigation mainly using the two structural genes, VP1 and VP4. Based on these, EV-71 is divided into three genotypes, A, B and C, and within the genotypes B and C, there are further subgenotypes, B1-B5 and C1-C5. Classification using these genes is useful but gives incomplete phylogenetic information. In the present study, the phylogenetic relationships amongst all the known EV-71 and human enterovirus A (HEV-A) isolates with complete genome sequences were examined. A different tree topology involving EV-71 isolates of subgenotypes, C4 and B5 was obtained in comparison to that drawn using VP1. The nucleotide sequence divergence of the C4 isolates was 18.11% (17-20%) when compared to other isolates of subgenotype C. However, this positions the C4 isolates within the cut-off divergence value of 17-22% used to designate the virus genotypes. Hence, it is proposed here that C4 should be designated as a new genotype D. In addition, the subgenotype B5 isolates had an average nucleotide divergence of only 6.14% (4-8%) when compared to other subgenotype B4 isolates. This places the B5 isolates within the subgenotype B4. It is proposed here that the B5 isolates to be redesignated as B4. With the newly proposed genotype D and inclusion of subgenotype B5 within B4, the average nucleotide divergence between genotypes was 18.99% (17-22%). Inter- and intra-subgenotype average divergences were 12.02% (10-14%) and 3.92% (1-10%), respectively. A phylogenetic tree built using the full genome sequences is robust as it takes into consideration changes in the sequences of both the structural and non-structural genes. Similar nucleotide similarities, however, were obtained if only VP1 and 3D RNA polymerase genes were used. Furthermore, addition of 3D RNA polymerase sequences will also show recombination events. Hence, in the absence of full genome sequences, it is proposed here that a combination of VP1 and 3D RNA polymerase gene sequences be used for initial genotyping of EV-71 isolates.
Several anuran groups of Laurasian origin are each co-distributed in four isolated regions of the Northern Hemisphere: central/southern Europe and adjacent areas, Korean Peninsula and adjacent areas, Indo-Malaya, and southern North America. Similar distribution patterns have been observed in diverse animal and plant groups. Savage [Savage, J.M., 1973. The geographic distribution of frogs: patterns and predictions. In: Vial, J.L. (Ed.), Evolutionary Biology of the Anurans. University of Missouri Press, Columbia, pp. 351-445] hypothesized that the Miocene global cooling and increasing aridities in interiors of Eurasia and North America caused a southward displacement and range contraction of Laurasian frogs (and other groups). We use the frog genus Bombina to test Savage's biogeographical hypothesis. A phylogeny of Bombina is reconstructed based on three mitochondrial and two nuclear gene fragments. The genus is divided into three major clades: an Indo-Malaya clade includes B. fortinuptialis, B. lichuanensis, B. maxima, and B. microdeladigitora; a European clade includes B. bombina, B. pachypus, and B. variegata; and a Korean clade contains B. orientalis. The European and Korean clades form sister-group relationship. Molecular dating of the phylogenetic tree using the penalized likelihood and Bayesian analyses suggests that the divergence between the Indo-Malaya clade and other Bombina species occurred 5.9-28.6 million years ago. The split time between the European clade and the Korean clade is estimated at 5.1-20.9 million years ago. The divergence times of these clades are not significantly later than the timing of Miocene cooling and drying, and therefore can not reject Savage's hypothesis. Some other aspects of biogeography of Bombina also are discussed. The Korean Peninsula and the Shandong Peninsula might have supplied distinct southern refugia for B. orientalis during the Pleistocene glacial maxima. In the Indo-Malaya clade, the uplift of the Tibetan Plateau might have promoted the split between B. maxima and the other species.