High-stakes species delimitation in eyeless cave spiders (Cicurina, Dictynidae, Araneae) from central Texas

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.

Cryptic diversity begets challenges and opportunities in biodiversity research

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

The phylogenetic structure and coalescent species delimitation of an endemic trapdoor spider genus, Stasimopus (Araneae, Mygalomorphae, Stasimopidae) in the Karoo region of South Africa

The Karoo region of South Africa is a unique and sensitive ecosystem which is facing pressure for development due to economic incentives such as mining, farming and shale gas exploration. The species diversity of many taxa in the area is largely unknown. A phylogenetic analysis of the cork-lid trapdoor spider genus, Stasimopus (Stasimopidae) was undertaken in order to gain insight into the relationships between the species that may be present in the area. The species within Stasimopus are challenging to identify and define using traditional morphological methods due to a high degree of morphological conservatism within the genus. For this reason, multiple coalescent based species delimitation methods were used to attempt to determine the species present for Stasimopus in the region which was tested against the morphological identifications and genetic clades (based on CO1, 16S and EF-1ɣ). We tested single-locus methods Automatic Barcode Gap Discovery (ABGD), Bayesian implementation of Poisson Tree Processes (bPTP) and General Mixed Yule- Coalescent (GMYC), as well as multi-locus Brownie. The phylogenetic analysis of Stasimopus in the Karoo showed that there is a high degree of genetic diversity within the genus. The species delimitation results proved unfruitful for the genus, as they appear to delimit population structure rather than species for most methods. Alternative methods should be investigated to aid in the identification of the species in order truly understand the species diversity of the genus.

Deep mtDNA Sequence Divergences and Possible Species Radiation of Whip Spiders (Arachnida, Amblypygi, Phrynidae, Phrynus/Paraphrynus) among Caribbean Oceanic and Cave Islands

Islands—whether classic oceanic islands or habitat islands such as isolated thermal vents, mountain tops, or caves—often promote the diversification of lineages that colonize them. We examined CO1 mtDNA sequence divergences within the tailless whip spider genus Phrynus Lamarck, 1809 (Amblypygi: Phrynidae) among oceanic islands and among cave ’islands´ distributed across the Caribbean archipelago and on the continental mainland. The significance of this study lies in the extensive taxon sampling of a supposedly depauperate lineage (considering its age), over a large proportion of its geographical range, and the discovery of deep mtDNA sequence divergences. We sampled thousands of specimens—and sequenced 544, including six outgroup species—across 173 localities on 17 islands (135 localities) and five countries on the North to South American mainland (38 localities), including a total of 63 caves. Classical taxonomy identified ten named Phrynus and two Paraphrynus Moreno, 1940 species. Paraphrynus seems to be paraphyletic and nested in Phrynus. Uncorrected genetic distances within named species and among morphological species ranged up to 15% and 19%, respectively. Geographic distances explained a significant portion of genetic distances on islands (19%, among both subterranean and epigean specimens), and for epigean specimens on the mainland (27%). Species delimitation analyses indicated that the 12 named species harbored from 66 to well over 100 putative species. The highest number of species was indicated by the GMYC method (114 species) while the Bayesian Poisson tree processes (bPTP) and the BP&P relaying on the Markov chain Monte Carlo Bayesian Phylogenetic model estimated an upper level of 110 species. On the other hand, the recently recommended and relatively conservative distance-based (phylogeny free) ASAP model has the greatest support for 73 species. In either case, nearly all putative species are tightly limited to a single locality, often a small cave system, and sometimes to the surrounding epigean area. Caribbean Phrynus diversity has likely been vastly underestimated, likely due to both morphological crypsis and the ignorance of Caribbean cave fauna. Although mtDNA sequences can suggest species limits, nuclear DNA sequencing and detailed morphological research are necessary to corroborate them and explore whether this phenomenon constitutes species radiation or perhaps just mtDNA divergences as a consequence of, for example, stationary females and actively dispersing males.

Biogeography and eye size evolution of the ogre-faced spiders

Net-casting spiders (Deinopidae) comprise a charismatic family with an enigmatic evolutionary history. There are 67 described species of deinopids, placed among three genera, Deinopis, Menneus, and Asianopis, that are distributed globally throughout the tropics and subtropics. Deinopis and Asianopis, the ogre-faced spiders, are best known for their giant light-capturing posterior median eyes (PME), whereas Menneus does not have enlarged PMEs. Molecular phylogenetic studies have revealed discordance between morphology and molecular data. We employed a character-rich ultra-conserved element (UCE) dataset and a taxon-rich cytochrome-oxidase I (COI) dataset to reconstruct a genus-level phylogeny of Deinopidae, aiming to investigate the group's historical biogeography, and examine PME size evolution. Although the phylogenetic results support the monophyly of Menneus and the single reduction of PME size in deinopids, these data also show that Deinopis is not monophyletic. Consequently, we formally transfer 24 Deinopis species to Asianopis; the transfers comprise all of the African, Australian, South Pacific, and a subset of Central American and Mexican species. Following the divergence of Eastern and Western deinopids in the Cretaceous, Deinopis/Asianopis dispersed from Africa, through Asia and into Australia with its biogeographic history reflecting separation of Western Gondwana as well as long-distance dispersal events.

