Pleistocene persistence and expansion in tarantulas on the Colorado Plateau and the effects of missing data on phylogeographical inferences from RADseq.

Montane species endemic to the "sky islands" of the North American southwest were significantly impacted by changing climates during the Pleistocene. We combined mitochondrial and genomic data with species distribution modelling to determine whether Aphonopelma marxi, a large tarantula from the nearby Colorado Plateau, was similarly impacted by glacial climates. Genetic analyses revealed that the species comprises three main clades that diverged in the Pleistocene. A clade distributed along the Mogollon Rim appears to have persisted in place during glacial conditions, whereas the other two clades probably colonized central and northeastern portions of the species' range from refugia in canyons. Climate models support this hypothesis for the Mogollon Rim, but late glacial climate data appear too coarse to detect suitable areas in canyons. Locations of canyon refugia could not be inferred from genomic analyses due to missing data, encouraging us to explore the effect of missing loci in phylogeographical inferences using RADseq. Results from analyses with varying amounts of missing data suggest that samples with large amounts of missing data can still improve inferences, and the specific loci that are missing matters more than the number of missing loci. This study highlights the profound impact of Pleistocene climates on tarantulas endemic to the Colorado Plateau, as well as the mixed nature of the region's fauna. Some animals recently colonized from nearby deserts as glacial climates receded, whereas others, like tarantulas, appear to have persisted on the Mogollon Rim and in refugia associated with the region's famous river-cut canyons.

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex.

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure, principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large-scale ecological data for niche-based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic-scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.

A new species of Aphonopelma (Araneae: Mygalomorphae: Theraphosidae) from the Madrean pine-oak woodlands of northeastern Sonora, Mexico.

The tarantula spider genus Aphonopelma Pocock, 1901 has received considerable attention in recent years but the group's diversity remains poorly understood in Mexico, particularly in the pine-oak woodlands of the Sierra Madre Occidental and associated Madrean "Sky Islands". A pair of tarantulas discovered from an unsampled region in the Sierra de Bacadéhuachi (the westernmost range of the Sierra Madre Occidental) in northeastern Sonora was found to be closely related to four species from the Madrean "Sky Islands" in Arizona and New Mexico. An integrative approach for delimiting species (incorporating data from molecular phylogenetics, morphology, distributions, and breeding periods) suggests that the specimens from Sierra de Bacadéhuachi belong to an undescribed species that is herein named Aphonopelma bacadehuachi sp. nov. This new species adds to our knowledge of an increasingly diverse assemblage of Aphonopelma from the Madrean Pine-Oak Woodlands Hotspot. Collaborations between Mexican and American researchers are needed to accelerate discovery and description of the group's remaining diversity, particularly in light of the many threats facing the ecoregion including habitat degradation and climate change.

Taxonomic revision of the tarantula genus Aphonopelma Pocock, 1901 (Araneae, Mygalomorphae, Theraphosidae) within the United States.

