Diversification of rhacophorid frogs provides evidence for accelerated faunal exchange between India and Eurasia during the Oligocene

A new species of Raorchestes (Anura, Rhacophoridae) from Yunnan Province, China

A new bush frog species is described from Yunnan, China, based on phylogenetic analyses, species delimitation analyses, and morphological comparisons. Raorchesteshekouensissp. nov. is distinguished from all other congeners by a combination of 11 morphological characters. The new species brings the current number of Raorchestes species in China to ten, nine of which are distributed in Yunnan. Molecular analyses supported an unnamed lineage previously recorded as "Raorchestesgryllus" in northern Vietnam. Further studies including additional samples are necessary to clarify the species diversity and boundaries of Raorchestes in China and Indochina.

Integrated phylogenetic analyses reveal the evolutionary, biogeographic, and diversification history of Asian warty treefrog genus Theloderma (Anura, Rhacophoridae)

Asian warty treefrogs, genus Theloderma, are morphologically variable arboreal frogs endemic to Southeast Asia and Southern China. However, integrated systematic studies are lacking, and knowledge of the genus in terms of diversity, origin, and historical diversification remains limited. To address these knowledge gaps, we used three mitochondrial and five nuclear gene fragments to reconstruct the Theloderma phylogeny, estimate divergence times, and examine the biogeography of the genus. Phylogenetic and species delimitation analyses suggest that the genus Theloderma comprises three major clades corresponding to two subgenera and seven species groups, and mPTP identified at least 12 putative cryptic species, suggesting that species diversity has been underestimated. Biogeographic analyses indicated that most recent common ancestor of Theloderma originated in the Indochina Peninsula during the Middle Oligocene (ca. 27.77 Ma) and the splitting of Clade A to C occurred in the Late Oligocene (ca. 23.55-25.57 Ma). Current biogeographic patterns result from two distinct processes: in situ diversification in the Indochina Peninsula and dispersal in multiple areas, namely southward dispersal to the Malay Peninsula and Borneo, northeastward dispersal to Southern China, northward dispersal to the Himalayas, and dispersal from Southern China to the Indochina Peninsula. Ancestral character reconstruction suggests that the ancestor of Theloderma may have possessed a small body size, rough dorsal skin, and absence of vomerine teeth and hand webbing, and that these four characters have undergone multiple evolutions. Principal component analysis based on eight bioclimatic variables did not clearly distinguish the three major clades of Theloderma, suggesting that species in these clades may occupy similar climatic ecological niches. Our research highlights the importance of orogeny and paleoclimatic changes, in shaping amphibian biodiversity in mountain ecosystems.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

An updated floristic map of the world

Floristic regions reflect the geographic organization of floras and provide essential tools for biological studies. Previous global floristic regions are generally based on floristic endemism, lacking a phylogenetic consideration that captures floristic evolution. Moreover, the contribution of tectonic dynamics and historical and current climate to the division of floristic regions remains unknown. Here, by integrating global distributions and a phylogeny of 12,664 angiosperm genera, we update global floristic regions and explore their temporal changes. Eight floristic realms and 16 nested sub-realms are identified. The previously-defined Holarctic, Neotropical and Australian realms are recognized, but Paleotropical, Antarctic and Cape realms are not. Most realms have formed since Paleogene. Geographic isolation induced by plate tectonics dominates the formation of floristic realms, while current/historical climate has little contribution. Our study demonstrates the necessity of integrating distributions and phylogenies in regionalizing floristic realms and the interplay of macroevolutionary and paleogeographic processes in shaping regional floras.

Divergence time of mites of the family Laelapidae based on mitochondrial barcoding region

Using the mitochondrial barcoding region to correlate research with 58 species in 19 genera of the family Laelapidae with the aim of determining the origin, phylogenetic relationships, and biogeographic historical distribution characteristics of mites in the family Laelapidae. Phylogenetic trees were obtained using Bayesian inference (BI) and Maximum-likelihood (ML) methods, based on three fossil records calibrated as molecular clock nodes, to estimate the divergence time of mites in the family Laelapidae as well as to apply Dispersal-Extinction-Cladogenesis (DEC) analyses to obtain biogeographic history inferences. The result showed species of the genera Hyperlaelaps and Haemolaelaps and some species of the genus Androlaelaps in the family Laelapidae were divided into clades of the genus Laelaps in both the BI and ML trees. Divergence time estimates and biogeographic history analysis revealed that the family Laelapidae likely diverged from other taxa during the Middle Jurassic (ca. 156.73 Mya), with Asia considered the most likely ancestral region for the family Laelapidae. Species of various genera began to undergo massive diversification events during the Cenozoic Tertiary. The results suggest that some genera in the family Laelapidae need to be re-defined or new genera need to be established; the Late Cretaceous to Late Neogene warm period would have promoted the divergence and expansion of species in the family Laelapidae. The divergence and dispersal of the family Laelapidae species is most likely a joint response to the continued northward drift of the Indian plate away from the Gondwana paleo-continent and gradually closer to Asia during the Late Cretaceous and the geological activity of the Tibetan Plateau during the Cenozoic Tertiary. The results strengthen our understanding of the origin and evolution of species in the family Laelapidae.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

First record of Theloderma khoii Ninh, Nguyen, Nguyen, Hoang, Siliyavong, Nguyen, Le, Le &amp; Ziegler, 2022 from China, with confirmation of Rhacophorus orlovi Ziegler &amp; Köhler, 2001 in China (Anura, Rhacophoridae)

We report the first country record of Theloderma khoii Ninh, Nguyen, Nguyen, Hoang, Siliyavong, Nguyen, Le, Le & Ziegler, 2022 from China based on a specimen collected from Wenshan Prefecture, southeastern Yunnan, China. Morphologically, the specimen from Wenshan Prefecture agrees well with the type specimens of T. khoii from Vietnam except for having a smaller body size, and phylogenetically clustered with the type specimens of T. khoii from Vietnam. In addition, we confirm the distribution of Rhacophorus orlovi Ziegler & Köhler, 2001 in China based on three specimens collected from Honghe Prefecture, southern Yunnan, China. Morphologically, the specimens from Honghe Prefecture agree well with the type specimens of R. orlovi from Vietnam except for having a slightly larger body size in adult females, and phylogenetically clustered with the specimens of R. orlovi from Vietnam (including the type locality of R. orlovi).

