Jason E. Bond

Phylogenomics Resolves a Spider Backbone Phylogeny and Rejects a Prevailing Paradigm for Orb Web Evolution

Spiders represent an ancient predatory lineage known for their extraordinary biomaterials, including venoms and silks. These adaptations make spiders key arthropod predators in most terrestrial ecosystems. Despite ecological, biomedical, and biomaterial importance, relationships among major spider lineages remain unresolved or poorly supported. Current working hypotheses for a spider backbone phylogeny are largely based on morphological evidence, as most molecular markers currently employed are generally inadequate for resolving deeper-level relationships. We present here a phylogenomic analysis of spiders including taxa representing all major spider lineages. Our robust phylogenetic hypothesis recovers some fundamental and uncontroversial spider clades, but rejects the prevailing paradigm of a monophyletic Orbiculariae, the most diverse lineage, containing orb-weaving spiders. Based on our results, the orb web either evolved much earlier than previously hypothesized and is ancestral for a majority of spiders or else it has multiple independent origins, as hypothesized by precladistic authors. Cribellate deinopoid orb weavers that use mechanically adhesive silk are more closely related to a diverse clade of mostly webless spiders than to the araneoid orb-weaving spiders that use adhesive droplet silks. The fundamental shift in our understanding of spider phylogeny proposed here has broad implications for interpreting the evolution of spiders, their remarkable biomaterials, and a key extended phenotype--the spider web.

Spider phylogenomics: untangling the Spider Tree of Life

Spiders (Order Araneae) are massively abundant generalist arthropod predators that are found in nearly every ecosystem on the planet and have persisted for over 380 million years. Spiders have long served as evolutionary models for studying complex mating and web spinning behaviors, key innovation and adaptive radiation hypotheses, and have been inspiration for important theories like sexual selection by female choice. Unfortunately, past major attempts to reconstruct spider phylogeny typically employing the “usual suspect” genes have been unable to produce a well-supported phylogenetic framework for the entire order. To further resolve spider evolutionary relationships we have assembled a transcriptome-based data set comprising 70 ingroup spider taxa. Using maximum likelihood and shortcut coalescence-based approaches, we analyze eight data sets, the largest of which contains 3,398 gene regions and 696,652 amino acid sites forming the largest phylogenomic analysis of spider relationships produced to date. Contrary to long held beliefs that the orb web is the crowning achievement of spider evolution, ancestral state reconstructions of web type support a phylogenetically ancient origin of the orb web, and diversification analyses show that the mostly ground-dwelling, web-less RTA clade diversified faster than orb weavers. Consistent with molecular dating estimates we report herein, this may reflect a major increase in biomass of non-flying insects during the Cretaceous Terrestrial Revolution 125–90 million years ago favoring diversification of spiders that feed on cursorial rather than flying prey. Our results also have major implications for our understanding of spider systematics. Phylogenomic analyses corroborate several well-accepted high level groupings: Opisthothele, Mygalomorphae, Atypoidina, Avicularoidea, Theraphosoidina, Araneomorphae, Entelegynae, Araneoidea, the RTA clade, Dionycha and the Lycosoidea. Alternatively, our results challenge the monophyly of Eresoidea, Orbiculariae, and Deinopoidea. The composition of the major paleocribellate and neocribellate clades, the basal divisions of Araneomorphae, appear to be falsified. Traditional Haplogynae is in need of revision, as our findings appear to support the newly conceived concept of Synspermiata. The sister pairing of filistatids with hypochilids implies that some peculiar features of each family may in fact be synapomorphic for the pair. Leptonetids now are seen as a possible sister group to the Entelegynae, illustrating possible intermediates in the evolution of the more complex entelegyne genitalic condition, spinning organs and respiratory organs.

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).

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.

Species Delimitation and Phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): Cryptic Diversity in North American Tarantulas

Background The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal “barcoding gene”. An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame. Methods and Findings A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a “barcode gap” using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups. Conclusions The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a “species”.

