Search

Studying rapid biological changes accompanying the introduction of alien organisms into native ecosystems can provide insights into fundamental ecological and evolutionary theory. While powerful, this quasi-experimental approach is difficult to implement because the timing of invasions and their consequences are hard to predict, meaning that baseline pre-invasion data are often missing. Exceptionally, the eventual arrival of Varroa destructor (hereafter Varroa) in Australia has been predicted for decades. Varroa is a major driver of honeybee declines worldwide, particularly as vectors of diverse RNA viruses. The detection of Varroa in 2022 at over a hundred sites poses a risk of further spread across the continent. At the same time, careful study of Varroa's spread, if it does become established, can provide a wealth of information that can fill knowledge gaps about its effects worldwide. This includes how Varroa affects honeybee populations and pollination. Even more generally, Varroa invasion can serve as a model for evolution, virology and ecological interactions between the parasite, the host and other organisms.

Novel transmission routes can allow infectious diseases to spread, often with devastating consequences. Ectoparasitic varroa mites vector a diversity of RNA viruses and, having switched hosts from the eastern to western honey bees (Apis ceranatoApis mellifera). They provide an opportunity to explore how novel transmission routes shape disease epidemiology. As the principal driver of the spread of Deformed Wing Viruses (mainly DWV-A and DWV-B), varroa infestation has also driven global honey bee health declines. The more virulent DWV-B strain has been replacing the original DWV-A strain in many regions over the past two decades. Yet, how these viruses originated and spread remains poorly understood. Here we use a phylogeographic analysis based on whole genome data to reconstruct the origins and demography of DWV spread. We found that, rather than reemerging in western honey bees after varroa switched hosts, as suggested by previous work, DWV-A most likely originated in Asia and spread in the mid-20thcentury. It also showed a massive population size expansion following the varroa host switch. By contrast, DWV-B was most likely acquired more recently from a source outside Asia, and appears absent from eastern honey bees, the original varroa host. These results highlight the dynamic nature of viral adaptation, whereby a vector’s host switch can give rise to competing and increasingly virulent disease pandemics. The evolutionary novelty and rapid global spread of these host-virus interactions, together with observed spillover into other species, illustrate how increasing globalisation poses urgent threats to biodiversity and food security.

Host switching allows parasites to expand their niches. However, successful switching may require suites of adaptations and may decrease performance on the old host. As a result, reductions in gene flow accompany many host switches, driving speciation. Because host switches tend to be rapid, it is difficult to study them in real time and their demographic parameters remain poorly understood. As a result, fundamental factors that control subsequent parasite evolution, such as the size of the switching population or the extent of immigration from the original host, remain largely unknown. To shed light on the host switching process, we explored how host switches occur in independent host shifts by two ectoparasitic honey bee mites (Varroa destructor and V. jacobsoni). Both switched to the western honey bee (Apis mellifera) after it was brought into contact with their ancestral host (Apis cerana), ~70 and ~12 years ago, respectively. Varroa destructor subsequently caused worldwide collapses of honey bee populations. Using whole-genome sequencing on 63 mites collected in their native ranges from both the ancestral and novel hosts, we were able to reconstruct the known temporal dynamics of the switch. We further found multiple previously undiscovered mitochondrial lineages on the novel host, along with genetic equivalent of tens of individuals that were involved in the initial host switch. Despite being greatly reduced, some gene flow remains between mites adapted to different hosts. Our findings suggest that while reproductive isolation may facilitate fixation of traits beneficial for exploitation of the new host, ongoing genetic exchange may allow genetic amelioration of inbreeding effects.

