Genealogical exclusivity in geographically proximate populations of Hypochilus thorelli Marx (Araneae, Hypochilidae) on the Cumberland Plateau of North America

Comparative Population Genomic Diversity and Differentiation in Trapdoor Spiders and Relatives (Araneae, Mygalomorphae)

Although patterns of population genomic variation are well-studied in animals, there remains room for studies that focus on non-model taxa with unique biologies. Here we characterise and attempt to explain such patterns in mygalomorph spiders, which are generally sedentary, often occur as spatially clustered demes and show remarkable longevity. Genome-wide single nucleotide polymorphism (SNP) data were collected for 500 individuals across a phylogenetically representative sample of taxa. We inferred genetic populations within focal taxa using a phylogenetically informed clustering approach, and characterised patterns of diversity and differentiation within- and among these genetic populations, respectively. Using phylogenetic comparative methods we asked whether geographical range sizes and ecomorphological variables (behavioural niche and body size) significantly explain patterns of diversity and differentiation. Specifically, we predicted higher genetic diversity in genetic populations with larger geographical ranges, and in small-bodied taxa. We also predicted greater genetic differentiation in small-bodied taxa, and in burrowing taxa. We recovered several significant predictors of genetic diversity, but not genetic differentiation. However, we found generally high differentiation across genetic populations for all focal taxa, and a consistent signal for isolation-by-distance irrespective of behavioural niche or body size. We hypothesise that high population genetic structuring, likely reflecting combined dispersal limitation and microhabitat specificity, is a shared trait for all mygalomorphs. Few studies have found ubiquitous genetic structuring for an entire ancient and species-rich animal clade.

Mitochondrial DNA and Population Genomics Reveal Additional Cryptic Diversity in the Green Salamander (Subgenus Castaneides) Species Complex

Cryptic species present particular challenges to biodiversity conservation, as true species diversity and distributional boundaries remain obscured. However, modern molecular tools have afforded unparalleled opportunities to elucidate cryptic species, define their distributions, and, ultimately, develop conservation interventions to extend their evolutionary trajectories into the future. The Green Salamander (Aneides aeneus) complex provides an evolutionary focal point and the Appalachian Highlands an ecological context for the exploration of cryptic speciation in an imperiled taxon. A recent study uncovered significant levels of genetic and genomic variation geographically structured across the Appalachian Highlands, including up to four lineages, one of which (A. caryaensis) was described therein. Here we extend the genetic and genomic examination of the Castaneides species complex by intensive sampling of additional populations along Cumberland Plateau and Appalachian Valley and Ridge of Alabama and Tennessee, employing both mtDNA and RADseq species delimitation approaches to delineate cryptic diversity and boundaries in this region. Analyses of two mitochondrial loci (nd4 and cytb) identified two reciprocally monophyletic lineages, which are also supported by population clustering and phylogenetic analyses of SNPs, that identified two population clusters with no evidence of gene flow. Our genetic and genomic results support the recognition of two additional cryptic lineages in the Castaneides species complex. Ultimately, this information is critical in developing successful adaptive management strategies for this important and endemic component of Appalachian Highland biodiversity.

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.

The role of glacial-interglacial climate change in shaping the genetic structure of eastern subterranean termites in the southern Appalachian Mountains, USA

The eastern subterranean termite, Reticulitermes flavipes, currently inhabits previously glaciated regions of the northeastern U.S., as well as the unglaciated southern Appalachian Mountains and surrounding areas. We hypothesized that Pleistocene climatic fluctuations have influenced the distribution of R. flavipes, and thus the evolutionary history of the species. We estimated contemporary and historical geographic distributions of R. flavipes by constructing Species Distribution Models (SDM). We also inferred the evolutionary and demographic history of the species using mitochondrial (cytochrome oxidase I and II) and nuclear (endo-beta-1,4-glucanase) DNA sequence data. To do this, genetic populations were delineated using Bayesian spatial-genetic clustering, competing hypotheses about population divergence were assessed using approximate Bayesian computation (ABC), and changes in population size were estimated using Bayesian skyline plots. SDMs identified areas in the north with suitable habitat during the transition from the Last Interglacial to the Last Glacial Maximum, as well as an expanding distribution from the mid-Holocene to the present. Genetic analyses identified three geographically cohesive populations, corresponding with northern, central, and southern portions of the study region. Based on ABC analyses, divergence between the Northern and Southern populations was the oldest, estimated to have occurred 64.80 thousand years ago (kya), which corresponds with the timing of available habitat in the north. The Central and Northern populations diverged in the mid-Holocene, 8.63 kya, after which the Central population continued to expand. Accordingly, phylogeographic patterns of R. flavipes in the southern Appalachians appear to have been strongly influenced by glacial-interglacial climate change.

