Nuclear translocation of vitellogenin in the honey bee (Apis mellifera)

Bixafen, Prothioconazole, and Trifloxystrobin Alone or in Combination Have a Greater Effect on Health Related Gene Expression in Honey Bees from Nutritionally Deprived than from Protein Supplemented Colonies

The aim of this study was to evaluate whether alterations in food availability compromise the metabolic homeostasis of honey bees exposed to three fungicides alone or together. Ten honey bee colonies were used, with half receiving carbohydrate-protein supplementation for 15 weeks while another five colonies had their protein supply reduced with pollen traps. Subsequently, forager bees were collected and exposed by contact to 1 or 7 µg of bixafen, prothioconazole, or trifloxystrobin, either individually or in combination. After 48 h, bee abdomens without the intestine were used for the analysis of expression of antioxidant genes (SOD-1, CAT, and GPX-1), detoxification genes (GST-1 and CYP306A1), the storage protein gene vitellogenin, and immune system antimicrobial peptide genes (defensin-1, abaecin, hymenoptaecin, and apidaecin), through real-time PCR. All fungicide treatments induced changes in gene expression, with bixafen showing the most prominent upregulation. Exposure to 1 µg of each of the three pesticides resulted in upregulation of genes associated with detoxification and nutrition processes, and downregulation of immune system genes. When the three pesticides were combined at a dose of 7 µg each, there was a pronounced downregulation of all genes. Food availability in the colonies affected the impact of fungicides on the expression of the studied genes in forager bees.

Transcriptional signature of host shift in the seed beetle Zabrotes subfasciatus

In phytophagous insects, adaptation to a new host is a dynamic process, in which early and later steps may be underpinned by different features of the insect genome. Here, we tested the hypothesis that early steps of this process are underpinned by a shift in gene expression patterns. We set up a short-term artificial selection experiment (10 generations) for the use of an alternative host (Cicer arietinum) on populations of the bean beetle Zabrotes subfasciatus. Using Illumina sequencing on young adult females, we show the selected populations differ in the expression of genes associated to stimuli, signalling, and developmental processes. Particularly, the "C. arietinum" population shows upregulation of histone methylation genes, which may constitute a strategy for fine-tuning the insect global gene expression network. Using qPCR on body regions, we demonstrated that the "Phaseolus vulgaris" population upregulates the genes polygalacturonase and egalitarian and that the expression of an odorant receptor transcript variant changes over generations. Moreover, in this population we detected the existence of vitellogenin (Vg) variants in both males and females, possibly harbouring canonical reproductive function in females and extracellular unknown functions in males. This study provides the basis for future genomic investigations seeking to shed light on the nature of the proximate mechanisms involved in promoting differential gene expression associated to insect development and adaptation to new hosts.

Identifying a developmental transition in honey bees using gene expression data

In many organisms, interactions among genes lead to multiple functional states, and changes to interactions can lead to transitions into new states. These transitions can be related to bifurcations (or critical points) in dynamical systems theory. Characterizing these collective transitions is a major challenge for systems biology. Here, we develop a statistical method for identifying bistability near a continuous transition directly from high-dimensional gene expression data. We apply the method to data from honey bees, where a known developmental transition occurs between bees performing tasks in the nest and leaving the nest to forage. Our method, which makes use of the expected shape of the distribution of gene expression levels near a transition, successfully identifies the emergence of bistability and links it to genes that are known to be involved in the behavioral transition. This proof of concept demonstrates that going beyond correlative analysis to infer the shape of gene expression distributions might be used more generally to identify collective transitions from gene expression data.

Genome-wide analysis of the chromatin sites targeted by HEX 70a storage protein in the honeybee brain and fat body

