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.

Genomic architecture and sexually dimorphic expression underlying immunity in the red mason bee, Osmia bicornis.

Insect pollinators provide crucial ecosystem services yet face increasing environmental pressures. The challenges posed by novel and reemerging pathogens on bee health means we need to improve our understanding of the immune system, an important barrier to infections and disease. Despite the importance of solitary bees, which are ecologically relevant, our understanding of the genomic basis and molecular mechanisms underlying their immune potential, and how intrinsic and extrinsic factors may influence it is limited. To improve our understanding of the genomic architecture underlying immunity of a key solitary bee pollinator, we characterized putative immune genes of the red mason bee, Osmia bicornis. In addition, we used publicly available RNA-seq datasets to determine how sexes differ in immune gene expression and splicing but also how pesticide exposure may affect immune gene expression in females. Through comparative genomics, we reveal an evolutionarily conserved set of more than 500 putative immune-related genes. We found genome-wide patterns of sex-biased gene expression, with greater enrichment of immune-related processes among genes with higher constitutive expression in males than females. Our results also suggest an up-regulation of immune-related genes in response to exposure to two common neonicotinoids, thiacloprid and imidacloprid. Collectively, our study provides important insights into the gene repertoire, regulation and expression differences in the sexes of O. bicornis, as well as providing additional support for how neonicotinoids can affect immune gene expression, which may affect the capacity of solitary bees to respond to pathogenic threats.

Mating precedes selective immune priming which is maintained throughout bumblebee queen diapause.

BACKGROUND: Understanding the mechanisms by which organisms adapt to unfavourable conditions is a fundamental question in ecology and evolutionary biology. One such mechanism is diapause, a period of dormancy typically found in nematodes, fish, crustaceans and insects. This state is a key life-history event characterised by arrested development, suppressed metabolism and increased stress tolerance and allows an organism to avoid prolonged periods of harsh and inhospitable environmental conditions. For some species, diapause is preceded by mating which can have a profound effect on female behaviour, physiology and key biological processes, including immunity. However, our understanding of how mating impacts long-term immunity and whether these effects persist throughout diapause is currently limited. To address this, we explored molecular changes in the haemolymph of the ecologically important pollinator, the buff-tailed bumblebee Bombus terrestris. B. terrestris queens mate prior to entering diapause, a non-feeding period of arrested development that can last 6-9 months. Using mass-spectrometry-based proteomics, we quantified changes in the pre-diapause queen haemolymph after mating, as well as the subsequent protein expression of mated queens during and post-diapause. RESULTS: Our analysis identified distinct proteome profiles associated with diapause preparation, maintenance and termination. More specifically, mating pre-diapause was followed by an increase in the abundance of antimicrobial peptides, key effectors of the immune system. Furthermore, we identified the elevated abundance of these proteins to be maintained throughout diapause. This finding was in contrast to the general reduction observed in immune proteins during diapause suggestive of selective immune priming and expression during diapause. Diapause also affected the expression of proteins involved in cuticular maintenance, olfaction, as well as proteins of unknown function, which may have roles in diapause regulation. CONCLUSIONS: Our results provide clear molecular evidence for the consequences and benefits of mating at the immune level as it precedes the selective increased abundance of antimicrobial peptides that are sustained throughout diapause. In addition, our results provide novel insights into the molecular mechanisms by which bumblebees prepare for, survive, and recover from diapause, insights that may have implications for our general understanding of these processes in other insect groups.

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.

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.

Signatures of Adaptation, Constraints, and Potential Redundancy in the Canonical Immune Genes of a Key Pollinator.

