Comprehensive DNA methylation analysis of the Aedes aegypti genome

Multigenerational expression of antimicrobial peptides in Aedes aegypti exposed to Metarhizium anisopliae: Is trans-generational immune priming involved?

We assessed, for the first time, a multigenerational expression of antimicrobial peptides (AMPs) in Aedes aegypti larvae exposed to the entomopathogenic fungus, Metarhizium anisopliae, and correlated it with a possible involvement in trans-generational immune priming (TGIP). Aedes aegypti larvae were first exposed to blastospores or conidia of M. anisopliae CG 489 for 24 and 48 h, and the relative expression of AMPs were measured using quantitative Real-Time PCR. A suspension of conidia was prepared, and two different survival tests were conducted with different larval generations (F0, F1, and F2). In the first bioassay, the survival curves of the three generations were conducted separately and compared with their respective control groups. In the other bioassay, the survival curves of the F0, F1, and F2 generations were compared simultaneously against a naïve group exposed to Tween 80. In both survival tests, the F0 generation was more susceptible to M. anisopliae than subsequent generations. For molecular analyses related to TGIP, F0, F1, and F2 larvae were exposed to conidia, and their expression of AMPs was compared with their control groups and a naïve group. There was no differential expression of cecropin, defensin A or cathepsin B between generations. Lysozyme C, however, showed an increase in expression across generations, suggesting a role in TGIP. These discoveries may help us develop biological insecticides against mosquito larvae based on entomopathogenic fungi.

The mechanisms and factors that induce trained immunity in arthropods and mollusks

Besides dividing the organism's immune system into adaptive and innate immunity, it has long been thought that only adaptive immunity can establish immune memory. However, many studies have shown that innate immunity can also build immunological memory through epigenetic reprogramming and modifications to resist pathogens' reinfection, known as trained immunity. This paper reviews the role of mitochondrial metabolism and epigenetic modifications and describes the molecular foundation in the trained immunity of arthropods and mollusks. Mitochondrial metabolism and epigenetic modifications complement each other and play a key role in trained immunity.

