Honeybee Genome Sequencing Consortium. Insights into social insects from the genome of the honeybee Apis mellifera

Effects of Functional Depletion of Doublesex on Male Development in the Sawfly, Athalia rosae

Simple Summary

The sawfly, Athalia rosae, exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. The doublesex (dsx) gene is a well-conserved transcription factor that regulates sexual differentiation in insects. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. As a result, knockdown of Ardsx in haploid males caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. The same knockdown approach using diploid males caused complete male-to-female sex reversal; they were able to produce fertile eggs and exhibited female behaviors. The same RNAi treatment did not affect female differentiation. These results demonstrated that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Abstract

The doublesex (dsx) gene, which encodes a transcription factor, regulates sexual differentiation in insects. Sex-specific splicing of dsx occurs to yield male- and female-specific isoforms, which promote male and female development, respectively. Thus, functional disruption of dsx leads to an intersexual phenotype in both sexes. We previously identified a dsx ortholog in the sawfly, Athalia rosae. Similar to dsx in other insects, dsx in the sawfly yields different isoforms in males and females as a result of alternative splicing. The sawfly exploits a haplodiploid mode of reproduction, in which fertilized eggs develop into diploid females, whereas unfertilized eggs parthenogenetically develop into haploid males. In the present study, we knocked down the A. rosae ortholog of dsx (Ardsx) during several developmental stages with repeated double-stranded RNA (dsRNA) injections. Knockdown of Ardsx via parental RNA interference (RNAi), which enables knockdown of genes in offspring embryos, led to a lack of internal and external genitalia in haploid male progeny. Additional injection of dsRNA targeting Ardsx in these animals caused almost complete male-to-female sex reversal, but the resulting eggs were infertile. Notably, the same knockdown approach using diploid males obtained by sib-crossing caused complete male-to-female sex reversal; they were morphologically and behaviorally females. The same RNAi treatment did not affect female differentiation. These results indicate that dsx in the sawfly is essential for male development and its depletion caused complete male-to-female sex reversal. This is the first demonstration of functional depletion of dsx not causing intersexuality but inducing total sex reversal in males instead.

Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes

Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.

Recombination mapping of the Brazilian stingless bee Frieseomelitta varia confirms high recombination rates in social hymenoptera

Background

Meiotic recombination is a fundamental genetic process that shuffles allele combinations and promotes accurate segregation of chromosomes. Analyses of the ubiquitous variation of recombination rates within and across species suggest that recombination is evolving adaptively. All studied insects with advanced eusociality have shown exceptionally high recombination rates, which may represent a prominent case of adaptive evolution of recombination. However, our understanding of the relationship between social evolution and recombination rates is incomplete, partly due to lacking empirical data. Here, we present a linkage map of the monandrous, advanced eusocial Brazilian stingless bee, Frieseomelitta varia, providing the first recombination analysis in the diverse Meliponini (Hymenoptera, Apidae).

Results

Our linkage map includes 1417 markers in 19 linkage groups. This map spans approximately 2580 centimorgans, and comparisons to the physical genome assembly indicate that it covers more than 75 % of the 275 Megabasepairs (Mbp) F. varia genome. Thus, our study results in a genome-wide recombination rate estimate of 9.3–12.5 centimorgan per Mbp. This value is higher than estimates from nonsocial insects and comparable to other highly social species, although it does not support our prediction that monandry and strong queen-worker caste divergence of F. varia lead to even higher recombination rates than other advanced eusocial species.

Conclusions

Our study expands the association between elevated recombination and sociality in the order Hymenoptera and strengthens the support for the hypothesis that advanced social evolution in hymenopteran insects invariably selects for high genomic recombination rates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07987-3.

Neofunctionalization of an ancient domain allows parasites to avoid intraspecific competition by manipulating host behaviour

Intraspecific competition is a major force in mediating population dynamics, fuelling adaptation, and potentially leading to evolutionary diversification. Among the evolutionary arms races between parasites, one of the most fundamental and intriguing behavioural adaptations and counter-adaptations are superparasitism and superparasitism avoidance. However, the underlying mechanisms and ecological contexts of these phenomena remain underexplored. Here, we apply the Drosophila parasite Leptopilina boulardi as a study system and find that this solitary endoparasitic wasp provokes a host escape response for superparasitism avoidance. We combine multi-omics and in vivo functional studies to characterize a small set of RhoGAP domain-containing genes that mediate the parasite's manipulation of host escape behaviour by inducing reactive oxygen species in the host central nervous system. We further uncover an evolutionary scenario in which neofunctionalization and specialization gave rise to the novel role of RhoGAP domain in avoiding superparasitism, with an ancestral origin prior to the divergence between Leptopilina specialist and generalist species. Our study suggests that superparasitism avoidance is adaptive for a parasite and adds to our understanding of how the molecular manipulation of host behaviour has evolved in this system.

Hermetia illucens (L.) (Diptera: Stratiomyidae) Odorant Binding Proteins and Their Interactions with Selected Volatile Organic Compounds: An In Silico Approach

The black soldier fly (BSF), Hermetia illucens (Diptera: Stratiomyidae), has considerable global interest due to its outstanding capacity in bioconverting organic waste to insect biomass, which can be used for livestock, poultry, and aquaculture feed. Mass production of this insect in colonies requires the development of methods concentrating oviposition in specific collection devices, while the mass production of larvae and disposing of waste may require substrates that are more palatable and more attractive to the insects. In insects, chemoreception plays an essential role throughout their life cycle, responding to an array of chemical, biological and environmental signals to locate and select food, mates, oviposition sites and avoid predators. To interpret these signals, insects use an arsenal of molecular components, including small proteins called odorant binding proteins (OBPs). Next generation sequencing was used to identify genes involved in chemoreception during the larval and adult stage of BSF, with particular attention to OBPs. The analysis of the de novo adult and larval transcriptome led to the identification of 27 and 31 OBPs for adults and larvae, respectively. Among these OBPs, 15 were common in larval and adult transcriptomes and the tertiary structures of 8 selected OBPs were modelled. In silico docking of ligands confirms the potential interaction with VOCs of interest. Starting from the information about the growth performance of H. illucens on different organic substrates from the agri-food sector, the present work demonstrates a possible correlation between a pool of selected VOCs, emitted by those substrates that are attractive for H. illucens females when searching for oviposition sites, as well as phagostimulants for larvae. The binding affinities between OBPs and selected ligands calculated by in silico modelling may indicate a correlation among OBPs, VOCs and behavioural preferences that will be the basis for further analysis.

Improved Apis mellifera reference genome based on the alternative long-read-based assemblies

