Comparative phylogeography reveals deep lineages and regional evolutionary hotspots in the M ojave and S onoran Deserts

Filling knowledge gaps in insect conservation by leveraging genetic data from public archives

Insect decline has become a growing concern in recent years, with studies showing alarming declines in populations of several taxa. Our knowledge about genetic spatial patterns and evolutionary history of insects still exhibits significant gaps hindering our ability to effectively conserve and manage insect populations and species. Genetic data may provide valuable insights into the diversity and the evolutionary relationships of insects' species and populations. Public repositories, such as GenBank and BOLD, containing vast archives of genetic data with associated metadata, offer an irreplaceable resource for researchers contributing to our understanding of species diversity, population structure and evolutionary relationships. However, there are some issues in using these data, as they are often scattered and may lack accuracy due to inconsistent sampling protocols and incomplete information. In this paper we describe a curated georeferenced database of genetic data collected in GenBank and BOLD, for insects listed in the International Union for Conservation of Nature (IUCN) Italian Red Lists (dragonflies, bees, saproxylic beetles and butterflies). After querying these repositories, we performed quality control and data standardization steps. We created a dataset containing approximately 33 000 mitochondrial sequences and associated metadata about taxonomy, collection localities, geographic coordinates and IUCN Red List status for 1466 species across the four insect lists. We describe the current state of geographical metadata in queried repositories for species listed under different conservation status in the Italian Red Lists to quantify data gaps posing barriers to prioritization of conservation actions. Our curated dataset is available for data repurposing and analysis, enabling researchers to conduct comparative studies. We emphasize the importance of filling knowledge gaps in insect diversity and distribution and highlight the potential of this dataset for promoting other research fields like phylogeography, macrogenetics and conservation strategies. Our database can be downloaded through the Zenodo repository in SQL format. Database URL:  https://zenodo.org/records/8375181.

Diversification and post-glacial range expansion of giant North American camel spiders in genus Eremocosta (Solifugae: Eremobatidae)

Species of camel spiders in the family Eremobatidae are an important component of arthropod communities in arid ecosystems throughout North America. Recently, research demonstrated that the evolutionary history and biogeography of the family are poorly understood. Herein we explore the biogeographic history of this group of arachnids using genome-wide single nucleotide polymorphism (SNP) data, morphology, and distribution modelling to study the eremobatid genus Eremocosta, which contains exceptionally large species distributed throughout North American deserts. Relationships among sampled species were resolved with strong support and they appear to have diversified within distinct desert regions along an east-to-west progression beginning in the Chihuahuan Desert. The unexpected phylogenetic position of some samples suggests that the genus may contain additional, morphologically cryptic species. Geometric morphometric analyses reveal a largely conserved cheliceral morphology among Eremocosta spp. Phylogeographic analyses indicate that the distribution of E. titania was substantially reduced during the last glacial maximum and the species only recently colonized much of the Mojave Desert. Results from this study underscore the power of genome-wide data for unlocking the genetic potential of museum specimens, which is especially promising for organisms like camel spiders that are notoriously difficult to collect.

Genetic landscape clustering of a large DNA barcoding data set reveals shared patterns of genetic divergence among freshwater fishes of the Maroni Basin

The Maroni is one of the most speciose basins of the Guianas and hosts a megadiverse freshwater fish community. Although taxonomic references based on morphological identification exist for both the Surinamese and Guianese parts of the basin, there are still taxonomic uncertainties concerning the status of several species. We used COI sequences of 1284 fish in conjunction with morphological and biogeographical evidence to assist with species delineation and discovery in order to validate and standardize the current taxonomy. This resulted in a final DNA barcode data set of 199 fish species (125 genera, 36 families and eight orders; 68.86% of strictly freshwater fishes from the basin), among which 25 are new putative candidate species flagged as requiring taxonomic update. DNA barcoding delineation through Barcode Index Numbers (BINs) revealed further cryptic diversity (230 BINs in total). To explore global genetic patterns across the basin, genetic divergence landscapes were computed for 128 species, showing a global trend of high genetic divergence between the Surinamese southwest (Tapanahony and Paloemeu), the Guianese southeast (Marouini, Litany, Tampok, etc.), and the river outlet in the north. This could be explained by lower levels of connectivity between these three main areas and/or the exchange of individuals between these areas and the neighbouring basins. A new method of ordination of genetic landscapes successfully assigned species into cluster groups based on their respective pattern of genetic divergence across the Maroni Basin: genetically homogeneous species were effectively discriminated from species showing high spatial genetic fragmentation and possible lower capacity for dispersal.

