The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics

BACKGROUND

The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí).

RESULTS

A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions.

CONCLUSIONS

Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

Transcriptional control of a metabolic switch regulating cellular methylation reactions is part of a common response to stress in divergent bee species

Recent global declines in bee health have elevated the need for a more complete understanding of the cellular stress mechanisms employed by diverse bee species. We recently uncovered the biomarker lethal (2) essential for life [l(2)efl] genes as part of a shared transcriptional program in response to a number of cell stressors in the western honey bee (Apis mellifera). Here, we describe another shared stress-responsive gene, glycine N-methyltransferase (Gnmt), which is known as a key metabolic switch controlling cellular methylation reactions. We observed Gnmt induction by both abiotic and biotic stressors. We also found increased levels of the GNMT reaction product sarcosine in the midgut after stress, linking metabolic changes with the observed changes in gene regulation. Prior to this study, Gnmt upregulation had not been associated with cellular stress responses in other organisms. To determine whether this novel stress-responsive gene would behave similarly in other bee species, we first characterized the cellular response to endoplasmic reticulum (ER) stress in lab-reared adults of the solitary alfalfa leafcutting bee (Megachile rotundata) and compared this with age-matched honey bees. The novel stress gene Gnmt was induced in addition to a number of canonical gene targets induced in both bee species upon unfolded protein response (UPR) activation, suggesting that stress-induced regulation of cellular methylation reactions is a common feature of bees. Therefore, this study suggests that the honey bee can serve as an important model for bee biology more broadly, although studies on diverse bee species will be required to fully understand global declines in bee populations.

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

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

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

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

The Bee Hemolymph Metabolome: A Window into the Impact of Viruses on Bumble Bees

State-of-the-art virus detection technology has advanced a lot, yet technology to evaluate the impacts of viruses on bee physiology and health is basically lacking. However, such technology is sorely needed to understand how multi-host viruses can impact the composition of the bee community. Here, we evaluated the potential of hemolymph metabolites as biomarkers to identify the viral infection status in bees. A metabolomics strategy based on ultra-high-performance liquid chromatography coupled to high-resolution mass spectrometry was implemented. First, we constructed a predictive model for standardized bumble bees, in which non-infected bees were metabolically differentiated from an overt Israeli acute paralysis virus (IAPV) infection (R²Y = 0.993; Q² = 0.906), as well as a covert slow bee paralysis virus (SBPV) infection (R²Y = 0.999; Q² = 0.875). Second, two sets of potential biomarkers were identified, being descriptors for the metabolomic changes in the bee's hemolymph following viral infection. Third, the biomarker sets were evaluated in a new dataset only containing wild bees and successfully discriminated virus infection versus non-virus infection with an AUC of 0.985. We concluded that screening hemolymph metabolite markers can underpin physiological changes linked to virus infection dynamics, opening promising avenues to identify, monitor, and predict the effects of virus infection in a bee community within a specific environment.

Using Manual and Computer-Based Text-Mining to Uncover Research Trends for Apis mellifera

Honey bee research is believed to be influenced dramatically by colony collapse disorder (CCD) and the sequenced genome release in 2006, but this assertion has never been tested. By employing text-mining approaches, research trends were tested by analyzing over 14,000 publications during the period of 1957 to 2017. Quantitatively, the data revealed an exponential growth until 2010 when the number of articles published per year ceased following the trend. Analysis of author-assigned keywords revealed that changes in keywords occurred roughly every decade with the most fundamental change in 1991-1992, instead of 2006. This change might be due to several factors including the research intensification on the Varroa mite. The genome release and CCD had quantitively only minor effects, mainly on honey bee health-related topics post-2006. Further analysis revealed that computational topic modeling can provide potentially hidden information and connections between some topics that might be ignored in author-assigned keywords.

Honey bee survival mechanisms against the parasite Varroa destructor: a systematic review of phenotypic and genomic research efforts

The ectoparasitic mite Varroa destructor is the most significant pathological threat to the western honey bee, Apis mellifera, leading to the death of most colonies if left untreated. An alternative approach to chemical treatments is to selectively enhance heritable honey bee traits of resistance or tolerance to the mite through breeding programs, or select for naturally surviving untreated colonies. We conducted a literature review of all studies documenting traits of A. mellifera populations either selectively bred or naturally selected for resistance and tolerance to mite parasitism. This allowed us to conduct an analysis of the diversity, distribution and importance of the traits in different honey bee populations that can survive V. destructor globally. In a second analysis, we investigated the genetic bases of these different phenotypes by comparing 'omics studies (genomics, transcriptomics, and proteomics) of A. mellifera resistance and tolerance to the parasite. Altogether, this review provides a detailed overview of the current state of the research projects and breeding efforts against the most devastating parasite of A. mellifera. By highlighting the most promising traits of Varroa-surviving bees and our current knowledge on their genetic bases, this work will help direct future research efforts and selection programs to control this pest. Additionally, by comparing the diverse populations of honey bees that exhibit those traits, this review highlights the consequences of anthropogenic and natural selection in the interactions between hosts and parasites.