Apis mellifera L., the western honey bee is a major crop pollinator that plays a key role in beekeeping and serves as an important model organism in social behavior studies. Recent efforts have improved on the quality of the honey bee reference genome and developed a chromosome-level assembly of 16 chromosomes, two of which are gapless. However, the rest suffer from 51 gaps, 160 unplaced/unlocalized scaffolds, and the lack of 2 distal telomeres. The gaps are located at the hard-to-assemble extended highly repetitive chromosomal regions that may contain functional genomic elements. Here, we use de novo re-assemblies from the most recent reference genome Amel_HAv_3.1 raw reads and other long-read-based assemblies (INRA_AMelMel_1.0, ASM1384120v1, and ASM1384124v1) of the honey bee genome to resolve 13 gaps, five unplaced/unlocalized scaffolds and, the lacking telomeres of the Amel_HAv_3.1. The total length of the resolved gaps is 848,747 bp. The accuracy of the corrected assembly was validated by mapping PacBio reads and performing gene annotation assessment. Comparative analysis suggests that the PacBio-reads-based assemblies of the honey bee genomes failed in the same highly repetitive extended regions of the chromosomes, especially on chromosome 10. To fully resolve these extended repetitive regions, further work using ultra-long Nanopore sequencing would be needed. Our updated assembly facilitates more accurate reference-guided scaffolding and marker/sequence mapping in honey bee genomics studies.

Societies to genes: can we get there from here?

Understanding the organization and evolution of social complexity is a major task because it requires building an understanding of mechanisms operating at different levels of biological organization from genes to social interactions. I discuss here, a unique forward genetic approach spanning more than 30 years beginning with human-assisted colony-level selection for a single social trait, the amount of pollen honey bees (Apis mellifera L.) store. The goal was to understand a complex social trait from the social phenotype to genes responsible for observed trait variation. The approach combined the results of colony-level selection with detailed studies of individual behavior and physiology resulting in a mapped, integrated phenotypic architecture composed of correlative relationships between traits spanning anatomy, physiology, sensory response systems, and individual behavior that affect individual foraging decisions. Colony-level selection reverse engineered the architecture of an integrated phenotype of individuals resulting in changes in the social trait. Quantitative trait locus (QTL) studies combined with an exceptionally high recombination rate (60 kb/cM), and a phenotypic map, provided a genotype-phenotype map of high complexity demonstrating broad QTL pleiotropy, epistasis, and epistatic pleiotropy suggesting that gene pleiotropy or tight linkage of genes within QTL integrated the phenotype. Gene expression and knockdown of identified positional candidates revealed genes affecting foraging behavior and confirmed one pleiotropic gene, a tyramine receptor, as a target for colony-level selection that was under selection in two different tissues in two different life stages. The approach presented here has resulted in a comprehensive understanding of the structure and evolution of honey bee social organization.

Sex-Specific Development in Haplodiploid Honeybee Is Controlled by the Female-Embryo-Specific Activation of Thousands of Intronic LncRNAs

Embryonic development depends on a highly coordinated shift in transcription programs known as the maternal-to-zygotic transition (MZT). It remains unclear how haploid and diploid embryo coordinate their genomic activation and embryonic development during MZT in haplodiploid animals. Here, we applied a single-embryo RNA-seq approach to characterize the embryonic transcriptome dynamics in haploid males vs. diploid females of the haplodiploid insect honeybee (Apis mellifera). We observed typical zygotic genome activation (ZGA) occurred in three major waves specifically in female honeybee embryos; haploid genome activation was much weaker and occurred later. Strikingly, we also observed three waves of transcriptional activation for thousands of long non-coding transcripts (lncRNA), 73% of which are transcribed from intronic regions and 65% were specific to female honeybee embryos. These findings support a model in which introns encode thousands of lncRNAs that are expressed in a diploid-embryo-specific and ZGA-triggered manner that may have potential functions to regulate gene expression during early embryonic development in the haplodiploid insect honeybee.

Elementary calcium release events in the skeletal muscle cells of the honey bee Apis mellifera

Calcium sparks are involved in major physiological and pathological processes in vertebrate muscles but have never been characterized in invertebrates. Here, dynamic confocal imaging on intact skeletal muscle cells isolated enzymatically from the adult honey bee legs allowed the first spatio-temporal characterization of subcellular calcium release events (CREs) in an insect species. The frequency of CREs, measured in x–y time lapse series, was higher than frequencies usually described in vertebrates. Honey bee CREs had a larger spatial spread at half maximum than their vertebrate counterparts and a slightly ellipsoidal shape, two characteristics that may be related to ultrastructural features specific to invertebrate cells. In line-scan experiments, the histogram of CREs’ duration followed a bimodal distribution, supporting the existence of both sparks and embers. Unlike in vertebrates, embers and sparks had similar amplitudes, a difference that could be related to genomic differences and/or excitation–contraction coupling specificities in honey bee skeletal muscle fibres. The first characterization of CREs from an arthropod which shows strong genomic, ultrastructural and physiological differences with vertebrates may help in improving the research field of sparkology and more generally the knowledge in invertebrates cell Ca²⁺ homeostasis, eventually leading to a better understanding of their roles and regulations in muscles but also the myotoxicity of new insecticides targeting ryanodine receptors.

BREC: an R package/Shiny app for automatically identifying heterochromatin boundaries and estimating local recombination rates along chromosomes

Background

Meiotic recombination is a vital biological process playing an essential role in genome's structural and functional dynamics. Genomes exhibit highly various recombination profiles along chromosomes associated with several chromatin states. However, eu-heterochromatin boundaries are not available nor easily provided for non-model organisms, especially for newly sequenced ones. Hence, we miss accurate local recombination rates necessary to address evolutionary questions.

Results

Here, we propose an automated computational tool, based on the Marey maps method, allowing to identify heterochromatin boundaries along chromosomes and estimating local recombination rates. Our method, called BREC (heterochromatin Boundaries and RECombination rate estimates) is non-genome-specific, running even on non-model genomes as long as genetic and physical maps are available. BREC is based on pure statistics and is data-driven, implying that good input data quality remains a strong requirement. Therefore, a data pre-processing module (data quality control and cleaning) is provided. Experiments show that BREC handles different markers' density and distribution issues.

Conclusions

BREC's heterochromatin boundaries have been validated with cytological equivalents experimentally generated on the fruit fly Drosophila melanogaster genome, for which BREC returns congruent corresponding values. Also, BREC's recombination rates have been compared with previously reported estimates. Based on the promising results, we believe our tool has the potential to help bring data science into the service of genome biology and evolution. We introduce BREC within an R-package and a Shiny web-based user-friendly application yielding a fast, easy-to-use, and broadly accessible resource. The BREC R-package is available at the GitHub repository https://github.com/GenomeStructureOrganization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04233-1.

Rapid Genomic Evolution Drives the Diversification of Male Reproductive Genes in Dung Beetles