The stepwise Indian-Eurasian collision and uplift of the Himalayan-Tibetan plateau drove the diversification of high-elevation Scytodes spiders

The Cenozoic Indian-Eurasian collision and uplift of the Himalayan-Tibetan Plateau (HTP) are among the most important geological events in the world. They have affected the diversification of regional biota of many taxonomic groups on and around the HTP. However, the exact timing and model of the collision and uplift events and speciation on and around the HTP are still in debate. The Himalayas group of Scytodes spitting spiders (Araneae: Scytodidae) are distributed at high elevations of the HTP and northern Indochina. Here, we reconstruct a dated molecular phylogeny for pan-Himalayan Scytodes spiders, including the Himalayas group, with full geographical sampling of the species from the HTP and Indochina. We test a hypothesis to explain that the rich montane biodiversity of the region is uplift-driven diversification-that orogeny drives rapid in situ speciation of the resident Scytodes lineages. Our findings revealed that the separation of the Himalayas clade from the Myanmar clade took place during the middle Oligocene, reflecting the final collision of India with Eurasia. The deep divergences among three clades (the Himalayas, the Myanmar and the Indochina clades) occurred from the middle Eocene to the middle Oligocene, corresponding to two early uplift events of the HTP. The evolutionary split between the Himalayas + Myanmar and Indochina clades were simultaneous with the rapid lateral extrusion of Indochina by the initial Himalayan uplift around the Eocene. This study highlights the importance of the diversification of dispersal-limited, high-elevation invertebrates as independent lines of evidence to reflect key tectonic events in the Himalayan-Indochina region, supporting the stepwise model for the Indian-Eurasian collision and uplift of HTP.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

Simple Summary

Duplicated control regions have been reported several times in the tree frog family Rhacophoridae, and previous studies have mostly relied on sequence analysis to reconstruct their evolution. This is the first study to employ a phylogenetic method to demonstrate the existence of concerted and parallel evolution succinctly and intuitively in the duplicated control regions of the family Rhacophoridae. Phylogenetic relationships were also used to illustrate the parallel evolution of ATP8 loss of function in the genus Polypedates. In general, this study elucidated the evolutionary patterns and pathways of mitochondrial gene rearrangement of the family Rhacophoridae from a phylogenetic perspective, which aids in understanding the evolutionary history of this fascinating tree frog taxon from a molecular evolution standpoint.

Abstract

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

The role of climate and islands in species diversification and reproductive-mode evolution of Old World tree frogs

Large diversifications of species are known to occur unevenly across space and evolutionary lineages, but the relative importance of their driving mechanisms, such as climate, ecological opportunity and key evolutionary innovations (KEI), remains poorly understood. Here, we explore the remarkable diversification of rhacophorid frogs, which represent six percent of global amphibian diversity, utilize four distinct reproductive modes, and span a climatically variable area across mainland Asia, associated continental islands, and Africa. Using a complete species-level phylogeny, we find near-constant diversification rates but a highly uneven distribution of species richness. Montane regions on islands and some mainland regions have higher phylogenetic diversity and unique assemblages of taxa; we identify these as cool-wet refugia. Starting from a centre of origin, rhacophorids reached these distant refugia by adapting to new climatic conditions (‘niche evolution’-dominant), especially following the origin of KEIs such as terrestrial reproduction (in the Late Eocene) or by dispersal during periods of favourable climate (‘niche conservatism’-dominant).

By examining climate, geographical and phylogenetic data, the diversification and evolution of rhacophorid frogs is explored

The biogeography of bent-toed geckos, Cyrtodactylus (Squamata: Gekkonidae)

The gekkonid genus Cyrtodactylus is the third largest vertebrate genus on the planet with well over 300 species that range across at least eight biogeographic regions from South Asia to Melanesia. The ecological and morphological plasticity within the genus, has contributed to its ability to disperse across ephemeral seaways, river systems, basins, land bridges, and mountain ranges-followed by in situ diversification within specific geographic areas. Ancestral ranges were reconstructed on a mitochondrial phylogeny with 346 described and undescribed species from which it was inferred that Cyrtodactylus evolved in a proto-Himalaya region during the early Eocene. From there, it dispersed to what is currently Indoburma and Indochina during the mid-Eocene-the latter becoming the first major center of origin for the remainder of the genus that seeded dispersals to the Indian subcontinent, Papua, and Sundaland. Sundaland became a second major center of radiation during the Oligocene and gave rise to a large number of species that radiated further within Sundaland and dispersed to Wallacea, the Philippines, and back to Indochina. One Papuan lineage dispersed west to recolonize and radiate in Sundaland. Currently, Indochina and Sundaland still harbor the vast majority of species of Cyrtodactylus.

Diversity, distribution and molecular species delimitation in frogs and toads from the Eastern Palaearctic

Biodiversity analyses can greatly benefit from coherent species delimitation schemes and up-to-date distribution data. In this article, we have made the daring attempt to delimit and map described and undescribed lineages of anuran amphibians in the Eastern Palaearctic (EP) region in its broad sense. Through a literature review, we have evaluated the species status considering reproductive isolation and genetic divergence, combined with an extensive occurrence dataset (nearly 85k localities). Altogether 274 native species from 46 genera and ten families were retrieved, plus eight additional species introduced from other realms. Independent hotspots of species richness were concentrated in southern Tibet (Medog County), the circum-Sichuan Basin region, Taiwan, the Korean Peninsula and the main Japanese islands. Phylogeographic breaks responsible for recent in situ speciation events were shared around the Sichuan Mountains, across Honshu and between the Ryukyu Island groups, but not across shallow water bodies like the Yellow Sea and the Taiwan Strait. Anuran compositions suggested to restrict the zoogeographical limits of the EP to East Asia. In a rapidly evolving field, our study provides a checkpoint to appreciate patterns of species diversity in the EP under a single, spatially explicit, species delimitation framework that integrates phylogeographic data in taxonomic research.

