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

Hydroxylation patterns associated with pheromone synthesis and composition in two honey bee subspecies Apis mellifera scutellata and A. m. capensis laying workers

Colony losses due to social parasitism in the form of reproductive workers of the Apis mellifera capensis clones results from the production of queen-like pheromonal signals coupled with ovarian activation in these socially parasitic honey bees. While the behavioral attributes of these social parasites have been described, their genetic attributes require more detailed exploration. Here, we investigate the production of mandibular gland pheromones in queenless workers of two sub-species of African honey bees; A. m. scutellata (low reproductive potential) and A. m. capensis clones (high reproductive potential). We used standard techniques in gas chromatography to assess the amounts of various pheromone components present, and qPCR to assess the expression of cytochrome P450 genes cyp6bd1 and cyp6as8, thought to be involved in the caste-dependent hydroxylation of acylated stearic acid in queens and workers, respectively. We found that, for both subspecies, the quality and quantity of the individual pheromone components vary with age, and that from the onset, A. m. capensis parasites make use of gene pathways typically upregulated in queens in achieving reproductive dominance. Due to the high production of 9-hydroxy-decenoic acid (9-HDA) the precursor to the queen substance 9-oxo-decenoic acid (9-ODA) in newly emerged capensis clones, we argue that clones are primed for parasitism upon emergence and develop into fully fledged parasites depending on the colony's social environment.

New approaches to old problems: Removal of phospholipase A2 results in highly active microsomal membranes from the honey bee, Apis mellifera

Over the last 50 years numerous studies were published by insect toxicologists using native microsomal membrane preparations in order to investigate in vitro cytochrome P450-(P450) mediated oxidative metabolism of xenobiotics, including insecticides. Whereas the preparation of active microsomal membranes from many pest insect species is straightforward, their isolation from honey bees, Apis mellifera (Hymenoptera: Apidae) remained difficult, if not impossible, due to the presence of a yet unidentified endogenous inhibitory factor released during abdominal gut membrane isolation. Thus hampering in vitro toxicological studies on microsomal oxidative phase 1 metabolism of xenobiotics, including compounds of ecotoxicological concern. The use of microsomal membranes rather than individually expressed P450s offers advantages and allows to develop a better understanding of phase 1 driven metabolic fate of foreign compounds. Here we biochemically investigated the problems associated with the isolation of active honey bee microsomes and developed a method resulting in highly active native microsomal preparations from adult female worker abdomens. This was achieved by removal of the abdominal venom gland sting complex prior to microsomal membrane preparation. Molecular sieve chromatography of the venom sac content leads to the identification of phospholipase A₂ as the enzyme responsible for the immediate inhibition of cytochrome P450 activity in microsomal preparations. The substrate specificity of functional honey bee microsomes was investigated with different fluorogenic substrates, and revealed a strong preference for coumarin over resorufin derivatives. Furthermore we were able to demonstrate the metabolism of insecticides by honey bee microsomes using an approach coupled to LC-MS/MS analysis of hydroxylated metabolites. Our work provides access to a new and simple in vitro tool to study honey bee phase 1 metabolism of xenobiotics utilising the entire range of microsomal cytochrome P450s.

Comparative Analysis of Intra- and Inter-Specific Genomic Variability in the Peach Potato Aphid, Myzus persicae

The availability of genomic data in the last decade relating to different aphid species has allowed the analysis of the genomic variability occurring among such species, whereas intra-specific variability has hitherto very largely been neglected. In order to analyse the intra-genomic variability in the peach potato aphid, Myzus persicae, comparative analyses were performed revealing several clone-specific gene duplications, together with numerous deletions/rearrangements. Our comparative approach also allowed us to evaluate the synteny existing between the two M. persicae clones tested and between the peach potato aphid and the pea aphid, Acyrthosiphon pisum. Even if part of the observed rearrangements are related to a low quality of some assembled contigs and/or to the high number of contigs present in these aphid genomes, our evidence reveals that aphid clones are genetically more different than expected. These results suggest that the choice of performing genomes sequencing combining different biotypes/populations, as revealed in the case of the soybean aphid, Aphis glycines, is unlikely to be very informative in aphids. Interestingly, it is possible that the holocentric nature of aphid chromosomes favours genome rearrangements that can be successively inherited transgenerationally via the aphid's apomictic (parthenogenetic) mode of reproduction. Lastly, we evaluated the structure of the cluster of genes coding for the five histones (H1, H2A, H2B, H3 and H4) in order to better understand the quality of the two M. persicae genomes and thereby to improve our knowledge of this functionally important gene family.

Honey Bees and Environmental Stress: Toxicologic Pathology of a Superorganism

As a eusocial species, Apis mellifera, the European honey bee, is effectively a superorganism-a group of genetically related individuals functioning as a collective unit. Because the unit of selection is the colony and not the individual, standard methods for assessing toxicologic pathology can miss colony-level responses to stress. For over a decade, US populations of honeybees have experienced severe annual losses attributed to a variety of environmental stressors varying temporally and geographically; differentiating among those stressors is accordingly a high priority. Social interactions among individuals in this social species, however, mean that the "footprint" of stressors such as pesticides, phytochemicals, pathogens, and parasites may be most discernible in individuals that did not themselves directly encounter the stressor. For example, neurotoxic effects of pesticides on nurse bees may impair their behavioral responses to queen-destined larvae, which may then emerge as adults with altered anatomy or physiology. Similarly, pesticide-induced size alterations in nurse hypopharyngeal glands, which produce royal jelly, the exclusive food of larval and adult queens, may disproportionately affect the queen's (and thus colony) health. Thus, evaluating toxicologic pathology in the honeybee requires a new perspective and development of assays that preserve the social context that ultimately determines colony health.

Analysis of the extent of synteny and conservation in the gene order in aphids: A first glimpse from the Aphis glycines genome

In the last decade several insect genomes have been sequenced, but for most the chromosomal mapping of the identified scaffolds/annotated genes is not available. The lack of this information makes it difficult to analyse various genetic aspects, including the presence of genome rearrangements and the extent of synteny within and across species. We mapped five multigenic DNA families (major and minor rDNAs, histone gene cluster, esterases and carotenoid desaturases) and seven scaffolds corresponding to 9 Mb of the soybean aphid, Aphis glycines, genome and identified loci spanning the four soybean aphid chromosomes. A comparative analysis of the localization of the annotated A. glycines genes with respect to the peach potato aphid, Myzus persicae, and the fly, Drosophila melanogaster, evidenced a lower degree of synteny between the two aphid species than in the aphid-fly comparison. Only 1.4 genes per syntenic block were observed in aphids in contrast to 2.3 genes per block in flies. This higher chromosomal rearrangement rate in aphids could be explained considering that they possess holocentric chromosomes that can favour the stabilization and inheritance of chromosomal rearrangements. Lastly, our experiments did not detect the presence of chimeric assemblies in the newly available A. glycines biotype 1 genome, differently from what reported in assembled genome of other aphid species, suggesting that chromosomal mapping can be used to ascertain the quality of assembled genomes.