Comparative phylogeography clarifies the complexity and problems of continental distribution that drove A. R. Wallace to favor islands

Deciphering the geographic context of diversification and distributional dynamics in continental biotas has long been an interest of biogeographers, ecologists, and evolutionary biologists. Thirty years ago, the approach now known as comparative phylogeography was introduced in a landmark study of a continental biota. Here, I use a set of 455 studies to explore the current scope of continental comparative phylogeography, including geographic, conceptual, temporal, ecological, and genomic attributes. Geographically, studies are more frequent in the northern hemisphere, but the south is catching up. Most studies focus on a Quaternary timeframe, but the Neogene is well represented. As such, explanations for geographic structure and history include geological and climatic events in Earth history, and responses include vicariance, dispersal, and range contraction-expansion into and out of refugia. Focal taxa are biased toward terrestrial or semiterrestrial vertebrates, although plants and invertebrates are well represented in some regions. The use of various kinds of nuclear DNA markers is increasing, as are multiple locus studies, but use of organelle DNA is not decreasing. Species distribution models are not yet widely incorporated into studies. In the future, continental comparative phylogeographers will continue to contribute to erosion of the simple vicariance vs. dispersal paradigm, including exposure of the widespread nature of temporal pseudocongruence and its implications for models of diversification; provide new templates for addressing a variety of ecological and evolutionary traits; and develop closer working relationships with earth scientists and biologists in a variety of disciplines.

multi-dice: r package for comparative population genomic inference under hierarchical co-demographic models of independent single-population size changes

Population genetic data from multiple taxa can address comparative phylogeographic questions about community-scale response to environmental shifts, and a useful strategy to this end is to employ hierarchical co-demographic models that directly test multi-taxa hypotheses within a single, unified analysis. This approach has been applied to classical phylogeographic data sets such as mitochondrial barcodes as well as reduced-genome polymorphism data sets that can yield 10,000s of SNPs, produced by emergent technologies such as RAD-seq and GBS. A strategy for the latter had been accomplished by adapting the site frequency spectrum to a novel summarization of population genomic data across multiple taxa called the aggregate site frequency spectrum (aSFS), which potentially can be deployed under various inferential frameworks including approximate Bayesian computation, random forest and composite likelihood optimization. Here, we introduce the r package multi-dice, a wrapper program that exploits existing simulation software for flexible execution of hierarchical model-based inference using the aSFS, which is derived from reduced genome data, as well as mitochondrial data. We validate several novel software features such as applying alternative inferential frameworks, enforcing a minimal threshold of time surrounding co-demographic pulses and specifying flexible hyperprior distributions. In sum, multi-dice provides comparative analysis within the familiar R environment while allowing a high degree of user customization, and will thus serve as a tool for comparative phylogeography and population genomics.

Phylogeography of a widespread small carnivore, the western spotted skunk (Spilogale gracilis) reveals temporally variable signatures of isolation across western North America