Hexamerins, the proteins massively stored in the larval haemolymph of insects, are gradually used throughout metamorphosis as a source of raw material and energy for the development of adult tissues. Such behaviour defined hexamerins as storage proteins. Immunofluorescence experiments coupled with confocal microscopy show a hexamerin, HEX 70a, in the nucleus of the brain and fat body cells from honeybee workers, an unexpected localization for a storage protein. HEX 70a colocalizes with fibrillarin, a nucleolar-specific protein and H3 histone, thus suggesting a potential role as a chromatin-binding protein. This was investigated through chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq). The significant HEX 70a-DNA binding sites were mainly localized at the intergenic, promoter and intronic regions. HEX 70a targeted DNA stretches mapped to the genomic regions encompassing genes with relevant functional attributes. Several HEX 70a targeted genes were associated with H3K27ac or/and H3K27me3, known as active and repressive histone marks. Brain and fat body tissues shared a fraction of the HEX 70 targeted genes, and tissue-specific targets were also detected. The presence of overrepresented DNA motifs in the binding sites is consistent with specific HEX 70a-chromatin association. In addition, a search for HEX 70a targets in RNA-seq public libraries of fat bodies from nurses and foragers revealed differentially expressed targets displaying hex 70a-correlated developmental expression, thus supporting a regulatory activity for HEX 70a. Our results support the premise that HEX 70a is a moonlighting protein that binds chromatin and has roles in the brain and fat body cell nuclei, apart from its canonical role as a storage protein.

Mixture toxicities of tetrachlorantraniliprole and tebuconazole to honey bees (Apis mellifera L.) and the potential mechanism

The extensive use of pesticides has negative effects on the health of insect pollinators. Although pollinators in the field are seldom exposed to individual pesticides, few reports have assessed the toxic impacts of pesticide combinations on them. In this work, we purposed to reveal the combined impacts of tetrachlorantraniliprole (TET) and tebuconazole (TEB) on honey bees (Apis mellifera L.). Our data exhibited that TET had greater toxicity to A. mellifera (96-h LC₅₀ value of 298.2 mg a.i. L^-1) than TEB (96-h LC₅₀ value of 1,841 mg a.i. L^-1). The mixture of TET and TEB displayed acute synergistic toxicity to the pollinators. Meanwhile, the activities of CarE, CYP450, trypsin, and sucrase, as well as the expressions of five genes (ppo, abaecin, cat, CYP4G11, and CYP6AS14) associated with immune response, oxidative stress, and detoxification metabolism, were conspicuously altered when exposed to the mixture relative to the individual exposures. These results provided an overall comprehension of honey bees upon the challenge of sublethal toxicity between neonicotinoid insecticides and triazole fungicides and could be used to assess the intricate toxic mechanisms in honey bees when exposed to pesticide mixtures. Additionally, these results might guide pesticide regulation strategies to enhance the honey bee populations.

Joint toxic effects of thiamethoxam and flusilazole on the adult worker honey bees (Apis mellifera L.)

Insect pollinators are routinely exposed to a complex mixture of many pesticides. However, traditional environmental risk assessment is only carried out based on ecotoxicological data of single substances. In this context, we aimed to explore the potential effects when worker honey bees (Apis mellifera L.) were simultaneously challenged by thiamethoxam (TMX) and flusilazole (FSZ). Results displayed that TMX possessed higher toxicity to A. mellifera (96-h LC₅₀ value of 0.11 mg a. i. L^-1) than FSZ (96-h LC₅₀ value of 738 mg a. i. L^-1). Furthermore, the mixture of TMX and FSZ exhibited an acute synergistic impact on the pollinators. Meanwhile, the activities of SOD, caspase 3, caspase 9, and PPO, as well as the expressions of six genes (abaecin, dorsal-2, defensin-2, vtg, caspase-1, and CYP6AS14) associated with oxidative stress, immune response, lifespan, cell apoptosis, and detoxification metabolism were noteworthily varied in the individual and mixture challenges than at the baseline level. These data revealed that it is imminently essential to investigate the combined toxicity of pesticides since the toxicity evaluation from individual compounds toward honey bees may underestimate the toxicity in realistic conditions. Overall, the present results could help understand the potential contribution of pesticide mixtures to the decline of bee populations.

