Candidate-species delimitation in Desmognathus salamanders reveals gene flow across lineage boundaries, confounding phylogenetic estimation and clarifying hybrid zones

Phylogenomic Discordance is Driven by Wide-Spread Introgression and Incomplete Lineage Sorting During Rapid Species Diversification Within Rattlesnakes (Viperidae: Crotalus and Sistrurus)

-Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here, we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like those seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

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.

Novel phylogenomic inference and 'Out of Asia' biogeography of cobras, coral snakes and their allies

Estimation of evolutionary relationships among lineages that rapidly diversified can be challenging, and, in such instances, inaccurate or unresolved phylogenetic estimates can lead to erroneous conclusions regarding historical geographical ranges of lineages. One example underscoring this issue has been the historical challenge posed by untangling the biogeographic origin of elapoid snakes, which includes numerous dangerously venomous species as well as species not known to be dangerous to humans. The worldwide distribution of this lineage makes it an ideal group for testing hypotheses related to historical faunal exchanges among the many continents and other landmasses occupied by contemporary elapoid species. We developed a novel suite of genomic resources, included worldwide sampling, and inferred a robust estimate of evolutionary relationships, which we leveraged to quantitatively estimate geographical range evolution through the deep-time history of this remarkable radiation. Our phylogenetic and biogeographical estimates of historical ranges definitively reject a lingering former 'Out of Africa' hypothesis and support an 'Out of Asia' scenario involving multiple faunal exchanges between Asia, Africa, Australasia, the Americas and Europe.

River Drainage Reorganization and Reticulate Evolution in the Two-Lined Salamander (Eurycea bislineata) Species Complex

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge." Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for > 100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical. [Amphibian; hybridization; introgression; Plethodontidae; stream capture.].

The natural and human-mediated expansion of a human-commensal lizard into the fringes of Southeast Asia

BACKGROUND

Human-commensal species often display deep ancestral genetic structure within their native range and founder-effects and/or evidence of multiple introductions and admixture in newly established areas. We investigated the phylogeography of Eutropis multifasciata, an abundant human-commensal scincid lizard that occurs across Southeast Asia, to determine the extent of its native range and to assess the sources and signatures of human introduction outside of the native range. We sequenced over 350 samples of E. multifasciata for the mitochondrial ND2 gene and reanalyzed a previous RADseq population genetic dataset in a phylogenetic framework.

RESULTS

Nuclear and mitochondrial trees are concordant and show that E. multifasciata has retained high levels of genetic structure across Southeast Asia despite being frequently moved by humans. Lineage boundaries in the native range roughly correspond to several major biogeographic barriers, including Wallace's Line and the Isthmus of Kra. Islands at the outer fringe of the range show evidence of founder-effects and multiple introductions.

CONCLUSIONS

Most of enormous range of E. multifasciata across Southeast Asia is native and it only displays signs of human-introduction or recent expansion along the eastern and northern fringe of its range. There were at least three events of human-introductions to Taiwan and offshore islands, and several oceanic islands in eastern Indonesia show a similar pattern. In Myanmar and Hainan, there is a founder-effect consistent with post-warming expansion after the last glacial maxima or human introduction.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

The Artefactual Branch Effect and Phylogenetic Conflict: Species Delimitation with Gene Flow in Mangrove Pit Vipers (Trimeresurus purpureomaculatus-erythrurus Complex)

Mangrove pit vipers of the Trimeresurus purpureomaculatus-erythrurus complex are the only species of viper known to naturally inhabit mangroves. Despite serving integral ecological functions in mangrove ecosystems, the evolutionary history, distribution, and species boundaries of mangrove pit vipers remain poorly understood, partly due to overlapping distributions, confusing phenotypic variations, and the lack of focused studies. Here, we present the first genomic study on mangrove pit vipers and introduce a robust hypothesis-driven species delimitation framework that considers gene flow and phylogenetic uncertainty in conjunction with a novel application of a new class of speciation-based delimitation model implemented through the program Delineate. Our results showed that gene flow produced phylogenetic conflict in our focal species and substantiates the artefactual branch effect where highly admixed populations appear as divergent nonmonophyletic lineages arranged in a stepwise manner at the basal position of clades. Despite the confounding effects of gene flow, we were able to obtain unequivocal support for the recognition of a new species based on the intersection and congruence of multiple lines of evidence. This study demonstrates that an integrative hypothesis-driven approach predicated on the consideration of multiple plausible evolutionary histories, population structure/differentiation, gene flow, and the implementation of a speciation-based delimitation model can effectively delimit species in the presence of gene flow and phylogenetic conflict.

