Population genomics of the invasive Northern Giant Hornet Vespa mandarinia in North America and across its native range.

The northern giant hornet Vespa mandarinia (NGH) is a voracious predator of other insect species, including honey bees. NGH's native range spans subtropical and temperate regions across much of east and southeast Asia and, in 2019, exotic populations of the species were discovered in North America. Despite this broad range and invasive potential, investigation of the population genomic structure of NGH across its native and introduced ranges has thus far been limited to a small number of mitochondrial samples. Here, we present analyses of genomic data from NGH individuals collected across the species' native range and from exotic individuals collected in North America. We provide the first survey of whole-genome population variation for any hornet species, covering this species' native and invasive ranges, and in doing so confirm likely origins in Japan and South Korea for the two introductions. We additionally show that, while this introduced population exhibited strongly elevated levels of inbreeding, these signatures of inbreeding are also present in some long-standing native populations, which may indicate that inbreeding depression alone is insufficient to prevent the persistence of NGH populations. As well as highlighting the importance of ongoing monitoring and eradication efforts to limit the spread of this species outside of its natural range, our data will serve as a foundational database for future genomic studies into introduced hornet populations.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality.

BACKGROUND: Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS: Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS: Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Rapid evolution, rearrangements and whole mitogenome duplication in the Australian stingless bees Tetragonula (Hymenoptera: Apidae): A steppingstone towards understanding mitochondrial function and evolution.

The extreme conservation of mitochondrial genomes in metazoans poses a significant challenge to understanding mitogenome evolution. However, the presence of variation in gene order or genome structure, found in a small number of taxa, can provide unique insights into this evolution. Previous work on two stingless bees in the genus Tetragonula (T. carbonaria and T. hockingsi) revealed highly divergent CO1 regions between them and when compared to the bees from the same tribe (Meliponini), indicating rapid evolution. Using mtDNA isolation and Illumina sequencing, we elucidated the mitogenomes of both species. In both species, there has been a duplication of the whole mitogenome to give a total genome size of 30,666 bp in T. carbonaria; and 30,662 bp in T. hockingsi. These duplicated genomes present a circular structure with two identical and mirrored copies of all 13 protein coding genes and 22 tRNAs, with the exception of a few tRNAs that are present as single copies. In addition, the mitogenomes are characterized by rearrangements of two block of genes. We believe that rapid evolution is present in the whole Indo-Malay/Australasian group of Meliponini but is extraordinarily elevated in T. carbonaria and T. hockingsi, probably due to founder effect, low effective population size and the mitogenome duplication. All these features - rapid evolution, rearrangements, and duplication - deviate significantly from the vast majority of the mitogenomes described so far, making the mitogenomes of Tetragonula unique opportunities to address fundamental questions of mitogenome function and evolution.

CRISPR-induced DNA reorganization for multiplexed nucleic acid detection.

Nucleic acid sensing powered by the sequence recognition of CRIPSR technologies has enabled major advancement toward rapid, accurate and deployable diagnostics. While exciting, there are still many challenges facing their practical implementation, such as the widespread need for a PAM sequence in the targeted nucleic acid, labile RNA inputs, and limited multiplexing. Here we report FACT (Functionalized Amplification CRISPR Tracing), a CRISPR-based nucleic acid barcoding technology compatible with Cas12a and Cas13a, enabling diagnostic outputs based on cis- and trans-cleavage from any sequence. Furthermore, we link the activation of CRISPR-Cas12a to the expression of proteins through a Reprogrammable PAIRing system (RePAIR). We then combine FACT and RePAIR to create FACTOR (FACT on RePAIR), a CRISPR-based diagnostic, that we use to detect infectious disease in an agricultural use case: honey bee viral infection. With high specificity and accuracy, we demonstrate the potential of FACTOR to be applied to the sensing of any nucleic acid of interest.

The suppressive potential of a gene drive in populations of invasive social wasps is currently limited.

