Active genic machinery for epigenetic RNA modifications in bees

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

N6-methyladenosine modification of RNA controls dopamine synthesis to influence labour division in ants

The N6-methyladenosine (m6A) modification of RNA has been reported to remodel gene expression in response to environmental conditions; however, the biological role of m6A in social insects remains largely unknown. In this study, we explored the role of m6A in the division of labour by worker ants (Solenopsis invicta). We first determined the presence of m6A in RNAs from the brains of worker ants and found that m6A methylation dynamics differed between foragers and nurses. Depletion of m6A methyltransferase or chemical suppression of m6A methylation in foragers resulted in a shift to 'nurse-like' behaviours. Specifically, mRNAs of dopamine receptor 1 (Dop1) and dopamine transporter (DAT) were modified by m6A, and their expression increased dopamine levels to promote the behavioural transition from foragers to nurses. The abundance of Dop1 and DAT mRNAs and their stability were reduced by the inhibition of m6A modification caused by the silencing of Mettl3, suggesting that m6A modification in worker ants modulates dopamine synthesis, which regulates labour division. Collectively, our results provide the first example of the epitranscriptomic regulation of labour division in social insects and implicate m6A regulatory mechanism as a potential novel target for controlling red imported fire ants.

The current and emerging Klotho-enhancement strategies

Klotho is well known as a gene with antiaging properties. It has membrane and soluble forms, providing a unique system that controls various metabolic processes essential to health and disease. Klotho deficiency has been revealed to be associated with various aging-related disorders. Based on its various known and unknown protective properties, upregulating the Klotho gene may be a possible therapeutic and/or preventive approach in aging-related complications. Some agents, such as hormonal compounds, renin-angiotensin system inhibitors, antioxidants, peroxisome proliferator-activated receptor gamma (PPAR-γ) agonists, statins, vitamin D receptor agonists, antioxidants, anti-inflammatory agents, mammalian target of rapamycin (mTOR) signaling inhibitors, and receptor-interacting serine/threonine-protein kinase 1 (RIPK1) inhibitors, can possibly lead to the upregulation and elevation of Klotho levels. Demethylation and deacetylation of the Klotho gene can also be considered other possible Klotho-enhancement methods. Some emerging techniques, such as RNA modifications, gene therapy, gene editing, and exosome therapy, probably have the potential to be applied for increasing Klotho. In the present study, these current and emerging Klotho-enhancement strategies and their underlying mechanisms were comprehensively reviewed, which could highlight some potential avenues for future research.

ALKBH8 as a potential N6 -methyladenosine (m6 A) eraser in insects

The N6 -methyladenosine (m6 A) machinery functions through three groups of proteins in eukaryotic cells, including m6 A writers, erasers and readers. The m6 A cellular machinery has mostly been characterised in mammalian species, and the relevant literature on insects is currently scant. While homologues of m6 A writers and readers have been reported from insects, no erasers have been described so far. Here, using BLAST search, we searched for potential erasers in insects. While we found homologues of human m6 A eraser ALKBH5 in termites, beetles and true bugs, they could not be found in representative dipteran and lepidopteran species. However, a potential m6 A eraser, ALKBH8, was identified and experimentally investigated. Our results showed that ALKBH8 can reduce the m6 A levels of Aedes aegypti and Drosophila melanogaster RNAs, suggesting that AeALKBH8 could be a candidate m6 A eraser in insects.

N6-adenosine (m6A) mRNA methylation is required for Tribolium castaneum development and reproduction