A phylogeographical framework for Zhangixalus gliding frogs, with insight on their plasticity of nesting behaviour

Old World tree frogs from the family Rhacophoridae, one of the most species-rich groups of amphibians worldwide, are becoming a model in ecological and evolutionary research, notably for their tremendous diversity of breeding systems. In this study, we provide the most comprehensive temporally and spatially explicit phylogeographical framework for the tropical neo-genus Zhangixalus, with a specific focus on the evolution of inter- and intraspecific polymorphism in their egg-laying strategies. We combined ~500 sequences spanning ~4 kb of mitochondrial genes, representative of nearly all known Zhangixalus species from 373 localities. We recovered a robust timetree featuring ~60 phylogeographical lineages assigned to 36 taxa, illustrating a vast under-appreciation of the diversity of the genus. Patterns of diversification were intimately linked to climatic trends since the Mid-Miocene, leading to a remarkable biogeographical scenario for the onset of Taiwanese diversity. Finally, in contrast to the macro-evolution of conserved breeding strategies in rhacophorids, the striking variation of nesting behaviour documented across Zhangixalus species, whereby foamy egg nests are alternatively fossorial, deposited on the edge of the water or suspended on the overhanging vegetation, seems to have evolved rapidly from arboreal ancestors in all major clades, implying micro-evolutionary processes and/or phenotypic plasticity.

Molecular species delimitation of shrub frogs of the genus Pseudophilautus (Anura, Rhacophoridae)

Sri Lanka is an amphibian hotspot of global significance. Its anuran fauna is dominated by the shrub frogs of the genus Pseudophilautus. Except for one small clade of four species in Peninsular India, these cool-wet adapted frogs, numbering some 59 extant species, are distributed mainly across the montane and lowland rain forests of the island. With species described primarily by morphological means, the diversification has never yet been subjected to a molecular species delimitation analysis, a procedure now routinely applied in taxonomy. Here we test the species boundaries of Pseudophilautus in the context of the phylogenetic species concept (PSC). We use all the putative species for which credible molecular data are available (nDNA-Rag-1; mt-DNA- 12S rRNA, 16S rRNA) to build a well resolved phylogeny, which is subjected to species delimitation analyses. The ABGD, bPTP, mPTP and bGMYC species delimitation methods applied to the 16S rRNA frog barcoding gene (for all species), 12S rRNA and Rag-1 nDNA grouped P. procax and P. abundus; P. hallidayi and P. fergusonianus; P. reticulatus and P. pappilosus; P. pleurotaenia and P. hoipolloi; P. hoffmani and P. asankai; P. silvaticus and P. limbus; P. dilmah and P. hankeni; P. fulvus and P. silus.. Surprisingly, all analyses recovered 14 unidentified potential new species as well. The geophylogeny affirms a distribution across the island's aseasonal 'wet zone' and its three principal hill ranges, suggestive of allopatric speciation playing a dominant role, especially between mountain masses. Among the species that are merged by the delimitation analyses, a pattern leading towards a model of parapatric speciation emerges-ongoing speciation in the presence of gene flow. This delimitation analysis reinforces the species hypotheses, paving the way to a reasonable understanding of Sri Lankan Pseudophilautus, enabling both deeper analyses and conservation efforts of this remarkable diversification. http://zoobank.org/urn:lsid:zoobank.org:pub:DA869B6B-870A-4ED3-BF5D-5AA3F69DDD27.

A tree of leaves: Phylogeny and historical biogeography of the leaf insects (Phasmatodea: Phylliidae)

The insect order Phasmatodea is known for large slender insects masquerading as twigs or bark. In contrast to these so-called stick insects, the subordinated clade of leaf insects (Phylliidae) are dorso-ventrally flattened and therefore resemble leaves in a unique way. Here we show that the origin of extant leaf insects lies in the Australasian/Pacific region with subsequent dispersal westwards to mainland Asia and colonisation of most Southeast Asian landmasses. We further hypothesise that the clade originated in the Early Eocene after the emergence of angiosperm-dominated rainforests. The genus Phyllium to which most of the ~100 described species pertain is recovered as paraphyletic and its three non-nominate subgenera are recovered as distinct, monophyletic groups and are consequently elevated to genus rank. This first phylogeny covering all major phylliid groups provides the basis for future studies on their taxonomy and a framework to unveil more of their cryptic and underestimated diversity.

An integrative approach to infer systematic relationships and define species groups in the shrub frog genus Raorchestes, with description of five new species from the Western Ghats, India

The genus Raorchestes is a large radiation of Old World tree frogs for which the Western Ghats in Peninsular India is the major center for origin and diversification. Extensive studies on this group during the past two decades have resolved long-standing taxonomic confusions and uncovered several new species, resulting in a four-fold increase in the number of known Raorchestes frogs from this region. Our ongoing research has revealed another five new species in the genus, formally described as Raorchestes drutaahu sp. nov., Raorchestes kakkayamensis sp. nov., Raorchestes keirasabinae sp. nov., Raorchestes sanjappai sp. nov., and Raorchestes vellikkannan sp. nov., all from the State of Kerala in southern Western Ghats. Based on new collections, we also provide insights on the taxonomic identity of three previously known taxa. Furthermore, since attempts for an up-to-date comprehensive study of this taxonomically challenging genus using multiple integrative taxonomic approaches have been lacking, here we review the systematic affinities of all known Raorchestes species and define 16 species groups based on evidence from multi-gene (2,327 bp) phylogenetic analyses, several morphological characters (including eye colouration and pattern), and acoustic parameters (temporal and spectral properties, as well as calling height). The results of our study present novel insights to facilitate a better working taxonomy for this rather speciose and morphologically conserved radiation of shrub frogs. This will further enable proper field identification, provide momentum for multi-disciplinary studies, as well as assist conservation of one of the most colourful and acoustically diverse frog groups of the Western Ghats biodiversity hotspot.

