Genomic Signatures of Recent Adaptation in a Wild Bumblebee.

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.

Transcriptomic responses to location learning by honeybee dancers are partly mirrored in the brains of dance-followers.

The waggle dances of honeybees are a strikingly complex form of animal communication that underlie the collective foraging behaviour of colonies. The mechanisms by which bees assess the locations of forage sites that they have visited for representation on the dancefloor are now well-understood, but few studies have considered the remarkable backward translation of such information into flight vectors by dance-followers. Here, we explore whether the gene expression patterns that are induced through individual learning about foraging locations are mirrored when bees learn about those same locations from their nest-mates. We first confirmed that the mushroom bodies of honeybee dancers show a specific transcriptomic response to learning about distance, and then showed that approximately 5% of those genes were also differentially expressed by bees that follow dances for the same foraging sites, but had never visited them. A subset of these genes were also differentially expressed when we manipulated distance perception through an optic flow paradigm, and responses to learning about target direction were also in part mirrored in the brains of dance followers. Our findings show a molecular footprint of the transfer of learnt information from one animal to another through this extraordinary communication system, highlighting the dynamic role of the genome in mediating even very short-term behavioural changes.

Expression of subunits of an insecticide target receptor varies across tissues, life stages, castes, and species of social bees.

Global losses of insects jeopardize ecosystem stability and crop pollination. Robust evidence indicates that insecticides have contributed to these losses. Notably, insecticides targeting nicotinic acetylcholine receptors (nAChRs) have neurotoxic effects on beneficial insects. Because each nAChR consists of five subunits, the alternative arrangements of subunits could create a multitude of receptors differing in structure and function. Therefore, understanding whether the use of subunits varies is essential for evaluating and predicting the effects of insecticides targeting such receptors. To better understand how the use and composition of nAChRs differ within and between insect pollinators, we analysed RNA-seq gene expression data from tissues and castes of Apis mellifera honey bees and life stages and castes of the Bombus terrestris bumble bees. We reveal that all analysed tissues express nAChRs and that relative expression levels of nAChR subunits vary widely across almost all comparisons. Our work thus shows fine-tuned spatial and temporal expression of nAChRs. Given that coexpression of subunits underpins the compositional diversity of functional receptors and that the affinities of insecticides depend on nAChR composition, our findings provide a likely mechanism for the various damaging effects of nAChR-targeting insecticides on insects. Furthermore, our results indicate that the appraisal of insecticide risks should carefully consider variation in molecular targets.

Larger, more connected societies of ants have a higher prevalence of viruses.

The benefits of cooperative living for foraging, nesting, defence and buffering environmental challenges lead animals with the most highly social lifestyles to dominate many ecosystems. However, living in larger, more highly connected groups should also increase the risks of pathogen exposure and transmission. While over long timescales selective responses could buffer the impacts of potential higher pathogen prevalence, similar processes are unlikely over short timescales. The red fire ant Solenopsis invicta is ideal for measuring the effects of group size on pathogen prevalence because two types of society coexist in this species: smaller single-nest single-queen colonies that are highly aggressive to their neighbours and larger multiple-queen colonies that exchange resources with neighbouring nests. We compare the presence of viruses between these two colony types using metagenomic sequence classification of RNA-sequencing reads. We find that queens from multiple-queen colonies have 8.3-times higher viral load and 1.5-times higher viral diversity than queens from single-queen colonies. This finding characterizes a rarely considered cost of transitions to more highly social living. Furthermore, our results show that highly social invertebrates can harbour many viruses.

Degenerative Expansion of a Young Supergene.