Gene expression is regulated at various levels, including post-transcriptional mRNA modifications, where m6A methylation is the most common modification of mRNA. The m6A methylation regulates multiple stages of mRNA processing, including splicing, export, decay, and translation. How m6A modification is involved in insect development is not well known. We used the red flour beetle, Tribolium castaneum, as a model insect to identify the role of m6A modification in insect development. RNA interference (RNAi)-mediated knockdown of genes coding for m6A writers (m6A methyltransferase complex, depositing m6A to mRNA) and readers (YTH-domain proteins, recognizing and executing the function of m6A) was conducted. Knockdown of most writers during the larval stage caused a failure of ecdysis during eclosion. The loss of m6A machinery sterilized both females and males by interfering with the functioning of reproductive systems. Females treated with dsMettl3, the main m6A methyltransferase, laid significantly fewer and reduced-size eggs than the control insects. In addition, the embryonic development in eggs laid by dsMettl3 injected females was terminated in the early stages. Knockdown studies also showed that the cytosol m6A reader, YTHDF, is likely responsible for executing the function of m6A modifications during insect development. These data suggest that m6A modifications are critical for T. castaneum development and reproduction.

Exploring RNA methylation as a promising biomarker for assessing sublethal effects of fipronil on honeybees (Apis mellifera L.)

Honeybees play a crucial role as pollinators for crops and are regarded as sensitive bioindicators of environmental health. The widespread use of pesticides poses a severe threat to honeybee survival. However, there is limited information available on the specific risks associated with fipronil exposure in honeybees, particularly concerning the impact on RNA methylation throughout their lifespan. This study aimed to evaluate the effects of sublethal concentrations of fipronil on RNA m⁶A and m⁵C methylations, along with the associated genes in honeybee larvae and newly emerged adults. LC-MS/MS analysis revealed a notable hypomethylation of m⁵C in larvae, while hypermethylation of m⁶A was observed in the adult brain. Significant changes in the expression of genes such as AmWTAP, AmYTHDF, AmALKBH4, AmALKBH6, AmALKBH8, AmNSUN5, AmNOP2, AmTET1, and AmYBX1 were observed in the adult brain, whereas alterations in the expression of AmNSUN2, AmMETTL14, AmALKBH1, AmALKBH4, AmALKBH6 AmALYREF, AmTET1, and AmYBX1 were observed in the larvae. Notably, the expression of AmALKBH1 was not detected in any fipronil-treated larvae, suggesting its potential as an early risk indicator for honeybee larvae in future assessments. This pioneering study provides insights into the effects of fipronil on RNA methylations in honeybees and explores the possibility of employing RNA methylation as a tool for assessing pesticide risks in this important pollinator species. These findings offer new perspectives on honeybee protection and the development of toxicity evaluation systems for pesticides.

Examining parent-of-origin effects on transcription and RNA methylation in mediating aggressive behavior in honey bees (Apis mellifera)

Conflict between genes inherited from the mother (matrigenes) and the father (patrigenes) is predicted to arise during social interactions among offspring if these genes are not evenly distributed among offspring genotypes. This intragenomic conflict drives parent-specific transcription patterns in offspring resulting from parent-specific epigenetic modifications. Previous tests of the kinship theory of intragenomic conflict in honey bees (Apis mellifera) provided evidence in support of theoretical predictions for variation in worker reproduction, which is associated with extreme variation in morphology and behavior. However, more subtle behaviors - such as aggression - have not been extensively studied. Additionally, the canonical epigenetic mark (DNA methylation) associated with parent-specific transcription in plant and mammalian model species does not appear to play the same role as in honey bees, and thus the molecular mechanisms underlying intragenomic conflict in this species is an open area of investigation. Here, we examined the role of intragenomic conflict in shaping aggression in honey bee workers through a reciprocal cross design and Oxford Nanopore direct RNA sequencing. We attempted to probe the underlying regulatory basis of this conflict through analyses of parent-specific RNA m6A and alternative splicing patterns. We report evidence that intragenomic conflict occurs in the context of honey bee aggression, with increased paternal and maternal allele-biased transcription in aggressive compared to non-aggressive bees, and higher paternal allele-biased transcription overall. However, we found no evidence to suggest that RNA m6A or alternative splicing mediate intragenomic conflict in this species.

Putting hornets on the genomic map