Single nucleotide polymorphism (SNP) in the doublesex ( dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex ( dsx) gene, involved in mosquitos' somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene ( Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, involved in mosquitos’ somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

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.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: The malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, encoding somatic sexually dimorphic traits in mosquitoes. However, the mechanism that regulates the expression of this gene in anopheline mosquitoes is poorly understood. This study aimed to screen the An. gambiae dsx gene splice site sequences for single nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Evidence of Adaptive Evolution in Wolbachia-Regulated Gene DNMT2 and Its Role in the Dipteran Immune Response and Pathogen Blocking

Eukaryotic nucleic acid methyltransferase (MTase) proteins are essential mediators of epigenetic and epitranscriptomic regulation. DNMT2 belongs to a large, conserved family of DNA MTases found in many organisms, including holometabolous insects such as fruit flies and mosquitoes, where it is the lone MTase. Interestingly, despite its nomenclature, DNMT2 is not a DNA MTase, but instead targets and methylates RNA species. A growing body of literature suggests that DNMT2 mediates the host immune response against a wide range of pathogens, including RNA viruses. Curiously, although DNMT2 is antiviral in Drosophila, its expression promotes virus replication in mosquito species. We, therefore, sought to understand the divergent regulation, function, and evolution of these orthologs. We describe the role of the Drosophila-specific host protein IPOD in regulating the expression and function of fruit fly DNMT2. Heterologous expression of these orthologs suggests that DNMT2's role as an antiviral is host-dependent, indicating a requirement for additional host-specific factors. Finally, we identify and describe potential evidence of positive selection at different times throughout DNMT2 evolution within dipteran insects. We identify specific codons within each ortholog that are under positive selection and find that they are restricted to four distinct protein domains, which likely influence substrate binding, target recognition, and adaptation of unique intermolecular interactions. Collectively, our findings highlight the evolution of DNMT2 in Dipteran insects and point to structural, regulatory, and functional differences between mosquito and fruit fly homologs.

Mechanisms of Attenuation by Genetic Recoding of Viruses

The development of safe and effective vaccines against viruses is central to disease control. With advancements in DNA synthesis technology, the production of synthetic viral genomes has fueled many research efforts that aim to generate attenuated viruses by introducing synonymous mutations. Elucidation of the mechanisms underlying virus attenuation through synonymous mutagenesis is revealing interesting new biology that can be exploited for vaccine development. Here, we review recent advancements in this field of synthetic virology and focus on the molecular mechanisms of attenuation by genetic recoding of viruses. We highlight the action of the zinc finger antiviral protein (ZAP) and RNase L, two proteins involved in the inhibition of viruses enriched for CpG and UpA dinucleotides, that are often the products of virus recoding algorithms. Additionally, we discuss current challenges in the field as well as studies that may illuminate how other host functions, such as translation, are potentially involved in the attenuation of recoded viruses.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

DNA Methylation in Anopheles albimanus Modulates the Midgut Immune Response Against Plasmodium berghei

Epigenetic mechanisms such as DNA methylation and histone post-translational modifications are fundamental for the phenotypic plasticity of insects during their interaction with the environment. In response to environmental cues, the methylation pattern in DNA is dynamically remodeled to achieve an epigenetic control of gene expression. DNA methylation is the focus of study in insects for its evolutionarily conserved character; however, there is scant knowledge about the epigenetic regulation in vector mosquitoes, especially during their infection by parasites. The aim of the present study was to evaluate the participation of DNA methylation in the immune response of Anopheles albimanus to a Plasmodium infection. For this, we first investigated the presence of a fully functional DNA methylation system in A. albimanus by assessing its potential role in larval development. Subsequently, we evaluated the transcriptional response to Plasmodium berghei of two mosquito phenotypes with different degrees of susceptibility to the parasite, in a scenario where their global DNA methylation had been pharmacologically inhibited. Our study revealed that A. albimanus has a functional DNA methylation system that is essential to larval viability, and that is also responsive to feeding and parasite challenges. The pharmacological erasure of the methylome with azacytidine or decitabine abolished the divergent responses of both mosquito phenotypes, leading to a transcriptionally similar response upon parasite challenge. This response was more specific, and the infection load in both phenotypes was lowered. Our findings suggest that DNA methylation may constitute a key factor in vector competence, and a promising target for preventing malaria transmission.

Induction of DNA methyltransferase genes in Helicoverpa armigera following injection of pathogenic bacteria modulates expression of antimicrobial peptides and affects bacterial proliferation

Following pathogen attack in a host, widespread changes are induced in the host's gene expression, in particular those involved in the immune system, growth and survival. Epigenetic mechanisms have been suggested to be involved in the regulation of these changes through a number of mechanisms. DNA methylation is one of the important epigenetic processes that is carried out by DNA (cytosine-5) methyltransferase (DNMT) and alters expression of target genes. Here, we identified two putative sequences of DNMT (i.e. DNMT1 and DNMT2) from the transcriptome dataset of Helicoverpa armigera that showed high similarity to the homologous sequences in Bombyx mori. Domain architectures of DNMT1 and DNMT2 exhibit the unique pattern of DNMTs that highlights conserved function of these genes in different insects. To see if these genes play any role in bacterial infection, we challenged the fifth instar larvae of H. armigera by injecting Bacillus thuringiensis and Serratia marcescens cells into the hemolymph. Transcript levels of the DNMTs were analyzed by RT-qPCR. The results showed that the expression levels of DNMT1 and DNMT2 increased in the bacteria-injected larvae. Injection of the heat-killed bacteria also induced the expression of the DNMTs, but lower than that of the live bacteria. To determine whether these genes function during bacterial infection, we injected the inhibitor of DNMTs, 5-azacytidine (5-AZA), into the larvae and 24 h later, the bacterial cells were also injected into the larvae. Bacterial replication and larval mortality were analyzed in the treated and control insects. We found that 5-AZA reduced bacterial replication and also mortality of the bacterial-injected larvae regardless of the pathogenic bacterial species. Interestingly, the expression levels of antimicrobial peptides (AMPs) were also modulated following 5-AZA treatment. In conclusion, we showed that upregulation of the DNMTs in H. armigera following bacterial infections modulates AMPs and thereby affects the insect-bacteria interactions.

Quantitative methods for assessing local and bodywide contributions to Wolbachia titer in maternal germline cells of Drosophila

BACKGROUND

Little is known about how bacterial endosymbionts colonize host tissues. Because many insect endosymbionts are maternally transmitted, egg colonization is critical for endosymbiont success. Wolbachia bacteria, carried by approximately half of all insect species, provide an excellent model for characterizing endosymbiont infection dynamics. To date, technical limitations have precluded stepwise analysis of germline colonization by Wolbachia. It is not clear to what extent titer-altering effects are primarily mediated by growth rates of Wolbachia within cell lineages or migration of Wolbachia between cells.

RESULTS

The objective of this work is to inform mechanisms of germline colonization through use of optimized methodology. The approaches are framed in terms of nutritional impacts on Wolbachia. Yeast-rich diets in particular have been shown to suppress Wolbachia titer in the Drosophila melanogaster germline. To determine the extent of Wolbachia sensitivity to diet, we optimized 3-dimensional, multi-stage quantification of Wolbachia titer in maternal germline cells. Technical and statistical validation confirmed the identity of Wolbachia in vivo, the reproducibility of Wolbachia quantification and the statistical power to detect these effects. The data from adult feeding experiments demonstrated that germline Wolbachia titer is distinctly sensitive to yeast-rich host diets in late oogenesis. To investigate the physiological basis for these nutritional impacts, we optimized methodology for absolute Wolbachia quantification by real-time qPCR. We found that yeast-rich diets exerted no significant effect on bodywide Wolbachia titer, although ovarian titers were significantly reduced. This suggests that host diets affects Wolbachia distribution between the soma and late stage germline cells. Notably, relative qPCR methods distorted apparent wsp abundance, due to altered host DNA copy number in yeast-rich conditions. This highlights the importance of absolute quantification data for testing mechanistic hypotheses.

CONCLUSIONS

We demonstrate that absolute quantification of Wolbachia, using well-controlled cytological and qPCR-based methods, creates new opportunities to determine how bacterial abundance within the germline relates to bacterial distribution within the body. This methodology can be applied to further test germline infection dynamics in response to chemical treatments, genetic conditions, new host/endosymbiont combinations, or potentially adapted to analyze other cell and tissue types.

Changes in gene DNA methylation and expression networks accompany caste specialization and age-related physiological changes in a social insect

Social insects provide systems for studying epigenetic regulation of phenotypes, particularly with respect to differentiation of reproductive and worker castes, which typically arise from a common genetic background. The role of gene expression in caste specialization has been extensively studied, but the role of DNA methylation remains controversial. Here, we perform well replicated, integrated analyses of DNA methylation and gene expression in brains of an ant (Formica exsecta) with distinct female castes using traditional approaches (tests of differential methylation) combined with a novel approach (analysis of co-expression and co-methylation networks). We found differences in expression and methylation profiles between workers and queens at different life stages, as well as some overlap between DNA methylation and expression at the functional level. Large portions of the transcriptome and methylome are organized into "modules" of genes, some significantly associated with phenotypic traits of castes and developmental stages. Several gene co-expression modules are preserved in co-methylation networks, consistent with possible regulation of caste-specific gene expression by DNA methylation. Surprisingly, brain co-expression modules were highly preserved when compared with a previous study that examined whole-body co-expression patterns in 16 ant species, suggesting that these modules are evolutionarily conserved and for specific functions in various tissues. Altogether, these results suggest that DNA methylation participates in regulation of caste specialization and age-related physiological changes in social insects.

Whole-Genome Bisulfite Sequencing for the Methylation Analysis of Insect Genomes

DNA methylation is a conserved epigenetic modification of animal genomes, but genome methylation patterns appear surprisingly diverse in insects. Whole-genome bisulfite sequencing (WGBS) represents a sensitive and robust method for the characterization of genome-wide methylation patterns at single-base resolution. Here, we describe a step-by-step protocol for the generation and analysis of WGBS datasets using standard Illumina sequencing platforms. In comparison to whole-genome sequencing, WGBS has additional caveats that require particular attention and are highlighted in this chapter.

Signatures of DNA Methylation across Insects Suggest Reduced DNA Methylation Levels in Holometabola

It has been experimentally shown that DNA methylation is involved in the regulation of gene expression and the silencing of transposable element activity in eukaryotes. The variable levels of DNA methylation among different insect species indicate an evolutionarily flexible role of DNA methylation in insects, which due to a lack of comparative data is not yet well-substantiated. Here, we use computational methods to trace signatures of DNA methylation across insects by analyzing transcriptomic and genomic sequence data from all currently recognized insect orders. We conclude that: 1) a functional methylation system relying exclusively on DNA methyltransferase 1 is widespread across insects. 2) DNA methylation has potentially been lost or extremely reduced in species belonging to springtails (Collembola), flies and relatives (Diptera), and twisted-winged parasites (Strepsiptera). 3) Holometabolous insects display signs of reduced DNA methylation levels in protein-coding sequences compared with hemimetabolous insects. 4) Evolutionarily conserved insect genes associated with housekeeping functions tend to display signs of heavier DNA methylation in comparison to the genomic/transcriptomic background. With this comparative study, we provide the much needed basis for experimental and detailed comparative analyses required to gain a deeper understanding on the evolution and function of DNA methylation in insects.

