Phylogenomics resolves a spider backbone phylogeny and rejects a prevailing paradigm for orb web evolution

Phylogenomic data reveal three new families of poorly studied Solifugae (camel spiders)

The systematics of the arachnid order Solifugae have been an enigma, owing to challenges in interpreting morphology, a paucity of molecular phylogenetic studies sampling across the group, and a dearth of taxonomic attention for many lineages. Recent work has suggested that solifuge families largely exhibit contiguous distributions and reflect patterns of vicariance, with the exception of three families: Melanoblossidae, Daesiidae and Gylippidae. Morphological studies have cast doubt on their existing circumscriptions and the present composition of these taxa renders their distributions as disjunct. We leveraged ultraconserved elements (UCEs) to test the phylogenetic placement of three key lineages of Solifugae that cause these anomalous distributions: Dinorhax rostrumpsittaci (putative melanoblossid), Namibesia (putative daesiid), and Trichotoma (putative gylippid). Phylogenetic placement of these three genera based on UCEs rendered the families that harbor them as para- or polyphyletic, recovering instead relationships that better accord with a biogeographic history driven by vicariance. Toward a stable and phylogenetically informed classification of Solifugae, we establish three new families, Dinorhaxidae new rank, Namibesiidae new rank and Lipophagidae new rank.

Chromosome-level genome of spider Pardosa pseudoannulata and cuticle protein genes in environmental stresses

Spiders are representative arthropods of adaptive radiation. The high-quality genomes have only been reported in several web weaver spider species, leaving the wandering spiders' genomic information scarce. The pond wolf spider, Pardosa pseudoannulata, is a representative species in the retrolateral titial apophysis (RTA) clade. We present a chromosome-level P. pseusoannulata genome assembly of 2.42 Gb in size with a scaffold N50 of 169.99 Mb. Hi-C scaffolding assigns 94.83% of the bases to 15 pseudo-chromosomes. The repeats account for 52.79% of the assembly. The assembly includes 96.2% of the complete arthropod universal single-copy orthologs. Gene annotation predicted 24,530 protein-coding genes with a BUSCO score of 95.8% complete. We identified duplicate clusters of Hox genes and an expanded cuticle protein gene family with 243 genes. The expression patterns of CPR genes change in response to environmental stresses such as coldness and insecticide exposure. The high-quality P. pseudoannulata genome provides valuable information for functional and comparative studies in spiders.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Advances in the reconstruction of the spider tree of life: A roadmap for spider systematics and comparative studies

In the last decade and a half, advances in genetic sequencing technologies have revolutionized systematics, transforming the field from studying morphological characters or a few genetic markers, to genomic datasets in the phylogenomic era. A plethora of molecular phylogenetic studies on many taxonomic groups have come about, converging on, or refuting prevailing morphology or legacy-marker-based hypotheses about evolutionary affinities. Spider systematics has been no exception to this transformation and the inter-relationships of several groups have now been studied using genomic data. About 51 500 extant spider species have been described, all with a conservative body plan, but innumerable morphological and behavioural peculiarities. Inferring the spider tree of life using morphological data has been a challenging task. Molecular data have corroborated many hypotheses of higher-level relationships, but also resulted in new groups that refute previous hypotheses. In this review, we discuss recent advances in the reconstruction of the spider tree of life and highlight areas where additional effort is needed with potential solutions. We base this review on the most comprehensive spider phylogeny to date, representing 131 of the 132 spider families. To achieve this sampling, we combined six Sanger-based markers with newly generated and publicly available genome-scale datasets. We find that some inferred relationships between major lineages of spiders (such as Austrochiloidea, Palpimanoidea and Synspermiata) are robust across different classes of data. However, several new hypotheses have emerged with different classes of molecular data. We identify and discuss the robust and controversial hypotheses and compile this blueprint to design future studies targeting systematic revisions of these problematic groups. We offer an evolutionary framework to explore comparative questions such as evolution of venoms, silk, webs, morphological traits and reproductive strategies.

Neglected no longer: Phylogenomic resolution of higher-level relationships in Solifugae

Advanced sequencing technologies have expedited resolution of higher-level arthropod relationships. Yet, dark branches persist, principally among groups occurring in cryptic habitats. Among chelicerates, Solifugae ("camel spiders") is the last order lacking a higher-level phylogeny and have thus been historically characterized as "neglected [arachnid] cousins". Though renowned for aggression, remarkable running speed, and xeric adaptation, inferring solifuge relationships has been hindered by inaccessibility of diagnostic morphological characters, whereas molecular investigations have been limited to one of 12 recognized families. Our phylogenomic dataset via capture of ultraconserved elements sampling all extant families recovered a well-resolved phylogeny, with two distinct groups of New World taxa nested within a broader Paleotropical radiation. Divergence times using fossil calibrations inferred that Solifugae radiated by the Permian, and most families diverged prior to the Paleogene-Cretaceous extinction, likely driven by continental breakup. We establish Boreosolifugae new suborder uniting five Laurasian families, and Australosolifugae new suborder uniting seven Gondwanan families using morphological and biogeographic signal.

Lampshade web spider Ectatosticta davidi chromosome-level genome assembly provides evidence for its phylogenetic position

The spider of Ectatosticta davidi, belonging to the lamp-shade web spider family, Hypochilidae, which is closely related to Hypochilidae and Filistatidae and recovered as sister of the rest Araneomorphs spiders. Here we show the final assembled genome of E. davidi with 2.16 Gb in 15 chromosomes. Then we confirm the evolutionary position of Hypochilidae. Moreover, we find that the GMC gene family exhibit high conservation throughout the evolution of true spiders. We also find that the MaSp genes of E. davidi may represent an early stage of MaSp and MiSp genes in other true spiders, while CrSp shares a common origin with AgSp and PySp but differ from MaSp. Altogether, this study contributes to addressing the limited availability of genomic sequences from Hypochilidae spiders, and provides a valuable resource for investigating the genomic evolution of spiders.

A molecular phylogeny of the European nesticid spiders (Nesticidae, Araneae): Implications for their systematics and biogeography

Nesticidae is a small family of spiders with a worldwide distribution that includes 15 genera and 272 described species. Seven genera and 56 species are known from Europe, distributed from the Iberian Peninsula to the Caucasus and the Ural Mountains. Most of these European species are cave dwellers and many of them are troglobites. In this study we present the first molecular phylogeny of the family Nesticidae in Europe with a wide geographical sampling across the continent. In our analysis the European nesticid fauna is well represented, including six genera and 40 of the 56 currently accepted species including the type species of all sampled genera. We have included in the analysis representatives of the North American and Asian fauna to test the monophyly of the European species and the phylogenetic relationships of European lineages. Phylogenetic relationships were reconstructed using maximum likelihood and Bayesian inference. As part of our Bayesian analyses, we also dated the phylogeny using two approaches, one based only on fossil calibrations and one that included an additional biogeographical constraint. Our results show paraphyly of the European nesticids with respect to the Asian and North American taxa. We recover four main lineages within Europe. These four European lineages and all European genera have 100% bootstrap support and high posterior probability support in the BEAST2 analysis. The Typhlonesticus lineage is the earliest branching clade present in Europe and includes seven species, the five currently accepted species plus T. parvus from Bosnia and Herzegovina and T. silvestrii from western North America. The Eastern lineage includes the genus Aituaria and is the sister group of the Asian genera Nesticella and Wraios. The Domitius lineage is likely the sister group of the Central European lineage and spreads over the Iberian and Italian peninsulas. Finally, the Central European lineage includes three genera: Kryptonesticus, distributed from the karstic massifs of the Balkan Peninsula to Turkey, Nesticus with a single synanthropic species N. cellulanus and Carpathonesticus, exclusive to the Carpathian Mountains. With the exception of the genus Typhlonesticus, all European genera show an allopatric distribution (except for the two European synanthropic species). The results obtained in this study together with the revision of the original descriptions, redescriptions, and illustrations, lead us to propose 11 nomenclatural changes (new combinations) concerning the genera Typhlonesticus, Nesticus and Carpathonesticus.