RNA-seq reveals disruption of gene regulation when honey bees are caged and deprived of hive conditions

In this study, we present phenotypic and genetic data characterizing the impact of imidacloprid and caging stress on honey bee Apis mellifera physiological responses and regulation of 45 genes using targeted-RNA seq. The term 'caging stress' characterizes the effects of depriving honey bees of all hive aspects and conditions. Two cohorts of 1 day old sister bees were subjected to different conditions. One cohort was caged and fed different imidacloprid-tainted sugar solutions and the second was marked and introduced back to its natal hive. Physiological bee parameters and diet behavior were monitored daily for caged bees over several weeks. Bee samples from both cohorts were sampled weekly for RNA sequencing and oxidative stress analyses. Imidacloprid induced significant protein damage and post-ingestive aversion responses in caged bees, leading to lower tainted syrup consumption and higher water intake compared with the controls. No differentially expressed genes were observed among caged bees in regards to imidacloprid treatment. However, significant upregulation in antioxidant genes was recorded in caged bees as compared with hive bees, with overwhelming downregulation in all gene categories in caged bees at week 4. We identified two sets of genes that were constantly regulated in caged bees, including Rsod with unknown function in insects that could potentially characterize caging stress in honey bees.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

Targeting the potassium ion channel genes SK and SH as a novel approach for control of insect pests: efficacy and biosafety

BACKGROUND

Potassium ion channels play a critical role in the generation of electrical signals and thus provide potential targets for control of insect pests by RNA interference.

RESULTS

Genes encoding the small conductance calcium-activated potassium channel (SK) and the voltage-gated potassium channel (SH) were knocked down in Tribolium castaneum by injection and oral delivery of dsRNA (dsTcSK and dsTcSH, respectively). Irrespective of the delivery mechanism a dose-dependent effect was observed for knockdown (KD) of gene expression and insect mortality for both genes. Larvae fed a 400 ng dsRNA mg^-1 diet showed significant gene (P < 0.05) knockdown (98% and 83%) for SK and SH, respectively, with corresponding mortalities of 100% and 98% after 7 days. When injected (248.4 ng larva^-1 ), gene KD was 99% and 98% for SK and SH, causing 100% and 73.4% mortality, respectively. All developmental stages tested (larvae, early- and late-stage pupae and adults) showed an RNAi-sensitive response for both genes. LC50 values were lower for SK than SH, irrespective of delivery method, demonstrating that the knockdown of SK had a greater effect on larval mortality. Biosafety studies using adult honeybee Apis mellifera showed that there were no significant differences either in expression levels or mortality of honeybees orally dosed with dsTcSK and dsTcSH compared to control-fed bees. Similarly, there was no significant difference in the titre of deformed wing virus, used as a measure of immune suppression, between experimental and control bees.

CONCLUSION

This study demonstrates the potential of using RNAi targeting neural receptors as a technology for the control of T. castaneum. © 2019 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

The Emerging Proteomic Research Facilitates in-Depth Understanding of the Biology of Honeybees

Advances in instrumentation and computational analysis in proteomics have opened new doors for honeybee biological research at the molecular and biochemical levels. Proteomics has greatly expanded the understanding of honeybee biology since its introduction in 2005, through which key signaling pathways and proteins that drive honeybee development and behavioral physiology have been identified. This is critical for downstream mechanistic investigation by knocking a gene down/out or overexpressing it and being able to attribute a specific phenotype/biochemical change to that gene. Here, we review how emerging proteome research has contributed to the new understanding of honeybee biology. A systematic and comprehensive analysis of global scientific progress in honeybee proteome research is essential for a better understanding of research topics and trends, and is potentially useful for future research directions.

Characterization and its implication of a novel taste receptor detecting nutrients in the honey bee, Apis mellifera

Umami taste perception indicates the presence of amino acids, which are essential nutrients. Although the physiology of umami perception has been described in mammals, how insects detect amino acids remains unknown except in Drosophila melanogaster. We functionally characterized a gustatory receptor responding to L-amino acids in the western honey bee, Apis mellifera. Using a calcium-imaging assay and two-voltage clamp recording, we found that one of the honey bee's gustatory receptors, AmGr10, functions as a broadly tuned amino acid receptor responding to glutamate, aspartate, asparagine, arginine, lysine, and glutamine, but not to other sweet or bitter compounds. Furthermore, the sensitivity of AmGr10 to these L-amino acids was dramatically enhanced by purine ribonucleotides, like inosine-5'-monophosphate (IMP). Contact sensory hairs in the mouthpart of the honey bee responded strongly to glutamate and aspartate, which house gustatory receptor neurons expressing AmGr10. Interestingly, AmGr10 protein is highly conserved among hymenopterans but not other insects, implying unique functions in eusocial insects.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Metabolomics-based biomarker discovery for bee health monitoring: A proof of concept study concerning nutritional stress in Bombus terrestris

Bee pollinators are exposed to multiple natural and anthropogenic stressors. Understanding the effects of a single stressor in the complex environmental context of antagonistic/synergistic interactions is critical to pollinator monitoring and may serve as early warning system before a pollination crisis. This study aimed to methodically improve the diagnosis of bee stressors using a simultaneous untargeted and targeted metabolomics-based approach. Analysis of 84 Bombus terrestris hemolymph samples found 8 metabolites retained as potential biomarkers that showed excellent discrimination for nutritional stress. In parallel, 8 significantly altered metabolites, as revealed by targeted profiling, were also assigned as candidate biomarkers. Furthermore, machine learning algorithms were applied to the above-described two biomarker sets, whereby the untargeted eight components showed the best classification performance with sensitivity and specificity up to 99% and 100%, respectively. Based on pathway and biochemistry analysis, we propose that gluconeogenesis contributed significantly to blood sugar stability in bumblebees maintained on a low carbohydrate diet. Taken together, this study demonstrates that metabolomics-based biomarker discovery holds promising potential for improving bee health monitoring and to identify stressor related to energy intake and other environmental stressors.

Neuropeptides and peptide hormones identified in codling moth, Cydia pomonella (Lepidoptera: Tortricidae)

