Comparative species delimitation of a biological conservation icon

The United States Endangered Species Act (ESA) of 1973 set a precedent for biodiversity conservation across the globe.1 A key requirement of protections afforded by the ESA is the accurate delimitation of imperiled species. We present a comparative reference-based taxonomic approach to species delimitation that integrates genomic and morphological data for objectively assessing the distinctiveness of species targeted for protection by governmental agencies. We apply this protocol to the Snail Darter (Percina tanasi), a freshwater fish from the Tennessee River that was discovered in 1973 and declared an endangered species under the ESA in 1975.2 Concurrently, the Snail Darter's habitat was slated to be destroyed through the construction of the Tellico Dam by the Tennessee Valley Authority (TVA),3,4 inspiring nationwide protests advocating for the suspension of the federal project. This David versus Goliath struggle between supporters of the 3-inch fish and the TVA culminated in the first major legal conflict over protections afforded by the ESA, the US Supreme Court case Hill v. TVA, 437 U.S. 153 (1978), with a 6 to 3 ruling in favor of protecting the Snail Darter and interrupting the completion of the Tellico Dam. Here, we integrate multiple lines of evidence in a comparative framework to demonstrate that despite its legacy, the Snail Darter is not a distinct species but is a population of the Stargazing Darter (Percina uranidea) described in 1887. These results illustrate how a reference-based framework for species delimitation dramatically aids the proper direction of efforts toward protecting biodiversity.

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.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Diverge and Conquer: Phylogenomics of southern Wallacean forest skinks (Genus: Sphenomorphus) and their colonization of the Lesser Sunda Archipelago

The archipelagos of Wallacea extend between the Sunda and Sahul Shelves, serving as a semipermeable two-way filter influencing faunal exchange between Asia and Australo-Papua. Forest skinks (Genus Sphenomorphus) are widespread throughout southern Wallacea and exhibit complex clinal, ontogenetic, sexual, and seasonal morphological variation, rendering species delimitation difficult. We screened a mitochondrial marker for 245 Sphenomorphus specimens from this area to inform the selection of 104 samples from which we used targeted sequence capture to generate a dataset of 1154 nuclear genes (∼1.8 Mb) plus complete mitochondrial genomes. Phylogenomic analyses recovered many deeply divergent lineages, three pairs of which are now sympatric, that began to diversify in the late Miocene shortly after the oldest islands are thought to have become emergent. We infer a complex and nonstepping-stone pattern of island colonization, with the group having originated in the Sunda Arc islands before using Sumba as a springboard for colonization of the Banda Arcs. Estimates of population structure and gene flow across the region suggest total isolation except between two Pleistocene Aggregate Island Complexes that become episodically land-bridged during glacial maxima. These historical processes have resulted in at least 11 Sphenomorphus species in the region, nine of which require formal description. This fine-scale geographic partitioning of undescribed species highlights the importance of utilizing comprehensive genomic studies for defining biodiversity hotspots to be considered for conservation protection.

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

With limited sampling, geographic variation within a single species can be difficult to distinguish from interspecific variation, confounding our ability to draw accurate species boundaries. We argue that thorough sampling and analysis of contact zones between putative taxa can determine if assortative mating or selection against hybrids exists (supporting the presence of two distinct species), or alternatively if mating is random among genotypes and admixture among adjacent populations is gradual and continuous (supporting geographic variation within a single species). Here, we test two alternative hypotheses for two pairs of named taxa at contact zones within the American milksnake (Lampropeltis triangulum) complex. A prior morphological analysis found areas of gradual intergradation among named taxa, and concluded that the taxa represented geographical races of a single polytypic species. In contrast, a subsequent analysis of gene sequence data, but with limited sampling near the contact zones, hypothesized distinct boundaries between species at the contact zones. At the contact zone between proposed species L. triangulum and L. gentilis, we examined a ∼700 km-wide transect across the states of Kansas and Missouri, with thorough sampling and reduced-representation genomic-level sequencing, to test the two opposing taxonomic hypotheses. Our transect analyses included examinations of population structure, fixed differences, cline-fitting, and an admixture index analysis. These analyses all supported a gradual and continuous geographic cline across a broad intergrade zone between two geographic forms of L. triangulum, thus providing strong support for a single species in this region (and no support for the recognition of L. gentilis as a distinct species). At a second contact zone between proposed species L. triangulum and L. elapsoides (but variously treated as species or subspecies by different researchers) in Kentucky and Tennessee, we re-evaluated morphological data. In this case, the contact zone analysis indicated sympatry and reproductive isolation of the two taxa, and thus strongly supported L. triangulum and L. elapsoides as distinct species. We conclude that detailed studies of contact zones, based on either genetic or morphological data, are essential for distinguishing intraspecific from interspecific variation in the case of widely and continuously distributed taxa.

Using natural history to guide supervised machine learning for cryptic species delimitation with genetic data

The diversity of biological and ecological characteristics of organisms, and the underlying genetic patterns and processes of speciation, makes the development of universally applicable genetic species delimitation methods challenging. Many approaches, like those incorporating the multispecies coalescent, sometimes delimit populations and overestimate species numbers. This issue is exacerbated in taxa with inherently high population structure due to low dispersal ability, and in cryptic species resulting from nonecological speciation. These taxa present a conundrum when delimiting species: analyses rely heavily, if not entirely, on genetic data which over split species, while other lines of evidence lump. We showcase this conundrum in the harvester Theromaster brunneus, a low dispersal taxon with a wide geographic distribution and high potential for cryptic species. Integrating morphology, mitochondrial, and sub-genomic (double-digest RADSeq and ultraconserved elements) data, we find high discordance across analyses and data types in the number of inferred species, with further evidence that multispecies coalescent approaches over split. We demonstrate the power of a supervised machine learning approach in effectively delimiting cryptic species by creating a "custom" training data set derived from a well-studied lineage with similar biological characteristics as Theromaster. This novel approach uses known taxa with particular biological characteristics to inform unknown taxa with similar characteristics, using modern computational tools ideally suited for species delimitation. The approach also considers the natural history of organisms to make more biologically informed species delimitation decisions, and in principle is broadly applicable for taxa across the tree of life.