Assessing the relative performance of fast molecular dating methods for phylogenomic data

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

The evolutionary history of cribellate orb-weaver capture thread spidroins

Background

Spiders have evolved two types of sticky capture threads: one with wet adhesive spun by ecribellate orb-weavers and another with dry adhesive spun by cribellate spiders. The evolutionary history of cribellate capture threads is especially poorly understood. Here, we use genomic approaches to catalog the spider-specific silk gene family (spidroins) for the cribellate orb-weaver Uloborus diversus.

Results

We show that the cribellar spidroin, which forms the puffy fibrils of cribellate threads, has three distinct repeat units, one of which is conserved across cribellate taxa separated by ~ 250 Mya. We also propose candidates for a new silk type, paracribellar spidroins, which connect the puffy fibrils to pseudoflagelliform support lines. Moreover, we describe the complete repeat architecture for the pseudoflagelliform spidroin (Pflag), which contributes to extensibility of pseudoflagelliform axial fibers.

Conclusions

Our finding that Pflag is closely related to Flag, supports homology of the support lines of cribellate and ecribellate capture threads. It further suggests an evolutionary phase following gene duplication, in which both Flag and Pflag were incorporated into the axial lines, with subsequent loss of Flag in uloborids, and increase in expression of Flag in ecribellate orb-weavers, explaining the distinct mechanical properties of the axial lines of these two groups.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02042-5.

Classification and functional characterization of spidroin genes in a wandering spider, Pardosa pseudoannulata

Spiders impress us with their sophisticated use of silk and the stunningly distinct silk proteins (spidroins) in each spider species. Understanding how silks and spidroins function and evolve within the spider world is one profound interest to expand our knowledge on spider evolution. Spidroins are characterized with the divergent repeat core region flanked with the relatively conserved N- and C-terminus. The structure and number of the repeats contribute to the unique mechanical properties of the spidroin and the silk. Spidroins have been intensively studied in web-weaver spiders, but information regarding their diversity in wandering spiders remains scarce. Here, twenty spidroin genes were identified in the pond wolf spider, Pardosa pseudoannulata, belonging to the retrolateral tibial apophysis (RTA) clade. These spidroins were categorized into four classes, including twelve ampullate spidroin (AmpSp), four aciniform spidroin (AcSp), one tubuliform spidroin (TuSp), one pyriform spidroin (PiSp), and two spidroin-like proteins. Multiple copies of the AmpSp and AcSp genes were tandemly arranged in a cluster within the genome, and the N-terminal domains and repetitive sequences of the proximately located spidroins were highly similar, suggesting that the spidroin genes diversified via tandem duplication. Only four types of morphologically distinct silk glands were found in P. pseudoannulata, namely Ma, Mi, Ac, and Pi glands, consistent with the glandular affiliation hypothesis that spidroins co-evolved with glandular specialization to fit species-specific needs. Expression profiling revealed that the single tubuliform spidroin (TuSp) gene was highly expressed in gravid females and two AcSp genes displayed synchronous expression. Knock-down of the TuSp gene via RNAi resulted in fragile and cracked eggsacs and prolonged the female pre-oviposition period, validating its importance in spider reproduction. The genome-scale characterization and functional study of spidroin genes allows associating the presence of specific spidroins with silk utility in P. pseudoannulata and will expand our knowledge of spider evolution.

1000 spider silkomes: Linking sequences to silk physical properties

Spider silks are among the toughest known materials and thus provide models for renewable, biodegradable, and sustainable biopolymers. However, the entirety of their diversity still remains elusive, and silks that exceed the performance limits of industrial fibers are constantly being found. We obtained transcriptome assemblies from 1098 species of spiders to comprehensively catalog silk gene sequences and measured the mechanical, thermal, structural, and hydration properties of the dragline silks of 446 species. The combination of these silk protein genotype-phenotype data revealed essential contributions of multicomponent structures with major ampullate spidroin 1 to 3 paralogs in high-performance dragline silks and numerous amino acid motifs contributing to each of the measured properties. We hope that our global sampling, comprehensive testing, integrated analysis, and open data will provide a solid starting point for future biomaterial designs.

Phylogenomics and classification of Notropis and related shiners (Cypriniformes: Leuciscidae) and the utility of exon capture on lower taxonomic groups

North American minnows of the Shiner Clade, within the family Leuciscidae, represent one of the most taxonomically complex clades of the order Cypriniformes due to the large number of taxa coupled with conserved morphologies. Species within this clade were moved between genera and subgenera until the community decided to lump many of the unclassified taxa with similar morphologies into one genus, Notropis, which has held up to 325 species. Despite phylogentic studies that began to re-elevate some genera merged into Notropis, such as Cyprinella, Luxilus, Lythrurus, and Pteronotropis, the large genus Notropis remained as a taxonomic repository for many shiners of uncertain placement. Recent molecular advances in sequencing technologies have provided the opportunity to re-examine the Shiner Clade using phylogenomic markers. Using a fish probe kit, we sequenced 90 specimens in 87 species representing 16 genera included in the Shiner Clade, with a resulting dataset of 1,004 loci and 286,455 base pairs. Despite the large dataset, only 32,349 bp (11.29%) were phylogenetically informative. In our maximum likelihood tree, 78% of nodes are 100% bootstrap supported demonstrating the utility of the phylogenomic markers at lower taxonomic levels. Unsurprisingly, species within Notropis as well as Hudsonius, Luxilus, and Alburnops are not resolved as monophyletic groups. Cyprinella is monophyletic if Cyprinella callistia is excluded, and Pteronotropis is monophyletic if it includes Hudsonius cummingsae. Taxonomic changes we propose are: restriction of species included in Alburnops and Notropis, elevation of the subgenus Hydrophlox, expansion of species included in Miniellus, movement of Hudsonius cummingsae to Pteronotropis, and resurrection of the genera Coccotis and Paranotropis. We additionally had two specimens of three species, Notropis atherinoides, Ericymba amplamala, and Pimephales vigilax and found signficant differences between the localities (1,086, 1,424, and 845 nucleotides respectively).

Biological challenges to conclusions from molecular phylogenies: behaviour strongly favours orb web monophyly, contradicting molecular analyses

This first-ever extensive review of the construction behaviour of orb webs, of webs secondarily derived from orbs, and of non-orbs shows that the evidence favouring monophyly over convergent evolution of orbs is stronger than previously appreciated. The two major orb-weaving groups, Uloboridae and Araneoidea, share 31 construction behaviour traits, 20 of which are likely to be both derived and to have feasible alternatives, making convergence an unlikely explanation. Convergence in two lineages seems unlikely, and convergence in five different lineages, as proposed in some recent molecular studies of phylogeny, is even less credible. A further set of seven shared responses in orb design to experimentally constrained spaces also supports orb monophyly. Finally, a ‘control’ case of confirmed convergence on similar ‘pseudo-orbs’ in a taxonomically distant group also supports this argument, as it shows a low frequency of behavioural similarities. I argue that the omission of behavioural data from recent molecular studies of orb web evolution represents a failure of the analytic techniques, not the data, and increases the risk of making mistakes. In general, phylogenetic studies that aim to understand the evolution of particular phenotypes can benefit from including careful study of the phenotypes themselves.

