Phylogeography of the Mekong mud snake (Enhydris subtaeniata): the biogeographic importance of dynamic river drainages and fluctuating sea levels for semiaquatic taxa in Indochina

Integrative methods reveal multiple drivers of diversification in rice paddy snakes

Divergence dating analyses in systematics provide a framework to develop and test biogeographic hypotheses regarding speciation. However, as molecular datasets grow from multilocus to genomic, sample sizes decrease due to computational burdens, and the testing of fine-scale biogeographic hypotheses becomes difficult. In this study, we use coalescent demographic models to investigate the diversification of poorly known rice paddy snakes from Southeast Asia (Homalopsidae: Hypsiscopus), which have conflicting dates of origin based on previous studies. We use coalescent modeling to test the hypothesis that Hypsiscopus diversified 2.5 mya during the Khorat Plateau uplift in Thailand. Additionally, we use ecological niche analyses to identify potential differences in the niche space of the two most widely distributed species in the past and present. Our results suggest Hypsiscopus diversified ~ 2.4 mya, supporting that the Khorat Plateau may have initiated the diversification of rice paddy snakes. We also find significant niche differentiation and shifts between species of Hypsiscopus, indicating that environmental differences may have sustained differentiation of this genus after the Khorat Plateau uplift. Our study expands on the diversification history of snakes in Southeast Asia, and highlights how results from smaller multilocus datasets can be useful in developing and testing biogeographic hypotheses alongside genomic datasets.

East‒West genetic differentiation across the Indo-Burma hotspot: evidence from two closely related dioecious figs

BACKGROUND

Understanding biodiversity patterns and their underlying mechanisms is of interest to ecologists, biogeographers and conservationists and is critically important for conservation efforts. The Indo-Burma hotspot features high species diversity and endemism, yet it also faces significant threats and biodiversity losses; however, few studies have explored the genetic structure and underlying mechanisms of Indo-Burmese species. Here, we conducted a comparative phylogeographic analysis of two closely related dioecious Ficus species, F. hispida and F. heterostyla, based on wide and intensive population sampling across Indo-Burma ranges, using chloroplast (psbA-trnH, trnS-trnG) and nuclear microsatellite (nSSR) markers, as well as ecological niche modeling.

RESULTS

The results indicated large numbers of population-specific cpDNA haplotypes and nSSR alleles in the two species. F. hispida showed slightly higher chloroplast diversity but lower nuclear diversity than F. heterostyla. Low-altitude mountainous areas of northern Indo-Burma were revealed to have high genetic diversity and high habitat suitability, suggesting potential climate refugia and conservation priority areas. Strong phylogeographic structure and a marked east‒west differentiation pattern were observed in both species, due to the interactions between biotic and abiotic factors. Interspecific dissimilarities at fine-scale genetic structure and asynchronized historical dynamics of east‒west differentiation between species were also detected, which were attributed to different species-specific traits.

CONCLUSIONS

We confirm hypothesized predictions that interactions between biotic and abiotic factors largely determine the patterns of genetic diversity and phylogeographic structure of Indo-Burmese plants. The east‒west genetic differentiation pattern observed in two targeted figs can be generalized to some other Indo-Burmese plants. The results and findings of this work will contribute to the conservation of Indo-Burmese biodiversity and facilitate targeted conservation efforts for different species.

Phylogeographic inference of Sumatran ranids bearing gastromyzophorous tadpoles with regard to the Pleistocene drainage systems of Sundaland

Rivers are known to act as biogeographic barriers in several strictly terrestrial taxa, while possibly serving as conduits of dispersal for freshwater-tolerant or -dependent species. However, the influence of river systems on genetic diversity depends on taxa-specific life history traits as well as other geographic factors. In amphibians, several studies have demonstrated that river systems have only minor influence on their divergence. Here, we assess the role of the paleodrainage systems of the Sunda region (with a focus on the island of Sumatra) in shaping the evolutionary history of two genera of frogs (Sumaterana and Wijayarana) whose tadpoles are highly dependent on cascading stream habitats. Our phylogenetic results show no clear association between the genetic diversification patterns of both anurans genera and the existence of paleodrainage systems. Time-calibrated phylogenies and biogeographical models suggest that these frogs colonized Sumatra and diversified on the island before the occurrence of the Pleistocene drainage systems. Both genera demonstrate phylogenetic structuring along a north–south geographic axis, the temporal dynamics of which coincide with the geological chronology of proto Sumatran and -Javan volcanic islands. Our results also highlight the chronic underestimation of Sumatran biodiversity and call for more intense sampling efforts on the island.

