Species delimitation in endangered groundwater salamanders: Implications for aquifer management and biodiversity conservation

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Groundwater is a hidden global keystone ecosystem

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Regional climate contributes more than geographic distance to beta diversity of copepods (Crustacea Copepoda) between caves of Italy

Despite the study of subterranean biodiversity facing harsh sampling and mapping challenges, the huge diversity of taxa, ecological adaptations and evolutionary trajectories in subterranean environments is gaining increasing attention. Yet, the spatial and environmental factors driving the composition of groundwater communities are still poorly understood. To partially fill this knowledge gap, we collected copepod crustaceans from 12 caves along the Italian peninsula between 2019 and 2022, sampling each cave twice. The resulting presence-absence data were analysed to assess: (i) between-cave taxonomic beta diversity, also partitioning between turnover and nestedness-resultant dissimilarity; (ii) the relative weight of geographic distance and climatic differences in shaping observed beta diversity. Seventy-one species of copepods were collected overall. Pairwise beta diversity was high for most pairs of caves, with turnover being the major component. Geographic distance-decay models partially explained total beta diversity and turnover patterns. However, in Generalized Dissimilarity Models (GDM), including surface climatic conditions as predictors, the contribution of seasonal temperature averages was generally higher than that of geographic distance. Further, the explanatory and predictive performance of the GDMs notably increased, along with temperature contribution, when widening the spatial extent from which climate data were gathered. Our results confirmed a high spatial turnover in groundwater copepods' assemblages and strengthened the link between regional climate and subterranean biodiversity.

Development and validation of an environmental DNA assay to detect federally threatened groundwater salamanders in central Texas

The molecular detection of DNA fragments that are shed into the environment (eDNA) has become an increasingly applied tool used to inventory biological communities and to perform targeted species surveys. This method is particularly useful in habitats where it is difficult or not practical to visually detect or trap the target organisms. Central Texas Eurycea salamanders inhabit both surface and subterranean aquatic environments. Subterranean surveys are challenging or infeasible, and the detection of salamander eDNA in water samples is an appealing survey technique for these situations. Here, we develop and validate an eDNA assay using quantitative PCR for E. chisholmensis, E. naufragia, and E. tonkawae. These three species are federally threatened and constitute the Septentriomolge clade that occurs in the northern segment of the Edwards Aquifer. First, we validated the specificity of the assay in silico and with DNA extracted from tissue samples of both target Septentriomolge and non-target amphibians that overlap in distribution. Then, we evaluated the sensitivity of the assay in two controls, one with salamander-positive water and one at field sites known to be occupied by Septentriomolge. For the salamander-positive control, the estimated probability of eDNA occurrence (ψ) was 0.981 (SE = 0.019), and the estimated probability of detecting eDNA in a qPCR replicate (p) was 0.981 (SE = 0.011). For the field control, the estimated probability of eDNA occurring at a site (ψ) was 0.938 (95% CRI: 0.714-0.998). The estimated probability of collecting eDNA in a water sample (θ) was positively correlated with salamander relative density and ranged from 0.371 (95% CRI: 0.201-0.561) to 0.999 (95% CRI: 0.850- > 0.999) among sampled sites. Therefore, sites with low salamander density require more water samples for eDNA evaluation, and we determined that our site with the lowest estimated θ would require seven water samples for the cumulative collection probability to exceed 0.95. The estimated probability of detecting eDNA in a qPCR replicate (p) was 0.882 (95% CRI: 0.807-0.936), and our assay required two qPCR replicates for the cumulative detection probability to exceed 0.95. In complementary visual encounter surveys, the estimated probability of salamanders occurring at a known-occupied site was 0.905 (SE = 0.096), and the estimated probability of detecting salamanders in a visual encounter survey was 0.925 (SE = 0.052). We additionally discuss future research needed to refine this method and understand its limitations before practical application and incorporation into formal survey protocols for these taxa.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

Gene flow assessment helps to distinguish strong genomic structure from speciation in an Iberian ant-eating spider

Although genomic data is boosting our understanding of evolution, we still lack a solid framework to perform reliable genome-based species delineation. This problem is especially critical in the case of phylogeographically structured organisms, with allopatric populations showing similar divergence patterns as species. Here, we assess the species limits and phylogeography of Zodarion alacre, an ant-eating spider widely distributed across the Iberian Peninsula. We first performed species delimitation based on genome-wide data and then validated these results using additional evidence. A commonly employed species delimitation strategy detected four distinct lineages with almost no admixture, which present allopatric distributions. These lineages showed ecological differentiation but no clear morphological differentiation, and evidence of introgression in a mitochondrial barcode. Phylogenomic networks found evidence of substantial gene flow between lineages. Finally, phylogeographic methods highlighted remarkable isolation by distance and detected evidence of range expansion from south-central Portugal to central-north Spain. We conclude that despite their deep genomic differentiation, the lineages of Z. alacre do not show evidence of complete speciation. Our results likely shed light on why Zodarion is among the most diversified spider genera despite its limited distribution and support the use of gene flow evidence to inform species boundaries.

Evolution in the dark: Unexpected genetic diversity and morphological stasis in the blind, aquifer-dwelling catfish Horaglanis