Long-term suppression of recombination ultimately leads to gene loss, as demonstrated by the depauperate Y and W chromosomes of long-established pairs of XY and ZW chromosomes. The young social supergene of the Solenopsis invicta red fire ant provides a powerful system to examine the effects of suppressed recombination over a shorter timescale. The two variants of this supergene are carried by a pair of heteromorphic chromosomes, referred to as the social B and social b (SB and Sb) chromosomes. The Sb variant of this supergene changes colony social organization and has an inheritance pattern similar to a Y or W chromosome because it is unable to recombine. We used high-resolution optical mapping, k-mer distribution analysis, and quantification of repetitive elements on haploid ants carrying alternate variants of this young supergene region. We find that instead of shrinking, the Sb variant of the supergene has increased in length by more than 30%. Surprisingly, only a portion of this length increase is due to consistent increases in the frequency of particular classes of repetitive elements. Instead, haplotypes of this supergene variant differ dramatically in the amounts of other repetitive elements, indicating that the accumulation of repetitive elements is a heterogeneous and dynamic process. This is the first comprehensive demonstration of degenerative expansion in an animal and shows that it occurs through nonlinear processes during the early evolution of a region of suppressed recombination.

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases.

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

The first draft genomes of the ant Formica exsecta, and its Wolbachia endosymbiont reveal extensive gene transfer from endosymbiont to host.

BACKGROUND: Adapting to changes in the environment is the foundation of species survival, and is usually thought to be a gradual process. However, transposable elements (TEs), epigenetic modifications, and/or genetic material acquired from other organisms by means of horizontal gene transfer (HGTs), can also lead to novel adaptive traits. Social insects form dense societies, which attract and maintain extra- and intracellular accessory inhabitants, which may facilitate gene transfer between species. The wood ant Formica exsecta (Formicidae; Hymenoptera), is a common ant species throughout the Palearctic region. The species is a well-established model for studies of ecological characteristics and evolutionary conflict. RESULTS: In this study, we sequenced and assembled draft genomes for F. exsecta and its endosymbiont Wolbachia. The F. exsecta draft genome is 277.7 Mb long; we identify 13,767 protein coding genes, for which we provide gene ontology and protein domain annotations. This is also the first report of a Wolbachia genome from ants, and provides insights into the phylogenetic position of this endosymbiont. We also identified multiple horizontal gene transfer events (HGTs) from Wolbachia to F. exsecta. Some of these HGTs have also occurred in parallel in multiple other insect genomes, highlighting the extent of HGTs in eukaryotes. CONCLUSION: We present the first draft genome of ant F. exsecta, and its endosymbiont Wolbachia (wFex), and show considerable rates of gene transfer from the symbiont to the host. We expect that especially the F. exsecta genome will be valuable resource in further exploration of the molecular basis of the evolution of social organization.

Caste- and pesticide-specific effects of neonicotinoid pesticide exposure on gene expression in bumblebees.

Social bees are important insect pollinators of wildflowers and agricultural crops, making their reported declines a global concern. A major factor implicated in these declines is the widespread use of neonicotinoid pesticides. Indeed, recent research has demonstrated that exposure to low doses of these neurotoxic pesticides impairs bee behaviours important for colony function and survival. However, our understanding of the molecular-genetic pathways that lead to such effects is limited, as is our knowledge of how effects may differ between colony members. To understand what genes and pathways are affected by exposure of bumblebee workers and queens to neonicotinoid pesticides, we implemented a transcriptome-wide gene expression study. We chronically exposed Bombus terrestriscolonies to either clothianidin or imidacloprid at field-realistic concentrations while controlling for factors including colony social environment and worker age. We reveal that genes involved in important biological processes including mitochondrial function are differentially expressed in response to neonicotinoid exposure. Additionally, clothianidin exposure had stronger effects on gene expression amplitude and alternative splicing than imidacloprid. Finally, exposure affected workers more strongly than queens. Our work demonstrates how RNA-Seq transcriptome profiling can provide detailed novel insight on the mechanisms mediating pesticide toxicity to a key insect pollinator.

A Y-like social chromosome causes alternative colony organization in fire ants.