OPEN RESEARCH BADGES

This article has been awarded Open Materials, Open Data Badges. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.5061/dryad.5hr7f31.

Desert salt flats as oases for the spider Saltonia incerta Banks (Araneae: Dictynidae)

The deserts of southwestern North America have undergone dramatic changes over their recent geological history including large changes in size and connectivity during the Pleistocene glaciopluvial cycles. This study examines the population history of the rare spider Saltonia incerta, once thought to be extinct, to determine the role of past climatological events in shaping the structure of the species. This species is restricted to salt crusts of intermittent or dry lakes, streams or rivers in the desert southwest, a region that was much wetter during glacial periods. We examine the distribution and genetic variability of populations to test whether there is recent dispersal throughout the range of the species. Analyses of mitochondrial and nuclear DNA indicate significant population structure, with one major clade comprising New Mexico localities and one comprising California-northern Baja California localities. Finer-scale structure is evident within the California clade, although not all of the subclades are reciprocally monophyletic. However, isolation with migration analysis suggests that migration is very low to non-existent. These results extend the known distribution of Saltonia, provide genetic evidence of strong isolation among localities within drainage basins and between drainage basins and provide a mechanistic understanding of population connectivity after the aridification of the American southwest. The implication is that although the species' distribution has been fragmented, populations have persisted throughout this area, suggesting that desert salt flats may have served as refugia for at least some terrestrial species.

Accelerated species inventory on Madagascar using coalescent-based models of species delineation

High-throughput DNA sequencing has the potential to accelerate species discovery if it is able to recognize evolutionary entities from sequence data that are comparable to species. The general mixed Yule-coalescent (GMYC) model estimates the species boundary from DNA surveys by identifying independently evolving lineages as a transition from coalescent to speciation branching patterns on a phylogenetic tree. Applied here to 12 families from 4 orders of insects in Madagascar, we used the model to delineate 370 putative species from mitochondrial DNA sequence variation among 1614 individuals. These were compared with data from the nuclear genome and morphological identification and found to be highly congruent (98% and 94%). We developed a modified GMYC that allows for a variable transition from coalescent to speciation among lineages. This revised model increased the congruence with morphology (97%), suggesting that a variable threshold better reflects the clustering of sequence data into biological species. Local endemism was pronounced in all 5 insect groups. Most species (60-91%) and haplotypes (88-99%) were found at only 1 of the 5 study sites (40-1000 km apart). This pronounced endemism resulted in a 37% increase in species numbers using diagnostic nucleotides in a population aggregation analysis. Sample sizes between 7 and 10 individuals represented a threshold above which there was minimal increase in genetic diversity, broadly agreeing with coalescent theory and other empirical studies. Our results from > 1.4 Mb of empirical data suggest that the GMYC model captures species boundaries comparable to those from traditional methods without the need for prior hypotheses of population coherence. This provides a method of species discovery and biodiversity assessment using single-locus data from mixed or environmental samples while building a globally available taxonomic database for future identifications.

Species status and conservation issues of New Zealand's endemic Latrodectus spider species (Araneae : Theridiidae)