Autosomal genetic variation is associated with DNA methylation in regions variably escaping X-chromosome inactivation

X-chromosome inactivation (XCI), i.e., the inactivation of one of the female X chromosomes, restores equal expression of X-chromosomal genes between females and males. However, ~10% of genes show variable degrees of escape from XCI between females, although little is known about the causes of variable XCI. Using a discovery data-set of 1867 females and 1398 males and a replication sample of 3351 females, we show that genetic variation at three autosomal loci is associated with female-specific changes in X-chromosome methylation. Through cis-eQTL expression analysis, we map these loci to the genes SMCHD1/METTL4, TRIM6/HBG2, and ZSCAN9. Low-expression alleles of the loci are predominantly associated with mild hypomethylation of CpG islands near genes known to variably escape XCI, implicating the autosomal genes in variable XCI. Together, these results suggest a genetic basis for variable escape from XCI and highlight the potential of a population genomics approach to identify genes involved in XCI.

Wolbachia wStri Blocks Zika Virus Growth at Two Independent Stages of Viral Replication

Mosquito-transmitted viruses are spread globally and present a great risk to human health. Among the many approaches investigated to limit the diseases caused by these viruses are attempts to make mosquitos resistant to virus infection. Coinfection of mosquitos with the bacterium Wolbachia pipientis from supergroup A is a recent strategy employed to reduce the capacity for major vectors in the Aedes mosquito genus to transmit viruses, including dengue virus (DENV), Chikungunya virus (CHIKV), and Zika virus (ZIKV). Recently, a supergroup B Wolbachia wStri, isolated from Laodelphax striatellus, was shown to inhibit multiple lineages of ZIKV in Aedes albopictus cells. Here, we show that wStri blocks the growth of positive-sense RNA viruses DENV, CHIKV, ZIKV, and yellow fever virus by greater than 99.9%. wStri presence did not affect the growth of the negative-sense RNA viruses LaCrosse virus or vesicular stomatitis virus. Investigation of the stages of the ZIKV life cycle inhibited by wStri identified two distinct blocks in viral replication. We found a reduction of ZIKV entry into wStri-infected cells. This was partially rescued by the addition of a cholesterol-lipid supplement. Independent of entry, transfected viral genome was unable to replicate in Wolbachia-infected cells. RNA transfection and metabolic labeling studies suggested that this replication defect is at the level of RNA translation, where we saw a 66% reduction in mosquito protein synthesis in wStri-infected cells. This study’s findings increase the potential for application of wStri to block additional arboviruses and also identify specific blocks in viral infection caused by Wolbachia coinfection.

