The Eurasian invasion: phylogenomic data reveal multiple Southeast Asian origins for Indian Dragon Lizards

Fossil-informed biogeographic analysis suggests Eurasian regionalization in crown Squamata during the early Jurassic

Background

Squamata (lizards, snakes, and amphisbaenians) is a Triassic lineage with an extensive and complex biogeographic history, yet no large-scale study has reconstructed the ancestral range of early squamate lineages. The fossil record indicates a broadly Pangaean distribution by the end- Cretaceous, though many lineages (e.g., Paramacellodidae, Mosasauria, Polyglyphanodontia) subsequently went extinct. Thus, the origin and occupancy of extant radiations is unclear and may have been localized within Pangaea to specific plates, with potential regionalization to distinct Laurasian and Gondwanan landmasses during the Mesozoic in some groups.

Methods

We used recent tectonic models to code extant and fossil squamate distributions occurring on nine discrete plates for 9,755 species, with Jurassic and Cretaceous fossil constraints from three extinct lineages. We modeled ancestral ranges for crown Squamata from an extant-only molecular phylogeny using a suite of biogeographic models accommodating different evolutionary processes and fossil-based node constraints from known Jurassic and Cretaceous localities. We hypothesized that the best-fit models would not support a full Pangaean distribution (i.e., including all areas) for the origin of crown Squamata, but would instead show regionalization to specific areas within the fragmenting supercontinent, likely in the Northern Hemisphere where most early squamate fossils have been found.

Results

Incorporating fossil data reconstructs a localized origin within Pangaea, with early regionalization of extant lineages to Eurasia and Laurasia, while Gondwanan regionalization did not occur until the middle Cretaceous for Alethinophidia, Scolecophidia, and some crown Gekkotan lineages. While the Mesozoic history of extant squamate biogeography can be summarized as a Eurasian origin with dispersal out of Laurasia into Gondwana, their Cenozoic history is complex with multiple events (including secondary and tertiary recolonizations) in several directions. As noted by previous authors, squamates have likely utilized over-land range expansion, land-bridge colonization, and trans-oceanic dispersal. Tropical Gondwana and Eurasia hold more ancient lineages than the Holarctic (Rhineuridae being a major exception), and some asymmetries in colonization (e.g., to North America from Eurasia during the Cenozoic through Beringia) deserve additional study. Future studies that incorporate fossil branches, rather than as node constraints, into the reconstruction can be used to explore this history further.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Hiding on jagged karst pinnacles: A new microendemic genus and species of a limestone-dwelling agamid lizard (Squamata: Agamidae: Draconinae) from Khammouan Province, Laos

We describe a unique new species and genus of agamid lizard from the karstic massifs of Khammouan Province, central Laos. Laodracon carsticola Gen. et sp. nov. is an elusive medium-sized lizard (maximum snout-vent length 101 mm) specifically adapted to life on limestone rocks and pinnacles. To assess the phylogenetic position of the new genus amongst other agamids, we generated DNA sequences from two mitochondrial gene fragments (16S rRNA and ND2) and three nuclear loci ( BDNF, RAG1 and c-mos), with a final alignment comprising 7 418 base pairs for 64 agamid species. Phylogenetic analyses unambiguously place the new genus in the mainland Asia subfamily Draconinae, where it forms a clade sister to the genus Diploderma from East Asia and the northern part of Southeast Asia. Morphologically, the new genus is distinguished from all other genera in Draconinae by possessing a notably swollen tail base with enlarged scales on its dorsal and ventral surfaces. Our work provides further evidence that limestone regions of Indochina represent unique "arks of biodiversity" and harbor numerous relict lineages. To date, Laodracon carsticola Gen. et sp. nov. is known from only two adult male specimens and its distribution seems to be restricted to a narrow limestone massif on the border of Khammouan and Bolikhamxai provinces of Laos. Additional studies are required to understand its life history, distribution, and conservation status.

Their fates intertwined: diversification patterns of the Asian gliding vertebrates may have been forged by dipterocarp trees

