Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida, Opiliones, Sclerobunus robustus ) from the southwestern United States

Species delimitation, biogeography, and natural history of dwarf funnel web spiders (Mygalomorphae, Hexurellidae, Hexurella) from the United States / Mexico borderlands

The rarely encountered spider genus Hexurella Gertsch & Platnick, 1979 includes some of the smallest mygalomorph spiders in the world, with four poorly known taxa from central and southeastern montane Arizona, southern California, and northern Baja California Norte. At time of description the genus was known from fewer than 20 individuals, with sparse natural history information suggesting a vagrant, web-building, litter-dwelling natural history. Here the first published taxonomic and natural history information for this taxon is provided in more than 50 years, working from extensive new geographic sampling, consideration of male and female morphology, and sequence capture-based nuclear phylogenomics and mitogenomics. Several new species are easily diagnosed based on distinctive male morphologies, while a complex of populations from central and northern Arizona required an integrative combination of genomic algorithmic species delimitation analyses and morphological study. Four new species are described, including H.ephedrasp. nov., H.uwiiltilsp. nov., H.xericasp. nov., and H.zassp. nov. Females of H.encina Gertsch & Platnick, 1979 are also described for the first time. It is predicted that additional new species will ultimately be found in the mountains of central and northwestern Arizona, northern mainland Mexico, and the Mojave Desert of California.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

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

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

Holocene Population Decline and Conservation Implication for the Western Hercules Beetle, Dynastes grantii (Coleoptera, Scarabaeidae)

The Western Hercules beetle (Dynastes grantii) is endemic to the highland forest habitats of southwestern United States and northern Mexico. The habitats harbor many endemic species, but are being threatened by rapid climate change and urban development. In this study, the genetic structure of D. grantii populations from southwestern United States was investigated. Specifically, genomic data from double-digest restriction-site-associated DNA sequencing libraries were utilized to test whether geographically distant populations from the Mogollon Rim (Arizona [N = 12 individuals] and New Mexico [N = 10 individuals]) are genetically structured. The study also estimated the effective population size of the Mogollon Rim populations based on genetic diversity. The results indicated that the 2 geographic populations from the Mogollon Rim were not genetically structured. A population size reduction was detected since the end of the last glacial period, which coincided with a reduction of forest habitat in the study area. The results implied that the connectivity and the size of highland forest habitats in the Mogollon Rim could have been the major factors shaping the population genetic structure and demographic history of D. grantii. The Western Hercules beetle could be a useful flagship species for local natural history education and to promote the conservation of highland forest habitats.

Reference-based RADseq resolves robust relationships among closely related species of lichen-forming fungi using metagenomic DNA

Despite increasing availability of phylogenomic datasets, strategies to generate genome-scale data from organisms involved in symbiotic relationships remains challenging. Restriction site-associated DNA sequencing (RADseq) can effectively generated reduced representation genomic loci. However, when using metagenomic DNA from inseparable symbiotic organisms, RADseq loci may belong to any number of the organisms involved in these intimate associations. In this study, we explored the potential for a reference-based RADseq approach to generate data for lichen-forming fungi from metagenomic DNA extracted from intact lichens. We simulated RAD data from draft genomes of closely related lichenized fungi to test if RADseq can reconstruct robust evolutionary relationships. Subsequently, we generated empirical RADseq data from metagenomic lichen DNA, with RADseq loci mapped back to a reference genome to exclude loci from other lichen symbionts that are represented in metagenomic libraries. In all cases, phylogenetic reconstructions using RADseq loci recovered diversification histories consistent with a previous study based on more comprehensive genome sampling. Furthermore, RADseq loci were found to resolve relationships among closely related species, which were otherwise indistinguishable using a phylogenetic species recognition criterion. Our studies revealed that a modified, reference-based RADseq approach can successfully be implemented to generate symbiont-specific phylogenomic data from metagenomic reads.

High phylogenetic utility of an ultraconserved element probe set designed for Arachnida

Arachnida is an ancient, diverse and ecologically important animal group that contains a number of species of interest for medical, agricultural and engineering applications. Despite their importance, many aspects of the arachnid tree of life remain unresolved, hindering comparative approaches to arachnid biology. Biologists have made considerable efforts to resolve the arachnid phylogeny; yet, limited and challenging morphological characters, as well as a dearth of genetic resources, have hindered progress. Here, we present a genomic toolkit for arachnids featuring hundreds of conserved DNA regions (ultraconserved elements or UCEs) that allow targeted sequencing of any species in the arachnid tree of life. We used recently developed capture probes designed from conserved regions of available arachnid genomes to enrich a sample of loci from 32 diverse arachnids. Sequence capture returned an average of 487 UCE loci for all species, with a range from 170 to 722. Phylogenetic analysis of these UCEs produced a highly resolved arachnid tree with relationships largely consistent with recent transcriptome-based phylogenies. We also tested the phylogenetic informativeness of UCE probes within the spider, scorpion and harvestman orders, demonstrating the utility of these markers at shallower taxonomic scales and suggesting that these loci will be useful for species-level differences. This probe set will open the door to phylogenomic and population genomic studies across the arachnid tree of life, enabling systematics, species delimitation, species discovery and conservation of these diverse arthropods.

New species in the Sitalcina sura species group (Opiliones, Laniatores, Phalangodidae), with evidence for a biogeographic link between California desert canyons and Arizona sky islands

The western United States is home to numerous narrowly endemic harvestman taxa (Arachnida, Opiliones), including members of the genus Sitalcina Banks, 1911. Sitalcina is comprised of three species groups, including the monospecific Sitalcina californica and Sitalcina lobata groups, and the Sitalcina sura group with eight described species. All species in the Sitalcina sura group have very small geographic distributions, with group members distributed like disjunct "beads on a string" from Monterey south to southern California and southeast to the sky-island mountain ranges of southern Arizona. Here, molecular phylogenetic and species delimitation analyses were conducted for all described species in the Sitalcina sura group, plus several newly discovered populations. Species trees were reconstructed using multispecies coalescent methods implemented in *BEAST, and species delimitation was accomplished using Bayes Factor Delimitation (BFD). Based on quantitative species delimitation results supported by consideration of morphological characters, two new species (Sitalcina oasiensis sp. n., Sitalcina ubicki sp. n.) are described. We also provide a description of the previously unknown male of Sitalcina borregoensis Briggs, 1968. Molecular phylogenetic evidence strongly supports distinctive desert versus coastal clades, with desert canyon taxa from southern California more closely related to Arizona taxa than to geographically proximate California coastal taxa. We hypothesize that southern ancestry and plate tectonics have played a role in the diversification history of this animal lineage, similar to sclerophyllous plant taxa of the Madro-Tertiary Geoflora. Molecular clock analyses for the Sitalcina sura group are generally consistent with these hypotheses. We also propose that additional Sitalcina species await discovery in the desert canyons of southern California and northern Baja, and the mountains of northwestern mainland Mexico.