Phylogenomic data reveal reticulation and incongruence among mitochondrial candidate species in Dusky Salamanders (Desmognathus)

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.].

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.

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.

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.

The emergence of a cryptic lineage and cytonuclear discordance through past hybridization in the Japanese fire-bellied newt, Cynops pyrrhogaster (Amphibia: Urodela)

Discrepancies in geographic variation patterns between nuclear DNA and mitochondrial DNA (mtDNA) are the result of the complicated differentiation processes in organisms and are the key to understanding their true evolutionary processes. The genetic differentiation of the northern and Southern-Izu lineages of the Japanese newt, Cynops pyrrhogaster, was investigated through their single nucleotide polymorphism variations obtained via multiplexed ISSR genotyping by sequencing (MIG-seq). We found three genetic groups (Tohoku, N-Kanto and S-Kanto), that were not detected by mtDNA variations, in the northern lineage. N-Kanto has intermediate genetic characteristics between Tohoku and S-Kanto. The genetic groups are now moderately isolated from each other and have unique genetic characteristics. An estimation of the evolutionary history using the approximate Bayesian computation (ABC) approach suggested that Tohoku diverged from the common ancestor of S-Kanto and S-Izu. Then, S-Kanto and S-Izu split, and the recent hybridization between Tohoku and S-Kanto gave rise to N-Kanto. The origin of N-Kanto through the hybridization is relatively young and seems to be related to changes in the distributions of Tohoku and S-Kanto as a result of climatic oscillation in the Pleistocene. We conclude that the mitochondrial genome of S-Kanto was captured in Tohoku and that the original mitochondrial genome of Tohoku was entirely removed through hybridization.

The Syngameon Enigma

Despite their evolutionary relevance, multispecies networks or syngameons are rarely reported in the literature. Discovering how syngameons form and how they are maintained can give insight into processes such as adaptive radiations, island colonizations, and the creation of new hybrid lineages. Understanding these complex hybridization networks is even more pressing with anthropogenic climate change, as syngameons may have unique synergistic properties that will allow participating species to persist. The formation of a syngameon is not insurmountable, as several ways for a syngameon to form have been proposed, depending mostly on the magnitude and frequency of gene flow events, as well as the relatedness of its participants. Episodic hybridization with small amounts of introgression may keep syngameons stable and protect their participants from any detrimental effects of gene flow. As genomic sequencing becomes cheaper and more species are included in studies, the number of known syngameons is expected to increase. Syngameons must be considered in conservation efforts as the extinction of one participating species may have detrimental effects on the survival of all other species in the network.

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

Dusky Salamanders (genus Desmognathus) currently comprise only 22 described, extant species. However, recent mitochondrial and nuclear estimates indicate the presence of up to 49 candidate species based on ecogeographic sampling. Previous studies also suggest a complex history of hybridization between these lineages. Studies in other groups suggest that disregarding admixture may affect both phylogenetic inference and clustering-based species delimitation. With a dataset comprising 233 Anchored Hybrid Enrichment (AHE) loci sequenced for 896 Desmognathus specimens from all 49 candidate species, we test three hypotheses regarding (i) species-level diversity, (ii) hybridization and admixture, and (iii) misleading phylogenetic inference. Using phylogenetic and population-clustering analyses considering gene flow, we find support for at least 47 candidate species in the phylogenomic dataset, some of which are newly characterized here while others represent combinations of previously named lineages that are collapsed in the current dataset. Within these, we observe significant phylogeographic structure, with up to 64 total geographic genetic lineages, many of which hybridize either narrowly at contact zones or extensively across ecological gradients. We find strong support for both recent admixture between terminal lineages and ancient hybridization across internal branches. This signal appears to distort concatenated phylogenetic inference, wherein more heavily admixed terminal specimens occupy apparently artifactual early-diverging topological positions, occasionally to the extent of forming false clades of intermediate hybrids. Additional geographic and genetic sampling and more robust computational approaches will be needed to clarify taxonomy, and to reconstruct a network topology to display evolutionary relationships in a manner that is consistent with their complex history of reticulation.

Gene flow in phylogenomics: Sequence capture resolves species limits and biogeography of Afromontane forest endemic frogs from the Cameroon Highlands