The lateritic aquifers of the southern Indian state of Kerala harbour a unique assemblage of enigmatic stygobitic fishes which are encountered very rarely, only when they surface during the digging and cleaning of homestead wells. Here, we focus on one of the most unusual members of this group, the catfish Horaglanis, a genus of rarely-collected, tiny, blind, pigment less, and strictly aquifer-residing species. A six-year exploratory and citizen-science backed survey supported by molecular phylogenetic analysis reveals novel insights into the diversity, distribution and population structure of Horaglanis. The genus is characterized by high levels of intraspecific and interspecific genetic divergence, with phylogenetically distinct species recovered above a 7.0% genetic-distance threshold in the mitochondrial cytochrome oxidase subunit 1 gene. Contrasting with this deep genetic divergence, however, is a remarkable stasis in external morphology. We identify and describe a new cryptic species, Horaglanis populi, a lineage that is the sister group of all currently known species. All four species are represented by multiple haplotypes. Mismatch distribution reveals that populations have not experienced recent expansions.

Fine-scale monitoring and mapping of biodiversity and ecosystem services reveals multiple synergies and few tradeoffs in urban green space management

Urban watersheds can play a critical role in supporting biodiversity and ecosystem services in a rapidly changing world. However, managing for multiple environmental and social objectives in urban landscapes is challenging, especially if the optimization of one ecosystem service conflicts with another. Urban ecology research has frequently been limited to a few indicators - typically either biodiversity or ecosystem service indices - making tradeoffs and synergies difficult to assess. Through a recently established watershed-scale monitoring network in Central Texas, we address this gap by evaluating biodiversity (flora and fauna), habitat quality, and ecosystem service indices of urban green spaces across the watershed. Our results reveal substantial heterogeneity in biodiversity and ecosystem service levels and multiple synergies (stacked benefits or "win-wins"). For example, we found that carbon sequestration positively correlated with tree species richness and the proportion of native trees in a green space, indicating that biodiversity goals for increased tree diversity can also provide carbon sequestration benefits. We also documented correlations between green spaces with greater riparian forest cover and lower particulate matter (PM_2.5) concentrations and cooler temperatures. In addition, we found that bee and wasp species richness was positively correlated with carbon sequestration and human visitation rates, meaning that urban green spaces can optimize carbon sequestration goals without losing pollinator habitat or access opportunities for city residents. Overall, our results indicate that many aspects of habitat quality, biodiversity, and ecosystem services can be simultaneously supported in urban green spaces. We conclude that urban design and management can optimize nature-based solutions and strategies to have distinct positive impacts on both people and nature.

Species boundary and phylogeographical pattern provide new insights into the management efforts of Megacopta cribraria (Hemiptera: Plataspidae), a bean bug invading North America

BACKGROUND

Correct identification together with information on distribution range, geographical origin and evolutionary history are the necessary basis for the management and control of invasive species. The bean bug Megacopta cribraria is a crucial agricultural pest of soybean. Recently, M. cribraria has invaded the United States and spread rapidly, causing severe reductions in soybean yields. However, the species boundary and phylogeographical pattern of this invasive bean bug are still unclear.

RESULTS

The results of different species delimitation methods (Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Bayesian Poisson Tree Processes and Bayesian Phylogenetics and Phylogeography) strongly demonstrated that M. cribraria and Megacopta punctatissima represent the same species. M. punctatissima should not be considered a distinct species but rather a variety of M. cribraria. Phylogenetic analyses revealed three well-supported clades (Southeast Asia [SEA], East Asia continent [EAC] and Japan [JA]) with distinct geographical structures in the M. cribraria-M. punctatissima complex. The SEA clade was at the base of the phylogenetic tree, and the sister relationship between the EAC clade and JA clade was strongly supported. The split between the EAC clade and JA clade occurred at approximately 0.71 Ma, corresponding to the submergence period of the East China Sea land bridge.

CONCLUSION

This study clarified the species boundary between M. cribraria and its closely related species and revealed the phylogeographical pattern and evolutionary history of M. cribraria. The species delimitation and phylogeography results achieved in this study could provide new insights into the monitoring and management of this agricultural pest. © 2022 Society of Chemical Industry.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

PREMISE

Understanding evolutionary history and classifying discrete units of organisms remain overwhelming tasks, and lags in this workload concomitantly impede an accurate documentation of biodiversity and conservation management. Rapid advances and improved accessibility of sensitive high-throughput sequencing tools are fortunately quickening the resolution of morphological complexes and   thereby improving the estimation of species diversity. The recently described and critically endangered Banksia vincentia is morphologically similar to the hairpin banksia complex (B. spinulosa s.l.), a group of eastern Australian flowering shrubs whose continuum of morphological diversity has been responsible for taxonomic controversy and possibly questionable conservation initiatives.

METHODS

To assist conservation while testing the current taxonomy of this group, we used high-throughput sequencing to infer a population-scale evolutionary scenario for a sample set that is comprehensive in its representation of morphological diversity and a 2500-km distribution.

RESULTS

Banksia spinulosa s.l. represents two clades, each with an internal genetic structure shaped through historical separation by biogeographic barriers. This structure conflicts with the existing taxonomy for the group. Corroboration between phylogeny and population statistics aligns with the hypothesis that B. collina, B. neoanglica, and B. vincentia should not be classified as species.

CONCLUSIONS