The molecular basis for the evolution of novel phenotypes is a central question in evolutionary biology. In recent years, dung beetles have emerged as models for novel trait evolution as they possess distinct precopulatory traits such as sexually dimorphic horns on their head and thorax. Here, we use functional and evolutionary genomics to investigate the origins and the evolution of postcopulatory reproductive traits in male dung beetles. Male ejaculates that underlie postcopulatory sexual selection are excellent candidates to study novel trait evolution as they are complex, fast evolving, and often highly divergent in insects. We assemble de novo transcriptomes of male accessory glands and testes of a widespread dung beetle, Catharsius molossus, and we perform an evolutionary analysis of closely and distantly related insect genomes. Our results show there is rapid innovation at the genomic level even among closely related dung beetles. Genomic expansion and contraction drive the divergence of male reproductive traits and their functions. The birth of scores of completely novel reproductive genes is reinforced by the recruitment of these genes for high expression in male reproductive tissues, especially in the accessory glands. We find that male accessory glands of C. molossus are specialized for secretory function and express female, egg, and embryo-related genes as well as serine protease inhibitors, whilst the testes are specialized for spermatogenesis and sperm function. Finally, we touch upon putative functions of these evolutionary novelties using structure-function analysis as these proteins bear no homology to any other known proteins.

Identification and characterization of odorant binding proteins in the diamondback moth, Plutella xylostella

Odorant binding proteins (OBPs) are a group of soluble proteins functioning as odorant carriers in insect antennae, mouth parts and other chemosensory organs. However, multiple insect OBPs have been detected in other tissues and various functions have been proposed. Therefore, a detailed expression profile including stages, tissues and sexes where OBPs are expressed will assist in building the links to their potential functions, enhancing the functional studies of insect OBPs. Here, we identified 39 putative OBP genes from its genome and transcriptome sequences of diamondback moth (DBM), Plutella xylostella. The expression patterns of identified PxylOBPs were further investigated from eggs, larvae, pupae, virgin adults, mated adults, larval midgut, larval heads, adult antennae, adult heads and adult tarsi. Moreover, P. xylostella larvae and adults with and without host plants for 5 h were utilized to study the interactions between OBP expression and host plants. The results showed that most PxylOBPs were highly expressed in male and female adult antennae. The expression levels of certain PxyOBPs could be regulated by mating activities and feeding host plants. This study advances our knowledge of P. xylostella OBPs, which may help develop new strategies for more environmentally sustainable management of P. xylostella.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

A SNP assay for assessing diversity in immune genes in the honey bee (Apis mellifera L.)

With a growing number of parasites and pathogens experiencing large-scale range expansions, monitoring diversity in immune genes of host populations has never been so important because it can inform on the adaptive potential to resist the invaders. Population surveys of immune genes are becoming common in many organisms, yet they are missing in the honey bee (Apis mellifera L.), a key managed pollinator species that has been severely affected by biological invasions. To fill the gap, here we identified single nucleotide polymorphisms (SNPs) in a wide range of honey bee immune genes and developed a medium-density assay targeting a subset of these genes. Using a discovery panel of 123 whole-genomes, representing seven A. mellifera subspecies and three evolutionary lineages, 180 immune genes were scanned for SNPs in exons, introns (< 4 bp from exons), 3' and 5´UTR, and < 1 kb upstream of the transcription start site. After application of multiple filtering criteria and validation, the final medium-density assay combines 91 quality-proved functional SNPs marking 89 innate immune genes and these can be readily typed using the high-sample-throughput iPLEX MassARRAY system. This medium-density-SNP assay was applied to 156 samples from four countries and the admixture analysis clustered the samples according to their lineage and subspecies, suggesting that honey bee ancestry can be delineated from functional variation. In addition to allowing analysis of immunogenetic variation, this newly-developed SNP assay can be used for inferring genetic structure and admixture in the honey bee.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Neural and Molecular Mechanisms of Biological Embedding of Social Interactions

Animals operate in complex environments, and salient social information is encoded in the nervous system and then processed to initiate adaptive behavior. This encoding involves biological embedding, the process by which social experience affects the brain to influence future behavior. Biological embedding is an important conceptual framework for understanding social decision-making in the brain, as it encompasses multiple levels of organization that regulate how information is encoded and used to modify behavior. The framework we emphasize here is that social stimuli provoke short-term changes in neural activity that lead to changes in gene expression on longer timescales. This process, simplified-neurons are for today and genes are for tomorrow-enables the assessment of the valence of a social interaction, an appropriate and rapid response, and subsequent modification of neural circuitry to change future behavioral inclinations in anticipation of environmental changes. We review recent research on the neural and molecular basis of biological embedding in the context of social interactions, with a special focus on the honeybee.

Parity and time reversal elucidate both decision-making in empirical models and attractor scaling in critical Boolean networks

We present new applications of parity inversion and time reversal to the emergence of complex behavior from simple dynamical rules in stochastic discrete models. Our parity-based encoding of causal relationships and time-reversal construction efficiently reveal discrete analogs of stable and unstable manifolds. We demonstrate their predictive power by studying decision-making in systems biology and statistical physics models. These applications underpin a novel attractor identification algorithm implemented for Boolean networks under stochastic dynamics. Its speed enables resolving a long-standing open question of how attractor count in critical random Boolean networks scales with network size and whether the scaling matches biological observations. Via 80-fold improvement in probed network size (N = 16,384), we find the unexpectedly low scaling exponent of 0.12 ± 0.05, approximately one-tenth the analytical upper bound. We demonstrate a general principle: A system's relationship to its time reversal and state-space inversion constrains its repertoire of emergent behaviors.

The genomic basis of evolutionary differentiation among honey bees

In contrast to the western honey bee, Apis mellifera, other honey bee species have been largely neglected despite their importance and diversity. The genetic basis of the evolutionary diversification of honey bees remains largely unknown. Here, we provide a genome-wide comparison of three honey bee species, each representing one of the three subgenera of honey bees, namely the dwarf (Apis florea), giant (A. dorsata), and cavity-nesting (A. mellifera) honey bees with bumblebees as an outgroup. Our analyses resolve the phylogeny of honey bees with the dwarf honey bees diverging first. We find that evolution of increased eusocial complexity in Apis proceeds via increases in the complexity of gene regulation, which is in agreement with previous studies. However, this process seems to be related to pathways other than transcriptional control. Positive selection patterns across Apis reveal a trade-off between maintaining genome stability and generating genetic diversity, with a rapidly evolving piRNA pathway leading to genomes depleted of transposable elements, and a rapidly evolving DNA repair pathway associated with high recombination rates in all Apis species. Diversification within Apis is accompanied by positive selection in several genes whose putative functions present candidate mechanisms for lineage-specific adaptations, such as migration, immunity, and nesting behavior.

Conservation genomics reveals pesticide and pathogen exposure in the declining bumble bee Bombus terricola

In recent years, many pollinators have experienced large population declines, which threaten food security and the stability of natural ecosystems. Bumble bees are particularly important because their ability to “buzz” pollinate and tolerate cooler conditions make them critical pollinators for certain plants and regions. Here, we apply a conservation genomics approach to study the vulnerable Bombus terricola. We sequenced RNA from 30 worker abdomens, 18 of which were collected from agricultural sites and 12 of which were collected from nonagricultural sites. We found transcriptional signatures associated with exposure to insecticides, with gene expression patterns suggesting that bumble bees were exposed to neonicotinoids and/or fipronil—two compounds known to negatively impact bees. We also found transcriptional signatures associated with pathogen infections. In addition to the transcriptomic analysis, we carried out a metatranscriptomic analysis and detected five pathogens in the abdomens of workers, three of which are common in managed honey bee and bumble bee colonies. Our conservation genomics study provides functional support for the role of pesticides and pathogen spillover in the decline of B. terricola. We demonstrate that conservation genomics is an invaluable tool which allows researchers to quantify the effects of multiple stressors that impact pollinator populations in the wild.

