Finding the missing honey bee genes: lessons learned from a genome upgrade

Reference genome for the Mojave poppy bee (Perdita meconis), a specialist pollinator of conservation concern

The Mojave poppy bee, Perdita meconis Griswold (Hymenoptera: Anthophila: Andrenidae), is a species of conservation concern that is restricted to the eastern Mojave Desert of North America. It is a specialist pollinator of two poppy genera, Arctomecon and Argemone (Papaveraceae), and is being considered for listing under the US Endangered Species Act along with one of its pollinator hosts, the Las Vegas bearpoppy (Arctomecon californica). Here, we present a near chromosome-level genome of the Mojave poppy bee to provide a genomic resource that will aid conservation efforts and future research. We isolated DNA from a single, small (<7 mm), male specimen collected using non-ideal preservation methods and then performed whole-genome sequencing using PacBio HiFi technology. After quality and contaminant filtering, the final draft genome assembly is 327 Mb, with an N50 length of 17.5 Mb. Annotated repetitive elements compose 37.3% of the genome, although a large proportion (24.87%) of those are unclassified repeats. Additionally, we annotated 18,245 protein-coding genes and 19,433 transcripts. This genome represents one of only a few genomes from the large bee family Andrenidae and one of only a few genomes for pollinator specialists. We highlight both the potential of this genome as a resource for future research, and how high-quality genomes generated from small, non-ideal (in terms of preservation) specimens could facilitate biodiversity genomics.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

BACKGROUND

The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses.

RESULTS

Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants.

CONCLUSIONS

Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Early life exposure to queen mandibular pheromone mediates persistent transcriptional changes in the brain of honey bee foragers

Behavioural regulation in insect societies remains a fundamental question in sociobiology. In hymenopteran societies, the queen plays a crucial role in regulating group behaviour by affecting individual behaviour and physiology through modulation of worker gene expression. Honey bee (Apis mellifera) queens signal their presence via queen mandibular pheromone (QMP). While QMP has been shown to influence behaviour and gene expression of young workers, we know little about how these changes translate in older workers. The effects of the queen pheromone could have prolonged molecular impacts on workers that depend on an early sensitive period. We demonstrate that removal of QMP impacts long-term gene expression in the brain and antennae in foragers that were treated early in life (1 day post emergence), but not when treated later in life. Genes important for division of labour, learning, chemosensory perception and ageing were among those differentially expressed in the antennae and brain tissues, suggesting that QMP influences diverse physiological and behavioural processes in workers. Surprisingly, removal of QMP did not have an impact on foraging behaviour. Overall, our study suggests a sensitive period early in the life of workers, where the presence or absence of a queen has potentially life-long effects on transcriptional activity.

Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.

Effects of Arboviral Infections on Transposable Element Transcript Levels in Aedes aegypti

Transposable elements are mobile repeated sequences found in all genomes. Transposable elements are controlled by RNA interference pathways in most organisms, and this control involves the PIWI-interacting RNA pathway and the small interfering RNA pathway, which is also known to be the first line of antiviral defense in invertebrates. Using Drosophila, we recently showed that viral infections result in the modulation of transposable element transcript levels through modulation of the small RNA repertoire. The Aedes aegypti mosquito is of particular interest because almost half of its genome is made of transposable elements, and it is described as a major vector of viruses (such as the dengue [DENV], Zika [ZIKV], and chikungunya [CHIKV] arboviruses). Moreover, Aedes mosquitoes are unique among insects in that the PIWI-interacting RNA pathway is also involved in the somatic antiviral response, in addition to the transposable element control and PIWI-interacting RNA pathway genes expanded in the mosquito genome. For these reasons, we studied the impacts of viral infections on transposable element transcript levels in A. aegypti samples. We retrieved public datasets corresponding to RNA-seq data obtained from viral infections by DENV, ZIKV, and CHIKV in various tissues. We found that transposable element transcripts are moderately modulated following viral infection and that the direction of the modulation varies greatly across tissues and viruses. These results highlight the need for an in-depth investigation of the tightly intertwined interactions between transposable elements and viruses.

Transcriptional dynamics and regulatory function of milRNAs in Ascosphaera apis invading Apis mellifera larvae

In the present study, small RNA (sRNA) data from Ascosphaera apis were filtered from sRNA-seq datasets from the gut tissues of A. apis-infected Apis mellifera ligustica worker larvae, which were combined with the previously gained sRNA-seq data from A. apis spores to screen differentially expressed milRNAs (DEmilRNAs), followed by trend analysis and investigation of the DEmilRNAs in relation to significant trends. Additionally, the interactions between the DEmilRNAs and their target mRNAs were verified using a dual-luciferase reporter assay. In total, 974 A. apis milRNAs were identified. The first base of these milRNAs was biased toward U. The expression of six milRNAs was confirmed by stem-loop RT-PCR, and the sequences of milR-3245-y and milR-10285-y were validated using Sanger sequencing. These miRNAs grouped into four significant trends, with the target mRNAs of DEmilRNAs involving 42 GO terms and 120 KEGG pathways, such as the fungal-type cell wall and biosynthesis of secondary metabolites. Further investigation demonstrated that 299 DEmilRNAs (novel-m0011-3p, milR-10048-y, bantam-y, etc.) potentially targeted nine genes encoding secondary metabolite-associated enzymes, while 258 (milR-25-y, milR-14-y, milR-932-x, etc.) and 419 (milR-4561-y, milR-10125-y, let-7-x, etc.) DEmilRNAs putatively targeted virulence factor-encoded genes and nine genes involved in the MAPK signaling pathway, respectively. Additionally, the interaction between ADM-B and milR-6882-x, as well as between PKIA and milR-7009-x were verified. Together, these results not only offer a basis for clarifying the mechanisms underlying DEmilRNA-regulated pathogenesis of A. apis and a novel insight into the interaction between A. apis and honey bee larvae, but also provide candidate DEmilRNA-gene axis for further investigation.

GC Content Across Insect Genomes: Phylogenetic Patterns, Causes and Consequences

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Diptera, Coleoptera and Lepidoptera. Two species of Diptera, Bombylius major and B. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Diptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Helenus and Ajax, Two Groups of Non-Autonomous LTR Retrotransposons, Represent a New Type of Small RNA Gene-Derived Mobile Elements