The pattern here supports how morphological diversity can be indicative of a locally expressed suite of traits rather than relationship. Oversplitting in the hairpin banksias is atypical since genomic analyses often reveal that species diversity is underestimated. However, we show that erring on overestimation can yield negative consequences, such as the disproportionate prioritization of a geographically anomalous population.

Towards evidence-based conservation of subterranean ecosystems

Subterranean ecosystems are among the most widespread environments on Earth, yet we still have poor knowledge of their biodiversity. To raise awareness of subterranean ecosystems, the essential services they provide, and their unique conservation challenges, 2021 and 2022 were designated International Years of Caves and Karst. As these ecosystems have traditionally been overlooked in global conservation agendas and multilateral agreements, a quantitative assessment of solution-based approaches to safeguard subterranean biota and associated habitats is timely. This assessment allows researchers and practitioners to understand the progress made and research needs in subterranean ecology and management. We conducted a systematic review of peer-reviewed and grey literature focused on subterranean ecosystems globally (terrestrial, freshwater, and saltwater systems), to quantify the available evidence-base for the effectiveness of conservation interventions. We selected 708 publications from the years 1964 to 2021 that discussed, recommended, or implemented 1,954 conservation interventions in subterranean ecosystems. We noted a steep increase in the number of studies from the 2000s while, surprisingly, the proportion of studies quantifying the impact of conservation interventions has steadily and significantly decreased in recent years. The effectiveness of 31% of conservation interventions has been tested statistically. We further highlight that 64% of the reported research occurred in the Palearctic and Nearctic biogeographic regions. Assessments of the effectiveness of conservation interventions were heavily biased towards indirect measures (monitoring and risk assessment), a limited sample of organisms (mostly arthropods and bats), and more accessible systems (terrestrial caves). Our results indicate that most conservation science in the field of subterranean biology does not apply a rigorous quantitative approach, resulting in sparse evidence for the effectiveness of interventions. This raises the important question of how to make conservation efforts more feasible to implement, cost-effective, and long-lasting. Although there is no single remedy, we propose a suite of potential solutions to focus our efforts better towards increasing statistical testing and stress the importance of standardising study reporting to facilitate meta-analytical exercises. We also provide a database summarising the available literature, which will help to build quantitative knowledge about interventions likely to yield the greatest impacts depending upon the subterranean species and habitats of interest. We view this as a starting point to shift away from the widespread tendency of recommending conservation interventions based on anecdotal and expert-based information rather than scientific evidence, without quantitatively testing their effectiveness.

Conservation genomics of urban populations of Streamside Salamander (Ambystoma barbouri)

In Tennessee, populations of the state endangered Streamside Salamander (Ambystoma barbouri) are in decline as their distribution lies mostly within rapidly developing areas in the Nashville Basin. Information regarding the partitioning of genetic variation among populations of A. barbouri and the taxonomic status of these populations relative to northern populations and their congener, the Small-mouthed Salamander (A. texanum), have important implications for management and conservation of this species. Here we combined mitochondrial sequencing and genome-wide single nucleotide polymorphism (SNP) data generated using Genotyping-by-Sequencing (GBS) to investigate patterns of genetic variation within Tennessee populations of A. barbouri, to assess their relationship to populations in Kentucky, and to examine their phylogenetic relationship to the closely related A. texanum. Results from phylogenetic reconstructions reveal a complex history of Tennessee A. barbouri populations with regards to northern populations, unisexual A. barbouri, and A. texanum. Patterns of mitochondrial sequence variation suggest that A. barbouri may have originated within Tennessee and expanded north multiple times into Kentucky, Ohio, Indiana, and West Virginia. Phylogenetic reconstructions based on genome-wide SNP data contradict results based on mitochondrial DNA and correspond to geographic and taxonomic boundaries. Variation in allele frequencies at SNP genotypes, as identified by multivariate analyses and Bayesian assignment tests, identified three evolutionary significant units (ESUs) for A. barbouri within Tennessee. Collectively, these results emphasize the need for prioritizing conservation needs for Tennessee populations of A. barbouri to ensure the long-term persistence of this species.

Landscape connectivity among coastal giant salamander (Dicamptodon tenebrosus) populations shows no association with land use, fire frequency, or river drainage but exhibits genetic signatures of potential conservation concern

Determining the genetic consequences of both historical and contemporary events can clarify the effects of the environment on population connectivity and inform conservation decisions. Historical events (like glaciations) and contemporary factors (like logging) can disrupt gene flow between populations. This is especially true among species with specialized ecological requirements and low dispersal ability, like amphibians. We test for the genetic consequences of historical and contemporary disturbances in the coastal giant salamander (Dicamptodon tenebrosus) in the Pacific Northwest of the United States. We consider predictions based on the contemporary landscape (habitat connectivity, logging, forest fires, and topography), in addition to relatively ancient post-Pleistocene range expansion (following the last glacial retreat). To assess local versus larger-scale effects, we sampled 318 individuals across 23 sites, which were clustered in five sampling regions. Genetic variation was assessed using five microsatellite markers. We found evidence of (i) historical regional isolation, with decreased genetic diversity among more recently colonized northern sites, as well as (ii) high levels of inbreeding and loss of heterozygosity at local scales, despite relatively low overall population differentiation (FST) or strong evidence for population bottlenecks. Genetic diversity was not associated with contemporary disturbances (logging or fire), and there were no detectable effects on the genetic connectivity of populations based on intervening landscape features (habitat fragmentation and topography). However, lower genetic diversity in more northern regions indicates a lag in recovery of genetic diversity following post-Pleistocene expansion. Additionally, some populations had evidence of having undergone a recent genetic bottleneck or had high inbreeding (FIS) values. Lower genetic diversity in more northern sites means populations may be more vulnerable to future environmental changes, and managing for connectivity alone may not be sufficient given low mobility. Recent apparent reductions in some populations were not clearly linked to anthropogenic disturbances we examined. This suggests the type of disturbances this species is sensitive to may not be well understood.