The codling moth, Cydia pomonella, is a worldwide pest of pome fruits. Neuropeptides regulate most physiological functions in insects and represent new targets for the development of control agents. The only neuropeptides reported from the codling moth to date are the allatostatin A family peptides. To identify other neuropeptides and peptide hormones from codling moth, we analyzed head transcriptomes, identified 50 transcripts, and predicted 120 prepropeptides for the codling moth neuropeptides and peptide hormones. All transcripts were amplified, and these sequences were verified. One of the notable findings in this study is that diapause hormones (DHs) reported from Tortricid moths, including the codling moth, do not have the WFGPRL sequence in C-terminal ends in the pban genes. The C-terminal motif is critical to characterize insect DH peptides, and always conserved in pban/dh genes in Lepidoptera and many insect orders. Interestingly, the WFGPRL sequence was produced only from the capa gene in the codling moth. The allatostatin A-family encoding transcript predicted nine peptides, seven of which, as expected, are identical to those previously isolated from the moth. We also identified new codling moth orthologs of insect neuropeptides including CCHamides, allatostatin CC, RYamides, and natalisins. The information provided in this study will benefit future codling moth investigations using peptidoproteomics to determine peptide presence and functions.

Rate variation in the evolution of non-coding DNA associated with social evolution in bees

The evolutionary origins of eusociality represent increases in complexity from individual to caste-based, group reproduction. These behavioural transitions have been hypothesized to go hand in hand with an increased ability to regulate when and where genes are expressed. Bees have convergently evolved eusociality up to five times, providing a framework to test this hypothesis. To examine potential links between putative gene regulatory elements and social evolution, we compare alignable, non-coding sequences in 11 diverse bee species, encompassing three independent origins of reproductive division of labour and two elaborations of eusocial complexity. We find that rates of evolution in a number of non-coding sequences correlate with key social transitions in bees. Interestingly, while we find little evidence for convergent rate changes associated with independent origins of social behaviour, a number of molecular pathways exhibit convergent rate changes in conjunction with subsequent elaborations of social organization. We also present evidence that many novel non-coding regions may have been recruited alongside the origin of sociality in corbiculate bees; these loci could represent gene regulatory elements associated with division of labour within this group. Thus, our findings are consistent with the hypothesis that gene regulatory innovations are associated with the evolution of eusociality and illustrate how a thorough examination of both coding and non-coding sequence can provide a more complete understanding of the molecular mechanisms underlying behavioural evolution. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Royal jelly-derived proteins enhance proliferation and migration of human epidermal keratinocytes in an in vitro scratch wound model

Background

Skin injury is inevitable in daily life. In recent years, with the increasing morbidity of diseases such as diabetes and metabolic disorders, chronic wounds have become a considerable challenge in clinical practice. Royal jelly, reported to have multifarious biological and physiological properties, has been used as a remedy for a variety of wounds since ancient times. However, the active components and mechanisms underlying the wound-healing properties of royal jelly are still largely unknown.

Methods

Water-soluble proteins of royal jelly were fractionated and investigated for the proliferative and migratory effects on human epidermal keratinocytes (HaCaT) in an in vitro wound healing model. The proteins present in bioactive fractions were characterised and quantified using Label-free protein quantification method. The potential functions of these proteins in biological systems were further analysed using bioinformatic tools.

Results

A protein fraction, mainly containing major royal jelly proteins 2 (MRJP2), MRJP3 and MRJP7, stimulated proliferative and migratory activities in HaCaT cells without visible cytotoxicity. It exerted the greatest effects on the growth of HaCaT cells in the first 48 h. Furthermore, when treated with this protein fraction, the closure rates of the in vitro scratch wound were significantly increased. Functional analysis indicated that MRJP2, MRJP3 and MRJP7 were associated with carbohydrate transport and metabolism.

Conclusions

We fractionated the water-soluble proteins of royal jelly and identified one fraction (Fraction 2) that induced both proliferative and migratory effects on a human epidermal keratinocyte cell line. Major royal jelly proteins (MRJP2, MRJP3 and/or MRJP7) were speculated to possess potential wound-healing bioactivity. This is the first report that royal jelly may improve wound closure via MRJP-induced cellular proliferation and migration. These proteins may be valuable lead compounds for the development of novel wound healing medications. Our findings would facilitate better understanding of the wound repair mechanisms of royal jelly.

Transcriptomic investigation of immune responses of the Apis cerana cerana larval gut infected by Ascosphaera apis

Chalkbrood is the most common fungal disease in honeybees. The objective of this study was to reveal immune responses in the Apis cerana cerana larval gut following Ascosphaera apis invasion. Combining a previously assembled transcriptome of A. c. cerana larval gut and the high-throughput sequencing data obtained in this study, 6152 differentially expressed genes (DEGs) were clustered into eight profiles. Trend analysis showed three significant up-regulated profiles (p ≤ 0.05) and three down-regulated profiles. Gene Ontology (GO) term analysis suggested that DEGs within significant up-regulated and down-regulated clusters were enriched in 46 and 38 functional groups, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis indicated a majority of DEGs were involved in ribosome structure or function, carbon metabolism, biosynthesis of amino acids, and oxidative phosphorylation. In addition, 142 and 14 DEGs were annotated in the cellular immune- and humoral immune-related pathways, respectively. Further investigation indicated that DEGs up-regulated in cellular immune and humoral immune pathways outnumbered those that were down-regulated. Moreover, immune responses of A. c. cerana and Apis mellifera ligustica larvae were compared and studied to decipher resistance of eastern honeybee larvae to A. apis. These results demonstrated that a large number of genes involved in immunity-related pathways were activated by A. apis. Our findings provided valuable information for elucidating the molecular mechanisms underlying immune responses of A. c. cerana larvae to A. apis infection and pathogen-host interactions during chalkbrood infection.

Long-Distance Transportation Causes Temperature Stress in the Honey Bee, Apis mellifera (Hymenoptera: Apidae)

Pollination services provided by the honey bee, Apis mellifera (Hymenoptera: Apidae, Linnaeus, 1758) have broad economic impacts and are necessary for production of a diversity of important crops. Hives may be transported multiple times per year to provide pollination. To test how temperature may contribute to transportation stress, temperature sensors were placed in hives in different locations and orientations on the trailer during shipping. Colony size prior to shipping significantly contributed to loss of population immediately after shipping which contributed to colony failure with smaller colonies more likely to fail and fail faster. Colony size also affects thermoregulation and temperature stress. Internal hive temperature varies significantly based on location and orientation. While colonies near the front and rear of the trailer and those oriented toward the center aisle had significantly different average internal temperatures, colony size best predicts loss of thermoregulation. Additionally, we profiled gene expression at departure, on arrival, and after a recovery period to identify transcriptional responses to transportation. Functional and enrichment analysis identified increased methylation and decreased ribosomal and protein-folding activity. Pheromone and odorant-binding transcripts were up-regulated after transportation. After recovery, transcripts associated with defense response, immune activity, and heat shock decreased, while production of antibiotic peptides increased. We conclude that hives experience considerable temperature stress possibly caused by turbulent airflow in exposed locations. Transportation stress should be considered an important component of annual colony losses which can be mitigated with improved management strategies.