The repeated evolution of gliding in diverse Asian vertebrate lineages is hypothesized to have been triggered by the dominance of tall dipterocarp trees in the tropical forests of Southeast Asia. These dipterocarp forests have acted as both centres of diversification and climatic refugia for gliding vertebrates, and support most of their extant diversity. We predict similarities in the diversification patterns of dipterocarp trees and gliding vertebrates, and specifically test whether episodic diversification events such as rate shifts and/or mass extinctions were temporally congruent in these groups. We analysed diversification patterns in reconstructed timetrees of Asian dipterocarps, the most speciose gliding vertebrates from different classes (Draco lizards, gliding frogs and Pteromyini squirrels) and compared them with similar-sized clades of non-gliding relatives (Diploderma lizards, Philautus frogs and Callosciurinae squirrels) from Southeast Asia. We found significant declines in net-diversification rates of dipterocarps and the gliding vertebrates during the Pliocene-Pleistocene, but not in the non-gliding groups. We conclude that the homogeneity and temporal coincidence of these rate declines point to a viable ecological correlation between dipterocarps and the gliding vertebrates. Further, we suggest that while the diversification decay in dipterocarps was precipitated by post-Miocene aridification of Asia, the crises in the gliding vertebrates were induced by both events concomitantly.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Taxonomic status of the Lesser Woolly Horseshoe bats (Chiroptera, Rhinolophidae, Rhinolophus beddomei) in peninsular India and Sri Lanka

Insularity provides ample opportunities for species diversification. Sri Lanka is home to a large diversity of species, many of which are endemic but morphologically similar to species found in southern India, due to recent speciation events, suggesting a complex evolutionary history. However, in some taxa although morphological diversity has been noted, the genetic level variations are minimal. Among the wide-ranging horseshoe bats such a phenomenon is noted. In this study, we used bioacoustics, morphometric and molecular data to evaluate the relationships between the taxa of lesser woolly horseshoe bats in the India and Sri Lanka. Our study reveals that the two taxa-Rhinolophus beddomei Andersen, 1905 and here we have validated the existing subspecies from peninsular India and R. sobrinus Andersen, 1918 from Sri Lanka are genetically very close to R. perniger Hodgson, 1843. Currently the taxa-beddomei and sobrinus are recognized as subspecies of Rhinolophus beddomei Andersen, 1905. We provide a detailed description of the taxa beddomei and sobrinus as the original descriptions are limited in nature.

Phylogenomics and historical biogeography of West Indian Rock Iguanas (genus Cyclura)

The genus Cyclura includes nine extant species and six subspecies of West Indian Rock Iguanas and is one of the most imperiled genera of squamate reptiles globally. An understanding of species diversity, evolutionary relationships, diversification, and historical biogeography in this group is crucial for implementing sound long-term conservation strategies. We collected DNA samples from 1 to 10 individuals per taxon from all Cyclura taxa (n = 70 ingroup individuals), focusing where possible on incorporating individuals from different populations of each species. We also collected 1-2 individuals from each of seven outgroup species of iguanas (Iguana delicatissima; five Ctenosaura species) and Anolis sagrei (n = 12 total outgroup individuals). We used targeted genomic sequence capture to isolate and to sequence 1,872 loci comprising of 687,308 base pairs (bp) from each of the 82 individuals from across the nuclear genome. We extracted mitochondrial reads and assembled and annotated mitogenomes for all Cyclura taxa plus outgroup species. We present well-supported phylogenomic gene tree/species tree analyses for all extant species of Cyclura using ASTRAL-III, SVDQuartets, and StarBEAST2 methods, and discuss the taxonomic, biogeographic, and conservation implications of these data. We find a most recent common ancestor of the genus 9.91 million years ago. The earliest divergence within Cyclura separates C. pinguis from a clade comprising all other Cyclura. Within the latter group, a clade comprising C. carinata from the southern Lucayan Islands and C. ricordii from Hispaniola is the sister taxon to a clade comprising the other Cyclura. Among the other Cyclura, the species C. cornuta and C. stejnegeri (from Hispaniola and Isla Mona) form the sister taxon to a clade of species from Jamaica (C. collei), Cuba and Cayman Islands (C. nubila and C. lewisi), and the eastern (C. rileyi) and western (C. cychlura) Lucayan Islands. Cyclura cychlura and C. rileyi form a clade whose sister taxa are C. nubila and C. lewisi. Cyclura collei is the sister taxon to these four species combined.

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.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