Insects have been key players in the assessments of biodiversity impacts of anthropogenically driven environmental change, including the evolutionary and ecological impacts of climate change. Populations of Edith’s Checkerspot Butterfly (Euphydryas editha) adapt rapidly to diverse environmental conditions, with numerous high-impact studies documenting these dynamics over several decades. However, studies of the underlying genetic bases of these responses have been hampered by missing genomic resources, limiting the ability to connect genomic responses to environmental change. Using a combination of Oxford Nanopore long reads, haplotype merging, HiC scaffolding followed by Illumina polishing, we generated a highly contiguous and complete assembly (contigs n = 142, N50 = 21.2 Mb, total length = 607.8 Mb; BUSCOs n = 5,286, single copy complete = 97.8%, duplicated = 0.9%, fragmented = 0.3%, missing = 1.0%). A total of 98% of the assembled genome was placed into 31 chromosomes, which displayed large-scale synteny with other well-characterized lepidopteran genomes. The E. editha genome, annotation, and functional descriptions now fill a missing gap for one of the leading field-based ecological model systems in North America.

Background The honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity. Results For both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4 °C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms. Conclusions Ethanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

The study of island community assembly has been fertile ground for developing and testing theoretical ideas in ecology and evolution. The ecoevolutionary trajectory of lineages after colonization has been a particular interest, as this is a key component of understanding community assembly. In this system, existing ideas, such as the taxon cycle, posit that lineages pass through a regular sequence of ecoevolutionary changes after colonization, with lineages shifting toward reduced dispersal ability, increased ecological specialization, and declines in abundance. However, these predictions have historically been difficult to test. Here, we integrate phylogenomics, population genomics, and X‐ray microtomography/3D morphometrics, to test hypotheses for whether the ecomorphological diversity of trap‐jaw ants (Strumigenys) in the Fijian archipelago is assembled primarily through colonization or postcolonization radiation, and whether species show ecological shifts toward niche specialization, toward upland habitats, and decline in abundance after colonization. We infer that most Fijian endemic Strumigenys evolved in situ from a single colonization and have diversified to fill a large fraction of global morphospace occupied by the genus. Within this adaptive radiation, lineages trend to different degrees toward high elevation, reduced dispersal ability, and demographic decline, and we find no evidence of repeated colonization that displaces the initial radiation. Overall these results are only partially consistent with taxon cycle and associated ideas, while highlighting the potential role of priority effects in assembling island communities.

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.

Herbivores use a wide range of factors to choose their host, including their own physiological states, physical characteristics of plants and the degree of competition. Field observations of herbivores in their native habitats provide a means for simultaneously estimating the relative importance of these factors, but statistical analysis of all these factors may be challenging. Here we used a 7-week dataset of leaf-cutting ant (Atta cephalotes) foraging in a diverse Neotropical arboretum containing 193 tree species (822 trees) to examine the relative role of tree phylogeny, territoriality and tree functional characteristics using a phylogenetic generalized least squares (PGLS) model. We observed that 54 tree species (117 trees) were foraged by the ants. This pattern was not random, but reflected known features of leaf-cutting ant foraging behaviour, such as a preference for larger trees and the decreased likelihood of foraging at the periphery of a colony’s territory. However, random effects such as tree phylogeny, the identity of individual trees and colony-specific effects explained most of the variation in foraging data. A significant phylogenetic effect on foraging likelihood (λ = 0.28), together with repeated measures of foraging on the same tree species, allowed estimation of relative palatability for each plant species. PGLS models can be flexibly scaled to include other covariates for even more complex investigation of foraging behaviour, and the link function can be modified to include the amount of plants foraged. As a result, PGLS can be used as a flexible framework for the study of LCA foraging.

While ectoparasitic Varroa mites cause minimal damage to their co-evolved ancestral host, the eastern honey bee (Apis cerana), they devastate their novel host, the western honey bee (Apis mellifera). Over several decades, the host switch caused worldwide population collapses, threatening global food security. Varroa management strategies have focused on breeding bees for tolerance. But, can Varroa overcome these counter-adaptations in a classic coevolutionary arms race? Despite increasing evidence for Varroa genetic diversity and evolvability, this eventuality has largely been neglected. We therefore suggest a more holistic paradigm for studying this host-parasite interaction, one in which ‘Varroa-tolerant’ bee traits should be viewed as a shared phenotype resulting from Varroa and honey bee interaction.

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities., Maeva Techer et al. report genome assemblies of two honeybee parasitic mites, Varroa destructor and V. jacobsoni. They find that 4.8% of orthologous genes show evidence of positive selection in at least one species, though the genes under selection are distinct between species indicating divergent evolution.