The Progress of Multi-Omics Technologies: Determining Function in Lactic Acid Bacteria Using a Systems Level Approach

Lactic Acid Bacteria (LAB) have long been recognized as having a significant impact ranging from commercial to health domains. A vast amount of research has been carried out on these microbes, deciphering many of the pathways and components responsible for these desirable effects. However, a large proportion of this functional information has been derived from a reductionist approach working with pure culture strains. This provides limited insight into understanding the impact of LAB within intricate systems such as the gut microbiome or multi strain starter cultures. Whole genome sequencing of strains and shotgun metagenomics of entire systems are powerful techniques that are currently widely used to decipher function in microbes, but they also have their limitations. An available genome or metagenome can provide an image of what a strain or microbiome, respectively, is potentially capable of and the functions that they may carry out. A top-down, multi-omics approach has the power to resolve the functional potential of an ecosystem into an image of what is being expressed, translated and produced. With this image, it is possible to see the real functions that members of a system are performing and allow more accurate and impactful predictions of the effects of these microorganisms. This review will discuss how technological advances have the potential to increase the yield of information from genomics, transcriptomics, proteomics and metabolomics. The potential for integrated omics to resolve the role of LAB in complex systems will also be assessed. Finally, the current software approaches for managing these omics data sets will be discussed.

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

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

Identification of Proteins and Genes Expressed by Methylophaga thiooxydans During Growth on Dimethylsulfide and Their Presence in Other Members of the Genus

Dimethylsulfide is a volatile organic sulfur compound that provides the largest input of biogenic sulfur from the oceans to the atmosphere, and thence back to land, constituting an important link in the global sulfur cycle. Microorganisms degrading DMS affect fluxes of DMS in the environment, but the underlying metabolic pathways are still poorly understood. Methylophaga thiooxydans is a marine methylotrophic bacterium capable of growth on DMS as sole source of carbon and energy. Using proteomics and transcriptomics we identified genes expressed during growth on dimethylsulfide and methanol to refine our knowledge of the metabolic pathways that are involved in DMS and methanol degradation in this strain. Amongst the most highly expressed genes on DMS were the two methanethiol oxidases driving the oxidation of this reactive and toxic intermediate of DMS metabolism. Growth on DMS also increased expression of the enzymes of the tetrahydrofolate linked pathway of formaldehyde oxidation, in addition to the tetrahydromethanopterin linked pathway. Key enzymes of the inorganic sulfur oxidation pathway included flavocytochrome c sulfide dehydrogenase, sulfide quinone oxidoreductase, and persulfide dioxygenases. A sulP permease was also expressed during growth on DMS. Proteomics and transcriptomics also identified a number of highly expressed proteins and gene products whose function is currently not understood. As the identity of some enzymes of organic and inorganic sulfur metabolism previously detected in Methylophaga has not been characterized at the genetic level yet, highly expressed uncharacterized genes provide new targets for further biochemical and genetic analysis. A pan-genome analysis of six available Methylophaga genomes showed that only two of the six investigated strains, M. thiooxydans and M. sulfidovorans have the gene encoding methanethiol oxidase, suggesting that growth on methylated sulfur compounds of M. aminisulfidivorans is likely to involve different enzymes and metabolic intermediates. Hence, the pathways of DMS-utilization and subsequent C1 and sulfur oxidation are not conserved across Methylophaga isolates that degrade methylated sulfur compounds.

Genomes of the Hymenoptera

Hymenoptera is the second-most sequenced arthropod order, with 52 publically archived genomes (71 with ants, reviewed elsewhere), however these genomes do not capture the breadth of this very diverse order (Figure 1, Table 1). These sequenced genomes represent only 15 of the 97 extant families. Although at least 55 other genomes are in progress in an additional 11 families (see Table 2), stinging wasps represent 35 (67%) of the available and 42 (76%) of the in progress genomes. A more comprehensive catalog of hymenopteran genomes is needed for research into the evolutionary processes underlying the expansive diversity in terms of ecology, behavior, and physiological traits within this group. Additional sequencing is needed to generate an assembly for even 0.05% of the estimated 1 million hymenopteran species, and we recommend premier level assemblies for at least 0.1% of the >150,000 named species dispersed across the order. Given the haplodiploid sex determination in Hymenoptera, haploid male sequencing will help minimize genome assembly issues to enable higher quality genome assemblies.