Midgut transcriptome assessment of the cockroach-hunting wasp Ampulex compressa (Apoidea: Ampulicidae)

The emerald jewel wasp Ampulex compressa (Hymenoptera: Ampulicidae) is a solitary wasp that is widely known for its specialized hunting of cockroaches as larvae provision. Adult wasps mainly feed on pollen and nectar, while their larvae feed on the cockroachs’ body, first as ecto- and later as endoparsitoids. Little is known about the expression of digestive, detoxification and stress-response-related genes in the midgut of A. compressa, or about its transcriptional versatility between life stages. To identify gut-biased genes related to digestion, detoxification, and stress response, we explored the midgut transcriptome of lab-reared A. compressa, for both adults and larvae, by focusing on the top 100 significantly up- and down-regulated genes. From the top 100 significantly differentially expressed genes (DEGs), we identified 39 and 36 DEGs putatively related to digestion and detoxification in the adult wasps and larvae, respectively. The two carbohydrases alpha-glucosidase (containing an alpha-amylase domain) and glycosyl hydrolase family 31, as well as the two proteinases chymotrypsin and trypsin, revealed the highest gene diversity. We identified six significant DEGs related to detoxification, which comprise glutathione S-transferase, cytochrome P450s and UDP-glucuronosyltransferase. The gene expression levels that were significantly expressed in both life stages vary strongly between life stages, as found in genes encoding for chymotrypsin and trypsin or glycosyl hydrolases family 31. The number of genes related to alpha-glucosidase, glycosyl hydrolase family 31, and cytochrome P450s was found to be similar across nine reference hymenopteran species, except for the identified glycosyl hydrolase family 31 gene, which was absent in all reference bee species. Phylogenetic analyses of the latter candidate genes revealed that they cluster together with their homologous genes found in the reference hymenopteran species. These identified candidate genes provide a basis for future comparative genomic and proteomic studies on (ontogenetic) dietary transitions in Hymenoptera.

The foraging gene as a modulator of division of labour in social insects

The social ants, bees, wasps, and termites include some of the most ecologically-successful groups of animal species. Their dominance in most terrestrial environments is attributed to their social lifestyle, which enable their colonies to exploit environmental resources with remarkable efficiency. One key attribute of social insect colonies is the division of labour that emerges among the sterile workers, which represent the majority of colony members. Studies of the mechanisms that drive division of labour systems across diverse social species have provided fundamental insights into the developmental, physiological, molecular, and genomic processes that regulate sociality, and the possible genetic routes that may have led to its evolution from a solitary ancestor. Here we specifically discuss the conserved role of the foraging gene, which encodes a cGMP-dependent protein kinase (PKG). Originally identified as a behaviourally polymorphic gene that drives alternative foraging strategies in the fruit fly Drosophila melanogaster, changes in foraging expression and kinase activity were later shown to play a key role in the division of labour in diverse social insect species as well. In particular, foraging appears to regulate worker transitions between behavioural tasks and specific behavioural traits associated with morphological castes. Although the specific neuroethological role of foraging in the insect brain remains mostly unknown, studies in genetically tractable insect species indicate that PKG signalling plays a conserved role in the neuronal plasticity of sensory, cognitive and motor functions, which underlie behaviours relevant to division of labour, including appetitive learning, aggression, stress response, phototaxis, and the response to pheromones.

Genus-Wide Characterization of Bumblebee Genomes Provides Insights into Their Evolution and Variation in Ecological and Behavioral Traits

Bumblebees are a diverse group of globally important pollinators in natural ecosystems and for agricultural food production. With both eusocial and solitary life-cycle phases, and some social parasite species, they are especially interesting models to understand social evolution, behavior, and ecology. Reports of many species in decline point to pathogen transmission, habitat loss, pesticide usage, and global climate change, as interconnected causes. These threats to bumblebee diversity make our reliance on a handful of well-studied species for agricultural pollination particularly precarious. To broadly sample bumblebee genomic and phenotypic diversity, we de novo sequenced and assembled the genomes of 17 species, representing all 15 subgenera, producing the first genus-wide quantification of genetic and genomic variation potentially underlying key ecological and behavioral traits. The species phylogeny resolves subgenera relationships, whereas incomplete lineage sorting likely drives high levels of gene tree discordance. Five chromosome-level assemblies show a stable 18-chromosome karyotype, with major rearrangements creating 25 chromosomes in social parasites. Differential transposable element activity drives changes in genome sizes, with putative domestications of repetitive sequences influencing gene coding and regulatory potential. Dynamically evolving gene families and signatures of positive selection point to genus-wide variation in processes linked to foraging, diet and metabolism, immunity and detoxification, as well as adaptations for life at high altitudes. Our study reveals how bumblebee genes and genomes have evolved across the Bombus phylogeny and identifies variations potentially linked to key ecological and behavioral traits of these important pollinators.

Probing the evolutionary history of epigenetic mechanisms: what can we learn from marine diatoms

Recent progress made on epigenetic studies revealed the conservation of epigenetic features in deep diverse branching species including Stramenopiles, plants and animals. This suggests their fundamental role in shaping species genomes across different evolutionary time scales. Diatoms are a highly successful and diverse group of phytoplankton with a fossil record of about 190 million years ago. They are distantly related from other super-groups of Eukaryotes and have retained some of the epigenetic features found in mammals and plants suggesting their ancient origin. Phaeodactylum tricornutum and Thalassiosira pseudonana, pennate and centric diatoms, respectively, emerged as model species to address questions on the evolution of epigenetic phenomena such as what has been lost, retained or has evolved in contemporary species. In the present work, we will discuss how the study of non-model or emerging model organisms, such as diatoms, helps understand the evolutionary history of epigenetic mechanisms with a particular focus on DNA methylation and histone modifications.

Identification of quantitative trait loci associated with calmness and gentleness in honey bees using whole-genome sequences