Evolution of the Torso activation cassette, a pathway required for terminal patterning and moulting

Embryonic terminal patterning and moulting are critical developmental processes in insects. In Drosophila and Tribolium both of these processes are regulated by the Torso-activation cassette (TAC). The TAC consists of a common receptor, Torso, ligands Trunk and prothoracicotropic hormone (PTTH), and the spatially restricted protein Torso-like, with combinations of these elements acting mechanistically to activate the receptor in different developmental contexts. In order to trace the evolutionary history of the TAC we determined the presence or absence of TAC components in the genomes of arthropods. Our analyses reveal that Torso, Trunk and PTTH are evolutionarily labile components of the TAC with multiple individual or combined losses occurring in the arthropod lineages leading to and within the insects. These losses are often correlated, with both ligands and receptor missing from the genome of the same species. We determine that the PTTH gene evolved in the common ancestor of Hemiptera and Holometabola, and is missing from the genomes of a number of species with experimentally demonstrated PTTH activity, implying another molecule may be involved in ecdysis in these species. In contrast, the torso-like gene is a common component of pancrustacean genomes.

lncRNA profile of Apis mellifera and its possible role in behavioural transition from nurses to foragers

BACKGROUND

The behavioural transition from nurses to foragers in honey bees is known to be affected by intrinsic and extrinsic factors, including colony demography, hormone levels, brain chemistry and structure, and gene expression in the brain. However, the molecular mechanism underlying this behavioural transition of honey bees is still obscure.

RESULTS

Through RNA sequencing, we performed a comprehensive analysis of lncRNAs and mRNAs in honey bee nurses and foragers. Nurses and foragers from both typical colonies and single-cohort colonies were used to prepare six libraries to generate 49 to 100 million clear reads per sample. We obtained 6863 novel lncRNAs, 1480 differentially expressed lncRNAs between nurses and foragers, and 9308 mRNAs. Consistent with previous studies, lncRNAs showed features distinct from mRNAs, such as shorter lengths, lower exon numbers, and lower expression levels compared to mRNAs. Bioinformatic analysis showed that differentially expressed genes were mostly involved in the regulation of sensory-related events, such as olfactory receptor activity and odorant binding, and enriched Wnt and FoxO signaling pathways. Moreover, we found that lncRNAs TCONS_00356023, TCONS_00357367, TCONS_00159909 and mRNAs dop1, Kr-h1 and HR38 may play important roles in behavioural transition in honey bees.

CONCLUSION

This study characterized the expression profile of lncRNAs in nurses and foragers and provided a framework for further study of the role of lncRNAs in honey bee behavioural transition.

A gene for resistance to the Varroa mite (Acari) in honey bee (Apis mellifera) pupae

Social insect colonies possess a range of defences which protect them against highly virulent parasites and colony collapse. The host-parasite interaction between honey bees (Apis mellifera) and the mite Varroa destructor is unusual, as honey bee colonies are relatively poorly defended against this parasite. The interaction has existed since the mid-20th Century, when Varroa switched host to parasitize A. mellifera. The combination of a virulent parasite and relatively naïve host means that, without acaricides, honey bee colonies typically die within 3 years of Varroa infestation. A consequence of acaricide use has been a reduced selective pressure for the evolution of Varroa resistance in honey bee colonies. However, in the past 20 years, several natural-selection-based breeding programmes have resulted in the evolution of Varroa-resistant populations. In these populations, the inhibition of Varroa's reproduction is a common trait. Using a high-density genome-wide association analysis in a Varroa-resistant honey bee population, we identify an ecdysone-induced gene significantly linked to resistance. Ecdysone both initiates metamorphosis in insects and reproduction in Varroa. Previously, using a less dense genetic map and a quantitative trait loci analysis, we have identified Ecdysone-related genes at resistance loci in an independently evolved resistant population. Varroa cannot biosynthesize ecdysone but can acquire it from its diet. Using qPCR, we are able to link the expression of ecdysone-linked resistance genes to Varroa's meals and reproduction. If Varroa co-opts pupal compounds to initiate and time its own reproduction, mutations in the host's ecdysone pathway may represent a key selection tool for honey bee resistance and breeding.

Telomerase Reverse Transcriptase Gene from Greasyback Shrimp Metapenaeus ensis: Isolation, Genomic Organization and Expression Pattern Analysis

Telomerase reverse transcriptase (TERT), a catalytic subunit of telomerase, plays a key role in the activity and biological functions of telomerase. In the present study we isolated and characterized the full-length cDNA and DNA sequences of the TERT gene (MeTERT) from Metapenaeus ensis. MeTERT cDNA was 4239 bp in length, which consisted of a 369 bp 5'UTR, a 3231 bp open reading frame encoding 1076 amino acids, and a 639 bp 3'UTR. The genomic DNA of MeTERT had only two introns, similar to beetle (two introns) and silkworm (intronless). The MeTERT protein showed only 5.2-7.9% identity with other known TERTs but contained all the four primary TERT domains of the N-terminal TEN, RNA binding domain (TRBD), reverse transcriptase (RT) and C-terminus CTE. Expression pattern analysis by RT-qPCR showed that, the MeTERT mRNA transcripts could be detected in all the tested samples, with relatively higher expression level in the gill, mysis, Oka organ and egg, but lower level in muscle, ovary, in vitro cultured 3-d Oka organ cells and heart. The significant decrease of MeTERT expression in the in vitro cultured 3-d Oka organ primary cells compared with their source tissue of Oka organ may have contributed to the cellular mitosisarrest. Thus trans-activation of TERT gene may be a candidate in attempts to immortalize in vitro cultured shrimp cells. This work will lay a solid foundation for future studies of the biological functions of telomerase in crustaceans.

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

Background

Adapting to changes in the environment is the foundation of species survival, and is usually thought to be a gradual process. However, transposable elements (TEs), epigenetic modifications, and/or genetic material acquired from other organisms by means of horizontal gene transfer (HGTs), can also lead to novel adaptive traits. Social insects form dense societies, which attract and maintain extra- and intracellular accessory inhabitants, which may facilitate gene transfer between species. The wood ant Formica exsecta (Formicidae; Hymenoptera), is a common ant species throughout the Palearctic region. The species is a well-established model for studies of ecological characteristics and evolutionary conflict.