Overview of Gene Expression Analysis: Transcriptomics

Currently, the study of the transcriptome is widely used to interpret the functional elements of the genome and molecular constituents of cells and tissues in an effort to unravel biological pathways associated with development and disease. The advent of technologies is now enabling the study of such comprehensive transcriptional characterization of mRNA, miRNA, lncRNA, and small RNA in a robust and successful manner. Transcriptomic strategies are gaining momentum across diverse areas of biological, plant sciences, medical, clinical, and pharmaceutical research for biomarker discovery, and disease diagnosis and prognosis.

Comparative proteomics and codon substitution analysis reveal mechanisms of differential resistance to hypoxia in congeneric snails

Although high-throughput proteomics has been widely applied to study mechanisms of environmental adaptation, the conclusions from studies that are based on one species can be confounded by phylogeny. We compare the freshwater snail Pomacea canaliculata (a notorious invasive species) and its congener Pomacea diffusa (a non-invasive species) to understand the molecular mechanisms of their differential resistance to hypoxia. A 72-h acute exposure experiment showed that P. canaliculata is more tolerant to hypoxia than P. diffusa. The two species were then exposed to three levels of dissolved oxygen (6.7, 2.0 and 1.0mgL^-1) for 8h, and their gill proteins were analyzed using iTRAQ-coupled LC-MS/MS. The two species showed striking differences in protein expression profiles, with the more hypoxia tolerant P. canaliculata having more up-regulated proteins in signal transduction and down-regulated proteins in glycolysis and the tricarboxylic acid cycle. Evolutionary analysis revealed five orthologous genes encoding differentially expressed proteins having clear signal of positive selection, indicating selection has acted on some of the hypoxia responsive genes. Our case study has highlighted the potential of integrated proteomics and comparative evolutionary analysis for understanding the genetic basis of adaptation to global environmental change in non-model species.

SIGNIFICANCE

Rapid globalization in recent decades has greatly facilitated species introduction around the world. Successfully established introduced species, so-called invasive species, have threatened the invaded ecosystems. There has been substantial interest in studying how invasive species respond to extreme environmental conditions because the results can help not only predict their range of expansion and manage their impact, but also may reveal the adaptive mechanisms underlying their invasiveness. Our study has adopted a comparative approach to study the differential physiological and proteomic responses of two congeneric snails to various hypoxic conditions, as well as codon substitution analysis at transcriptomic level to detect signals of positive selection in hypoxia-responsive genes. The integrated physiological, proteomic and transcriptomic approach can be applied in other non-model species to understand the molecular mechanisms of adaptation to global environmental change.

Genetics and developmental biology of cooperation

Despite essential progress towards understanding the evolution of cooperative behaviour, we still lack detailed knowledge about its underlying molecular mechanisms, genetic basis, evolutionary dynamics and ontogeny. An international workshop "Genetics and Development of Cooperation," organized by the University of Bern (Switzerland), aimed at discussing the current progress in this research field and suggesting avenues for future research. This review uses the major themes of the meeting as a springboard to synthesize the concepts of genetic and nongenetic inheritance of cooperation, and to review a quantitative genetic framework that allows for the inclusion of indirect genetic effects. Furthermore, we argue that including nongenetic inheritance, such as transgenerational epigenetic effects, parental effects, ecological and cultural inheritance, provides a more nuanced view of the evolution of cooperation. We summarize those genes and molecular pathways in a range of species that seem promising candidates for mechanisms underlying cooperative behaviours. Concerning the neurobiological substrate of cooperation, we suggest three cognitive skills necessary for the ability to cooperate: (i) event memory, (ii) synchrony with others and (iii) responsiveness to others. Taking a closer look at the developmental trajectories that lead to the expression of cooperative behaviours, we discuss the dichotomy between early morphological specialization in social insects and more flexible behavioural specialization in cooperatively breeding vertebrates. Finally, we provide recommendations for which biological systems and species may be particularly suitable, which specific traits and parameters should be measured, what type of approaches should be followed, and which methods should be employed in studies of cooperation to better understand how cooperation evolves and manifests in nature.