The identification of quantitative trait loci (QTL) through genome-wide association studies (GWAS) is a powerful method for unravelling the genetic background of selected traits and improving early-stage predictions. In honey bees (Apis mellifera), past genetic analyses have particularly focused on individual queens and workers. In this study, we used pooled whole-genome sequences to ascertain the genetic variation of the entire colony. In total, we sampled 216 Apis mellifera mellifera and 28 Apis mellifera carnica colonies. Different experts subjectively assessed the gentleness and calmness of the colonies using a standardised protocol. Conducting a GWAS for calmness on 211 purebred A. m. mellifera colonies, we identified three QTL, on chromosomes 8, 6, and 12. The two first QTL correspond to LOC409692 gene, coding for a disintegrin and metalloproteinase domain-containing protein 10, and to Abscam gene, coding for a Dscam family member Abscam protein, respectively. The last gene has been reported to be involved in the domestication of A. mellifera. The third QTL is located 13 kb upstream of LOC102655631, coding for a trehalose transporter. For gentleness, two QTL were identified on chromosomes 4 and 3. They are located within gene LOC413669, coding for a lap4 protein, and gene LOC413416, coding for a bicaudal C homolog 1-B protein, respectively. The identified positional candidate genes of both traits mainly affect the olfaction and nervous system of honey bees. Further research is needed to confirm the results and to better understand the genetic and phenotypic basis of calmness and gentleness.

Peripheral taste detection in honey bees: What do taste receptors respond to?

Understanding the neural principles governing taste perception in species that bear economic importance or serve as research models for other sensory modalities constitutes a strategic goal. Such is the case of the honey bee (Apis mellifera), which is environmentally and socioeconomically important, given its crucial role as pollinator agent in agricultural landscapes and which has served as a traditional model for visual and olfactory neurosciences and for research on communication, navigation, and learning and memory. Here we review the current knowledge on honey bee gustatory receptors to provide an integrative view of peripheral taste detection in this insect, highlighting specificities and commonalities with other insect species. We describe behavioral and electrophysiological responses to several tastant categories and relate these responses, whenever possible, to known molecular receptor mechanisms. Overall, we adopted an evolutionary and comparative perspective to understand the neural principles of honey bee taste and define key questions that should be answered in future gustatory research centered on this insect.

qPCR assays with dual-labeled probes for genotyping honey bee variants associated with varroa resistance

Background

The varroa mite is one of the main causes of honey bee mortality. An important mechanism by which honey bees increase their resistance against this mite is the expression of suppressed mite reproduction. This trait describes the physiological inability of mites to produce viable offspring and was found associated with eight genomic variants in previous research.

Results

This paper presents the development and validation of high-throughput qPCR assays with dual-labeled probes for discriminating these eight single-nucleotide variants. Amplicon sequences used for assay validation revealed additional variants in the primer/probe binding sites in four out of the eight assays. As for two of these the additional variants interfered with the genotyping outcome supplementary primers and/or probes were developed. Inclusion of these primers and probes in the assay mixes allowed for the correct genotyping of all eight variants of interest within our bee population.

Conclusion

These outcomes underline the importance of checking for interfering variants in designing qPCR assays. Ultimately, the availability of this assay allows genotyping for the suppressed mite reproduction trait and paves the way for marker assisted selection in breeding programs.

Oxidative stress and senescence in social insects: a significant but inconsistent link?

The life-prolonging effects of antioxidants have long entered popular culture, but the scientific community still debates whether free radicals and the resulting oxidative stress negatively affect longevity. Social insects are intriguing models for analysing the relationship between oxidative stress and senescence because life histories differ vastly between long-lived reproductives and the genetically similar but short-lived workers. Here, we present the results of an experiment on the accumulation of oxidative damage to proteins, and a comparative analysis of the expression of 20 selected genes commonly involved in managing oxidative damage, across four species of social insects: a termite, two bees and an ant. Although the source of analysed tissue varied across the four species, our results suggest that oxidative stress is a significant factor in senescence and that its manifestation and antioxidant defenses differ among species, making it difficult to find general patterns. More detailed and controlled investigations on why responses to oxidative stress may differ across social species may lead to a better understanding of the relations between oxidative stress, antioxidants, social life history and senescence. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

The Power of Drosophila melanogaster for Modeling Neonicotinoid Effects on Pollinators and Identifying Novel Mechanisms

Neonicotinoids are the most widely used insecticides in the world and are implicated in the widespread population declines of insects including pollinators. Neonicotinoids target nicotinic acetylcholine receptors which are expressed throughout the insect central nervous system, causing a wide range of sub-lethal effects on non-target insects. Here, we review the potential of the fruit fly Drosophila melanogaster to model the sub-lethal effects of neonicotinoids on pollinators, by utilizing its well-established assays that allow rapid identification and mechanistic characterization of these effects. We compare studies on the effects of neonicotinoids on lethality, reproduction, locomotion, immunity, learning, circadian rhythms and sleep in D. melanogaster and a range of pollinators. We also highlight how the genetic tools available in D. melanogaster, such as GAL4/UAS targeted transgene expression system combined with RNAi lines to any gene in the genome including the different nicotinic acetylcholine receptor subunit genes, are set to elucidate the mechanisms that underlie the sub-lethal effects of these common pesticides. We argue that studying pollinators and D. melanogaster in tandem allows rapid elucidation of mechanisms of action, which translate well from D. melanogaster to pollinators. We focus on the recent identification of novel and important sublethal effects of neonicotinoids on circadian rhythms and sleep. The comparison of effects between D. melanogaster and pollinators and the use of genetic tools to identify mechanisms make a powerful partnership for the future discovery and testing of more specific insecticides.

Honey Bee Diversity Is Swayed by Migratory Beekeeping and Trade Despite Conservation Practices: Genetic Evidence for the Impact of Anthropogenic Factors on Population Structure

The intense admixture of honey bee (Apis melliferaL.) populations at a global scale is mostly attributed to the widespread migratory beekeeping practices and replacement of queens and colonies with non-native races or hybrids of different subspecies. These practices are also common in Anatolia and Thrace, but their influence on the genetic make-up of the five native subspecies of honey bees has not been explored. Here, we present an analysis of 30 microsatellite markers from honey bees from six different regions in Anatolia and Thrace (N= 250 samples), with the aim of comparing the impact of: (1) migratory beekeeping, (2) queen and colony trade, and (3) conservation efforts on the genetic structure of native populations. Populations exposed to migratory beekeeping showed less allegiance than stationary ones. We found genetic evidence for migratory colonies, acting as a hybrid zone mobile in space and time, becoming vectors of otherwise local gene combinations. The effect of honey bee trade leaves very high introgression levels in native honey bees. Despite their narrow geographic range, introgression occurs mainly with the highly commercial Caucasian bees. We also measured the direction and magnitude of gene flow associated with bee trade. A comparison between regions that are open and those closed to migratory beekeeping allowed the evaluation of conservation sites as centers with limited gene flow and demonstrated the importance of establishing such isolated regions. Despite evidence of gene flow, our findings confirm high levels of geographically structured genetic diversity in four subspecies of honey bees in Turkey and emphasize the need to develop policies to maintain this diversity. Our overall results are of interest to the wider scientific community studying anthropogenic effects on the population diversity of honey bees and other insects. Our findings on the effects of migratory beekeeping, replacement of queens and colonies have implications for the conservation of honey bees, other pollinators, and invertebrates, in general, and are informative for policy-makers and other stakeholders in Europe and beyond.

Long Non-Coding RNAs in Insects