Social insects are very successful invasive species, and the continued increase of global trade and transportation has exacerbated this problem. The yellow-legged hornet, Vespa velutina nigrithorax (henceforth Asian hornet), is drastically expanding its range in Western Europe. As an apex insect predator, this hornet poses a serious threat to the honey bee industry and endemic pollinators. Current suppression methods have proven too inefficient and expensive to limit its spread. Gene drives might be an effective tool to control this species, but their use has not yet been thoroughly investigated in social insects. Here, we built a model that matches the hornet's life history and modelled the effect of different gene drive scenarios on an established invasive population. To test the broader applicability and sensitivity of the model, we also incorporated the invasive European paper wasp Polistes dominula. We find that, due to the haplodiploidy of social hymenopterans, only a gene drive targeting female fertility is promising for population control. Our results show that although a gene drive can suppress a social wasp population, it can only do so under fairly stringent gene drive-specific conditions. This is due to a combination of two factors: first, the large number of surviving offspring that social wasp colonies produce make it possible that, even with very limited formation of resistance alleles, such alleles can quickly spread and rescue the population. Second, due to social wasp life history, infertile individuals do not compete with fertile ones, allowing fertile individuals to maintain a large population size even when drive alleles are widespread. Nevertheless, continued improvements in gene drive technology may make it a promising method for the control of invasive social insects in the future.

Breeding honey bees (Apis mellifera L.) for low and high Varroa destructor population growth: Gene expression of bees performing grooming behavior.

Introduction: Social organisms, including honey bees (Apis mellifera L.), have defense mechanisms to control the multiplication and transmission of parasites and pathogens within their colonies. Self-grooming, a mechanism of behavioral immunity, seems to contribute to restrain the population growth of the ectoparasitic mite Varroa destructor in honey bee colonies. Because V. destructor is the most damaging parasite of honey bees, breeding them for resistance against the mite is a high priority of the beekeeping industry. Methods: A bidirectional breeding program to select honey bee colonies with low and high V. destructor population growth (LVG and HVG, respectively) was conducted. Having high and low lines of bees allowed the study of genetic mechanisms underlying self-grooming behavior between the extreme genotypes. Worker bees were classified into two categories: 'light groomers' and 'intense groomers'. The brains of bees from the different categories (LVG-intense, LVG-light, HVG-intense, and HVG-light) were used for gene expression and viral quantification analyses. Differentially expressed genes (DEGs) associated with the LVG and HVG lines were identified. Results: Four odorant-binding proteins and a gustatory receptor were identified as differentially expressed genes. A functional enrichment analysis showed 19 enriched pathways from a list of 219 down-regulated DEGs in HVG bees, including the Kyoto Encyclopedia of Genes and Genomes (KEGG) term of oxidative phosphorylation. Additionally, bees from the LVG line showed lower levels of Apis rhabdovirus 1 and 2, Varroa destructor virus -1 (VDV-1/DWV-B), and Deformed wing virus-A (DWV-A) compared to bees of the HVG line. The difference in expression of odorant-binding protein genes and a gustatory receptor between bee lines suggests a possible link between them and the perception of irritants to trigger rapid self-grooming instances that require the activation of energy metabolic pathways. Discussion: These results provide new insights on the molecular mechanisms involved in honey bee grooming behavior. Differences in viral levels in the brains of LVG and HVG bees showed the importance of investigating the pathogenicity and potential impacts of neurotropic viruses on behavioral immunity. The results of this study advance the understanding of a trait used for selective breeding, self-grooming, and the potential of using genomic assisted selection to improve breeding programs.

Unraveling the Mechanisms of the Apis mellifera Honeycomb Construction by 4D X-ray Microscopy.

Honeycomb is one of nature's best engineered structures. Even though it has inspired several modern engineering structures, an understanding of the process by which the hexagonal cells are formed in 3D space is lacking. Previous studies on the structure of the honeycomb are based on either 2D microscopy or by direct visual observations. As a result, several critical features of its microstructure and the precise mechanisms of its growth are not well understood. Using 4D X-ray microscopy, this study shows how individual and groups of honeycomb cells are formed. Cells grow additively from a corrugated central spine in a dynamic manner. The previously undocumented, corrugated spine contributes significantly to the comb's robust mechanical properties in all three dimensions. As cells grow, honey bees create a "coping," which this study shows to be the location where new wax material is deposited behind where compaction and densification take place. This is exemplified by pores in the wax observed at the coping and alternating rear junctions between the comb cells that arise from the additive building technique and the highly efficient cell packing methodology, respectively. Additional mechanisms for growth and formation are discussed and described.

Haploid and Sexual Selection Shape the Rate of Evolution of Genes across the Honey Bee (Apis mellifera L.) Genome.