Intraspecific variability in social organization is common, yet the underlying causes are rarely known. In the fire ant Solenopsis invicta, the existence of two divergent forms of social organization is under the control of a single Mendelian genomic element marked by two variants of an odorant-binding protein gene. Here we characterize the genomic region responsible for this important social polymorphism, and show that it is part of a pair of heteromorphic chromosomes that have many of the key properties of sex chromosomes. The two variants, hereafter referred to as the social B and social b (SB and Sb) chromosomes, are characterized by a large region of approximately 13 megabases (55% of the chromosome) in which recombination is completely suppressed between SB and Sb. Recombination seems to occur normally between the SB chromosomes but not between Sb chromosomes because Sb/Sb individuals are non-viable. Genomic comparisons revealed limited differentiation between SB and Sb, and the vast majority of the 616 genes identified in the non-recombining region are present in the two variants. The lack of recombination over more than half of the two heteromorphic social chromosomes can be explained by at least one large inversion of around 9 megabases, and this absence of recombination has led to the accumulation of deleterious mutations, including repetitive elements in the non-recombining region of Sb compared with the homologous region of SB. Importantly, most of the genes with demonstrated expression differences between individuals of the two social forms reside in the non-recombining region. These findings highlight how genomic rearrangements can maintain divergent adaptive social phenotypes involving many genes acting together by locally limiting recombination.

Foraging bumblebees acquire a preference for neonicotinoid-treated food with prolonged exposure.

Social bees represent an important group of pollinating insects that can be exposed to potentially harmful pesticides when foraging on treated or contaminated flowering plants. To investigate if such exposure is detrimental to bees, many studies have exclusively fed individuals with pesticide-spiked food, informing us about the hazard but not necessarily the risk of exposure. While such studies are important to establish the physiological and behavioural effects on individuals, they do not consider the possibility that the risk of exposure may change over time. For example, many pesticide assays exclude potential behavioural adaptations to novel toxins, such as rejection of harmful compounds by choosing to feed on an uncontaminated food source, thus behaviourally lowering the risk of exposure. In this paper, we conducted an experiment over 10 days in which bumblebees could forage on an array of sucrose feeders containing 0, 2 and 11 parts per billion of the neonicotinoid pesticide thiamethoxam. This more closely mimics pesticide exposure in the wild by allowing foraging bees to (i) experience a field realistic range of pesticide concentrations across a chronic exposure period, (ii) have repeated interactions with the pesticide in their environment, and (iii) retain the social cues associated with foraging by using whole colonies. We found that the proportion of visits to pesticide-laced feeders increased over time, resulting in greater consumption of pesticide-laced sucrose relative to untreated sucrose. After changing the spatial position of each feeder, foragers continued to preferentially visit the pesticide-laced feeders which indicates that workers can detect thiamethoxam and alter their behaviour to continue feeding on it. The increasing preference for consuming the neonicotinoid-treated food therefore increases the risk of exposure for the colony during prolonged pesticide exposure. Our results highlight the need to incorporate attractiveness of pesticides to foraging bees (and potentially other insect pollinators) in addition to simply considering the proportion of pesticide-contaminated floral resources within the foraging landscape.

GeneValidator: identify problems with protein-coding gene predictions.

UNLABELLED: : Genomes of emerging model organisms are now being sequenced at very low cost. However, obtaining accurate gene predictions remains challenging: even the best gene prediction algorithms make substantial errors and can jeopardize subsequent analyses. Therefore, many predicted genes must be time-consumingly visually inspected and manually curated. We developed GeneValidator (GV) to automatically identify problematic gene predictions and to aid manual curation. For each gene, GV performs multiple analyses based on comparisons to gene sequences from large databases. The resulting report identifies problematic gene predictions and includes extensive statistics and graphs for each prediction to guide manual curation efforts. GV thus accelerates and enhances the work of biocurators and researchers who need accurate gene predictions from newly sequenced genomes. AVAILABILITY AND IMPLEMENTATION: GV can be used through a web interface or in the command-line. GV is open-source (AGPL), available at https://wurmlab.github.io/tools/genevalidator CONTACT: : y.wurm@qmul.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The fire ant social chromosome supergene variant Sb shows low diversity but high divergence from SB.