Only a small subset of all the transcribed RNAs are used as a template for protein translation, whereas RNA molecules that are not translated play a very important role as regulatory non-coding RNAs (ncRNAs). Besides traditionally known RNAs (ribosomal and transfer RNAs), ncRNAs also include small non-coding RNAs (sncRNAs) and long non-coding RNAs (lncRNAs). The lncRNAs, which were initially thought to be junk, have gained a great deal attention because of their regulatory roles in diverse biological processes in animals and plants. Insects are the most abundant and diverse group of animals on this planet. Recent studies have demonstrated the role of lncRNAs in almost all aspects of insect development, reproduction, and genetic plasticity. In this review, we describe the function and molecular mechanisms of the mode of action of different insect lncRNAs discovered up to date.

Regulators and signalling in insect antimicrobial innate immunity: Functional molecules and cellular pathways

Insects possess an immune system that protects them from attacks by various pathogenic microorganisms that would otherwise threaten their survival. Immune mechanisms may deal directly with the pathogens by eliminating them from the host organism or disarm them by suppressing the synthesis of toxins and virulence factors that promote the invasion and destructive action of the intruder within the host. Insects have been established as outstanding models for studying immune system regulation because innate immunity can be explored as an integrated system at the level of the whole organism. Innate immunity in insects consists of basal immunity that controls the constitutive synthesis of effector molecules such as antimicrobial peptides, and inducible immunity that is activated after detection of a microbe or its product(s). Activation and coordination of innate immune defenses in insects involve evolutionary conserved immune factors. Previous research in insects has led to the identification and characterization of distinct immune signalling pathways that modulate the response to microbial infections. This work has not only advanced the field of insect immunology, but it has also rekindled interest in the innate immune system of mammals. Here we review the current knowledge on key molecular components of insect immunity and discuss the opportunities they present for confronting infectious diseases in humans.

The trophocytes and oenocytes of worker and queen honey bees (Apis mellifera) exhibit distinct age-associated transcriptome profiles

Despite the identical genomic context, trophocytes and oenocytes in worker bees exhibit aging-related phenotypes, in contrast to the longevity phenotypes in queen bees. To explore this phenomenon at the molecular level, we evaluated the age-associated transcriptomes of trophocytes and oenocytes in worker bees and queen bees using high-throughput RNA-sequencing technology (RNA-seq). The results showed that (i) while gene expression profiles were different between worker and queen bees, they remained similar between young and old counterparts; (ii) worker bees express a high proportion of low-abundance genes, whereas queen bee transcriptomes display a high proportion of moderate-expression genes; (iii) genes were upregulated to a greater extent in queen bees vs. worker bees; and (iv) distinct aging-related and longevity-related candidate genes were found in worker and queen bees. These results provide new insights into the cellular aging and longevity of trophocytes and oenocytes in honey bees. Identification of aging-associated biomarker genes also constitutes a basis for translational research of aging in higher organisms.

Novel structure in the nuclei of honey bee brain neurons revealed by immunostaining

In the course of a screen designed to produce antibodies (ABs) with affinity to proteins in the honey bee brain we found an interesting AB that detects a highly specific epitope predominantly in the nuclei of Kenyon cells (KCs). The observed staining pattern is unique, and its unfamiliarity indicates a novel previously unseen nuclear structure that does not colocalize with the cytoskeletal protein f-actin. A single rod-like assembly, 3.7-4.1 µm long, is present in each nucleus of KCs in adult brains of worker bees and drones with the strongest immuno-labelling found in foraging bees. In brains of young queens, the labelling is more sporadic, and the rod-like structure appears to be shorter (~ 2.1 µm). No immunostaining is detectable in worker larvae. In pupal stage 5 during a peak of brain development only some occasional staining was identified. Although the cellular function of this unexpected structure has not been determined, the unusual distinctiveness of the revealed pattern suggests an unknown and potentially important protein assembly. One possibility is that this nuclear assembly is part of the KCs plasticity underlying the brain maturation in adult honey bees. Because no labelling with this AB is detectable in brains of the fly Drosophila melanogaster and the ant Camponotus floridanus, we tentatively named this antibody AmBNSab (Apis mellifera Brain Neurons Specific antibody). Here we report our results to make them accessible to a broader community and invite further research to unravel the biological role of this curious nuclear structure in the honey bee central brain.

RNA Interference-Mediated Knockdown of Genes Encoding Spore Wall Proteins Confers Protection against Nosema ceranae Infection in the European Honey Bee, Apis mellifera

Nosema ceranae (Opisthosporidia: Microsporidia) is an emergent intracellular parasite of the European honey bee (Apis mellifera) and causes serious Nosema disease which has been associated with worldwide honey bee colony losses. The only registered treatment for Nosema disease is fumagillin-b, and this has raised concerns about resistance and off-target effects. Fumagillin-B is banned from use in honey bee colonies in many countries, particularly in Europe. As a result, there is an urgent need for new and effective therapeutic options to treat Nosema disease in honey bees. An RNA interference (RNAi)-based approach can be a potent strategy for controlling diseases in honey bees. We explored the therapeutic potential of silencing the sequences of two N. ceranae encoded spore wall protein (SWP) genes by means of the RNAi-based methodology. Our study revealed that the oral ingestion of dsRNAs corresponding to SWP8 and SWP12 used separately or in combination could lead to a significant reduction in spore load, improve immunity, and extend the lifespan of N. ceranae-infected bees. The results from the work completed here enhance our understanding of honey bee host responses to microsporidia infection and highlight that RNAi-based therapeutics are a promising treatment for honey bee diseases.

Disentangling Ethiopian Honey Bee (Apis mellifera) Populations Based on Standard Morphometric and Genetic Analyses

Simple Summary

We conducted this population study of Ethiopian honey bees, using morphometric and genetic methods, to decipher their controversial classification. These honey bees are highly diverse and showed differentiation based on size and genetic information according to prevailing agro-ecological conditions, demonstrating morphological and molecular signatures of local adaptation. The results of both morphometric and genetic analyses suggest that Ethiopian honey bees differ from populations in the neighboring geographic regions and are characterized by extensive gene flow within the country, enhanced by honey bee colony trade. Consequently, future research that includes studying traits of vitality, behavior and colony performance of honey bees in remaining pocket areas of highland agro-ecological zones could contribute to the development of appropriate conservation management.

Abstract