Mesozoic origin and 'out-of-India' radiation of ricefishes (Adrianichthyidae)

The Indian subcontinent has an origin geologically different from Eurasia, but many terrestrial animal and plant species on it have congeneric or sister species in other parts of Asia, especially in the Southeast. This faunal and floral similarity between India and Southeast Asia is explained by either of the two biogeographic scenarios, 'into-India' or 'out-of-India'. Phylogenies based on complete mitochondrial genomes and five nuclear genes were undertaken for ricefishes (Adrianichthyidae) to examine which of these two biogeographic scenarios fits better. We found that Oryzias setnai, the only adrianichthyid distributed in and endemic to the Western Ghats, a mountain range running parallel to the western coast of the Indian subcontinent, is sister to all other adrianichthyids from eastern India and Southeast-East Asia. Divergence time estimates and ancestral area reconstructions reveal that this western Indian species diverged in the late Mesozoic during the northward drift of the Indian subcontinent. These findings indicate that adrianichthyids dispersed eastward 'out-of-India' after the collision of the Indian subcontinent with Eurasia, and subsequently diversified in Southeast-East Asia. A review of geographic distributions of 'out-of-India' taxa reveals that they may have largely fuelled or modified the biodiversity of Eurasia.

On the systematics and the phylogenetic position of the poorly known, montane dragon-lizard species Pseudocalotes austeniana (Annandale, 1908) (Squamata, Agamidae, Draconinae)

The montane agamid species Pseudocalotes austeniana has had a complicated taxonomic history, as the species was initially described as a member of the genus Salea Gray, 1845. Later, the species was placed in a monotypic genus Mictopholis Smith, 1935, which was erected only to include this species; however, the species was later on transferred to the genus Pseudocalotes Fitzinger, 1843, owing to the morphological similarities, and lack of strong characters to diagnose the genus Mictopholis. Nonetheless, its precise phylogenetic and systematic position has remained unresolved due to the lack of molecular sequence data. During a herpetological expedition to Arunachal Pradesh, specimens of P. austeniana were collected from the hills near the type locality. The mitochondrial 16S rRNA, ND2 and ND4, and the nuclear RAG1 regions were subjected to molecular phylogenetics. Maximum Likelihood and Bayesian Inference gene trees revealed that P. austeniana is a member of the subfamily Draconinae. The analyses showed that the genus Pseudocalotes is polyphyletic, and P. austeniana was embedded within the genus Japalura Gray, 1853 sensu stricto. We here, thus, propose to transfer the species P. austeniana to the genus Japlaura, as Japalura austenianacomb. nov. Biogeographic and evolutionary significance of the findings are discussed.

Phylogenomic Analysis Reveals Dispersal-Driven Speciation and Divergence with Gene Flow in Lesser Sunda Flying Lizards (Genus Draco)

The Lesser Sunda Archipelago offers exceptional potential as a model system for studying the dynamics of dispersal-driven diversification. The geographic proximity of the islands suggests the possibility for successful dispersal, but this is countered by the permanence of the marine barriers and extreme intervening currents that are expected to hinder gene flow. Phylogenetic and species delimitation analyses of flying lizards (genus Draco) using single mitochondrial genes, complete mitochondrial genomes, and exome-capture data sets identified 9-11 deeply divergent lineages including single-island endemics, lineages that span multiple islands, and parapatrically-distributed non-sister lineages on the larger islands. Population clustering and PCA confirmed these genetic boundaries with isolation-by-distance playing a role in some islands or island sets. While gdi estimates place most candidate species comparisons in the ambiguous zone, migration estimates suggest 9 or 10 species exist with nuclear introgression detected across some intra-island contact zones. Initial entry of Draco into the archipelago occurred at 5.5-7.5 Ma, with most inter-island colonization events having occurred between 1-3 Ma. Biogeographical model testing favors scenarios integrating geographic distance and historical island connectivity, including an initial stepping-stone dispersal process from the Greater Sunda Shelf through the Sunda Arc as far eastward as Lembata Island. However, rather than reaching the adjacent island of Pantar by dispersing over the 15-km wide Alor Strait, Draco ultimately reached Pantar (and much of the rest of the archipelago) by way of a circuitous route involving at least five over-water dispersal events. These findings suggest that historical geological and oceanographic conditions heavily influenced dispersal pathways and gene flow, which in turn drove species formation and shaped species boundaries.