Phylogenomics and species delimitation of the economically important Black Basses (Micropterus)

Informed management and conservation efforts are vital to sustainable recreational fishing and biodiversity conservation. Because the taxonomic rank of species is important in conservation and management strategies, success of these efforts depends on accurate species delimitation. The Black Basses (Micropterus) are an iconic lineage of freshwater fishes that include some of the world’s most popular species for recreational fishing and world's most invasive species. Despite their popularity, previous studies to delimit species and lineages in Micropterus suffer from insufficient geographic coverage and uninformative molecular markers. Our phylogenomic analyses of ddRAD data result in the delimitation of 19 species of Micropterus, which includes 14 described species, the undescribed but well-known Altamaha, Bartram’s, and Choctaw basses, and two additional undescribed species currently classified as Smallmouth Bass (M. dolomieu). We provide a revised delimitation of species in the Largemouth Bass complex that necessitates a change in scientific nomenclature: Micropterus salmoides is retained for the Florida Bass and Micropterus nigricans is elevated from synonymy for the Largemouth Bass. The new understanding of diversity, distribution, and systematics of Black Basses will serve as important basis for the management and conservation of this charismatic and economically important clade of fishes.

Predicting surface abundance of federally threatened Jollyville Plateau Salamanders (Eurycea tonkawae) to inform management activities at a highly modified urban spring

Urban expansion has contributed to the loss of habitat for range restricted species across the globe. Managing wildlife populations within these urban settings presents the challenge of balancing human and wildlife needs. Jollyville Plateau Salamanders (Eurycea tonkawae) are a range restricted, federally threatened, species of neotenic brook salamander endemic to central Texas. Almost the entire geographic range of E. tonkawae is embedded in the Austin, Cedar Park, and Round Rock metropolitan areas of Travis and Williamson counties, Texas. Among E. tonkawae occupied sites, Brushy Creek Spring has experienced some of the most extensive anthropogenic disturbance. Today the site consists of small groundwater outlets that emerge in the seams within a concrete culvert underlying a highway. Salamanders persist within this system though they are rarely detected. Here, we model the occurrence of salamanders within the surface habitat of Brushy Creek Spring using generalized linear models. In the absence of available data regarding the amount of water that is discharged from the spring, we use accumulated rainfall as a proxy for discharge to estimate salamander abundance. Additionally, we present evidence of reproduction, recruitment, and subterranean movement by E. tonkawae throughout this site. Infrastructure maintenance is inevitable at Brushy Creek Spring. We intend for our results to inform when maintenance should occur, i.e., during environmental conditions when salamanders are less likely to be observed in the surface habitat, to avoid unnecessary impacts to this federally threatened species.

An integrative re-evaluation of Typhlatya shrimp within the karst aquifer of the Yucatán Peninsula, Mexico

The Yucatán Peninsula, Mexico is a carbonate platform well-known for extensive karst networks of densely stratified aquifer ecosystems. This aquifer supports diverse anchialine fauna, including species of the globally distributed anchialine shrimp genus Typhlatya (Atyidae). Four species (T. campecheae, T. pearsei, T. dzilamensis and T. mitchelli) are endemic to the Peninsula, of which three are federally listed in Mexico. This first integrative evaluation (i.e., molecular, morphological, broad geographic and type locality sampling, and environmental data) of Yucatán Typhlatya reveals considerable species identity conflict in prior phylogenetic assessments, broad species ranges, syntopy within cave systems and five genetic lineages (of which two are new to science). Despite sampling from the type locality of endangered T. campecheae, specimens (and molecular data) were indistinguishable from vulnerable T. pearsei. Ancestral/divergence reconstructions support convergent evolution of a low-salinity ancestor for a post-Paleogene arc Yucatán + Cuba Typhlatya clade within the anchialine Atyidae clade. A secondary adaptation for the coastal-restricted euryhaline (2–37 psu), Typhlatya dzilamensis (unknown conservation status) was identified, while remaining species lineages were low-salinity (< 5 psu) adapted and found within the meteoric lens of inland and coastal caves. This study demonstrates the need for integrative/interdisciplinary approaches when conducting biodiversity assessments in complex and poorly studied aquifers.

SPEDE‐sampler: an R Shiny application to assess how methodological choices and taxon‐sampling can affect Generalised Mixed Yule Coalescent (GMYC) output and interpretation