Terminal repeat retrotransposons in miniature (TRIMs) are short non-autonomous long terminal repeat (LTR) retrotransposons found from various eukaryotes. Cassandra is a unique TRIM lineage which contains a 5S rRNA-derived sequence in its LTRs. Here, two new groups of TRIMs, designated Helenus and Ajax, are reported based on bioinformatics analysis and the usage of Repbase. Helenus is found from fungi, animals, and plants, and its LTRs contain a tRNA-like sequence. It includes two LTRs and between them, a primer-binding site (PBS) and polypurine tract (PPT) exist. Fungal and plant Helenus generate 5 bp target site duplications (TSDs) upon integration, while animal Helenus generates 4 bp TSDs. Ajax includes a 5S rRNA-derived sequence in its LTR and is found from two nemertean genomes. Ajax generates 5 bp TSDs upon integration. These results suggest that despite their unique promoters, Helenus and Ajax are TRIMs whose transposition is dependent on autonomous LTR retrotransposon. These TRIMs can originate through an insertion of SINE in an LTR of TRIM. The discovery of Helenus and Ajax suggests the presence of TRIMs with a promoter for RNA polymerase III derived from a small RNA gene, which is here collectively termed TRIMp3.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

Reminiscences on the honeybee genome project and the rise of epigenetic concepts in insect science

The sequencing of the honeybee genome in 2006 was an important technological and logistic achievement experience. But what benefits have flown from the honeybee genome project? What does the annotated genomic assembly mean for the study of behavioural complexity and organismal function in honeybees? Here, I discuss several lines of research that have arisen from this project and highlight the rapidly expanding studies on insect epigenomics, emergent properties of royal jelly, the mechanism of nutritional control of development and the contribution of epigenomic regulation to the evolution of sociality. I also argue that the term 'insect epigenetics' needs to be carefully redefined to reflect the diversity of epigenomic toolkits in insects and the impact of lineage-specific innovations on organismal outcomes. The honeybee genome project helped pioneer advances in social insect molecular biology, and fuelled breakthrough research into the role of flexible epigenomic control systems in linking genotype to phenotype.

The response of heat shock proteins in honey bees to abiotic and biotic stressors

Honey bees, Apis mellifera, are the most important managed pollinators worldwide. They are highly impacted by various abiotic and biotic stressors, especially temperature extremes, which can lead to cellular damage and death. The induction of heat shock proteins (HSPs) has been recorded in honey bees as a response to various types of stressors. HSPs are classified into different gene families according to their molecular weights. HSPs play an important role in maintaining cellular protein homeostasis due to their contribution as molecular chaperones or co-chaperones. HSPs in honey bees have complex functions with induction even under normal colony conditions. Previous studies have suggested various functions of HSPs to protect cells from damage under exposure to environmental stressors, pollutants, and pathogens. Surprisingly, HSPs have also been found to play roles in larval development and age-related tasks. The expression of HSPs varies depending on tissue type, developmental stage, age, and stress period. This article reviews studies on HSPs (sHSPs, HSP40, HSP60, HSP70, and HSP90) in honey bees and highlights gaps in the available knowledge. This review is crucial for honey bee research, particularly in the face of climate change challenges.

Comparative Genome Annotation

Newly sequenced genomes are being added to the tree of life at an unprecedented fast pace. A large proportion of such new genomes are phylogenetically close to previously sequenced and annotated genomes. In other cases, whole clades of closely related species or strains ought to be annotated simultaneously. Often, in subsequent studies, differences between the closely related species or strains are in the focus of research when the shared gene structures prevail. We here review methods for comparative structural genome annotation. The reviewed methods include classical approaches such as the alignment of protein sequences or protein profiles against the genome and comparative gene prediction methods that exploit a genome alignment to annotate either a single target genome or all input genomes simultaneously. We discuss how the methods depend on the phylogenetic placement of genomes, give advice on the choice of methods, and examine the consistency between gene structure annotations in an example. Furthermore, we provide practical advice on genome annotation in general.

Toll receptor ligand Spätzle 4 responses to the highly pathogenic Enterococcus faecalis from Varroa mites in honeybees

Honeybees play a major role in crop pollination, which supports the agricultural economy and international food supply. The colony health of honeybees is threatened by the parasitic mite Varroa destructor, which inflicts physical injury on the hosts and serves as the vector for variable viruses. Recently, it shows that V. destructor may also transmit bacteria through the feeding wound, yet it remains unclear whether the invading bacteria can exhibit pathogenicity to the honeybees. Here, we incidentally isolate Enterococcus faecalis, one of the most abundant bacteria in Varroa mites, from dead bees during our routine generation of microbiota-free bees in the lab. In vivo tests show that E. faecalis is only pathogenic in Apis mellifera but not in Apis cerana. The expression of antimicrobial peptide genes is elevated following infection in A. cerana. The gene-based molecular evolution analysis identifies positive selection of genes encoding Späetzle 4 (Spz4) in A. cerana, a signaling protein in the Toll pathway. The amino acid sites under positive selection are related to structural changes in Spz4 protein, suggesting improvement of immunity in A. cerana. The knock-down of Spz4 in A. cerana significantly reduces the survival rates under E. faecalis challenge and the expression of antimicrobial peptide genes. Our results indicate that bacteria associated with Varroa mites are pathogenic to adult bees, and the positively selected gene Spz4 in A. cerana is crucial in response to this mite-related pathogen.

A postmeiotically bifurcated roadmap of honeybee spermatogenesis marked by phylogenetically restricted genes

Haploid males of hymenopteran species produce gametes through an abortive meiosis I followed by meiosis II that can either be symmetric or asymmetric in different species. Thus, one spermatocyte could give rise to two spermatids with either equal or unequal amounts of cytoplasm. It is currently unknown what molecular features accompany these postmeiotic sperm cells especially in species with asymmetric meiosis II such as bees. Here we present testis single-cell RNA sequencing datasets from the honeybee (Apis mellifera) drones of 3 and 14 days after emergence (3d and 14d). We show that, while 3d testes exhibit active, ongoing spermatogenesis, 14d testes only have late-stage spermatids. We identify a postmeiotic bifurcation in the transcriptional roadmap during spermatogenesis, with cells progressing toward the annotated spermatids (SPT) and small spermatids (sSPT), respectively. Despite an overall similarity in their transcriptomic profiles, sSPTs express the fewest genes and the least RNA content among all the sperm cell types. Intriguingly, sSPTs exhibit a relatively high expression level for Hymenoptera-restricted genes and a high mutation load, suggesting that the special meiosis II during spermatogenesis in the honeybee is accompanied by phylogenetically young gene activities.

Novel Insights into the circRNA-Modulated Developmental Mechanism of Western Honey Bee Larval Guts