Background Vector-borne viral diseases threaten human and wildlife worldwide. Vectors are often viewed as a passive syringe injecting the virus. However, to survive, replicate and spread, viruses must manipulate vector biology. While most vector-borne viral research focuses on vectors transmitting a single virus, in reality, vectors often carry diverse viruses. Yet how viruses affect the vectors remains poorly understood. Here, we focused on the varroa mite (Varroa destructor), an emergent parasite that can carry over 20 honey bee viruses, and has been responsible for colony collapses worldwide, as well as changes in global viral populations. Co-evolution of the varroa and the viral community makes it possible to investigate whether viruses affect vector gene expression and whether these interactions affect viral epidemiology. Results Using a large set of available varroa transcriptomes, we identified how abundances of individual viruses affect the vector’s transcriptional network. We found no evidence of competition between viruses, but rather that some virus abundances are positively correlated. Furthermore, viruses that are found together interact with the vector’s gene co-expression modules in similar ways, suggesting that interactions with the vector affect viral epidemiology. We experimentally validated this observation by silencing candidate genes using RNAi and found that the reduction in varroa gene expression was accompanied by a change in viral load. Conclusions Combined, the meta-transcriptomic analysis and experimental results shed light on the mechanism by which viruses interact with each other and with their vector to shape the disease course.

Vector-borne viral diseases threaten human and wildlife worldwide. Vectors are often viewed as a passive syringe injecting the virus, however to survive, replicate and spread, viruses must manipulate vector biology. While most vector-borne viral research focuses on vectors transmitting a single virus, in reality vectors often carry diverse viruses. Yet how viruses affect the vectors remains poorly understood. Here we focused on the varroa mite, an emergent parasite that vectors over 20 honey bee viruses, and has been responsible for colony collapses worldwide, as well as changes in global viral populations. Co-evolution of the varroa and the viral community makes it possible to investigate whether viruses affect vector gene expression, and whether these interactions affect viral epidemiology. Using a large set of available varroa transcriptomes we identified how abundances of individual viruses affect the vector’s transcriptional network. Perhaps surprisingly, we found no evidence of competition between viruses, but rather that some virus abundances are positively correlated. Furthermore, we found a strong correlation between the extent to which a virus interacts with the vector’s gene expression and co-occurrence with other viruses, suggesting that interactions with the vector affect epidemiology. We experimentally validated this observation by silencing candidate genes using RNAi and found that reduction in varroa gene expression was accompanied by a change in viral load. Combined, the meta-transcriptomic analysis and experimental results shed light on the mechanism by which viruses interact with each other and with their vector to shape the disease course.

Evolution can occur with surprising predictability when organisms face similar ecological challenges. How and why traits arise repeatedly remains a central question in evolutionary biology, but the complexity of most traits makes it challenging to answer. Reptiles and mammals independently evolved oral venoms consisting of proteinaceous cocktails that allow mapping between genotype and phenotype. Although biochemically similar toxins can occur as major venom components across many taxa, whether these toxins evolved via convergent or parallel processes remains unknown. Most notable among them are kallikrein-like serine proteases that form the core of most vertebrate oral venoms. We used a combination of comparative genomics and phylogenetics to investigate whether serine protease recruitment into the venom occurred independently or in parallel in mammalian and reptilian lineages. Using syntenic relationships between genes flanking known toxins, we traced the origin of kallikreins to a single locus containing one or more nearby paralogous kallikrein-like clusters. Independently, phylogenetic analysis of vertebrate serine proteases revealed that kallikrein-like toxins in mammals and reptiles are homologous. Given the shared regulatory and genetic machinery, these findings suggest a unified model underlying vertebrate venom evolution. Namely, the common tetrapod ancestor possessed salivary genes that were biochemically suitable for envenomation. We term such genes ‘toxipotent’ – in the case of salivary kallikreins they already had potent vasodilatory activity that was weaponized by venomous lineages. Furthermore, the ubiquitous distribution of kallikreins across vertebrates suggests that the evolution of envenomation may be more common than previously recognized, blurring the line between venomous and non-venomous animals.