Species delimitation tools are vital to taxonomy and the discovery of new species. These tools can make use of genetic data to estimate species boundaries, where one of the most widely used methods is the Generalized Mixed Yule Coalescent (GMYC) model. Despite its popularity, a number of factors are known to influence the performance and resulting inferences of the GMYC. Moreover, the few studies that have assessed model performance to date have been predominantly based on simulated data sets, where model assumptions are not violated. Here, we present a user‐friendly R Shiny application, ‘SPEDE‐sampler’ (SPEcies DElimitation sampler), that assesses the effect of computational and methodological choices, in combination with sampling effects, on the GMYC model. Output phylogenies are used to test the effect that (1) sample size, (2) BEAST and GMYC parameters (e.g. prior settings, single vs multiple threshold, clock model), and (3) singletons have on GMYC output. Optional predefined grouping information (e.g. morphospecies/ecotypes) can be uploaded in order to compare it with GMYC species and estimate percentage match scores. Additionally, predefined groups that contribute to inflated species richness estimates are identified by SPEDE‐sampler, allowing for the further investigation of potential cryptic species or geographical substructuring in those groups. Merging by the GMYC is also recorded to identify where traditional taxonomy has overestimated species numbers. Four worked examples are provided to illustrate the functionality of the program's workflow, and the variation that can arise when applying the GMYC model to empirical data sets. The R Shiny program is available for download at https://github.com/clarkevansteenderen/spede_sampler_R.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Mammoth Cave: A Hotspot of Subterranean Biodiversity in the United States

The Mammoth Cave System in the Interior Low Plateau karst region in central Kentucky, USA is a global hotspot of cave-limited biodiversity, particularly terrestrial species. We searched the literature, museum accessions, and database records to compile an updated list of troglobiotic and stygobiotic species for the Mammoth Cave System and compare our list with previously published checklists. Our list of cave-limited fauna totals 49 species, with 32 troglobionts and 17 stygobionts. Seven species are endemic to the Mammoth Cave System and other small caves in Mammoth Cave National Park. The Mammoth Cave System is the type locality for 33 cave-limited species. The exceptional diversity at Mammoth Cave is likely related to several factors, such as the high dispersal potential of cave fauna associated with expansive karst exposures, high surface productivity, and a long history of exploration and study. Nearly 80% of the cave-limited fauna is of conservation concern, many of which are at an elevated risk of extinction because of small ranges, few occurrences, and several potential threats.

Linking Hydrogeology and Ecology in Karst Landscapes: The Response of Epigean and Obligate Groundwater Copepods (Crustacea: Copepoda)

Groundwater invertebrate communities in karst landscapes are known to vary in response to multiple environmental factors. This study aims to explore the invertebrate assemblages’ composition of an Apennine karst system in Italy mainly described by the Rio Gamberale surface stream and the Stiffe Cave. The stream sinks into the carbonate rock and predominantly feeds the saturated karst into the cave. For a minor portion, groundwater flows from the epikarst and the perched aquifer within it. The spatial distribution of the species belonging to the selected target group of the Crustacea Copepoda between the surface stream and the groundwater habitats inside the cave highlighted a different response of surface-water species and obligate groundwater dwellers to the hydrogeological traits of the karst unit. Our results suggest that fast endorheic infiltration routes promoted the drift of epigean species from the surface to groundwater via the sinking stream while most of the obligate groundwater dwellers come from the perched aquifer in the epikarst from diffuse infiltration pathways.

Comparative development and ocular histology between epigean and subterranean salamanders (Eurycea) from central Texas

The salamander clade Eurycea from the karst regions of central Texas provides an ideal platform for comparing divergent nervous and sensory systems since some species exhibit extreme phenotypes thought to be associated with inhabiting a subterranean environment, including highly reduced eyes, while others retain an ancestral ocular phenotype appropriate for life above ground. We describe ocular morphology, comparing three salamander species representing two phenotypes-the surface-dwelling Barton Springs salamander (E. sosorum) and San Marcos salamander (E. nana) and the obligate subterranean Texas blind salamander (E. rathbuni) - in terms of structure and size of their eyes. Eyes were examined using confocal microscopy and measurements were made using ImageJ. Statistical analysis of data was carried out using R. We also provide a developmental series and track eye development and immunolocalization of Pax6 in E. sosorum and E. rathbuni. Adult histology of the surface-dwelling San Marcos salamander (E. nana) shows similarities to E. sosorum. The eyes of adults of the epigean species E. nana and E. sosorum appear fully developed with all the histological features of a fully functional eye. In contrast, the eyes of E. rathbuni adults have fewer layers, lack lenses and other features associated with vision as has been reported previously. However, in early developmental stages eye morphology did not differ significantly between E. rathbuni and E. sosorum. Parallel development is observed between the two phenotypes in terms of morphology; however, Pax6 labeling seems to decrease in the latter stages of development in E.rathbuni. We test for immunolabeling of the visual pigment proteins opsin and rhodopsin and observe immunolocalization around photoreceptor disks in E. nana and E. sosorum, but not in the subterranean E. rathbuni. Our results from examining developing salamanders suggest a combination of underdevelopment and degeneration contribute to the reduced eyes of adult E. rathbuni.

Phylogenomic species delimitation dramatically reduces species diversity in an Antarctic adaptive radiation