Resolving the zinc binding capacity of honey bee vitellogenin and locating its putative binding sites

The protein vitellogenin (Vg) plays a central role in lipid transportation in most egg-laying animals. High Vg levels correlate with stress resistance and lifespan potential in honey bees (Apis mellifera). Vg is the primary circulating zinc-carrying protein in honey bees. Zinc is an essential metal ion in numerous biological processes, including the function and structure of many proteins. Measurements of Zn²⁺ suggest a variable number of ions per Vg molecule in different animal species, but the molecular implications of zinc-binding by this protein are not well-understood. We used inductively coupled plasma mass spectrometry to determine that, on average, each honey bee Vg molecule binds 3 Zn²⁺ -ions. Our full-length protein structure and sequence analysis revealed seven potential zinc-binding sites. These are located in the β-barrel and α-helical subdomains of the N-terminal domain, the lipid binding site, and the cysteine-rich C-terminal region of unknown function. Interestingly, two potential zinc-binding sites in the β-barrel can support a proposed role for this structure in DNA-binding. Overall, our findings suggest that honey bee Vg bind zinc at several functional regions, indicating that Zn²⁺ -ions are important for many of the activities of this protein. In addition to being potentially relevant for other egg-laying species, these insights provide a platform for studies of metal ions in bee health, which is of global interest due to recent declines in pollinator numbers.

Stabilized Double-Stranded RNA Strategy Improves Cotton Resistance to CBW (Anthonomus grandis)

Cotton is the most important crop for fiber production worldwide. However, the cotton boll weevil (CBW) is an insect pest that causes significant economic losses in infested areas. Current control methods are costly, inefficient, and environmentally hazardous. Herein, we generated transgenic cotton lines expressing double-stranded RNA (dsRNA) molecules to trigger RNA interference-mediated gene silencing in CBW. Thus, we targeted three essential genes coding for chitin synthase 2, vitellogenin, and ecdysis-triggering hormone receptor. The stability of expressed dsRNAs was improved by designing a structured RNA based on a viroid genome architecture. We transformed cotton embryos by inserting a promoter-driven expression cassette that overexpressed the dsRNA into flower buds. The transgenic cotton plants were characterized, and positive PCR transformed events were detected with an average heritability of 80%. Expression of dsRNAs was confirmed in floral buds by RT-qPCR, and the T1 cotton plant generation was challenged with fertilized CBW females. After 30 days, data showed high mortality (around 70%) in oviposited yolks. In adult insects fed on transgenic lines, chitin synthase II and vitellogenin showed reduced expression in larvae and adults, respectively. Developmental delays and abnormalities were also observed in these individuals. Our data remark on the potential of transgenic cotton based on a viroid-structured dsRNA to control CBW.

Identification of 121 variants of honey bee Vitellogenin protein sequences with structural differences at functional sites

Proteins are under selection to maintain central functions and to accommodate needs that arise in ever-changing environments. The positive selection and neutral drift that preserve functions result in a diversity of protein variants. The amount of diversity differs between proteins: multifunctional or disease-related proteins tend to have fewer variants than proteins involved in some aspects of immunity. Our work focuses on the extensively studied protein Vitellogenin (Vg), which in honey bees (Apis mellifera) is multifunctional and highly expressed and plays roles in immunity. Yet, almost nothing is known about the natural variation in the coding sequences of this protein or how amino acid-altering variants might impact structure-function relationships. Here, we map out allelic variation in honey bee Vg using biological samples from 15 countries. The successful barcoded amplicon Nanopore sequencing of 543 bees revealed 121 protein variants, indicating a high level of diversity in Vg. We find that the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) differs between protein regions with different functions; domains involved in DNA and protein-protein interactions contain fewer nsSNPs than the protein's lipid binding cavities. We outline how the central functions of the protein can be maintained in different variants and how the variation pattern may inform about selection from pathogens and nutrition.