Homogenized Phylogeographic Structure across the Indo-Burma Ranges of a Large Monoecious Fig, Ficus altissima Blume

As well as bountiful natural resources, the Indo-Burma biodiversity hotspot features high rates of habitat destruction and fragmentation due to increasing human activity; however, most of the Indo-Burma species are poorly studied. The exploration of plants closely associated with human activity will further assist us to understand our influence in the context of the ongoing extinction events in the Anthropocene. This study, based on widely and intensively sampled F. altissima across Indo-Burma and the adjacent south China ranges, using both the chloroplast psbA-trnH spacer and sixteen newly developed nuclear microsatellite markers (nSSRs), aims to explore its spatial genetic structure. The results indicated low chloroplast haplotype diversity and a moderate level of nuclear genetic diversity. Although limited seed flow was revealed by psbA-trnH, no discernible phylogeographic structure was shown due to the low resolution of cpDNA markers and dominance of an ancestral haplotype. From the nSSRs data set, phylogeographic structure was homogenized, most likely due to extensive pollen flow mediated by pollinating fig wasps. Additionally, human cultivation and human-mediated transplanting further confounded the analyses of population structure. No geographic barriers are evident across the large study range, with F. altissima constituting a single population, and extensive human cultivation is likely to have had beneficial consequences for protecting the genetic diversity of F. altissima.

Isolation by geographical distance after release from Pleistocene refugia explains genetic and phenotypic variation in Xylotrupes siamensis (Coleoptera: Scarabaeidae)

Consistent and objective species delimitation is crucial to biodiversity studies, but challenges remain when conflicting taxonomic decisions have been made because different data sets and analytical methods were used to delineate species. In the rhinoceros beetle, Xylotrupes siamensis, the use of different morphological characters has resulted in taxonomic disagreement between studies. We used three molecular loci (mitochondrial CO1 and nuclear ITS2 and H3) to investigate the genetic divergence between populations exhibiting different male horn phenotypes. We also applied an approximate Bayesian computation approach to test alternative historical hypotheses that might explain the present genetic diversity among geographical populations. Furthermore, we used species distribution models to estimate the temporal variation in the geographical distribution of suitable habitats. The results show that the two phenotypic taxa within X. siamensis are not genetically structured and that their genetic structure can be explained using isolation by geographical distance. The emergence of the two phenotypic taxa might have been associated with historical isolation in separate refugia. However, spatial expansion and genetic interchange between populations might have gradually eroded the spatial genetic structure. We demonstrate that understanding the historical processes responsible for phenotypic divergence and genetic diversity among current populations could help with making evolutionarily coherent taxonomic decisions.

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.

Rivers of Indochina as potential drivers of lineage diversification in the spotted flying lizard (Draco maculatus) species complex

Southeast Asia hosts a rich concentration of biodiversity within multiple biodiversity hotspots. Indochina, a region with remarkably high levels of in situ diversification, possesses five major rivers (Ayeyarwady, Chiang Mai, Mekong, Red, and Salween), several of which coincide with phylogenetic breaks of terrestrial taxa. Draco maculatus possesses a range that stretches across Indochina, which widespread geographic distribution along with potential discrete variation within subspecies alludes to the possibility of this taxon constituting multiple divergent lineages. Using sequence data from three mitochondrial (12S, 16S, and ND2) and three nuclear (BDNF, CMOS, and PNN) genes, we provide the first estimated phylogeny of this hypothesized species complex and examine its phylogeographic architecture with maximum likelihood and Bayes factor delimitation (BFD) approaches. Our results support multiple divergent lineages with phylogenetic breaks coincident with rivers, indicating that river barriers may be contributing to the elevated levels of in situ diversification of Indochina.