Circular RNAs (circRNAs) are a class of novel non-coding RNAs (ncRNAs) that play essential roles in the development and growth of vertebrates through multiple manners. However, the mechanism by which circRNAs modulate the honey bee gut development is currently poorly understood. Utilizing the transcriptome data we obtained earlier, the highly expressed circRNAs in the Apis mellifera worker 4-, 5-, and 6-day-old larval guts were analyzed, which was followed by an in-depth investigation of the expression pattern of circRNAs during the process of larval guts development and the potential regulatory roles of differentially expressed circRNAs (DEcircRNAs). In total, 1728 expressed circRNAs were detected in the A. mellifera larval guts. Among the most highly expressed 10 circRNAs, seven (novel_circ_000069, novel_circ_000027, novel_circ_000438, etc.) were shared by the 4-, 5-, and 6-day-old larval guts. In addition, 21 (46) up-regulated and 22 (27) down-regulated circRNAs were, respectively, screened in the Am4 vs. Am5 (Am5 vs. Am6) comparison groups. Additionally, nine DEcircRNAs, such as novel_circ_000340, novel_circ_000758 and novel_circ_001116, were shared by these two comparison groups. These DEcircRNAs were predicted to be transcribed from 14 and 29 parental genes; these were respectively annotated to 15 and 22 GO terms such as biological regulation and catalytic activity as well as 16 and 21 KEGG pathways such as dorsoventral axis formation and apoptosis. Moreover, a complicated competing endogenous RNA (ceRNA) network was observed; novel_circ_000838 in the Am4 vs. Am5 comparison group potentially targeted ame-miR-6000a-3p, further targeting 518 mRNAs engaged in several developmental signaling pathways (e.g., TGF-beta, hedgehog, and wnt signaling pathway) and immune pathways (e.g., phagosome, lysosome, and MAPK signaling pathway). The results demonstrated that the novel_circ_000838-ame-miR-6000a-3p axis may plays a critical regulatory part in the larval gut development and immunity. Furthermore, back-splicing sites of six randomly selected DEcircRNAs were amplified and verified by PCR; an RT-qPCR assay of these six DEcircRNAs confirmed the reliability of the used high-throughput sequencing data. Our findings provide a novel insight into the honey bee gut development and pave a way for illustration of the circRNA-modulated developmental mechanisms underlying the A. mellifera worker larval guts.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Cellular heterogeneity of the developing worker honey bee (Apis mellifera) pupa: a single cell transcriptomics analysis

It is estimated that animals pollinate 87.5% of flowering plants worldwide and that managed honey bees (Apis mellifera) account for 30-50% of this ecosystem service to agriculture. In addition to their important role as pollinators, honey bees are well-established insect models for studying learning and memory, behavior, caste differentiation, epigenetic mechanisms, olfactory biology, sex determination, and eusociality. Despite their importance to agriculture, knowledge of honey bee biology lags behind many other livestock species. In this study, we have used scRNA-Seq to map cell types to different developmental stages of the worker honey bee (prepupa at day 11 and pupa at day 15) and sought to determine their gene expression signatures. To identify cell-type populations, we examined the cell-to-cell network based on the similarity of the single-cells transcriptomic profiles. Grouping similar cells together we identified 63 different cell clusters of which 17 clusters were identifiable at both stages. To determine genes associated with specific cell populations or with a particular biological process involved in honey bee development, we used gene coexpression analysis. We combined this analysis with literature mining, the honey bee protein atlas, and gene ontology analysis to determine cell cluster identity. Of the cell clusters identified, 17 were related to the nervous system and sensory organs, 7 to the fat body, 19 to the cuticle, 5 to muscle, 4 to compound eye, 2 to midgut, 2 to hemocytes, and 1 to malpighian tubule/pericardial nephrocyte. To our knowledge, this is the first whole single-cell atlas of honey bees at any stage of development and demonstrates the potential for further work to investigate their biology at the cellular level.

The NIH Comparative Genomics Resource: addressing the promises and challenges of comparative genomics on human health

Comparative genomics is the comparison of genetic information within and across organisms to understand the evolution, structure, and function of genes, proteins, and non-coding regions (Sivashankari and Shanmughavel, Bioinformation 1:376-8, 2007). Advances in sequencing technology and assembly algorithms have resulted in the ability to sequence large genomes and provided a wealth of data that are being used in comparative genomic analyses. Comparative analysis can be leveraged to systematically explore and evaluate the biological relationships and evolution between species, aid in understanding the structure and function of genes, and gain a better understanding of disease and potential drug targets. As our knowledge of genetics expands, comparative genomics can help identify emerging model organisms among a broader span of the tree of life, positively impacting human health. This impact includes, but is not limited to, zoonotic disease research, therapeutics development, microbiome research, xenotransplantation, oncology, and toxicology. Despite advancements in comparative genomics, new challenges have arisen around the quantity, quality assurance, annotation, and interoperability of genomic data and metadata. New tools and approaches are required to meet these challenges and fulfill the needs of researchers. This paper focuses on how the National Institutes of Health (NIH) Comparative Genomics Resource (CGR) can address both the opportunities for comparative genomics to further impact human health and confront an increasingly complex set of challenges facing researchers.

The evolutionary history of bees in time and space

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

Apis mellifera Solinvivirus-1, a Novel Honey Bee Virus That Remained Undetected for over a Decade, Is Widespread in the USA

A metagenomic analysis of the virome of honey bees (Apis mellifera) from an apiary with high rates of unexplained colony losses identified a novel RNA virus. The virus, which was named Apis mellifera solinvivirus 1 (AmSV1), contains a 10.6 kb positive-strand genomic RNA with a single ORF coding for a polyprotein with the protease, helicase, and RNA-dependent RNA polymerase domains, as well as a single jelly-roll structural protein domain, showing highest similarity with viruses in the family Solinviviridae. The injection of honey bee pupae with AmSV1 preparation showed an increase in virus titer and the accumulation of the negative-strand of AmSV1 RNA 3 days after injection, indicating the replication of AmSV1. In the infected worker bees, AmSV1 was present in heads, thoraxes, and abdomens, indicating that this virus causes systemic infection. An analysis of the geographic and historic distribution of AmSV1, using over 900 apiary samples collected across the United States, showed AmSV1 presence since at least 2010. In the year 2021, AmSV1 was detected in 10.45% of apiaries (95%CI: 8.41-12.79%), mostly sampled in June and July in Northwestern and Northeastern United States. The diagnostic methods and information on the AmSV1 distribution will be used to investigate the connection of AmSV1 to honey bee colony losses.