All organisms require an immune system to recognize, differentiate, and defend against pathogens. From an evolutionary perspective, immune systems evolve under strong selective pressures exerted by fast-evolving pathogens. However, the functional diversity of the immune system means that different immune components and their associated genes may evolve under varying forms of selection. Insect pollinators, which provide essential ecosystem services, are an important system in which to understand how selection has shaped immune gene evolution as their populations are experiencing declines with pathogens highlighted as a potential contributing factor. To improve our understanding of the genetic variation found in the immune genes of an essential pollinator, we performed whole-genome resequencing of wild-caught Bombus terrestris males. We first assessed nucleotide diversity and extended haplotype homozygosity for canonical immune genes finding the strongest signatures of positive selection acting on genes involved in pathogen recognition and antiviral defense, possibly driven by growing pathogen spread in wild populations. We also identified immune genes evolving under strong purifying selection, highlighting potential constraints on the bumblebee immune system. Lastly, we highlight the potential loss of function alleles present in the immune genes of wild-caught haploid males, suggesting that such genes are potentially less essential for development and survival and represent redundancy in the gene repertoire of the bumblebee immune system. Collectively, our analysis provides novel insights into the recent evolutionary history of the immune system of a key pollinator, highlighting targets of selection, constraints to adaptation, and potential redundancy.

Unexpected worker mating and colony-founding in a superorganism.

The emergence of caste-differentiated colonies, which have been defined as 'superorganisms', in ants, bees, and wasps represents a major transition in evolution. Lifetime mating commitment by queens, pre-imaginal caste determination and lifetime unmatedness of workers are key features of these animal societies. Workers in superorganismal species like honey bees and many ants have consequently lost, or retain only vestigial spermathecal structures. However, bumble bee workers retain complete spermathecae despite 25-40 million years since their origin of superorganismality, which remains an evolutionary mystery. Here, we show (i) that bumble bee workers retain queen-like reproductive traits, being able to mate and produce colonies, underlain by queen-like gene expression, (ii) the social conditions required for worker mating, and (iii) that these abilities may be selected for by early queen-loss in these annual species. These results challenge the idea of lifetime worker unmatedness in superorganisms, and provide an exciting new tool for the conservation of endangered bumble bee species.

Polyphenism in social insects: insights from a transcriptome-wide analysis of gene expression in the life stages of the key pollinator, Bombus terrestris.

BACKGROUND: Understanding polyphenism, the ability of a single genome to express multiple morphologically and behaviourally distinct phenotypes, is an important goal for evolutionary and developmental biology. Polyphenism has been key to the evolution of the Hymenoptera, and particularly the social Hymenoptera where the genome of a single species regulates distinct larval stages, sexual dimorphism and physical castes within the female sex. Transcriptomic analyses of social Hymenoptera will therefore provide unique insights into how changes in gene expression underlie such complexity. Here we describe gene expression in individual specimens of the pre-adult stages, sexes and castes of the key pollinator, the buff-tailed bumblebee Bombus terrestris. RESULTS: cDNA was prepared from mRNA from five life cycle stages (one larva, one pupa, one male, one gyne and two workers) and a total of 1,610,742 expressed sequence tags (ESTs) were generated using Roche 454 technology, substantially increasing the sequence data available for this important species. Overlapping ESTs were assembled into 36,354 B. terrestris putative transcripts, and functionally annotated. A preliminary assessment of differences in gene expression across non-replicated specimens from the pre-adult stages, castes and sexes was performed using R-STAT analysis. Individual samples from the life cycle stages of the bumblebee differed in the expression of a wide array of genes, including genes involved in amino acid storage, metabolism, immunity and olfaction. CONCLUSIONS: Detailed analyses of immune and olfaction gene expression across phenotypes demonstrated how transcriptomic analyses can inform our understanding of processes central to the biology of B. terrestris and the social Hymenoptera in general. For example, examination of immunity-related genes identified high conservation of important immunity pathway components across individual specimens from the life cycle stages while olfactory-related genes exhibited differential expression with a wider repertoire of gene expression within adults, especially sexuals, in comparison to immature stages. As there is an absence of replication across the samples, the results of this study are preliminary but provide a number of candidate genes which may be related to distinct phenotypic stage expression. This comprehensive transcriptome catalogue will provide an important gene discovery resource for directed programmes in ecology, evolution and conservation of a key pollinator.

Evolution and Expression of the Immune System of a Facultatively Anadromous Salmonid.