Phylogenomics of paleoendemic lampshade spiders (Araneae, Hypochilidae, Hypochilus), with the description of a new species from montane California

Hypochilus is a relictual lineage of Nearctic spiders distributed disjunctly across the United States in three montane regions (California, southern Rocky Mountains, southern Appalachia). Phylogenetic resolution of species relationships in Hypochilus has been challenging, and conserved morphology coupled with extreme genetic divergence has led to uncertain species limits in some complexes. Here, Hypochilus interspecies relationships have been reconstructed and cryptic speciation more critically evaluated using a combination of ultraconserved elements, mitochondrial CO1 by-catch, and morphology. Phylogenomic data strongly support the monophyly of regional clades and support a ((California, Appalachia), southern Rocky Mountains) topology. In Appalachia, five species are resolved as four lineages (H.thorelli Marx, 1888 and H.coylei Platnick, 1987 are clearly sister taxa), but the interrelationships of these four lineages remain unresolved. The Appalachian species H.pococki Platnick, 1987 is recovered as monophyletic but is highly genetically structured at the nuclear level. While algorithmic analyses of nuclear data indicate many species (e.g., all H.pococki populations as species), male morphology instead reveals striking stasis. Within the California clade, nuclear and mitochondrial lineages of H.petrunkevitchi Gertsch, 1958 correspond directly to drainage basins of the southern Sierra Nevada, with H.bernardino Catley, 1994 nested within H.petrunkevitchi and sister to the southernmost basin populations. Combining nuclear, mitochondrial, geographical, and morphological evidence a new species from the Tule River and Cedar Creek drainages is described, Hypochilusxomotesp. nov. We also emphasize the conservation issues that face several microendemic, habitat-specialized species in this remarkable genus.

Improving Taxonomic Practices and Enhancing Its Extensibility—An Example from Araneology

Planetary extinction of biodiversity underscores the need for taxonomy. Here, we scrutinize spider taxonomy over the last decade (2008–2018), compiling 2083 published accounts of newly described species. We evaluated what type of data were used to delineate species, whether data were made freely available, whether an explicit species hypothesis was stated, what types of media were used, the sample sizes, and the degree to which species constructs were integrative. The findings we report reveal that taxonomy remains largely descriptive, not integrative, and provides no explicit conceptual framework. Less than 4% of accounts explicitly stated a species concept and over one-third of all new species described were based on 1–2 specimens or only one sex. Only ~5% of studies made data freely available, and only ~14% of all newly described species employed more than one line of evidence, with molecular data used in ~6% of the studies. These same trends have been discovered in other animal groups, and therefore we find it logical that taxonomists face an uphill challenge when justifying the scientific rigor of their field and securing the needed resources. To move taxonomy forward, we make recommendations that, if implemented, will enhance its rigor, repeatability, and scientific standards.

Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web

High throughput sequencing and phylogenomic analyses focusing on relationships among spiders have both reinforced and upturned long-standing hypotheses. Likewise, the evolution of spider webs-perhaps their most emblematic attribute-is being understood in new ways. With a matrix including 272 spider species and close arachnid relatives, we analyze and evaluate the relationships among these lineages using a variety of orthology assessment methods, occupancy thresholds, tree inference methods and support metrics. Our analyses include families not previously sampled in transcriptomic analyses, such as Symphytognathidae, the only araneoid family absent in such prior works. We find support for the major established spider lineages, including Mygalomorphae, Araneomorphae, Synspermiata, Palpimanoidea, Araneoidea and the Retrolateral Tibial Apophysis Clade, as well as the uloborids, deinopids, oecobiids and hersiliids Grade. Resulting trees are evaluated using bootstrapping, Shimodaira-Hasegawa approximate likelihood ratio test, local posterior probabilities and concordance factors. Using structured Markov models to assess the evolution of spider webs while accounting for hierarchically nested traits, we find multiple convergent occurrences of the orb web across the spider tree-of-life. Overall, we provide the most comprehensive spider tree-of-life to date using transcriptomic data and use new methods to explore controversial issues of web evolution, including the origins and multiple losses of the orb web.

Incorporating the speciation process into species delimitation

The "multispecies" coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.

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.

A review of the taxonomy of spiny-backed orb-weaving spiders of the subfamily Gasteracanthinae (Araneae, Araneidae) in Thailand

Spiny-backed orb-weaving spiders of the subfamily Gasteracanthinae are broadly distributed in the Old World. Despite their use as a model species in biology, evolution, and behavior because of their extraordinary characteristics, the systematics of this group of spiders are still poorly understood. This study elucidates the systematics of Gasteracanthinae in Thailand based on morphological and molecular-based analyses. In total, seven species from three genera, namely Gasteracantha, Macracantha, and Thelacantha, were recorded in Thailand. Shape of abdominal spines, pattern of sigilla, and female genitalia are significant characters for species identification. In contrast, coloration shows highly intraspecific variation in most species within Gasteracanthinae. A phylogenetic tree based on partial sequences of COI, 16S, and H3 genes recovered Gasteracanthinae as a monophyletic group and supports the existence of three clades. Gasteracantha hasselti is placed as a sister taxon to Macracantha arcuata. Hence, we propose to transfer G. hasselti to Macracantha. Moreover, molecular species delimitation analyses (ABGD, bPTP, and GMYC) using 675 bp of COI gene support all nominal species, with evidence of possible additional cryptic species.