Genetic diversity and population structure of blue-crested lizard, Calotes mystaceus Duméril & Bibron, 1837 (Squamata: Agamidae) in Thailand

The blue-crested lizard, Calotes mystaceus Duméril & Bibron, 1837, is listed as a protected wild animal in Thailand. Its population is likely to be dramatically reduced due to massive hunting in several areas in this country. Basic information on its population genetics is therefore needed to facilitate its conservation. Thus, in this study we investigated the mitochondrial cytochrome c oxidase subunit 1 (CO1) sequence variation of 238 individualC.mystaceus from 42 different geographical localities of the north, west, central, east and northeast regions of Thailand. High genetic diversity and genetic differentiation at the intrapopulation and interpopulation levels was observed.We detected 63 unique haplotypes and 12 common/shared haplotypes. The phylogenetic analysis reveals two major lineages, I and II. These two lineages are separated by mountain ranges, which play an important role as natural barriers blocking gene flow. Our finding reveal at least two cryptic lineages represented in C. mystaceus populations in Thailand. However, a comprehensive investigation of the morphology, biology, ecology and genetic diversity of C. mystaceus in other regions within its area of distribution is needed.

Integrative Taxonomy of Southeast Asian Snail-Eating Turtles (Geoemydidae: Malayemys) Reveals a New Species and Mitochondrial Introgression

Based on an integrative taxonomic approach, we examine the differentiation of Southeast Asian snail-eating turtles using information from 1863 bp of mitochondrial DNA, 12 microsatellite loci, morphology and a correlative species distribution model. Our analyses reveal three genetically distinct groups with limited mitochondrial introgression in one group. All three groups exhibit distinct nuclear gene pools and distinct morphology. Two of these groups correspond to the previously recognized species Malayemys macrocephala (Chao Phraya Basin) and M. subtrijuga (Lower Mekong Basin). The third and genetically most divergent group from the Khorat Basin represents a previously unrecognized species, which is described herein. Although Malayemys are extensively traded and used for religious release, only few studied turtles appear to be translocated by humans. Historic fluctuations in potential distributions were assessed using species distribution models (SDMs). The Last Glacial Maximum (LGM) projection of the predictive SDMs suggests two distinct glacial distribution ranges, implying that the divergence of M. macrocephala and M. subtrijuga occurred in allopatry and was triggered by Pleistocene climate fluctuations. Only the projection derived from the global circulation model MIROC reveals a distinct third glacial distribution range for the newly discovered Malayemys species.

The Centipede Genus Scolopendra in Mainland Southeast Asia: Molecular Phylogenetics, Geometric Morphometrics and External Morphology as Tools for Species Delimitation

Seven Scolopendra species from the Southeast Asian mainland delimited based on standard external morphological characters represent monophyletic groups in phylogenetic trees inferred from concatenated sequences of three gene fragments (cytochrome c oxidase subunit 1, 16S rRNA and 28S rRNA) using Maximum likelihood and Bayesian inference. Geometric morphometric description of shape variation in the cephalic plate, forcipular coxosternite, and tergite of the ultimate leg-bearing segment provides additional criteria for distinguishing species. Colouration patterns in some Scolopendra species show a high degree of fit to phylogenetic trees at the population level. The most densely sampled species, Scolopendra dehaani Brandt, 1840, has three subclades with allopatric distributions in mainland SE Asia. The molecular phylogeny of S. pinguis Pocock, 1891, indicated ontogenetic colour variation among its populations. The taxonomic validation of S. dawydoffi Kronmüller, 2012, S. japonica Koch, 1878, and S. dehaani Brandt, 1840, each a former subspecies of S. subspinipes Leach, 1814 sensu Lewis, 2010, as full species was supported by molecular information and additional morphological data. Species delimitation in these taxonomically challenging animals is facilitated by an integrative approach that draws on both morphology and molecular phylogeny.