The Importance of Contact Zones for Distinguishing Interspecific from Intraspecific Geographic Variation

Distinguishing species boundaries from geographic variation

In an era of unprecedented biodiversity loss, the need for standardized practices to describe biological variation is becoming increasingly important. As with all scientific endeavors, species delimitation needs to be explicit, testable, and refutable. A fundamental task in species delimitation is distinguishing within-species variation from among-species variation. Many species that are distributed across large geographic areas exhibit levels of genetic variation that are as great or greater than those that exist between well-defined sympatric species. Here, we provide a workflow to distinguish between intra- and interspecific genetic variation and apply the workflow to a taxonomically problematic group of frogs (the Rana pipiens complex, or leopard frogs) that are widely distributed across Mexico and Central America. Our workflow makes use of recent advancements that pair genome-scale datasets with model-based species delimitation methods, while emphasizing the need for positive evidence of reproductive isolation to confirm the validity of geographically contiguous species boundaries. We find that intraspecific geographic variation in widespread leopard frog species has resulted in considerable taxonomic inflation of species. Ten currently recognized species are not supported in our analyses, and we here synonymize them with previously named taxa. Furthermore, we find positive evidence for the presence of three undescribed species. In addition to proposing these taxonomic changes, we provide descriptions of the data or analyses that would be needed to refute and overturn our recommendations. We recommend that all species delimitation studies (especially of geographically variable groups) clarify what new evidence would be sufficient to change the taxonomic recommendations.

Deforestation-induced Hybridization in Philippine Frogs Creates a Distinct Phenotype With an Inviable Genotype

Hybridization plays a major role in the evolutionary history of many taxa and can generate confounding patterns affecting many downstream applications. In this study, we empirically demonstrate how hybridization obfuscates phylogenetic inference (via the artefactual branch effect), species boundaries, and taxonomy in an adaptive radiation of frogs. Philippine narrow-mouthed frogs of the genus Kaloula exhibit a wide range of phenotypic and ecological adaptations but their evolutionary history and taxonomy remain poorly understood. In particular, the Kaloula conjuncta complex contains numerous subspecies with unresolved taxonomic boundaries and unclear evolutionary relationships. Within this complex, Kaloula conjuncta stickeli, until now was considered a rare, enigmatic, and phenotypically distinct subspecies that had not been encountered since its original description nearly 80 years ago. Here, we show that K. c. stickeli shares alleles with K. conjuncta meridionalis and another species outside the conjuncta group, K. picta. Using target-capture sequencing and a robust analytical framework, we show that despite having a unique phenotype, K. c. stickeli is likely an inviable F1 hybrid between K. c. meridionalis and K. picta and thus, does not warrant taxonomic recognition. Our results show how industry-standard approaches in systematic inference and integrative taxonomy-morphological, phylogenomic, clustering, and distance-based methods-can generate misleading results for identifying and understanding affinities of hybrids. In contrast, we demonstrate how network multispecies coalescent and population genetic approaches are more effective at accurately inferring reticulated evolutionary history. We also propose a rare phenomenon of deforestation-induced hybridization, which could have important consequences in light of large-scale Southeast Asian forest destruction.

Genomic, Phenotypic and Environmental Correlates of Speciation in the Midwife Toads (Alytes)

Speciation, i.e., the formation of new species, implies that diverging populations evolve genetic, phenotypic or ecological factors that promote reproductive isolation (RI), but the relative contributions of these factors remain elusive. Here we test which of genomic, bioacoustic, morphological, and environmental differences best predicts RI across a continuum of divergence in the midwife toads (genus Alytes), a group of Western Mediterranean amphibians, using a total evidence approach. We found that, without strong geographic barriers to dispersal, the extent of introgression across hybrid zones between phylogeographic lineages, which should reflect the strength of RI, predominantly covaries with genomic divergence. Overall phenotypic differentiation becomes substantial only between well established, fully isolated species. These results suggest that speciation in midwife toads initially involve cryptic lineages, which probably evolve RI through intrinsic (genetic) hybrid incompatibilities. As they continue to diverge, these nascent species eventually differentiate externally, which potentially enforces pre-mating barriers and facilitates sympatry. This speciation scenario has practical implications for species delimitation, notably when using hybrid zones and divergence thresholds as proxies for reproductive isolation.