Expression of spidroin proteins in the silk glands of golden orb-weaver spiders

The expression of spidroins in the major ampullate, minor ampullate, flagelliform, and tubuliform silk glands of Trichonephila clavipes spiders was analyzed using proteomics analysis techniques. Spidroin peptides were identified and assigned to different gene products based on sequence concurrence when compared with the whole genome of the spider. It was found that only a relatively low proportion of the spidroin genes are expressed as proteins in any of the studied glands. In addition, the expression of spidroin genes in different glands presents a wide range of patterns, with some spidroins being found in a single gland exclusively, while others appear in the content of several glands. The combination of precise genomics, proteomics, microstructural, and mechanical data provides new insights both on the design principles of these materials and how these principles might be translated for the production of high-performance bioinspired artificial fibers.

Embryonic expression patterns of Wnt genes in the RTA-clade spider Cupiennius salei

Spiders represent widely used model organisms for chelicerate and even arthropod development and evolution. Wnt genes are important and evolutionary conserved factors that control and regulate numerous developmental processes. Recent studies comprehensively investigated the complement and expression of spider Wnt genes revealing conserved as well as diverged aspects of their expression and thus (likely) function among different groups of spiders representing Mygalomorphae (tarantulas), and both main groups of Araneae (true spiders) (Haplogynae/Synspermiata and Entelegynae). The allegedly most modern/derived group of entelegyne spiders is represented by the RTA-clade of which no comprehensive data on Wnt expression were available prior to this study. Here, we investigated the embryonic expression of all Wnt genes of the RTA-clade spider Cupiennius salei. We found that most of the Wnt expression patterns are conserved between Cupiennius and other spiders, especially more basally branching species. Surprisingly, most differences in Wnt gene expression are seen in the common model spider Parasteatoda tepidariorum (a non-RTA clade entelegyne species). These results show that data and conclusions drawn from research on one member of a group of animals (or any other organism) cannot necessarily be extrapolated to the group as a whole, and instead highlight the need for comprehensive taxon sampling.

Improving Taxonomic Practices and Enhancing Its Extensibility—An Example from Araneology

Planetary extinction of biodiversity underscores the need for taxonomy. Here, we scrutinize spider taxonomy over the last decade (2008–2018), compiling 2083 published accounts of newly described species. We evaluated what type of data were used to delineate species, whether data were made freely available, whether an explicit species hypothesis was stated, what types of media were used, the sample sizes, and the degree to which species constructs were integrative. The findings we report reveal that taxonomy remains largely descriptive, not integrative, and provides no explicit conceptual framework. Less than 4% of accounts explicitly stated a species concept and over one-third of all new species described were based on 1–2 specimens or only one sex. Only ~5% of studies made data freely available, and only ~14% of all newly described species employed more than one line of evidence, with molecular data used in ~6% of the studies. These same trends have been discovered in other animal groups, and therefore we find it logical that taxonomists face an uphill challenge when justifying the scientific rigor of their field and securing the needed resources. To move taxonomy forward, we make recommendations that, if implemented, will enhance its rigor, repeatability, and scientific standards.

Biogeography of Long-Jawed Spiders Reveals Multiple Colonization of the Caribbean

Dispersal ability can affect levels of gene flow thereby shaping species distributions and richness patterns. The intermediate dispersal model of biogeography (IDM) predicts that in island systems, species diversity of those lineages with an intermediate dispersal potential is the highest. Here, we tested this prediction on long-jawed spiders (Tetragnatha) of the Caribbean archipelago using phylogenies from a total of 318 individuals delineated into 54 putative species. Our results support a Tetragnatha monophyly (within our sampling) but reject the monophyly of the Caribbean lineages, where we found low endemism yet high diversity. The reconstructed biogeographic history detects a potential early overwater colonization of the Caribbean, refuting an ancient vicariant origin of the Caribbean Tetragnatha as well as the GAARlandia land-bridge scenario. Instead, the results imply multiple colonization events to and from the Caribbean from the mid-Eocene to late-Miocene. Among arachnids, Tetragnatha uniquely comprises both excellently and poorly dispersing species. A direct test of the IDM would require consideration of three categories of dispersers; however, long-jawed spiders do not fit one of these three a priori definitions, but rather represent a more complex combination of attributes. A taxon such as Tetragnatha, one that readily undergoes evolutionary changes in dispersal propensity, can be referred to as a ‘dynamic disperser’.

Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genus-level spider phylogenetics

RAD sequencing yields large amounts of genome-wide data at a relatively low cost and without requiring previous taxon-specific information, making it ideal for evolutionary studies of highly diversified and neglected organisms. However, concerns about information decay with phylogenetic distance have discouraged its use for assessing supraspecific relationships. Here, using Double Digest Restriction Associated DNA (ddRAD) data, we perform the first deep-level approach to the phylogeny of Zodarion, a highly diversified spider genus. We explore the impact of loci and taxon filtering across concatenated and multispecies coalescent reconstruction methods and investigate the patterns of information dropout in reference to both the time of divergence and the mitochondrial divergence between taxa. We found that relaxed loci-filtering and nested taxon-filtering strategies maximized the amount of molecular information and improved phylogenetic inference. As expected, there was a clear pattern of allele dropout towards deeper time and mitochondrial divergences, but the phylogenetic signal remained strong throughout the phylogeny. Therefore, we inferred topologies that were almost fully resolved, highly supported, and noticeably congruent between setups and inference methods, which highlights overall inconsistency in the taxonomy of Zodarion. Because Zodarion appears to be among the oldest and most mitochondrially diversified spider genera, our results suggest that ddRAD data show high potential for inferring intra-generic relationships across spiders and probably also in other taxonomic groups.

Drivers of taxonomic, functional and phylogenetic diversities in dominant ground-dwelling arthropods of coastal heathlands

Although functional and phylogenetic diversities are increasingly used in ecology for a variety of purposes, their relationship remains unclear, and this relationship likely differs among taxa, yet most recent studies focused on plants. We hypothesize that communities may be diverse in functional traits due to presence of: many phylogenetic lineages, trait divergence within lineages, many species and random functional variation among species, weak filtering of traits in favorable environments, or strong trait divergence in unfavorable environments. We tested these predictions for taxa showing higher (ants), or lower (spiders, ground beetles) degrees of competition and niche construction, both of which might decouple functional traits from phylogenetic position or from the environment. Studying > 11,000 individuals and 216 species from coastal heathlands, we estimated functional as minimum spanning trees using traits related to the morphology, feeding habits and dispersal, respectively. Relationships between functional and phylogenetic diversities were overall positive and strong. In ants, this relationship disappeared after accounting for taxonomic diversities and environments, whereas in beetles and spiders taxonomic diversity is related to functional diversity only via increasing phylogenetic diversity. Environmental constraints reduced functional diversity in ants, but affected functional diversity only indirectly via phylogenetic diversity (ground beetles) and taxonomic and then phylogenetic diversity (spiders and ground beetles). Results are consistent with phylogenetic conservatism in traits in spiders and ground beetles. In ants, in contrast, traits appear more phylogenetically neutral with any new species potentially representing a new trait state, tentatively suggesting that competition or niche construction might decouple phylogenetics from trait diversity.