Paraphyly and evolutionary independent lineages in Gymnotus pantherinus (Gymnotiformes: Gymnotidae) in the Brazilian Atlantic Forest Coastal Streams

In the Brazilian Atlantic Forest (AF), many terrestrial species with broad geographical distributions show high diversity and endemism of intraspecific lineages, as revealed by molecular genetic data. This pattern, however, is less explored in freshwater fishes. Gymnotus pantherinus is an electric fish endemic to the Brazilian coastal drainages that shows a wide distribution, ranging from the states of Bahia to Santa Catarina, an unusual pattern for AF fishes. It has been hypothesized that G. pantherinus is a species complex because distinct morphotypes were described for the species based on morphometric and meristic features. We used mitochondrial and nuclear data to test this hypothesis. Based on phylogenetic inference and multi-locus, multispecies coalescent methods, we identified six independent lineages, flagging them as candidate species. One such lineage is the recently described species G. refugio that is nested within G. pantherinus and renders it paraphyletic, showing it is a species complex. We named G. pantherinus stricto sensu the lineage that includes samples from the type locality (Santos, SP). Our results show that genetic lineages correspond only partially and far exceed the number of previously reported morphotypes. Genetic breaks in the group correspond to landscape features associated with the Serra do Mar mountain range and with riverine dynamics caused by sea level changes during the last glacial maximum. Moreover, we found evidence of river capture events affecting phylogeographic structure in the group. We uncovered an important dimension of diversity in the group and encourage further integration of genetic and phenotypic data. Such integration is a fruitful approach not only to reduce the gap between taxonomy and evolutionary history in Gymnotidae, but also to uncover the real AF biodiversity.

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.

Local- versus broad-scale environmental drivers of continental β-diversity patterns in subterranean spider communities across Europe

Macroecologists seek to identify drivers of community turnover (β-diversity) through broad spatial scales. However, the influence of local habitat features in driving broad-scale β-diversity patterns remains largely untested, owing to the objective challenges of associating local-scale variables to continental-framed datasets. We examined the relative contribution of local- versus broad-scale drivers of continental β-diversity patterns, using a uniquely suited dataset of cave-dwelling spider communities across Europe (35-70° latitude). Generalized dissimilarity modelling showed that geographical distance, mean annual temperature and size of the karst area in which caves occurred drove most of β-diversity, with differential contributions of each factor according to the level of subterranean specialization. Highly specialized communities were mostly influenced by geographical distance, while less specialized communities were mostly driven by mean annual temperature. Conversely, local-scale habitat features turned out to be meaningless predictors of community change, which emphasizes the idea of caves as the human accessible fraction of the extended network of fissures that more properly represents the elective habitat of the subterranean fauna. To the extent that the effect of local features turned to be inconspicuous, caves emerge as experimental model systems in which to study broad biological patterns without the confounding effect of local habitat features.

Molecular phylogeny and species delimitation of the genus Dicerapanorpa (Mecoptera: Panorpidae)

Given that species is the fundamental unit in systematic biology, rigorous species delimitation is crucial for taxonomic studies, yet routine species delimitation remains an ongoing challenge in the taxonomic practice of insects. The two-horned scorpionfly Dicerapanorpa is a small genus in Panorpidae (Mecoptera) endemic to the Qinling-Bashan and Hengduan mountains, a biodiversity hotspot. However, species of Dicerapanorpa are difficult to delineate owing to marked intraspecific variation and interspecific similarity. Here, we investigate the diversity and species boundaries of Dicerapanorpa using an integrative approach based on DNA barcoding, morphological, geometric morphometric and molecular phylogenetic analyses. This integrative analyses confirmed the 13 described species of Dicerapanorpa and revealed three new species: Dicerapanorpa lativalva sp. nov., Dicerapanorpa hualongshana sp. nov. and Dicerapanorpa minshana sp. nov. Most molecular operational taxonomic units are in congruence with morphological clusters. Possible reasons for several discordances in Dicerapanorpa are tentatively discussed.

Continental data on cave-dwelling spider communities across Europe (Arachnida: Araneae)

BACKGROUND

Spiders (Arachnida: Araneae) are widespread in subterranean ecosystems worldwide and represent an important component of subterranean trophic webs. Yet, global-scale diversity patterns of subterranean spiders are still mostly unknown. In the frame of the CAWEB project, a European joint network of cave arachnologists, we collected data on cave-dwelling spider communities across Europe in order to explore their continental diversity patterns. Two main datasets were compiled: one listing all subterranean spider species recorded in numerous subterranean localities across Europe and another with high resolution data about the subterranean habitat in which they were collected. From these two datasets, we further generated a third dataset with individual geo-referenced occurrence records for all these species.

NEW INFORMATION