Phylogenomics yields new systematic and taxonomical insights for Lissotriton newts, a genus with a strong legacy of introgressive hybridization

The ease with which genome-wide data can nowadays be collected allows complicated phylogenetic questions to be re-evaluated. Phylogenetic relationships among newts have often proven difficult to resolve due to the prevalence of incomplete lineage sorting and introgressive hybridization. For the newt genus Lissotriton, phylogenetic relationships are not settled and there is controversy surrounding the species status of several taxa. We obtain c. 7 k nuclear DNA markers with target enrichment by sequence capture and conduct a concatenated analysis with RAxML, gene-tree summarization with ASTRAL, and species tree estimation with SNAPPER. We explore introgression between evolutionary lineages with TreeMix and Dsuite and compare how introgression events influence the different phylogenetic tools employed. We retrieve tree topologies that are discordant with previous mtDNA-based attempts, particularly concerning the phylogenetic placement of L. italicus and the L. vulgaris species complex. Yet, we also observe deviations between the phylogenetic hypotheses resulting from the different analyses. We interpret the placement of L. montandoni deep within the L. vulgaris species complex by SNAPPER, rather than as the sister taxon to the remainder of the L. vulgaris species complex according to RAxML and ASTRAL, as an artifact of introgression - well-documented in previous work and backed up by TreeMix and Dsuite analyses. Our analyses allow us to make some taxonomical recommendations: we confirm the recently proposed species status of L. kosswigi and L. graecus and propose that L. v. lantzi and L. v. schmidtleri had better be treated as subspecies. Our work also highlights areas for further taxonomic research: range-wide phylogenomic data are required to disentangle the L. boscai - L. maltzani species complex and the northern and southern lineages of L. v. vulgaris. Our study illustrates the power of target enrichment by sequence capture in tackling longstanding questions in taxa with an extensive history of hybridization and introgression.

Cryptic speciation in arid mountains: An integrative revision of the Pristurus rupestris species complex (Squamata, Sphaerodactylidae) from Arabia based on morphological, genetic and genomic data, with the description of four new species

In the arid landscapes of the Arabian Peninsula, high levels of cryptic diversity among reptiles, and especially in geckos, have recently been revealed. Mountain ranges within the peninsula were shown to contain the highest richness of reptile endemicity, serving as refugia to species less adapted to the hyper-arid conditions of the lowlands. With up to 19 endemic reptile species, the Hajar Mountains of southeastern Arabia are a clear example of this pattern. Owing to its old geological history, complex topography and geographic isolation from the rest of the peninsula, this mountain range rises as a hotspot of reptile biodiversity and endemicity in Arabia, and provides the perfect scenario to study the processes of evolution and diversification of reptiles in arid mountain ranges. In the present study we investigate the systematics of the Pristurus rupestris species complex, a group of geckos exhibiting cryptic morphological traits along with a remarkably deep evolutionary history. Initially considered a single species distributed throughout coastal Arabia, and with some scattered populations at the Horn of Africa, several recent studies have shown that Pristurus rupestris actually comprises a species complex restricted to the Hajar Mountains of southeastern Arabia. Here, we utilize an integrative approach assembling several morphological, genetic, genomic, and ecological datasets to resolve this long-standing systematic challenge. Results support the existence of four new cryptic Pristurus species in the Hajar Mountains with three new Oman endemics. While no unique diagnostic morphological characters were identified, some slight morphological differences occur between species, especially among high-elevation species relative to the rest. Despite the lack of clear morphological differentiation, extreme levels of genetic variation were found between species with genetic distances of up to 24% in the 12S mitochondrial marker, resulting from deep divergence times of up to 10 mya. Moreover, all species have been found in sympatry with at least another representative of the species complex and without any signs of apparent and ongoing gene flow among them. These findings yield profound implications for conservation efforts, as one of these newly described species presents an extremely restricted distribution (only known from a single locality and three individuals), requiring immediate attention for protection. Overall, this study sheds light on the hidden diversity within the P. rupestris species complex, emphasizing the importance of preserving biodiversity in the face of ongoing environmental changes, while highlighting, once again, the Hajar Mountains of southeastern Arabia as a cradle of reptile biodiversity.