Variation in social behaviour is common, yet little is known about the genetic architectures underpinning its evolution. A rare exception is in the fire ant Solenopsis invicta: Alternative variants of a supergene region determine whether a colony will have exactly one or up to dozens of queens. The two variants of this region are carried by a pair of 'social chromosomes', SB and Sb, which resemble a pair of sex chromosomes. Recombination is suppressed between the two chromosomes in the supergene region. While the X-like SB can recombine with itself in SB/SB queens, recombination is effectively absent in the Y-like Sb because Sb/Sb queens die before reproducing. Here, we analyse whole-genome sequences of eight haploid SB males and eight haploid Sb males. We find extensive SB-Sb differentiation throughout the >19-Mb-long supergene region. We find no evidence of 'evolutionary strata' with different levels of divergence comparable to those reported in several sex chromosomes. A high proportion of substitutions between the SB and Sb haplotypes are nonsynonymous, suggesting inefficacy of purifying selection in Sb sequences, similar to that for Y-linked sequences in XY systems. Finally, we show that the Sb haplotype of the supergene region has 635-fold less nucleotide diversity than the rest of the genome. We discuss how this reduction could be due to a recent selective sweep affecting Sb specifically or associated with a population bottleneck during the invasion of North America by the sampled population.

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Vitellogenin underwent subfunctionalization to acquire caste and behavioral specific expression in the harvester ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

Sociogenomics of cooperation and conflict during colony founding in the fire ant Solenopsis invicta.

One of the fundamental questions in biology is how cooperative and altruistic behaviors evolved. The majority of studies seeking to identify the genes regulating these behaviors have been performed in systems where behavioral and physiological differences are relatively fixed, such as in the honey bee. During colony founding in the monogyne (one queen per colony) social form of the fire ant Solenopsis invicta, newly-mated queens may start new colonies either individually (haplometrosis) or in groups (pleometrosis). However, only one queen (the "winner") in pleometrotic associations survives and takes the lead of the young colony while the others (the "losers") are executed. Thus, colony founding in fire ants provides an excellent system in which to examine the genes underpinning cooperative behavior and how the social environment shapes the expression of these genes. We developed a new whole genome microarray platform for S. invicta to characterize the gene expression patterns associated with colony founding behavior. First, we compared haplometrotic queens, pleometrotic winners and pleometrotic losers. Second, we manipulated pleometrotic couples in order to switch or maintain the social ranks of the two cofoundresses. Haplometrotic and pleometrotic queens differed in the expression of genes involved in stress response, aging, immunity, reproduction and lipid biosynthesis. Smaller sets of genes were differentially expressed between winners and losers. In the second experiment, switching social rank had a much greater impact on gene expression patterns than the initial/final rank. Expression differences for several candidate genes involved in key biological processes were confirmed using qRT-PCR. Our findings indicate that, in S. invicta, social environment plays a major role in the determination of the patterns of gene expression, while the queen's physiological state is secondary. These results highlight the powerful influence of social environment on regulation of the genomic state, physiology and ultimately, social behavior of animals.

The genomic impact of 100 million years of social evolution in seven ant species.

Ants (Hymenoptera, Formicidae) represent one of the most successful eusocial taxa in terms of both their geographic distribution and species number. The publication of seven ant genomes within the past year was a quantum leap for socio- and ant genomics. The diversity of social organization in ants makes them excellent model organisms to study the evolution of social systems. Comparing the ant genomes with those of the honeybee, a lineage that evolved eusociality independently from ants, and solitary insects suggests that there are significant differences in key aspects of genome organization between social and solitary insects, as well as among ant species. Altogether, these seven ant genomes open exciting new research avenues and opportunities for understanding the genetic basis and regulation of social species, and adaptive complex systems in general.