New Zealand has two endemic widow spiders, Latrodectus katipo Powell, 1871 and L. atritus Urquhart, 1890. Both species face many conservation threats and are actively managed. The species status of the Latrodectus spiders of New Zealand was assessed using molecular (COI, ITS1, ITS2) and morphological methods and with cross-breeding experiments. Latrodectus katipo and L. atritus were not found to be reciprocally monophyletic for any of the gene regions or morphological traits. Other than colour, which is variable, there were no morphological characters that separated the two species, which cross-bred in the laboratory and produced fertile eggsacs. Colour variation is clinal over latitude and correlates significantly with mean annual temperature. We conclude that L. atritus is a junior synonym of L. katipo. An example of introgression from the Australian species L. hasseltii Thorell, 1870 was also detected and its conservation implications are discussed.

Delimiting cohesion species: extreme population structuring and the role of ecological interchangeability

Species exhibiting morphological homogeneity and strong population structuring present challenging taxonomic problems: morphology-based approaches infer few species, whereas genetic approaches often indicate more. Morphologically cryptic, yet genetically divergent species groups require alternative approaches to delimiting species that assess adaptive divergence and ecological interchangeability of lineages. We apply such an approach to Promyrmekiaphila, a small genus (three nominal taxa) of trapdoor spiders endemic to northern California to define cohesion species (lineages that are genetically exchangeable and ecologically interchangeable). Genetic exchangeability is evaluated using standard phylogeographical techniques (e.g. nested clade analysis); ecological interchangeability is assessed using two GIS-based approaches. First, climatic values are extracted from layer data for each locality point and utilized in a principal components analysis followed by MANOVA. Second, niche-based distribution models of genetically divergent lineages are created using a maximum-entropy modelling approach; the amount of overlap among lineages is calculated and evaluated against a probability distribution of null overlap. Lineages that have significant amounts of predicted overlap are considered ecologically interchangeable. Based on a synthetic evaluation of ecological interchangeability, geographical concordance, and morphological differentiation, we conclude that Promyrmekiaphila comprises six cohesion species, five of which are cryptic (i.e. undetectable by conventional means).

Multi-locus DNA sequence data reveal a history of deep cryptic vicariance and habitat-driven convergence in the desert night lizardXantusia vigilisspecies complex (Squamata: Xantusiidae)

The lizard genus Xantusia of southwestern North America has received recent attention in relation to delimiting species. Using more than 500 lizards from 156 localities, we further test hypothesized species boundaries and clarify phylogeographical patterns, particularly in regions of potential secondary contact. We sequenced the entire mitochondrial cytochrome b gene for every lizard in the study, plus a second mitochondrial DNA (mtDNA) region and two nuclear introns for subsets of the total sample. Phylogenetic analyses of the mtDNA recover a well-resolved, novel hypothesis for species in the Xantusia vigilis complex. The nuclear DNA (nDNA) data provide independent support for the recognition of X. arizonae, X. bezyi and X. wigginsi. Differences between the respective mtDNA and nDNA topologies result from either the effects of lineage sorting or ancient introgression. Nuclear data confirm the inference that some populations of X. vigilis in northwestern Arizona converged on rock-crevice-dwelling morphology and are not X. arizonae with an introgressed X. vigilis mtDNA genome. The historical independence of ancient cryptic lineages of Xantusia in southern California is also corroborated, though limited introgression is detected. Our proposed biogeographical scenario indicates that diversification of this group was driven by vicariance beginning in the late Miocene. Additionally, Pleistocene climatical changes influenced Xantusia distribution, and the now inhospitable Colorado Desert previously supported night lizard presence. The current taxonomy of the group likely underestimates species diversity within the group, and our results collectively show that while convergence on the rock-crevice-dwelling morphology is one hallmark of Xantusia evolution, morphological stasis is paradoxically another.

The age and evolution of sociality in Stegodyphus spiders: a molecular phylogenetic perspective

Social, cooperative breeding behaviour is rare in spiders and generally characterized by inbreeding, skewed sex ratios and high rates of colony turnover, processes that when combined may reduce genetic variation and lower individual fitness quickly. On these grounds, social spider species have been suggested to be unstable in evolutionary time, and hence sociality a rare phenomenon in spiders. Based on a partial molecular phylogeny of the genus Stegodyphus, we address the hypothesis that social spiders in this genus are evolutionary transient. We estimate the age of the three social species, test whether they represent an ancestral or derived state and assess diversification relative to subsocial congeners. Intraspecific sequence divergence was high in all of the social species, lending no support for the idea that they are young, transient species. The age of the social lineages, constant lineage branching and the likelihood that social species are independently derived suggest that either the social species are 'caught in sociality' or they have evolved into cryptic species.