Data from 475 geo-referenced subterranean localities (caves, mines and other artificial subterranean sites, interstitial habitats) are herein made available. For each subterranean locality, information about the composition of the spider community is provided, along with local geomorphological and habitat features. Altogether, these communities account for > 300 unique taxonomic entities and 2,091 unique geo-referenced occurrence records, that are made available via the Global Biodiversity Information Facility (GBIF) (Mammola and Cardoso 2019). This dataset is unique in that it covers both a large geographic extent (from 35° south to 67° north) and contains high-resolution local data on geomorphological and habitat features. Given that this kind of high-resolution data are rarely associated with broad-scale datasets used in macroecology, this dataset has high potential for helping researchers in tackling a range of biogeographical and macroecological questions, not necessarily uniquely related to arachnology or subterranean biology.

A demonstration of unsupervised machine learning in species delimitation

One major challenge to delimiting species with genetic data is successfully differentiating population structure from species-level divergence, an issue exacerbated in taxa inhabiting naturally fragmented habitats. Many fields of science are now using machine learning, and in evolutionary biology supervised machine learning has recently been used to infer species boundaries. These supervised methods require training data with associated labels. Conversely, unsupervised machine learning (UML) uses inherent data structure and does not require user-specified training labels, potentially providing more objectivity in species delimitation. In the context of integrative taxonomy, we demonstrate the utility of three UML approaches (random forests, variational autoencoders, t-distributed stochastic neighbor embedding) for species delimitation in an arachnid taxon with high population genetic structure (Opiliones, Laniatores, Metanonychus). We find that UML approaches successfully cluster samples according to species-level divergences and not high levels of population structure, while model-based validation methods severely over-split putative species. UML offers intuitive data visualization in two-dimensional space, the ability to accommodate various data types, and has potential in many areas of systematic and evolutionary biology. We argue that machine learning methods are ideally suited for species delimitation and may perform well in many natural systems and across taxa with diverse biological characteristics.

'At the crossroads towards polyploidy': genomic divergence and extent of homoploid hybridization are drivers for the formation of the ox-eye daisy polyploid complex (Leucanthemum, Compositae-Anthemideae)

Polyploidy plays a paramount role in phytodiversity, but the causes of this evolutionary pathway require further study. Here, we use phylogenetic methods to examine possible polyploidy-promoting factors by comparing diploid representatives of the comprehensive European polyploid complex Leucanthemum with members of its strictly diploid North African counterpart Rhodanthemum. We investigate genetic divergence and gene flow among all diploid lineages of both genera to evaluate the role of genomic differentiation and hybridization for polyploid speciation. To test whether hybridization in Leucanthemum has been triggered by the geological conditions during its diversification, we additionally generate a time-calibrated phylogeny of 46 species of the subtribe Leucantheminae. Leucanthemum shows a significantly higher genetic divergence and hybridization signal among diploid lineages compared with Rhodanthemum, in spite of a similar crown age and diversification pattern during the Quaternary. Our study demonstrates the importance of genetic differentiation among diploid progenitors and their concurrent affinity for natural hybridization for the formation of a polyploid complex. Furthermore, the role of climate-induced range overlaps on hybridization and polyploid speciation during the Quaternary is discussed.

The Kingdom of the Blind: Disentangling Fundamental Drivers in the Evolution of Eye Loss

Light is a fundamentally important biological cue used by almost every animal on earth, to maintain daily rhythms, navigate, forage, find mates, or avoid predators. But an enormous number of species live in darkness: in subterranean caves, deep oceans, underground burrows, and within parasitic host bodies, and the loss of eyes appears consistently across these ecosystems. However, the evolutionary mechanisms that lead to the reduction of the visual system remain the subject of great interest and debate more than 150 years after Darwin tackled the issue. Studies of model taxa have discovered significant roles for natural selection, neutral evolution, and pleiotropy, but the interplay between them remains unclear. To nail down unifying concepts surrounding the evolution of eye loss, we must embrace the enormous range of affected animals and habitats. The fine developmental details of model systems such as the Mexican cave tetra Astyanax mexicanus have transformed and enriched the field, but these should be complemented by wider studies to identify truly overarching patterns that apply throughout animals. Here, the major evolutionary drivers are placed within a conceptual cost-benefit framework that incorporates the fundamental constraints and forces that influence evolution in the dark. Major physiological, ecological, and environmental factors are considered within the context of this framework, which appears faithful to observed patterns in deep-sea and cavernicolous animals. To test evolutionary hypotheses, a comparative phylogenetic approach is recommended, with the goal of studying large groups exhibiting repeated reduction, and then comparing these across habitats, taxa, and lifestyles. Currently, developmental and physiological methods cannot feasibly be used on such large scales, but penetrative imaging techniques could provide detailed morphological data non-invasively and economically for large numbers of species. Comprehensive structural datasets can then be contextualized phylogenetically to examine recurrent trends and associations, and to reconstruct character histories through multiple independent transitions into darkness. By assessing these evolutionary trajectories within an energetic cost-benefit framework, the relationships between fundamental influences can be inferred and compared across different biological and physical parameters. However, substantial numbers of biological and environmental factors affect the evolutionary trajectory of loss, and it is critical that researchers make fair and reasonable comparisons between objectively similar groups.

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

