Phylogeography, geographic structure, genetic variation, and potential species boundaries in Philippine slender toads

Philippine herpetology (Amphibia, Reptilia), 20 years on: two decades of progress towards an increasingly collaborative, equitable, and inclusive approach to the study of the archipelago's amphibians and reptiles

A first review of the history, status, and prospects for Philippine herpetology conducted more than two decades ago (2002) summarized the diverse topics studied and highlighted the development and achievements in research up to the year 2000. This study revisits and re-assesses what Philippine herpetology has accomplished, both as a discipline and a community, during the last two decades (2002-2022). A total of 423 herpetological publications was collated, revealing a substantial increase in annual publications, rising from approximately four per year during 2002-2008 to around 28 per year in 2009-2022. Half of the published studies focused on squamate reptiles (lizards 30.5%, snakes 21%) and 28.4% on amphibians, 5.9% on turtles, and 2.6% on crocodiles. The remaining 11.6% of studies focused simultaneously on multiple taxa (i.e., faunal inventories). Diversity and distribution (35.2%) and ecological (26.5%) studies remained popular, while studies on taxonomy (14.9%), phylogenetics and biogeography (11.8%), and conservation (11.6%) all increased. However, geographical gaps persist urging immediate surveys in many understudied regions of the country. Finally, we found a balanced representation between Filipino and foreign first authors (1.0:1.1), yet a substantial gender gap exists between male and female first authors (7.1:1.0). Nonetheless, the steep increase in publications and the diversity of people engaged in Philippine herpetology is a remarkable positive finding compared to the 20 years preceding the last review (1980-2000). Our hope is that the next decades will bring increasingly equitable, internationally collaborative, and broadly inclusive engagement in the study of amphibians and reptiles in the Philippines.

Unexpectedly high levels of lineage diversity in Sundaland puddle frogs (Dicroglossidae: Occidozyga Kuhl and van Hasselt, 1822)

One of the most urgent contemporary tasks for taxonomists and evolutionary biologists is to estimate the number of species on earth. Recording alpha diversity is crucial for protecting biodiversity, especially in areas of elevated species richness, which coincide geographically with increased anthropogenic environmental pressures - the world's so-called biodiversity hotspots. Although the distribution of Puddle frogs of the genus Occidozyga in South and Southeast Asia includes five biodiversity hotspots, the available data on phylogeny, species diversity, and biogeography are surprisingly patchy. Samples analyzed in this study were collected throughout Southeast Asia, with a primary focus on Sundaland and the Philippines. A mitochondrial gene region comprising ~ 2000 bp of 12S and 16S rRNA with intervening tRNA Valine and three nuclear loci (BDNF, NTF3, POMC) were analyzed to obtain a robust, time-calibrated phylogenetic hypothesis. We found a surprisingly high level of genetic diversity within Occidozyga, based on uncorrected p-distance values corroborated by species delimitation analyses. This extensive genetic diversity revealed 29 evolutionary lineages, defined by the > 5% uncorrected p-distance criterion for the 16S rRNA gene, suggesting that species diversity in this clade of phenotypically homogeneous forms probably has been underestimated. The comparison with results of other anuran groups leads to the assumption that anuran species diversity could still be substantially underestimated in Southeast Asia in general. Many genetically divergent lineages of frogs are phenotypically similar, indicating a tendency towards extensive morphological conservatism. We present a biogeographic reconstruction of the colonization of Sundaland and nearby islands which, together with our temporal framework, suggests that lineage diversification centered on the landmasses of the northern Sunda Shelf. This remarkably genetically structured group of amphibians could represent an exceptional case for future studies of geographical structure and diversification in a widespread anuran clade spanning some of the most pronounced geographical barriers on the planet (e.g., Wallace's Line). Studies considering gene flow, morphology, ecological and bioacoustic data are needed to answer these questions and to test whether observed diversity of Puddle frog lineages warrants taxonomic recognition.

Gene Flow Increases Phylogenetic Structure and Inflates Cryptic Species Estimations: A Case Study on Widespread Philippine Puddle Frogs (Occidozyga laevis)