Converging on the orb: denser taxon sampling elucidates spider phylogeny and new analytical methods support repeated evolution of the orb web

High throughput sequencing and phylogenomic analyses focusing on relationships among spiders have both reinforced and upturned long-standing hypotheses. Likewise, the evolution of spider webs-perhaps their most emblematic attribute-is being understood in new ways. With a matrix including 272 spider species and close arachnid relatives, we analyze and evaluate the relationships among these lineages using a variety of orthology assessment methods, occupancy thresholds, tree inference methods and support metrics. Our analyses include families not previously sampled in transcriptomic analyses, such as Symphytognathidae, the only araneoid family absent in such prior works. We find support for the major established spider lineages, including Mygalomorphae, Araneomorphae, Synspermiata, Palpimanoidea, Araneoidea and the Retrolateral Tibial Apophysis Clade, as well as the uloborids, deinopids, oecobiids and hersiliids Grade. Resulting trees are evaluated using bootstrapping, Shimodaira-Hasegawa approximate likelihood ratio test, local posterior probabilities and concordance factors. Using structured Markov models to assess the evolution of spider webs while accounting for hierarchically nested traits, we find multiple convergent occurrences of the orb web across the spider tree-of-life. Overall, we provide the most comprehensive spider tree-of-life to date using transcriptomic data and use new methods to explore controversial issues of web evolution, including the origins and multiple losses of the orb web.

The biology and evolution of spider venoms

Spiders are diverse, predatory arthropods that have inhabited Earth for around 400 million years. They are well known for their complex venom systems that are used to overpower their prey. Spider venoms contain many proteins and peptides with highly specific and potent activities suitable for biomedical or agrochemical applications, but the key role of venoms as an evolutionary innovation is often overlooked, even though this has enabled spiders to emerge as one of the most successful animal lineages. In this review, we discuss these neglected biological aspects of spider venoms. We focus on the morphology of spider venom systems, their major components, biochemical and chemical plasticity, as well as ecological and evolutionary trends. We argue that the effectiveness of spider venoms is due to their unprecedented complexity, with diverse components working synergistically to increase the overall potency. The analysis of spider venoms is difficult to standardize because they are dynamic systems, fine-tuned and modified by factors such as sex, life-history stage and biological role. Finally, we summarize the mechanisms that drive spider venom evolution and highlight the need for genome-based studies to reconstruct the evolutionary history and physiological networks of spider venom compounds with more certainty.

New specimens from Mid-Cretaceous Myanmar amber illuminate the phylogenetic placement of Lagonomegopidae (Arachnida: Araneae)

New lagonomegopid spiders are described from Mid-Cretaceous Myanmar (Burmese) amber. Two new genera and species based on single specimens, Scopomegops fax gen. & sp. nov. and Hiatomegops spinalis gen. & sp. nov. are described. Two specimens belonging to Lineaburmops beigeli are further described. Additionally, after re-examining the holotype of Odontomegops titan, a detailed description of its basal ventral abdomen is added here. A phylogenetic analysis was performed to investigate the phylogenetic placement of Lagonomegopidae. A matrix of 79 morphological characters, scored for six lagonomegopid taxa and 26 non-lagonomegopid taxa, was analysed through parsimony and Bayesian phylogenetic inference. Our results recover extant Palpimanoidea as a monophyletic group and partly suggest that Lagonomegopidae is the sister-group to extant Palpimanoidea. The external sexual organs, retrolateral tibial apophysis on the male palp and tracheal spiracle in lagonomegopids are discussed.

Phylogeny of Micronesian emperor fishes and evolution of trophic types

Island communities that rely on reef fish are currently faced with declining marine resources due to unsustainable fishing and climate change. Identification of genetic stocks through phylogenetic analyses has become a growing field of study with conservation implications, but genetic information on reef fish in Micronesia is limited. In this study we focus on Lethrinidae, one of the most commonly fished reef fish families in Micronesia. Our main goal was to establish a phylogeny for Lethrinidae based on Micronesian data with the intent to help future conservation efforts and clarify the evolutionary history of trophic types in this family. Thirty-eight Lethrinidae specimens collected across five Micronesian islands were used to build a phylogeny with three mitochondrial and one nuclear gene. The phylogenetic analyses allowed us to clarify the identity and position of 11 commonly harvested species and provided a novel genetic identification for Monotaxis heterodon in Micronesia. Our improved and dated phylogeny supports a new hypothesis for the ancestral trophic type of emperor fishes: "stalkers" with low-bodies and conical teeth. We correlated the radiation of most Lethrinidae species with the radiation of major scleractinian coral lineages in the middle Miocene, highlighting the tight relationships between declining reefs and the survival of emperor fishes.

Phylogenomic analyses reveal a Gondwanan origin and repeated out of India colonizations into Asia by tarantulas (Araneae: Theraphosidae)

The study of biogeography seeks taxa that share a key set of characteristics, such as timescale of diversification, dispersal ability, and ecological lability. Tarantulas are ideal organisms for studying evolution over continental-scale biogeography given their time period of diversification, their mostly long-lived sedentary lives, low dispersal rate, and their nevertheless wide circumtropical distribution. In tandem with a time-calibrated transcriptome-based phylogeny generated by PhyloBayes, we estimate the ancestral ranges of ancient tarantulas using two methods, DEC+j and BBM, in the context of their evolution. We recover two ecologically distinct tarantula lineages that evolved on the Indian Plate before it collided with Asia, emphasizing the evolutionary significance of the region, and show that both lineages diversified across Asia at different times. The most ancestral tarantulas emerge on the Americas and Africa 120 Ma-105.5 Ma. We provide support for a dual colonization of Asia by two different tarantula lineages that occur at least 20 million years apart, as well as a Gondwanan origin for the group. We determine that their current distributions are attributable to a combination of Gondwanan vicariance, continental rafting, and geographic radiation. We also discuss emergent patterns in tarantula habitat preferences through time.

Relative Efficiencies of Simple and Complex Substitution Models in Estimating Divergence Times in Phylogenomics

The conventional wisdom in molecular evolution is to apply parameter-rich models of nucleotide and amino acid substitutions for estimating divergence times. However, the actual extent of the difference between time estimates produced by highly complex models compared with those from simple models is yet to be quantified for contemporary data sets that frequently contain sequences from many species and genes. In a reanalysis of many large multispecies alignments from diverse groups of taxa, we found that the use of the simplest models can produce divergence time estimates and credibility intervals similar to those obtained from the complex models applied in the original studies. This result is surprising because the use of simple models underestimates sequence divergence for all the data sets analyzed. We found three fundamental reasons for the observed robustness of time estimates to model complexity in many practical data sets. First, the estimates of branch lengths and node-to-tip distances under the simplest model show an approximately linear relationship with those produced by using the most complex models applied on data sets with many sequences. Second, relaxed clock methods automatically adjust rates on branches that experience considerable underestimation of sequence divergences, resulting in time estimates that are similar to those from complex models. And, third, the inclusion of even a few good calibrations in an analysis can reduce the difference in time estimates from simple and complex models. The robustness of time estimates to model complexity in these empirical data analyses is encouraging, because all phylogenomics studies use statistical models that are oversimplified descriptions of actual evolutionary substitution processes.

Interrogating Genomic-Scale Data to Resolve Recalcitrant Nodes in the Spider Tree of Life