Hornets are the largest of the social wasps, and are important regulators of insect populations in their native ranges. Hornets are also very successful as invasive species, with often devastating economic, ecological and societal effects. Understanding why these wasps are such successful invaders is critical to managing future introductions and minimising impact on native biodiversity. Critical to the management toolkit is a comprehensive genomic resource for these insects. Here we provide the annotated genomes for two hornets, Vespa crabro and Vespa velutina. We compare their genomes with those of other social Hymenoptera, including the northern giant hornet Vespa mandarinia. The three hornet genomes show evidence of selection pressure on genes associated with reproduction, which might facilitate the transition into invasive ranges. Vespa crabro has experienced positive selection on the highest number of genes, including those putatively associated with molecular binding and olfactory systems. Caste-specific brain transcriptomic analysis also revealed 133 differentially expressed genes, some of which are associated with olfactory functions. This report provides a spring-board for advancing our understanding of the evolution and ecology of hornets, and opens up opportunities for using molecular methods in the future management of both native and invasive populations of these over-looked insects.

Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers

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.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

Characterisation of the epigenetic architecture of the major malaria vector Anopheles arabiensis (Diptera: Culicidae) after treatment with epigenetic modulators and heavy metals

Anopheles arabiensis (a member of the An. gambiae species complex) is a major vector of malaria in sub-Saharan Africa. Despite its disease vector status, there is currently a paucity of epigenetic information for this species. The aim this study was therefore to analyse global epigenetic markers and their response to metal exposure in insecticide susceptible and resistant laboratory strains of An. arabiensis. This was done using commercially available epigenetic marker quantification kits. In order to validate the efficacy of the kits, several kits were assessed to determine whether changes induced by known epigenetic modulators were detectable using these platforms. The efficacy of the dosages used were determined by examining the effect of the dosages used on insecticide resistant phenotypes. Upon confirmation that the dosages used were sufficient to induce a phenotypic change, the effect on epigenetic markers was assessed. Commercial kits were used to quantify 5-methylcysteine (5-mC) and 5-hydroxymethylcysteine (5-hmC) methylation in DNA, m6A methylation in mRNA as well as Histone Acetyl Transferase (HAT) activity. There was a marked difference in the phenotypic response in adult mosquitoes of the insecticide susceptible strain compared to that of its' resistant counterpart. For males and females of the resistant strain, exposure to nucleic acid modifying drugs typically increased their tolerance to insecticides. The patterns of changes in 5-mC methylation by epigenetic modulators was congruent with previous studies which quantified by mass spectrometry. The two strains differed in methylation patterns under control conditions and responded differentially to larval metal exposure. In the resistant strain, which previously was demonstrated to show increased detoxification enzyme activity and insecticide tolerance after the same treatment, the potential increase in transcriptional activity appeared to be modulated by reduced methylation and increased HAT activity. This study suggests that the commercial epigenetic quantification kits can be used to characterise phenotypic changes in An. arabiensis, and also shows that epigenetic regulation of the response to metal exposure is regulated at the DNA as opposed to the RNA level.

Dual isoform sequencing reveals complex transcriptomic and epitranscriptomic landscapes of a prototype baculovirus

In this study, two long-read sequencing (LRS) techniques, MinION from Oxford Nanopore Technologies and Sequel from the Pacific Biosciences, were used for the transcriptional characterization of a prototype baculovirus, Autographa californica multiple nucleopolyhedrovirus. LRS is able to read full-length RNA molecules, and thereby distinguish between transcript isoforms, mono- and polycistronic RNAs, and overlapping transcripts. Altogether, we detected 875 transcript species, of which 759 were novel and 116 were annotated previously. These RNA molecules include 41 novel putative protein coding transcripts [each containing 5'-truncated in-frame open reading frames (ORFs), 14 monocistronic transcripts, 99 polygenic RNAs, 101 non-coding RNAs, and 504 untranslated region isoforms. This work also identified novel replication origin-associated transcripts, upstream ORFs, cis-regulatory sequences and poly(A) sites. We also detected RNA methylation in 99 viral genes and RNA hyper-editing in the longer 5'-UTR transcript isoform of the canonical ORF 19 transcript.