Accelerated loss of diastase in mānuka honey: Investigation of mānuka specific compounds

Diastase is used internationally as a quality monitor for excessive heat treatment and prolonged storage of honey; honey must contain an activity of at least 8 diastase numbers (DN) for it to be considered export quality. Freshly harvested mānuka honey can have diastase activity close to the export threshold of 8 DN without excess heating, increasing susceptibility for export failure. This research investigated the effect of compounds unique to or high in concentration in mānuka honey on diastase activity. Investigation of the effect of methylglyoxal, dihydroxyacetone, 2-methoxybenzoic acid, 3-phenyllatic acid, 4-hydroxyphenyllactic acid and 2'-methoxyacetophenone on diastase activity was carried out. Mānuka honey was stored at 20 and 27 °C and clover honey spiked with compounds of interest were stored at 20, 27 and 34 °C and monitored overtime. Methylglyoxal and 3-phenyllactic acid were found to accelerate the loss of diastase above the loss normally observed with time and elevated temperature.

Expression profile of the entire detoxification gene inventory of the western honeybee, Apis mellifera across life stages

The western honeybee, Apis mellifera, is a managed pollinator of many crops and potentially exposed to a wide range of foreign compounds, including pesticides throughout its life cycle. Honeybees as well as other insects recruit molecular defense mechanisms to facilitate the detoxification of xenobiotic compounds. The inventory of detoxification genes (DETOXome) is comprised of five protein superfamilies: cytochrome P450 monooxygenases (P450), carboxylesterases, glutathione S-transferases (GST), UDP-glycosyl transferases (UGT) and ATP-binding cassette (ABC) transporters. Here we characterized the gene expression profile of the entire honeybee DETOXome by analyzing 47 transcriptomes across the honeybee life cycle, including different larval instars, pupae, and adults. All life stages were well separated by principal component analysis, and K-means clustering revealed distinct temporal patterns of gene expression. Indeed, >50% of the honeybee detoxification gene inventory is found in one cluster and follows strikingly similar expression profiles, i.e., increased expression during larval development, followed by a sharp decline after pupation and a steep increase again in adults. This cluster includes 29 P450 genes dominated by CYP3 and CYP4 clan members, 15 ABC transporter genes mostly belonging to the ABCC subfamily and 13 carboxylesterase genes including almost all members involved in dietary/detox and hormone/semiochemical processing. RT-qPCR analysis of selected detoxification genes from all families revealed high expression levels in various tissues, especially Malpighian tubules, fatbody and midgut, supporting the view that these tissues are essential for metabolic clearance of environmental toxins and pollutants in honeybees. Our study is meant to spark further research on the molecular basis of detoxification in this critical pollinator to better understand and evaluate negative impacts from potentially toxic substances. Additionally, the entire gene set of 47 transcriptomes collected and analyzed provides a valuable resource for future honeybee research across different disciplines.

A consensus phosphoserine within the large cytoplasmic loop of insect nAChR α8 subunits modulated interaction between 14-3-3ε and nAChRs to regulate neonicotinoid efficacy

Neonicotinoids are insect-selective nicotinic acetylcholine receptors (nAChRs) agonists that are used extensively for plant protection and animal health care. Some chaperone proteins, such as 14-3-3 proteins, importantly modulate nAChRs to display the physiological and pharmacological properties. Here we found that there is a 14-3-3 binding motif RSPSTH within the cytoplasmic loop of most insect α8 subunits. In the motif, a potential phosphorylated serine residue, serine 337, was a putative protein kinase A (PKA) substrate. Using Locusta migratoria α8 subunit as a representative, here we demonstrated that Loc14-3-3ε interacted with the unique phosphoserine (α8^S337) of Locα8 subunit to regulate agonist efficacy on hybrid Locα8/β2 nAChRs in Xenopus oocytes. Co-expression of Loc14-3-3ε caused a dramatic rise of maximal inward currents (I_max) of Locα8/β2 for acetylcholine and imidacloprid to 2.9-fold and 3.1-fold of that of Locα8/β2 alone. The S337A substitution of Locα8 reduced the I_max rise when Locα8^S337A/β2 and Loc14-3-3ε were co-expressed. The increased agonist currents by exogenous Loc14-3-3ε on Locα8/β2 could be almost abolished by either PKA inhibitor KT5720 or 14-3-3 inhibitor difopein. The findings revealed that serine 337 within motif RSPSTH was important for the interaction between insect nAChRs and 14-3-3ε, and inhibiting the interaction would change the pharmacological property of insect nAChRs to agonist such as neonicotinoids which may provide insights to develop new targets for insecticide design.

Positive selection in cytochrome P450 genes is associated with gonad phenotype and mating strategy in social bees

The honey bee, Apis mellifera differs from all other social bees in its gonad phenotype and mating strategy. Honey bee queens and drones have tremendously enlarged gonads, and virgin queens mate with several males. In contrast, in all the other bees, the male and female gonads are small, and the females mate with only one or very few males, thus, suggesting an evolutionary and developmental link between gonad phenotype and mating strategy. RNA-seq comparisons of A. mellifera larval gonads revealed 870 genes as differentially expressed in queens versus workers and drones. Based on Gene Ontology enrichment we selected 45 genes for comparing the expression levels of their orthologs in the larval gonads of the bumble bee Bombus terrestris and the stingless bee, Melipona quadrifasciata, which revealed 24 genes as differentially represented. An evolutionary analysis of their orthologs in 13 solitary and social bee genomes revealed four genes with evidence of positive selection. Two of these encode cytochrome P450 proteins, and their gene trees indicated a lineage-specific evolution in the genus Apis, indicating that cytochrome P450 genes may be involved in the evolutionary association of polyandry and the exaggerated gonad phenotype in social bees.

A comparative transcriptome analysis of the head of 1 and 9 days old worker honeybees (Apis mellifera)