Genome-scale data sets are converging on robust, stable phylogenetic hypotheses for many lineages; however, some nodes have shown disagreement across classes of data. We use spiders (Araneae) as a system to identify the causes of incongruence in phylogenetic signal between three classes of data: exons (as in phylotranscriptomics), noncoding regions (included in ultraconserved elements [UCE] analyses), and a combination of both (as in UCE analyses). Gene orthologs, coded as amino acids and nucleotides (with and without third codon positions), were generated by querying published transcriptomes for UCEs, recovering 1,931 UCE loci (codingUCEs). We expected that congeners represented in the codingUCE and UCEs data would form clades in the presence of phylogenetic signal. Noncoding regions derived from UCE sequences were recovered to test the stability of relationships. Phylogenetic relationships resulting from all analyses were largely congruent. All nucleotide data sets from transcriptomes, UCEs, or a combination of both recovered similar topologies in contrast with results from transcriptomes analyzed as amino acids. Most relationships inferred from low-occupancy data sets, containing several hundreds of loci, were congruent across Araneae, as opposed to high occupancy data matrices with fewer loci, which showed more variation. Furthermore, we found that low-occupancy data sets analyzed as nucleotides (as is typical of UCE data sets) can result in more congruent relationships than high occupancy data sets analyzed as amino acids (as in phylotranscriptomics). Thus, omitting data, through amino acid translation or via retention of only high occupancy loci, may have a deleterious effect in phylogenetic reconstruction.

First chromosomal analysis in Deinopidae (Araneae) reveals Sex Chromosome System X1X2X3X4, B chromosomes and polymorphism for centric fusion

Here we present the first cytogentic study concerning Deinopidae and their controversial phylogenetic position. This study karyologically analyzed one population of Deinopis biaculeata Simon, 1906 and five populations of Deinopis plurituberculata Mello-Leitão, 1925. The majority of specimens of D. plurituberculata exhibited 2n♂ = 40 and 2n♀ = 44 telocentric chromosomes (however some of them showed B chromosomes, belongs to Aquidauana and Botucatu population). The Deinopis biaculeata and D. plurituberculata meiosis of males showed 18 autosomal bivalents + X₁X₂X₃X₄, n = 22 and n = 18, a rare sex chromosome system (SCS) in spiders. Some individuals of D. plurituberculata from the Campo Grande population exhibited 2n♂ = 39 and 2n♀ = 43, with a metacentric chromosome (heterozygotes for centric fusion). The D. plurituberculata males with the rearrangement exhibit diplotenes with 16 autosomal bivalents + 1 autosomal trivalent + X₁X₂X₃X₄ and metaphases II with n = 22 (18 telocentric autosomes + X₁X₂X₃X₄), n = 21 (16 telocentric autosomes + a metacentric autosome + X₁X₂X₃X₄), n = 18 (18 telocentric autosomes) and n = 17 (16 telocentric autosomes + a metacentric autosome). The Ag-NORs (silver impregnation) are terminally located in a pair, coinciding with secondary constriction, which is the most common configuration for Araneae. The relatively high diploid number in Deinopis corroborates phylogenies that place it in a basal position among Entelegynes, in the UDOH grade (Uloboridae, Deinopidae, Oecobiidae and Hersiliidae). Centric fusion in only one population of D. plurituberculata suggests low dispersion capacity of this species and an absence of homozygotes for fusion suggests their low viability or a need to increase the population sampling of D. plurituberculata exhibiting the rearrangement. B chromosomes were detected in D. plurituberculata, with interpopulacional, intrapopulacional and intraindividual numerical variation, with cells presenting 0 - 3 and 0 - 6 B chromosomes in populations of Aquidauana and Botucatu, respectively.

Unique behavioural modifications in the web structure of the cave orb spider Meta menardi (Araneae, Tetragnathidae)

In the last decade there has been a renewed interest in the study of behavioural adaptations to environmental constraints with a focus on adaptations to challenging habitats due to their reduced ecological complexity. However, behavioural studies on organisms adapted to nutrient poor subterranean habitats are few and far between. Here, we compared both morphological traits, in terms of relative leg lengths, and behavioural traits, captured in the geometry of the spider web, between the cave-dwelling spider, Meta menardi, and two aboveground species from the same family (Tetragnathidae); Metellina mengei and Tetragnatha montana. We found that the webs of the cave spider differed significantly from the two surface-dwelling species. The most dramatic difference was the lack of frame threads with the radii in the webs instead attaching directly to the surrounding rock, but other differences in relative web size, web asymmetry and number of capture spiral threads were also found. We argue that these modifications are likely to be adaptations to allow for a novel foraging behaviour to additionally capture walking prey within the vicinity of the web. We found only limited evidence for morphological adaptations and suggest that the cave orb spider could act as a model organism for studies of behaviour in energy-poor environments.

Spider Diversification Through Space and Time

Spiders (Araneae) make up a remarkably diverse lineage of predators that have successfully colonized most terrestrial ecosystems. All spiders produce silk, and many species use it to build capture webs with an extraordinary diversity of forms. Spider diversity is distributed in a highly uneven fashion across lineages. This strong imbalance in species richness has led to several causal hypotheses, such as codiversification with insects, key innovations in silk structure and web architecture, and loss of foraging webs. Recent advances in spider phylogenetics have allowed testing of some of these hypotheses, but results are often contradictory, highlighting the need to consider additional drivers of spider diversification. The spatial and historical patterns of diversity and diversification remain contentious. Comparative analyses of spider diversification will advance only if we continue to make progress with studies of species diversity, distribution, and phenotypic traits, together with finer-scale phylogenies and genomic data.

Global Diversification of Anelosimus Spiders Driven by Long-Distance Overwater Dispersal and Neogene Climate Oscillations

Vicariance and dispersal events, combined with intricate global climatic history, have left an imprint on the spatiotemporal distribution and diversity of many organisms. Anelosimus cobweb spiders (Theridiidae), are organisms ranging in behavior from solitary to highly social, with a cosmopolitan distribution in temperate to tropical areas. Their evolutionary history and the discontinuous distribution of species richness suggest that 1) long-distance overwater dispersal and 2) climate change during the Neogene (23-2.6 Ma), may be major factors in explaining their distribution and diversification. Here, we test these hypotheses, and explicitly test if global Miocene/Pliocene climatic cooling in the last 8 Ma affected Anelosimus radiation in parallel in South America and Madagascar. To do so, we investigate the phylogeny and spatiotemporal biogeography of Anelosimus through a culmination of a 20-year comprehensive global sampling at the species level (69 species, including 84% of the known 75 species worldwide, represented by 268 individuals) using nucleotide data from seven loci (5.5 kb). Our results strongly support the monophyly of Anelosimus with an Oligocene ($\sim $30 Ma) South American origin. Major clades on other continents originate via multiple, long-distance dispersal events, of solitary or subsocial-but not social-lineages, from the Americas. These intercontinental dispersals were to Africa, Madagascar (twice), and SE Asia/Australasia. The early diversification of Anelosimus spiders coincides with a sudden thermal increase in the late Oligocene ($\sim $27-25 Ma), though no causal connection can be made. Our results, however, strongly support the hypothesis that global Neogene climatic cooling in the last 8 Ma drove Anelosimus radiation in parallel in South America and Madagascar, offering a rare empirical evidence for diversification of a socially diverse group driven by an interplay between long-distance dispersal and global Neogene climatic changes. [Cobweb spiders; diversification; global biogeography; long-distance dispersal; molecular phylogenetics; neogene climate changes; sociality; vicariance.].