Application of genetic data to species delimitation often builds confidence in delimitations previously hypothesized using morphological, ecological, and geographic data and frequently yields recognition of previously-undescribed cryptic diversity. However, a recent critique of genomic data-based species delimitation approaches is that they have the potential to conflate population structure with species diversity, resulting in taxonomic oversplitting. The need for an integrative approach to species delimitation, in which molecular, morphological, ecological, and geographic lines of evidence are evaluated together, is becoming increasingly apparent. Here, we integrate phylogenetic, population genetic, and coalescent analyses of genome-wide sequence data with investigation of variation in multiple morphological traits to delimit species within the Antarctic barbeled plunderfishes (Artedidraconidae: Pogonophryne). Pogonophryne currently comprises 29 valid species, most of which are distinguished solely by variation in ornamentation of the mental barbel that projects from the lower jaw, a structure previously shown to vary widely within a single species. However, our genomic and phenotypic analyses result in a dramatic reduction in the number of distinct species recognized within the clade, providing evidence to support the recognition of no more than six species. We propose to synonymize 24 of the currently recognized species with five species of Pogonophryne. We find genomic and phenotypic evidence for a new species of Pogonophryne from specimens collected in the Ross Sea. Our findings represent a rare example in which application of molecular data provides evidence of taxonomic oversplitting on the basis of morphology, clearly demonstrating the utility of an integrative species delimitation framework.

wQFM: Highly Accurate Genome-scale Species Tree Estimation from Weighted Quartets

MOTIVATION

Species tree estimation from genes sampled from throughout the whole genome is complicated due to the gene tree-species tree discordance. Incomplete lineage sorting (ILS) is one of the most frequent causes for this discordance, where alleles can coexist in populations for periods that may span several speciation events. Quartet-based summary methods for estimating species trees from a collection of gene trees are becoming popular due to their high accuracy and statistical guarantee under ILS. Generating quartets with appropriate weights, where weights correspond to the relative importance of quartets, and subsequently amalgamating the weighted quartets to infer a single coherent species tree can allow for a statistically consistent way of estimating species trees. However, handling weighted quartets is challenging.

RESULTS

We propose wQFM, a highly accurate method for species tree estimation from multi-locus data, by extending the quartet FM (QFM) algorithm to a weighted setting. wQFM was assessed on a collection of simulated and real biological datasets, including the avian phylogenomic dataset which is one of the largest phylogenomic datasets to date. We compared wQFM with wQMC, which is the best alternate method for weighted quartet amalgamation, and with ASTRAL, which is one of the most accurate and widely used coalescent-based species tree estimation methods. Our results suggest that wQFM matches or improves upon the accuracy of wQMC and ASTRAL.

AVAILABILITY

wQFM is available in open source form at https://github.com/Mahim1997/wQFM-2020.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Subterranean movement inferred by temporary emigration in Barton Springs salamanders (Eurycea sosorum)

Movement behavior is an important aspect of animal ecology but is challenging to study in species that are unobservable for some portion of their lives, such as those inhabiting subterranean environments. Using four years of robust-design capture-recapture data, we examined the probability of movement into subterranean habitat by a population of endangered Barton Springs salamanders (Eurycea sosorum), a species that inhabits both surface and subterranean groundwater habitats. We tested the effects of environmental variables and body size on survival and temporary emigration, using the latter as a measure of subterranean habitat use. Based on 2,046 observations of 1,578 individuals, we found that temporary emigration was higher for larger salamanders, 79% of which temporarily emigrated into subterranean habitat between primary sampling intervals, on average. Body size was a better predictor of temporary emigration and survival compared to environmental covariates, although coefficients from lower ranked models suggested turbidity and dissolved oxygen may influence salamander movement between the surface and subsurface. Surface population dynamics are partly driven by movement below ground and therefore surface abundance estimates represent a fraction of the superpopulation. As such, while surface habitat management remains an important conservation strategy for this species, periodic declines in apparent surface abundance do not necessarily indicate declines of the superpopulation associated with the spring habitat.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

Phylogeography, historical demography and systematics of the world's smallest pythons (Pythonidae, Antaresia)

Advances from empirical studies in phylogeography, systematics and species delimitation highlight the importance of integrative approaches for quantifying taxonomic diversity. Genomic data have greatly improved our ability to discern both systematic diversity and evolutionary history. Here we combine analyses of mitochondrial DNA sequences, thousands of genome-wide SNPs and linear and geometric morphometrics on Antaresia, a clade of four currently recognised dwarf pythons from Australia and New Guinea (Antaresia childreni, A. stimsoni, A. maculosa and A. perthensis). Our integrative analyses of phylogenetics, population structure, species delimitation, historical demography and morphometrics revealed that the true evolutionary diversity is not well reflected in the current appraisal of the diversity of the group. We find that Antaresia childreni and A. stimsoni comprise a widespread network of populations connected by gene flow and without evidence of species-level divergence among them. However, A. maculosa shows considerable genetic structuring which leads us to recognise two subspecies in northeastern Australia and a new species in Torres Strait and New Guinea. These two contrasting cases of over and under estimation of diversity, respectively, illustrate the power of thorough integrative approaches into understanding evolution of biodiversity. Furthermore, our analyses of historical demographic patterns highlight the importance of the Kimberley, Pilbara and Cape York as origins of biodiversity in Australia.

Osmoregulatory ability predicts geographical range size in marine amniotes