The role of bees in the environment, economic, biodiversity and pharmaceutical industries is due to its social behavior, which is oriented from the brain and hypopharyngeal gland that is the center of royal jelly (RJ) production. Limited studies have been performed on the head gene expression profile at the RJ production stage. The aim of this study was to compare the gene expressions in 9 and 1-day-old (DO) honeybee workers in order to achieve better understanding about head gene expression pattern. After sequencing of RNAs, transcriptome and their networks were compared. The head expression profile undergoes various changes. 1662 gene transcripts had differential expressions which 1125 and 537 were up and down regulated, respectively, in 9_DO compared with 1_DO honey bees. The day 1th had more significant role in the expression of genes related to RJ production as major RJ protein 1, 2, 3, 5, 6 and 9 encoding genes, but their maximum secretion occurred at day 9th. All process related to hypopharyngeal glands activities as CYP450 gene, fatty acid synthase gene, vitamin B6 metabolism and some of genes involved in fatty acid elongation and degradation process had an upward trend from 1_DO and were age-dependent. By increasing the age, the activity of pathways related to immune system increased for keeping the health of bees against the chemical compound. The expression of aromatic amino acid genes involved in Phenylalanine, tyrosine and tryptophan biosynthesis pathway are essential for early stage of life. In 9_DO honeybees, the energy supplying, reducing stress, protein production and export pathways have a crucial role for support the body development and the social duties. It can be stated that the activity of honeybee head is focused on energy supply instead of storage, while actively trying to improve the level of cell dynamics for increasing the immunity and reducing stress. Results of current study identified key genes of certain behaviors of honeybee workers. Deeper considering of some pathways will be evaluated in future studies.

A Preliminary Study on “Personalised Treatment” against Varroa destructor Infestations in Honey Bee (Apis mellifera) Colonies

The requirement for the protection of bee colonies against Varroa destructor invasions has been noted by many breeders and is included as an aspect of the development of beekeeping. This research aimed to check the effect of the development of a colony exposed to laying eggs (brood surface) by queen bees with similar chemical potential (sisters) on the effect of a preparation combating V. destructor depending on the number of mites found in a given colony. We chose this as a standard model of conduct that treats each bee colony as one organism subjected to individual parasite control. For this purpose, we created a bee colony with a mother-of-one breeding line and fertilised drones from one colony. Infection with V. destructor occurred naturally and uncontrollably. Without interfering with the colony’s development, the frame insulator helped each colony’s brood (mothers’ reproductive potential) and the initial and final individuals from the mites themselves. The study was carried out in four species (two control species and two species with up to 20 and over 21 mites, respectively). Treatments with amitraz to combat damage were divided into four treatment subgroups: two treatments every four days or four treatments every two days. We observed the number of individuals that were protected in all subgroups in the average brood area. The reproductive potential of the sisters’ mothers did not change after the treatments with amitraz, which indicated that amitraz did not affect the delegation of egg laying. The invasion rate was also tracked relative to the control group, which allowed us to conclude that a two-time treatment with amitraz reduced the frequency of mites and a four-time treatment checked the effectiveness. Tailoring the control of V. destructor in bee colonies may be an effective measure in the fight against this parasite.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

Analysis of the Toll and Spaetzle Genes Involved in Toll Pathway-Dependent Antimicrobial Gene Induction in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)

Insects rely only on their innate immune system to protect themselves from pathogens. Antimicrobial peptide (AMP) production is the main immune reaction in insects. In Drosophila melanogaster, the reaction is regulated mainly by the Toll and immune deficiency (IMD) pathways. Spaetzle proteins, activated by immune signals from upstream components, bind to Toll proteins, thus, activating the Toll pathway, which in turn, induces AMP genes. Previous studies have shown the difference in immune systems related to Toll and IMD pathways between D. melanogaster and Tribolium castaneum. In T. castaneum, nine Toll and seven spaetzle (spz) genes were identified. To extend our understanding of AMP production by T. castaneum, we conducted functional assays of Toll and spaetzle genes related to Toll-pathway-dependent AMP gene expression in T. castaneum under challenge with bacteria or budding yeast. The results revealed that Toll3 and Toll4 double-knockdown and spz7 knockdown strongly and moderately reduced the Toll-pathway-dependent expression of AMP genes, respectively. Moreover, Toll3 and Toll4 double-knockdown pupae more rapidly succumbed to entomopathogenic bacteria than the control pupae, but spz7 knockdown pupae did not. The results suggest that Toll3 and Toll4 play a large role in Toll-pathway-dependent immune reactions, whereas spz7 plays a small part.

Current Knowledge on Bee Innate Immunity Based on Genomics and Transcriptomics

As important pollinators, bees play a critical role in maintaining the balance of the ecosystem and improving the yield and quality of crops. However, in recent years, the bee population has significantly declined due to various pathogens and environmental stressors including viruses, bacteria, parasites, and increased pesticide application. The above threats trigger or suppress the innate immunity of bees, their only immune defense system, which is essential to maintaining individual health and that of the colony. In addition, bees can be divided into solitary and eusocial bees based on their life traits, and eusocial bees possess special social immunities, such as grooming behavior, which cooperate with innate immunity to maintain the health of the colony. The omics approach gives us an opportunity to recognize the distinctive innate immunity of bees. In this regard, we summarize innate bee immunity from a genomic and transcriptomic perspective. The genetic characteristics of innate immunity were revealed by the multiple genomes of bees with different kinds of sociality, including honeybees, bumblebees, wasps, leaf-cutter bees, and so on. Further substantial transcriptomic data of different tissues from diverse bees directly present the activation or suppression of immune genes under the infestation of pathogens or toxicity of pesticides.

A genotyping by sequencing approach can disclose Apis mellifera population genomic information contained in honey environmental DNA

Awareness has been raised over the last years on the genetic integrity of autochthonous honey bee subspecies. Genomic tools available in Apis mellifera can make it possible to measure this information by targeting individual honey bee DNA. Honey contains DNA traces from all organisms that contributed or were involved in its production steps, including the honey bees of the colony. In this study, we designed and tested a genotyping by sequencing (GBS) assay to analyse single nucleotide polymorphisms (SNPs) of A. mellifera nuclear genome using environmental DNA extracted from honey. A total of 121 SNPs (97 SNPs informative for honey bee subspecies identification and 24 SNPs associated with relevant traits of the colonies) were used in the assay to genotype honey DNA, which derives from thousands of honey bees. Results were integrated with information derived from previous studies and whole genome resequencing datasets. This GBS method is highly reliable in estimating honey bee SNP allele frequencies of the whole colony from which the honey derived. This assay can be used to identify the honey bee subspecies of the colony that produced the honey and, in turn, to authenticate the entomological origin of the honey.

Transcriptome characterization and SSR discovery in the giant spiny frog Quasipaa spinosa