The Philodromus pulchellus-group in the Mediterranean: taxonomic revision, phylogenetic analysis and biogeography (Araneae:Philodromidae)

The Philodromus pulchellus species-group is defined and diagnosed. Eleven species are included, described or redescribed, keyed and illustrated: P. afroglaucinus, sp. nov. from Algeria; P. bistigma Simon, 1870, P. glaucinus Simon, 1870, P. lamellipalpis, sp. nov. from Algeria; P. medius O. P.-Cambridge, 1872; P. pardalis, sp. nov. from northern Africa and the Iberian peninsula; P. pulchellus Lucas, 1846, P. punctigerus O. P.-Cambridge, 1908, P. ruficapillus Simon, 1885, P. simoni de Mello-Leitão, 1929, and P. wunderlichi, sp. nov. from the western Canary Islands. The validity of P. bistigma and P. medius is re-established (formerly in synonymy with P. pulchellus); neotype and lectotype, respectively are newly designated. The following new synonymies are proposed: P. torquatus O. P.- Cambridge, 1908 = P. pulchellus; P. salinarum Denis, 1939 = P. glaucinus; P. glaucinoides Wunderlich, 1987 = P. punctigerus; P. marionschmidti (Schmidt, 1990) = P. pulchellus. Philodromus albopictus Simon, 1875 and P. rubidus Simon, 1870 are considered nomina dubia. Determination of phylogenetic relationships within the group is difficult owing to continuous character variation, resulting in partially incongruent reconstructions using morphological and molecular data (partial mitochondrial cytochrome c oxidase I gene). Dispersal vicariance analysis provides support for a western Mediterranean origin of the group.

Molecular taxonomy in the dark: evolutionary history, phylogeography, and diversity of cave crayfish in the subgenus Aviticambarus, genus Cambarus

Freshwater crayfish species in the subgenus Aviticambarus (Cambaridae: Cambarus) are restricted to caves along the Cumberland Plateau, the Sequatchie Valley, and the Highland Rim which extend along the Tennessee River in southcentral Tennessee and northern Alabama. Historically, three stygobitic species, Cambarus jonesi, Cambarus hamulatus, and Cambarus veitchorum, comprise this subgenus. We examine species' boundaries and phylogeographic structure of this imperiled Southern Appalachian assemblage to shed light on patterns of cave colonization. We also provide estimates of genetic diversity for conservation status assessment. Using geologic evidence, phylogeographic analyses, and sequence data from five gene regions (two nuclear: Histone H3 and GAPDH and three mitochondrial: 12S, 16S, and CO1 totaling almost 2700 base pairs), we identify two well-supported cryptic species in addition to the three currently recognized taxa. Four of these taxa exhibit low levels of genetic variation both currently and historically, which may indicate local extirpation events associated with geological and river basin changes. Our results also support other recent findings that pre-Pleistocene paleodrainages may best explain the current patterns of aquatic faunal biodiversity in the Southern Appalachians.

Studies of morphological and molecular phylogenetic divergence in spiders (Araneae: Homalonychus) from the American southwest, including divergence along the Baja California Peninsula