Phylogenetics of mud snakes (Squamata: Serpentes: Homalopsidae): A paradox of both undescribed diversity and taxonomic inflation

Mud snakes (Serpentes: Homalopsidae) are a family of 55 described, mainly aquatic, species primarily distributed throughout mainland Southeast Asia and the Indo-Australian Archipelago. Although they have been the focus of prior research, the basic relationships amongst genera and species remain poorly known. We used a combined mitochondrial and nuclear gene dataset to infer their phylogenetic relationships, using the highest levels of taxon and geographic sampling for any homalopsid phylogeny to date (62% generic and 62% species coverage; 140 individuals). Our results recover two reciprocally monophyletic groups: the fangless Brachyorrhos and its sister clade comprised of all rear-fanged homalopsids. Most genera and interspecific relationships were monophyletic and strongly supported, but intergeneric relationships and intraspecific population structure lack support. We find evidence of both undescribed diversity as well as cases of taxonomic inflation within several species. Tree-based species delimitation approaches (mPTP) support potential new candidate species as distinct from their conspecifics and also suggest that many named taxa may not be distinct species. Divergence date estimation and lineage-through-time analyses indicate lower levels of speciation in the Eocene, with a subsequent burst in diversification in the Miocene. Homalopsids may have diversified most rapidly during the Pliocene and Pleistocene, possibly in relation to tectonic shifts and sea-level fluctuations that took place in Sundaland and the Sahul Shelf. Our analyses provide new insights on homalopsid taxonomy, a baseline phylogeny for the family, and further biogeographic implications demonstrating how dynamic tectonics and Quaternary sea level changes may have shaped a widespread, diverse family of snakes.

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.].

India's biogeographic history through the eyes of blindsnakes- filling the gaps in the global typhlopoid phylogeny

The Indian subcontinent's unique geological history is reflected in the diverse assemblage of its biota. The blindsnake superfamily Typhlopoidea, with its unique mix of ancient as well as younger lineages in Asia, provides an opportunity to understand the various biotic exchange scenarios proposed for the Indian landmass. In this study, we aim to understand the biogeographic origins of the four genera of typhlopoids found in India and to decipher their times and modes of arrival in the subcontinent. Five nuclear markers were sequenced for 12 samples collected from across India, encompassing all four genera under study. Published sequences of typhlopoid genera were compiled and combined with Indian sequences to generate a global dataset. Phylogenetic relationships were reconstructed using maximum likelihood and Bayesian inference methods. Divergence times were estimated using BEAST 1.8.2. Ancestral geographical ranges were estimated using DEC + J, implemented in BioGeoBEARS. Divergence time estimates suggest that Gerrhopilus is an ancient lineage, and the lineage leading to it was present on the Indian landmass since the last 100 million years. The other three genera are more recent dispersals into India, possibly trans-oceanic. Biogeographic reconstructions suggest an East Gondwanan origin for Typhlopoidea, an African origin for Grypotyphlops and an Asian origin for Indotyphlops and Argyrophis. It appears that India harbours a combination of ancient and more recently dispersed lineages of typhlopoids. The genus Gerrhopilus is of Gondwanan origin that likely dispersed out of India into Southeast Asia. The other genera are intrusive elements that dispersed into India from Africa (Grypotyphlops) and Asia (Indotyphlops and possibly Argyrophis) post break-up of Gondwana. Thus, our study provides further evidence on the ability of blindsnakes to undergo long distance trans-oceanic dispersal. Results also suggest an Asian origin for typhlopoids from Australasia, Philippines and Wallacea.