The giant spiny frog Quasipaa spinosa (Amphibia: Ranidae) is a large unique frog species found mainly in southern China with a low amount of fat and high protein, and it has become one of the most important aquaculture animal species in China. To better understand its genetic background and screen potential molecular markers for artificial breeding and species conservation, we constructed an expression profile of Q. spinosa with high-throughput RNA sequencing and acquired potential SSR markers. Approximately 81.7 Gb of data and 93,887 unigenes were generated. The transcriptome contains 2085 (80.7 %) complete BUSCOs, suggesting that our assembly methods were effective and accurate.These unigenes were functionally classified using 7 functional databases, yielding 17,482 Pfam-, 12,752 Sting-, 17,526 KEGG-, 24,341 Swiss-Prot-, 28,604 Nr-, 16,287 GO- and 12,752 COG-annotated unigenes. Among several amphibian species, Q. spinosa unigenes had the highest number of hits to Xenopus tropicalis (35.25 %), followed by Xenopus laevis (12.68 %). 1417 unigenes were assigned to the immune system. In addition, a total of 33,019 candidate SSR markers were identified from the constructed library. Further tests with 20 loci and 118 large-scale breeding specimens gathered from four culture farms in China showed that 15 (75 %) loci were polymorphic, with the number of alleles per locus varying from 3 to 9 (mean of 4.3). The PIC values for the SSR markers ranged from 0.19 to 0.82, with an average value of 0.43, indicating moderate polymorphism in Q. spinosa. The transcriptomic profile and SSR repertoire obtained in the present study will facilitate population genetic studies and the selective breeding of amphibian species.

Identification of sex-biased and neurodevelopment genes via brain transcriptome in Ostrinia furnacalis

Insect brains play important roles in the regulation of sex-biased behaviors such as mating and oviposition. The neural structure and function of brain differences between males and females have been identified, in which the antenna lobes (AL) showed the most discrepancy, however, the whole repertoire of the genes expressed in the brains and the molecular mechanism of neural signaling and structural development are still unclear. In this study, high-throughput transcriptome analysis of male and female brains was carried on in the Asia corn borer, Ostrinia furnacalis, and a total of 39.23 Gb data and 34,092 unigenes were obtained. Among them, 276 genes displayed sex-biased expression by DEG analysis, of which 125 genes were highly expressed in the males and 151 genes were highly expressed in the females. Besides, by homology analysis against genes that have been confirmed to be related to brain neurodevelopment, a total of 24 candidate genes were identified in O. furnacalis. In addition, to further screen the core genes that may be important for sex-biased nerve signaling and neurodevelopment, protein-protein interaction networks were constructed for the sex-biased genes and neurodevelopment genes. We identified 10 (Mhc, Mlc1, Mlc2, Prm, Mf, wupA, TpnC25D, fln, l(2)efl, and Act5C), 11 (PPO2, GNBP3, Spn77Ba, Ppn, yellow-d2, PGRP-LB, PGRP-SD, PGRP-SC2, Hml, Cg25C, and vkg) and 8 (dac, wg, hh, ci, run, Lim1, Rbp9, and Bx) core hub genes that may be related to brain neural development from male-biased, female-biased, and neurodevelopment gene groups. Our results provide a reference for further analysis of the dimorphism of male and female brain structures in agricultural pests.

Revealing Landscapes of Transposable Elements in Apis Species by Meta-Analysis

Transposable elements (TEs) are grouped into several families with diverse sequences. Owing to their diversity, studies involving the detection, classification, and annotation of TEs are difficult tasks. Moreover, simple comparisons of TEs among different species with different methods can lead to misinterpretations. The genome data of several honey bee (Apis) species are available in public databases. Therefore, we conducted a meta-analysis of TEs, using 11 sets of genome data for Apis species, in order to establish data of “landscape of TEs”. Consensus TE sequences were constructed and their distributions in the Apis genomes were determined. Our results showed that TEs belonged to four to seven TE families among 13 and 15 families of TEs detected in classes I and II respectively mainly consisted of Apis TEs and that more DNA/TcMar-Mariner consensus sequences and copies were present in all Apis genomes tested. In addition, more consensus sequences and copy numbers of DNA/TcMar-Mariner were detected in Apis mellifera than in other Apis species. These results suggest that TcMar-Mariner might exert A. mellifera-specific effects on the host A. mellifera species. In conclusion, our unified approach enabled comparison of Apis genome sequences to determine the TE landscape, which provide novel evolutionary insights into Apis species.

The genome of the rice planthopper egg parasitoid wasps Anagrus nilaparvatae casts light on the chemo- and mechanosensation in parasitism

BACKGROUND

Mymaridae is an ancient insect group and is a basal lineage of the superfamily Chalcidoidea. Species of Mymaridae have great potential for biological control. Anagrus nilaparvatae, a representative species of Mymaridae, is ideal for controlling rice planthopper due to its high rate of parasitism and ability to find hosts efficiently in paddy ridges and fields.

RESULTS

Using both PacBio single-molecule real-time and Illumina sequencing, we sequenced and assembled the whole genome of A. nilaparvatae, a first for the family Mymaridae. The assembly consists of 394 scaffolds, totaling 488.8 Mb. The assembly is of high continuity and completeness, indicated by the N50 value of 25.4 Mb and 98.2% mapping rate of Benchmarking Universal Single-Copy Orthologs. In total, 16,894 protein-coding genes in the genome were annotated. A phylogenomic tree constructed for A. nilaparvatae and other 12 species of Hymenoptera confirmed that the family Mymaridae is sister to all remaining chalcidoids. The divergence time between A. nilaparvatae and the other seven Chalcidoidea species was dated at ~ 126.9 Mya. Chemoreceptor and mechanoreceptor genes are important in explaining parasitic behavior. We identified 17 odorant binding proteins, 11 chemosensory proteins, four Niemann-Pick type C2 proteins, 88 olfactory receptors, 12 gustatory receptors, 22 ionotropic receptors and 13 sensory neuron membrane proteins in the genome of A. nilaparvatae, which are associated with the chemosensory functions. Strikingly, there is only one pickpocket receptors and nine transient receptor potential genes in the genome that have a mechanosensory function.

CONCLUSIONS

We obtained a high-quality genome assembly for A. nilaparvatae using PacBio single-molecule real-time sequencing, which provides phylogenomic insights for its evolutionary history. The small numbers of chemo- and mechanosensory genes in A. nilaparvatae indicate the species-specific host detection and oviposition behavior of A. nilaparvatae might be regulated by relatively simple molecular pathways.

Behavioral genetics and genomics: Mendel's peas, mice, and bees