Many species have separate haploid and diploid phases. Theory predicts that each phase should experience the effects of evolutionary forces (like selection) differently. In the haploid phase, all fitness-affecting alleles are exposed to selection, whereas in the diploid phase, those same alleles can be masked by homologous alleles. This predicts that selection acting on genes expressed in haploids should be more effective than diploid-biased genes. Unfortunately, in arrhenotokous species, this prediction can be confounded with the effects of sex-specific expression, as haploids are usually reproductive males. Theory posits that, when accounting for ploidal- and sex-specific expression, selection should be equally efficient on haploid- and diploid-biased genes relative to constitutive genes. Here, we used a multiomic approach in honey bees to quantify the evolutionary rates of haploid-biased genes and test the relative effects of sexual- and haploid-expression on molecular evolution. We found that 16% of the honey bee's protein-coding genome is highly expressed in haploid tissue. When accounting for ploidy and sex, haploid- and diploid-biased genes evolve at a lower rate than expected, indicating that they experience strong negative selection. However, the rate of molecular evolution of haploid-biased genes was higher than diploid-based genes. Genes associated with sperm storage are a clear exception to this trend with evidence of strong positive selection. Our results provide an important empirical test of theory outlining how selection acts on genes expressed in arrhenotokous species. We propose the haploid life history stage affects genome-wide patterns of diversity and divergence because of both sexual and haploid selection.

Population Genetics of Reproductive Genes in Haplodiploid Species.

Many animal species are haplodiploid: their fertilized eggs develop into diploid females and their unfertilized eggs develop into haploid males. The unique genetic features of haplodiploidy raise the prospect that these systems can be used to disentangle the population genetic consequences of haploid and diploid selection. To this end, sex-specific reproductive genes are of particular interest because, while they are shared within the same genome, they consistently experience selection in different ploidal environments. However, other features of these genes, including sex-specific expression and putative involvement in postcopulatory sexual selection, are potentially confounding factors because they may also impact the efficacy of selection asymmetrically between the sexes. Thus, to properly interpret evolutionary genomic patterns, it is necessary to generate a null expectation for the relative amount of polymorphism and divergence we expect to observe among sex-specific genes in haplodiploid species, given differences in ploidal environment, sex-limited expression, and their potential role in sexual selection. Here, we derive the theoretical expectation for the rate of evolution of sex-specific genes in haplodiploid species, under the assumption that they experience the same selective environment as genes expressed in both sexes. We find that the null expectation is that reproductive genes evolve more rapidly than constitutively expressed genes in haplodiploid genomes. However, despite the aforementioned differences, the null expectation does not differ between male- and female-specific reproductive genes, when assuming additivity. Our theoretical results provide an important baseline expectation that should be used in molecular evolution studies comparing rates of evolution among classes of genes in haplodiploid species.

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

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

Prospects in Connecting Genetic Variation to Variation in Fertility in Male Bees.

Bees are economically and ecologically important pollinating species. Managed and native bee species face increasing pressures from human-created stressors such as habitat loss, pesticide use, and introduced pathogens. There has been increasing attention towards how each of these factors impacts fertility, especially sperm production and maintenance in males. Here, we turn our attention towards another important factor impacting phenotypic variation: genetics. Using honey bees as a model, we explore the current understanding of how genetic variation within and between populations contributes to variation in sperm production, sperm maintenance, and insemination success among males. We conclude with perspectives and future directions in the study of male fertility in honey bees and non-Apis pollinators more broadly, which still remain largely understudied.

A gene drive does not spread easily in populations of the honey bee parasite Varroa destructor.

Varroa mites (Varroa destructor) are the most significant threat to beekeeping worldwide. They are directly or indirectly responsible for millions of colony losses each year. Beekeepers are somewhat able to control varroa populations through the use of physical and chemical treatments. However, these methods range in effectiveness, can harm honey bees, can be physically demanding on the beekeeper, and do not always provide complete protection from varroa. More importantly, in some populations varroa mites have developed resistance to available acaricides. Overcoming the varroa mite problem will require novel and targeted treatment options. Here, we explore the potential of gene drive technology to control varroa. We show that spreading a neutral gene drive in varroa is possible but requires specific colony-level management practices to overcome the challenges of both inbreeding and haplodiploidy. Furthermore, continued treatment with acaricides is necessary to give a gene drive time to fix in the varroa population. Unfortunately, a gene drive that impacts female or male fertility does not spread in varroa. Therefore, we suggest that the most promising way forward is to use a gene drive which carries a toxin precursor or removes acaricide resistance alleles. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13592-021-00891-5.