Tropical Niche Conservatism Explains the Eocene Migration from India to Southeast Asia in Ochyroceratid Spiders

Biological migrations between India and Southeast (SE) Asia provide an ideal system for exploring the effects of geology and climate on species ranges. Geologists have confirmed that the direct collision between India and Eurasia occurred in the Early Eocene, but most migrations occurred between the Indian subcontinent and SE Asia rather than the former and the southern margin of Eurasia. To explain this seemingly paradoxical disconnect between the routes of plate movement and biological migration, we studied the evolutionary history of the tropical spider family Ochyroceratidae based on 101 globally distributed species. We infer a robust dated phylogeny using both transcriptomic data and a data set of classical markers and relate these to biogeographic and climatic analyses. Our results indicate that the monophyly of Ochyroceratidae is strongly supported, and the divergence times suggest a Cretaceous Gondwanan origin of the family. Reconstructed biogeographic histories support a dispersal event from the Indian subcontinent to islands of SE Asia 55-38 Ma. Climatic analyses and the fossil record reveal that ochyroceratids are characterized by a high degree of tropical niche conservatism, and that the ancestor of the Indian and SE Asian clades originated in very warm, wet environments. Early Eocene tropical, perhumid climates in India, and SE Asia may have facilitated ochyroceratid migration, whereas the dry or seasonal climate extending from the eastern coast of China to Central Asia may have acted as a barrier, preventing dispersal. Our analyses suggest that climate plays a more important role than geology in biological migration from the Indian subcontinent to SE Asia, providing new insights into the Indian-Asian biogeographic link. [Biogeography; ecology; geological connections; macroevolution; paleoclimate.].

Phylogenetic Systematics and Evolution of the Spider Infraorder Mygalomorphae Using Genomic Scale Data

The infraorder Mygalomorphae is one of the three main lineages of spiders comprising over 3000 nominal species. This ancient group has a worldwide distribution that includes among its ranks large and charismatic taxa such as tarantulas, trapdoor spiders, and highly venomous funnel-web spiders. Based on past molecular studies using Sanger-sequencing approaches, numerous mygalomorph families (e.g., Hexathelidae, Ctenizidae, Cyrtaucheniidae, Dipluridae, and Nemesiidae) have been identified as non-monophyletic. However, these data were unable to sufficiently resolve the higher-level (intra- and interfamilial) relationships such that the necessary changes in classification could be made with confidence. Here, we present a comprehensive phylogenomic treatment of the spider infraorder Mygalomorphae. We employ 472 loci obtained through anchored hybrid enrichment to reconstruct relationships among all the mygalomorph spider families and estimate the timeframe of their diversification. We sampled nearly all currently recognized families, which has allowed us to assess their status, and as a result, propose a new classification scheme. Our generic-level sampling has also provided an evolutionary framework for revisiting questions regarding silk use in mygalomorph spiders. The first such analysis for the group within a strict phylogenetic framework shows that a sheet web is likely the plesiomorphic condition for mygalomorphs, as well as providing insights to the ancestral foraging behavior for all spiders. Our divergence time estimates, concomitant with detailed biogeographic analysis, suggest that both ancient continental-level vicariance and more recent dispersal events have played an important role in shaping modern day distributional patterns. Based on our results, we relimit the generic composition of the Ctenizidae, Cyrtaucheniidae, Dipluridae, and Nemesiidae. We also elevate five subfamilies to family rank: Anamidae (NEW RANK), Euagridae (NEW RANK), Ischnothelidae (NEW RANK), Pycnothelidae (NEW RANK), and Bemmeridae (NEW RANK). Three families Entypesidae (NEW FAMILY), Microhexuridae (NEW FAMILY), and Stasimopidae (NEW FAMILY), and one subfamily Australothelinae (NEW SUBFAMILY) are newly proposed. Such a major rearrangement in classification, recognizing nine newly established family-level rank taxa, is the largest the group has seen in over three decades. [Biogeography; molecular clocks; phylogenomics; spider web foraging; taxonomy.].

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

Spatial distribution and sample size to estimate Lycosa erythrognatha (Araneae: Lycosidae) population density overwintering

Lycosa erythrognatha Lucas, 1833 (Araneae: Lycosidae) is a predatory arthropod with potential for conservation biological control. In addition to being considered a bioindicator of environmental quality, this arthropod provides an important service for agriculture by reducing insect-pest populations. In this work we seek to understand how the plants Andropogon bicornis L., Saccharum angustifolium Nees and Eustachys retusa Lag (Poales: Poaceae) and their different clump sizes affect the population density, spatial distribution and determination of the minimum number of samples to estimate its population density during the winter. Among the evaluated host plants, S. angustifolium and A. bicornis presented higher population density than E. retusa, but we observed that the clump diameter significantly influences the population density and the minimum number of samples. We observed a gregarious behavior in plants of A. bicornis and E. retusa. For S. angustifolium, a uniform distribution was observed.

Spidroin profiling of cribellate spiders provides insight into the evolution of spider prey capture strategies

Orb-weaving spiders have two main methods of prey capture: cribellate spiders use dry, sticky capture threads, and ecribellate spiders use viscid glue droplets. Predation behaviour is a major evolutionary driving force, and it is important on spider phylogeny whether the cribellate and ecribellate spiders each evolved the orb architecture independently or both strategies were derived from an ancient orb web. These hypotheses have been discussed based on behavioural and morphological characteristics, with little discussion on this subject from the perspective of molecular materials of orb web, since there is little information about cribellate spider-associated spidroin genes. Here, we present in detail a spidroin catalogue of six uloborid species of cribellate orb-weaving spiders, including cribellate and pseudoflagelliform spidroins, with transcriptome assembly complemented with long read sequencing, where silk composition is confirmed by proteomics. Comparative analysis across families (Araneidae and Uloboridae) shows that the gene architecture, repetitive domains, and amino acid frequencies of the orb web constituting silk proteins are similar among orb-weaving spiders regardless of the prey capture strategy. Notably, the fact that there is a difference only in the prey capture thread proteins strongly supports the monophyletic origin of the orb web.

Molecular dating for phylogenies containing a mix of populations and species by using Bayesian and RelTime approaches

Simultaneous molecular dating of population and species divergences is essential in many biological investigations, including phylogeography, phylodynamics and species delimitation studies. In these investigations, multiple sequence alignments consist of both intra- and interspecies samples (mixed samples). As a result, the phylogenetic trees contain interspecies, interpopulation and within-population divergences. Bayesian relaxed clock methods are often employed in these analyses, but they assume the same tree prior for both inter- and intraspecies branching processes and require specification of a clock model for branch rates (independent vs. autocorrelated rates models). We evaluated the impact of a single tree prior on Bayesian divergence time estimates by analysing computer-simulated data sets. We also examined the effect of the assumption of independence of evolutionary rate variation among branches when the branch rates are autocorrelated. Bayesian approach with coalescent tree priors generally produced excellent molecular dates and highest posterior densities with high coverage probabilities. We also evaluated the performance of a non-Bayesian method, RelTime, which does not require the specification of a tree prior or a clock model. RelTime's performance was similar to that of the Bayesian approach, suggesting that it is also suitable to analyse data sets containing both populations and species variation when its computational efficiency is needed.

Machine learning approaches identify male body size as the most accurate predictor of species richness

BACKGROUND