In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs-a group that is widespread, phenotypically conserved, and exhibits high levels of geographically-based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: (1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and (2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear p-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely-adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow.

Dynamics behind disjunct distribution, hotspot-edge refugia, and discordant RADseq/mtDNA variability: insights from the Emei mustache toad

Background

The distribution of genetic diversity and the underlying processes are important for conservation planning but are unknown for most species and have not been well studied in many regions. In East Asia, the Sichuan Basin and surrounding mountains constitute an understudied region that exhibits a “ring” of high species richness overlapping the eastern edge of the global biodiversity hotspot Mountains of Southwest China. We examine the distributional history and genetic diversification of the Emei mustache toad Leptobrachium boringii, a typical “ring” element characterized by disjunct ranges in the mountains, by integrating time-calibrated gene tree, genetic variability, individual-level clustering, inference of population splitting and mixing from allele frequencies, and paleoclimatic suitability modeling.

Results

The results reveal extensive range dynamics, including secondary contact after long-term isolation via westward dispersal accompanied by variability loss. They allow the proposal of a model that combines recurrent contractions caused by Quaternary climatic changes and some failed expansions under suitable conditions for explaining the shared disjunct distribution pattern. Providing exceptional low-elevation habitats in the hotspot area, the eastern edge harbors both long-term refugial and young immigrant populations. This finding and a synthesis of evidence from other taxa demonstrate that a certain contributor to biodiversity, one that preserves and receives low-elevation elements of the east in this case, can be significant for only a particular part of a hotspot. By clarifying the low variability of these refugial populations, we show that discordant mitochondrial estimates of diversity can be obtained for populations that experienced admixture, which would have unlikely left proportional immigrant alleles for each locus.

Conclusions

Dispersal after long-term isolation can explain much of the spatial distribution of genetic diversity in this species, while secondary contact and long-term persistence do not guarantee a large variation. The model for the formation of disjunct ranges may apply to many other taxa isolated in the mountains surrounding the Sichuan Basin. Furthermore, this study provides insights into the heterogeneous nature of hotspots and discordant variability obtained from genome-wide and mitochondrial data.

From mainland to islands: colonization history in the tree frog Kurixalus (Anura: Rhacophoridae)

The origin and colonization history of Kurixalus, a genus of small arboreal tree frogs breeding exclusively in shallow swamps, is under disputed. On the basis of comprehensive sampling program, the evolutionary history of Kurixalus is investigated based on 3 mitochondrial genes. Our results indicate that the genus Kurixalus originated in the Asian mainland and subsequently arrived at its current distribution in Borneo, Taiwan, Ryukyu, and Hainan islands by a series of dispersal events. Moreover, the colonization of Taiwan from mainland Asia has occurred 2 times. The initial colonization of Taiwan occurred at 3.46–8.68 Mya (95% highest posterior density), which rejects the hypothesis that Kurixalus probably originated from Taiwan during the early Oligocene and favors the model of Neogene-origin rather than the model of Quaternary-origin for Taiwanese Kurixalus. Kurixalus eiffingeri has dispersed from Taiwan to the Ryukyus once or 2 times pending more data. Both transoceanic dispersal and landbridge dispersal have played a role in the colonization process; the former resulted in the colonization of Taiwan and the Ryukyus and the latter led to the colonization of Borneo and Hainan.

Examining the interglacial high-elevation refugia scenario in East Asian subtropical mountain systems with the frog species Leptobrachium liui

The effects of Quaternary climatic oscillations on the distributions of organisms in different parts of the world are not equally well understood, limiting the ability to understand the determinants of biodiversity. Compared with the mountain regions in southern Europe and southwestern North America, such effects on high-elevation species in the East Asian subtropical mountain systems located in southern and southeastern China have seldom been addressed. In this study, using Leptobrachium liui (Megophryidae), we made one of the earliest attempts to examine the interglacial high-elevation refugia scenario in these Asian mountains. Based on our current understanding of the study system, we formulated a hypothesis that these frogs of western origin were distributed more widely and continuously during glacial phases, allowing eastward dispersal, and that they are currently isolated in interglacial refugia at higher elevations. Microsatellite data and mitochondrial and nuclear sequence data were obtained with extensive sampling followed by the synthesis of phylogeographic and population genetic analyses and modeling of the species distribution. The analyses revealed a sequential eastward divergence of microsatellite clusters and gene lineages accompanied by a decline in genetic diversity. Molecular dating estimates revealed divergence events in the Pleistocene, and a reduction in local populations was inferred to have occurred at a time comparable to the end of the last glacial. Strong genetic isolation by distance reflecting a more continuous historical distribution was detected. Furthermore, environmental niche models inferred a wide planar distribution during the last glacial maximum, providing further support for the hypothesis.

Colonization and diversification of the white-browed shortwing (Aves: Muscicapidae: Brachypteryx montana) in the Philippines

Molecular phylogenetic approaches have greatly improved our knowledge of the pattern and process of biological diversification across the globe; however, many regions remain poorly documented, even for well-studied vertebrate taxa. The Philippine archipelago, one of the least-studied 'biodiversity hotspots', is an ideal natural laboratory for investigating the factors driving diversification in an insular and geologically dynamic setting. We investigated the history and geography of diversification of the Philippine populations of a widespread montane bird, the White-browed Shortwing (Brachypteryx montana). Leveraging dense archipelago-wide sampling, we generated a multi-locus genetic dataset (one nuclear and two mtDNA markers), which we analyzed using phylogenetic, population genetic, and coalescent-based methods. Our results demonstrate that Philippine shortwings (1) likely colonized the Philippines from the Sunda Shelf to Mindanao in the late Miocene or Pliocene, (2) diversified across inter-island barriers into three divergent lineages during the Pliocene and early Pleistocene, (3) have not diversified within the largest island, Luzon, contrary to patterns observed in other montane taxa, and (4) colonized Palawan from the oceanic Philippines rather than from Borneo, challenging the assumption of Palawan functioning exclusively as a biogeographic extension of the Sunda Shelf. Additionally, our finding that divergent (c. 4.0 mya) lineages are coexisting in secondary sympatry on Mindanao without apparent gene flow suggests that the speciation process is likely complete for these shortwing lineages. Overall, these investigations provide insight into how topography and island boundaries influence diversification within remote oceanic archipelagos and echo the results of many other studies in demonstrating that taxonomic diversity continues to be underestimated in the Philippines.

The amphibians and reptiles of Mindanao Island, southern Philippines, II: the herpetofauna of northeast Mindanao and adjacent islands

We summarize all available amphibian and reptile species distribution data from the northeast Mindanao faunal region, including small islands associated with this subcenter of endemic vertebrate biodiversity. Together with all publicly available historical information from biodiversity repositories, we present new data from several major herpetological surveys, including recently conducted inventories on four major mountains of northeast Mindanao, and adjacent islands of Camiguin Sur, Dinagat, and Siargao. We present species accounts for all taxa, comment on unresolved taxonomic problems, and provide revisions to outdated IUCN conservation status assessments in cases where our new data significantly alter earlier classification status summaries. Together, our comprehensive analysis of this fauna suggests that the greater Mindanao faunal region possesses distinct subcenters of amphibian and reptile species diversity, and that until this area is revisited and its fauna and actually studied, with on-the-ground field work including targeted surveys of species distributions coupled to the study their natural history, our understanding of the diversity and conservation status of southern Philippine herpetological fauna will remain incomplete. Nevertheless, the northeast Mindanao geographical area (Caraga Region) appears to have the highest herpetological species diversity (at least 126 species) of any comparably-sized Philippine faunal subregion.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.

Genetic Differentiation in Insular Lowland Rainforests: Insights from Historical Demographic Patterns in Philippine Birds

Phylogeographic studies of Philippine birds support that deep genetic structure occurs across continuous lowland forests within islands, despite the lack of obvious contemporary isolation mechanisms. To examine the pattern and tempo of diversification within Philippine island forests, and test if common mechanisms are responsible for observed differentiation, we focused on three co-distributed lowland bird taxa endemic to Greater Luzon and Greater Negros-Panay: Blue-headed Fantail (Rhipidura cyaniceps), White-browed Shama (Copsychus luzoniensis), and Lemon-throated Leaf-Warbler (Phylloscopus cebuensis). Each species has two described subspecies within Greater Luzon, and a single described subspecies on Greater Negros/Panay. Each of the three focal species showed a common geographic pattern of two monophyletic groups in Greater Luzon sister to a third monophyletic group found in Greater Negros-Panay, suggesting that common or similar biogeographic processes may have produced similar distributions. However, studied species displayed variable levels of mitochondrial DNA differentiation between clades, and genetic differentiation within Luzon was not necessarily concordant with described subspecies boundaries. Population genetic parameters for the three species suggested both rapid population growth from small numbers and geographic expansion across Luzon Island. Estimates of the timing of population expansion further supported that these events occurred asynchronously throughout the Pleistocene in the focal species, demanding particular explanations for differentiation, and support that co-distribution may be secondarily congruent.

Archipelago colonization by ecologically dissimilar amphibians: evaluating the expectation of common evolutionary history of geographical diffusion in co-distributed rainforest tree frogs in islands of Southeast Asia

Widespread, co-distributed species with limited relative dispersal abilities represent compelling focal taxa for comparative phylogeography. Forest vertebrates in island archipelagos often exhibit pronounced population structure resulting from limited dispersal abilities or capacity to overcome marine barriers to dispersal. The exceptionally diverse Old World tree frogs of the family Rhacophoridae have colonized the forested island archipelagos of Southeast Asia on multiple occasions, entering the islands of Indonesia and the Philippines via a "stepping stone" mode of dispersal along elongate island chains, separated by a series of marine channels. Here we evaluate the prediction that two tightly co-distributed Philippine rhacophorids colonized the archipelago during concomitant timescales and in the same, linear, "island-hopping" progression. We use a new multilocus dataset, utilize dense genetic sampling from the eastern arc of the Philippines, and we take a model-based phylogeographic approach to examining the two species for similar topological patterns of diversification, genetic structure, and timescales of diversification. Our results support some common mechanistic predictions (a general south-to-north polarity of colonization) but not others (timescale for colonization and manner and degree of lineage diversification), suggesting differing biogeographic scenarios of geographical diffusion through the archipelago and unique and idiosyncratic ecological capacities and evolutionary histories of each species.

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae)

New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490 base pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.