Genealogy and Demographic History of a Widespread Amphibian throughout Indochina

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.

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.

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.

Spatial dynamics of Chinese Muntjac related to past and future climate fluctuations

Climate fluctuations in the past and in the future are likely to result in population expansions, shifts, or the contraction of the ecological niche of many species, and potentially leading to the changes in their geographical distributions. Prediction of suitable habitats has been developed as a useful tool for the assessment of habitat suitability and resource conservation to protect wildlife. Here, we model the ancestral demographic history of the extant modern Chinese Muntjac Muntiacus reevesi populations using approximate Bayesian computation (ABC) and used the maximum entropy model to simulate the past and predict the future spatial dynamics of the species under climate oscillations. Our results indicated that the suitable habitats for the M. reevesi shifted to the Southeast and contracted during the Last Glacial Maximum, whereas they covered a broader and more northern position in the Middle Holocene. The ABC analyses revealed that the modern M. reevesi populations diverged in the Middle Holocene coinciding with the significant contraction of the highly suitable habitat areas. Furthermore, our predictions suggest that the potentially suitable environment distribution for the species will expand under all future climate scenarios. These results indicated that the M. reevesi diverged in the recent time after the glacial period and simultaneously as its habitat’s expanded in the Middle Holocene. Furthermore, the past and future climate fluctuation triggered the change of Chinese muntjac spatial distribution, which has great influence on the Chinese muntjac’s population demographic history.

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.

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.

Molecular phylogeny, species delimitation and biogeographic history of the Stegana (Steganina) shirozui species group (Diptera: Drosophilidae) from East Asia

The Stegana (Steganina) shirozui species group is mainly distributed in East Asia. In the present study, the molecular phylogeny of the S. shirozui group was investigated based on mitochondrial (COI and ND2) and nuclear (28S rRNA) markers. The resulting trees support the S. shirozui group as monophyletic and indicate that in this group, species associated with closer affinities show higher structural homogeneity in male genitalia. Molecular species delimitation assess most species limits and recognize four new species in the S. shirozui group from south-west China: S. alianya sp. nov., S. diodonta sp. nov., S. zebromyia sp. nov. and S. zopheria sp. nov. One new synonym was also recognized. Additionally, three typical male genital characters of the S. shirozui group were placed on the molecular phylogenetic framework. The outcome of both divergence-time estimation and ancestral area reconstruction suggests that the S. shirozui group likely originated in south-west China in the Middle Miocene.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Phylogenetic relationships and genetic diversity of the Polypedates leucomystax complex in Thailand

Taxonomic uncertainty of the Asian tree frog Polypedates leucomystax complex presents the challenging task of inferring its biogeographical history. Here, we describe its dispersion and the genetic relationships among different populations in Thailand, where we connect the population of the P. leucomystax complex of the Sunda Islands to the Indochina (mainland) population based on analyses of 266 sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Our maternal genealogy implies that there are four well-supported lineages in Thailand, consisting of Northern A (clade A: Polypedates sp.), Nan (clade B: P. cf. impresus), Southern (clade C: P. cf. leucomystax) and Northern D (clade D: P. cf. megacephalus), with Bayesian posterior probability >0.9. Phylogeny and haplotype networks indicate that clades A, B and D are sympatric. In contrast, clade C (P. cf. leucomystax) and clade D (P. cf. megacephalus) are genetically divergent due to the geographical barrier of the Isthmus of Kra, resulting in an allopatric distribution. Climatic conditions, in particular differences in rainfall on each side of the Isthmus of Kra, may play an important role in limiting the immigration of both clades. For the within-populations of either clades C or D, there was no significant correlation between geographic and genetic distance by the isolation-by-distance test, indicating intraspecific-dispersal of each clade. Population expansion occurred in clade C, whereas clade D showed a constant population. Taken together, the P. leucomystax complex in South East Asia may have diversified under climatic pressure, leading to allopatric and/or sympatric speciation.

Population genetic structure of the endemic rosewoods Dalbergia cochinchinensis and D. oliveri at a regional scale reflects the Indochinese landscape and life-history traits