Comparative phylogenetic and phylogeographic analyses have revealed a pervasive midpeninsular divergence in the mitochondrial genealogies of numerous vertebrate taxa distributed on the Baja California Peninsula. In this study, we extend the investigation of regional vicariance in Baja California to an arthropod taxon by examining patterns of phylogenetic and morphological divergence in the spider genus Homalonychus (Araneae, Homalonychidae). We analyzed data from two mtDNA genes (16S rRNA and NADH dehydrogenase subunit (1) and a nuclear gene (28S rRNA) using maximum parsimony and Bayesian phylogenetic analyses, and also conducted geometric morphometric analyses employing landmark data on male and female genitalia. Genes and morphology both reveal a deep split across the Colorado River and Gulf of California, separating Homalonychus selenopoides on the east side of river from its congener Homalonychus theologus on the west side of the river, including the Baja California Peninsula. Along the north-south axis of the Baja Peninsula, an apparently more recent midpeninsular phylogenetic break is evident within H. theologus in the mitochondrial genome and in female genitalia. However, there is no measurable divergence between northern and southern populations in either nuclear DNA or male genitalia. We suggest that this discordance between datasets reflects either a difference in rates of evolution between male versus female systems, or that male-based nuclear gene flow is obscuring a phylogenetic split that is fixed in the female-based systems. Our findings provide additional support for a midpeninsular Baja divergence event, although the timing and geological evidence for such an event remain elusive.

The power and perils of 'molecular taxonomy': a case study of eyeless and endangered Cicurina (Araneae: Dictynidae) from Texas caves

Rapid development in karst-rich regions of the US state of Texas has prompted the listing of four Cicurina species (Araneae, Dictynidae) as US Federally Endangered. A major constraint in the management of these taxa is the extreme rarity of adult specimens, which are required for accurate species identification. We report a first attempt at using mitochondrial DNA (mtDNA) sequences to accurately identify immature Cicurina specimens. This identification is founded on a phylogenetic framework that is anchored by identified adult and/or topotypic specimens. Analysis of approximately 1 kb of cytochrome oxidase subunit I (CO1) mtDNA data for over 100 samples results in a phylogenetic tree that includes a large number of distinctive, easily recognizable, tip clades. These tip clades almost always correspond to a priori species hypotheses, and show nonoverlapping patterns of sequence divergence, making it possible to place species names on a number of immature specimens. Three cases of inconsistency between recovered tip clades and a priori species hypotheses suggest possible introgression between cave-dwelling Cicurina, or alternatively, species synonymy. Although species determination is not possible in these instances, the inconsistencies point to areas of taxonomic ambiguity that require further study. Our molecular phylogenetic sample is largest for the Federally Endangered C. madla. These data suggest that C. madla occurs in more than twice the number of caves as previously reported, and indicate the possible synonymy of C. madla with C. vespera, which is also Federally Endangered. Network analyses reveal considerable genetic divergence and structuring across caves in this species. Although the use of DNA sequences to identify previously 'unidentifiable' specimens illustrates the potential power of molecular data in taxonomy, many other aspects of the same dataset speak to the necessity of a balanced taxonomic approach.

Testing Species Boundaries in an Ancient Species Complex with Deep Phylogeographic History: GenusXantusia(Squamata: Xantusiidae)

Identification of species in natural populations has recently received increased attention with a number of investigators proposing rigorous methods for species delimitation. Morphologically conservative species (or species complexes) with deep phylogenetic histories (and limited gene flow) are likely to pose particular problems when attempting to delimit species, yet this is crucial to comparative studies of the geography of speciation. We apply two methods of species delimitation to an ancient group of lizards (genus Xantusia) that occur throughout southwestern North America. Mitochondrial cytochrome b and nicotinamide adenine dinucleotide dehydrogenase subunit 4 gene sequences were generated from samples taken throughout the geographic range of Xantusia. Maximum likelihood, Bayesian, and nested cladogram analyses were used to estimate relationships among haplotypes and to infer evolutionary processes. We found multiple well-supported independent lineages within Xantusia, for which there is considerable discordance with the currently recognized taxonomy. High levels of sequence divergence (21.3%) suggest that the pattern in Xantusia may predate the vicariant events usually hypothesized for the fauna of the Baja California peninsula, and the existence of deeply divergent clades (18.8%-26.9%) elsewhere in the complex indicates the occurrence of ancient sundering events whose genetic signatures were not erased by the late Wisconsin vegetation changes. We present a revised taxonomic arrangement for this genus consistent with the distinct mtDNA lineages and discuss the phylogeographic history of this genus as a model system for studies of speciation in North American deserts.