Results

In this study, we sequenced and assembled draft genomes for F. exsecta and its endosymbiont Wolbachia. The F. exsecta draft genome is 277.7 Mb long; we identify 13,767 protein coding genes, for which we provide gene ontology and protein domain annotations. This is also the first report of a Wolbachia genome from ants, and provides insights into the phylogenetic position of this endosymbiont. We also identified multiple horizontal gene transfer events (HGTs) from Wolbachia to F. exsecta. Some of these HGTs have also occurred in parallel in multiple other insect genomes, highlighting the extent of HGTs in eukaryotes.

Conclusion

We present the first draft genome of ant F. exsecta, and its endosymbiont Wolbachia (wFex), and show considerable rates of gene transfer from the symbiont to the host. We expect that especially the F. exsecta genome will be valuable resource in further exploration of the molecular basis of the evolution of social organization.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5665-6) contains supplementary material, which is available to authorized users.

Comparative Analysis of Brain and Fat Body Gene Splicing Patterns in the Honey Bee, Apis mellifera

RNA-seq has proven to be a powerful tool to unravel various aspects of the transcriptome, especially the quantification of alternative splicing (AS) that leads to isoform diversity. The honey bee (Apis mellifera) is an important model organism for studying the molecular underpinnings of behavioral plasticity and social behavior, and recent RNA-seq studies of honey bees have revealed AS patterns and their regulation by DNA methylation. However, tissue-specific AS patterns have not been fully explored. In this paper, we characterized AS patterns in two different honey bee tissue types, and also explored their conservation and regulation. We used the RNA-seq data from brain and fat body to improve the existing models of honey bee genes and identified tissue-specific AS patterns. We found that AS genes show high conservation between honey bee and Drosophila melanogaster. We also confirmed and extended previous findings of a correlation between gene body DNA methylation and AS patterns, providing further support for the role of DNA methylation in regulating AS. In addition, our analysis suggests distinct functional roles for tissue-specific alternatively spliced genes. Taken together, our work provides new insights into the conservation and dynamics of AS patterns across different tissue types.

Comparative transcriptomics of social insect queen pheromones

Queen pheromones are chemical signals that mediate reproductive division of labor in eusocial animals. Remarkably, queen pheromones are composed of identical or chemically similar compounds in some ants, wasps and bees, even though these taxa diverged >150MYA and evolved queens and workers independently. Here, we measure the transcriptomic consequences of experimental exposure to queen pheromones in workers from two ant and two bee species (genera: Lasius, Apis, Bombus), and test whether they are similar across species. Queen pheromone exposure affected transcription and splicing at many loci. Many genes responded consistently in multiple species, and the set of pheromone-sensitive genes was enriched for functions relating to lipid biosynthesis and transport, olfaction, production of cuticle, oogenesis, and histone (de)acetylation. Pheromone-sensitive genes tend to be evolutionarily ancient, positively selected, peripheral in the gene coexpression network, hypomethylated, and caste-specific in their expression. Our results reveal how queen pheromones achieve their effects, and suggest that ants and bees use similar genetic modules to achieve reproductive division of labor.

A hybrid de novo genome assembly of the honeybee, Apis mellifera, with chromosome-length scaffolds

Background

The ability to generate long sequencing reads and access long-range linkage information is revolutionizing the quality and completeness of genome assemblies. Here we use a hybrid approach that combines data from four genome sequencing and mapping technologies to generate a new genome assembly of the honeybee Apis mellifera. We first generated contigs based on PacBio sequencing libraries, which were then merged with linked-read 10x Chromium data followed by scaffolding using a BioNano optical genome map and a Hi-C chromatin interaction map, complemented by a genetic linkage map.

Results

Each of the assembly steps reduced the number of gaps and incorporated a substantial amount of additional sequence into scaffolds. The new assembly (Amel_HAv3) is significantly more contiguous and complete than the previous one (Amel_4.5), based mainly on Sanger sequencing reads. N50 of contigs is 120-fold higher (5.381 Mbp compared to 0.053 Mbp) and we anchor > 98% of the sequence to chromosomes. All of the 16 chromosomes are represented as single scaffolds with an average of three sequence gaps per chromosome. The improvements are largely due to the inclusion of repetitive sequence that was unplaced in previous assemblies. In particular, our assembly is highly contiguous across centromeres and telomeres and includes hundreds of AvaI and AluI repeats associated with these features.

Conclusions

The improved assembly will be of utility for refining gene models, studying genome function, mapping functional genetic variation, identification of structural variants, and comparative genomics.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5642-0) contains supplementary material, which is available to authorized users.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Characterization of locomotor response to psychostimulants in the parthenogenetic marbled crayfish (Procambarus fallax forma virginalis): A promising model for studying the neural and molecular mechanisms of drug addiction

Although scientific research using mammalian models has made great strides in uncovering the enigmatic neural and molecular mechanisms orchestrating the state of drug addiction, a complete understanding has thus far eluded researchers. The complexity of the task has led to the use of invertebrate model systems to complement the research of drug-induced reward in mammalian systems. Invertebrates, such as crayfish, offer excellent model systems to help reveal the underlying mechanisms of drug addiction as they retain the ancestral neural reward circuit that is evolutionarily conserved across taxa, and they possess relatively few, large neurons, laid out in an accessible, modularly organized nervous system. Crayfish offer the benefits of delineated developmental life stages, a large body size suitable for a variety of experimental methods, and stereotyped behaviors. Unique among crayfish is the parthenogenetic marbled crayfish (Procambarus fallax forma virginalis), a species of asexually reproducing, genetically identical clones. With the benefits of reduced individual variation, high fecundity, and easy lab husbandry, the marbled crayfish would make a particularly powerful addition to the animal model repertoire. Here we characterize the locomotor response of juvenile P. f. f. virginalis exposed to the psychostimulant, d-amphetamine sulfate. Custom video-tracking software was used to record the movement patterns of juveniles exposed to water infused with varying concentrations of d-amphetamine sulfate. ANOVA demonstrated that crayfish locomotion was significantly impacted by drug concentration. These psychostimulant effects provide the foundation of P. f. f. virginalis as a model for parsing the neural and molecular mechanisms of drug addiction.

Revisiting Polybia paulista wasp venom using shotgun proteomics - Insights into the N-linked glycosylated venom proteins