Vertebrates have evolved a complex immune system required for the identification of and coordinated response to harmful pathogens. Migratory species spend periods of their life-cycle in more than one environment, and their immune system consequently faces a greater diversity of pathogens residing in different environments. In facultatively anadromous salmonids, individuals may spend parts of their life-cycle in freshwater and marine environments. For species such as the brown trout Salmo trutta, sexes differ in their life-histories with females more likely to migrate to sea while males are more likely to stay and complete their life-cycle in their natal river. Salmonids have also undergone a lineage-specific whole genome duplication event, which may provide novel immune innovations but our current understanding of the differences in salmonid immune expression between the sexes is limited. We characterized the brown trout immune gene repertoire, identifying a number of canonical immune genes in non-salmonid teleosts to be duplicated in S. trutta, with genes involved in innate and adaptive immunity. Through genome-wide transcriptional profiling ("RNA-seq") of male and female livers to investigate sex differences in gene expression amplitude and alternative splicing, we identified immune genes as being generally male-biased in expression. Our study provides important insights into the evolutionary consequences of whole genome duplication events on the salmonid immune gene repertoire and how the sexes differ in constitutive immune expression.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes.

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Exposure to microplastics reduces attachment strength and alters the haemolymph proteome of blue mussels (Mytilus edulis).

The contamination of marine ecosystems with microplastics, such as the polymer polyethylene, a commonly used component of single-use packaging, is of global concern. Although it has been suggested that biodegradable polymers, such as polylactic acid, may be used to replace some polyethylene packaging, little is known about their effects on marine organisms. Blue mussels, Mytilus edulis, have become a "model organism" for investigating the effects of microplastics in marine ecosystems. We show here that repeated exposure, over a period of 52 days in an outdoor mesocosm setting, of M. edulis to polyethylene microplastics reduced the number of byssal threads produced and the attachment strength (tenacity) by ∼50%. Exposure to either type of microplastic altered the haemolymph proteome and, although a conserved response to microplastic exposure was observed, overall polyethylene resulted in more changes to protein abundances than polylactic acid. Many of the proteins affected are involved in vital biological processes, such as immune regulation, detoxification, metabolism and structural development. Our study highlights the utility of mass spectrometry-based proteomics to assess the health of key marine organisms and identifies the potential mechanisms by which microplastics, both conventional and biodegradable, could affect their ability to form and maintain reefs.

The liver proteome in a mouse model for Ascaris suum resistance and susceptibility: evidence for an altered innate immune response.

BACKGROUND: Ascariasis is a neglected tropical disease that affects 800 million people worldwide. Whereas most people only experience light worm burden, some people experience heavy worm burdens even after several rounds of chemotherapy, a phenomenon known as predisposition. Such heavy infections are associated with more severe symptoms and increased chronic morbidity. METHODS: In order to investigate potential mechanisms that may explain the observed predisposition, we infected mice with the porcine ascarid Ascaris suum using an established mouse model with two different mouse strains, where the C57BL/6J strain is more susceptible to infection and therefore a model for heavy infection and the CBA/Ca strain is more resistant and thus a model for light infection. At day 7 post-infection we investigated the liver proteome, using shotgun mass spectrometry, of both infected and control mice of each strain. RESULTS: We identified intrinsic differences, between the two mouse strains, in both oxidative phosphorylation proteins and proteins involved in retinol metabolism. Additionally, we found differences between the two mouse strains in activation of the complement system, where the CBA/Ca strain has higher protein abundances for lectin pathway proteins and the C57BL/6J strain has higher protein abundances for complement inhibiting proteins. The CBA/Ca strain had a higher abundance of proteins involved in the activation of the complement cascade via the lectin pathway. In contrast, the C57BL/6J strain demonstrated a higher abundance of proteins involved in arresting the complement pathway. CONCLUSIONS: We observed clear differences between the two mouse strains both intrinsically and under infection.

A Proteomic Investigation of Hepatic Resistance to Ascaris in a Murine Model.