The diversity and local differentiation of honey bees are subjects of broad general interest. In particular, the classification of Ethiopian honey bees has been a subject of debate for decades. Here, we conducted an integrated analysis based on classical morphometrics and a putative nuclear marker (denoted r7-frag) for elevational adaptation to classify and characterize these honey bees. Therefore, 660 worker bees were collected out of 66 colonies from highland, midland and lowland agro-ecological zones (AEZs) and were analyzed in reference to populations from neighboring countries. Multivariate morphometric analyses show that our Ethiopian samples are separate from Apis mellifera scutellata, A. m. jemenitica, A. m. litorea and A. m. monticola, but are closely related to A. m. simensis reference. Linear discriminant analysis showed differentiation according to AEZs in the form of highland, midland and lowland ecotypes. Moreover, size was positively correlated with elevation. Similarly, our Ethiopian samples were differentiated from A. m. monticola and A. m. scutellata based on r7-frag. There was a low tendency towards genetic differentiation between the Ethiopian samples, likely impacted by increased gene flow. However, the differentiation slightly increased with increasing elevational differences, demonstrated by the highland bees that showed higher differentiation from the lowland bees (F_ST = 0.024) compared to the midland bees (F_ST = 0.015). An allelic length polymorphism was detected (denoted as d) within r7-frag, showing a patterned distribution strongly associated with AEZ (X² = 11.84, p < 0.01) and found predominantly in highland and midland bees of some pocket areas. In conclusion, the Ethiopian honey bees represented in this study are characterized by high gene flow that suppresses differentiation.

Multifunctionality and intrinsic disorder of royal jelly proteome

Royal Jelly (RJ) is a gelatinous white-yellowish fluid, possessing a sour taste and a slight phenolic smell that is secreted by the hypopharyngeal and mandibular salivary glands of the nurse honeybees, and is used in nutrition of larvae and adult queens. Similar to other substances associated with the activities of honeybees, RJ not only contains nutritive components, such as carbohydrates, proteins, peptides, lipids, vitamins, and mineral salts, but also represents a natural ingredient with cosmetic and health-promoting properties. RJ is characterized by remarkable multifunctionality, possessing numerous biological activities. Although this multifunctionality of RJ can be considered as a consequence of its complex nature, many proteins and peptides in RJ are polyfunctional entities themselves. In this article, we show that RJ proteins contain different levels of intrinsic disorder, have sites of post-translational modifications, can be found in multiple isoforms, and many of them possess disorder-based binding sites, suggesting that the conformational ensembles of the RJ proteins might undergo change as a result of their interaction with specific binding partners. All these observations suggest that the multifunctionality of proteins and peptides from RJ is determined by their structural heterogeneity and polymorphism, and serve as an illustration of the protein structure-function continuum concept.

Rational design of Lepidoptera-specific insecticidal inhibitors targeting farnesyl diphosphate synthase, a key enzyme of the juvenile hormone biosynthetic pathway

Reducing the use of broad-spectrum insecticides is one of the many challenges currently faced by insect pest management practitioners. For this reason, efforts are being made to develop environmentally benign pest-control products through bio-rational approaches that aim at disrupting physiological processes unique to specific groups of pests. Perturbation of hormonal regulation of insect development and reproduction is one such strategy. It has long been hypothesized that some enzymes in the juvenile hormone biosynthetic pathway of moths, butterflies and caterpillars (order Lepidoptera) display unique structural features that could be targeted for the development of Lepidoptera-specific insecticides, a promising avenue given the numerous agricultural and forest pests belonging to this order. Farnesyl diphosphate synthase, FPPS, is one such enzyme, with recent work suggesting that it has structural characteristics that may enable its selective inhibition. This review synthesizes current knowledge on FPPS and summarizes recent advances in its use as a target for insecticide development.

Assessment of risk to honey bees and honey consumers resulting from the insect exposure to captan, thiacloprid, penthiopyrad, and λ-cyhalothrin used in a commercial apple orchard

Samples of leaves, flowers, soil, pollen, bee workers, bee brood, honey, and beeswax were collected to assess the possibility of a transfer of captan, thiacloprid, penthiopyrad, and λ-cyhalothrin from apple trees of Idared variety to honey bee (Apis mellifera) hives. Chemical analyses were performed using the Agilent 7890 Gas Chromatograph equipped with the Micro-cell Electron Capture Detector. It was found that significant amounts of penthiopyrad, the active ingredient of Fontelis 200 SC, were present in leaves, flowers, pollen, bee workers, and beeswax. Simultaneously, captan was present in the brood, worker bees, and honey samples. Significant levels of the captan residues were also detected on the soil surface. In honey samples, captan residue levels exceeded the acceptable standard, reaching 160% of its maximum residue level. However, in no case the amounts of captan, thiacloprid, penthiopyrad, and λ-cyhalothrin ingested with honey by an adult consumer exceeded the level of 0.02% of the acceptable daily intake. Despite the trace amounts of pesticide residues in honey samples collected during the field trial, bee honey consumption can be considered safe. An adult consumer can safely consume about 16 kg of honey.

Metatranscriptome Analysis of Sympatric Bee Species Identifies Bee Virus Variants and a New Virus, Andrena-Associated Bee Virus-1

Bees are important plant pollinators in agricultural and natural ecosystems. High average annual losses of honey bee (Apis mellifera) colonies in some parts of the world, and regional population declines of some mining bee species (Andrena spp.), are attributed to multiple factors including habitat loss, lack of quality forage, insecticide exposure, and pathogens, including viruses. While research has primarily focused on viruses in honey bees, many of these viruses have a broad host range. It is therefore important to apply a community level approach in studying the epidemiology of bee viruses. We utilized high-throughput sequencing to evaluate viral diversity and viral sharing in sympatric, co-foraging bees in the context of habitat type. Variants of four common viruses (i.e., black queen cell virus, deformed wing virus, Lake Sinai virus 2, and Lake Sinai virus NE) were identified in honey bee and mining bee samples, and the high degree of nucleotide identity in the virus consensus sequences obtained from both taxa indicates virus sharing. We discovered a unique bipartite + ssRNA Tombo-like virus, Andrena-associated bee virus-1 (AnBV-1). AnBV-1 infects mining bees, honey bees, and primary honey bee pupal cells maintained in culture. AnBV-1 prevalence and abundance was greater in mining bees than in honey bees. Statistical modeling that examined the roles of ecological factors, including floral diversity and abundance, indicated that AnBV-1 infection prevalence in honey bees was greater in habitats with low floral diversity and abundance, and that interspecific virus transmission is strongly modulated by the floral community in the habitat. These results suggest that land management strategies that aim to enhance floral diversity and abundance may reduce AnBV-1 spread between co-foraging bees.

Transcriptomic Responses of the Honey Bee Brain to Infection with Deformed Wing Virus

Managed colonies of European honey bees (Apis mellifera) are under threat from Varroa destructor mite infestation and infection with viruses vectored by mites. In particular, deformed wing virus (DWV) is a common viral pathogen infecting honey bees worldwide that has been shown to induce behavioral changes including precocious foraging and reduced associative learning. We investigated how DWV infection of bees affects the transcriptomic response of the brain. The transcriptomes of individual brains were analyzed using RNA-Seq after experimental infection of newly emerged adult bees with DWV. Two analytical methods were used to identify differentially expressed genes from the ~15,000 genes in the Apis mellifera genome. The 269 genes that had increased expression in DWV infected brains included genes involved in innate immunity such as antimicrobial peptides (AMPs), Ago2, and Dicer. Single bee brain NMR metabolomics methodology was developed for this work and indicates that proline is strongly elevated in DWV infected brains, consistent with the increased presence of the AMPs abaecin and apidaecin. The 1361 genes with reduced expression levels includes genes involved in cellular communication including G-protein coupled, tyrosine kinase, and ion-channel regulated signaling pathways. The number and function of the downregulated genes suggest that DWV has a major impact on neuron signaling that could explain DWV related behavioral changes.