Groundwater-dependent species are among the least-known components of global biodiversity, as well as some of the most vulnerable because of rapid groundwater depletion at regional and global scales. The karstic Edwards-Trinity aquifer system of west-central Texas is one of the most species-rich groundwater systems in the world, represented by dozens of endemic groundwater-obligate species with narrow, naturally fragmented distributions. Here, we examine how geomorphological and hydrogeological processes have driven population divergence and speciation in a radiation of salamanders (Eurycea) endemic to the Edwards-Trinity system using phylogenetic and population genetic analysis of genome-wide DNA sequence data. Results revealed complex patterns of isolation and reconnection driven by surface and subsurface hydrology, resulting in both adaptive and nonadaptive population divergence and speciation. Our results uncover cryptic species diversity and refine the borders of several threatened and endangered species. The US Endangered Species Act has been used to bring state regulation to unrestricted groundwater withdrawals in the Edwards (Balcones Fault Zone) Aquifer, where listed species are found. However, the Trinity and Edwards-Trinity (Plateau) aquifers harbor additional species with similarly small ranges that currently receive no protection from regulatory programs designed to prevent groundwater depletion. Based on regional climate models that predict increased air temperature, together with hydrologic models that project decreased springflow, we conclude that Edwards-Trinity salamanders and other codistributed groundwater-dependent organisms are highly vulnerable to extinction within the next century.

Mitochondrial sequence data indicate “Vicariance by Erosion” as a mechanism of species diversification in North American Ptomaphagus (Coleoptera, Leiodidae, Cholevinae) cave beetles

Small carrion beetles (Coleoptera: Leiodidae: Cholevinae) are members of cave communities around the world and important models for understanding the colonization of caves, adaptation to cave life, and the diversification of cave-adapted lineages. We developed a molecular phylogeny to examine the diversification of the hirtus-group of the small carrion beetle genus Ptomaphagus. The hirtus-group has no surface-dwelling members; it consists of 19 short-range endemic cave- and soil-dwelling species in the central and southeastern United States of America. Taxonomic, phylogenetic and biogeographic data were previously interpreted to suggest the hirtus-group diversified within the past 350,000 years through a series of cave colonization and speciation events related to Pleistocene climate fluctuations. However, our time-calibrated molecular phylogeny resulting from the analysis of 2,300 nucleotides from five genes across three mitochondrial regions (cox1, cytb, rrnL-trnL-nad1) for all members of the clade paints a different picture. We identify three stages of diversification in the hirtus-group: (1) ~10 million years ago (mya), the lineage that develops into P.shapardi, a soil-dwelling species from the Ozarks, diverged from the lineage that gives rise to the 18 cave-obligate members of the group; (2) between 8.5 mya and 6 mya, seven geographically distinct lineages diverged across Kentucky, Tennessee, Alabama and Georgia; six of these lineages represent a single species today, whereas (3) the ‘South Cumberlands’ lineage in Tennessee and Alabama diversified into 12 species over the past ~6 my. While the events triggering diversification during the first two stages remain to be determined, the distributions, phylogenetic relationships and divergence times in the South Cumberlands lineage are consistent with populations being isolated by vicariant events as the southern Cumberland Plateau eroded and fragmented over millions of years.

Integrative taxonomy resolves taxonomic uncertainty for freshwater mussels being considered for protection under the U.S. Endangered Species Act

Objectively delimiting species boundaries remains an important challenge in systematics and becomes urgent when unresolved taxonomy complicates conservation and recovery efforts. We examined species boundaries in the imperiled freshwater mussel genus Cyclonaias (Bivalvia: Unionidae) using morphometrics, molecular phylogenetics, and multispecies coalescent models to help guide pending conservation assessments and legislative decisions. Congruence across multiple lines of evidence indicated that current taxonomy overestimates diversity in the C. pustulosa species complex. The only genetically and morphologically diagnosable species in the C. pustulosa species complex were C. pustulosa and C. succissa and we consider C. aurea, C. houstonensis, C. mortoni, and C. refulgens to be synonyms of C. pustulosa. In contrast, all three species in the C. nodulata complex (C. necki, C. nodulata, and C. petrina) were genetically, geographically, and morphologically diagnosable. Our findings have important conservation and management implications, as three nominal species (C. aurea, C. houstonensis, and C. petrina) are being considered for protection under the Endangered Species Act.

Sequence capture phylogenomics of eyeless Cicurina spiders from Texas caves, with emphasis on US federally-endangered species from Bexar County (Araneae, Hahniidae)