Primeval populations replicating at high error rates required a mechanism to overcome the accumulation of mutations and information deterioration. Known strategies to overcome mutation pressures include RNA processivity, epistasis, selection, and quasispecies. We investigated the mechanism by which small molecular ribozyme populations can survive under high error rates by propagating several lineages under different mutagen concentrations. We found that every population that evolved without mutagen went extinct, while those subjected to mutagenic evolution survived. To understand how they survived, we characterized the evolved genotypic diversity, the formation of genotype-genotype interaction networks, the fitness of the most common mutants for each enzymatic step, and changes in population size along the course of evolution. We found that the elevated mutation rate was necessary for the populations to survive in the novel environment, in which all the steps of the metabolism worked to promote the survival of even less catalytically efficient ligases. Besides, an increase in population size and the mutational coupling of genotypes in close-knit networks, which helped maintain or recover lost genotypes making their disappearance transient, prevented Muller's ratchet and extinction.

In this study, we report the first complete mitochondrial genome sequence of the Aquatic Coralsnake Micrurus surinamensis. The mitochondrial genome lengthis 17,375 bp, comprising 13 protein-coding genes, 2 rRNA (12S and 16S) and 22 tRNA, as well as two typical control regions. Phylogenetic analysis based upon 13 protein-coding genes showed clusters based on terrestrial and marine species.

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of host-symbiont interactions between the northernmost fungus-farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome-wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant-fungus genome-genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population-genetic structure was concordant between the ants and one cultivar type (M-fungi, concordant clines) but discordant for the other cultivar type (T-fungi). Discordance in population-genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between-nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant-fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.

