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.
The sequences of the mitochondrial ND4 gene (1339 bp) and the ND4L gene (290 bp) were determined for all the 14 extant taxa of the Drosophila nasuta subgroup. The average A + T content of ND4 genes is 76.5% and that of ND4L genes is 83.5%. A total of 114 variable sites were scored. The ND4 gene sequence divergence ranged from 0 to 5.4% within the subgroup. The substitution rate of the ND4 gene is about 1.25% per million years. The base substitution of the genes is strongly transition biased. Neighbor-joining and parsimony were used to construct a phylogeny based on the resultant sequence data set. According to these trees, five distinct mtDNA clades can be identified. D. niveifrons represents the most diverged lineage. D. sulfurigaster bilimbata and D. kepulauana form two independent lineages. The other two clades are the kohkoa complex and the albomicans complex. The kohkoa complex consists of D. sulfurigaster sulfurigaster, D. pulaua, D. kohkoa, and Taxon-F. The albomicans complex can be divided into two groups: D. nasuta, D. sulfurigaster neonasuta, D. sulfurigaster albostrigata, and D. albomicans from Chiangmai form one group; and D. pallidifrons, Taxon-I, Taxon-J, and D. albomicans from China form the other group. High genetic differentiation was found among D. albomicans populations. Based on our phylogenetic results, we hypothesize that D. niveifrons diverged first from the D. nasuta subgroup in Papua New Guinea about 3.5 Mya. The ancestral population spread to the north and when it reached Borneo, it diversified sequentially into the kohkoa complex, D. s. bilimbata, and D. kepulauana. About 1 Mya, another radiation occurred when the ancestral populations reached the Indo-China Peninsula, forming the albomicans complex. Discrepancy between morphological groupings and phylogenetic results suggests that the male morphological traits may not be orthologous.
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.
Platanthera is one of the largest genera of temperate orchids in the Holarctic and exemplifies a lineage that has adaptively radiated into diverse habitats within North America, Asia, Europe, North Africa, Borneo, and Sarawak. Major centers of diversity in this genus are North America and eastern Asia. Despite its diversity, a thorough phylogenetic hypothesis for the genus is lacking because no studies have yet sampled taxa exhaustively or developed a robust molecular toolkit. While there is strong evidence that suggests monophyly of subgenus Limnorchis, most taxa in this group have not been included in a phylogenetic analysis. In this study, we developed a new toolkit for Platanthera consisting of genomic information from 617 low-copy nuclear loci. Using a targeted enrichment approach, we collected high-throughput sequence data in 23 accessions of nine of the 12 diploid species of subgenus Limnorchis and outgroup species across Platanthera. A maximum likelihood analysis resolved a strongly supported monophyletic clade for subgenus Limnorchis. Ancestral biogeographic reconstruction indicated that subgenus Limnorchis originated in western North America ca. 3-4.5 Mya from an ancestor that was widespread in western North America and eastern Asia and subsequently diversified in western North America, followed by dispersal of some species to eastern North America. Our results indicate complex biogeographic connections between Asia and North America, and therefore it suggests that Platanthera is a suitable system to test biogeographic hypotheses over time and space in the Holarctic. Our results are also expected to facilitate further study of diversification and biogeographic spread across Platanthera and lay the groundwork for understanding independent origins, biogeography, and morphological diversification of polyploid species within subgenus Limnorchis.
Myrteae (c. 2500 species; 51 genera) is the largest tribe of Myrtaceae and an ecologically important groups of angiosperms in the Neotropics. Systematic relationships in Myrteae are complex, hindering conservation initiatives and jeopardizing evolutionary modelling. A well-supported and robust phylogenetic hypothesis was here targeted towards a comprehensive understanding of the relationships within the tribe. The resultant topology was used as a base for key evolutionary analyses such as age estimation, historical biogeography and diversification rate patterns. One nuclear (ITS) and seven chloroplast (psbA-trnH, matK, ndhF, trnl-trnF, trnQ-rps16, rpl16 and rpl32-trnL) DNA regions for 115 taxa representing 46 out of the 51 genera in the tribe were accessed and analysed using maximum likelihood and Bayesian inference tools for phylogenetic reconstruction. Dates of diversification events were estimated and contrasted using two distinct fossil sets (macro and pollen) in BEAST. The subsequent dated phylogenies were compared and analysed for biogeographical patterns using BioGeoBEARS and diversification rates using BAMM. Myrteae phylogeny presents strong statistical support for three major clades within the tribe: Australasian group, Myrtus group and Main Neotropical Lineage. Dating results from calibration using macrofossil are an average of 20 million years older and show an early Paleocene origin of Myrteae, against a mid-Eocene one from the pollen fossil calibration. Biogeographic analysis shows the origin of Myrteae in Zealandia in both calibration approaches, followed by a widespread distribution throughout the still-linked Gondwana continents and diversification of Neotropical endemic lineages by later vicariance. Best configuration shift indicates three points of acceleration in diversification rates, all of them occurring in the Main Neotropical Lineage. Based on the reconstructed topology, several new taxonomic placements were recovered, including: the relative position of Myrtus communis, the placement of the Blepharocalyx group, the absence of generic endemism in the Caribbean, and the paraphyletism of the former Pimenta group. Distinct calibration approaches affect biogeography interpretation, increasing the number of necessary long distance dispersal events in the topology with older nodes. It is hypothesised that biological intrinsic factors such as modifications of embryo type and polyploidy might have played a role in accelerating shifts of diversification rates in Neotropical lineages. Future perspectives include formal subtribal classification, standardization of fossil calibration approaches and better links between diversification shifts and trait evolution.
The increased rate at which complete mitogenomes are being sequenced and their increasing use for phylogenetic studies have resulted in a bioinformatic bottleneck in preparing and utilising such data for phylogenetic analysis. Hence, we present MitoPhAST, an automated tool that (1) identifies annotated protein-coding gene features and generates a standardised, concatenated and partitioned amino acid alignment directly from complete/partial GenBank/EMBL-format mitogenome flat files, (2) generates a maximum likelihood phylogenetic tree using optimised protein models and (3) reports various mitochondrial genes and sequence information in a table format. To demonstrate the capacity of MitoPhAST in handling a large dataset, we used 81 publicly available decapod mitogenomes, together with eight new complete mitogenomes of Australian freshwater crayfishes, including the first for the genus Gramastacus, to undertake an updated test of the monophyly of the major groups of the order Decapoda and their phylogenetic relationships. The recovered phylogenetic trees using both Bayesian and ML methods support the results of studies using fragments of mtDNA and nuclear markers and other smaller-scale studies using whole mitogenomes. In comparison to the fragment-based phylogenies, nodal support values are generally higher despite reduced taxon sampling suggesting there is value in utilising more fully mitogenomic data. Additionally, the simple table output from MitoPhAST provides an efficient summary and statistical overview of the mitogenomes under study at the gene level, allowing the identification of missing or duplicated genes and gene rearrangements. The finding of new mtDNA gene rearrangements in several genera of Australian freshwater crayfishes indicates that this group has undergone an unusually high rate of evolutionary change for this organelle compared to other major families of decapod crustaceans. As a result, freshwater crayfishes are likely to be a useful model for studies designed to understand the evolution of mtDNA rearrangements. We anticipate that our bioinformatics pipeline will substantially help mitogenome-based studies increase the speed, accuracy and efficiency of phylogenetic studies utilising mitogenome information. MitoPhAST is available for download at https://github.com/mht85/MitoPhAST.
The infraorder Anomura consists of a morphologically and ecologically heterogeneous group of decapod crustaceans, and has attracted interest from taxonomists for decades attempting to find some order out of the seemingly chaotic diversity within the group. Species-level diversity within the Anomura runs the gamut from the "hairy" spindly-legged yeti crab found in deep-sea hydrothermal vent environments to the largest known terrestrial invertebrate, the robust coconut or robber crab. Owing to a well-developed capacity for parallel evolution, as evidenced by the occurrence of multiple independent carcinization events, Anomura has long tested the patience and skill of both taxonomists attempting to find order, and phylogeneticists trying to establish stable hypotheses of evolutionary inter-relationships. In this study, we performed genome skimming to recover the mitogenome sequences of 12 anomuran species including the world's largest extant invertebrate, the robber crab (Birgus latro), thereby over doubling these resources for this group, together with 8 new brachyuran mitogenomes. Maximum-likelihood (ML) and Bayesian-inferred (BI) phylogenetic reconstructions based on amino acid sequences from mitogenome protein-coding genes provided strong support for the monophyly of the Anomura and Brachyura and their sister relationship, consistent with previous studies. The majority of relationships within families were supported and were largely consistent with current taxonomic classifications, whereas many relationships at higher taxonomic levels were unresolved. Nevertheless, we have strong support for a polyphyletic Paguroidea and recovered a well-supported clade of a subset of paguroids (Diogenidae + Coenobitidae) basal to all other anomurans, though this requires further testing with greater taxonomic sampling. We also introduce a new feature to the MitoPhAST bioinformatics pipeline (https://github.com/mht85/MitoPhAST) that enables the extraction of mitochondrial gene order (MGO) information directly from GenBank files and clusters groups based on common MGOs. Using this tool, we compared MGOs across the Anomura and Brachyura, identifying Anomura as a taxonomic "hot spot" with high variability in MGOs among congeneric species from multiple families while noting the broad association of highly-rearranged MGOs with several anomuran lineages inhabiting extreme niches. We also demonstrate the value of MGOs as a source of novel synapomorphies for independently reinforcing tree-based relationships and for shedding light on relationships among challenging groups such as the Aegloidea and Lomisoidea that were unresolved in phylogenetic reconstructions. Overall, this study contributes a substantial amount of new genetic material for Anomura and attempts to further resolve anomuran evolutionary relationships where possible based on a combination of sequence and MGO information. The new feature in MitoPhAST adds to the relatively limited number of bioinformatics tools available for MGO analyses, which can be utilized widely across animal groups.
In the species groups related to Diphasiastrum multispicatum and D. veitchii, hybridization was investigated in samples from northern and southern Vietnam and the island of Taiwan, including available herbarium specimens from southeast Asia. The accessions were analyzed using flow cytometry (living material only), Sanger sequencing and multiplexed inter-simple sequence repeat genotyping by sequencing. We detected two cases of ancient hybridization involving different combinations of parental species; both led via subsequent duplication to tetraploid taxa. A cross D. multispicatum × D. veitchii from Malaysia represents D. wightianum, a tetraploid taxon according to reported DNA measurements of dried material (genome formulas MM, VV and MMVV, respectively). The second case involves D. veitchii (genome formula VV) and an unknown diploid parent (genome formula XX). Three hybridogenous taxa (genome formulas VVX, VVXX, VVVX) were discernable by a combination of flow cytometry and molecular data. Taxon I (VVX, three clones found in Taiwan island) is apparently triploid. Taxon II represents another genetically diverse and sexual tetraploid species (VVXX) and can be assigned to D. yueshanense, described from Taiwan island but occurring as well in mainland China and Vietnam. Taxon III (VVVX) is as well most likely tetraploid and represented by at least one, more likely two, clones from Taiwan island. Taxa I and III are presumably asexual and new to science. Two independently inherited nuclear markers recombine only within, not between these hybrids, pointing towards reproductive isolation. We present an evolutionary scheme which explains the origin of the hybrids and the evolution of new and fully sexual species by hybridization and subsequent allopolyploidization in flat-branched clubmosses.
Family I.3 lipase is distinguished from other families by the amino acid sequence and secretion mechanism. Little is known about the evolutionary process driving these differences. This study attempt to understand how the diverse temperature stabilities of bacterial lipases from family I.3 evolved. To achieve that, eighty-three protein sequences sharing a minimum 30% sequence identity with Antarctic Pseudomonas sp. AMS8 lipase were used to infer phylogenetic tree. Using ancestral sequence reconstruction (ASR) technique, the last universal common ancestor (LUCA) sequence of family I.3 was reconstructed. A gene encoding LUCA was synthesized, cloned and expressed as inclusion bodies in E. coli system. Insoluble form of LUCA was refolded using urea dilution method and then purified using affinity chromatography. The purified LUCA exhibited an optimum temperature and pH at 70 ℃ and 10 respectively. Various metal ions increased or retained the activity of LUCA. LUCA also demonstrated tolerance towards various organic solvents in 25% v/v concentration. The finding from this study could support the understanding on temperature and environment during ancient time. In overall, reconstructed ancestral enzymes have improved physicochemical properties that make them suitable for industrial applications and ASR technique can be employed as a general technique for enzyme engineering.
Inferring interfamilial relationships within the eudicot order Ericales has remained one of the more recalcitrant problems in angiosperm phylogenetics, likely due to a rapid, ancient radiation. As a result, no comprehensive time-calibrated tree or biogeographical analysis of the order has been published. Here, we elucidate phylogenetic relationships within the order and then conduct time-dependent biogeographical and diversification analyses by using a taxon and locus-rich supermatrix approach on one-third of the extant species diversity calibrated with 23 macrofossils and two secondary calibration points. Our results corroborate previous studies and also suggest several new but poorly supported relationships. Newly suggested relationships are: (1) holoparasitic Mitrastemonaceae is sister to Lecythidaceae, (2) the clade formed by Mitrastemonaceae + Lecythidaceae is sister to Ericales excluding balsaminoids, (3) Theaceae is sister to the styracoids + sarracenioids + ericoids, and (4) subfamilial relationships with Ericaceae suggest that Arbutoideae is sister to Monotropoideae and Pyroloideae is sister to all subfamilies excluding Arbutoideae, Enkianthoideae, and Monotropoideae. Our results indicate Ericales began to diversify 110 Mya, within Indo-Malaysia and the Neotropics, with exchange between the two areas and expansion out of Indo-Malaysia becoming an important area in shaping the extant diversity of many families. Rapid cladogenesis occurred along the backbone of the order between 104 and 106 Mya. Jump dispersal is important within the order in the last 30 My, but vicariance is the most important cladogenetic driver of disjunctions at deeper levels of the phylogeny. We detect between 69 and 81 shifts in speciation rate throughout the order, the vast majority of which occurred within the last 30 My. We propose that range shifting may be responsible for older shifts in speciation rate, but more recent shifts may be better explained by morphological innovation.
This study provides a first description of the phylogeographic patterns and evolutionary history of two species of the mudskipper genus Periophthalmus. These amphibious gobies are distributed throughout the whole Indo-Pacific region and Atlantic coast of Africa, in peritidal habitats of soft-bottom coastal ecosystems. Three sequence datasets of two widely distributed species, Periophthalmus argentilineatus and P. kalolo, were obtained by amplifying and sequencing two mtDNA markers (D-loop and 16S rDNA) and the nDNA rag1 region. The three datasets were then used to perform phylogeographic, demographic and population genetic analyses. Our results indicate that tectonic events and past climatic oscillations strongly contributed to shape present genetic differentiation, phylogeographic and demographic patterns. We found support for the monophyly of P. kalolo, and only shallow genetic differentiation between East-African and Indo-Malayan populations of this species. However, our collections of the morphospecies P. argentilineatus include three molecularly distinct lineages, one of them more closely related to P. kalolo. The presence of Miocenic timings for the most recent common ancestors of some of these morphologically similar clades, suggests the presence of strong stabilising selection in mudskippers' habitats. At population level, demographic analyses and palaeoecological records of mangrove ecosystems suggest that Pleistocene bottlenecks and expansion plus secondary contact events of the studied species were associated with recurrent sea transgressions during interglacials, and sea regressions or stable regimes during glacials, respectively.
Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.
As opposed to angiosperms, moss species richness is similar among tropical regions of the world, in line with the hypothesis that tropical bryophytes are extremely good dispersers. Here, we reconstructed the phylogeny of the pantropical moss genus Pelekium to test the hypothesis that high migration rates erase any difference in species richness among tropical regions. In contrast with this hypothesis, several species considered to have a pantropical range were resolved as a complex of species with a strong geographic structure. Consequently, a significant phylogeographical signal was found in the data, evidencing that cladogenetic diversification within regions takes place at a faster rate than intercontinental migration. The shape of the Pelekium phylogeny, along with the selection of a constant-rate model of diversification among species in the genus, suggests, however, that the cladogenetic speciation patterns observed in Pelekium are not comparable to some of the spectacular examples of tropical radiations reported in angiosperms. Rather, the results presented here point to the constant accumulation of diversity through time in Pelekium. This, combined with evidence for long-distance dispersal limitations in the genus, suggests that the similar patterns of species richness among tropical areas are better explained in terms of comparable rates of diversification across tropical regions than by the homogenization of species richness by recurrent migrations.
Nepenthaceae is one of the largest carnivorous plant families and features ecological and morphological adaptations indicating an impressive adaptive radiation. However, investigation of evolutionary and taxonomic questions is hindered by poor phylogenetic understanding, with previous molecular studies based on limited loci and taxa. We use high-throughput sequencing with a target-capture methodology based on a 353-loci, probe set to recover sequences for 197 samples, representing 151 described or putative Nepenthes species. Phylogenetic analyses were performed using supermatrix and maximum quartet species tree approaches. Our analyses confirm five Western outlier taxa, followed by N. danseri, as successively sister to the remainder of the group. We also find mostly consistent recovery of two major Southeast Asian clades. The first contains common or widespread lowland species plus a Wallacean-New Guinean clade. Within the second clade, sects. Insignes and Tentaculatae are well supported, while geographically defined clades representing Sumatra, Indochina, Peninsular Malaysia, Palawan, Mindanao and Borneo are also consistently recovered. However, we find considerable conflicting signal at the site and locus level, and often unstable backbone relationships. A handful of Bornean taxa are inconsistently placed and require further investigation. We make further suggestions for a modified infra-generic classification of genus Nepenthes.
Phylogeographic inference has provided extensive insight into the relative roles of geographical isolation and ecological processes during evolutionary radiations. However, the importance of cross-lineage admixture in facilitating adaptive radiations is increasingly being recognised, and suggested as a main cause of phylogenetic uncertainty. In this study, we used a double digest RADseq protocol to provide a high resolution (~4 Million bp) nuclear phylogeny of the Delphininae. Phylogenetic resolution of this group has been especially intractable, likely because it has experienced a recent species radiation. We carried out cross-lineage reticulation analyses, and tested for several sources of potential bias in determining phylogenies from genome sampling data. We assessed the divergence time and historical demography of T. truncatus and T. aduncus by sequencing the T. aduncus genome and comparing it with the T. truncatus reference genome. Our results suggest monophyly for the genus Tursiops, with the recently proposed T. australis species falling within the T. aduncus lineage. We also show the presence of extensive cross-lineage gene flow between pelagic and European coastal ecotypes of T. truncatus, as well as in the early stages of diversification between spotted (Stenella frontalis; Stenella attenuata), spinner (Stenella longirostris), striped (Stenella coeruleoalba), common (Delphinus delphis), and Fraser's (Lagenodelphis hosei) dolphins. Our study suggests that cross-lineage gene flow in this group has been more extensive and complex than previously thought. In the context of biogeography and local habitat dependence, these results improve our understanding of the evolutionary processes determining the history of this lineage.
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.
The genus Rana, notably diversified in Oriental regions from China to Southeast Asia, includes a group of cascade frogs assigned to subgenera Odorrana and Eburana. Among them, R. ishikawae and the R. narina complex represent the northernmost members occurring from Taiwan to the Ryukyu Archipelago of Japan. Relationships of these frogs with the continental members, as well as the history of their invasions to islands, have been unclear. The taxonomic status of Odorrana and related genera varies among authors and no phylogenetic reassessment has been done. Using partial sequences of mitochondrial 12S and 16S rRNA genes, we estimated phylogenetic relationships among 17 species of the section Hylarana including Odorrana and Eburana, and related species from the Ryukyus, Taiwan, China, Thailand, Malaysia, and Indonesia. We estimate that (1) Odorrana is monophyletic and encompasses species of Eburana and R. hosii, which is now placed in Chalcorana, (2) the ancestor of R. ishikawae separated from other Rana in the middle to late Miocene prior to its entry to the Ryukyu Archipelago, (3) the ancestor of the R. narina complex later diversified in continental Asia, and invaded the Ryukyu Archipelago through Taiwan, (4) the R. narina complex attained its current distribution within the Ryukyus through niche segregations, and (5) vicariance of R. hosii between Malay Peninsula and Borneo occurred much later than the divergence events in the R. narina complex. Current subgeneric classification of Rana, at least of Southeast Asian members, requires full reassessment in the light of phylogenetic relationships.
A fanged frog Limnonectes kuhlii was once thought to be wide-ranging in Southeast Asia, but is now confined to its type locality Java through recent phylogenetic studies, which clarified heterospecific status of non-Javanese populations, and monophyly of Bornean populations. However, large genetic differences among Bornean populations suggest occurrence of cryptic species, which we test using dense geographic sampling. We estimated the phylogenetic relationships among samples of Bornean populations together with their putative relatives from the continental Southeast Asia, using 2517bp sequences of the 12S rRNA, tRNA(val), and 16S rRNA of mitochondrial DNA, and 2367bp sequences of the NCX1, POMC, and RAG1 of nuclear genes. In the mtDNA trees, Bornean L. kuhlii-like frogs formed a monophyletic group split into 18 species lineages including L. hikidai, with the deepest phylogenetic split separating L. cintalubang from the remaining species. Almost all of these lineages co-occur geographically, and two to three lineages were found syntopically in each locality. Co-occurrence of more than one lineage may be maintained by differential morphology and microhabitat selection. These syntopic lineages should be regarded as distinct species. Our results clearly indicate that taxonomic revision is urgent to clarify many evolutionary problems of Bornean L. kuhlii-like frogs.