Morphological, mitochondrial, and nuclear phylogenomic data were combined to address phylogenetic and species delimitation questions in cave-limited Cicurina spiders from central Texas. Special effort was focused on specimens and cave locations in the San Antonio region (Bexar County), home to four eyeless species listed as US Federally Endangered. Sequence capture experiments resulted in the recovery of ~200-400 homologous ultra-conserved element (UCE) nuclear loci across taxa, and nearly complete COI mitochondrial DNA sequences from the same set of individuals. Some of these nuclear and mitochondrial sequences were recovered from "standard" museum specimens without special preservation of DNA material, including museum specimens preserved in the 1990s. Multiple phylogenetic analyses of the UCE data agree in the recovery of two major lineages of eyeless Cicurina in Texas. These lineages also differ in mitochondrial clade membership, female genitalic morphology, degree of troglomorphy (as measured by relative leg length), and are mostly allopatric across much of Texas. Rare sympatry was confirmed in Bexar County, where members of the two major clades sometimes co-exist in the same karst feature. Both nuclear phylogenomic and mitochondrial data indicate the existence of undescribed species from the San Antonio region, although further sampling and collection of adult specimens is needed to explicitly test these hypotheses. Our data support the two following species synonymies (Cicurina venii Gertsch, 1992 = Cicurina madla Gertsch, 1992; Cicurina loftini Cokendolpher, 2004 = Cicurina vespera Gertsch, 1992), formally proposed here. Overall, our taxonomy-focused research has many important conservation implications, and again highlights the fundamental importance of robust taxonomy in conservation research.

Prevalence of cryptic species in morphologically uniform taxa - Fast speciation and evolutionary radiation in Asian frogs

Diversity and distributions of cryptic species have long been a vexing issue. Identification of species boundaries is made difficult by the lack of obvious morphological differences. Here, we investigate the cryptic diversity and evolutionary history of an underappreciated group of Asian frog species (Megophrys) to explore the pattern and dynamic of amphibian cryptic species. We sequenced four mitochondrial genes and five nuclear genes and delineated species using multiple approaches, combining DNA and mating-call data. A Bayesian species tree was generated to estimate divergence times and to reconstruct ancestral ranges. Macroevolutionary analyses and hybridization tests were conducted to explore the evolutionary dynamics of this cryptic group. Our phylogenies support the current subgenera. We revealed 43 cryptic species, 158% higher than previously thought. The species-delimitation results were further confirmed by mating-call data and morphological divergence. We found that these Asian frogs entered China from the Sunda Shelf 48 Mya, followed by an ancient radiation event during middle Miocene. We confirmed the efficiency of the multispecies coalescent model for delimitation of species with low morphological diversity. Species diversity of Megophrys is severely underappreciated, and species distributions have been misestimated as a result.

Molecular and morphological data reveal non-monophyly and speciation in imperiled freshwater mussels (Anodontoides and Strophitus)

Accurate taxonomic placement is vital to conservation efforts considering many intrinsic biological characteristics of understudied species are inferred from closely related taxa. The rayed creekshell, Anodontoides radiatus (Conrad, 1834), exists in the Gulf of Mexico drainages from western Florida to Louisiana and has been petitioned for listing under the Endangered Species Act. We set out to resolve the evolutionary history of A. radiatus, primarily generic placement and species boundaries, using phylogenetic, morphometric, and geographic information. Our molecular matrix contained 3 loci: cytochrome c oxidase subunit I, NADH dehydrogenase subunit I, and the nuclear-encoded ribosomal internal transcribed spacer I. We employed maximum likelihood and Bayesian inference to estimate a phylogeny and test the monophyly of Anodontoides and Strophitus. We implemented two coalescent-based species delimitation models to test seven species models and evaluate species boundaries within A. radiatus. Concomitant to molecular data, we also employed linear morphometrics and geographic information to further evaluate species boundaries. Molecular and morphological evidence supports the inclusion of A. radiatus in the genus Strophitus, and we resurrect the binomial Strophitus radiatus to reflect their shared common ancestry. We also found strong support for polyphyly in Strophitus and advocate the resurrection of the genus Pseudodontoideus to represent 'Strophitus' connasaugaensis and 'Strophitus' subvexus. Strophitus radiatus exists in six well-supported clades that were distinguished as evolutionary independent lineages using Bayesian inference, maximum likelihood, and coalescent-based species delimitation models. Our integrative approach found evidence for as many as 4 evolutionary divergent clades within S. radiatus. Therefore, we formally describe two new species from the S. radiatus species complex (Strophitus williamsi and Strophitus pascagoulaensis) and recognize the potential for a third putative species (Strophitus sp. cf. pascagoulaensis). Our findings aid stakeholders in establishing conservation and management strategies for the members of Anodontoides, Strophitus, and Pseudodontoideus.

"Fix Me Another Marguerite!": Species delimitation in a group of intensively hybridizing lineages of ox-eye daisies (Leucanthemum Mill., Compositae-Anthemideae)

Delineating species boundaries in the framework of the multi-species coalescent (MSC) proves to be a reliable, objective, and reproducible method in an increasing number of studies. However, the underlying model assumes the lack of gene flow after speciation; an assumption which may be frequently violated in plant evolution. This study evaluates the robustness of currently available species delimitation methods implemented in beast (BFD, BFD*, and dissect) in the closely-knit ox-eye daisy group around Leucanthemum ageratifolium Pau. Comprising five taxa being allopatrically distributed between northern Spain and southern Italy this study group shows signs of hybridization with the widespread and codistributed species Leucanthemum vulgare (Vaill.) Lam. to various extent. As expected, our empirical analyses based on both AFLP fingerprinting and sequence data demonstrate that the robustness of species delimitation results is considerably influenced by the intensity of hybridization among species and the number of hybrid individuals included. Therefore, we set up a methodological pipeline with a first step of identification and subsequent removal of individuals showing admixed genetic patterns caused by actual interbreeding using AFLP-fingerprint and morphometric data, followed by application of different Bayesian MSC species delimitation methods based on the remnant individuals using both AFLP-fingerprint and sequence data (four nuclear markers, five concatenated intergenic spacer regions of the plastid genome). The results argue for acknowledgement of Leucanthemum laciniatum, L. legraeanum, and L. ligusticum as independent species, show the close relationship of L. ageratifolium, L. monspeliense, and L. vulgare, and give rise to the description of three nothospecies new to science.