Species that are distributed over wide geographical ranges are likely to encounter a greater diversity of environmental conditions than do narrowly distributed taxa, and thus we expect a correlation between size of geographical range and breadth of physiological tolerances to abiotic challenges. That correlation could arise either because higher physiological capacity enables range expansion, or because widely distributed taxa experience more intense (but spatially variable) selection on physiological tolerances. The invasion of oceanic habitats by amniotic vertebrates provides an ideal system with which to test the predicted correlation between range size and physiological tolerances, because all three lineages that have secondarily moved into marine habitats (mammals, birds, reptiles) exhibit morphological and physiological adaptations to excrete excess salt. Our analyses of data on 62 species (19 mammals, 18 birds, 24 reptiles) confirm that more-widely distributed taxa encounter habitats with a wider range of salinities, and that they have higher osmoregulatory ability as determined by sodium concentrations in fluids expelled from salt-excreting organs. This result remains highly significant even in models that incorporate additional explanatory variables such as metabolic mode, body size and dietary habits. Physiological data thus may help to predict potential range size and perhaps a species' vulnerability to anthropogenic disturbance.

Stygobromus bakeri, a new species of groundwater amphipod (Amphipoda, Crangonyctidae) associated with the Trinity and Edwards aquifers of central Texas, USA

A new stygobitic groundwater amphipod species, Stygobromus bakerisp. nov., is described from 4 central Texas limestone karst springs; John Knox Spring (Comal County), Jacob’s Well (Hays County), Mormon Spring (Travis County) and Salado Springs (Bell County). This species belongs to the predominately western Nearctic hubbsi species group of Stygobromus and differs from other species in the group by gnathopod 1 with row of 4 to 5 setae posterior to the defining angle, pereopods 6 and 7 with broadly expanded bases and distinct distoposterior lobes, gnathopods 1 and 2 with 2 rows of 3 singly inserted setae on the inner palm, and uropod 3 with single, slightly distal peduncular seta. Habitat, sympatric groundwater species, and conservation issues are discussed.

Reliable Evapotranspiration Predictions with a Probabilistic Machine Learning Framework

Evapotranspiration is often expressed in terms of reference crop evapotranspiration (ETo), actual evapotranspiration (ETa), or surface water evaporation (Esw), and their reliable predictions are critical for groundwater, irrigation, and aquatic ecosystem management in semi-arid regions. We demonstrated that a newly developed probabilistic machine learning (ML) model, using a hybridized “boosting” framework, can simultaneously predict the daily ETo, Esw, & ETa from local hydroclimate data with high accuracy. The probabilistic approach exhibited great potential to overcome data uncertainties, in which 100% of the ETo, 89.9% of the Esw, and 93% of the ETa test data at three watersheds were within the models’ 95% prediction intervals. The modeling results revealed that the hybrid boosting framework can be used as a reliable computational tool to predict ETo while bypassing net solar radiation calculations, estimate Esw while overcoming uncertainties associated with pan evaporation & pan coefficients, and predict ETa while offsetting high capital & operational costs of EC towers. In addition, using the Shapley analysis built on a coalition game theory, we identified the order of importance and interactions between the hydroclimatic variables to enhance the models’ transparency and trustworthiness.

Gene flow creates a mirage of cryptic species in a Southeast Asian spotted stream frog complex

Most new cryptic species are described using conventional tree- and distance-based species delimitation methods (SDMs), which rely on phylogenetic arrangements and measures of genetic divergence. However, although numerous factors such as population structure and gene flow are known to confound phylogenetic inference and species delimitation, the influence of these processes is not frequently evaluated. Using large numbers of exons, introns, and ultraconserved elements obtained using the FrogCap sequence-capture protocol, we compared conventional SDMs with more robust genomic analyses that assess population structure and gene flow to characterize species boundaries in a Southeast Asian frog complex (Pulchrana picturata). Our results showed that gene flow and introgression can produce phylogenetic patterns and levels of divergence that resemble distinct species (up to 10% divergence in mitochondrial DNA). Hybrid populations were inferred as independent (singleton) clades that were highly divergent from adjacent populations (7%-10%) and unusually similar (<3%) to allopatric populations. Such anomalous patterns are not uncommon in Southeast Asian amphibians, which brings into question whether the high levels of cryptic diversity observed in other amphibian groups reflect distinct cryptic species-or, instead, highly admixed and structured metapopulation lineages. Our results also provide an alternative explanation to the conundrum of divergent (sometimes nonsister) sympatric lineages-a pattern that has been celebrated as indicative of true cryptic speciation. Based on these findings, we recommend that species delimitation of continuously distributed "cryptic" groups should not rely solely on conventional SDMs, but should necessarily examine population structure and gene flow to avoid taxonomic inflation.

Defining relictual biodiversity: Conservation units in speckled dace (Leuciscidae: Rhinichthys osculus) of the Greater Death Valley ecosystem