This systematic study documents the taxonomy, diversity, and distribution of the tarantula spider genus Aphonopelma Pocock, 1901 within the United States. By employing phylogenomic, morphological, and geospatial data, we evaluated all 55 nominal species in the United States to examine the evolutionary history of Aphonopelma and the group's taxonomy by implementing an integrative approach to species delimitation. Based on our analyses, we now recognize only 29 distinct species in the United States. We propose 33 new synonymies (Aphonopelma apacheum, Aphonopelma minchi, Aphonopelma rothi, Aphonopelma schmidti, Aphonopelma stahnkei = Aphonopelma chalcodes; Aphonopelma arnoldi = Aphonopelma armada; Aphonopelma behlei, Aphonopelma vogelae = Aphonopelma marxi; Aphonopelma breenei = Aphonopelma anax; Aphonopelma chambersi, Aphonopelma clarum, Aphonopelma cryptethum, Aphonopelma sandersoni, Aphonopelma sullivani = Aphonopelma eutylenum; Aphonopelma clarki, Aphonopelma coloradanum, Aphonopelma echinum, Aphonopelma gurleyi, Aphonopelma harlingenum, Aphonopelma odelli, Aphonopelma waconum, Aphonopelma wichitanum = Aphonopelma hentzi; Aphonopelma heterops = Aphonopelma moderatum; Aphonopelma jungi, Aphonopelma punzoi = Aphonopelma vorhiesi; Aphonopelma brunnius, Aphonopelma chamberlini, Aphonopelma iviei, Aphonopelma lithodomum, Aphonopelma smithi, Aphonopelma zionis = Aphonopelma iodius; Aphonopelma phanum, Aphonopelma reversum = Aphonopelma steindachneri), 14 new species (Aphonopelma atomicum sp. n., Aphonopelma catalina sp. n., Aphonopelma chiricahua sp. n., Aphonopelma icenoglei sp. n., Aphonopelma johnnycashi sp. n., Aphonopelma madera sp. n., Aphonopelma mareki sp. n., Aphonopelma moellendorfi sp. n., Aphonopelma parvum sp. n., Aphonopelma peloncillo sp. n., Aphonopelma prenticei sp. n., Aphonopelma saguaro sp. n., Aphonopelma superstitionense sp. n., and Aphonopelma xwalxwal sp. n.), and seven nomina dubia (Aphonopelma baergi, Aphonopelma cratium, Aphonopelma hollyi, Aphonopelma mordax, Aphonopelma radinum, Aphonopelma rusticum, Aphonopelma texense). Our proposed species tree based on Anchored Enrichment data delimits five major lineages: a monotypic group confined to California, a western group, an eastern group, a group primarily distributed in high-elevation areas, and a group that comprises several miniaturized species. Multiple species are distributed throughout two biodiversity hotspots in the United States (i.e., California Floristic Province and Madrean Pine-Oak Woodlands). Keys are provided for identification of both males and females. By conducting the most comprehensive sampling of a single theraphosid genus to date, this research significantly broadens the scope of prior molecular and morphological investigations, finally bringing a modern understanding of species delimitation in this dynamic and charismatic group of spiders.

An evaluation of sampling effects on multiple DNA barcoding methods leads to an integrative approach for delimiting species: a case study of the North American tarantula genus Aphonopelma (Araneae, Mygalomorphae, Theraphosidae).

The North American tarantula genus Aphonopelma provides one of the greatest challenges to species delimitation and downstream identification in spiders because traditional morphological characters appear ineffective for evaluating limits of intra- and interspecific variation in the group. We evaluated the efficacy of numerous molecular-based approaches to species delimitation within Aphonopelma based upon the most extensive sampling of theraphosids to date, while also investigating the sensitivity of randomized taxon sampling on the reproducibility of species boundaries. Mitochondrial DNA (cytochrome c oxidase subunit I) sequences were sampled from 682 specimens spanning the genetic, taxonomic, and geographic breadth of the genus within the United States. The effects of random taxon sampling compared traditional Neighbor-Joining with three modern quantitative species delimitation approaches (ABGD, P ID(Liberal), and GMYC). Our findings reveal remarkable consistency and congruence across various approaches and sampling regimes, while highlighting highly divergent outcomes in GMYC. Our investigation allowed us to integrate methodologies into an efficient, consistent, and more effective general methodological workflow for estimating species boundaries within the mygalomorph spider genus Aphonopelma. Taken alone, these approaches are not particularly useful - especially in the absence of prior knowledge of the focal taxa. Only through the incorporation of multiple lines of evidence, employed in a hypothesis-testing framework, can the identification and delimitation of confident species boundaries be determined. A key point in studying closely related species, and perhaps one of the most important aspects of DNA barcoding, is to combine a sampling strategy that broadly identifies the extent of genetic diversity across the distributions of the species of interest and incorporates previous knowledge into the "species equation" (morphology, molecules, and natural history).

An exploration of species boundaries in turret-building tarantulas of the Mojave Desert (Araneae, Mygalomorphae, Theraphosidae, Aphonopelma).

Tarantulas in the North American genus Aphonopelma are poorly known due to their challenging patterns of morphological variation and questionable taxonomy; few specimens can be confidently identified using existing keys or comparisons to original descriptions. In an effort to identify new strategies for resolving what has been characterized as a "taxonomic and nomenclatural nightmare", we employed five different approaches for delimiting species in a group of closely related tarantulas from the Mojave Desert in the southwestern United States. These methods included the application of single techniques (morphology, DNA barcoding, shared genealogical exclusivity among independent loci, and generalized mixed Yule coalescent) and an integrative approach that incorporates genealogical and ecological information. Results demonstrate that the taxonomy of these spiders as presently defined underestimates actual species-level diversity and the group is in need of revision. The number of species delimited by each approach, however, was variable and we argue that it is this discordance that emphasizes the importance of incorporating multiple lines of evidence into an integrative taxonomic framework that can be used for constructing robust taxonomic hypotheses for Aphonopelma species.