Global phylogeny of the family Gomphillaceae (Ascomycota, Graphidales) sheds light on the origin, diversification and endemism in foliicolous lineages

Foliicolous lichens grow on living leaves of vascular plants. They are mostly found in tropical to subtropical or temperate rainforests. Many phenotype-based species are considered as pantropical or even sub-cosmopolitan, either attributed to old ages, having existed prior to continental breakups or long-distance dispersal. We built a much expanded, global phylogeny of Gomphillaceae, the most diverse group of leaf-dwelling lichenised fungi. Our sampling encompassed six major biodiversity hotspots: MIOI (Madagascar and the Indian Ocean Islands), the Caribbean, New Caledonia, the Colombian Chocó, Mesoamerica and the Atlantic coast of Brazil. It was based on multilocus sequence data (mtSSU rDNA, nuLSU rDNA and RPB1), including 2207 sequences of 1256 specimens. Species delimitation methods combined with a phenotype matrix identified 473 putative species. Amongst these, 104 are confirmed as described, 213 are classified as cryptic or near cryptic (hidden diversity), 100 represent new species to science (identified on the basis of phenotype) and 56 remain unidentified. Amongst the 104 species with a valid name, 40.5% are distributed across 2-5 continents (lichenogeographical regions) by applying the phenotype-based species concept. However, using the integrative approach to delineate species, this estimate is reduced to 9%. We estimate the global species richness of Gomphillaceae at 1,861-2,356 species. The timing of species-level divergences suggests that the current distribution of foliicolous lichens is shaped more by long-distance dispersal and rapid diversification than by vicariance. The origin of the family and major clades appears to be in the Neotropics, with subsequent numerous dispersal events. Our results support the separation of three major lineages, corresponding to the former families Asterothyriaceae, Gomphillaceae s.str. and Solorinellaceae, which should be recognised at the subfamily level.

The Limits of the Metapopulation: Lineage Fragmentation in a Widespread Terrestrial Salamander (Plethodon cinereus)

In vicariant species formation, divergence results primarily from periods of allopatry and restricted gene flow. Widespread species harboring differentiated, geographically distinct sublineages offer a window into what may be a common mode of species formation, whereby a species originates, spreads across the landscape, then fragments into multiple units. However, incipient lineages usually lack reproductive barriers that prevent their fusion upon secondary contact, blurring the boundaries between a single, large metapopulation-level lineage and multiple independent species. Here, we explore this model of species formation in the Eastern Red-backed Salamander (Plethodon cinereus), a widespread terrestrial vertebrate with at least 6 divergent mitochondrial clades throughout its range. Using anchored hybrid enrichment data, we applied phylogenomic and population genomic approaches to investigate patterns of divergence, gene flow, and secondary contact. Genomic data broadly match most mitochondrial groups but reveal mitochondrial introgression and extensive admixture at several contact zones. While species delimitation analyses in Bayesian Phylogenetics and Phylogeography supported 5 lineages of P. cinereus, genealogical divergence indices (gdi) were highly sensitive to the inclusion of admixed samples and the geographic representation of candidate species, with increasing support for multiple species when removing admixed samples or limiting sampling to a single locality per group. An analysis of morphometric data revealed differences in body size and limb proportions among groups, with a reduction of forelimb length among warmer and drier localities consistent with increased fossoriality. We conclude that P. cinereus is a single species, but one with highly structured component lineages of various degrees of independence.

Secondary contact erodes Pleistocene diversification in a wide-ranging freshwater mussel (Quadrula)