We analyzed phylogeographic patterns in the western spotted skunk, Spilogale gracilis Merriam, 1890 (Carnivora: Mephitidae) in relation to historical events associated with Pre-Pleistocene Divergence (PPD) and Quaternary climate change (QCC) using mitochondrial DNA from 97 individuals distributed across Western North America. Divergence times were generated using BEAST to estimate when isolation in putative refugia occurred. Patterns and timing of demographic expansion was performed using Bayesian skyline plot. Putative climatic refugia resulting from Quaternary climate change were identified using paleoecological niche modeling and divergence dates compared to major vicariant events associated with Pre-Pleistocene conditions. We recovered three major mitochondrial clades corresponding to western North America (California, Baja, and across the Great Basin), east-central North America (Texas, central Mexico, New Mexico), and southwestern Arizona/northwestern Mexico. Time to most recent common ancestor for S. gracilis occurred ~1.36 Ma. Divergence times for each major clade occurred between 0.25 and 0.12 Ma, with signature of population expansion occurring 0.15 and 0.10 Ma. Ecological niche models identified three potential climatic refugia during the Last Interglacial, (1) west coast of California and Oregon, (2) northwestern Mexico, and (3) southern Texas/northeastern Mexico as well as two refugia during the Last Glacial Maximum, (1) western USA and (2) southern Texas/northeastern Mexico. This study supports PPD in shaping species-level diversity compared to QCC-driven changes at the intraspecific level for Spilogale, similar to the patterns reported for other small mammals (e.g., rodents and bats). Phylogeographic patterns also appear to have been shaped by both habitat and river vicariance, especially across the desert southwest. Further, continuing climate change during the Holocene coupled with anthropogenic modifications during the Anthropocene appears to be removing both of these barriers to current dispersal of western spotted skunks.

High gene flow in the American badger overrides habitat preferences and limits broadscale genetic structure

Habitat associations are a function of habitat preferences and dispersal capabilities, both of which can influence how species responded to Quaternary climatic changes and contemporary habitat heterogeneity. Predicting resultant genetic structure is not always straightforward, especially in species where high dispersal potential and habitat preferences yield opposing predictions. The American badger has high dispersal capabilities that predict widespread panmixia, but avoids closed-canopy forests and clay soils, which could restrict gene flow and create ecologically based population genetic structure. We used mitochondrial sequence and microsatellite data sets to characterize how these opposing forces contribute to genetic structure in badgers at a continent-wide scale. Our data revealed an overall lack of ecologically based population genetic structure, suggesting that high dispersal capabilities were sufficiently realized to overcome most habitat-based genetic structure. At a broadscale, badger gene flow is limited only by geographic distance (isolation by distance) and large water barriers (Lake Michigan and the Mississippi River). The absence of genetic structure in a species with strong avoidance of unsuitable habitats advances our understanding of when and how genetic structure emerges in widespread, highly mobile species.

Testing Taxon Tenacity of Tortoises: evidence for a geographical selection gradient at a secondary contact zone

We examined a secondary contact zone between two species of desert tortoise, Gopherus agassizii and G. morafkai. The taxa were isolated from a common ancestor during the formation of the Colorado River (4–8 mya) and are a classic example of allopatric speciation. However, an anomalous population of G. agassizii comes into secondary contact with G. morafkai east of the Colorado River in the Black Mountains of Arizona and provides an opportunity to examine reinforcement of species' boundaries under natural conditions. We sampled 234 tortoises representing G. agassizii in California (n - 103), G. morafkai in Arizona (n - 78), and 53 individuals of undetermined assignment in the contact zone including and surrounding the Black Mountains. We genotyped individuals for 25 STR loci and determined maternal lineage using mtDNA sequence data. We performed multilocus genetic clustering analyses and used multiple statistical methods to detect levels of hybridization. We tested hypotheses about habitat use between G. agassizii and G. morafkai in the region where they co-occur using habitat suitability models. Gopherus agassizii and G. morafkai maintain independent taxonomic identities likely due to ecological niche partitioning, and the maintenance of the hybrid zone is best described by a geographical selection gradient model.

Detecting concerted demographic response across community assemblages using hierarchical approximate Bayesian computation