Spiders in caves

World experts of different disciplines, from molecular biology to macro-ecology, recognize the value of cave ecosystems as ideal ecological and evolutionary laboratories. Among other subterranean taxa, spiders stand out as intriguing model organisms for their ecological role of top predators, their unique adaptations to the hypogean medium and their sensitivity to anthropogenic disturbance. As the description of the first eyeless spider (Stalita taenaria), an array of papers on subterranean spider biology, ecology and evolution has been published, but a comprehensive review on these topics is still lacking. We provide a general overview of the spider families recorded in hypogean habitats worldwide, we review the different adaptations of hypogean spiders to subterranean life, and we summarize the information gathered so far about their origin, population structure, ecology and conservation status. Finally, we point out the limits of the knowledge we currently have regarding hypogean spiders, aiming to stimulate future research.

Comparative transcriptomics of Entelegyne spiders (Araneae, Entelegynae), with emphasis on molecular evolution of orphan genes

Next-generation sequencing technology is rapidly transforming the landscape of evolutionary biology, and has become a cost-effective and efficient means of collecting exome information for non-model organisms. Due to their taxonomic diversity, production of interesting venom and silk proteins, and the relative scarcity of existing genomic resources, spiders in particular are excellent targets for next-generation sequencing (NGS) methods. In this study, the transcriptomes of six entelegyne spider species from three genera (Cicurina travisae, C. vibora, Habronattus signatus, H. ustulatus, Nesticus bishopi, and N. cooperi) were sequenced and de novo assembled. Each assembly was assessed for quality and completeness and functionally annotated using gene ontology information. Approximately 100 transcripts with evidence of homology to venom proteins were discovered. After identifying more than 3,000 putatively orthologous genes across all six taxa, we used comparative analyses to identify 24 instances of positively selected genes. In addition, between ~ 550 and 1,100 unique orphan genes were found in each genus. These unique, uncharacterized genes exhibited elevated rates of amino acid substitution, potentially consistent with lineage-specific adaptive evolution. The data generated for this study represent a valuable resource for future phylogenetic and molecular evolutionary research, and our results provide new insight into the forces driving genome evolution in taxa that span the root of entelegyne spider phylogeny.

Cryptic diversity in Tranzscheliella spp. (Ustilaginales) is driven by host switches

Species of Tranzscheliella have been reported as pathogens of more than 30 genera of grasses (Poaceae). In this study, a combined morphological and molecular phylogenetic approach was used to examine 33 specimens provisionally identified as belonging to the T. hypodytes species complex. The phylogenetic analysis resolved several well-supported clades that corresponded to known and novel species of Tranzscheliella. Four new species are described and illustrated. In addition, a new combination in Tranzscheliella is proposed for Sorosporium reverdattoanum. Cophylogenetic analyses assessed by distance-based and event-cost based methods, indicated host switches are likely the prominent force driving speciation in Tranzscheliella.

Rapid dissemination of taxonomic discoveries based on DNA barcoding and morphology

The taxonomic impediment is characterized by dwindling classical taxonomic expertise, and slow pace of revisionary work, thus more rapid taxonomic assessments are needed. Here we pair rapid DNA barcoding methods with swift assessment of morphology in an effort to gauge diversity, establish species limits, and rapidly disseminate taxonomic information prior to completion of formal taxonomic revisions. We focus on a poorly studied, but diverse spider genus, Pseudopoda, from East Asia. We augmented the standard barcoding locus (COI) with nuclear DNA sequence data (ITS2) and analyzed congruence among datasets and species delimitation methods for a total of 572 individuals representing 23 described species and many potentially new species. Our results suggest that a combination of CO1 + ITS2 fragments identify and diagnose species better than the mitochondrial barcodes alone, and that certain tree based methods yield considerably higher diversity estimates than the distance-based approaches and morphology. Combined, through an extensive field survey, we detect a twofold increase in species diversity in the surveyed area, at 42–45, with most species representing short range endemics. Our study demonstrates the power of biodiversity assessments and swift dissemination of taxonomic data through rapid inventory, and through a combination of morphological and multi-locus DNA barcoding diagnoses of diverse arthropod lineages.

Catalogue of Texas spiders

