Luana Bataglia
Zilá Luz Paulino Simões
Francis Morais Franco Nunes- 1Departamento de Genética, Faculdade de Medicina de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, Brazil
- 2Departamento de Biologia, Faculdade de Filosofia Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, Brazil
- 3Departamento de Genética e Evolução, Centro de Ciências Biológicas e da Saúde, Universidade Federal de São Carlos, São Carlos, Brazil
Honey bee (Apis mellifera) adult workers change behaviors and nutrition according to age progression. Young workers, such as nurses, perform in-hive tasks and consume protein-rich pollen, while older workers (foragers) leave the colony to search for food, and consume carbohydrate-rich nectar. These environmentally stimulated events involve transcriptional and DNA epigenetic marks alterations in worker tissues. However, post-transcriptional RNA modifications (epitranscriptomics) are still poorly explored in bees. We investigated the transcriptional profiles of m6A and m5C RNA methyltransferases in the brain and fat body of adult workers of 1) different ages and performing different tasks [nurses of 8 days-old (N-8D) and foragers of 29 days-old (F-29D), sampled from wild-type colonies], and 2) same-aged young workers caged in an incubator and treated with a pollen-rich [PR] or a pollen-deprived [PD] diet for 8 days. In the brain, METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 genes increased expression during adulthood (from N-8D to F-29D), while the opposite pattern was observed in the fat body for METTL3, DNMT2, and NSUN2 genes. Regarding diet treatments, high expression levels were observed in the brains of the pollen-deprived group (DNMT2, NOP2, and NSUN2 genes) and the fat bodies of the pollen-rich group (NOP2, NSUN4, and NSUN5 genes) compared to the brains of the PR group and the fat bodies of the PD group, respectively. Our data indicate that RNA epigenetics may be an important regulatory layer in the development of adult workers, presenting tissue-specific signatures of RNA methyltransferases expression in response to age, behavior, and diet content.
Introduction
The adult life of honey bee (Apis mellifera) workers is characterized by a temporal polyethism, in which individuals progressively experience different tasks and behaviors during their lifetime. This age-related polyethism is accompanied by alteration of the feeding regimes. Young workers feed on pollen, a predominant source of protein, lipids, vitamins, and minerals (Brodschneider and Crailsheim, 2010). During this period, they perform in-hive tasks, such as feeding the developing larvae, a behavior that characterizes them as nurses (Winston, 1987; Page and Peng, 2001). Around the third week of adulthood, workers start to go outside the colony to collect pollen, nectar, resins, and water, a behavior that characterizes them as foragers (Winston, 1987; Page and Peng, 2001). The foragers consume nectar, a diet rich in carbohydrates and low in proteins, which increases the bee’s energy metabolism, necessary for flight activity and resource gathering (Winston, 1987). This behavioral transition in honey bees involves transcriptional and metabolic changes in tissues, especially in the brain and the fat body (Ament et al., 2011). Foragers have an increase in cognitive capacity (memory and learning) compared to young workers performing in-hive tasks (Zayed and Robinson, 2012), and their brain transcriptomes differ significantly (Whitfield et al., 2003). The fat body is a metabolic center and during the behavior transition, the foragers lose half of their lipid stores and the capacity to metabolize these molecules (Toth and Robinson, 2005).
Epigenetics is proposed to be a plastic mechanism for regulating gene expression in response to environmental conditions (Vaiserman, 2014; Abdul et al., 2017). Diet is one of the external factors that alters gene expression through epigenetic modifications of the genomic DNA (epigenomics). The nutrients and bioactive components present in food can alter the activity of genes involved in DNA methylation or histone modifications or can change the availability of substrates required for these modifications’ reactions (reviewed by Abdul et al., 2017). In honey bees, the different diets consumed by workers and queens’ larvae trigger distinct abundance and patterns of DNA methylation between castes (Kucharski et al., 2008). In workers, several DNA CpG sites are differentially methylated comparing the brains of nurses and foragers (Herb et al., 2012; Lockett et al., 2012). Although knowledge in bee epigenomics has advanced in recent years, little has been explored about how epitranscriptomics, an emerging field of investigation, acts in bee development.
Epitranscriptomics is the study of RNA post-transcriptional modifications. More than 170 chemical modifications present in RNAs are currently known (Boccaletto et al., 2017; Huang et al., 2020). Two of them, the methylation of the nitrogen N6 of adenosine (N6-methyladenosine, m6A) and the methylation of the carbon C5 of cytosine (5-methylcytosine, m5C) are the most studied and regulate several biological processes such as development in mouse (Tuorto et al., 2012) and Drosophila melanogaster (Lence et al., 2016), behavioral adaptation in mammals (Meyer et al., 2012), aging and response to environmental (feeding, temperature) stimuli in D. melanogaster, Caenorhabditis elegans and Saccharomyces cerevisiae (Schosserer et al., 2015). Interestingly, these biological processes are closely related to those observed in the context of age polyethism such that RNA modifications have become excellent candidates to be explored in honey bees.
To date, over 1,000 RNA epigenetic studies were published, however, only 3% of them investigated insects (source: PubMed, June 2022). Most of them used D. melanogaster, with RNA methylation being related to sex determination (Haussmann, et al., 2016; Kan et al., 2017), fertility (Hongay and Orr-Weaver, 2011), and nervous system development (Lence et al., 2016). In Bombyx mori, m6A was reported to control gene expression and cell cycle progression (Li et al., 2019). In honey bees, two studies revealed RNA modifications acting in caste differentiation (Bataglia et al., 2021; Wang et al., 2021). Our group recently identified the orthologs of RNA methyltransferase genes in the A. mellifera genome, being METTL3 (methyltransferase like 3) and METTL14 (methyltransferase like 14) for m6A methylation, and DNMT2 (DNA methyltransferase 2) and the NSUN (NOP2/Sun RNA methyltransferase member) family encoding genes (NOP2, NSUN2, NSUN4, NSUN5, and NSUN7) for m5C methylation (Bataglia et al., 2021).
Here, we verified if the expression of m6A and m5C RNA methyltransferase genes vary in the development of A. mellifera adult workers according to age-related behavior (nurses of 8 days-old versus foragers of 29 days-old), diet content (pollen-rich or pollen-deprived), or tissue (brain or fat body). We found that the expression of some RNA methyltransferase transcripts is modulated during temporal polyethism in a tissue-specific and nutrition-dependent manner.
Materials and methods
Samples
Apis mellifera workers were collected from three queenright colonies maintained at the experimental apiary of the Universidade de São Paulo in Ribeirão Preto, Brazil. Three independent biological experiments were performed, each one starting with 1 day of difference, on three consecutive days, each day with bees from a different colony. For each experiment, combs containing pharate adults close to emergence (Pbd stage) were placed in an incubator at 34°C and ∼80% relative humidity for 8 h to obtain at least 400 newly-emerged adult workers.
Each experiment was performed according to the following design. A total of 100 newly-emerged workers were divided (day 1) into two cages (n = 50 per cage of 20 cm × 13 cm × 13 cm) and kept in the incubator, with food and water ad libitum for 8 days. Every day, the food and the water were replaced, and the dead bees were removed. The bees of one cage received a pollen-deprived (PD) diet (a mix of 20% honey plus 80% powdered sugar, named candy), and the bees of the other cage received a pollen-rich (PR) diet (a mix of 70% candy and 30% of fresh poly floral pollen collected from the same apiary). On day 8, five individuals per cage were collected. The newly-emerged workers that were not placed in cages (at least 300 workers), were then marked on the thorax with uni POSCA non-toxic pens (1.8–2.5 mm), for age control, and returned to their original colony. We collected five nurses of 8 days-old visiting larval cells and five foragers of 29 days-old returning to the colony (N-8D and F-29D, respectively).
The collected workers were anesthetized at 4°C for 5 min for tissue dissection. Then, we quickly dissected the brain and the fat body of each worker. Brains were dissected using a Leica M125 stereomicroscope and all glands and adjacent tissues were removed. For sampling the fat body, we pulled the sting along with the last abdominal segment, to remove the venom gland, the intestine, and adjacent tissues (such as the ovaries), leaving only the abdominal carcass. The abdominal carcass is internally lined by the fat body (Ament et al., 2011). Each individual brain or abdominal fat body was placed separately in tubes containing TRIzol® (Invitrogen) for RNA extraction.
Total RNA extraction, cDNA synthesis, and real-time PCR
The total RNA was extracted from individual brains and fat bodies according to Trizol® protocol. To avoid genomic DNA contamination, the total RNA was treated with DNase (DNase I Amplification Grade—Invitrogen). Next, cDNA was synthesized by reverse transcription catalyzed by SuperScript II reverse transcriptase (200 U/μL, Invitrogen) and oligo (dT)12–18 primers.
Real-time PCR was performed with 60 brain samples and 60 fat body samples. Among these 60 samples of each tissue, 15 samples represent each condition (PD, PR, N-8D, F-29D), being five samples of each experiment. The expression of the following methyltransferase genes was evaluated: for m6A, METTL3 (NCBI Gene ID: 551911) and METTL14 (409900), and for m5C, DNMT2 (410512), NOP2 (726212), NSUN2 (411579), NSUN4 (102655259), NSUN5 (410262), and NSUN7 (413907).
We perform the real-time PCR following the PowerUP™ SYBR™ Green Master Mix (2X) protocol and cycling (Applied Biosystems). For each sample, reactions were assembled in triplicate, and every single reaction consisted of a final volume of 20 μl, containing: 10 μl of PowerUP™ SYBR™ Green Master Mix (2X) solution (Applied Biosystems), 6.4 μl of water, 0.8 μl of each primer, forward and reverse (10 pmol/μL), and 2 μl of the cDNA solution (10 ng/μl). The primer pairs used are listed in Supplementary Table S1, including the pair used for amplification of the ribosomal protein L32 (RpL32) encoding gene, used as a reference for data normalization (Grozinger et al., 2003; Lourenço et al., 2008). The amplification reactions were conducted on a 7,500 Real-Time PCR System (Applied Biosystems). For primer pairs with optimal efficiency at 60°C, the following cycling conditions were used: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. For primer pairs with annealing temperature lower than 60°C, the following cycle was used: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 15 s, primer annealing temperature for 15 s, and 72°C for 1 min.
The relative expressions were calculated using the 2-∆∆CT method (Livak and Schmittgen, 2001). The statistical analyses were performed in R (version 4.1.3). We applied the Generalized Linear Mixed (GLM) models (package lme4 version 1.1–29) followed by the Anova function (package car version 3.0.12) to compare the different diets or ages with gene expression, using colony or cage replicate as a random effect. Significant outliers were removed [Grubbs’ test (alpha = 0.05)] when detected.
Results
The expression of METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 methyltransferases statistically vary in the brain of workers of different ages, performing different tasks (N-8D versus F-29D comparison). Also, DNMT2, NOP2, and NSUN2 were differentially expressed in the brain of same-aged young workers that consumed distinct diets (PR versus PD comparison). In the fat body, the expression of METTL3, DNMT2, and NSUN2 genes were influenced by age/task (N-8D versus F-29D comparison), while the expression of NOP2, NSUN4, and NSUN5 was influenced by diet (PR versus PD comparison). All these results are graphically represented in Figure 1 (and the results of the statistical analyses are in Supplementary Table S2).
FIGURE 1. Relative expression of m6A (METTL3 and METTL14) and m5C (DNMT2, NOP2, NSUN2, NSUN4, NSUN5, and NSUN7) RNA methyltransferase genes, measured by real-time PCR in the brain and fat body samples from workers of different ages/behaviors (N-8D for 8 days-old nurses–blue—and F-29D for 29 days-old foragers–orange) and under different dietary regimes (PR for pollen-rich diet–green—and PD for pollen-deprived–purple), n = 15 per condition. The boxplots represent the minimum, maximum, and median values. The presence of asterisks indicates a significant difference in expression, *p < 0.05, *p < 0.01, **p < 0.001.
A clear expression pattern for methyltransferase genes should be highlighted. Regarding the gene set that presented significant transcriptional differences in the brain, we observed lower levels in both N-8D and PR groups when compared to F-29D and PD groups. In addition, an opposite pattern was observed in the fat body, with increased expression in N-8D and PR groups in relation to F-29D and PD groups.
Discussion
The adult life progression of A. mellifera workers is marked by changes ranging from behavior to nutrition, mainly in the brain and the fat body (Ament et al., 2011). Despite the age-related division of labor in workers (temporal polyethism) being well-documented (Page et al., 2006), the participation of RNA modifications in these processes remains largely unknown (Bataglia et al., 2021). To advance this knowledge, we have investigated the transcriptional expression of m6A and m5C RNA methyltransferases in the brain and fat body of workers of different ages, performing different tasks, or same-aged young workers that consumed a pollen-deprived or a pollen-rich diet.
Foraging tasks of workers require well-developed brain abilities such as, recognizing the environment, identifying food resources, and communicating this information with nestmates. During the transition from in-hive to outside tasks, structural modifications such as the increase of the mushroom bodies (Withers et al., 1993; Durst et al., 1994) and the neuronal dendritic arborization (Farris et al., 2001) are observed. Higher expression levels of METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 transcripts were found in the brain of foragers compared to nurses. All these genes are reported in the literature as being involved in brain development, memory formation, or increase of cognitive capacity. The overexpression of METTL3 enhances hippocampus-dependent long-term memory in mice (Zhang et al., 2018) while its depletion perturbs short-term memory formation in flies (Kan et al., 2021). High abundance of mRNAs of DNMT2 (Rai et al., 2007) and NOP2 (Kosi et al., 2015) are related to neurogenesis in zebrafish brains and neural stem cell proliferation in mammals, respectively. The lack of NSUN2 expression inhibits normal brain development by reducing the number of differentiated neurons in mice (Flores et al., 2017) and short-term memory deficits in Drosophila (Abbasi-Moheb et al., 2012). NSUN5 mutant mice showed deficits in spatial cognitive abilities (Zhang et al., 2019) and in cerebral cortex development (Chen et al., 2019). Transcriptomic evidence suggests that NSUN7 may act on the cognitive processes of mammals (Chi and Delgado-Olguin. 2013; Tang et al., 2019). It is increasingly evident that the development of the nervous system, cognitive and behavioral aspects, and learning and memory abilities, in different organisms, are associated with a reprogramming (spatial and temporal) of the epitranscriptome (Jung and Goldman, 2018; Leighton et al., 2018). Our study supports and extends such evidence.
The main function of fat body of insects is the metabolism of macromolecules (reviewed by Skowronek et al., 2021). RNA methylation acts on nutritional metabolism, mainly on glucose and lipid metabolism (Wu et al., 2020). For example, METTL3 knockdown decreased lipid accumulation in human liver cell culture (Zhong et al., 2018). The difference in the expression of RNA methyltransferases in the fat body of nurses and foragers found here may be integrated into the gene regulatory network associated with the conversion of a protein-lipid metabolism in younger workers to a carbohydrate metabolism in older ones (Ament et al., 2011). Our data showed that not only age, but diet also influenced the expression of some methyltransferases. NOP2, NSUN4, and NSUN5 were differentially expressed in the fat body of same-aged workers that fed on different diets under caged conditions and did not perform behavioral tasks. Thus, these genes are excellent candidates for further studies about nutritional metabolism uncoupled to social interactions in honey bees. Moreover, in the worker brains, DNMT2, NOP2, and NSUN2 expression levels were diet-responsive. Our data are further corroborated by a study performed with rats fed on high or low-calorie diets, that showed the nutritional influence on the differential expression of RNA methyltransferases depending on the tissue analyzed (D’Aquila et al., 2022). It is important to highlight that the pollen-rich (PR) diet mimics that consumed by nurses while the pollen-deprived (PD) diet is similar to the nutritional content of foragers. Recently, Martelli et al. (2022) showed that young honey bee workers fed on PD diet age precociously and present biological parameters similar to foragers. Here, we observed corresponding expression profiles of some m5C genes (DNMT2, NOP2, and NSUN2) between N-8D + PR groups, which contrast with data from F-29D + PD groups. These genes will be tested in the future as potential candidates for aging markers in adult worker brains.
The fat body is also important for the production of energy and for the synthesis of immune system proteins (reviewed by Skowronek et al., 2021). Changes in the activity of methyltransferases have already been observed in response to stressful conditions. METTL3 plays a crucial role in oxidative stress, glycolytic stress, and DNA damage, and promotes an increased lipid metabolism in different mammalian cell lines (reviewed by Wilkinson et al., 2021). DNMT2 integrates stress response pathways by protecting tRNAs against ribonuclease cleavage induced by heat shock conditions in Drosophila (Schaefer et al., 2010), while NOP2 is involved with resistance to heat stress in C. elegans (Heissenberger et al., 2020). In mice, NSUN2 is required in cellular adaptation in response to toxic compounds and UV radiation (Gkatza et al., 2019). In honey bees, nurses have a more robust immune response to stress conditions compared to foragers (Amdam et al., 2004; Lourenço et al., 2019). In addition, our findings on the higher expression of some methyltransferases in the fat body of young bees compared to old ones seem somewhat coherent with this immunological scenario. Further research is needed to check if this suggestive crosstalk between immunity and epitranscriptomics is also true for bees.
This scientific topic is promising and some studies already involve epigenetic mechanisms to the xenobiotic response in insects. Epigenomic marks are associated with modulation of the expression of pesticide resistant genes in the beetle Meligethes aeneus (Erban et al., 2017). Also, strong evidence suggests that the expression of METTL3, METTL14 and a set of detoxification P450 genes along with increased levels of m6A marks confer thiamethoxam resistance to a strain of whitefly crop pest, Bemisia tabaci (Yang et al., 2021). It is important to consider that pollen is a common source of pesticides in agricultural sites and the contamination can be the cause of poisoning incidents of honey bees (Kiljanek et al., 2016). Although we did not screen the pollen for pesticides, our apiary is located within a non-agricultural area of a University with six million m2 containing an abundance of flowering plants representing nearly 300 species, and no poisoning effects were observed on bees we used as a source for the experiments. Furthermore, it is known that pollen consumption reduces the toxicity of pesticides in Apis mellifera (Barascou et al., 2021) and that diets influence the expression of genes of the epitranscriptomic machinery (this study). Perspectives are open to investigate whether RNA modifications respond to the detrimental impacts of insecticides and/or whether they confer some resistance to these pollinators.
Here we observed co-expression of RNA methyltransferase genes in all tested contexts in worker bees irrespective of their nutrition or age. Does this co-expression indicate functional redundancy? Do RNA methyltransferases have specific or integrated actions? It is widely accepted that all RNA methyltransferases have the canonical function of methylating RNAs, and RNA modifications are found in viruses (Kennedy et al., 2017), prokaryotes (Deng et al., 2015), yeasts (Clancy et al., 2002), plants (Fray & Simpson 2015), invertebrates (Hongay & Orr-Weaver, 2011), and vertebrates (D’Aquila et al., 2022). Despite METTL3 gene silencing in human embryonic stem cells is not sufficient for a complete reduction of m6A methylation levels, this knockdown at least removes the methylation marks of a specific subset of pluripotency-related genes, influencing the process of cell differentiation (Batista et al., 2014). This suggests a co-working, but it is also true that an individual RNA methyltransferase may have preferential RNA species as targets. NSUN2 preferentially methylates mRNA, tRNAs, and mitochondrial tRNAs (Tuorto et al., 2012; Shinoda et al., 2019), whereas NSUN5 preferentially methylates rRNA (Schosserer et al., 2015). Also, m6A sites are found mainly in mRNAs and affect the physical interaction of RNA with proteins or regulatory factors, resulting in altered translation or changes in mRNA stability (Meyer et al., 2015). This methylation also participates in alternative splicing (Zhao et al., 2014) and miRNA processing (Alarcón et al., 2015) events. The m5C marks are more frequent in non-protein-coding RNAs (ncRNAs, mainly tRNAs, and rRNAs) and maintain the stability of their secondary structures for the correct translation events (Micura et al., 2001; Schaefer et al., 2010). We believe that new findings should emerge in the coming years to clarify both species-specific or conserved functional issues related to epitranscriptomics. This is the case of a recent article that studied the m5C methylation patterns of maternal mRNAs in six species (five vertebrates: zebrafish, Xenopus laevis, Xenopus tropicalis, mouse, human, and one invertebrate: D. melanogaster) spanning ∼800 million years of evolution. For all tested species, they found a massive methylation of maternal mRNAs occurring in early embryonic development, and the rates drop dramatically after the mater-zygotic transition. This data shows an extremely conserved pattern in animals. On the other hand, they also observed regulatory innovations, such as the gain of m5C sites in mRNAs during evolution (Liu et al., 2022).
Concerning the differences we found in the expression of m6A and m5C methyltransferases between honey bee tissues, we speculate that they may reflect in the regulation of tissue-specific splicing variants previously found in the brain and fat body (Kannan et al., 2019), and/or several differentially expressed coding-genes and miRNA from workers performing different behaviors (Whitfield et al., 2003; Behura and Whitfield, 2010).
Bataglia et al. (2021) reported that global levels of m6A and m5C methylation quantified in the total RNA from brain and fat body samples increases in age progression. In this study we observed an increase in the expression of RNA methyltransferases in brain which is in accordance with the overall increase in global m6A and m5C methylation levels (Bataglia et al., 2021). However, in the fat body, although the higher expression of RNA methyltransferases is found in younger bees, the global m6A and m5C methylation levels are lower in young bees when compared with older bees. This point may sound contradictory, but it is important to emphasize that epitranscriptomic machinery is not only composed of RNA methyltransferases (named as writers) but also of readers and demethylases (erasers), and these last two sets of enzymes act in the recognition and removal of methyl groups, respectively (Schaefer et al., 2017). So, when the global methylation levels are accessed, we are measuring the resulting action of the whole machinery (writers, readers, and erasers) on the analyzed transcriptome.
We conclude that the adult life progression of A. mellifera workers is accompanied by changes in the expression of m6A and m5C RNA methyltransferases. It is possible that the mRNA methylome is dynamically regulated all through the lifespan of workers. We found that brain and fat body influence the transcriptional levels of these methyltransferases in opposite ways. Further the transcript levels of some of these genes are impacted by the protein content of food consumed as well as by the age-related behavior performed. Our results point to a tissue-specific transcriptional signature of RNA methyltransferases during the progression of adulthood.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
FMFN conceived and supervised the study. LB performed sampling, experiments, data collection, and analysis. ZLPS provided materials and reagents. FMFN and ZLPS acquired funding. FMFN and LB wrote the manuscript. All authors performed data interpretation, reviewed, and edited the manuscript, and approved its final version.
Funding
This research was financially supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, universal grant, process number 461711/2014-1), the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, processes number 16/06657-0 and 19/02374-1), the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES, Finance Code 001), and a fellowship for LB (CNPq, process number 130285/2016-1).
Acknowledgments
We are very grateful to Marcela Laure, Luís Aguiar, and Vera Figueiredo for their valuable technical assistance, to Kate Ihle, Felipe Martelli and Carlos Cardoso Jr for statistical advice, and to Navdeep Mutti for critical reading.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.921503/full#supplementary-material
References
Abbasi-Moheb, L., Mertel, S., Gonsior, M., Nouri-Vahid, L., Kahrizi, K., Cirak, S., et al. (2012). Mutations in NSUN2 cause autosomal-recessive intellectual disability. Am. J. Hum. Genet. 90 (5), 847–855. doi:10.1016/j.ajhg.2012.03.021
Abdul, Q. A., Yu, B. P., Chung, H. Y., Jung, H. A., and Choi, J. S. (2017). Epigenetic modifications of gene expression by lifestyle and environment. Arch. Pharm. Res. 40 (11), 1219–1237. doi:10.1007/s12272-017-0973-3
Alarcón, C. R., Lee, H., Goodarzi, H., Halberg, N., and Tavazoie, S. F. (2015). N6-methyladenosine marks primary microRNAs for processing. Nature 519 (7544), 482–485. doi:10.1038/nature14281
Amdam, G. V., Simões, Z. L., Hagen, A., Norberg, K., Schrøder, K., Mikkelsen, Ø., et al. (2004). Hormonal control of the yolk precursor vitellogenin regulates immune function and longevity in honeybees. Exp. Gerontol. 39 (5), 767–773. doi:10.1016/j.exger.2004.02.010
Ament, S. A., Chan, Q. W., Wheeler, M. M., Nixon, S. E., Johnson, S. P., Rodriguez-Zas, S. L., et al. (2011). Mechanisms of stable lipid loss in a social insect. J. Exp. Biol. 214 (22), 3808–3821. doi:10.1242/jeb.060244
Barascou, L., Sene, D., Barraud, A., Michez, D., Lefebvre, V., Medrzycki, P., et al. (2021). Pollen nutrition fosters honeybee tolerance to pesticides. R. Soc. Open Sci. 18 (9), 210818. doi:10.1098/rsos.210818
Bataglia, L., Simões, Z. L. P., and Nunes, F. M. F. (2021). Active genic machinery for epigenetic RNA modifications in bees. Insect Mol. Biol. 30 (6), 566–579. doi:10.1111/imb.12726
Batista, P. J., Molinie, B., Wang, J., Qu, K., Zhang, J., Li, L., et al. (2014). m(6)A RNA modification controls cell fate transition in mammalian embryonic stem cells. Cell Stem Cell 15 (6), 707–719. doi:10.1016/j.stem.2014.09.019
Behura, S. K., and Whitfield, C. W. (2010). Correlated expression patterns of microRNA genes with age-dependent behavioural changes in honeybee. Insect Mol. Biol. 19 (4), 431–439. doi:10.1111/j.1365-2583.2010.01010.x
Boccaletto, P., Machnicka, M. A., Purta, E., Piatkowski, P., Baginski, B., Wirecki, T. K., et al. (2017). Modomics: A database of RNA modification pathways. 2017 update. Nucleic Acids Res. 46 (1), D303–D307. doi:10.1093/nar/gkx1030
Brodschneider, R., and Crailsheim, K. (2010). Nutrition and health in honey bees. Apidologie 41 (3), 278–294. doi:10.1051/apido/2010012
Chen, P., Zhang, T., Yuan, Z., Shen, B., and Chen, L. (2019). Expression of the RNA methyltransferase Nsun5 is essential for developing cerebral cortex. Mol. Brain 12 (1), 74. doi:10.1186/s13041-019-0496-6
Chi, L., and Delgado-Olguín, P. (2013). Expression of NOL1/NOP2/sun domain (Nsun) RNA methyltransferase family genes in early mouse embryogenesis. Gene Expr. Patterns. 13 (8), 319–327. doi:10.1016/j.gep.2013.06.003
Clancy, M. J., Shambaugh, M. E., Timpte, C. S., and Bokar, J. A. (2002). Induction of sporulation in Saccharomyces cerevisiae leads to the formation of N6-methyladenosine in mRNA: A potential mechanism for the activity of the IME4 gene. Nucleic Acids Res. 30 (20), 4509–4518. doi:10.1093/nar/gkf573
D'Aquila, P., De Rango, F., Paparazzo, E., Mandalà, M., Bellizzi, D., and Passarino, G. (2022). Impact of nutrition on age-related epigenetic RNA modifications in rats. Nutrients 14 (6), 1232. doi:10.3390/nu14061232
Deng, X., Chen, K., Luo, G. Z., Weng, X., Ji, Q., Zhou, T., et al. (2015). Widespread occurrence of N6-methyladenosine in bacterial mRNA. Nucleic Acids Res. 43 (13), 6557–6567. doi:10.1093/nar/gkv596
Durst, C., Eichmüller, S., and Menzel, R. (1994). Development and experience lead to increased volume of subcompartments of the honeybee mushroom body. Behav. Neural Biol. 62 (3), 259–263. doi:10.1016/s0163-1047(05)80025-1
Erban, T., Harant, K., Chalupnikova, J., Kocourek, F., and Stara, J. (2017). Beyond the survival and death of the deltamethrin-threatened pollen beetle Meligethes aeneus: An in-depth proteomic study employing a transcriptome database. J. Proteomics 50, 281–289. doi:10.1016/j.jprot.2016.09.016
Farris, S. M., Robinson, G. E., and Fahrbach, S. E. (2001). Experience- and age-related outgrowth of intrinsic neurons in the mushroom bodies of the adult worker honeybee. J. Neurosci. 21 (16), 6395–6404. doi:10.1523/JNEUROSCI.21-16-06395.2001
Flores, J. V., Cordero-Espinoza, L., Oeztuerk-Winder, F., Andersson-Rolf, A., Selmi, T., Blanco, S., et al. (2017). Cytosine-5 RNA methylation regulates neural stem cell differentiation and motility. Stem Cell Rep. 8 (1), 112–124. doi:10.1016/j.stemcr.2016.11.014
Fray, R. G., and Simpson, G. G. (2015). The Arabidopsis epitranscriptome. Curr. Opin. Plant Biol. 27, 17–21. doi:10.1016/j.pbi.2015.05.015
Gkatza, N. A., Castro, C., Harvey, R. F., Heiß, M., Popis, M. C., Blanco, S., et al. (2019). Cytosine-5 RNA methylation links protein synthesis to cell metabolism. PLoS Biol. 17 (6), e3000297. doi:10.1371/journal.pbio.3000297
Grozinger, C. M., Sharabash, N. M., Whitfield, C. W., and Robinson, G. E. (2003). Pheromone-mediated gene expression in the honey bee brain. Proc. Natl. Acad. Sci. U. S. A. 100 (2), 14519–14525. doi:10.1073/pnas.2335884100
Haussmann, I. U., Bodi, Z., Sanchez-Moran, E., Mongan, N. P., Archer, N., Fray, R. G., et al. (2016). m6A potentiates Sxl alternative pre-mRNA splicing for robust Drosophila sex determination. Nature 540 (7632), 301–304. doi:10.1038/nature20577
Heissenberger, C., Rollins, J. A., Krammer, T. L., Nagelreiter, F., Stocker, I., Wacheul, L., et al. (2020). The ribosomal RNA m5C methyltransferase NSUN-1 modulates healthspan and oogenesis in Caenorhabditis elegans. Elife 9, e56205. doi:10.7554/eLife.56205
Herb, B. R., Wolschin, F., Hansen, K. D., Aryee, M. J., Langmead, B., Irizarry, R., et al. (2012). Reversible switching between epigenetic states in honeybee behavioral subcastes. Nat. Neurosci. 15 (10), 1371–1373. doi:10.1038/nn.3218
Hongay, C. F., and Orr-Weaver, T. L. (2011). Drosophila Inducer of MEiosis 4 (IME4) is required for Notch signaling during oogenesis. Proc. Natl. Acad. Sci. U. S. A. 108 (36), 14855–14860. doi:10.1073/pnas.1111577108
Huang, H., Weng, H., and Chen, J. (2020). The biogenesis and precise control of RNA m6A methylation. Trends Genet. 36 (1), 44–52. doi:10.1016/j.tig.2019.10.011
Jung, Y., and Goldman, D. (2018). Role of RNA modifications in brain and behavior. Genes Brain Behav. 17 (3), e12444. doi:10.1111/gbb.12444
Kan, L., Grozhik, A. V., Vedanayagam, J., Patil, D. P., Pang, N., Lim, K. S., et al. (2017). The m6A pathway facilitates sex determination in Drosophila. Nat. Commun. 8, 15737. doi:10.1038/ncomms15737
Kan, L., Ott, S., Joseph, B., Park, E. S., Dai, W., Kleiner, R. E., et al. (2021). A neural m6A/Ythdf pathway is required for learning and memory in Drosophila. Nat. Commun. 12 (1), 1458. doi:10.1038/s41467-021-21537-1
Kannan, K., Shook, M., Li, Y., Robinson, G. E., and Ma, J. (2019). Comparative analysis of brain and fat body gene splicing patterns in the honey bee, Apis mellifera. G3 (Bethesda) 9 (4), 1055–1063. doi:10.1534/g3.118.200857
Kennedy, E. M., Courtney, D. G., Tsai, K., and Cullen, B. R. (2017). Viral epitranscriptomics. J. Virol. 91 (9), e02263–16. doi:10.1128/JVI.02263-16
Kiljanek, T., Niewiadowska, A., and Posyniak, A. (2016). Pesticide poisoning of honeybees: A review of symptoms, incident classification, and causes of poisoning. J. Apic. Sci. 60 (2), 5–24. doi:10.1515/jas-2016-0024
Kosi, N., Alić, I., Kolačević, M., Vrsaljko, N., Jovanov Milošević, N., Sobol, M., et al. (2015). Nop2 is expressed during proliferation of neural stem cells and in adult mouse and human brain. Brain Res. 1597, 65–76. doi:10.1016/j.brainres.2014.11.040
Kucharski, R., Maleszka, J., Foret, S., and Maleszka, R. (2008). Nutritional control of reproductive status in honeybees via DNA methylation. Science 319 (5871), 1827–1830. doi:10.1126/science.1153069
Leighton, L. J., Ke, K., Zajaczkowski, E. L., Edmunds, J., Spitale, R. C., and Bredy, T. W. (2018). Experience-dependent neural plasticity, learning, and memory in the era of epitranscriptomics. Genes Brain Behav. 17 (3), e12426. doi:10.1111/gbb.12426
Lence, T., Akhtar, J., Bayer, M., Schmid, K., Spindler, L., Ho, C. H., et al. (2016). m6A modulates neuronal functions and sex determination in Drosophila. Nature 540 (7632), 242–247. doi:10.1038/nature20568
Li, B., Wang, X., Li, Z., Lu, C., Zhang, Q., Chang, L., et al. (2019). Transcriptome-wide analysis of N6-methyladenosine uncovers its regulatory role in gene expression in the lepidopteran Bombyx mori. Insect Mol. Biol. 28 (5), 703–715. doi:10.1111/imb.12584
Liu, J., Huang, T., Chen, W., Ding, C., Zhao, T., Zhao, X., et al. (2022). Developmental mRNA m5C landscape and regulatory innovations of massive m5C modification of maternal mRNAs in animals. Nat. Commun. 13 (1), 2484. doi:10.1038/s41467-022-30210-0
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262
Lockett, G. A., Kucharski, R., and Maleszka, R. (2012). DNA methylation changes elicited by social stimuli in the brains of worker honey bees. Genes Brain Behav. 11 (2), 235–242. doi:10.1111/j.1601-183X.2011.00751.x
Lourenço, A. P., Mackert, A., Cristino, Ad.S., and Simões, Z. L. P. (2008). Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR. Apidologie 39, 372–385. doi:10.1051/apido:2008015
Lourenço, A. P., Martins, J. R., Torres, F. A. S., Mackert, A., Aguiar, L. R., Hartfelder, K., et al. (2019). Immunosenescence in honey bees (Apis mellifera L.) is caused by intrinsic senescence and behavioral physiology. Exp. Gerontol. 119, 174–183. doi:10.1016/j.exger.2019.02.005
Martelli, F., Falcon, T., Pinheiro, D. G., Simões, Z. L. P., and Nunes, F. M. F. (2022). Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging. Insect Biochem. Mol. Biol. 146, 103774. doi:10.1016/j.ibmb.2022.103774
Meyer, K. D., Patil, D. P., Zhou, J., Zinoviev, A., Skabkin, M. A., Elemento, O., et al. (2015). 5' UTR m(6)A promotes cap-independent translation. Cell 163 (4), 999–1010. doi:10.1016/j.cell.2015.10.012
Meyer, K. D., Saletore, Y., Zumbo, P., Elemento, O., Mason, C. E., and Jaffrey, S. R. (2012). Comprehensive analysis of mRNA methylation reveals enrichment in 3' UTRs and near stop codons. Cell 149 (7), 1635–1646. doi:10.1016/j.cell.2012.05.003
Micura, R., Pils, W., Höbartner, C., Grubmayr, K., Ebert, M. O., and Jaun, B. (2001). Methylation of the nucleobases in RNA oligonucleotides mediates duplex-hairpin conversion. Nucleic Acids Res. 29 (19), 3997–4005. doi:10.1093/nar/29.19.3997
Page, R. E., and Peng, C. Y. (2001). Aging and development in social insects with emphasis on the honey bee, Apis mellifera L. Exp. Gerontol. 36 (4-6), 695–711. doi:10.1016/s0531-5565(00)00236-9
Page, R. E., Scheiner, R., Erber, J., and Amdam, G. V. (2006). The development and evolution of division of labor and foraging specialization in a social insect (Apis mellifera L.). Curr. Top. Dev. Biol. 74, 253–286. doi:10.1016/S0070-2153(06)74008-X
Rai, K., Chidester, S., Zavala, C. V., Manos, E. J., James, S. R., Karpf, A. R., et al. (2007). Dnmt2 functions in the cytoplasm to promote liver, brain, and retina development in zebrafish. Genes Dev. 21 (3), 261–266. doi:10.1101/gad.1472907
Schaefer, M., Kapoor, U., and Jantsch, M. F. (2017). Understanding RNA modifications: The promises and technological bottlenecks of the 'epitranscriptome. Open Biol. 7 (5), 170077. doi:10.1098/rsob.170077
Schaefer, M., Pollex, T., Hanna, K., Tuorto, F., Meusburger, M., Helm, M., et al. (2010). RNA methylation by Dnmt2 protects transfer RNAs against stress-induced cleavage. Genes Dev. 24 (15), 1590–1595. doi:10.1101/gad.586710
Schosserer, M., Minois, N., Angerer, T. B., Amring, M., Dellago, H., Harreither, E., et al. (2015). Methylation of ribosomal RNA by NSUN5 is a conserved mechanism modulating organismal lifespan. Nat. Commun. 6, 6158. doi:10.1038/ncomms7158
Shinoda, S., Kitagawa, S., Nakagawa, S., Wei, F. Y., Tomizawa, K., Araki, K., et al. (2019). Mammalian NSUN2 introduces 5-methylcytidines into mitochondrial tRNAs. Nucleic Acids Res. 47 (16), 8734–8745. doi:10.1093/nar/gkz575
Skowronek, P., Wójcik, Ł., and Strachecka, A. (2021). Fat body-multifunctional insect tissue. Insects 12 (6), 547. doi:10.3390/insects12060547
Tang, J., Chen, X., Cai, B., and Chen, G. (2019). A logical relationship for schizophrenia, bipolar, and major depressive disorder. Part 4: Evidence from chromosome 4 high-density association screen. J. Comp. Neurol. 527 (2), 392–405. doi:10.1002/cne.24543
Toth, A. L., and Robinson, G. E. (2005). Worker nutrition and division of labour in honeybees. Anim. Behav. 69 (2), 427–435. doi:10.1016/j.anbehav.2004.03.017
Tuorto, F., Liebers, R., Musch, T., Schaefer, M., Hofmann, S., Kellner, S., et al. (2012). RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat. Struct. Mol. Biol. 19 (9), 900–905. doi:10.1038/nsmb.2357
Vaiserman, A. (2014). Developmental epigenetic programming of caste-specific differences in social insects: An impact on longevity. Curr. Aging Sci. 7 (3), 176–186. doi:10.2174/1874609807666141129173749
Wang, M., Xiao, Y., Li, Y., Wang, X., Qi, S., Wang, Y., et al. (2021). RNA m6A modification functions in larval development and caste differentiation in honeybee (Apis mellifera). Cell Rep. 34 (1), 108580. doi:10.1016/j.celrep.2020.108580
Whitfield, C. W., Cziko, A. M., and Robinson, G. E. (2003). Gene expression profiles in the brain predict behavior in individual honey bees. Science 302 (5643), 296–299. doi:10.1126/science.1086807
Wilkinson, E., Cui, Y. H., and He, Y. Y. (2021). Context-dependent roles of RNA modifications in stress responses and diseases. Int. J. Mol. Sci. 22 (4), 1949. doi:10.3390/ijms22041949
Winston, M. L. (1987). “The honey bee colony: Life history,” in The hive and the honey bee. Editor J. M. Graham (Illinois: Dadant & Sons-Hamilton).
Withers, G. S., Fahrbach, S. E., and Robinson, G. E. (1993). Selective neuroanatomical plasticity and division of labour in the honeybee. Nature 364 (6434), 238–240. doi:10.1038/364238a0
Wu, J., Frazier, K., Zhang, J., Gan, Z., Wang, T., and Zhong, X. (2020). Emerging role of m6 A RNA methylation in nutritional physiology and metabolism. Obes. Rev. 21 (1), e12942. doi:10.1111/obr.12942
Yang, X., Wei, X., Yang, J., Du, T., Yin, C., Fu, B., et al. (2021). Epitranscriptomic regulation of insecticide resistance. Sci. Adv. 7 (19), eabe5903. doi:10.1126/sciadv.abe5903
Zayed, A., and Robinson, G. E. (2012). Understanding the relationship between brain gene expression and social behavior: Lessons from the honey bee. Annu. Rev. Genet. 46, 591–615. doi:10.1146/annurev-genet-110711-155517
Zhang, T., Chen, P., Li, W., Sha, S., Wang, Y., Yuan, Z., et al. (2019). Cognitive deficits in mice lacking Nsun5, a cytosine-5 RNA methyltransferase, with impairment of oligodendrocyte precursor cells. Glia 67 (4), 688–702. doi:10.1002/glia.23565
Zhang, Z., Wang, M., Xie, D., Huang, Z., Zhang, L., Yang, Y., et al. (2018). METTL3-mediated N6-methyladenosine mRNA modification enhances long-term memory consolidation. Cell Res. 28 (11), 1050–1061. doi:10.1038/s41422-018-0092-9
Zhao, X., Yang, Y., Sun, B. F., Shi, Y., Yang, X., Xiao, W., et al. (2014). FTO-dependent demethylation of N6-methyladenosine regulates mRNA splicing and is required for adipogenesis. Cell Res. 24 (12), 1403–1419. doi:10.1038/cr.2014.151
Keywords: RNA methylation, m6A, m5C, epitranscriptomics, bee, nutrition, aging, behavior
Citation: Bataglia L, Simões ZLP and Nunes FMF (2022) Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers. Front. Cell Dev. Biol. 10:921503. doi: 10.3389/fcell.2022.921503
Received: 15 April 2022; Accepted: 01 August 2022;
Published: 29 August 2022.
Edited by:
Juan Rafael Riesgo-Escovar, Universidad Nacional Autónoma de México, MexicoReviewed by:
Elizabeth Jenness Duncan, University of Leeds, United KingdomAlexandra Chittka, University College London, United Kingdom
Tomas Erban, Crop Research Institute (CRI), Czechia
Copyright © 2022 Bataglia, Simões and Nunes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Francis Morais Franco Nunes, francis.nunes@ufscar.br