The helminth Ascaris causes ascariasis in both humans and pigs. Humans, especially children, experience significant morbidity including respiratory complications, growth deficits and intestinal obstruction. Given that 800 million people worldwide are infected by Ascaris, this represents a significant global public health concern. The severity of the symptoms and associated morbidity are related to the parasite burden and not all hosts are infected equally. While the pathology of the disease has been extensively examined, our understanding of the molecular mechanisms underlying resistance and susceptibility to this nematode infection is poor. In order to investigate host differences associated with heavy and light parasite burden, an experimental murine model was developed utilising Ascaris-susceptible and -resistant mice strains, C57BL/6J and CBA/Ca, respectively, which experience differential burdens of migratory Ascaris larvae in the host lungs. Previous studies identified the liver as the site where this difference in susceptibility occurs. Using a label free quantitative proteomic approach, we analysed the hepatic proteomes of day four post infection C57BL/6J and CBA/Ca mice with and without Ascaris infection to identify proteins changes potentially linked to both resistance and susceptibility amongst the two strains, respectively. Over 3000 proteins were identified in total and clear intrinsic differences were elucidated between the two strains. These included a higher abundance of mitochondrial proteins, particularly those associated with the oxidative phosphorylation pathway and reactive oxygen species (ROS) production in the relatively resistant CBA/Ca mice. We hypothesise that the increased ROS levels associated with higher levels of mitochondrial activity results in a highly oxidative cellular environment that has a dramatic effect on the nematode's ability to successfully sustain a parasitic association with its resistant host. Under infection, both strains had increased abundances in proteins associated with the oxidative phosphorylation pathway, as well as the tricarboxylic acid cycle, with respect to their controls, indicating a general stress response to Ascaris infection. Despite the early stage of infection, some immune-associated proteins were identified to be differentially abundant, providing a novel insight into the host response to Ascaris. In general, the susceptible C57BL/6J mice displayed higher abundances in immune-associated proteins, most likely signifying a more active nematode cohort with respect to their CBA/Ca counterparts. The complement component C8a and S100 proteins, S100a8 and S100a9, were highly differentially abundant in both infected strains, signifying a potential innate immune response and the importance of the complement pathway in defence against macroparasite infection. In addition, the signatures of an early adaptive immune response were observed through the presence of proteins, such as plastin-2 and dipeptidyl peptidase 1. A marked decrease in proteins associated with translation was also observed in both C57BL/6J and CBA/Ca mice under infection, indicative of either a general response to Ascaris or a modulatory effect by the nematode itself. Our research provides novel insights into the in vivo host-Ascaris relationship on the molecular level and provides new research perspectives in the development of Ascaris control and treatment strategies.

A depauperate immune repertoire precedes evolution of sociality in bees.

BACKGROUND: Sociality has many rewards, but can also be dangerous, as high population density and low genetic diversity, common in social insects, is ideal for parasite transmission. Despite this risk, honeybees and other sequenced social insects have far fewer canonical immune genes relative to solitary insects. Social protection from infection, including behavioral responses, may explain this depauperate immune repertoire. Here, based on full genome sequences, we describe the immune repertoire of two ecologically and commercially important bumblebee species that diverged approximately 18 million years ago, the North American Bombus impatiens and European Bombus terrestris. RESULTS: We find that the immune systems of these bumblebees, two species of honeybee, and a solitary leafcutting bee, are strikingly similar. Transcriptional assays confirm the expression of many of these genes in an immunological context and more strongly in young queens than males, affirming Bateman's principle of greater investment in female immunity. We find evidence of positive selection in genes encoding antiviral responses, components of the Toll and JAK/STAT pathways, and serine protease inhibitors in both social and solitary bees. Finally, we detect many genes across pathways that differ in selection between bumblebees and honeybees, or between the social and solitary clades. CONCLUSIONS: The similarity in immune complement across a gradient of sociality suggests that a reduced immune repertoire predates the evolution of sociality in bees. The differences in selection on immune genes likely reflect divergent pressures exerted by parasites across social contexts.

The genomes of two key bumblebee species with primitive eusocial organization.

BACKGROUND: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. RESULTS: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. CONCLUSIONS: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.