Geographical range evolution of the genus Polypedates (Anura: Rhacophoridae) from the Oligocene to present

Currently, the genus Polypedates comprises 26 species distributed in South, Southeast, and East Asia. Because of their relatively low dispersal capability and intolerance to seawater, this genus is ideal for the study of terrestrial range evolution that extends into the island archipelagos of southeastern Asia. In this study, based on data compiled for Polypedatesfrom previous studies and partial mitochondrial and nuclear genes collected in this study, we performed systematic biogeographical analysis. We confirmed a Sundaland origin for the extant genus and showed northward dispersal into mainland Southeast Asia and Asia, which coincided with the timing of paleoclimatic change from the Oligocene to Middle Miocene. Climate fluctuations had a profound impact on species diversification within the genus Polypedates. Furthermore, the Red River did not mediate species exchange between Southeast Asia and mainland Asia until the end of the Miocene, with the sudden onset of northward dispersal in several clades independently at that time. Alternatively, the lineage of widespread insular P. leucomystaxstrongly supports the hypothesis of terrestrial connection between island archipelagos of Southeast Asia during the Mid-Pleistocene paleoclimate fluctuations. Our biogeographical analysis also supports the recent introduction of P. leucomystax to the Philippines and Ryukyus, as previously suggested.

Tropical Niche Conservatism Explains the Eocene Migration from India to Southeast Asia in Ochyroceratid Spiders

Biological migrations between India and Southeast (SE) Asia provide an ideal system for exploring the effects of geology and climate on species ranges. Geologists have confirmed that the direct collision between India and Eurasia occurred in the Early Eocene, but most migrations occurred between the Indian subcontinent and SE Asia rather than the former and the southern margin of Eurasia. To explain this seemingly paradoxical disconnect between the routes of plate movement and biological migration, we studied the evolutionary history of the tropical spider family Ochyroceratidae based on 101 globally distributed species. We infer a robust dated phylogeny using both transcriptomic data and a data set of classical markers and relate these to biogeographic and climatic analyses. Our results indicate that the monophyly of Ochyroceratidae is strongly supported, and the divergence times suggest a Cretaceous Gondwanan origin of the family. Reconstructed biogeographic histories support a dispersal event from the Indian subcontinent to islands of SE Asia 55-38 Ma. Climatic analyses and the fossil record reveal that ochyroceratids are characterized by a high degree of tropical niche conservatism, and that the ancestor of the Indian and SE Asian clades originated in very warm, wet environments. Early Eocene tropical, perhumid climates in India, and SE Asia may have facilitated ochyroceratid migration, whereas the dry or seasonal climate extending from the eastern coast of China to Central Asia may have acted as a barrier, preventing dispersal. Our analyses suggest that climate plays a more important role than geology in biological migration from the Indian subcontinent to SE Asia, providing new insights into the Indian-Asian biogeographic link. [Biogeography; ecology; geological connections; macroevolution; paleoclimate.].

Out of India, thrice: diversification of Asian forest scorpions reveals three colonizations of Southeast Asia

The 'Out of India' hypothesis is often invoked to explain patterns of distribution among Southeast Asian taxa. According to this hypothesis, Southeast Asian taxa originated in Gondwana, diverged from their Gondwanan relatives when the Indian subcontinent rifted from Gondwana in the Late Jurassic, and colonized Southeast Asia when it collided with Eurasia in the early Cenozoic. A growing body of evidence suggests these events were far more complex than previously understood, however. The first quantitative reconstruction of the biogeography of Asian forest scorpions (Scorpionidae Latreille, 1802: Heterometrinae Simon, 1879) is presented here. Divergence time estimation, ancestral range estimation, and diversification analyses are used to determine the origins, dispersal and diversification patterns of these scorpions, providing a timeline for their biogeographical history that can be summarized into four major events. (1) Heterometrinae diverged from other Scorpionidae on the African continent after the Indian subcontinent became separated in the Cretaceous. (2) Environmental stresses during the Cretaceous-Tertiary (KT) mass extinction caused range contraction, restricting one clade of Heterometrinae to refugia in southern India (the Western Ghats) and Sri Lanka (the Central Highlands). (3) Heterometrinae dispersed to Southeast Asia three times during India's collision with Eurasia, the first dispersal event occurring as the Indian subcontinent brushed up against the western side of Sumatra, and the other two events occurring as India moved closer to Eurasia. (4) Indian Heterometrinae, confined to southern India and Sri Lanka during the KT mass extinction, recolonized the Deccan Plateau and northern India, diversifying into new, more arid habitats after environmental conditions stabilized. These hypotheses, which are congruent with the geological literature and biogeographical analyses of other taxa from South and Southeast Asia, contribute to an improved understanding of the dispersal and diversification patterns of taxa in this biodiverse and geologically complex region.

Target-capture phylogenomics provide insights on gene and species tree discordances in Old World treefrogs (Anura: Rhacophoridae)

Genome-scale data have greatly facilitated the resolution of recalcitrant nodes that Sanger-based datasets have been unable to resolve. However, phylogenomic studies continue to use traditional methods such as bootstrapping to estimate branch support; and high bootstrap values are still interpreted as providing strong support for the correct topology. Furthermore, relatively little attention has been given to assessing discordances between gene and species trees, and the underlying processes that produce phylogenetic conflict. We generated novel genomic datasets to characterize and determine the causes of discordance in Old World treefrogs (Family: Rhacophoridae)-a group that is fraught with conflicting and poorly supported topologies among major clades. Additionally, a suite of data filtering strategies and analytical methods were applied to assess their impact on phylogenetic inference. We showed that incomplete lineage sorting was detected at all nodes that exhibited high levels of discordance. Those nodes were also associated with extremely short internal branches. We also clearly demonstrate that bootstrap values do not reflect uncertainty or confidence for the correct topology and, hence, should not be used as a measure of branch support in phylogenomic datasets. Overall, we showed that phylogenetic discordances in Old World treefrogs resulted from incomplete lineage sorting and that species tree inference can be improved using a multi-faceted, total-evidence approach, which uses the most amount of data and considers results from different analytical methods and datasets.