A reconsideration of the classification of the spider infraorder Mygalomorphae (Arachnida: Araneae) based on three nuclear genes and morphology.

BACKGROUND: The infraorder Mygalomorphae (i.e., trapdoor spiders, tarantulas, funnel web spiders, etc.) is one of three main lineages of spiders. Comprising 15 families, 325 genera, and over 2,600 species, the group is a diverse assemblage that has retained a number of features considered primitive for spiders. Despite an evolutionary history dating back to the lower Triassic, the group has received comparatively little attention with respect to its phylogeny and higher classification. The few phylogenies published all share the common thread that a stable classification scheme for the group remains unresolved. METHODS AND FINDINGS: We report here a reevaluation of mygalomorph phylogeny using the rRNA genes 18S and 28S, the nuclear protein-coding gene EF-1γ, and a morphological character matrix. Taxon sampling includes members of all 15 families representing 58 genera. The following results are supported in our phylogenetic analyses of the data: (1) the Atypoidea (i.e., antrodiaetids, atypids, and mecicobothriids) is a monophyletic group sister to all other mygalomorphs; and (2) the families Mecicobothriidae, Hexathelidae, Cyrtaucheniidae, Nemesiidae, Ctenizidae, and Dipluridae are not monophyletic. The Microstigmatidae is likely to be subsumed into Nemesiidae. Nearly half of all mygalomorph families require reevaluation of generic composition and placement. The polyphyletic family Cyrtaucheniidae is most problematic, representing no fewer than four unrelated lineages. CONCLUSIONS: Based on these analyses we propose the following nomenclatural changes: (1) the establishment of the family Euctenizidae (NEW RANK); (2) establishment of the subfamily Apomastinae within the Euctenizidae; and (3) the transfer of the cyrtaucheniid genus Kiama to Nemesiidae. Additional changes include relimitation of Domiothelina and Theraphosoidea, and the establishment of the Euctenizoidina clade (Idiopidae + Euctenizidae). In addition to these changes, we propose a "road map" for future sampling across the infraorder with the aim of solving many remaining questions that hinder mygalomorph systematics.

Phylogeny and classification of the trapdoor spider genus Myrmekiaphila: an integrative approach to evaluating taxonomic hypotheses.

BACKGROUND: Revised by Bond and Platnick in 2007, the trapdoor spider genus Myrmekiaphila comprises 11 species. Species delimitation and placement within one of three species groups was based on modifications of the male copulatory device. Because a phylogeny of the group was not available these species groups might not represent monophyletic lineages; species definitions likewise were untested hypotheses. The purpose of this study is to reconstruct the phylogeny of Myrmekiaphila species using molecular data to formally test the delimitation of species and species-groups. We seek to refine a set of established systematic hypotheses by integrating across molecular and morphological data sets. METHODS AND FINDINGS: Phylogenetic analyses comprising Bayesian searches were conducted for a mtDNA matrix composed of contiguous 12S rRNA, tRNA-val, and 16S rRNA genes and a nuclear DNA matrix comprising the glutamyl and prolyl tRNA synthetase gene each consisting of 1348 and 481 bp, respectively. Separate analyses of the mitochondrial and nuclear genome data and a concatenated data set yield M. torreya and M. millerae paraphyletic with respect to M. coreyi and M. howelli and polyphyletic fluviatilis and foliata species groups. CONCLUSIONS: Despite the perception that molecular data present a solution to a crisis in taxonomy, studies like this demonstrate the efficacy of an approach that considers data from multiple sources. A DNA barcoding approach during the species discovery process would fail to recognize at least two species (M. coreyi and M. howelli) whereas a combined approach more accurately assesses species diversity and illuminates speciation pattern and process. Concomitantly these data also demonstrate that morphological characters likewise fail in their ability to recover monophyletic species groups and result in an unnatural classification. Optimizations of these characters demonstrate a pattern of "Dollo evolution" wherein a complex character evolves only once but is lost multiple times throughout the group's history.