The partial proteome of Polybia paulista wasp venom was previously reported elsewhere using a gel-dependent approach and resulted in the identification of a limited number of venom toxins. Here, we reinvestigated the P. paulista venom using a gel-free shotgun proteomic approach; the highly dynamic range of this approach facilitated the detection and identification of 1673 proteins, of which 23 venom proteins presented N-linked glycosylation as a posttranslational modification. Three different molecular forms of PLA1 were identified as allergenic proteins, and two of these forms were modified by N-linked glycosylation. This study reveals an extensive repertoire of hitherto undescribed proteins that were classified into the following six different functional groups: (i) typical venom proteins; (ii) proteins related to the folding/conformation and PTMs of toxins; (iii) proteins that protect toxins from oxidative stress; (iv) proteins involved in chemical communication; (v) housekeeping proteins; and (vi) uncharacterized proteins. It was possible to identify venom toxin-like proteins that are commonly reported in other animal venoms, including arthropods such as spiders and scorpions. Thus, the findings reported here may contribute to improving our understanding of the composition of P. paulista venom, its envenoming mechanism and the pathologies experienced by the victim after the wasp stinging accident. BIOLOGICAL SIGNIFICANCE: The present study significantly expanded the number of proteins identified in P. paulista venom, contributing to improvements in our understanding of the envenoming mechanism produced by sting accidents caused by this wasp. For example, novel wasp venom neurotoxins have been identified, but no studies have assessed the presence of this type of toxin in social wasp venoms. In addition, 23 N-linked glycosylated venom proteins were identified in the P. paulista venom proteome, and some of these proteins might be relevant allergens that are immunoreactive to human IgE.

Exploring integument transcriptomes, cuticle ultrastructure, and cuticular hydrocarbons profiles in eusocial and solitary bee species displaying heterochronic adult cuticle maturation

Differences in the timing of exoskeleton melanization and sclerotization are evident when comparing eusocial and solitary bees. This cuticular maturation heterochrony may be associated with life style, considering that eusocial bees remain protected inside the nest for many days after emergence, while the solitary bees immediately start outside activities. To address this issue, we characterized gene expression using large-scale RNA sequencing (RNA-seq), and quantified cuticular hydrocarbon (CHC) through gas chromatography-mass spectrometry in comparative studies of the integument (cuticle plus its underlying epidermis) of two eusocial and a solitary bee species. In addition, we used transmission electron microscopy (TEM) for studying the developing cuticle of these and other three bee species also differing in life style. We found 13,200, 55,209 and 30,161 transcript types in the integument of the eusocial Apis mellifera and Frieseomelitta varia, and the solitary Centris analis, respectively. In general, structural cuticle proteins and chitin-related genes were upregulated in pharate-adults and newly-emerged bees whereas transcripts for odorant binding proteins, cytochrome P450 and antioxidant proteins were overrepresented in foragers. Consistent with our hypothesis, a distance correlation analysis based on the differentially expressed genes suggested delayed cuticle maturation in A. mellifera in comparison to the solitary bee. However, this was not confirmed in the comparison with F. varia. The expression profiles of 27 of 119 genes displaying functional attributes related to cuticle formation/differentiation were positively correlated between A. mellifera and F. varia, and negatively or non-correlated with C. analis, suggesting roles in cuticular maturation heterochrony. However, we also found transcript profiles positively correlated between each one of the eusocial species and C. analis. Gene co-expression networks greatly differed between the bee species, but we identified common gene interactions exclusively between the eusocial species. Except for F. varia, the TEM analysis is consistent with cuticle development timing adapted to the social or solitary life style. In support to our hypothesis, the absolute quantities of n-alkanes and unsaturated CHCs were significantly higher in foragers than in the earlier developmental phases of the eusocial bees, but did not discriminate newly-emerged from foragers in C. analis. By highlighting differences in integument gene expression, cuticle ultrastructure, and CHC profiles between eusocial and solitary bees, our data provided insights into the process of heterochronic cuticle maturation associated to the way of life.

Pigment-Dispersing Factor-expressing neurons convey circadian information in the honey bee brain

Pigment-Dispersing Factor (PDF) is an important neuropeptide in the brain circadian network of Drosophila and other insects, but its role in bees in which the circadian clock influences complex behaviour is not well understood. We combined high-resolution neuroanatomical characterizations, quantification of PDF levels over the day and brain injections of synthetic PDF peptide to study the role of PDF in the honey bee Apis mellifera We show that PDF co-localizes with the clock protein Period (PER) in a cluster of laterally located neurons and that the widespread arborizations of these PER/PDF neurons are in close vicinity to other PER-positive cells (neurons and glia). PDF-immunostaining intensity oscillates in a diurnal and circadian manner with possible influences for age or worker task on synchrony of oscillations in different brain areas. Finally, PDF injection into the area between optic lobes and the central brain at the end of the subjective day produced a consistent trend of phase-delayed circadian rhythms in locomotor activity. Altogether, these results are consistent with the hypothesis that PDF is a neuromodulator that conveys circadian information from pacemaker cells to brain centres involved in diverse functions including locomotion, time memory and sun-compass orientation.

Ras-Specific GTPase-Activating Proteins-Structures, Mechanisms, and Interactions

Ras-specific GTPase-activating proteins (RasGAPs) down-regulate the biological activity of Ras proteins by accelerating their intrinsic rate of GTP hydrolysis, basically by a transition state stabilizing mechanism. Oncogenic Ras is commonly not sensitive to RasGAPs caused by interference of mutants with the electronic or steric requirements of the transition state, resulting in up-regulation of activated Ras in respective cells. RasGAPs are modular proteins containing a helical catalytic RasGAP module surrounded by smaller domains that are frequently involved in the subcellular localization or contributing to regulatory features of their host proteins. In this review, we summarize current knowledge about RasGAP structure, mechanism, regulation, and dual-substrate specificity and discuss in some detail neurofibromin, one of the most important negative Ras regulators in cellular growth control and neuronal function.

Laboratory bioassays on the response of honey bee (Apis mellifera L.) glutathione S-transferase and acetylcholinesterase to the oral exposure to copper, cadmium, and lead

In the present study, the influence of cadmium, copper, and lead on two enzymes often used as biomarkers in toxicological analysis was investigated. Bees were fed with 1 M sucrose solution containing 10-fold serial dilutions of CuCl₂ (1000 mg L^-1, 100 mg L^-1, and 10 mg L^-1), CdCl₂ (0.1 mg L^-1, 0.01 mg L^-1, and 0.001 mg L^-1), or PbCl₂ (10 mg L^-1, 1 mg L^-1, and 0.1 mg L^-1) during 48 h. Our results showed that the total glutathione S-transferase activity was not changed under the influence of cadmium and lead, and it was decreased with the highest concentration of copper. The level of gene expression of the three analyzed classes of glutathione S-transferase was significantly increased with increasing concentrations of copper and cadmium. Lead did not cause significant changes in glutathione S-transferase activity and gene expression, while it showed biphasic effect on acetylcholinesterase activity: lower concentration of lead, 0.1 mg L^-1 inhibited and higher dose, 10 mg L^-1 induced acetylcholinesterase activity in honey bees. Furthermore, our results showed a significant decrease of the acetylcholinesterase activity in honey bees treated with 0.001 and 0.01 mg L^-1 CdCl₂. Our results indicate the influence of cadmium, copper, and lead on GST and AChE in the honey bees. These results form the basis for future research on the impact of metallic trace element pollution on honey bees.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.

Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