The isolated river drainages of eastern North America serve as a natural laboratory to investigate the roles of allopatry and secondary contact in the evolutionary trajectories of recently diverged lineages. Drainage divides facilitate allopatric speciation, but due to their sensitivity to climatic and geomorphological changes, neighboring rivers frequently coalesce, creating recurrent opportunities of isolation and contact throughout the history of aquatic lineages. The freshwater mussel Quadrula quadrula is widely distributed across isolated rivers of eastern North America and possesses high phenotypic and molecular variation across its range. We integrate sequence data from three genomes, including female- and male-inherited mitochondrial markers and thousands of nuclear encoded SNPs with morphology and geography to illuminate the group's divergence history. Across contemporary isolated rivers, we found continuums of molecular and morphological variation, following a pattern of isolation by distance. In contact zones, hybridization was frequent with no apparent fitness consequences, as advanced hybrids were common. Accordingly, we recognize Q. quadrula as a single cohesive species with subspecific variation (Q. quadrula rumphiana). Demographic modeling and divergence dating supported a divergence history characterized by allopatric vicariance followed by secondary contact, likely driven by river rearrangements and Pleistocene glacial cycles. Despite clinal range-wide variation and hybridization in contact zones, the process-based species delimitation tool delimitR, which considers demographic scenarios like secondary contact, supported the delimitation of the maximum number of species tested. As such, when interpreting species delimitation results, we suggest careful consideration of spatial sampling and subsequent geographic patterns of biological variation, particularly for wide-ranging taxa.

Understanding Species Boundaries that Arise from Complex Histories: Gene Flow Across the Speciation Continuum in the Spotted Whiptail Lizards

-Gene flow between diverging lineages challenges the resolution of species boundaries and the understanding of evolutionary history in recent radiations. Here, we integrate phylogenetic and coalescent tools to resolve reticulate patterns of diversification and use a perspective focused on evolutionary mechanisms to distinguish interspecific and intraspecific taxonomic variation. We use this approach to resolve the systematics for one of the most intensively studied but difficult to understand groups of reptiles: the spotted whiptail lizards of the genus Aspidoscelis (A. gularis complex). Whiptails contain the largest number of unisexual species known within any vertebrate group and the spotted whiptail complex has played a key role in the generation of this diversity through hybrid speciation. Understanding lineage boundaries and the evolutionary history of divergence and reticulation within this group is therefore key to understanding the generation of unisexual diversity in whiptails. Despite this importance, long-standing confusion about their systematics has impeded understanding of which gonochoristic species have contributed to the formation of unisexual lineages. Using reduced representation genomic data, we resolve patterns of divergence and gene flow within the spotted whiptails and clarify patterns of hybrid speciation. We find evidence that biogeographically structured ecological and environmental variation has been important in morphological and genetic diversification, as well as the maintenance of species boundaries in this system. Our study elucidates how gene flow among lineages and the continuous nature of speciation can bias the practice of species delimitation and lead taxonomists operating under different frameworks to different conclusions (here we propose that a 2 species arrangement best reflects our current understanding). In doing so, this study provides conceptual and methodological insights into approaches to resolving diversification patterns and species boundaries in rapid radiations with complex histories, as well as long-standing taxonomic challenges in the field of systematic biology.

Understanding species limits through the formation of phylogeographic lineages

The outcomes of speciation across organismal dimensions (e.g., ecological, genetic, phenotypic) are often assessed using phylogeographic methods. At one extreme, reproductively isolated lineages represent easily delimitable species differing in many or all dimensions, and at the other, geographically distinct genetic segments introgress across broad environmental gradients with limited phenotypic disparity. In the ambiguous gray zone of speciation, where lineages are genetically delimitable but still interacting ecologically, it is expected that these lineages represent species in the context of ontology and the evolutionary species concept when they are maintained over time with geographically well-defined hybrid zones, particularly at the intersection of distinct environments. As a result, genetic structure is correlated with environmental differences and not space alone, and a subset of genes fail to introgress across these zones as underlying genomic differences accumulate. We present a set of tests that synthesize species delimitation with the speciation process. We can thereby assess historical demographics and diversification processes while understanding how lineages are maintained through space and time by exploring spatial and genome clines, genotype-environment interactions, and genome scans for selected loci. Employing these tests in eight lineage-pairs of snakes in North America, we show that six pairs represent 12 "good" species and that two pairs represent local adaptation and regional population structure. The distinct species pairs all have the signature of divergence before or near the mid-Pleistocene, often with low migration, stable hybrid zones of varying size, and a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between distinct ecoregions. Locally adapted populations are younger, exhibit higher migration, and less ecological differentiation. Our results demonstrate that interacting lineages can be delimited using phylogeographic and population genetic methods that properly integrate spatial, temporal, and environmental data.