The spider tree of life: phylogeny of Araneae based on target‐gene analyses from an extensive taxon sampling

We present a phylogenetic analysis of spiders using a dataset of 932 spider species, representing 115 families (only the family Synaphridae is unrepresented), 700 known genera, and additional representatives of 26 unidentified or undescribed genera. Eleven genera of the orders Amblypygi, Palpigradi, Schizomida and Uropygi are included as outgroups. The dataset includes six markers from the mitochondrial (12S, 16S, COI) and nuclear (histone H3, 18S, 28S) genomes, and was analysed by multiple methods, including constrained analyses using a highly supported backbone tree from transcriptomic data. We recover most of the higher‐level structure of the spider tree with good support, including Mesothelae, Opisthothelae, Mygalomorphae and Araneomorphae. Several of our analyses recover Hypochilidae and Filistatidae as sister groups, as suggested by previous transcriptomic analyses. The Synspermiata are robustly supported, and the families Trogloraptoridae and Caponiidae are found as sister to the Dysderoidea. Our results support the Lost Tracheae clade, including Pholcidae, Tetrablemmidae, Diguetidae, Plectreuridae and the family Pacullidae (restored status) separate from Tetrablemmidae. The Scytodoidea include Ochyroceratidae along with Sicariidae, Scytodidae, Drymusidae and Periegopidae; our results are inconclusive about the separation of these last two families. We did not recover monophyletic Austrochiloidea and Leptonetidae, but our data suggest that both groups are more closely related to the Cylindrical Gland Spigot clade rather than to Synspermiata. Palpimanoidea is not recovered by our analyses, but also not strongly contradicted. We find support for Entelegynae and Oecobioidea (Oecobiidae plus Hersiliidae), and ambiguous placement of cribellate orb‐weavers, compatible with their non‐monophyly. Nicodamoidea (Nicodamidae plus Megadictynidae) and Araneoidea composition and relationships are consistent with recent analyses. We did not obtain resolution for the titanoecoids (Titanoecidae and Phyxelididae), but the Retrolateral Tibial Apophysis clade is well supported. Penestomidae, and probably Homalonychidae, are part of Zodarioidea, although the latter family was set apart by recent transcriptomic analyses. Our data support a large group that we call the marronoid clade (including the families Amaurobiidae, Desidae, Dictynidae, Hahniidae, Stiphidiidae, Agelenidae and Toxopidae). The circumscription of most marronoid families is redefined here. Amaurobiidae include the Amaurobiinae and provisionally Macrobuninae. We transfer Malenellinae (Malenella, from Anyphaenidae), Chummidae (Chumma) (new syn.) and Tasmarubriinae (Tasmarubrius, Tasmabrochus and Teeatta, from Amphinectidae) to Macrobuninae. Cybaeidae are redefined to include Calymmaria, Cryphoeca, Ethobuella and Willisius (transferred from Hahniidae), and Blabomma and Yorima (transferred from Dictynidae). Cycloctenidae are redefined to include Orepukia (transferred from Agelenidae) and Pakeha and Paravoca (transferred from Amaurobiidae). Desidae are redefined to include five subfamilies: Amphinectinae, with Amphinecta, Mamoea, Maniho, Paramamoea and Rangitata (transferred from Amphinectidae); Ischaleinae, with Bakala and Manjala (transferred from Amaurobiidae) and Ischalea (transferred from Stiphidiidae); Metaltellinae, with Austmusia, Buyina, Calacadia, Cunnawarra, Jalkaraburra, Keera, Magua, Metaltella, Penaoola and Quemusia; Porteriinae (new rank), with Baiami, Cambridgea, Corasoides and Nanocambridgea (transferred from Stiphidiidae); and Desinae, with Desis, and provisionally Poaka (transferred from Amaurobiidae) and Barahna (transferred from Stiphidiidae). Argyroneta is transferred from Cybaeidae to Dictynidae. Cicurina is transferred from Dictynidae to Hahniidae. The genera Neoramia (from Agelenidae) and Aorangia, Marplesia and Neolana (from Amphinectidae) are transferred to Stiphidiidae. The family Toxopidae (restored status) includes two subfamilies: Myroinae, with Gasparia, Gohia, Hulua, Neomyro, Myro, Ommatauxesis and Otagoa (transferred from Desidae); and Toxopinae, with Midgee and Jamara, formerly Midgeeinae, new syn. (transferred from Amaurobiidae) and Hapona, Laestrygones, Lamina, Toxops and Toxopsoides (transferred from Desidae). We obtain a monophyletic Oval Calamistrum clade and Dionycha; Sparassidae, however, are not dionychans, but probably the sister group of those two clades. The composition of the Oval Calamistrum clade is confirmed (including Zoropsidae, Udubidae, Ctenidae, Oxyopidae, Senoculidae, Pisauridae, Trechaleidae, Lycosidae, Psechridae and Thomisidae), affirming previous findings on the uncertain relationships of the “ctenids” Ancylometes and Cupiennius, although a core group of Ctenidae are well supported. Our data were ambiguous as to the monophyly of Oxyopidae. In Dionycha, we found a first split of core Prodidomidae, excluding the Australian Molycriinae, which fall distantly from core prodidomids, among gnaphosoids. The rest of the dionychans form two main groups, Dionycha part A and part B. The former includes much of the Oblique Median Tapetum clade (Trochanteriidae, Gnaphosidae, Gallieniellidae, Phrurolithidae, Trachelidae, Gnaphosidae, Ammoxenidae, Lamponidae and the Molycriinae), and also Anyphaenidae and Clubionidae. Orthobula is transferred from Phrurolithidae to Trachelidae. Our data did not allow for complete resolution for the gnaphosoid families. Dionycha part B includes the families Salticidae, Eutichuridae, Miturgidae, Philodromidae, Viridasiidae, Selenopidae, Corinnidae and Xenoctenidae (new fam., including Xenoctenus, Paravulsor and Odo, transferred from Miturgidae, as well as Incasoctenus from Ctenidae). We confirm the inclusion of Zora (formerly Zoridae) within Miturgidae.

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.