Recent advances in population and quantitative genomics of honey bees

The increase in the availability of individual Apis mellifera genomes has resulted in significant progress toward understanding the evolution and adaptation of the honey bee. These efforts have identified new subspecies, evolutionary lineages, and a significant number of genes involved with adaptations and colony-level quantitative traits. Many studies have also developed genetic assays that are being used to monitor the movement and admixture of honey bee populations. These resources are valuable for conservation and breeding programs that seek to improve the economic value of colonies or preserve locally adapted populations and subspecies. This review provides a brief discussion on how population and quantitative genomic studies has improved our understanding of the honey bee.

Hormonal control and target genes of ftz-f1 expression in the honeybee Apis mellifera: a positive loop linking juvenile hormone, ftz-f1, and vitellogenin

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg.

Immediate early genes in social insects: a tool to identify brain regions involved in complex behaviors and molecular processes underlying neuroplasticity

Social insects show complex behaviors and master cognitive tasks. The underlying neuronal mechanisms, however, are in most cases only poorly understood due to challenges in monitoring brain activity in freely moving animals. Immediate early genes (IEGs) that get rapidly and transiently expressed following neuronal stimulation provide a powerful tool for detecting behavior-related neuronal activity in vertebrates. In social insects, like honey bees, and in insects in general, this approach is not yet routinely established, even though these genes are highly conserved. First studies revealed a vast potential of using IEGs as neuronal activity markers to analyze the localization, function, and plasticity of neuronal circuits underlying complex social behaviors. We summarize the current knowledge on IEGs in social insects and provide ideas for future research directions.

Draft genome of the cotton aphid Aphis gossypii

The cotton aphid Aphis gossypii Glover is a worldwide agricultural pest that feeds on cotton, melon, and other landscape plants, causing a high level of economic loss. In addition to the common characteristics shared with other aphids, Ap. gossypii has evolved multiple biotypes that present substantial differences in host adaption. These intriguing biological features are of interest from both a fundamental and applied perspective. However, the molecular studies of Ap. gossypii have been restrained by the lack of a reference genome. Furthermore, in order to establish a platform for the development of novel and sustainable control methods, it is necessary to generate genomic resources for Ap. gossypii. Here, we present a 294 Mb draft genome sequence of Ap. gossypii, which consists of 4,724 scaffolds with an N50 size of 438 kb. Compared to other aphid species with published genomes, Ap. gossypii presents the most compact genome size. A total of 14,694 protein-coding genes were predicted and annotated in the consensus gene set, 98.03% of CEGMA genes and 93.5% of BUSCO genes were captured respectively. Genome-wide selection analyses revealed that significantly evolving pathways in the genus Aphis are related to biological processes of detoxification, steroid biosynthesis, and ethylbenzene degradation. The acquisition of the genome of Ap. gossypii makes it possible to understand the molecular mechanism of intricate biological traits of this species, and will further facilitate the study of aphid evolution.

Diversity and evolution of the transposable element repertoire in arthropods with particular reference to insects

BACKGROUND

Transposable elements (TEs) are a major component of metazoan genomes and are associated with a variety of mechanisms that shape genome architecture and evolution. Despite the ever-growing number of insect genomes sequenced to date, our understanding of the diversity and evolution of insect TEs remains poor.

RESULTS

Here, we present a standardized characterization and an order-level comparison of arthropod TE repertoires, encompassing 62 insect and 11 outgroup species. The insect TE repertoire contains TEs of almost every class previously described, and in some cases even TEs previously reported only from vertebrates and plants. Additionally, we identified a large fraction of unclassifiable TEs. We found high variation in TE content, ranging from less than 6% in the antarctic midge (Diptera), the honey bee and the turnip sawfly (Hymenoptera) to more than 58% in the malaria mosquito (Diptera) and the migratory locust (Orthoptera), and a possible relationship between the content and diversity of TEs and the genome size.

CONCLUSION

While most insect orders exhibit a characteristic TE composition, we also observed intraordinal differences, e.g., in Diptera, Hymenoptera, and Hemiptera. Our findings shed light on common patterns and reveal lineage-specific differences in content and evolution of TEs in insects. We anticipate our study to provide the basis for future comparative research on the insect TE repertoire.

Parasitoid Jewel Wasp Mounts Multipronged Neurochemical Attack to Hijack a Host Brain

The parasitoid emerald jewel wasp Ampulex compressa induces a compliant state of hypokinesia in its host, the American cockroach Periplaneta americana through direct envenomation of the central nervous system (CNS). To elucidate the biochemical strategy underlying venom-induced hypokinesia, we subjected the venom apparatus and milked venom to RNAseq and proteomics analyses to construct a comprehensive "venome," consisting of 264 proteins. Abundant in the venome are enzymes endogenous to the host brain, including M13 family metalloproteases, phospholipases, adenosine deaminase, hyaluronidase, and neuropeptide precursors. The amphipathic, alpha-helical ampulexins are among the most abundant venom components. Also prominent are members of the Toll/NF-κB signaling pathway, including proteases Persephone, Snake, Easter, and the Toll receptor ligand Spätzle. We find evidence that venom components are processed following envenomation. The acidic (pH∼4) venom contains unprocessed neuropeptide tachykinin and corazonin precursors and is conspicuously devoid of the corresponding processed, biologically active peptides. Neutralization of venom leads to appearance of mature tachykinin and corazonin, suggesting that the wasp employs precursors as a prolonged time-release strategy within the host brain post-envenomation. Injection of fully processed tachykinin into host cephalic ganglia elicits short-term hypokinesia. Ion channel modifiers and cytolytic toxins are absent in A. compressa venom, which appears to hijack control of the host brain by introducing a "storm" of its own neurochemicals. Our findings deepen understanding of the chemical warfare underlying host-parasitoid interactions and in particular neuromodulatory mechanisms that enable manipulation of host behavior to suit the nutritional needs of opportunistic parasitoid progeny.

The Phosphofructokinase and Pyruvate Kinase Genes in Apis andreniformis and Apis cerana indica: Exon Intron Organisation and Evolution