Puddle frogs of the Phrynobatrachus steindachneri species complex are a useful group for investigating speciation and phylogeography in Afromontane forests of the Cameroon Volcanic Line, western Central Africa. The species complex is represented by six morphologically relatively cryptic mitochondrial DNA lineages, only two of which are distinguished at the species level - southern P. jimzimkusi and Lake Oku endemic P. njiomock, leaving the remaining four lineages identified as 'P. steindachneri'. In this study, the six mtDNA lineages are subjected to genomic sequence capture analyses and morphological examination to delimit species and to study biogeography. The nuclear DNA data (387 loci; 571,936 aligned base pairs) distinguished all six mtDNA lineages, but the topological pattern and divergence depths supported only four main clades: P. jimzimkusi, P. njiomock, and only two divergent evolutionary lineages within the four 'P. steindachneri' mtDNA lineages. One of the two lineages is herein described as a new species, P. amieti sp. nov. Reticulate evolution (hybridization) was detected within the species complex with morphologically intermediate hybrid individuals placed between the parental species in phylogenomic analyses, forming a ladder-like phylogenetic pattern. The presence of hybrids is undesirable in standard phylogenetic analyses but is essential and beneficial in the network multispecies coalescent. This latter approach provided insight into the reticulate evolutionary history of these endemic frogs. Introgressions likely occurred during the Middle and Late Pleistocene climatic oscillations, due to the cyclic connections (likely dominating during cold glacials) and separations (during warm interglacials) of montane forests. The genomic phylogeographic pattern supports the separation of the southern (Mt. Manengouba to Mt. Oku) and northern mountains at the onset of the Pleistocene. Further subdivisions occurred in the Early Pleistocene, separating populations from the northernmost (Tchabal Mbabo, Gotel Mts.) and middle mountains (Mt. Mbam, Mt. Oku, Mambilla Plateau), as well as the microendemic lineage restricted to Lake Oku (Mt. Oku). This unique model system is highly threatened as all the species within the complex have exhibited severe population declines in the past decade, placing them on the brink of extinction. In addition, Mount Oku is identified to be of particular conservation importance because it harbors three species of this complex. We, therefore, urge for conservation actions in the Cameroon Highlands to preserve their diversity before it is too late.

UCE Phylogenomics, detection of a putative hybrid population, and one older mitogenomic node age of Batrachuperus salamanders

The prevalence of incomplete lineage sorting complicates the examination of hybridization and species-level paraphyly with gene trees of a small number of loci. In Asian mountain salamanders of the genus Batrachuperus, possible hybridization and species paraphyly had been identified by utilizing mitochondrial genealogy and fixed allozyme differences. Here we sampled 2909 UCEs in 44 local populations from all six Batrachuperus species, inferred gene and species trees, compared them with mitochondrial and allozyme results, and examined the potential hybridization and species paraphyly. The clustering pattern of single-locus trees, increased proportion of heterozygous SNPs, allele frequency-based migration edge estimation, and intrapopulation long branches (as expected from an increase of genetic lineage and nucleotide diversity) support that an eastern B. karlschmidti population has experienced admixture with B. tibetanus. On the 2909-UCE concatenated and species trees, lower nodal supports were observed when similar proportions of loci agreed with alternative topologies, i.e., a reciprocal monophyly between a Pengxian lineage and the remainder of B. pinchonii (0.379) or a paraphyly of the latter with respect to Pengxian (0.362). The UCE phylogenomics agreed with the relatively recent groupings in the allozyme dendrogram. Despite incomplete lineage sorting, the mitochondrial trees were similar to the UCE trees for deeper relationships of the genus. However, one significant branch-length level discordance was identified. The branch between the common ancestor of B. daochengensis and B. yenyuanensis and common ancestor of the genus was approximately three times shorter on the mitochondrial tree than on the UCE tree, suggesting that the split of the mitochondrial lineages was likely a few million years earlier than the split of species. This finding supports considering possible ancestral polymorphism when interpreting different divergence dates estimated from mitochondrial and genome-wide data.

Preparing for invasion: Assessing risk of infection by chytrid fungi in southeastern plethodontid salamanders