Eusociality influences the strength of negative selection on insect genomes.

While much of the focus of sociobiology concerns identifying genomic changes that influence social behaviour, we know little about the consequences of social behaviour on genome evolution. It has been hypothesized that social evolution can influence the strength of negative selection via two mechanisms. First, division of labour can influence the efficiency of negative selection in a caste-specific manner; indirect negative selection on worker traits is theoretically expected to be weaker than direct selection on queen traits. Second, increasing social complexity is expected to lead to relaxed negative selection because of its influence on effective population size. We tested these two hypotheses by estimating the strength of negative selection in honeybees, bumblebees, paper wasps, fire ants and six other insects that span the range of social complexity. We found no consistent evidence that negative selection was significantly stronger on queen-biased genes relative to worker-biased genes. However, we found strong evidence that increased social complexity reduced the efficiency of negative selection. Our study clearly illustrates how changes in behaviour can influence patterns of genome evolution by modulating the strength of natural selection.

Defense Response in Brazilian Honey Bees (Apis mellifera scutellata × spp.) Is Underpinned by Complex Patterns of Admixture.

In 1957, an invasive and highly defensive honey bee began to spread across Brazil. In the previous year, Brazilian researchers hoped to produce a subtropical-adapted honey bee by crossing local commercial honey bees (of European origin) with a South African honey bee subspecies (Apis mellifera scutellata; an A-lineage honey bee subspecies). The resulting cross-African hybrid honey bees (AHBs)-escaped from their enclosure and spread through the Americas. Today, AHB is the most common honey bee from Northern Argentina to the Southern United States. AHBs are much more likely to sting nest intruders than managed European-derived honey bee colonies. Previous studies have explored how genetic variation contributes to differences in defense response between European-derived honey bee and AHB. Although this work demonstrated very strong genetic effects on defense response, they have yet to pinpoint which genes influence variation in defense response within AHBs, specifically. We quantified defense response for 116 colonies in Brazil and performed pooled sequencing on the most phenotypically divergent samples. We identified 65 loci containing 322 genes that were significantly associated with defense response. Loci were strongly associated with metabolic function, consistent with previous functional genomic analyses of this phenotype. Additionally, defense-associated loci had nonrandom and unexpected patterns of admixture. Defense response was not simply the product of more A-lineage honey bee ancestry as previously assumed, but rather an interaction between A-lineage and European alleles. Our results suggest that a combination of A-lineage and European alleles play roles in defensive behavior in AHBs.

Developmental plasticity shapes social traits and selection in a facultatively eusocial bee.

Developmental plasticity generates phenotypic variation, but how it contributes to evolutionary change is unclear. Phenotypes of individuals in caste-based (eusocial) societies are particularly sensitive to developmental processes, and the evolutionary origins of eusociality may be rooted in developmental plasticity of ancestral forms. We used an integrative genomics approach to evaluate the relationships among developmental plasticity, molecular evolution, and social behavior in a bee species (Megalopta genalis) that expresses flexible sociality, and thus provides a window into the factors that may have been important at the evolutionary origins of eusociality. We find that differences in social behavior are derived from genes that also regulate sex differentiation and metamorphosis. Positive selection on social traits is influenced by the function of these genes in development. We further identify evidence that social polyphenisms may become encoded in the genome via genetic changes in regulatory regions, specifically in transcription factor binding sites. Taken together, our results provide evidence that developmental plasticity provides the substrate for evolutionary novelty and shapes the selective landscape for molecular evolution in a major evolutionary innovation: Eusociality.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos.

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.

The Long Noncoding RNA Pnky Is a Trans-acting Regulator of Cortical Development In Vivo.