Phylogeny and biogeography of Sumatra´s cloud forest lizards of the genus Dendragama and status of Acanthosaura schneideri

Lizards of the genus Dendragama are endemic to the highland cloud forests of Sumatra’s Barisan Mountain Range in western Indonesia, and recent studies have uncovered widespread diversity within the genus. Here, a suite of morphological characters and mitochondrial DNA are used to compare three geographically isolated populations of D.boulengeri from (1) Mount Kerinci in Jambi province, (2) Mount Marapi of west Sumatra, and (3) the Karo Highlands of north Sumatra. Additional phylogeographic analyses with two recently described sister species, D.australis and D.dioidema were conducted. Five genetically distinct clades of Dendragama, all distributed allopatrically of one another were identified and some are suspected to inhabit small distributions. Morphological and genetic data confirm the Karo Highlands population D.schneideri (previously Acanthosauraschneideri Ahl, 1926) should be revalidated from the synonymy of D.boulengeri. Dendragamaschneideri is endemic to montane forests of the Karo Highlands surrounding Lake Toba in Sumatra Utara province. Pairwise genetic distances of 6–11% separate D.schneideri from congeners. Two distinct clades of D.boulengeri from Mount Kerinci and Mount Marapi were identified, which are 5.0% genetically distant from one another. Using morphological characters, we provide the first key for distinguishing between species of Dendragama. Based on biogeographic patterns and levels of genetic variation it is suspected that at least 18 other isolated cloud forest locations may hold new species or divergent populations of Dendragama but lack survey work. Collectively, these comparisons among populations of montane lizards further elucidate the complex biogeographic history of Sumatra’s montane forest species and the first phylogeny of the genus Dendragama.

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.

Phylogeography of montane dragons could shed light on the history of forests and diversification processes on Sumatra

Biogeographical evidence, both, for and against the "regional endemism paradigm" hypothesis has been uncovered across the Greater Sunda Region (Sundaland) of Southeast Asia. Additionally, there are competing hypotheses regarding how Pleistocene forests may have impacted biological patterns and processes in Sumatra. Using montane agamid lizards from Sumatra, we derived and analyzed a phylogenetic dataset, genetic divergence estimates, and contemporary distributional patterns among species. We tested whether (1) Sumatra's highland Draconinae diversification fits the regional endemism paradigm hypothesis and (2) Draconinae phylogeography provides biological evidence for Pleistocene forest extent at various points in history. Our results suggest in situ diversification was the main driver behind montane Draconinae lizard diversification in Sumatra, rejecting the "regional endemism paradigm". Contemporary distribution of endemic species and their genetic relationships may potentially provide biologicalevidence for determining more precise elevational lower limits of montane forests during the Pleistocene epoch. Our data suggests montane forests did not retreat more than 700-750 m during glacial maxima because lower retreating forests would have become interconnected, allowing for widespread dispersal, exchange of gene flow and sympatric distributions contemporarily. To the contrary, our divergence estimates show that cloud forest dragons have been isolated for millions of years, suggesting there may have been a continuous disconnect between some areas, predating the Pleistocene. There may also be other ecological and evolutionary factors that impacted Draconinae distributions, such as competition, making this an excellent system for testing questions regarding montane biogeography. Additionally, we provide the first phylogeny for a wide range of Sundaland agamid species and identify some biogeographic pressures that may have triggered montane Draconinae diversification in Sumatra.

Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements

Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

Ultraconserved element (UCE) probe set design: Base genome and initial design parameters critical for optimization