Dengue, Zika, and yellow fever viruses are mosquito-transmitted diseases that have spread throughout the world, causing millions of infections and thousands of deaths each year. Existing programs that seek to contain these diseases through elimination of the mosquito population have so far failed, making it crucial to explore new ways of limiting the spread of these viruses. Here, we show that introduction of an insect symbiont Wolbachia wStri, into mosquito cells is highly effective at reducing yellow fever virus, dengue virus, Zika virus, and Chikungunya virus production. Reduction of virus replication was attributable to decreases in entry and a strong block of virus gene expression at the translational level. These findings expand the potential use of Wolbachia wStri to block viruses and identify two separate steps for limiting virus replication in mosquitos that could be targeted via microbes or other means as an antiviral strategy.

Epigenetic modifications and their relation to caste and sex determination and adult division of labor in the stingless bee Melipona scutellaris

Stingless bees of the genus Melipona, have long been considered an enigmatic case among social insects for their mode of caste determination, where in addition to larval food type and quantity, the genotype also has a saying, as proposed over 50 years ago by Warwick E. Kerr. Several attempts have since tried to test his Mendelian two-loci/two-alleles segregation hypothesis, but only recently a single gene crucial for sex determination in bees was evidenced to be sex-specifically spliced and also caste-specifically expressed in a Melipona species. Since alternative splicing is frequently associated with epigenetic marks, and the epigenetic status plays a major role in setting the caste phenotype in the honey bee, we investigated here epigenetic chromatin modification in the stingless bee Melipona scutellaris. We used an ELISA-based methodology to quantify global methylation status and western blot assays to reveal histone modifications. The results evidenced DNA methylation/demethylation events in larvae and pupae, and significant differences in histone methylation and phosphorylation between newly emerged adult queens and workers. The epigenetic dynamics seen in this stingless bee species represent a new facet in the caste determination process in Melipona bees and suggest a possible mechanism that is likely to link a genotype component to the larval diet and adult social behavior of these bees.