Genes related to carbohydrate metabolism have evolved rapidly in eusocial bees, including honey bees. However, the characterisation of carbohydrate metabolism genes has not been reported in Apis andreniformis or Apis cerana indica. This study aimed to characterise phosphofructokinase (PFK) and pyruvate kinase (PK) genes in these honey bee species and to analyse the evolution of the genus Apis using these genes. This study found the first data regarding A. andreniformis PFK and PK-like nucleotide sequences. A BLAST-n algorithm-based search showed that A. andreniformis and A. c. indica PFK and PK genes were homologous with those of Apis florea and Apis cerana cerana from Korea, respectively. Multiple alignments of PFKs from five Apis species showed many exon gains and losses, but only one among the PKs. Thus, the exon-intron organisation of the PK genes may be more conserved compare with that of the PFKs. Another evolutionary pattern indicated that more nucleotide substitutions occurred in Apis' PK than PFK genes. Deduced PFK amino acid sequences revealed a PFK consensus pattern of 19 amino acids, while the deduced PK amino acid sequences were predicted to have barrel and alpha/beta domains. Based on these two metabolism-related genes, The Neighbour-joining and Maximum likelihood phylogenetic trees are congruent and revealed that the A. andreniformis and A. florea group were in the basal position. Apis mellifera, A. cerana, and Apis dorsata formed a monophyletic clade, although the positions of A. mellifera and A. dorsata were different in the nucleotide- and amino acid-based phylogenetic trees.

Genomic Signals of Adaptation towards Mutualism and Sociality in Two Ambrosia Beetle Complexes

Mutualistic symbiosis and eusociality have developed through gradual evolutionary processes at different times in specific lineages. Like some species of termites and ants, ambrosia beetles have independently evolved a mutualistic nutritional symbiosis with fungi, which has been associated with the evolution of complex social behaviors in some members of this group. We sequenced the transcriptomes of two ambrosia complexes (Euwallacea sp. near fornicatus⁻Fusarium euwallaceae and Xyleborus glabratus⁻Raffaelea lauricola) to find evolutionary signatures associated with mutualism and behavior evolution. We identified signatures of positive selection in genes related to nutrient homeostasis; regulation of gene expression; development and function of the nervous system, which may be involved in diet specialization; behavioral changes; and social evolution in this lineage. Finally, we found convergent changes in evolutionary rates of proteins across lineages with phylogenetically independent origins of sociality and mutualism, suggesting a constrained evolution of conserved genes in social species, and an evolutionary rate acceleration related to changes in selective pressures in mutualistic lineages.

Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism

Background

Trombidid mites have a unique life cycle in which only the larval stage is ectoparasitic. In the superfamily Trombiculoidea ("chiggers"), the larvae feed preferentially on vertebrates, including humans. Species in the genus Leptotrombidium are vectors of a potentially fatal bacterial infection, scrub typhus, that affects 1 million people annually. Moreover, chiggers can cause pruritic dermatitis (trombiculiasis) in humans and domesticated animals. In the Trombidioidea (velvet mites), the larvae feed on other arthropods and are potential biological control agents for agricultural pests. Here, we present the first trombidid mites genomes, obtained both for a chigger, Leptotrombidium deliense, and for a velvet mite, Dinothrombium tinctorium.

Results

Sequencing was performed using Illumina technology. A 180 Mb draft assembly for D. tinctorium was generated from two paired-end and one mate-pair library using a single adult specimen. For L. deliense, a lower-coverage draft assembly (117 Mb) was obtained using pooled, engorged larvae with a single paired-end library. Remarkably, both genomes exhibited evidence of ancient lateral gene transfer from soil-derived bacteria or fungi. The transferred genes confer functions that are rare in animals, including terpene and carotenoid synthesis. Thirty-seven allergenic protein families were predicted in the L. deliense genome, of which nine were unique. Preliminary proteomic analyses identified several of these putative allergens in larvae.

Conclusions

Trombidid mite genomes appear to be more dynamic than those of other acariform mites. A priority for future research is to determine the biological function of terpene synthesis in this taxon and its potential for exploitation in disease control.

Genome of the small hive beetle (Aethina tumida, Coleoptera: Nitidulidae), a worldwide parasite of social bee colonies, provides insights into detoxification and herbivory

Background

The small hive beetle (Aethina tumida; ATUMI) is an invasive parasite of bee colonies. ATUMI feeds on both fruits and bee nest products, facilitating its spread and increasing its impact on honey bees and other pollinators. We have sequenced and annotated the ATUMI genome, providing the first genomic resources for this species and for the Nitidulidae, a beetle family that is closely related to the extraordinarily species-rich clade of beetles known as the Phytophaga. ATUMI thus provides a contrasting view as a neighbor for one of the most successful known animal groups.

Results

We present a robust genome assembly and a gene set possessing 97.5% of the core proteins known from the holometabolous insects. The ATUMI genome encodes fewer enzymes for plant digestion than the genomes of wood-feeding beetles but nonetheless shows signs of broad metabolic plasticity. Gustatory receptors are few in number compared to other beetles, especially receptors with known sensitivity (in other beetles) to bitter substances. In contrast, several gene families implicated in detoxification of insecticides and adaptation to diverse dietary resources show increased copy numbers. The presence and diversity of homologs involved in detoxification differ substantially from the bee hosts of ATUMI.

Conclusions

Our results provide new insights into the genomic basis for local adaption and invasiveness in ATUMI and a blueprint for control strategies that target this pest without harming their honey bee hosts. A minimal set of gustatory receptors is consistent with the observation that, once a host colony is invaded, food resources are predictable. Unique detoxification pathways and pathway members can help identify which treatments might control this species even in the presence of honey bees, which are notoriously sensitive to pesticides.

Reconstructing the functions of endosymbiotic Mollicutes in fungus-growing ants

Mollicutes, a widespread class of bacteria associated with animals and plants, were recently identified as abundant abdominal endosymbionts in healthy workers of attine fungus-farming leaf-cutting ants. We obtained draft genomes of the two most common strains harbored by Panamanian fungus-growing ants. Reconstructions of their functional significance showed that they are independently acquired symbionts, most likely to decompose excess arginine consistent with the farmed fungal cultivars providing this nitrogen-rich amino-acid in variable quantities. Across the attine lineages, the relative abundances of the two Mollicutes strains are associated with the substrate types that foraging workers offer to fungus gardens. One of the symbionts is specific to the leaf-cutting ants and has special genomic machinery to catabolize citrate/glucose into acetate, which appears to deliver direct metabolic energy to the ant workers. Unlike other Mollicutes associated with insect hosts, both attine ant strains have complete phage-defense systems, underlining that they are actively maintained as mutualistic symbionts.