Evaluating the efficacy of continuous quantitative characters for reconstructing the phylogeny of a morphologically homogeneous spider taxon (Araneae, Mygalomorphae, Antrodiaetidae, Antrodiaetus).

The use of continuous quantitative characters for phylogenetic analyses has long been contentious in the systematics literature. Recent studies argue for and against their use, but there have been relatively few attempts to evaluate whether these characters provide an accurate estimate of phylogeny, despite the fact that a number of methods have been developed to analyze these types of data for phylogenetic inference. A tree topology will be produced for a given methodology and set of characters, but little can be concluded with regards to the accuracy of phylogenetic signal without an independent evaluation of those characters. We assess the performance of continuous quantitative characters for the mygalomorph spider genus Antrodiaetus, a group that is morphologically homogeneous and one for which few discrete (morphological) characters have been observed. Phylogenetic signal contained in continuous quantitative characters is compared to an independently derived phylogeny inferred on the basis of multiple nuclear and mitochondrial gene loci. Tree topology randomizations, regression techniques, and topological tests all demonstrate that continuous quantitative characters in Antrodiaetus conflict with the phylogenetic signal contained in the gene trees. Our results show that the use of continuous quantitative characters for phylogenetic reconstruction may be inappropriate for reconstructing Antrodiaetus phylogeny and indicate that due caution should be exercised before employing this character type in the absence of other independently derived sources of characters.

Molecular phylogeny and biogeography of an ancient Holarctic lineage of mygalomorph spiders (Araneae: Antrodiaetidae: Antrodiaetus).

The mygalomorph spider genera Antrodiaetus and Atypoides (Antrodiaetidae) belong to an ancient lineage that has persisted since at least the Cretaceous. These spiders display a classic disjunct Holarctic distribution with species in the eastern Palaearctic plus the western and eastern Nearctic. Prior phylogenetic analyses of this group have been proposed on the basis of morphology, but lack strong support and independent corroboration. Here we present the first phylogenetic analysis of species-level relationships based on molecular data obtained from the mitochondrial (cytochrome c oxidase subunit I) and nuclear (18S and 28S rRNA) genomes. Analyses corroborate earlier findings that Atypoides forms a paraphyletic grade with respect to Antrodiaetus, and consequently, that genus is formally synonymized under Antrodiaetus. In addition, our results support the relatively early divergence of Antrodiaetus roretzi. Antrodiaetus pacificus is "paraphyletic" with respect to the A. lincolnianus group and is likely an assemblage of numerous species. The final topology based on a combined molecular dataset, in conjunction with two different molecular dating techniques (penalized likelihood plus a Bayesian approach) and ancestral distribution reconstructions, was used to infer the historical biogeography of these spiders. Trans-Beringian and trans-Atlantic routes appear to account for the present-day distribution of Antrodiaetus in Japan and North America. Future studies on Antrodiaetus phylogeny will be used to address questions regarding morphological stasis and the evolution of quantitative morphological characters.

Testing species boundaries in the Antrodiaetus unicolor complex (Araneae: Mygalomorphae: Antrodiaetidae): "paraphyly" and cryptic diversity.

The inability to correctly identify species has far reaching implications in nearly all areas of biology, yet few studies investigate methods for delineating species boundaries. Moreover, once these boundaries have been hypothesized, little thought has been given to how these constructs can be further evaluated. We employ a molecular phylogenetic approach using nuclear 28S rRNA and mitochondrial cytochrome c oxidase subunit I genes to test the general efficacy of species boundaries in the Antrodiaetus unicolor spider species complex. Our analyses provide evidence that An. unicolor is "paraphyletic" with respect to An. microunicolor, indicating that morphological criteria used to delineate species boundaries undersplits actual species-level diversity in this group of spiders. These analyses also demonstrate that individuals from geographically proximate populations sometimes exhibit considerable molecular divergence, strongly suggesting that An. unicolor is a cryptic species complex. Finally, this molecular approach has provided the phylogenetic framework that is necessary to begin interpreting the vast amount of morphological variation observed in these spiders based upon findings from previous studies. Our approach using multiple genes appears to be a rigorous method to critically examine species boundaries originally based on traditional morphological approaches to spider taxonomy.