Evolutionary history of zoogeographical regions surrounding the Tibetan Plateau

The Tibetan Plateau (TP) and surrounding regions have one of the most complex biotas on Earth. However, the evolutionary history of these regions in deep time is poorly understood. Here, we quantify the temporal changes in beta dissimilarities among zoogeographical regions during the Cenozoic using 4,966 extant terrestrial vertebrates and 1,278 extinct mammal genera. We identify ten present-day zoogeographical regions and find that they underwent a striking change over time. Specifically, the fauna on the TP was close to the Oriental realm in deep time but became more similar to the Palearctic realms more recently. The present-day zoogeographical regions generally emerged during the Miocene/Pliocene boundary (ca. 5 Ma). These results indicate that geological events such as the Indo-Asian Collision, the TP uplift, and the aridification of the Asian interior underpinned the evolutionary history of the zoogeographical regions surrounding the TP over different time periods.

A little frog leaps a long way: compounded colonizations of the Indian Subcontinent discovered in the tiny Oriental frog genus Microhyla (Amphibia: Microhylidae)

Frogs of the genus Microhyla include some of the world's smallest amphibians and represent the largest radiation of Asian microhylids, currently encompassing 50 species, distributed across the Oriental biogeographic region. The genus Microhyla remains one of the taxonomically most challenging groups of Asian frogs and was found to be paraphyletic with respect to large-sized fossorial Glyphoglossus. In this study we present a time-calibrated phylogeny for frogs in the genus Microhyla, and discuss taxonomy, historical biogeography, and morphological evolution of these frogs. Our updated phylogeny of the genus with nearly complete taxon sampling includes 48 nominal Microhyla species and several undescribed candidate species. Phylogenetic analyses of 3,207 bp of combined mtDNA and nuDNA data recovered three well-supported groups: the Glyphoglossus clade, Southeast Asian Microhyla II clade (includes M. annectens species group), and a diverse Microhyla I clade including all other species. Within the largest major clade of Microhyla are seven well-supported subclades that we identify as the M. achatina, M. fissipes, M. berdmorei, M. superciliaris, M. ornata, M. butleri, and M. palmipes species groups. The phylogenetic position of 12 poorly known Microhyla species is clarified for the first time. These phylogenetic results, along with molecular clock and ancestral area analyses, show the Microhyla-Glyphoglossus assemblage to have originated in Southeast Asia in the middle Eocene just after the first hypothesized land connections between the Indian Plate and the Asian mainland. While Glyphoglossus and Microhyla II remained within their ancestral ranges, Microhyla I expanded its distribution generally east to west, colonizing and diversifying through the Cenozoic. The Indian Subcontinent was colonized by members of five Microhyla species groups independently, starting with the end Oligocene-early Miocene that coincides with an onset of seasonally dry climates in South Asia. Body size evolution modeling suggests that four groups of Microhyla have independently achieved extreme miniaturization with adult body size below 15 mm. Three of the five smallest Microhyla species are obligate phytotelm-breeders and we argue that their peculiar reproductive biology may be a factor involved in miniaturization. Body size increases in Microhyla-Glyphoglossus seem to be associated with a burrowing adaptation to seasonally dry habitats. Species delimitation analyses suggest a vast underestimation of species richness and diversity in Microhyla and reveal 15-33 undescribed species. We revalidate M. nepenthicola, synonymize M. pulverata with M. marmorata, and provide insights on taxonomic statuses of a number of poorly known species. Further integrative studies, combining evidence from phylogeny, morphology, advertisement calls, and behavior will result in a better systematic understanding of this morphologically cryptic radiation of Asian frogs.

An integrative phylogenomic approach illuminates the evolutionary history of Old World tree frogs (Anura: Rhacophoridae)

Rhacophoridae are one of the most speciose and ecologically diverse families of amphibians. Resolution of their evolutionary relationships is key to understanding the accumulation of biodiversity, yet previous hypotheses based on Sanger sequencing exhibit much discordance amongst generic relationships. This conflict precludes the making of sound macroevolutionary conclusions. Herein, we conduct the first phylogenomic study using broad-scale sampling and sequences of 352 nuclear DNA loci obtained using anchored hybrid enrichment targeted sequencing. The robust time-calibrated phylogenetic hypothesis clarifies several long-disputed relationships and facilitates the testing of evolutionary hypotheses on spatiotemporal diversification and reproductive modes. The major extant lineages of Rhacophoridae appear to have radiated in mainland Asia, and the spatiotemporal process corresponds with several common accumulations of biodiversity in Asia. Analyses do not detect any case of "Out of Himalaya" in Rhacophoridae. All transitions of reproductive modes appear to have evolved in an ordered, gradual sequence associated with gaining independence of standing water for larval development. The different reproductive modes are phylogenetically conserved and the completion of their transitions appear to have occurred over a period of ~30 Ma, which does not fit a pattern of a rapid burst of diversification. Innovations in reproductive modes associate statistically with the uneven distribution of species-richness between clades, where higher diversification is linked to increased terrestrial modes of reproduction. These results strengthen the hypothesis that breeding innovations drive diversification by providing new opportunities for ecological release and dispersion.

Indopolystoma n. gen. (Monogenea, Polystomatidae) with the description of three new species and reassignment of eight known Polystoma species from Asian frogs (Anura, Rhacophoridae)