Methods that integrate population-level sampling from multiple taxa into a single community-level analysis are an essential addition to the comparative phylogeographic toolkit. Detecting how species within communities have demographically tracked each other in space and time is important for understanding the effects of future climate and landscape changes and the resulting acceleration of extinctions, biological invasions, and potential surges in adaptive evolution. Here, we present a statistical framework for such an analysis based on hierarchical approximate Bayesian computation (hABC) with the goal of detecting concerted demographic histories across an ecological assemblage. Our method combines population genetic data sets from multiple taxa into a single analysis to estimate: 1) the proportion of a community sample that demographically expanded in a temporally clustered pulse and 2) when the pulse occurred. To validate the accuracy and utility of this new approach, we use simulation cross-validation experiments and subsequently analyze an empirical data set of 32 avian populations from Australia that are hypothesized to have expanded from smaller refugia populations in the late Pleistocene. The method can accommodate data set heterogeneity such as variability in effective population size, mutation rates, and sample sizes across species and exploits the statistical strength from the simultaneous analysis of multiple species. This hABC framework used in a multitaxa demographic context can increase our understanding of the impact of historical climate change by determining what proportion of the community responded in concert or independently and can be used with a wide variety of comparative phylogeographic data sets as biota-wide DNA barcoding data sets accumulate.

Speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a species of conservation concern

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation-with-migration (IM) models and a novel coalescent-based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.

Fuzzy boundaries: color and gene flow patterns among parapatric lineages of the western shovel-nosed snake and taxonomic implication

Accurate delineation of lineage diversity is increasingly important, as species distributions are becoming more reduced and threatened. During the last century, the subspecies category was often used to denote phenotypic variation within a species range and to provide a framework for understanding lineage differentiation, often considered incipient speciation. While this category has largely fallen into disuse, previously recognized subspecies often serve as important units for conservation policy and management when other information is lacking. In this study, we evaluated phenotypic subspecies hypotheses within shovel-nosed snakes on the basis of genetic data and considered how evolutionary processes such as gene flow influenced possible incongruence between phenotypic and genetic patterns. We used both traditional phylogenetic and Bayesian clustering analyses to infer range-wide genetic structure and spatially explicit analyses to detect possible boundary locations of lineage contact. Multilocus analyses supported three historically isolated groups with low to moderate levels of contemporary gene exchange. Genetic data did not support phenotypic subspecies as exclusive groups, and we detected patterns of discordance in areas where three subspecies are presumed to be in contact. Based on genetic and phenotypic evidence, we suggested that species-level diversity is underestimated in this group and we proposed that two species be recognized, Chionactis occipitalis and C. annulata. In addition, we recommend retention of two subspecific designations within C. annulata (C. a. annulata and C. a. klauberi) that reflect regional shifts in both genetic and phenotypic variation within the species. Our results highlight the difficultly in validating taxonomic boundaries within lineages that are evolving under a time-dependent, continuous process.

Cryptic genetic diversity and complex phylogeography of the boreal North American scorpion, Paruroctonus boreus (Vaejovidae)

Diverse studies in western North America have revealed the role of topography for dynamically shaping genetic diversity within species though vicariance, dispersal and range expansion. We examined patterns of phylogeographical diversity in the widespread but poorly studied North American vaejovid scorpion, Paruroctonus boreus Girard 1854. We used mitochondrial sequence data and parsimony, likelihood, and Bayesian inference to reconstruct phylogenetic relationships across the distributional range of P. boreus, focusing on intermontane western North America. Additionally, we developed a species distribution model to predict its present and historical distributions during the Last Glacial Maximum and the Last Interglacial Maximum. Our results documented complex phylogeographic relationships within P. boreus, with multiple, well-supported crown clades that are either geographically-circumscribed or widespread and separated by short, poorly supported internodes. We also observed subtle variation in predicted habitat suitability, especially at the northern, eastern and southern edges of the predicted distributional range under past climatic conditions. The complex phylogenetic relationships of P. boreus suggests that historical isolation and expansion of populations may have occurred. Variation in the predicted distributional range over time may implicate past climatic fluctuations in generating the patterns of genetic diversity observed in P. boreus. These findings highlight both the potential for cryptic biodiversity in widespread North American scorpion species and the importance of phylogeographical studies for understanding the factors responsible for generating the biodiversity of western North America.