Unsupervised machine learning for species delimitation, integrative taxonomy, and biodiversity conservation

Integrative taxonomy, combining data from multiple axes of biologically relevant variation, is a major goal of systematics. Ideally, such taxonomies will derive from similarly integrative species-delimitation analyses. Yet, most current methods rely solely or primarily on molecular data, with other layers often incorporated only in a post hoc qualitative or comparative manner. A major limitation is the difficulty of devising quantitative parametric models linking different datasets in a unified ecological and evolutionary framework. Machine Learning (ML) methods offer flexibility in this arena by easily learning high-dimensional associations between observations (e.g., individual specimens) across a wide array of input features (e.g., genetics, geography, environment, and phenotype) to delimit statistically meaningful clusters. Here, I implement an unsupervised method using Self-Organizing (or "Kohonen") Maps (SOMs) for such purposes. Recent extensions called "SuperSOMs" can integrate multiple layers, each of which exerts independent influence on a two-dimensional output grid via empirically estimated weights. The grid cells are then delimited into K distinct units that can be interpreted as species or other entities. I show empirical examples in salamanders (Desmognathus) and snakes (Storeria) with layers representing alleles, space, climate, and traits. Simulations reveal that the SuperSOM approach can detect K = 1, tends not to over-split, reflects contributions from all layers, and limits large layers (e.g., genetic matrices) from overwhelming other datasets, desirable properties addressing major concerns from previous studies. Finally, I suggest that these and similar methods could integrate conservation-relevant layers such as population trends and human encroachment to delimit management units from an explicitly quantitative framework grounded in the ecology and evolution of species limits and boundaries.

Different patterns of introgression in a three species hybrid zone among European cave salamanders

Hybrid zones occur where genetically distinct populations meet, mate and produce offspring with mixed ancestry. In Plethodontid salamanders, introgressive hybridization is a common phenomenon, where hybrids backcross with parental populations leading to the spread of new alleles into the parental genomes. Whereas many hybrid zones have been reported in American Plethodontid salamanders, only a single hybrid zone has been documented in European plethodontids so far, which is located at the Apuan Alps in the Italian Peninsula. Here, we describe a previously unreported hybrid zone in the Northern Apennines involving all the three Plethodontid salamander species inhabiting the Italian Peninsula. We found 21 new Speleomantes sites of occurrence, from a hitherto unexplored area located at the boundaries between three Speleomantes species ranges. Using mitochondrial (Cytb and ND2 genes) and nuclear markers (two diagnostic SNPs at the NCX1 gene), we revealed a three-way contact zone where all the three mainland species hybridize: S. strinatii, S. ambrosii and S. italicus. We observed a strong mitonuclear discordance, with mitochondrial markers showing a conspicuous geographic pattern, while diagnostic nuclear SNPs coexisted in both the same populations and individuals, providing evidence of hybridization in many possible combinations. The introgression is asymmetric, with S. italicus mitogenome usually associated with S. a. ambrosii and, to a lesser extent, to S. strinatii nuclear alleles. This finding confirms that Plethodontid are a group of choice to investigate hybridization mechanisms and suggests that behavioural, genetic and ecological components may concur in determining the direction and extent of introgression.

Systematic revision of the Spotted and Northern Dusky Salamanders (Plethodontidae: Desmognathus conanti and D. fuscus), with six new species from the eastern United States

Spotted and Northern Dusky Salamanders (Desmognathus conanti and D. fuscus) have a long and complex taxonomic history. At least 10 other currently recognized species in the genus were either described from populations previously considered D. fuscus, described as or later considered subspecies thereof, or later considered synonyms thereof, before ultimately being recognized as distinct. Recent molecular analyses have also revealed extensive cryptic diversity within both species, which are polyphyletic assemblages of 13 distinct mitochondrial lineages with 5.7-10.3% uncorrected 'p' distances in the COI barcode locus. Based on phylogenomic data and population-clustering analyses considering admixture between lineages, 11 candidate species were circumscribed by recent authors. Those within D. conanti are also ecomorphologically variable, comprising both large, robust, keel-tailed populations, and small, gracile, round-tailed forms. Evaluating their distinctiveness based on genetic, geographic, and morphological evidence, we conclude that six of the candidates represent new species: Desmognathus anicetus sp. nov., D. bairdi sp. nov., D. campi sp. nov., D. catahoula sp. nov., D. lycos sp. nov., and D. tilleyi sp. nov. Consequently, we recognize eight total species from populations formerly associated with the nominal species D. conanti and D. fuscus, the re-delimited concepts of which also contain additional phylogeographic lineage diversity that may represent further distinct species. In addition to existing mitochondrial and nuclear phylogenetic, network, and clustering results, we present preliminary analyses of linear morphometrics to bolster diagnostic specificity based on phenotypic characteristics. These changes stabilize the previously paraphyletic taxonomy of species-level lineages within Desmognathus, though additional cryptic diversity may exist both within the species considered here, and elsewhere in the genus.