The polystomes (Monogenea, Polystomatidae) radiated across semi-aquatic tetrapods including all three amphibian orders, freshwater turtles and the hippopotamus. Prior to this study, phylogenetic analyses revealed that the most diverse and widespread genus, Polystoma, was not monophyletic; a lineage comprising four undescribed species from the bladder of Zhangixalus spp. (Rhacophoridae) in Asia occupied a deep phylogenetic position. Regarding vicariance biogeography and molecular dating, the origin of this lineage is correlated with the breakup of Gondwanaland in the Mesozoic period. Based on a Bayesian analysis of four concatenated genes (18S, 28S, COI and 12S) and morphological evidence, one new genus, Indopolystoma n. gen., and three new species, sampled in Japan and China, are described here: Indopolystoma viridi n. sp. from Z. viridis of Japan, Indopolystoma elongatum n. sp. from Z. arboreus of Japan, and Indopolystoma parvum n. sp. from Z. omeimontis of China. Indopolystoma is unique amongst polystome genera infecting anurans by possessing a small haptor relative to the body size, posteriormost marginal hooklet C1 much bigger than hooklets C2–C8 with conspicuous broad blade and guard and a pair of hamuli lacking a deep notch. Eight species of Asian Polystoma, all from rhacophorids, are transferred as Indopolystoma carvirostris (Fan, Li & He, 2008) n. comb., I. hakgalense (Crusz & Ching, 1975) n. comb., I. indicum (Diengdoh & Tandon, 1991) n. comb., I. leucomystax (Zhang & Long, 1987) n. comb., I. mutus (Meng, Song & Ding, 2010) n. comb., I. pingbianensis (Fan, Wang & Li, 2004) n. comb., I. rhacophori (Yamaguti, 1936) n. comb., and I. zuoi (Shen, Wang & Fan, 2013) n. comb.

Phylogeography of Bellamya (Mollusca: Gastropoda: Viviparidae) snails on different continents: contrasting patterns of diversification in China and East Africa

BACKGROUND

Species diversity is determined by both local environmental conditions that control differentiation and extinction and the outcome of large-scale processes that affect migration. The latter primarily comprises climatic change and dynamic landscape alteration. In the past few million years, both Southeast Asia and Eastern Africa experienced drastic climatic and geological oscillations: in Southeast Asia, especially in China, the Tibetan Plateau significantly rose up, and the flow of the Yangtze River was reversed. In East Africa, lakes and rivers experienced frequent range expansions and regressions due to the African mega-droughts. To test how such climatic and geological histories of both regions relate to their respective regional species and genetic diversity, a large scale comparative phylogeographic study is essential. Bellamya, a species rich freshwater snail genus that is widely distributed across China and East Africa, represents a suitable model system to address this question. We sequenced mitochondrial and nuclear DNA for members of the genus from China and used published sequences from Africa and some other locations in Asia to investigate their phylogeny and distribution of genetic diversity.

RESULTS

Our phylogenetic analysis revealed two monophyletic groups, one in China and one in East Africa. Within the Chinese group, Bellamya species show little genetic differentiation. In contrast, we observe fairly deep divergence among the East African lakes with almost every lake possessing its unique clade. Our results show that strong divergence does not necessarily depend on intrinsic characteristics of a species, but rather is related to the landscape dynamics of a region.

CONCLUSION

Our phylogenetic results suggest that the Bellamya in China and East Africa are independent phylogenetic clades with different evolutionary trajectories. The different climate and geological histories likely contributed to the diverging evolutionary patterns. Repeated range expansions and regressions of lakes likely contributed to the great divergence of Bellamya in East Africa, while reversal of the river courses and intermingling of different lineages had an opposite effect on Bellamya diversification in China.

New microhylid frog genus from Peninsular India with Southeast Asian affinity suggests multiple Cenozoic biotic exchanges between India and Eurasia

Anurans in Peninsular India exhibit close biogeographical links with Gondwana as well as Laurasia, often explainable by the geological history of the Indian subcontinent; its breakup from Gondwanan landmasses followed by long isolation that resulted in diversification of endemic lineages, and subsequent land connections with Asia that enabled dispersal of widespread groups. Although widely distributed, the frog subfamily Microhylinae mostly comprises of geographically restricted genera found either in Southeast and East Asia or Peninsular India and Sri Lanka. Here we report a previously unknown microhylid from the Western Ghats in Peninsular India with closest relatives found over 2,000 km away in Southeast Asia. Based on integrated evidence from mitochondrial and nuclear DNA, adult and tadpole morphology, hand musculature, male advertisement call, and geographical distance, we recognize this enigmatic frog as a distinct new species and genus endemic to the Western Ghats. The discovery of Mysticellus franki gen. et sp. nov. and its close evolutionary relationship with the Southeast Asian genus Micryletta also provide insights on the biogeography of Microhylinae. Genus-level divergences within the subfamily suggest multiple Cenozoic biotic exchange events between India and Eurasia, particularly through postulated Eocene land bridges via Southeast Asia prior to accretion of the two landmasses.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

A new species of Kurixalus from western Yunnan, China (Anura, Rhacophoridae)

A new species of the genus Kurixalus (Anura: Rhacophoridae) is described from western Yunnan, China. Genetically the new species, Kurixalusyangisp. n., is closer to Kurixalusnaso than to other known congeners. Morphologically the new species is distinguished from all other known congeners by a combination of the following characters: smaller ratios of head, snout, limbs, IND, and UEW to body size; male body size larger than 30 mm; curved canthus rostralis; weak nuptial pad; brown dorsal color; absence of large dark spots on surface of upper-middle abdomen; presence of vomerine teeth; gold brown iris; single internal vocal sac; serrated dermal fringes along outer edge of limbs; granular throat and chest; rudimentary web between fingers; and presence of supernumerary tubercles and outer metacarpal tubercle.

A new species of the genus Theloderma Tschudi, 1838 (Amphibia: Anura: Rhacophoridae) from Tay Nguyen Plateau, central Vietnam

