Abstract: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. [Cryptic species; gene flow; introgression; isolation-by-distance; isolation-by-environment; phylogenetic network; species delimitation.].

The effect of gene flow on coalescent-based species-tree inference

Genealogical asymmetry under the IM model and a two-taxon test for gene flow

Identifiability of speciation times under the multispecies coalescent

Inference of continuous gene flow between species under misspecified models

The Inference of Gene Trees with Species Trees

The accuracy of species tree inference under gene tree dependence

From gene trees to species trees II: Species tree inference in the deep coalescence model

Coalescent-based species tree estimation: a stochastic Farris transform

Anomaly zones for uniformly sampled gene trees under the gene duplication and loss model

Species Delimitation with Gene Flow

Computing the probability of gene trees concordant with the species tree in the multispecies coalescent

Sympatric speciation based on pure assortative mating

Bayesian Inference of Species Trees from Multilocus Data

On the connections between the spatial Lambda–Fleming–Viot model and other processes for analysing geo-referenced genetic data

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

Inferring Species Trees Directly from Biallelic Genetic Markers: Bypassing Gene Trees in a Full Coalescent Analysis

Estimating Species Trees from Quartet Gene Tree Distributions under the Coalescent Model

Efficient Bayesian species tree inference under the multi-species coalescent

ASTRAL: genome-scale coalescent-based species tree estimation

What Predicts Gene Flow During Speciation? The Relative Roles of Time, Space, Morphology and Climate