Understanding the responses of naïve communities to the invasion of multihost pathogens requires accurate estimates of susceptibility across taxa. In the Americas, the likely emergence of a second amphibian pathogenic fungus (Batrachochytrium salamandrivorans, Bsal) calls for new ways of prioritizing disease mitigation among species due to the high diversity of naïve hosts with prior B. dendrobatidis (Bd) infections. Here, we applied the concept of pathogenic potential to quantify the virulence of chytrid fungi on naïve amphibians and evaluate species for conservation efforts in the event of an outbreak. The benefit of this measure is that it combines and summarizes the variation in disease effects into a single numerical index, allowing for comparisons across species, populations or groups of individuals that may inherently exhibit differences in susceptibility. As a proof of concept, we obtained standardized responses of disease severity by performing experimental infections with Bsal on five plethodontid salamanders from southeastern United States. Four out of five species carried natural infections of Bd at the start of the experiments. We showed that Bsal exhibited its highest value of pathogenic potential in a species that is already declining (Desmognathus auriculatus). We find that this index provides additional information beyond the standard measures of disease prevalence, intensity, and mortality, because it leveraged these disease parameters within each categorical group. Scientists and practitioners could use this measure to justify research, funding, trade, or conservation measures.

Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal

Phylogenomic data from a rapidly increasing number of studies provide new evidence for resolving relationships in recently radiated clades, but they also pose new challenges for inferring evolutionary histories. Most existing methods for reconstructing phylogenetic hypotheses rely solely on algorithms that only consider incomplete lineage sorting (ILS) as a cause of intra- or intergenomic discordance. Here, we utilize a variety of methods, including those to infer phylogenetic networks, to account for both ILS and introgression as a cause for nuclear and cytoplasmic-nuclear discordance using phylogenomic data from the recently radiated flowering plant genus Polemonium (Polemoniaceae), an ecologically diverse genus in Western North America with known and suspected gene flow between species. We find evidence for widespread discordance among nuclear loci that can be explained by both ILS and reticulate evolution in the evolutionary history of Polemonium. Furthermore, the histories of organellar genomes show strong discordance with the inferred species tree from the nuclear genome. Discordance between the nuclear and plastid genome is not completely explained by ILS, and only one case of discordance is explained by detected introgression events. Our results suggest that multiple processes have been involved in the evolutionary history of Polemonium and that the plastid genome does not accurately reflect species relationships. We discuss several potential causes for this cytoplasmic-nuclear discordance, which emerging evidence suggests is more widespread across the Tree of Life than previously thought. [Cyto-nuclear discordance, genomic discordance, phylogenetic networks, plastid capture, Polemoniaceae, Polemonium, reticulations.].

Developmental life history is associated with variation in rates of climatic niche evolution in a salamander adaptive radiation

Rates of climatic niche evolution vary widely across the tree of life and are strongly associated with rates of diversification among clades. However, why the climatic niche evolves more rapidly in some clades than others remains unclear. Variation in life history traits often plays a key role in determining the environmental conditions under which species can survive, and therefore, could impact the rate at which lineages can expand in available climatic niche space. Here, we explore the relationships among life-history variation, climatic niche breadth, and rates of climatic niche evolution. We reconstruct a phylogeny for the genus Desmognathus, an adaptive radiation of salamanders distributed across eastern North America, based on nuclear and mitochondrial genes. Using this phylogeny, we estimate rates of climatic niche evolution for species with long, short, and no aquatic larval stage. Rates of climatic niche evolution are unrelated to the mean climatic niche breadth of species with different life histories. Instead, we find that the evolution of a short larval period promotes greater exploration of climatic space, leading to increased rates of climatic niche evolution across species having this trait. We propose that morphological and physiological differences associated with variation in larval stage length underlie the heterogeneous ability of lineages to explore climatic niche space. Rapid rates of climatic niche evolution among species with short larval periods were an important dimension of the clade's adaptive radiation and likely contributed to the rapid rate of lineage accumulation following the evolution of an aquatic life history in this clade. Our results show how variation in a key life-history trait can constrain or promote divergence of the climatic niche, leading to variation in rates of climatic niche evolution among species.