A new species of small tree frog from a primary montane tropical forest of central Vietnam, Tay Nguyen Plateau, is described based on morphological, molecular, and acoustic evidence. The Golden Bug-Eyed Frog, Theloderma auratum sp. nov., is distinguishable from its congeners and other small rhacophorid species based on a combination of the following morphological attributes: (1) bony ridges on head absent; (2) smooth skin completely lacking calcified warts or asperities; (3) pointed elongated tapering snout; (4) vocal opening in males absent; (5) vomerine teeth absent; (6) males of small body size (SVL 21.8-26.4 mm); (7) head longer than wide; ED/SVL ratio 13%-15%; ESL/SVL ratio 16%-20%; (8) small tympanum (TD/EL ratio 50%-60%) with few tiny tubercles; (9) supratympanic fold absent; (10) ventral surfaces completely smooth; (11) webbing between fingers absent; (12) outer and inner metacarpal tubercles present, supernumerary metacarpal tubercle single, medial, oval in shape; (13) toes half-webbed: I 2-2¼ II 1½-2¾ III 2-3¼ IV 3-1½ V; (14) inner metatarsal tubercle present, oval; outer metatarsal tubercle absent; (15) iris bicolored; (16) dorsal surfaces golden-yellow with sparse golden-orange speckling or reticulations and few small dark-brown spots; (17) lateral sides of head and body with wide dark reddish-brown to black lateral stripes, clearly separated from lighter dorsal coloration by straight contrasting edge; (18) ventral surfaces of body, throat, and chest greyish-blue with indistinct brown confluent blotches; (19) upper eyelids with few (3-5) very small flat reddish superciliary tubercles; (20) limbs dorsally reddish-brown, ventrally brown with small bluish-white speckles. The new species is also distinct from all congeners in 12S rRNA to 16S rRNA mitochondrial DNA fragment sequences (uncorrected genetic distance P>8.9%). Advertisement call and tadpole morphology of the new species are described. Our molecular data showed Theloderma auratum sp. nov. to be a sister species of Th. palliatum from Langbian Plateau in southern Vietnam.

Geological data indicate that the interpretation for the age-calibrated phylogeny for the Kurixalus-genus frogs of South, South-east and East Asia (Lv et al., 2018) needs to be rethought

Recently, Lv et al. (2018) published an age-calibrated phylogenetic tree for the Kurixalus frogs, members of which occur across parts of South, South-east and East Asia. A clade on Taiwan, represented by Kurixalus idiootocus and the Kurixalus eiffingeri species complex, is deemed to have been resident since the middle Cenozoic; its closest congeners are in southern Indochina (not in the adjacent parts of south-east China), and the split between the two is dated at 32.8 Ma. Furthermore, a sub-population of Kurixalus eiffingeri is believed to have colonized islands in the western Ryukyus c. 13.5 Ma. There is, however, a problem with this scenario: the landmass regarded as modern-day Taiwan has existed only for 4-5 million years (it results from a young and ongoing tectonic-plate collision). Assuming the Kurixalus phylogeny and the dating of its branchings are correct, then a palaeobiogeographical scenario involving an older, alternative land surface with later transfer to Taiwan, possibly involving over-water dispersal, would reconcile the biology, but testing this may be difficult/impossible. If the ages of the nodes in the proposed tree are found to be significantly overestimated, the geology and biology might more easily be accommodated.

Coalescent-based delimitation outperforms distance-based methods for delineating less divergent species: the case of Kurixalus odontotarsus species group

Few empirical studies have compared coalescent-based methods to distance-based methods for delimitation of less divergent species. In this study, we used two coalescent-based (BFD and BPP) and two distance-based barcoding (ABGD and jMOTU) methods to delimit closely related species in the Kurixalus odontotarsus species group. Phylogenetic analyses revealed that the K. odontotarsus species group comprises 11 distinct maternal clades with strong support values. Based on the genetic and morphological evidences, we consider that species diversity in the K. odontotarsus species group was underestimated and the 11 clades represent 11 species, of which six are unnamed. The coalescent-based delimitations decisively supported the scenario of 11-species corresponding to the 11 clades. However, the distance-based ABGD only obtained 3-6 candidate species, which is not consistent with morphological evidence. These results indicate that BFD and BPP are more conservative than ABGD to false negatives (lumping). Method of fixed threshold (jMOTU) may obtain a resolution similar to that inferred by BFD and BPP, but it severely relies on subjective choice of the threshold and lacks statistical support. We consider that coalescent-based BFD and BPP approaches outperform distance-based methods for delineation of less divergent species.

The reanalysis of biogeography of the Asian tree frog, Rhacophorus (Anura: Rhacophoridae): geographic shifts and climatic change influenced the dispersal process and diversification

Rapid uplifts of the Tibetan Plateau and climate change in Asia are thought to have profoundly modulated the diversification of most of the species distributed throughout Asia. The ranoid tree frog genus Rhacophorus, the largest genus in the Rhacophoridae, is widely distributed in Asia and especially speciose in the areas south and east of the Tibetan Plateau. Here, we infer phylogenetic relationships among species and estimate divergence times, asking whether the spatiotemporal characteristics of diversification within Rhacophorus were related to rapid uplifts of the Tibetan Plateau and concomitant climate change. Phylogenetic analysis recovered distinct lineage structures in Rhacophorus, which indicated a clear distribution pattern from Southeast Asia toward East Asia and India. Molecular dating suggests that the first split within the genus date back to the Middle Oligocene (approx. 30 Ma). The Rhacophorus lineage through time (LTT) showed that there were periods of increased speciation rate: 14–12 Ma and 10–4 Ma. In addition, ancestral area reconstructions supported Southeast Asia as the ancestral area of Rhacophorus. According to the results of molecular dating, ancestral area reconstructions and LTT we think the geographic shifts, the staged rapid rises of the Tibetan Plateau with parallel climatic changes and reinforcement of the Asian monsoons (15 Ma, 8 Ma and 4–3 Ma), possibly prompted a burst of diversification in Rhacophorus.

A comprehensive phylogeny of the genus Kurixalus (Rhacophoridae, Anura) sheds light on the geographical range evolution of frilled swamp treefrogs