A major challenge in biodiversity science is to understand the factors contributing to the variability of species richness -the number of different species in a community or region - among comparable taxonomic lineages. Multiple biotic and abiotic factors have been hypothesized to have an effect on species richness and have been used as its predictors, but identifying accurate predictors is not straightforward. Spiders are a highly diverse group, with some 48,000 species in 120 families; yet nearly 75% of all species are found within just the ten most speciose families. Here we use a Random Forest machine learning algorithm to test the predictive power of different variables hypothesized to affect species richness of spider genera.

RESULTS

We test the predictive power of 22 variables from spiders' morphological, genetic, geographic, ecological and behavioral landscapes on species richness of 45 genera selected to represent the phylogenetic and biological breath of Araneae. Among the variables, Random Forest analyses find body size (specifically, minimum male body size) to best predict species richness. Multiple Correspondence analysis confirms this outcome through a negative relationship between male body size and species richness. Multiple Correspondence analyses furthermore establish that geographic distribution of congeneric species is positively associated with genus diversity, and that genera from phylogenetically older lineages are species poorer. Of the spider-specific traits, neither the presence of ballooning behavior, nor sexual size dimorphism, can predict species richness.

CONCLUSIONS

We show that machine learning analyses can be used in deciphering the factors associated with diversity patterns. Since no spider-specific biology could predict species richness, but the biologically universal body size did, we believe these conclusions are worthy of broader biological testing. Future work on other groups of organisms will establish whether the detected associations of species richness with small body size and wide geographic ranges hold more broadly.

An Economic Dilemma Between Molecular Weapon Systems May Explain an Arachno-atypical Venom in Wasp Spiders (Argiope bruennichi)

Spiders use venom to subdue their prey, but little is known about the diversity of venoms in different spider families. Given the limited data available for orb-weaver spiders (Araneidae), we selected the wasp spider Argiope bruennichi for detailed analysis. Our strategy combined a transcriptomics pipeline based on multiple assemblies with a dual proteomics workflow involving parallel mass spectrometry techniques and electrophoretic profiling. We found that the remarkably simple venom of A. bruennichi has an atypical composition compared to other spider venoms, prominently featuring members of the cysteine-rich secretory protein, antigen 5 and pathogenesis-related protein 1 (CAP) superfamily and other, mostly high-molecular-weight proteins. We also detected a subset of potentially novel toxins similar to neuropeptides. We discuss the potential function of these proteins in the context of the unique hunting behavior of wasp spiders, which rely mostly on silk to trap their prey. We propose that the simplicity of the venom evolved to solve an economic dilemma between two competing yet metabolically expensive weapon systems. This study emphasizes the importance of cutting-edge methods to encompass the lineages of smaller venomous species that have yet to be characterized in detail, allowing us to understand the biology of their venom systems and to mine this prolific resource for translational research.

Phylogeography of the widespread Caribbean spiny orb weaver Gasteracantha cancriformis

BACKGROUND

Modern molecular analyses are often inconsistent with pre-cladistic taxonomic hypotheses, frequently indicating higher richness than morphological taxonomy estimates. Among Caribbean spiders, widespread species are relatively few compared to the prevalence of single island endemics. The taxonomic hypothesis Gasteracantha cancriformis circumscribes a species with profuse variation in size, color and body form. Distributed throughout the Neotropics, G. cancriformis is the only morphological species of Gasteracantha in the New World in this globally distributed genus.

METHODS

We inferred phylogenetic relationships across Neotropical populations of Gasteracantha using three target genes. Within the Caribbean, we estimated genetic diversity, population structure, and gene flow among island populations.

RESULTS

Our findings revealed a single widespread species of Gasteracantha throughout the Caribbean, G. cancriformis, while suggesting two recently divergent mainland populations that may represent separate species, diverging linages, or geographically isolated demes. The concatenated and COI (Cytochrome c oxidase subunit 1) phylogeny supported a Caribbean clade nested within the New World. Genetic variability was high between island populations for our COI dataset; however, gene flow was also high, especially between large, adjacent islands. We found structured genetic and morphological variation within G. cancriformis island populations; however, this variation does not reflect genealogical relationships. Rather, isolation by distance and local morphological adaptation may explain the observed variation.

Cytogenetics of entelegyne spiders (Arachnida, Araneae) from southern Africa

Spiders represent one of the most studied arachnid orders. They are particularly intriguing from a cytogenetic point of view, due to their complex and dynamic sex chromosome determination systems. Despite intensive research on this group, cytogenetic data from African spiders are still mostly lacking. In this study, we describe the karyotypes of 38 species of spiders belonging to 16 entelegyne families from South Africa and Namibia. In the majority of analysed families, the observed chromosome numbers and morphology (mainly acrocentric) did not deviate from the family-level cytogenetic characteristics based on material from other continents: Tetragnathidae (2n♂ = 24), Ctenidae and Oxyopidae (2n♂ = 28), Sparassidae (2n♂ = 42), Gnaphosidae, Trachelidae and Trochanteriidae (2n♂ = 22), and Salticidae (2n♂ = 28). On the other hand, we identified interspecific variability within Hersiliidae (2n♂ = 33 and 35), Oecobiidae (2n♂ = 19 and 25), Selenopidae (2n♂ = 26 and 29) and Theridiidae (2n♂ = 21 and 22). We examined the karyotypes of Ammoxenidae and Gallieniellidae for the first time. Their diploid counts (2n♂ = 22) correspond to the superfamily Gnaphosoidea and support their placement in this lineage. On the other hand, the karyotypes of Prodidominae (2n♂ = 28 and 29) contrast with all other Gnaphosoidea. Similarly, the unusually high diploid number in Borboropactus sp. (2n♂ = 28) within the otherwise cytogenetically uniform family Thomisidae (mainly 2n♂ = 21-24) supports molecular data suggesting a basal position of the genus in the family. The implementation of FISH methods for visualisation of rDNA clusters facilitated the detection of complex dynamics of numbers of these loci. We identified up to five loci of the 18S rDNA clusters in our samples. Three different sex chromosome systems (X0, X1X20 and X1X2X30) were also detected among the studied taxa.

Reliable Confidence Intervals for RelTime Estimates of Evolutionary Divergence Times

Confidence intervals (CIs) depict the statistical uncertainty surrounding evolutionary divergence time estimates. They capture variance contributed by the finite number of sequences and sites used in the alignment, deviations of evolutionary rates from a strict molecular clock in a phylogeny, and uncertainty associated with clock calibrations. Reliable tests of biological hypotheses demand reliable CIs. However, current non-Bayesian methods may produce unreliable CIs because they do not incorporate rate variation among lineages and interactions among clock calibrations properly. Here, we present a new analytical method to calculate CIs of divergence times estimated using the RelTime method, along with an approach to utilize multiple calibration uncertainty densities in dating analyses. Empirical data analyses showed that the new methods produce CIs that overlap with Bayesian highest posterior density intervals. In the analysis of computer-simulated data, we found that RelTime CIs show excellent average coverage probabilities, that is, the actual time is contained within the CIs with a 94% probability. These developments will encourage broader use of computationally efficient RelTime approaches in molecular dating analyses and biological hypothesis testing.

The fossil record of spiders revisited: implications for calibrating trees and evidence for a major faunal turnover since the Mesozoic