While it is now appreciated that certain long noncoding RNAs (lncRNAs) have important functions in cell biology, relatively few have been shown to regulate development in vivo, particularly with genetic strategies that establish cis versus trans mechanisms. Pnky is a nuclear-enriched lncRNA that is transcribed divergently from the neighboring proneural transcription factor Pou3f2. Here, we show that conditional deletion of Pnky from the developing cortex regulates the production of projection neurons from neural stem cells (NSCs) in a cell-autonomous manner, altering postnatal cortical lamination. Surprisingly, Pou3f2 expression is not disrupted by deletion of the entire Pnky gene. Moreover, expression of Pnky from a BAC transgene rescues the differential gene expression and increased neurogenesis of Pnky-knockout NSCs, as well as the developmental phenotypes of Pnky-deletion in vivo. Thus, despite being transcribed divergently from a key developmental transcription factor, the lncRNA Pnky regulates development in trans.

Integrative Genomics Reveals the Genetics and Evolution of the Honey Bee's Social Immune System.

Social organisms combat pathogens through individual innate immune responses or through social immunity-behaviors among individuals that limit pathogen transmission within groups. Although we have a relatively detailed understanding of the genetics and evolution of the innate immune system of animals, we know little about social immunity. Addressing this knowledge gap is crucial for understanding how life-history traits influence immunity, and identifying if trade-offs exist between innate and social immunity. Hygienic behavior in the Western honey bee, Apis mellifera, provides an excellent model for investigating the genetics and evolution of social immunity in animals. This heritable, colony-level behavior is performed by nurse bees when they detect and remove infected or dead brood from the colony. We sequenced 125 haploid genomes from two artificially selected highly hygienic populations and a baseline unselected population. Genomic contrasts allowed us to identify a minimum of 73 genes tentatively associated with hygienic behavior. Many genes were within previously discovered QTLs associated with hygienic behavior and were predictive of hygienic behavior within the unselected population. These genes were often involved in neuronal development and sensory perception in solitary insects. We found that genes associated with hygienic behavior have evidence of positive selection within honey bees (Apis), supporting the hypothesis that social immunity contributes to fitness. Our results indicate that genes influencing developmental neurobiology and behavior in solitary insects may have been co-opted to give rise to a novel and adaptive social immune phenotype in honey bees.

Draft Genome Assembly and Population Genetics of an Agricultural Pollinator, the Solitary Alkali Bee (Halictidae: Nomia melanderi).

Alkali bees (Nomia melanderi) are solitary relatives of the halictine bees, which have become an important model for the evolution of social behavior, but for which few solitary comparisons exist. These ground-nesting bees defend their developing offspring against pathogens and predators, and thus exhibit some of the key traits that preceded insect sociality. Alkali bees are also efficient native pollinators of alfalfa seed, which is a crop of major economic value in the United States. We sequenced, assembled, and annotated a high-quality draft genome of 299.6 Mbp for this species. Repetitive content makes up more than one-third of this genome, and previously uncharacterized transposable elements are the most abundant type of repetitive DNA. We predicted 10,847 protein coding genes, and identify 479 of these undergoing positive directional selection with the use of population genetic analysis based on low-coverage whole genome sequencing of 19 individuals. We found evidence of recent population bottlenecks, but no significant evidence of population structure. We also identify 45 genes enriched for protein translation and folding, transcriptional regulation, and triglyceride metabolism evolving slower in alkali bees compared to other halictid bees. These resources will be useful for future studies of bee comparative genomics and pollinator health research.

Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee.

Inbreeding (the mating between closely related individuals) often has detrimental effects that are associated with loss of heterozygosity at overdominant loci, and the expression of deleterious recessive alleles. However, determining which loci are detrimental when homozygous, and the extent of their phenotypic effects, remains poorly understood. Here, we utilize a unique inbred population of clonal (thelytokous) honey bees, Apis mellifera capensis, to determine which loci reduce individual fitness when homozygous. This asexual population arose from a single worker ancestor approximately 20 years ago and has persisted for at least 100 generations. Thelytokous parthenogenesis results in a 1/3 of loss of heterozygosity with each generation. Yet, this population retains heterozygosity throughout its genome due to selection against homozygotes. Deep sequencing of one bee from each of the three known sub-lineages of the population revealed that 3,766 of 10,884 genes (34%) have retained heterozygosity across all sub-lineages, suggesting that these genes have heterozygote advantage. The maintenance of heterozygosity in the same genes and genomic regions in all three sub-lineages suggests that nearly every chromosome carries genes that show sufficient heterozygote advantage to be selectively detrimental when homozygous.