Targeted capture and enrichment approaches have proven effective for phylogenetic study. Ultraconserved elements (UCEs) in particular have exhibited great utility for phylogenomic analyses, with the software package phyluce being among the most utilized pipelines for UCE phylogenomics, including probe design. Despite the success of UCEs, it is becoming increasing apparent that diverse lineages require probe sets tailored to focal taxa in order to improve locus recovery. However, factors affecting probe design and methods for optimizing probe sets to focal taxa remain underexplored. Here, we use newly available beetle (Coleoptera) genomic resources to investigate factors affecting UCE probe set design using phyluce. In particular, we explore the effects of stringency during initial design steps, as well as base genome choice on resulting probe sets and locus recovery. We found that both base genome choice and initial bait design stringency parameters greatly alter the number of resultant probes included in final probe sets and strongly affect the number of loci detected and recovered during in silico testing of these probe sets. In addition, we identify attributes of base genomes that correlated with high performance in probe design. Ultimately, we provide a recommended workflow for using phyluce to design an optimized UCE probe set that will work across a targeted lineage, and use our findings to develop a new, open-source UCE probe set for beetles of the suborder Adephaga.

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.

Aridification driven diversification of fan-throated lizards from the Indian subcontinent

The establishment of monsoon climate and the consequent aridification has been one of the most important climate change episodes in the Indian subcontinent. However, little is known about how these events might have shaped the diversification patterns among the widely distributed taxa. Fan-throated lizards (FTL) (Genus: Sitana, Sarada) are widespread, diurnal and restricted to the semi-arid zones of the Indian subcontinent. We sampled FTL in 107 localities across its range. We used molecular species delimitation method and delineated 15 species including six putative species. Thirteen of them were distinguishable based on morphology but two sister species were indistinguishable and have minor overlaps in distribution. Five fossils were used to calibrate and date the phylogeny. Diversification of fan-throated lizards lineage started ~18 mya and higher lineage diversification was observed after 11 my. The initial diversification corresponds to the time when monsoon climate was established and the latter was a period of intensification of monsoon and initiation of aridification. Thirteen out of the fifteen FTL species delimited are from Peninsular India; this is probably due to the landscape heterogeneity in this region. The species poor sister genus Otocryptis is paraphyletic and probably represents relict lineages which are now confined to forested areas. Thus, the seasonality led changes in habitat, from forests to open habitats appear to have driven diversification of fan-throated lizards.

Phylogenomic analysis of the Chilean clade of Liolaemus lizards (Squamata: Liolaemidae) based on sequence capture data

The genus Liolaemus is one of the most ecologically diverse and species-rich genera of lizards worldwide. It currently includes more than 250 recognized species, which have been subject to many ecological and evolutionary studies. Nevertheless, Liolaemus lizards have a complex taxonomic history, mainly due to the incongruence between morphological and genetic data, incomplete taxon sampling, incomplete lineage sorting and hybridization. In addition, as many species have restricted and remote distributions, this has hampered their examination and inclusion in molecular systematic studies. The aims of this study are to infer a robust phylogeny for a subsample of lizards representing the Chilean clade (subgenus Liolaemus sensu stricto), and to test the monophyly of several of the major species groups. We use a phylogenomic approach, targeting 541 ultra-conserved elements (UCEs) and 44 protein-coding genes for 16 taxa. We conduct a comparison of phylogenetic analyses using maximum-likelihood and several species tree inference methods. The UCEs provide stronger support for phylogenetic relationships compared to the protein-coding genes; however, the UCEs outnumber the protein-coding genes by 10-fold. On average, the protein-coding genes contain over twice the number of informative sites. Based on our phylogenomic analyses, all the groups sampled are polyphyletic. Liolaemus tenuis tenuis is difficult to place in the phylogeny, because only a few loci (nine) were recovered for this species. Topologies or support values did not change dramatically upon exclusion of L. t. tenuis from analyses, suggesting that missing data did not had a significant impact on phylogenetic inference in this data set. The phylogenomic analyses provide strong support for sister group relationships between L. fuscus, L. monticola, L. nigroviridis and L. nitidus, and L. platei and L. velosoi. Despite our limited taxon sampling, we have provided a reliable starting hypothesis for the relationships among many major groups of the Chilean clade of Liolaemus that will help future work aimed at resolving the Liolaemus phylogeny.

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.