Studies in evolutionary biology and biogeography increasingly rely on the estimation of dated phylogenetic trees using molecular clocks. In turn, the calibration of such clocks is critically dependent on external evidence (i.e. fossils) anchoring the ages of particular nodes to known absolute ages. In recent years, a plethora of new fossil spiders, especially from the Mesozoic, have been described, while the number of studies presenting dated spider phylogenies based on fossil calibrations increased sharply. We critically evaluate 44 of these studies, which collectively employed 67 unique fossils in 180 calibrations. Approximately 54% of these calibrations are problematic, particularly regarding unsupported assignment of fossils to extant clades (44%) and crown (rather than stem) dating (9%). Most of these cases result from an assumed equivalence between taxonomic placement of fossils and their phylogenetic position. To overcome this limitation, we extensively review the literature on fossil spiders, with a special focus on putative synapomorphies and the phylogenetic placement of fossil species with regard to their importance for calibrating higher taxa (families and above) in the spider tree of life. We provide a curated list including 41 key fossils intended to be a basis for future estimations of dated spider phylogenies. In a second step, we use a revised set of 23 calibrations to estimate a new dated spider tree of life based on transcriptomic data. The revised placement of key fossils and the new calibrated tree are used to resolve a long-standing debate in spider evolution - we tested whether there has been a major turnover in the spider fauna between the Mesozoic and Cenozoic. At least 17 (out of 117) extant families have been recorded from the Cretaceous, implying that at least 41 spider lineages in the family level or above crossed the Cretaeous-Paleogene (K-Pg) boundary. The putative phylogenetic affinities of families known only from the Mesozoic suggest that at least seven Cretaceous families appear to have no close living relatives and might represent extinct lineages. There is no unambiguous fossil evidence of the retrolateral tibial apophysis clade (RTA-clade) in the Mesozoic, although molecular clock analyses estimated the major lineages within this clade to be at least ∼100 million years old. Our review of the fossil record supports a major turnover showing that the spider faunas in the Mesozoic and the Cenozoic are very distinct at high taxonomic levels, with the Mesozoic dominated by Palpimanoidea and Synspermiata, while the Cenozoic is dominated by Araneoidea and RTA-clade spiders.

A morphological and combined phylogenetic analysis of pirate spiders (Araneae, Mimetidae): evolutionary relationships, taxonomy and character evolution

Mimetidae is one of the three families within Araneoidea whose members do not spin foraging webs, but are unique in displaying a complex prey-capture behaviour known as aggressive mimicry. Mimetids are distributed worldwide and are most diverse in the tropics of Central and South America. Here we provide a comprehensive phylogeny of pirate spiders (Mimetidae) based on analyses that combine morphological and multigene nucleotide sequence data. We scored 147 morphological characters for 55 mimetids and 16 outgroup taxa and combined it in a total-evidence approach with the sequence data of Benavides et al. (2017) which included two nuclear ribosomal genes, 18S rRNA and 28S rRNA, two mitochondrial ribosomal genes, 12S rRNA and 16S rRNA, the nuclear protein-encoding gene histone H3 and the mitochondrial protein-encoding gene cytochromec oxidase subunitI. We analysed the combined dataset using parsimony, maximum-likelihood and Bayesian inference methods. Our results support the monophyly of Mimetidae and of the genera Gelanor, Ero, Anansi and Australomimetus. Mimetidae is sister to Arkyidae + Tetragnathidae. Mimetus as currently circumscribed is not monophyletic under any analytical approach used, although several lineages within the genus are consistently found in our analyses. We describe, illustrate and discuss the morphological synapomorphies that support the main clades of Mimetidae. The following nomenclatural changes are proposed: Ermetus koreanus (Paik, 1967), the sole species of the genus, is transferred to Ero C.L. Koch, 1836 and thus Ermetus Ponomarev, 2008 is a junior synonym of Ero C.L. Koch, 1836 (new synonymy) and Ero koreana Paik, 1967 becomes a revalidated combination. Phobetinus sagittifer Simon, 1895, the type species of the genus, is transferred to Mimetus Hentz, 1832 and thus Phobetinus Simon, 1895 is a junior synonym of Mimetus Hentz, 1832 (new synonymy), which results in two changes: Mimetus sagittifer (Simon, 1895), new combination and Mimetus investus (Simon, 1909), new combination. Reo latro Brignoli, 1979, the type species of the genus, is transferred to Mimetus and thus Reo Brignoli, 1979 is a junior synonym of Mimetus (new synonymy), which results in the following two changes: Mimetus latro Brignoli, 1979, new combination and Mimetus eutypus Chamberlin & Ivie, 1935, revalidated combination. Arochoides integrans Mello-Leitão, 1935 is transferred to Tetragnathidae (new family placement). The type specimen of Arochoides integrans, the only species in this genus, is a subadult male of Azilia (Tetragnathidae), most likely Azilia histrio Simon, 1895. Arochoides is a junior synonym of Azilia (new synonymy). http://zoobank.org/urn:lsid:zoobank.org:pub:90F6B3DA-232B-428C-BF38-AEA8953D7685

Spider-specific probe set for ultraconserved elements offers new perspectives on the evolutionary history of spiders (Arachnida, Araneae)

Phylogenomic methods have proven useful for resolving deep nodes and recalcitrant groups in the spider tree of life. Across arachnids, transcriptomic approaches may generate thousands of loci, and target-capture methods, using the previously designed arachnid-specific probe set, can target a maximum of about 1,000 loci. Here, we develop a specialized target-capture probe set for spiders that contains over 2,000 ultraconserved elements (UCEs) and then demonstrate the utility of this probe set through sequencing and phylogenetic analysis. We designed the 'spider-specific' probe set using three spider genomes (Loxosceles, Parasteatoda and Stegodyphus) and ensured that the newly designed probe set includes UCEs from the previously designed Arachnida probe set. The new 'spider-specific' probes were used to sequence UCE loci in 51 specimens. The remaining samples included five spider genomes and taxa that were enriched using Arachnida probe set. The 'spider-specific' probes were also used to gather loci from a total of 84 representative taxa across Araneae. On mapping these 84 taxa to the Arachnida probe set, we captured at most 710 UCE loci, while the spider-specific probe set captured up to 1,547 UCE loci from the same taxon sample. Phylogenetic analyses using maximum likelihood and coalescent methods corroborate most nodes resolved by recent transcriptomic analyses, but not all (e.g. UCE data suggest monophyly of 'symphytognathoids'). Our preferred hypothesis based on topology tests, suggests monophyly of the 'symphytognathoids' (the miniature orb weavers), which in previous studies has only been supported by a combination of morphological and behavioural characters.

Phylogeny-Guided Selection of Priority Groups for Venom Bioprospecting: Harvesting Toxin Sequences in Tarantulas as a Case Study

Animal venoms are promising sources of novel drug leads, but their translational potential is hampered by the low success rate of earlier biodiscovery programs, in part reflecting the narrow selection of targets for investigation. To increase the number of lead candidates, here we discuss a phylogeny-guided approach for the rational selection of venomous taxa, using tarantulas (family Theraphosidae) as a case study. We found that previous biodiscovery programs have prioritized the three subfamilies Ornithoctoninae, Selenocosmiinae, and Theraphosinae, which provide almost all of the toxin sequences currently available in public databases. The remaining subfamilies are poorly represented, if at all. These overlooked subfamilies include several that form entire clades of the theraphosid life tree, such as the subfamilies Eumenophorinae, Harpactirinae, and Stromatopelminae, indicating that biodiversity space has not been covered effectively for venom biodiscovery in Theraphosidae. Focusing on these underrepresented taxa will increase the likelihood that promising candidates with novel structures and mechanisms of action can be identified in future bioprospecting programs.