Expanding anchored hybrid enrichment to resolve both deep and shallow relationships within the spider tree of life.

BackgroundDespite considerable effort, progress in spider molecular systematics has lagged behind many other comparable arthropod groups, thereby hindering family-level resolution, classification, and testing of important macroevolutionary hypotheses. Recently, alternative targeted sequence capture techniques have provided molecular systematics a powerful tool for resolving relationships across the Tree of Life. One of these approaches, Anchored Hybrid Enrichment (AHE), is designed to recover hundreds of unique orthologous loci from across the genome, for resolving both shallow and deep-scale evolutionary relationships within non-model systems. Herein we present a modification of the AHE approach that expands its use for application in spiders, with a particular emphasis on the infraorder Mygalomorphae.ResultsOur aim was to design a set of probes that effectively capture loci informative at a diversity of phylogenetic timescales. Following identification of putative arthropod-wide loci, we utilized homologous transcriptome sequences from 17 species across all spiders to identify exon boundaries. Conserved regions with variable flanking regions were then sought across the tick genome, three published araneomorph spider genomes, and raw genomic reads of two mygalomorph taxa. Following development of the 585 target loci in the Spider Probe Kit, we applied AHE across three taxonomic depths to evaluate performance: deep-level spider family relationships (33 taxa, 327 loci); family and generic relationships within the mygalomorph family Euctenizidae (25 taxa, 403 loci); and species relationships in the North American tarantula genus Aphonopelma (83 taxa, 581 loci). At the deepest level, all three major spider lineages (the Mesothelae, Mygalomorphae, and Araneomorphae) were supported with high bootstrap support. Strong support was also found throughout the Euctenizidae, including generic relationships within the family and species relationships within the genus Aptostichus. As in the Euctenizidae, virtually identical topologies were inferred with high support throughout Aphonopelma.ConclusionsThe Spider Probe Kit, the first implementation of AHE methodology in Class Arachnida, holds great promise for gathering the types and quantities of molecular data needed to accelerate an understanding of the spider Tree of Life by providing a mechanism whereby different researchers can confidently and effectively use the same loci for independent projects, yet allowing synthesis of data across independent research groups.