The tips in the tree of life serve as foci for conservation and management, yet clear delimitations are masked by inherent variance at the species-population interface. Analyses using thousands of nuclear loci can potentially sort inconsistencies, yet standard categories applied to this parsing are themselves potentially conflicting and/or subjective [e.g., DPS (distinct population segments); DUs (Diagnosable Units-Canada); MUs (management units); SSP (subspecies); ESUs (Evolutionarily Significant Units); and UIEUs (uniquely identified evolutionary units)]. One potential solution for consistent categorization is to create a comparative framework by accumulating statistical results from independent studies and evaluating congruence among data sets. Our study illustrates this approach in speckled dace (Leuciscidae: Rhinichthys osculus) endemic to two basins (Owens and Amargosa) in the Death Valley ecosystem. These fish persist in the Mojave Desert as isolated Plio-Pleistocene relicts and are of conservation concern, but lack formal taxonomic descriptions/designations. Double digest RAD (ddRAD) methods identified 14,355 SNP loci across 10 populations (N = 140). Species delimitation analyses [multispecies coalescent (MSC) and unsupervised machine learning (UML)] delineated four putative ESUs. F _ST outlier loci (N = 106) were juxtaposed to uncover the potential for localized adaptations. We detected one hybrid population that resulted from upstream reconnection of habitat following contemporary pluvial periods, whereas remaining populations represent relics of ancient tectonism within geographically isolated springs and groundwater-fed streams. Our study offers three salient conclusions: a blueprint for a multifaceted delimitation of conservation units; a proposed mechanism by which criteria for intraspecific biodiversity can be potentially standardized; and a strong argument for the proactive management of critically endangered Death Valley ecosystem fishes.

Reticulate evolution as a management challenge: Patterns of admixture with phylogenetic distance in endemic fishes of western North America

Admixture in natural populations is a long-standing management challenge, with population genomic approaches offering means for adjudication. We now more clearly understand the permeability of species boundaries and the potential of admixture for promoting adaptive evolution. These issues particularly resonate in western North America, where tectonism and aridity have fragmented and reshuffled rivers over millennia, in turn promoting reticulation among endemic fishes, a situation compounded by anthropogenic habitat modifications and non-native introductions. The melding of historic and contemporary admixture has both confused and stymied management. We underscore this situation with a case study that quantifies basin-wide admixture among a group of native and introduced fishes by employing double-digest restriction site-associated DNA (ddRAD) sequencing. Our approach: (a) quantifies the admixed history of 343 suckers (10 species of Catostomidae) across the Colorado River Basin; (b) gauges admixture within the context of phylogenetic distance and "ecological specialization"; and (c) extrapolates potential drivers of introgression across hybrid crosses that involve endemic as well as invasive species. Our study extends across an entire freshwater basin and expands previous studies more limited in scope both geographically and taxonomically. Our results detected admixture involving all 10 species, with habitat alterations not only accelerating the breakdown of reproductive isolation, but also promoting introgression. Hybridization occurred across the genus despite phylogenetic distance, whereas introgression was only detected within subgenera, implicating phylogenetic distance and/or ecological specialization as drivers of reproductive isolation. Understanding the extent of admixture and reproductive isolation across multiple species serves to disentangle their reticulate evolutionary histories and provides a broadscale perspective for basin-wide conservation and management.

Evidence for rapid phenotypic and behavioural shifts in a recently established cavefish population

Cave colonization offers a natural laboratory to study an extreme environmental shift, and diverse cave species from around the world often have converged on robust morphological, physiological and behavioural traits. The Mexican tetra (Astyanax mexicanus) has repeatedly colonized caves in the Sierra de El Abra and Sierra de Guatemala regions of north-east Mexico ~0.20–1 Mya, indicating an ability to adapt to the cave environment. The time frame for the evolution of these traits in any cave animal, however, is poorly understood. Astyanax mexicanus from the Río Grande in South Texas were brought to Central Texas beginning in the early 1900s and colonized underground environments. Here, we investigate whether phenotypic and behavioural differences have occurred rapidly between a surface population and a geographically proximate cave population, probably of recent origin. Fish from the cave and surface populations differ significantly in morphological traits, including coloration, lateral line expansion and dorsal fin placement. Striking behavioural shifts in aggression, feeding and wall-following have also occurred. Together, our results suggest that morphological and behavioural changes accompanying cave colonization can be established rapidly, and this system offers an exciting and unique opportunity for isolating the genetic and environmental contributions to colonization of extreme environments.

Hydrologic variability contributes to reduced survival through metamorphosis in a stream salamander

Changes in the amount, intensity, and timing of precipitation are increasing hydrologic variability in many regions, but we have little understanding of how these changes are affecting freshwater species. Stream-breeding amphibians-a diverse group in North America-may be particularly sensitive to hydrologic variability during aquatic larval and metamorphic stages. Here, we tested the prediction that hydrologic variability in streams decreases survival through metamorphosis in the salamander Gyrinophilus porphyriticus, reducing recruitment to the adult stage. Using a 20-y dataset from Merrill Brook, a stream in northern New Hampshire, we show that abundance of G. porphyriticus adults has declined by ∼50% since 1999, but there has been no trend in larval abundance. We then tested whether hydrologic variability during summers influences survival through metamorphosis, using capture-mark-recapture data from Merrill Brook (1999 to 2004) and from 4 streams in the Hubbard Brook Experimental Forest (2012 to 2014), also in New Hampshire. At both sites, survival through metamorphosis declined with increasing variability of stream discharge. These results suggest that hydrologic variability reduces the demographic resilience and adaptive capacity of G. porphyriticus populations by decreasing recruitment of breeding adults. They also provide insight on how increasing hydrologic variability is affecting freshwater species, and on the broader effects of environmental variability on species with vulnerable metamorphic stages.

Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae)

Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.