This catalogue lists 1,084 species of spiders (three identified to genus only) in 311 genera from 53 families currently recorded from Texas and is based on the "Bibliography of Texas Spiders" published by Bea Vogel in 1970. The online list of species can be found at http://pecanspiders.tamu.edu/spidersoftexas.htm. Many taxonomic revisions have since been published, particularly in the families Araneidae, Gnaphosidae and Leptonetidae. Many genera in other families have been revised. The Anyphaenidae, Ctenidae, Hahniidae, Nesticidae, Sicariidae and Tetragnathidae were also revised. Several families have been added and others split up. Several genera of Corinnidae were transferred to Phrurolithidae and Trachelidae. Two genera from Miturgidae were transferred to Eutichuridae. Zoridae was synonymized under Miturgidae. A single species formerly in Amaurobiidae is now in the Family Amphinectidae. Some trapdoor spiders in the family Ctenizidae have been transferred to Euctenizidae. Gertsch and Mulaik started a list of Texas spiders in 1940. In a letter from Willis J. Gertsch dated October 20, 1982, he stated "Years ago a first listing of the Texas fauna was published by me based largely on Stanley Mulaik material, but it had to be abandoned because of other tasks." This paper is a compendium of the spiders of Texas with distribution, habitat, collecting method and other data available from revisions and collections. This includes many records and unpublished data (including data from three unpublished studies). One of these studies included 16,000 adult spiders belonging to 177 species in 29 families. All specimens in that study were measured and results are in the appendix. Hidalgo County has 340 species recorded with Brazos County at 323 and Travis County at 314 species. These reflect the amount of collecting in the area.

Taxonomist's Nightmare … Evolutionist's Delight : An Integrative Approach Resolves Species Limits in Jumping Bristletails Despite Widespread Hybridization and Parthenogenesis

Accurate species delimitation is fundamental to biology. Traditionally, species were delimited based on morphological characters, sometimes leading to taxonomic uncertainty in morphologically conserved taxa. Recently, multiple taxonomically challenging cases have benefited from integrative taxonomy-an approach that highlights congruence among different disciplines and invokes evolutionary explanations for incongruence, acknowledging that different methods can mirror different stages of the speciation continuum. Here, we used a cohesive protocol for integrative taxonomy to revise species limits in 20 nominal species and 4 morphospecies of an ancestrally wingless insect group, the jumping bristletail genus Machilis from the European Eastern Alps. Even though morphologically conserved, several small-scale endemic species have been described from the Eastern Alps based on variation in hypodermal pigmentation patterns-a highly questionable character. As valuable as these endemics are for conservation, they have never been verified by alternative methods. Using traditional morphometrics, mitochondrial DNA, ribosomal DNA, and amplified fragment-length polymorphism markers, we identify six nominal species as taxonomic junior synonyms (Machilis alpicola Janetschek, 1953 syn. n. under M. vagans Wygodzinsky, 1941; M. ladensis Janetschek, 1950 syn. n., M. robusta Wygodzinsky, 1941 syn. n., and M. vicina Wygodzinsky, 1941 syn. n. under M. inermis Wygodzinsky, 1941; M. aleamaculata Wygodzinsky, 1941 syn. n. under M. montana Wygodzinsky, 1941; M. pulchra Janetschek, 1950 syn. n. under M. helleri Verhoeff, 1910) and describe two new species (Machilis cryptoglacialis sp. n. and Machilis albida sp. n.), one uncovered from morphological crypsis and one never sampled before. Building on numerous cases of incongruence among data sources, we further shed light on complex evolutionary histories including hybrid speciation, historical and recent hybridization, and ongoing speciation. We hypothesize that an inherent affinity to hybridization, combined with parallel switches to parthenogenesis and repeated postglacial colonization events may have boosted endemicity in Eastern Alpine Machilis We thus emphasize the importance of integrative taxonomy for rigorous species delimitation and its implication for evolutionary research and conservation in taxonomically challenging taxa.

Understanding diploid diversity: A first step in unraveling polyploid, apomictic complexity in Amelanchier

PREMISE OF THE STUDY

Delimitation of Amelanchier species is difficult because of polyploidy and gametophytic apomixis. A first step in unraveling this species problem is understanding the diversity of the diploids that contributed genomes to polyploid apomicts. This research helps clarify challenging species-delimitation problems attending polyploid, apomictic complexity.

METHODS

We sampled 431 diploid accessions from 13 species, of which 10 are North American and three are Old World. Quantitative morphological analyses tested the null hypothesis of no discrete groups. Using three to nine diploid accessions per species, we constructed phylogenies with DNA sequences from ETS, ITS, the second intron of LEAFY, and chloroplast regions rpoB-trnC, rpl16, trnD-trnT, and ycf6-psbM.

KEY RESULTS

Most Amelanchier diploid taxa are morphologically and ecogeographically distinct and genetically exclusive lineages. They rarely hybridize with one another. Nuclear and chloroplast DNA sequences almost completely resolve the Amelanchier phylogeny. The backbone is the mostly western North American clade A, eastern North American clade B, and Old World clade O. DNA sequences and morphology support clades A and O as sister taxa. Despite extensive paralogy, our LEAFY data are phylogenetically informative and identify a clade (T) of three arborescent taxa within clade B.

CONCLUSIONS

Amelanchier diploids differ strikingly from polyploid apomicts, in that hybridization among them is rare, and they form taxa that would qualify as species by most species concepts. Knowledge of diploid morphology, phylogeny, and ecogeography provides a foundation for understanding the evolutionary history of polyploid apomicts, their patterns of diversification, and their species status.