The question of the heritability of behavior has been of long fascination to scientists and the broader public. It is now widely accepted that most behavioral variation has a genetic component, although the degree of genetic influence differs widely across behaviors. Starting with Mendel's remarkable discovery of "inheritance factors," it has become increasingly clear that specific genetic variants that influence behavior can be identified. This goal is not without its challenges: Unlike pea morphology, most natural behavioral variation has a complex genetic architecture. However, we can now apply powerful genome-wide approaches to connect variation in DNA to variation in behavior as well as analyses of behaviorally related variation in brain gene expression, which together have provided insights into both the genetic mechanisms underlying behavior and the dynamic relationship between genes and behavior, respectively, in a wide range of species and for a diversity of behaviors. Here, we focus on two systems to illustrate both of these approaches: the genetic basis of burrowing in deer mice and transcriptomic analyses of division of labor in honey bees. Finally, we discuss the troubled relationship between the field of behavioral genetics and eugenics, which reminds us that we must be cautious about how we discuss and contextualize the connections between genes and behavior, especially in humans.

Experimental disruption of social structure reveals totipotency in the orchid bee, Euglossa dilemma

Eusociality has evolved multiple times across the insect phylogeny. Social insects with greater levels of social complexity tend to exhibit specialized castes with low levels of individual phenotypic plasticity. In contrast, species with simple social groups may consist of totipotent individuals that transition among behavioral and reproductive states. However, recent work has shown that in simple social groups, there can still be constraint on individual plasticity, caused by differences in maternal nourishment or social interaction. It is not well understood how these constraints arise, ultimately leading to the evolution of nonreproductive workers. Some species of orchid bees form social groups of a dominant and-one to two subordinate helpers where all individuals are reproductive. Females can also disperse to start their own nest as a solitary foundress, which includes a nonreproductive phase characterized by ovary inactivation, not typically expressed by subordinates. Little is known about individual flexibility across these trajectories. Here, using the orchid bee Euglossa dilemma, we assess the plasticity of subordinate helpers, finding that they are capable of the same behavioral, physiological, transcriptomic, and chemical changes seen in foundresses. Our results suggest that the lack of nonreproductive workers in E. dilemma is not due to a lack of subordinate plasticity.

Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging

Pollinator populations, including bees, are in rapid decline in many parts of the world, raising concerns over the future of ecosystems and food production. Among the factors involved in these declines, poor nutrition deserves attention. The diet consumed by adult worker honeybees (Apis mellifera) is crucial for their behavioral maturation, i.e., the progressive division of labor they perform, such as nurse bees initially and later in life as foragers. Poor pollen nutrition is known to reduce the workers' lifespan, but the underlying physiological and genetic mechanisms are not fully understood. Here we investigate how the lack of pollen in the diet of workers during their first week of adult life can affect age-related phenotypes. During the first seven days of adult life, newly emerged workers were fed either a pollen-deprived (PD) diet mimicking that of an older bee, or a control pollen-rich (PR) diet, as typically consumed by young bees. The PD-fed bees showed alterations in their fat body transcriptome, such as a switch from a protein-lipid based metabolism to a carbohydrate-based metabolism, and a reduced expression of genes involved with immune response. The absence of pollen in the diet also led to an accumulation of oxidative stress markers in fat body tissue and alterations in the cuticular hydrocarbon profiles, which became similar to those of chronologically older bees. Together, our data indicate that the absence of pollen during first week of adulthood triggers the premature onset of an aging-related worker phenotype.

G protein-coupled receptors (GPCRs) in rotifers and cladocerans: Potential applications in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology

The G protein-coupled receptor (GPCR) superfamily plays a fundamental role in both sensory functions and the regulation of homeostasis, and is highly conserved across the eukaryote taxa. Its functional diversity is related to a conserved seven-transmembrane core and invariant set of intracellular signaling mechanisms. The interplay between these properties is key to the evolutionary success of GPCR. As this superfamily originated from a common ancestor, GPCR genes have evolved via lineage-specific duplications through the process of adaptation. Here we summarized information on GPCR gene families in rotifers and cladocerans based on their evolutionary position in aquatic invertebrates and their potential application in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology. Phylogenetic analyses were conducted to examine the evolutionary significance of GPCR gene families and to provide structural insight on their role in aquatic invertebrates. In particular, most GPCR gene families have undergone sporadic evolutionary processes, but some GPCRs are highly conserved across species despite the dynamics of GPCR evolution. Overall, this review provides a better understanding of GPCR evolution in aquatic invertebrates and expand our knowledge of the potential application of these receptors in various fields.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

Non-additive gene interactions underpin molecular and phenotypic responses in honey bee larvae exposed to imidacloprid and thymol

Understanding the cumulative risk of chemical mixtures at environmentally realistic concentrations is a key challenge in honey bee ecotoxicology. Ecotoxicogenomics, including transcriptomics, measures responses in individual organisms at the molecular level which can provide insights into the mechanisms underlying phenotypic responses induced by one or more stressors and link impacts on individuals to populations. Here, fifth instar honey bee larvae were sampled from a previously reported field experiment exploring the phenotypic impacts of environmentally realistic chronic exposures of the pesticide imidacloprid (5 μg.kg^-1 for six weeks) and the acaricide thymol (250 g.kg^-1 applied via Apiguard gel in-hive for four weeks), both separately and in combination. RNA-seq was used to discover individual and interactive chemical effects on larval gene expression and to uncover molecular mechanisms linked to reported adult and colony phenotypes. The separate and combined treatments had distinct gene expression profiles which represented differentially affected signaling and metabolic pathways. The molecular signature of the mixture was characterised by additive interactions in canonical stress responses associated with oxidative stress and detoxification, and non-additive interactions in secondary responses including developmental, neurological, and immune pathways. Novel emergent impacts on eye development genes correlated with long-term defects in visual learning performance as adults. This is consistent with these chemicals working through independent modes of action that combine to impact common downstream pathways, and highlights the importance of establishing mechanistic links between molecular and phenotypic responses when predicting effects of chemical mixtures on ecologically relevant population outcomes.

iLoci: robust evaluation of genome content and organization for provisional and mature genome assemblies