Caste-specific gene expression underlying the differential adult brain development in the honeybee Apis mellifera

Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynβ). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories.

Identification of a new P450s gene (AccCYP4AV1) and its roles in abiotic stress resistance in the Apis cerana cerana Fabricius

Cytochrome P450 monooxygenases (P450s) play significant roles in protecting organisms from abiotic stress damage. Here, we report the sequence and characterization of a P450s gene (AccCYP4AV1), isolated from Apis cerana cerana Fabricius. The open reading frame of AccCYP4AV1 is 1506 base pairs long and encodes a predicted protein of 501 amino acids and 57.84 kDa, with an isoelectric point of 8.67. Real-time quantitative polymerase chain reaction (RT-qPCR) analysis indicated that AccCYP4AV1 is more highly expressed in the midgut than in other tissues. In addition, the highest expression occurs in newly emerged adult workers, followed by the first instar of the larval stage. In addition, the expression of the AccCYP4AV1 was upregulated by low temperature (4 °C), ultraviolet radiation, hydrogen peroxide, paraquat, and dichlorvos treatments. In contrast, AccCYP4AV1 transcription was downregulated by other abiotic stress conditions: exposure to increased temperature (44 °C), deltamethrin, cadmium chloride, and mercury (II) chloride. Moreover, when AccCYP4AV1 was knocked-down by RNA interference, the results suggested that multiple antioxidant genes (AccsHSP22.6, AccSOD2, AccTpx1, and AccTpx4) were downregulated and antioxidant genes AccGSTO1 and AccTrx1 were upregulated. The activity levels of peroxidase and catalase were upregulated in the AccCYP4AV1-knocked-down samples, compared with those in the control groups. These findings suggest that the AccCYP4AV1 protein might be involved in the defense against abiotic stress damage.

Tissue-specific transcription patterns support the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Kin selection may act differently on genes inherited from parents (matrigenes and patrigenes), resulting in intragenomic conflict. This conflict can be observed as differential expression of matrigenes and patrigenes, or parent-specific gene expression (PSGE). In honey bees (Apis mellifera), intragenomic conflict is hypothesized to occur in multiple social contexts. Previously, we found that patrigene-biased expression in reproductive tissues was associated with increased reproductive potential in worker honey bees, consistent with the prediction that patrigenes are selected to promote selfish behaviour in this context. Here, we examined brain gene expression patterns to determine if PSGE is also found in other tissues. As before, the number of transcripts showing patrigene expression bias was significantly greater in the brains of reproductive vs. sterile workers, while the number of matrigene-biased transcripts was not significantly different. Twelve transcripts out of the 374 showing PSGE in either tissue showed PSGE in both brain and reproductive tissues; this overlap was significantly greater than expected by chance. However, the majority of transcripts show PSGE only in one tissue, suggesting the epigenetic mechanisms mediating PSGE exhibit plasticity between tissues. There was no significant overlap between transcripts that showed PSGE and transcripts that were significantly differentially expressed. Weighted gene correlation network analysis identified modules which were significantly enriched in both types of transcripts, suggesting that these genes may influence each other through gene networks. Our results provide further support for the kin selection theory of intragenomic conflict, and provide valuable insights into the mechanisms which may mediate this process.

Chromosome-level genome assembly of the aphid parasitoid Aphidius gifuensis using Oxford Nanopore sequencing and Hi-C technology

Aphidius gifuensis is a parasitoid wasp that has been commercially bred and released in large scale as a biocontrol agent for the management of aphid pests. As a highly efficient endoparasitoid, it is also an important model for exploring mechanisms of parasitism. Currently, artificially bred populations of this wasp are facing rapid decline with undetermined cause, and mechanisms underlying its parasitoid strategy remain poorly understood. Exploring the mechanism behind its population decline and the host–parasitoid relationship is impeded partly due to the lack of a comprehensive genome data for this species. In this study, we constructed a high‐quality reference genome of A. gifuensis using Oxford Nanopore sequencing and Hi‐C (proximity ligation chromatin conformation capture) technology. The final genomic assembly was 156.9 Mb, with a contig N50 length of 3.93 Mb, the longest contig length of 10.4 Mb and 28.89% repetitive sequences. 99.8% of genome sequences were anchored onto six linkage groups. A total of 11,535 genes were predicted, of which 90.53% were functionally annotated. Benchmarking Universal Single‐Copy Orthologs (BUSCO) analysis showed the completeness of assembled genome is 98.3%. We found significantly expanded gene families involved in metabolic processes, transmembrane transport, cell signal communication and oxidoreductase activity, in particular ATP‐binding cassette (ABC) transporter, Cytochrome P450 and venom proteins. The olfactory receptors (ORs) showed significant contraction, which may be associated with the decrease in host recognition. Our study provides a solid foundation for future studies on the molecular mechanisms of population decline as well as host–parasitoid relationship for parasitoid wasps.

Two novel venom proteins underlie divergent parasitic strategies between a generalist and a specialist parasite

Parasitoids are ubiquitous in natural ecosystems. Parasitic strategies are highly diverse among parasitoid species, yet their underlying genetic bases are poorly understood. Here, we focus on the divergent adaptation of a specialist and a generalist drosophilid parasitoids. We find that a novel protein (Lar) enables active immune suppression by lysing the host lymph glands, eventually leading to successful parasitism by the generalist. Meanwhile, another novel protein (Warm) contributes to a passive strategy by attaching the laid eggs to the gut and other organs of the host, leading to incomplete encapsulation and helping the specialist escape the host immune response. We find that these diverse parasitic strategies both originated from lateral gene transfer, followed with duplication and specialization, and that they might contribute to the shift in host ranges between parasitoids. Our results increase our understanding of how novel gene functions originate and how they contribute to host adaptation.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

A Combined Strategy of Additive Manufacturing to Support Multidisciplinary Education in Arts, Biology, and Engineering

This research presents results for the design and creation of supporting teaching materials using additive manufacturing. The materials are inspired by selected artwork of four animal species, which belong to a collection from the museum of the University of Antioquia. The topic selected was fauna in Colombia, and the animals in question were chosen based on important roles they have in areas like health, the environment, and food. These animals will complement science education given to several age groups visiting the museum. In addition to the 3D-manufactured objects, a study was conducted using several age groups that are very relevant to the museum: children, teenagers, and first year undergraduate students. A video showing technical information cards about the manufacturing process was also developed. This project was multidisciplinary, involving collaboration between the engineering school, the museum, and a local high school. The results showed that young visitors want complete information on the animals and to have interaction with the animal models, which is not always possible. This project serves as a local strategy not only for taking arts and knowledge out of the museum but also for planning first year courses in the university and thus reducing problems like school dropout, low motivation, and poor performance in national exams.