Indochina is a biodiversity hot spot and harbors a high number of endemic species, most of which are poorly studied. This study explores the genetic structure and reproductive system of the threatened endemic timber species Dalbergia cochinchinensis and Dalbergia oliveri using microsatellite data from populations across Indochina and relates it to landscape characteristics and life‐history traits. We found that the major water bodies in the region, Mekong and Tonle Sap, represented barriers to gene flow and that higher levels of genetic diversity were found in populations in the center of the distribution area, particularly in Cambodia. We suggest that this pattern is ancient, reflecting the demographic history of the species and possible location of refugia during earlier time periods with limited forest cover, which was supported by signs of old genetic bottlenecks. The D. oliveri populations had generally high levels of genetic diversity (mean H_e = 0.73), but also strong genetic differentiation among populations (global G_ST = 0.13), while D. cochinchinensis had a moderate level of genetic diversity (mean H_e = 0.55), and an even stronger level of differentiation (global G_ST = 0.25). These differences in genetic structure can be accounted for by a higher level of gene flow in D. oliveri due to a higher dispersal capacity, but also by the broader distribution area for D. oliveri, and the pioneer characteristics of D. cochinchinensis. This study represents the first detailed analysis of landscape genetics for tree species in Indochina, and the found patterns might be common for other species with similar ecology.

Phylogeographic and phylogenetic analysis for Tripterygium species delimitation

Tripterygium wilfordii (Celastraceae) is a traditional Chinese medicine; and the dried root and rhizome constitute the main officinal parts. Tripterygium wilfordii has been identified as a potential candidate for the treatment of systemic lupus erythematosus, rheumatoid arthritis, nephritis, asthma, leprosy, and cancer. The phylogenetic relationships within the Tripterygium genus are ambiguous; thus, our aim is to clarify the relationships within this genus using phylogeographic and phylogenetic analyses. Here, we first sequenced three plastid DNA regions (i.e., psbA‐trnH, rpl32‐trnL, and trnL‐trnF) and found that Tripterygium hypoglaucum and T. wilfordii were clustered together based on the strength of the topology in the phylogenetic analysis: T. hypoglaucum is polyphyletic, and T. wilfordii is paraphyletic. A spatial analysis of molecular variance showed that the best group value is 4, and the groups were almost consistent with the topology of in the phylogenetic analysis. The Mantel analyses of Tripterygium using IBD web showed statistically significant relationships between genetic and geographical distance distributions (r = .3479, p < .0001). The molecular dating using Fossil calibration indicated that the divergence in Tripterygium was approximately 8.13 Ma. Furthermore, we also analyzed four DNA regions (i.e., ITS2, psbA‐trnH, matK, and rbcL) that were obtained from the NCBI nucleotide database; these results showed that T. wilfordii and T. hypoglaucum clustered together, while Tripterygium regelii represented a separate cluster. Tripterygium hypoglaucum and T. wilfordii were never distinct lineages, and the species circumscriptions are artificial. We propose that T. wilfordii and T. hypoglaucum are conspecific, while T. regelii likely constitutes a separate species.

Recent range expansion of an intermediate host for animal schistosome parasites in the Indo-Australian Archipelago: phylogeography of the freshwater gastropod Indoplanorbis exustus in South and Southeast Asia

Background

The planorbid snail Indoplanorbis exustus is the sole intermediate host for the Schistosoma indicum species group, trematode parasites responsible for cattle schistosomiasis and human cercarial dermatitis. This freshwater snail is widely distributed in Southern Asia, ranging from Iran to China eastwards including India and from the southeastern Himalayas to Southeast Asia southwards. The veterinary and medical importance of this snail explains the interest in understanding its geographical distribution patterns and evolutionary history. In this study, we used a large and comprehensive sampling throughout Indo-Malaya, including specimens from South India and Indonesia, areas that have been formerly less studied.

Results

The phylogenetic inference revealed five highly divergent clades (genetic distances among clades: 4.4–13.9%) that are morphologically indistinguishable, supporting the assumption that this presumed nominal species may represent a cryptic species complex. The species group may have originated in the humid subtropical plains of Nepal or in southern adjacent regions in the Early Miocene. The major cladogenetic events leading to the fives clades occurred successively from the Early Miocene to the Early Pleistocene, coinciding with major periods of monsoonal intensification associated with major regional paleogeographic events in the Miocene and repeated climate changes due to the Plio-Pleistocene climatic oscillations. Our coverage of the Indo-Australian Archipelago (IAA) highlights the presence of a single clade there. Contrary to expectations, an AMOVA did not reveal any population genetic structure among islands or along a widely recognised zoogeographical regional barrier, suggesting a recent colonisation independent of natural biogeographical constraints. Neutrality tests and mismatch distributions suggested a sudden demographic and spatial population expansion that could have occurred naturally in the Pleistocene or may possibly result of a modern colonisation triggered by anthropogenic activities.

Conclusions