Millipede Taxonomy after 250 Years: Classification and Taxonomic Practices in a Mega-Diverse yet Understudied Arthropod Group

Background The arthropod class Diplopoda is a mega-diverse group comprising >12,000 described millipede species. The history of taxonomic research within the group is tumultuous and, consequently, has yielded a questionable higher-level classification. Few higher-taxa are defined using synapomorphies, and the practice of single taxon descriptions lacking a revisionary framework has produced many monotypic taxa. Additionally, taxonomic and geographic biases render global species diversity estimations unreliable. We test whether the ordinal taxa of the Diplopoda are consistent with regards to underlying taxonomic diversity, attempt to provide estimates for global species diversity, and examine millipede taxonomic effort at a global geographic scale. Methodology/Principal Findings A taxonomic distinctness metric was employed to assess uniformity of millipede ordinal taxa. We found that ordinal-level taxa are not uniform and are likely overinflated with higher-taxa when compared to related groups. Several methods of estimating global species richness were employed (Bayesian, variation in taxonomic productivity, extrapolation from nearly fully described taxa). Two of the three methods provided estimates ranging from 13,413–16,760 species. Variations in geographic diversity show biases to North America and Europe and a paucity of works on tropical taxa. Conclusions/Significance Before taxa can be used in an extensible way, they must be definable with respect to the diversity they contain and the diagnostic characters used to delineate them. The higher classification for millipedes is shown to be problematic from a number of perspectives. Namely, the ordinal taxa are not uniform in their underlying diversity, and millipedes appear to have a disproportionate number of higher-taxa. Species diversity estimates are unreliable due to inconsistent taxonomic effort at temporal, geographic, and phylogenetic scales. Lack of knowledge concerning many millipede groups compounds these issues. Diplopods are likely not unique in this regard as these issues may persist in many other diverse yet poorly studied groups.

Ordinal-Level Phylogenomics of the Arthropod Class Diplopoda (Millipedes) Based on an Analysis of 221 Nuclear Protein-Coding Loci Generated Using Next-Generation Sequence Analyses

Background The ancient and diverse, yet understudied arthropod class Diplopoda, the millipedes, has a muddled taxonomic history. Despite having a cosmopolitan distribution and a number of unique and interesting characteristics, the group has received relatively little attention; interest in millipede systematics is low compared to taxa of comparable diversity. The existing classification of the group comprises 16 orders. Past attempts to reconstruct millipede phylogenies have suffered from a paucity of characters and included too few taxa to confidently resolve relationships and make formal nomenclatural changes. Herein, we reconstruct an ordinal-level phylogeny for the class Diplopoda using the largest character set ever assembled for the group. Methods Transcriptomic sequences were obtained from exemplar taxa representing much of the diversity of millipede orders using second-generation (i.e., next-generation or high-throughput) sequencing. These data were subject to rigorous orthology selection and phylogenetic dataset optimization and then used to reconstruct phylogenies employing Bayesian inference and maximum likelihood optimality criteria. Ancestral reconstructions of sperm transfer appendage development (gonopods), presence of lateral defense secretion pores (ozopores), and presence of spinnerets were considered. The timings of major millipede lineage divergence points were estimated. Results The resulting phylogeny differed from the existing classifications in a number of fundamental ways. Our phylogeny includes a grouping that has never been described (Juliformia+Merocheta+Stemmiulida), and the ancestral reconstructions suggest caution with respect to using spinnerets as a unifying characteristic for the Nematophora. Our results are shown to have significantly stronger support than previous hypotheses given our data. Our efforts represent the first step toward obtaining a well-supported and robust phylogeny of the Diplopoda that can be used to answer many questions concerning the evolution of this ancient and diverse animal group.