A systematic revision of the Shovel-nosed Salamander (Plethodontidae: Desmognathus marmoratus), with re-description of the related D. aureatus and D. intermedius

Shovel-nosed Salamanders, Desmognathus marmoratus (Moore, 1899), were long thought to represent a single species from the southern Appalachian Mountains of the eastern United States, ranging from northeastern Georgia to extreme southwestern Virginia. These populations have a highly derived ecomorphology, being fully aquatic with a specialized flattened and elongated phenotype adapted to rocky riffle zones in fast-flowing, high-gradient mountain streams. Because of this, they were originally described in a separate genus, Leurognathus Moore, 1899. Four additional species or subspecies were described from 1928-1956 based on regional geographic variation in phenotype before being synonymized with L. marmoratus in 1962, which was reassigned to Desmognathus in 1996. Molecular analyses subsequently revealed four distinct candidate lineages in two distantly related clades, which were recently re-delimited into three species. These are D. aureatus (Martof, 1956) from northeastern Georgia, D. intermedius (Pope, 1928) from western North Carolina, and D. marmoratus from northwestern North Carolina. We provide a systematic revision of these taxa, which do not represent a natural group but instead exhibit convergent phenotypes across multiple species, potentially driven by ancient episodes of adaptive introgression between ancestral lineages. Our recent fieldwork revealed an astonishingly disjunct and morphologically distinct population of D. marmoratus in the New River Gorge of West Virginia, which were previously confused with D. kanawha Pyron and Beamer, 2022. This locality is ~120 airline km away from the nearest populations of D. marmoratus in Virginia. No Shovel-nosed Salamanders have ever been found in the New River drainage during our extensive previous explorations or credibly reported in museum specimens or the literature. Additional cryptic populations of these taxa may remain.

Gene drive in species complexes: defining target organisms

Engineered gene drives, which bias their own inheritance to increase in frequency in target populations, are being developed to control mosquito malaria vectors. Such mosquitoes can belong to complexes of both vector and nonvector species that can produce fertile interspecific hybrids, making vertical gene drive transfer (VGDT) to sibling species biologically plausible. While VGDT to other vectors could positively impact human health protection goals, VGDT to nonvectors might challenge biodiversity ones. Therefore, environmental risk assessment of gene drive use in species complexes invites more nuanced considerations of target organisms and nontarget organisms than for transgenes not intended to increase in frequency in target populations. Incorporating the concept of target species complexes offers more flexibility when assessing potential impacts from VGDT.

A two-species distribution model for parapatric newts, with inferences on their history of spatial replacement

Related species often engage in abutting or overlapping contact zones with various strengths of interspecific competition. Biotic interactions such as these preclude the registration of the full profile of environmental variables that would describe the otherwise larger species ranges. Here, I advocate to forego full range species distribution modelling and instead focus on the ecography of the contact zone, for example with ‘two-species distribution models’ (TSDMs), in which presence data are contrasted against the background of environmental data. The newts Triturus cristatus and Triturus marmoratus meet in the west of France. A countrywide TSDM suggests that the contact zone of the species is located at a climatic gradient, in line with their north-eastern vs. south-western ranges. The species are also ecologically segregated by elevation and forestation, which is in line with a documented movement of the contact zone caused by hedgerow removal in lowland areas. Hindcasts for the Holocene suggest that the species contact zone was positioned at either the same place as at present or more to the south, depending on the amount of forestation. A forecast under climate warming predicts a fast movement to the north, but this scenario is deemed unrealistic. One reason is that recent habitat loss compromises dispersal and range expansion. Other species pairs to which TSDMs have been applied are listed for comparison.