Currently, the genus Kurixalus comprises 14 species distributed in Southern, Southeast and East Asia. Because of their relatively low dispersal capability and intolerance of seawater, this group is ideal for the study of terrestrial range evolution, especially that portion of its range that extends into the island archipelagos of Southern Asia. We assembled a large dataset of mitochondrial and nuclear genes, and estimated phylogeny by maximum likelihood and Bayesian methods, and we explored the history of each species via divergence-time estimation based on fossil-calibrations. A variety of ancestral-area reconstruction strategies were employed to estimate past changes of the species' geographical range, and to evaluate the impact of different abiotic barriers on range evolution. We found that frilled swamp treefrogs probably originated in Taiwan or South Vietnam in the Oligocene. Alternatively, the lineage leading to Kurixalus appendiculatus strongly supports a hypothesis of terrestrial connection between the Indian and Asian continents in the Oligocene. The outcome of both our divergence-time estimates and ancestral-area reconstruction suggests that the divergence between species from Indochina and Taiwan can probably be attributed to the opening of the South China Sea, approximately 33 million years ago. We could not find evidence for dispersal between mainland China and Taiwan Island. Formation of both Mekong and Red River valleys did not have any impact on Kurixalus species diversification. However, coincidence in timing of climate change and availability of plausible dispersal routes from the Oligocene to the middle Miocene, plausibly implied that Kurixalus diversification in Asia resulted from contemporaneous, climate-induced environmental upheaval (Late Oligocene Warming at 29 Ma; Mi-1 glaciation since 24.4-21.5 Ma; Mid-Miocene Climatic Optimum at 14 Ma), which alternatively opened and closed dispersal routes.

Paleoclimate determines diversification patterns in the fossorial snake family Uropeltidae Cuvier, 1829

Understanding how and why diversification rates vary across evolutionary time is central to understanding how biodiversity is generated and maintained. Recent mathematical models that allow estimation of diversification rates across time from reconstructed phylogenies have enabled us to make inferences on how biodiversity copes with environmental change. Here, we explore patterns of temporal diversification in Uropeltidae, a diverse fossorial snake family. We generate a time-calibrated phylogenetic hypothesis for Uropeltidae and show a significant correlation between diversification rate and paleotemperature during the Cenozoic. We show that the temporal diversification pattern of this group is punctuated by one rate shift event with a decrease in diversification and turnover rate between ca. 11Ma to present, but there is no strong support for mass extinction events. The analysis indicates higher turnover during periods of drastic climatic fluctuations and reduced diversification rates associated with contraction and fragmentation of forest habitats during the late Miocene. Our study highlights the influence of environmental fluctuations on diversification rates in fossorial taxa such as uropeltids, and raises conservation concerns related to present rate of climate change.

Plio-Pleistocene phylogeography of the Southeast Asian Blue Panchax killifish, Aplocheilus panchax

The complex climatic and geological history of Southeast Asia has shaped this region’s high biodiversity. In particular, sea level fluctuations associated with repeated glacial cycles during the Pleistocene both facilitated, and limited, connectivity between populations. In this study, we used data from two mitochondrial and three anonymous nuclear markers to determine whether a fresh/brackish water killifish, Aplocheilus panchax, Hamilton, 1822, could be used to further understand how climatic oscillations and associated sea level fluctuations have shaped the distribution of biota within this region, and whether such patterns show evidence of isolation within palaeodrainage basins. Our analyses revealed three major mitochondrial clades within A. panchax. The basal divergence of A. panchax mitochondrial lineages was approximately 3.5 Ma, whilst the subsequent divergence timings of these clades occurred early Pleistocene (~2.6 Ma), proceeding through the Pleistocene. Continuous phylogeographic analysis showed a clear west-east dispersal followed by rapid radiation across Southeast Asia. Individuals from Krabi, just north of the Isthmus of Kra, were more closely related to the Indian lineages, providing further evidence for a freshwater faunal disjunction at the Isthmus of Kra biogeographic barrier. Our results suggest that Sulawesi, across the Wallace Line, was colonised relatively recently (~30 ka). Nuclear DNA is less geographically structured, although Mantel tests indicated that nuclear genetic distances were correlated with geographic proximity. Overall, these results imply that recent gene flow, as opposed to historical isolation, has been the key factor determining patterns of nuclear genetic variation in A. panchax, however, some evidence of historical isolation is retained within the mitochondrial genome. Our study further validates the existence of a major biogeographic boundary at the Kra Isthmus, and also demonstrates the use of widely distributed fresh/brackishwater species in phylogeographic studies, and their ability to disperse across major marine barriers in relatively recent time periods.

Phylogenetic relationships and divergence dates of softshell turtles (Testudines: Trionychidae) inferred from complete mitochondrial genomes

The softshell turtles (Trionychidae) are one of the most widely distributed reptile groups in the world, and fossils have been found on all continents except Antarctica. The phylogenetic relationships among members of this group have been previously studied; however, disagreements regarding its taxonomy, its phylogeography and divergence times are still poorly understood as well. Here, we present a comprehensive mitogenomic study of softshell turtles. We sequenced the complete mitochondrial genomes of 10 softshell turtles, in addition to the GenBank sequence of Dogania subplana, Lissemys punctata, Trionyx triunguis, which cover all extant genera within Trionychidae except for Cyclanorbis and Cycloderma. These data were combined with other mitogenomes of turtles for phylogenetic analyses. Divergence time calibration and ancestral reconstruction were calculated using BEAST and RASP software, respectively. Our phylogenetic analyses indicate that Trionychidae is the sister taxon of Carettochelyidae, and support the monophyly of Trionychinae and Cyclanorbinae, which is consistent with morphological data and molecular analysis. Our phylogenetic analyses have established a sister taxon relationship between the Asian Rafetus and the Asian Palea + Pelodiscus + Dogania + Nilssonia + Amyda, whereas a previous study grouped the Asian Rafetus with the American Apalone. The results of divergence time estimates and area ancestral reconstruction show that extant Trionychidae originated in Asia at around 108 million years ago (MA), and radiations mainly occurred during two warm periods, namely Late Cretaceous-Early Eocene and Oligocene. By combining the estimated divergence time and the reconstructed ancestral area of softshell turtles, we determined that the dispersal of softshell turtles out of Asia may have taken three routes. Furthermore, the times of dispersal seem to be in agreement with the time of the India-Asia collision and opening of the Bering Strait, which provide evidence for the accuracy of our estimation of divergence time. Overall, the mitogenomes of this group were used to explore the origin and dispersal route of Trionychidae and have provided new insights on the evolution of this group.