We introduce a new framework for genome analyses based on parsing an annotated genome assembly into distinct interval loci (iLoci), available as open-source software as part of the AEGeAn Toolkit (https://github.com/BrendelGroup/AEGeAn). We demonstrate that iLoci provide an alternative coordinate system that is robust to changes in assembly and annotation versions and facilitates granular quality control of genome data. We discuss how statistics computed on iLoci reflect various characteristics of genome content and organization and illustrate how these statistics can be used to establish a baseline for assessment of the completeness and accuracy of the data. We also introduce a well-defined measure of relative genome compactness and compute other iLocus statistics that reveal genome-wide characteristics of gene arrangements in the whole genome context. Given the fast pace of assembly/annotation updates, our AEGeAn Toolkit fills a niche in computational genomics based on deriving persistent and species-specific genome statistics. Gene structure model-centric iLoci provide a precisely defined coordinate system that can be used to store assembly/annotation updates that reflect either stable or changed assessments. Large-scale application of the approach revealed species- and clade-specific genome organization in precisely defined computational terms, promising intriguing forays into the forces of shaping genome structure as more and more genome assemblies are being deposited.

Annotation of glycolysis, gluconeogenesis, and trehaloneogenesis pathways provide insight into carbohydrate metabolism in the Asian citrus psyllid

Citrus greening disease is caused by the pathogen Candidatus Liberibacter asiaticus and transmitted by the Asian citrus psyllid, Diaphorina citri. No curative treatment or significant prevention mechanism exists for this disease, which causes economic losses from reduced citrus production. A high-quality genome of D. citri is being manually annotated to provide accurate gene models to identify novel control targets and increase understanding of this pest. Here, we annotated 25 D. citri genes involved in glycolysis and gluconeogenesis, and seven in trehaloneogenesis. Comparative analysis showed that glycolysis genes in D. citri are highly conserved but copy numbers vary. Analysis of expression levels revealed upregulation of several enzymes in the glycolysis pathway in the thorax, consistent with the primary use of glucose by thoracic flight muscles. Manually annotating these core metabolic pathways provides accurate genomic foundation for developing gene-targeting therapeutics to control D. citri.

Genomic identification, annotation, and comparative analysis of Vacuolar-type ATP synthase subunits in Diaphorina citri

The hemipteran insect Diaphorina citri, or Asian citrus psyllid, is a vector for Candidatus Liberibacter asiaticus (CLas), the bacterium causing citrus greening disease, or Huanglongbing (HLB). Millions of citrus trees have been destroyed, and every grove in Florida, USA, has been directly affected by this disease. In eukaryotes, vacuolar-type ATP synthase (V-ATPase) is an abundant heterodimeric enzyme that serves the cell with essential compartment acidification through the active processes that transport protons across the membrane. Fifteen putative V-ATPase genes in the D. citri genome were manually curated. Comparative genomic analysis revealed that D. citri V-ATPase subunits share domains and motifs with other insects, including the V-ATPase-A superfamily domain. Phylogenetic analysis separates D. citri V-ATPase subunits into expected clades with orthologous sequences. Annotation of the D. citri genome is a critical step towards developing directed pest management strategies to reduce the spread of HLB throughout the citrus industry.

Population Structure and Diversity in European Honey Bees (Apis mellifera L.)—An Empirical Comparison of Pool and Individual Whole-Genome Sequencing

Background: Whole-genome sequencing has become routine for population genetic studies. Sequencing of individuals provides maximal data but is rather expensive and fewer samples can be studied. In contrast, sequencing a pool of samples (pool-seq) can provide sufficient data, while presenting less of an economic challenge. Few studies have compared the two approaches to infer population genetic structure and diversity in real datasets. Here, we apply individual sequencing (ind-seq) and pool-seq to the study of Western honey bees (Apis mellifera). Methods: We collected honey bee workers that belonged to 14 populations, including 13 subspecies, totaling 1347 colonies, who were individually (139 individuals) and pool-sequenced (14 pools). We compared allele frequencies, genetic diversity estimates, and population structure as inferred by the two approaches. Results: Pool-seq and ind-seq revealed near identical population structure and genetic diversities, albeit at different costs. While pool-seq provides genome-wide polymorphism data at considerably lower costs, ind-seq can provide additional information, including the identification of population substructures, hybridization, or individual outliers. Conclusions: If costs are not the limiting factor, we recommend using ind-seq, as population genetic structure can be inferred similarly well, with the advantage gained from individual genetic information. Not least, it also significantly reduces the effort required for the collection of numerous samples and their further processing in the laboratory.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Thrice out of Asia and the adaptive radiation of the western honey bee

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world’s most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic “hotspots.” We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Tools and applications for integrative analysis of DNA methylation in social insects

DNA methylation is a common epigenetic signalling tool and an important biological process which is widely studied in a large array of species. The presence, level and function of DNA methylation vary greatly across species. In some insects, DNA methylation systems are minimal, and overall methylation rates tend to be low in all studied insect species. Low methylation levels probed by whole-genome bisulphite sequencing require great care with respect to data quality control and interpretation. Here, we introduce BWASP/R, a complete workflow that allows efficient, scalable and entirely reproducible analyses of raw DNA methylation sequencing data. Consistent application of quality control filters and analysis parameters provides fair comparisons among different studies and an integrated view of all experiments on one species. We describe the capabilities of the BWASP/R workflow by re-analysing several publicly available social insect WGBS data sets, comprising 70 samples and cumulatively 147 replicates from four different species. We show that the CpG methylome comprises only about 1.5% of CpG sites in the honeybee genome and that the cumulative data are consistent with genetic signatures of site accessibility and physiological control of methylation levels.

Halofuginone triggers a transcriptional program centered on ribosome biogenesis and function in honey bees

We previously found that pharmacological inhibition of prolyl-tRNA synthetase by halofuginone has potent activity against Nosema ceranae, an important pathogen of honey bees. However, we also observed that prolyl-tRNA synthetase inhibition is toxic to bees, suggesting further work is necessary to make this a feasible therapeutic strategy. As expected, we found that pharmacological inhibition of prolyl-tRNA synthetase activity resulted in robust induction of select canonical ATF4 target genes in honey bees. However, our understanding of this and other cellular stress responses in general in honey bees is incomplete. Thus, we used RNAseq to identify novel changes in gene expression after halofuginone treatment and observed induction of genes involved in ribosome biogenesis, translation, tRNA synthesis, and ribosome-associated quality control (RQC). These results suggest that halofuginone, potentially acting through the Integrated Stress Response (ISR), promotes a transcriptional response to ribosome functional impairment in honey bees rather than the response designed to oppose amino acid limitation, which has been observed in other organisms after ISR induction. In support of this idea, we found that cycloheximide (CHX) administration also induced all tested target genes, indicating that this gene expression program could be induced by ribosome stalling in addition to tRNA synthetase inhibition. Only a subset of halofuginone-induced genes was upregulated by Unfolded Protein Response (UPR) induction, suggesting that mode of activation and cross-talk with other cellular signaling pathways significantly influence ISR function and cellular response to its activation. Future work will focus on understanding how the apparently divergent transcriptional output of the ISR in honey bees impacts the health and disease of this important pollinator species.