Even though Indoplanorbis is the main focus of this study, our findings may also have important implications for fully understanding its role in hosting digenetic trematodes. The existence of a cryptic species complex, the historical phylogeographical patterns and the recent range expansion in the IAA provide meaningful insights to the understanding and monitoring of the parasites potential spread. It brings a substantial contribution to veterinary and public health issues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2043-6) contains supplementary material, which is available to authorized users.

Deep genetic structure and ecological divergence in a widespread human commensal toad

The Asian common toad (Duttaphrynus melanostictus) is a human commensal species that occupies a wide variety of habitats across tropical Southeast Asia. We test the hypothesis that genetic variation in D. melanostictus is weakly associated with geography owing to natural and human-mediated dispersal facilitated by its commensal nature. Phylogenetic and population genetic analyses of mitochondrial and nuclear DNA sequence variation, and predictive species distribution modelling, unexpectedly recovered three distinct evolutionary lineages that differ genetically and ecologically, corresponding to the Asian mainland, coastal Myanmar and the Sundaic islands. The persistence of these three divergent lineages, despite ample opportunities for recent human-mediated and geological dispersal, suggests that D. melanostictus actually consists of multiple species, each having narrower geographical ranges and ecological niches, and higher conservation value, than is currently recognized. These findings also have implications for the invasion potential of this human commensal elsewhere, such as in its recently introduced ranges on the islands of Borneo, Sulawesi, Seram and Madagascar.

Red River barrier and Pleistocene climatic fluctuations shaped the genetic structure of Microhyla fissipes complex (Anura: Microhylidae) in southern China and Indochina

South China and Indochina host striking species diversity and endemism. Complex tectonic and climatic evolutions appear to be the main drivers of the biogeographic patterns. In this study, based on the geologic history of this region, we test 2 hypotheses using the evolutionary history of Microhyla fissipes species complex. Using DNA sequence data from both mitochondrial and nuclear genes, we first test the hypothesis that the Red River is a barrier to gene flow and dispersal. Second, we test the hypothesis that Pleistocene climatic cycling affected the genetic structure and population history of these frogs. We detect 2 major genetic splits that associate with the Red River. Time estimation suggests that late Miocene tectonic movement associated with the Red River drove their diversification. Species distribution modeling (SDM) resolves significant ecological differences between sides of the Red River. Thus, ecological divergence also probably promoted and maintained the diversification. Genogeography, historical demography, and SDM associate patterns in southern China with climate changes of the last glacial maximum (LGM), but not Indochina. Differences in geography and climate between the 2 areas best explain the discovery. Responses to the Pleistocene glacial–interglacial cycling vary among species and regions.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Ecological genetics of Chinese rhesus macaque in response to mountain building: all things are not equal

BACKGROUND

Pliocene uplifting of the Qinghai-Tibetan Plateau (QTP) and Quaternary glaciation may have impacted the Asian biota more than any other events. Little is documented with respect to how the geological and climatological events influenced speciation as well as spatial and genetic structuring, especially in vertebrate endotherms. Macaca mulatta is the most widely distributed non-human primate. It may be the most suitable model to test hypotheses regarding the genetic consequences of orogenesis on an endotherm.

METHODOLOGY AND PRINCIPAL FINDINGS

Using a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, we discovered two maternal super-haplogroups exist, one in western China and the other in eastern China. M. mulatta formed around 2.31 Ma (1.51-3.15, 95%), and divergence of the two major matrilines was estimated at 1.15 Ma (0.78-1.55, 95%). The western super-haplogroup exhibits significant geographic structure. In contrast, the eastern super-haplogroup has far greater haplotypic variability with little structure based on analyses of six variable microsatellite loci using Structure and Geneland. Analysis using Migrate detected greater gene flow from WEST to EAST than vice versa. We did not detect signals of bottlenecking in most populations.

CONCLUSIONS

Analyses of the nuclear and mitochondrial datasets obtained large differences in genetic patterns for M. mulatta. The difference likely reflects inheritance mechanisms of the maternally inherited mtDNA genome versus nuclear biparentally inherited STRs and male-mediated gene flow. Dramatic environmental changes may be responsible for shaping the matrilineal history of macaques. The timing of events, the formation of M. mulatta, and the divergence of the super-haplogroups, corresponds to both the uplifting of the QTP and Quaternary climatic oscillations. Orogenesis likely drove divergence of western populations in China, and Pleistocene glaciations are likely responsible for genetic structuring in the eastern super-haplogroup via geographic isolation and secondary contact.