The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming.

We present a high-quality (>100× depth) Illumina genome sequence of the leaf-cutting ant Acromyrmex echinatior, a model species for symbiosis and reproductive conflict studies. We compare this genome with three previously sequenced genomes of ants from different subfamilies and focus our analyses on aspects of the genome likely to be associated with known evolutionary changes. The first is the specialized fungal diet of A. echinatior, where we find gene loss in the ant's arginine synthesis pathway, loss of detoxification genes, and expansion of a group of peptidase proteins. One of these is a unique ant-derived contribution to the fecal fluid, which otherwise consists of "garden manuring" fungal enzymes that are unaffected by ant digestion. The second is multiple mating of queens and ejaculate competition, which may be associated with a greatly expanded nardilysin-like peptidase gene family. The third is sex determination, where we could identify only a single homolog of the feminizer gene. As other ants and the honeybee have duplications of this gene, we hypothesize that this may partly explain the frequent production of diploid male larvae in A. echinatior. The fourth is the evolution of eusociality, where we find a highly conserved ant-specific profile of neuropeptide genes that may be related to caste determination. These first analyses of the A. echinatior genome indicate that considerable genetic changes are likely to have accompanied the transition from hunter-gathering to agricultural food production 50 million years ago, and the transition from single to multiple queen mating 10 million years ago.

Effects of ploidy and sex-locus genotype on gene expression patterns in the fire ant Solenopsis invicta.

Males in many animal species differ greatly from females in morphology, physiology and behaviour. Ants, bees and wasps have a haplodiploid mechanism of sex determination whereby unfertilized eggs become males while fertilized eggs become females. However, many species also have a low frequency of diploid males, which are thought to develop from diploid eggs when individuals are homozygous at one or more sex determination loci. Diploid males are morphologically similar to haploids, though often larger and typically sterile. To determine how ploidy level and sex-locus genotype affect gene expression during development, we compared expression patterns between diploid males, haploid males and females (queens) at three developmental timepoints in Solenopsis invicta. In pupae, gene expression profiles of diploid males were very different from those of haploid males but nearly identical to those of queens. An unexpected shift in expression patterns emerged soon after adult eclosion, with diploid male patterns diverging from those of queens to resemble those of haploid males, a pattern retained in older adults. The finding that ploidy level effects on early gene expression override sex effects (including genes implicated in sperm production and pheromone production/perception) may explain diploid male sterility and lack of worker discrimination against them during development.

Relaxed selection is a precursor to the evolution of phenotypic plasticity.

Phenotypic plasticity allows organisms to produce alternative phenotypes under different conditions and represents one of the most important ways by which organisms adaptively respond to the environment. However, the relationship between phenotypic plasticity and molecular evolution remains poorly understood. We addressed this issue by investigating the evolution of genes associated with phenotypically plastic castes, sexes, and developmental stages of the fire ant Solenopsis invicta. We first determined if genes associated with phenotypic plasticity in S. invicta evolved at a rapid rate, as predicted under theoretical models. We found that genes differentially expressed between S. invicta castes, sexes, and developmental stages all exhibited elevated rates of evolution compared with ubiquitously expressed genes. We next investigated the evolutionary history of genes associated with the production of castes. Surprisingly, we found that orthologs of caste-biased genes in S. invicta and the social bee Apis mellifera evolved rapidly in lineages without castes. Thus, in contrast to some theoretical predictions, our results suggest that rapid rates of molecular evolution may not arise primarily as a consequence of phenotypic plasticity. Instead, genes evolving under relaxed purifying selection may more readily adopt new forms of biased expression during the evolution of alternate phenotypes. These results suggest that relaxed selective constraint on protein-coding genes is an important and underappreciated element in the evolutionary origin of phenotypic plasticity.