Establishing species boundaries in Bornean geckos

Species delimitation using mitochondrial DNA (mtDNA) remains an important and accessible approach for discovering and delimiting species. However, delimiting species with a single locus (e.g. DNA barcoding) is biased towards overestimating species diversity. The highly diverse gecko genus Cyrtodactylus is one such group where delimitation using mtDNA remains the paradigm. In this study, we use genomic data to test putative species boundaries established using mtDNA within three recognized species of Cyrtodactylus on the island of Borneo. We predict that multi-locus genomic data will estimate fewer species than mtDNA, which could have important ramifications for the species diversity within the genus. We aim to (i) investigate the correspondence between species delimitations using mtDNA and genomic data, (ii) infer species trees for each target species, and (iii) quantify gene flow and identify migration patterns to assess population connectivity. We find that species diversity is overestimated and that species boundaries differ between mtDNA and nuclear data. This underscores the value of using genomic data to reassess mtDNA-based species delimitations for taxa lacking clear species boundaries. We expect the number of recognized species within Cyrtodactylus to continue increasing, but, when possible, genomic data should be included to inform more accurate species boundaries.

A cryptic syngameon within Betula shrubs revealed: Implications for conservation in changing subarctic environments

Arctic and subarctic ecosystems are rapidly transforming due to global warming, emphasizing the need to understand the genetic diversity and adaptive strategies of northern plant species for effective conservation. This study focuses on Betula glandulosa, a native North American tundra shrub known as dwarf birch, which demonstrates an apparent capacity to adapt to changing climate conditions. To address the taxonomic challenges associated with shrub birches and logistical complexities of sampling in the northernmost areas where species' ranges overlap, we adopted a multicriteria approach. Incorporating molecular data, ploidy level assessment and leaf morphology, we aimed to distinguish B. glandulosa individuals from other shrub birch species sampled. Our results revealed three distinct species and their hybrids within the 537 collected samples, suggesting the existence of a shrub birch syngameon, a reproductive network of interconnected species. Additionally, we identified two discrete genetic clusters within the core species, B. glandulosa, that likely correspond to two different glacial lineages. A comparison between the nuclear and chloroplast SNP data emphasizes a long history of gene exchange between different birch species and genetic clusters. Furthermore, our results highlight the significance of incorporating interfertile congeneric species in conservation strategies and underscores the need for a holistic approach to conservation in the context of climate change, considering the complex dynamics of species interactions. While further research will be needed to describe this shrub birches syngameon and its constituents, this study is a first step in recognizing its existence and disseminating awareness among ecologists and conservation practitioners. This biological phenomenon, which offers evolutionary flexibility and resilience beyond what its constituent species can achieve individually, may have significant ecological implications.

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Microgeographic population structuring in a genus of California trapdoor spiders and discovery of an enigmatic new species (Euctenizidae: Promyrmekiaphila korematsui sp. nov.)

The recognition and delineation of cryptic species remains a perplexing problem in systematics, evolution, and species delimitation. Once recognized as such, cryptic species complexes provide fertile ground for studying genetic divergence within the context of phenotypic and ecological divergence (or lack thereof). Herein we document the discovery of a new cryptic species of trapdoor spider, Promyrmekiaphila korematsui sp. nov. Using subgenomic data obtained via target enrichment, we document the phylogeography of the California endemic genus Promyrmekiaphila and its constituent species, which also includes P. clathrata and P. winnemem. Based on these data we show a pattern of strong geographic structuring among populations but cannot entirely discount recent gene flow among populations that are parapatric, particularly for deeply diverged lineages within P. clathrata. The genetic data, in addition to revealing a new undescribed species, also allude to a pattern of potential phenotypic differentiation where species likely come into close contact. Alternatively, phenotypic cohesion among genetically divergent P. clathrata lineages suggests that some level of gene flow is ongoing or occurred in the recent past. Despite considerable field collection efforts over many years, additional sampling in potential zones of contact for both species and lineages is needed to completely resolve the dynamics of divergence in Promyrmekiaphila at the population-species interface.

Reproductive isolation in a three-way contact zone

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.