Traditional static publication formats make visualization, exploration, and sharing of massive phylogenetic trees difficult. A phylogenetic study often involves hundreds of taxa, and the resulting tree has to be split across multiple journal pages, or be shrunk onto one, which jeopardizes legibility. Furthermore, additional data layers, such as species-specific information or time calibrations are often displayed in separate figures, making the entire picture difficult for readers to grasp. Web-based technologies, such as the Data Driven Document (D3) JavaScript library, were created to overcome such challenges by allowing interactive displays of complex data sets. The new phylogeny.IO web server (https://phylogeny.io) overcomes this issue by allowing users to easily import, annotate, and share interactive phylogenetic trees. It allows a range of static (e.g. such as shapes and colors) and dynamic (e.g. pop-up text and images) annotations. Annotated trees can be saved on the server for subsequent modification or they may be shared as IFrame HTML objects, easily embeddable in any web page. The principal goal of phylogeny.IO is not to produce publication-ready figures, but rather to provide a simple and intuitive annotation interface that allows easy and rapid sharing of figures in blogs, lecture notes, press releases, etc.

Gene expression changes contribute to complex trait variations in both individuals and populations. However, the evolution of gene expression underlying complex traits over macroevolutionary timescales remains poorly understood. Snake venoms are proteinaceous cocktails where the expression of each toxin can be quantified and mapped to a distinct genomic locus and traced for millions of years. Using a phylogenetic generalized linear mixed model, we analyzed expression data of toxin genes from 52 snake species spanning the 3 venomous snake families and estimated phylogenetic covariance, which acts as a measure of evolutionary constraint. We find that evolution of toxin combinations is not constrained. However, although all combinations are in principle possible, the actual dimensionality of phylomorphic space is low, with envenomation strategies focused around only four major toxin families: metalloproteases, three-finger toxins, serine proteases, and phospholipases A2. Although most extant snakes prioritize either a single or a combination of major toxin families, they are repeatedly recruited and lost. We find that over macroevolutionary timescales, the venom phenotypes were not shaped by phylogenetic constraints, which include important microevolutionary constraints such as epistasis and pleiotropy, but more likely by ecological filtering that permits a small number of optimal solutions. As a result, phenotypic optima were repeatedly attained by distantly related species. These results indicate that venoms evolve by selection on biochemistry of prey envenomation, which permit diversity through parallelism, and impose strong limits, since only a few of the theoretically possible strategies seem to work well and are observed in extant snakes.

BackgroundThe ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybeeApis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.ResultsEach of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor >98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds ofAvaIandAluIrepeats associated with these features.ConclusionsThe improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Oral venom systems evolved multiple times in numerous vertebrates, enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on toxins. However, using toxins present in modern-day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with toxins and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though they still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Evolutionary innovations underlie the rise of diversity and complexity—the 2 long-term trends in the history of life. How does natural selection redesign multiple interacting parts to achieve a new emergent function? We investigated the evolution of a biomechanical innovation, the latch-spring mechanism of trap-jaw ants, to address 2 outstanding evolutionary problems: how form and function change in a system during the evolution of new complex traits, and whether such innovations and the diversity they beget are repeatable in time and space. Using a new phylogenetic reconstruction of 470 species, and X-ray microtomography and high-speed videography of representative taxa, we found the trap-jaw mechanism evolved independently 7 to 10 times in a single ant genus (Strumigenys), resulting in the repeated evolution of diverse forms on different continents. The trap mechanism facilitates a 6 to 7 order of magnitude greater mandible acceleration relative to simpler ancestors, currently the fastest recorded acceleration of a resettable animal movement. We found that most morphological diversification occurred after evolution of latch-spring mechanisms, which evolved via minor realignments of mouthpart structures. This finding, whereby incremental changes in form lead to a change of function, followed by large morphological reorganization around the new function, provides a model for understanding the evolution of complex biomechanical traits, as well as insights into why such innovations often happen repeatedly., The mousetrap-like mandibles of trap-jaw ants are a functional innovation capable of record-breaking acceleration. This study uses phylogenomics, X-ray microtomography and high-speed videography to reveal that they evolved multiple times around the world through a pathway of incremental steps.

Sociality is a striking phenotypic innovation that independently evolved dozens of times across animals. Sociogenomic approaches have begun to elucidate the molecular underpinnings of social life in social insects and vertebrates, but the degree to which the convergent evolution of sociality involves convergent molecular evolution remains controversial and largely unknown. Spiders are a powerful system for identifying the genomic causes and consequences of social life because sociality is estimated to have independently evolved 15 times, and each origin likely occurred recently, within the past few million years. To determine if there are statistically supported genomic signatures of protein-coding sequence evolution associated with the convergent evolution of sociality in spiders, we compared the genomes or transcriptomes of 24 spider species that vary in social organization and represent at least seven independent origins of sociality. We identified hundreds of genes that experienced shifts in patterns of molecular evolution during the convergent evolution of sociality, and these genes were enriched for several annotated functions, including neural function, neurogenesis, and behavior, as well as immune function, growth, and metabolism. We also found evidence that directional selection for specific substitutions repeatedly occurred in social species for several genes, in particular the calcium channel gene TPC1. Finally, supporting previous results, we found elevated genome-wide rates of molecular evolution in social species, resulting mainly from relaxation of selection. Altogether, we identify genome-wide, genic, and site-specific changes that repeatedly occurred during the evolution of sociality in spiders. Significance StatementThe transition from solitary to social life is a major transition in evolution that repeatedly occurred, but the genetic underpinnings are largely unknown. To identify genomic changes associated with sociality, we compared the genomes of 24 spider species, representing seven recent and independent origins of sociality. We identified hundreds of genes and many functional classes of genes that tended to experience shifts in molecular evolution, including genes that have previously been implicated in animal social behavior and human behavioral disorders. Our study shows that while the precise genetic details vary, the repeated evolution of sociality in spiders predictably leaves genome-wide signatures and involves sets of genes with conserved functions that may play general roles in the evolution of social behavior.