Function and evolution of DNA methylation in Nasonia vitripennis

DNA methylase 1 influences temperature responses and development in the invasive pest Tuta absoluta

DNA methylase 1 (Dnmt1) is an important regulatory factor associated with biochemical signals required for insect development. It responds to changes in the environment and triggers phenotypic plasticity. Meanwhile, Tuta absoluta Meyrick (Lepidoptera: Gelechiidae)-a destructive invasive pest-can rapidly invade and adapt to different habitats; however, the role of Dnmt1 in this organism has not been elucidated. Accordingly, this study investigates the mechanism(s) underlying the rapid adaptation of Tuta absoluta to temperature stress. Potential regulatory genes were screened via RNAi (RNA interference), and the DNA methylase in Tuta absoluta was cloned by RACE (Rapid amplification of cDNA ends). TaDnmt1 was identified as a potential regulatory gene via bioinformatics; its expression was evaluated in response to temperature stress and during different development stages using real-time polymerase chain reaction. Results revealed that TaDnmt1 participates in hot/cold tolerance, temperature preference and larval development. The full-length cDNA sequence of TaDnmt1 is 3765 bp and encodes a 1254 kDa protein with typical Dnmt1 node-conserved structural features and six conserved DNA-binding active motifs. Moreover, TaDnmt1 expression is significantly altered by temperature stress treatments and within different development stages. Hence, TaDnmt1 likely contributes to temperature responses and organismal development. Furthermore, after treating with double-stranded RNA and exposing Tuta absoluta to 35°C heat shock or -12°C cold shock for 1 h, the survival rate significantly decreases; the preferred temperature is 2°C lower than that of the control group. In addition, the epidermal segments become enlarged and irregularly folded while the surface dries up. This results in a significant increase in larval mortality (57%) and a decrease in pupation (49.3%) and eclosion (50.9%) rates. Hence, TaDnmt1 contributes to temperature stress responses and temperature perception, as well as organismal growth and development, via DNA methylation regulation. These findings suggest that the rapid geographic expansion of T absoluta has been closely associated with TaDnmt1-mediated temperature tolerance. This study advances the research on 'thermos Dnmt' and provides a potential target for RNAi-driven regulation of Tuta absoluta.

Exploring the ageing methylome in the model insect, Nasonia vitripennis

BACKGROUND

The ageing process is a multifaceted phenomenon marked by the gradual deterioration of cellular and organismal functions, accompanied by an elevated susceptibility to diseases. The intricate interplay between genetic and environmental factors complicates research, particularly in complex mammalian models. In this context, simple invertebrate organisms have been pivotal, but the current models lack detectable DNA methylation limiting the exploration of this critical epigenetic ageing mechanism. This study introduces Nasonia vitripennis, the jewel wasp, as an innovative invertebrate model for investigating the epigenetics of ageing. Leveraging its advantages as a model organism and possessing a functional DNA methylation system, Nasonia emerges as a valuable addition to ageing research.

RESULTS

Whole-genome bisulfite sequencing unveiled dynamic alterations in DNA methylation, with differentially methylated CpGs between distinct time points in both male and female wasps. These changes were associated with numerous genes, enriching for functions related to telomere maintenance, histone methylation, and mRNA catabolic processes. Additionally, other CpGs were found to be variably methylated at each timepoint. Sex-specific effects on epigenetic entropy were observed, indicating differential patterns in the loss of epigenetic stability over time. Constructing an epigenetic clock containing 19 CpGs revealed a robust correlation between epigenetic age and chronological age.

CONCLUSIONS

Nasonia vitripennis emerges as a promising model for investigating the epigenetics of ageing, shedding light on the intricate dynamics of DNA methylation and their implications for age-related processes. This research not only expands the repertoire of ageing models but also opens avenues for deeper exploration of epigenetic mechanisms in the context of ageing.

Modifications to gene body methylation do not alter gene expression plasticity in a reef-building coral

As coral reefs continue to decline due to climate change, the role of coral epigenetics (specifically, gene body methylation, GBM) in coral acclimatization warrants investigation. The evidence is currently conflicting. In diverse animal phyla, the baseline GBM level is associated with gene function: continuously expressed "housekeeping" genes are typically highly methylated, while inducible context-dependent genes have low or no methylation at all. Some authors report an association between GBM and the environment and interpret this observation as evidence of the GBM's role in acclimatization. Yet, others argue that the correlation between GBM change and gene expression change is typically absent or negligible. Here, we used the reef-building coral, Acropora millepora, to test whether environmentally driven changes in GBM are associated with a gene's ability to respond to environmental changes (plasticity) rather than expression level. We analyzed two cases of modified gene expression plasticity observed in a 3-week-long heat acclimatization experiment. The first one was a group of heat-induced genes that failed to revert their expression after the coral was translocated back to the control tank. The second case involved genes that changed the magnitude of their response to the daily temperature fluctuations over the course of the experiment. In both cases, we found negligible or no association with GBM change. We conclude that although both gene expression plasticity and GBM can change during acclimatization, there is no direct association between the two. This adds to the increasing volume of evidence that the function of GBM in invertebrates is unrelated to acclimatization on physiological timescales.

Metabolic and immune costs balance during natural acclimation of corals in fluctuating environments

Epigenetic modifications based on DNA methylation can rapidly improve the potential of corals to adapt to environmental pressures by increasing their phenotypic plasticity, a factor important for scleractinian corals to adapt to future global warming. However, the extent to which corals develop similar adaptive mechanisms and their specific adaptation processes remain unclear. Here, to reveal the regulatory mechanism by which DNA methylation improves thermal tolerance in Pocillopora damicornis under fluctuating environments, we analyzed genome-wide DNA methylation signatures in P. damicornis and compared the differences in the methylation and transcriptional responses of P. damicornis from fluctuating and stable environments using whole-genome bisulfite sequencing and nanopore-based RNA sequencingtranscriptome sequencing. We discovered low methylation levels in P. damicornis (average methylation 4.14%), with CpG accounting for 74.88%, CHH for 13.27%, and CHG for 11.85% of this methylation. However, methylation levels did not change between coral samples from the fluctuating and stable environments. The varied methylation levels in different regions of the gene revealed that the overall methylation level of the gene body was relatively high and showed a bimodal methylation pattern. Methylation occurs primarily in exons rather than introns within the gene body In P. damicornis, there was only a weak correlation between methylation and transcriptional changes at the individual gene level, and the methylation and gene expression levels generally exhibited a bell-shaped relationship, which we speculate may be due to the specificity of cnidarian species. Correlation analysis between methylation levels and the transcriptome revealed that the highest proportion of the top 20 enriched KEGG pathways was related to immunity. Additionally, P. damicornis collected from a high-temperature pool had a lower metabolic rate than those collected from a low-temperature pool. We hypothesize that the dynamic balance of energy-expenditure costs between immunity and metabolism is an important strategy for increasing P. damicornis tolerance. The fluctuating environment of high-temperature pools may increase the heat tolerance in corals by increasing their immunity and thus lowering their metabolism.

Coevolution of the CDCA7-HELLS ICF-related nucleosome remodeling complex and DNA methyltransferases

5-Methylcytosine (5mC) and DNA methyltransferases (DNMTs) are broadly conserved in eukaryotes but are also frequently lost during evolution. The mammalian SNF2 family ATPase HELLS and its plant ortholog DDM1 are critical for maintaining 5mC. Mutations in HELLS, its activator CDCA7, and the de novo DNA methyltransferase DNMT3B, cause immunodeficiency-centromeric instability-facial anomalies (ICF) syndrome, a genetic disorder associated with the loss of DNA methylation. We here examine the coevolution of CDCA7, HELLS and DNMTs. While DNMT3, the maintenance DNA methyltransferase DNMT1, HELLS, and CDCA7 are all highly conserved in vertebrates and green plants, they are frequently co-lost in other evolutionary clades. The presence-absence patterns of these genes are not random; almost all CDCA7 harboring eukaryote species also have HELLS and DNMT1 (or another maintenance methyltransferase, DNMT5). Coevolution of presence-absence patterns (CoPAP) analysis in Ecdysozoa further indicates coevolutionary linkages among CDCA7, HELLS, DNMT1 and its activator UHRF1. We hypothesize that CDCA7 becomes dispensable in species that lost HELLS or DNA methylation, and/or the loss of CDCA7 triggers the replacement of DNA methylation by other chromatin regulation mechanisms. Our study suggests that a unique specialized role of CDCA7 in HELLS-dependent DNA methylation maintenance is broadly inherited from the last eukaryotic common ancestor.

Genome-wide DNA methylation patterns in bumble bee (Bombus vosnesenskii) populations from spatial-environmental range extremes

Unraveling molecular mechanisms of adaptation to complex environments is crucial to understanding tolerance of abiotic pressures and responses to climatic change. Epigenetic variation is increasingly recognized as a mechanism that can facilitate rapid responses to changing environmental cues. To investigate variation in genetic and epigenetic diversity at spatial and thermal extremes, we use whole genome and methylome sequencing to generate a high-resolution map of DNA methylation in the bumble bee Bombus vosnesenskii. We sample two populations representing spatial and environmental range extremes (a warm southern low-elevation site and a cold northern high-elevation site) previously shown to exhibit differences in thermal tolerance and determine positions in the genome that are consistently and variably methylated across samples. Bisulfite sequencing reveals methylation characteristics similar to other arthropods, with low global CpG methylation but high methylation concentrated in gene bodies and in genome regions with low nucleotide diversity. Differentially methylated sites (n = 2066) were largely hypomethylated in the northern high-elevation population but not related to local sequence differentiation. The concentration of methylated and differentially methylated sites in exons and putative promoter regions suggests a possible role in gene regulation, and this high-resolution analysis of intraspecific epigenetic variation in wild Bombus suggests that the function of methylation in niche adaptation would be worth further investigation.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Genomic signatures associated with maintenance of genome stability and venom turnover in two parasitoid wasps

Parasitoid wasps are rapidly developing as a model for evolutionary biology. Here we present chromosomal genomes of two Anastatus wasps, A. japonicus and A. fulloi, and leverage these genomes to study two fundamental questions-genome size evolution and venom evolution. Anastatus shows a much larger genome than is known among other wasps, with unexpectedly recent bursts of LTR retrotransposons. Importantly, several genomic innovations, including Piwi gene family expansion, ubiquitous Piwi expression profiles, as well as transposable element-piRNA coevolution, have likely emerged for transposable element silencing to maintain genomic stability. Additionally, we show that the co-option evolution arose by expression shifts in the venom gland plays a dominant role in venom turnover. We also highlight the potential importance of non-venom genes that are coexpressed with venom genes during venom evolution. Our findings greatly advance the current understanding of genome size evolution and venom evolution, and these genomic resources will facilitate comparative genomics studies of insects in the future.

Dnmt1a is essential for gene body methylation and the regulation of the zygotic genome in a wasp

Gene body methylation (GBM) is an ancestral mode of DNA methylation whose role in development has been obscured by the more prominent roles of promoter and CpG island methylation. The wasp Nasonia vitripennis has little promoter and CpG island methylation, yet retains strong GBM, making it an excellent model for elucidating the roles of GBM. Here we show that N. vitripennis DNA methyltransferase 1a (Nv-Dnmt1a) knockdown leads to failures in cellularization and gastrulation of the embryo. Both of these disrupted events are hallmarks of the maternal-zygotic transition (MZT) in insects. Analysis of the embryonic transcriptome and methylome revealed strong reduction of GBM and widespread disruption of gene expression during embryogenesis after Nv-Dnmt1a knockdown. Strikingly, there was a strong correlation between loss of GBM and reduced gene expression in thousands of methylated loci, consistent with the hypothesis that GBM directly facilitates high levels of transcription. We propose that lower expression levels of methylated genes due to reduced GBM is the crucial direct effect of Nv-Dnmt1 knockdown. Subsequently, the disruption of methylated genes leads to downstream dysregulation of the MZT, culminating in developmental failure at gastrulation.

DNA methylation across the tree of life, from micro to macro-organism

DNA methylation is a process in which methyl (CH3) groups are added to the DNA molecule. The DNA segment does not change in the sequence, but DNA methylation could alter the action of DNA. Different enzymes like DNA methyltransferases (DNMTs) take part in methylation of cytosine/adenine nucleosides in DNA. In prokaryotes, DNA methylation is performed to prevent the attack of phage and also plays a role in the chromosome replication and repair. In fungi, DNA methylation is studied to see the transcriptional changes, as in insects, the DNA methylation is not that well-known, it plays a different role like other organisms. In mammals, the DNA methylation is related to different types of cancers and plays the most important role in the placental development and abnormal DNA methylation connected with diseases like cancer, autoimmune diseases, and rheumatoid arthritis.

Ecotoxicological epigenetics in invertebrates: Emerging tool for the evaluation of present and past pollution burden

The effect of environmental pollution on epigenetic changes and their heredity in affected organisms is of major concern as such changes can play a significant role in adaptation to changing environmental conditions. Changes of epigenetic marks including DNA methylation, histone modifications, and non-coding RNA's can induce changes in gene transcription leading to physiological long-term changes or even transgenerational inheritance. Such mechanisms have until recently been scarcely studied in invertebrate organisms, mainly focusing on model species including Caenorhabditis elegans and Daphnia magna. However, more data are becoming available, particularly focused on DNA methylation changes caused by anthropogenic pollutants in a wide range of invertebrates. This review examines the literature from field and laboratory studies utilising invertebrate species exposed to environmental pollutants and their effect on DNA methylation. Possible mechanisms of epigenetic modifications and their role on physiology and adaptation as well as the incidence of intergenerational and transgenerational inheritance are discussed. Furthermore, critical research challenges are defined and the way forward is proposed. Future studies should focus on the use of next generation sequencing tools to define invertebrate methylomes under environmental stress in higher resolution, those data should further be linked to gene expression patterns and phenotypes and detailed studies focusing on transgenerational effects are encouraged. Moreover, studies of other epigenetic mechanisms in various invertebrate species, apart from DNA methylation would provide better understanding of interconnected cross-talk between epigenetic marks. Taken together incorporating epigenetic studies in ecotoxicology context presents a promising tool for development of sensitive biomarkers for environmental stress assessment.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

Lineage and Parent-of-Origin Effects in DNA Methylation of Honey Bees (Apis mellifera) Revealed by Reciprocal Crosses and Whole-Genome Bisulfite Sequencing

Parent-of-origin methylation arises when the methylation patterns of a particular allele are dependent on the parent it was inherited from. Previous work in honey bees has shown evidence of parent-of-origin-specific expression, yet the mechanisms regulating such pattern remain unknown in honey bees. In mammals and plants, DNA methylation is known to regulate parent-of-origin effects such as genomic imprinting. Here, we utilize genotyping of reciprocal European and Africanized honey bee crosses to study genome-wide allele-specific methylation patterns in sterile and reproductive individuals. Our data confirm the presence of allele-specific methylation in honey bees in lineage-specific contexts but also importantly, though to a lesser degree, parent-of-origin contexts. We show that the majority of allele-specific methylation occurs due to lineage rather than parent-of-origin factors, regardless of the reproductive state. Interestingly, genes affected by allele-specific DNA methylation often exhibit both lineage and parent-of-origin effects, indicating that they are particularly labile in terms of DNA methylation patterns. Additionally, we re-analyzed our previous study on parent-of-origin-specific expression in honey bees and found little association with parent-of-origin-specific methylation. These results indicate strong genetic background effects on allelic DNA methylation and suggest that although parent-of-origin effects are manifested in both DNA methylation and gene expression, they are not directly associated with each other.

A review on the DNA methyltransferase family of insects: Aspect and prospects

The DNA methyltransferase family contains a conserved set of DNA-modifying enzymatic proteins. They are responsible for epigenetic gene modulation, such as transcriptional silencing, transcription activation, and post-transcriptional modulation. Recent research has revealed that the canonical DNA methyltransferases (DNMTs) biological roles go beyond their traditional functions of establishing and maintaining DNA methylation patterns. Although a complete DNA methylation toolkit is absent in most insect orders, recent evidence indicates the de novo DNA methylation and maintenance function remain conserved. Studies using various molecular approaches provided evidence that DNMTs are multi-functional proteins. However, still in-depth studies on their biological role lack due to the least studied area in insects. Here, we review the DNA methylation toolkit of insects, focusing on recent research on various insect orders, which exhibit DNA methylation at different levels, and for which DNMTs functional studies have become available in recent years. We survey research on the potential roles of DNMTs in the regulation of gene transcription in insect species. DNMTs participate in different physiological processes by interacting with other epigenetic factors. Future studies on insect's DNMTs will benefit to understand developmental processes, responses to various stimuli, and adaptability of insects to different environmental conditions.

Male-Biased microRNA Discovery in the Pea Aphid

Epigenetic mechanisms modulate gene expression levels during development, shaping how a single genome produces a diversity of phenotypes. Here, we begin to explore the epigenetic regulation of sexual dimorphism in pea aphids (Acyrthosiphon pisum) by focusing on microRNAs. Previous analyses of microRNAs in aphids have focused solely on females, so we performed deep sequencing of a sample containing early-stage males. We used this sample, plus samples from Genbank, to find 207 novel pea aphid microRNA coding loci. We localized microRNA loci to a chromosome-level assembly of the pea aphid genome and found that those on the X chromosome have lower overall expression compared to those on autosomes. We then identified a set of 19 putative male-biased microRNAs and found them enriched on the X chromosome. Finally, we performed protein-coding RNA-Seq of first instar female and male pea aphids to identify genes with lower expression in males. 10 of these genes were predicted targets of the 19 male-biased microRNAs. Our study provides the most complete set of microRNAs in the pea aphid to date and serves as foundational work for future studies on the epigenetic control of sexual dimorphism.

Inference of DNA methylation patterns in molluscs

Mollusca are the second largest and arguably most diverse phylum of the animal kingdom. This is in sharp contrast to our very limited knowledge concerning epigenetic mechanisms including DNA methylation in this invertebrate group. Here, we inferred DNA methylation patterns by analysing the normalized dinucleotide CG content in protein-coding sequences and identified DNA methyltransferases (DNMT1 and 3) in published transcriptomes and genomes of 140 species across all eight classes of molluscs. Given the evolutionary age and morphological diversity of molluscs, we expected to find evidence for diverse methylation patterns. Our inferences suggest that molluscs possess substantial levels of DNA methylation in gene bodies as a rule. Yet, we found deviations from this general picture with regard to (i) the CpG observed/expected distributions indicating a reduction in DNA methylation in certain groups and (ii) the completeness of the DNMT toolkit. Reductions were evident in Caudofoveata, Solenogastres, Polyplacophora, Monoplacophora, as well as Scaphopoda. Heterobranchia and Oegopsida were remarkable as they lacked DNMT3, usually responsible for de novo methylation, yet showed signs of DNA methylation. Our survey may serve as guidance for direct empirical analyses of DNA methylation in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Complete Metamorphosis in Manduca sexta Involves Specific Changes in DNA Methylation Patterns

The transition between morphologically distinct phenotypes during complete metamorphosis in holometabolous insects is accompanied by fundamental transcriptional reprogramming. Using the tobacco hornworm (Manduca sexta), a powerful model for the analysis of insect evolution and development, we conducted a genome-wide comparative analysis of gene expression and DNA methylation in caterpillars and adults to determine whether complete metamorphosis has an epigenetic basis in this species. Bisulfite sequencing indicated a generally low level of DNA methylation with a unimodal CpGO/E distribution. Expression analysis revealed that 24 % of all known M. sexta genes (3.729) were upregulated in last-instar larvae relative to the adult moth, whereas 26 % (4.077) were downregulated. We also identified 4.946 loci and 4.960 regions showing stage-specific differential methylation. Interestingly, genes encoding histone acetyltransferases and histone deacetylases were differentially methylated in the larvae and adults, indicating there is crosstalk between different epigenetic mechanisms. The distinct sets of methylated genes in M. sexta larvae and adults suggest that complete metamorphosis involves epigenetic modifications associated with profound transcriptional reprogramming, involving approximately half of all the genes in this species.

Sex-specific expression and DNA methylation in a species with extreme sexual dimorphism and paternal genome elimination

Phenotypic differences between sexes are often mediated by differential expression and alternative splicing of genes. However, the mechanisms that regulate these expression and splicing patterns remain poorly understood. The mealybug, Planococcus citri, displays extreme sexual dimorphism and exhibits an unusual instance of sex-specific genomic imprinting, paternal genome elimination (PGE), in which the paternal chromosomes in males are highly condensed and eliminated from the sperm. Planococcus citri has no sex chromosomes and both sexual dimorphism and PGE are predicted to be under epigenetic control. We recently showed that P. citri females display a highly unusual DNA methylation profile for an insect species, with the presence of promoter methylation associated with lower levels of gene expression. Here, we therefore decided to explore genome-wide differences in DNA methylation between male and female P. citri using whole-genome bisulphite sequencing. We identified extreme differences in genome-wide levels and patterns between the sexes. Males display overall higher levels of DNA methylation which manifest as more uniform low levels across the genome. Whereas females display more targeted high levels of methylation. We suggest these unique sex-specific differences are due to chromosomal differences caused by PGE and may be linked to possible ploidy compensation. Using RNA-Seq, we identify extensive sex-specific gene expression and alternative splicing, but we find no correlation with cis-acting DNA methylation.

The effect of DNA methylation on bumblebee colony development

BACKGROUND

Although around 1% of cytosines in bees' genomes are known to be methylated, less is known about methylation's effect on bee behavior and fitness. Chemically altered DNA methylation levels have shown clear changes in the dominance and reproductive behavior of workers in queen-less colonies, but the global effect of DNA methylation on caste determination and colony development remains unclear, mainly because of difficulties in controlling for genetic differences among experimental subjects in the parental line. Here, we investigated the effect of the methylation altering agent decitabine on the developmental rate of full bumblebee colonies. Whole genome bisulfite sequencing was used to assess differences in methylation status.

RESULTS

Our results showed fewer methylated loci in the control group. A total of 22 CpG loci were identified as significantly differentially methylated between treated and control workers with a change in methylation levels of 10% or more. Loci that were methylated differentially between groups participated in pathways including neuron function, oocyte regulation and metabolic processes. Treated colonies tended to develop faster, and therefore more workers were found at a given developmental stage. However, male production followed the opposite trend and it tended to be higher in control colonies.

CONCLUSION

Overall, our results indicate that altered methylation patterns resulted in an improved cooperation between workers, while there were no signs of abnormal worker dominance or caste determination.

Whole-genome bisulfite sequencing of abdominal adipose reveals DNA methylation pattern variations in broiler lines divergently selected for fatness

The methylation status of pivotal genes involved in fat deposition in chickens has been extensively studied. However, the whole-genome DNA methylation profiles of broiler abdominal adipose tissue remain poorly understood. Using whole-genome bisulfite sequencing, we generated DNA methylation profiles of chicken abdominal adipose tissue from Northeast Agricultural University broiler lines divergently selected for abdominal fat content. We aimed to explore whether DNA methylation was associated with abdominal fat deposition in broilers. The whole-genome DNA methylation profiles of fat- and lean-line broilers abdominal adipose tissue were constructed. The DNA methylation levels of functional genomic regions in the fat broiler were higher than those in the lean broiler, especially in the 3' untranslated regions (UTRs) and exons in the non-CG contexts. Additionally, we identified 29,631 differentially methylated regions and, subsequently, annotated 6,484 and 2,016 differentially methylated genes (DMGs) in the gene body and promoter regions between the two lines, respectively. Functional annotation showed that the DMGs in promoter regions were significantly enriched mainly in the triglyceride catabolic process, lipid metabolism-related pathways, and extracellular matrix signal pathways. When the DMG in promoter regions and differentially expressed genes were integrated, we identified 30 genes with DNA methylation levels that negatively correlated with their messenger RNA (mRNA) expression, of which CMSS1 reached significant levels (false discovery rate < 0.05). These 30 genes were mainly involved in fatty acid metabolism, peroxisome-proliferator-activated receptor signaling, Wnt signaling pathways, transmembrane transport, RNA degradation, and glycosaminoglycan degradation. Comparing the DNA methylation profiles between fat- and lean-line broilers demonstrated that DNA methylation is involved in regulating broiler abdominal fat deposition. Our study offers a basis for further exploring the underlying mechanisms of abdominal adipose deposition in broilers.

DNMT1 Promotes Genome Methylation and Early Embryo Development in Cockroaches

The influence of DNA methylation on gene behavior and its consequent phenotypic effects appear to be very important, but the details are not well understood. Insects offer a diversity of DNA methylation modes, making them an excellent lineage for comparative analyses. However, functional studies have tended to focus on quite specialized holometabolan species, such as wasps, bees, beetles, and flies. Here, we have studied DNA methylation in the hemimetabolan insect Blattella germanica. In this cockroach, a gene involved in DNA methylation, DNA methyltransferase 1 (DNMT1), is expressed in early embryogenesis. In our experiments, RNAi of DNMT1 reduces DNA methylation and impairs blastoderm formation. Using reduced representation bisulfite sequencing and transcriptome analyses, we observed that methylated genes are associated with metabolism and are highly expressed, whereas unmethylated genes are related to signaling and show low expression. Moreover, methylated genes show greater expression change and less expression variability than unmethylated genes.

Comparative genome-wide DNA methylation analysis reveals epigenomic differences in response to heat-humidity stress in Bombyx mori

DNA methylation is an important epigenetic modification and has been shown to be involved in the response to abiotic stress. However, there are few studies on DNA methylation in insect response to environmental signals. In this study, we conducted a comprehensive comparative analysis of DNA methylation profiles between two silkworm strains with significantly different resistance to heat and humidity by whole-genome bisulfite sequencing (WGBS). We identified, in total, 2934 differentially methylated regions (DMRs) between RT_48h (resistant strain with high-temperature/humidity treatment for 48 h) and ST_48h (sensitive strain with high-temperature/humidity treatment for 48 h) under cytosine context (CG), which corresponded to 1230 DMR-related genes (DMGs), and the DMRs were primarily located in the gene body (exon and intron) region. Gene ontology (GO) and KEGG analysis showed that these DMGs were most significantly enriched in binding, cellular metabolic process, and RNA transport pathways. Moreover, 10 DMGs have been revealed to be involved in the heat-humidity stress response in the silkworm. The results of this study indicated that DNA methylation plays crucial roles in silkworm response to environmental stressors and provides important clues to identify key resistance genes in silkworm under high-temperature/humidity stress response.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

DNA methylation suppresses chitin degradation and promotes the wing development by inhibiting Bmara-mediated chitinase expression in the silkworm, Bombyx mori

BACKGROUND

DNA methylation, as an essential epigenetic modification found in mammals and plants, has been implicated to play an important role in insect reproduction. However, the functional role and the regulatory mechanism of DNA methylation during insect organ or tissue development are far from being clear.

RESULTS

Here, we found that DNA methylation inhibitor (5-aza-dC) treatment in newly molted pupae decreased the chitin content of pupal wing discs and adult wings and resulted in wing deformity of Bombyx mori. Transcriptome analysis revealed that the up-regulation of chitinase 10 (BmCHT10) gene might be related to the decrease of chitin content induced by 5-aza-dC treatment. Further, the luciferase activity assays demonstrated that DNA methylation suppressed the promoter activity of BmCHT10 by down-regulating the transcription factor, homeobox protein araucan (Bmara). Electrophoretic mobility shift assay, DNA pull-down and chromatin immunoprecipitation demonstrated that Bmara directly bound to the BmCHT10 promoter. Therefore, DNA methylation is involved in keeping the structural integrity of the silkworm wings from unwanted chitin degradation, as a consequence, it promotes the wing development of B. mori.

CONCLUSIONS

This study reveals that DNA methylation plays an important role in the wing development of B. mori. Our results support that the indirect transcriptional repression of a chitin degradation-related gene BmCHT10 by DNA methylation is necessary to keep the proper wing development in B. mori.

Comparative analysis of genome-wide DNA methylation in Neurospora

DNA methylation is an epigenetic mark that plays an important role in genetic regulation in eukaryotes. Major progress has been made in dissecting the molecular pathways that regulate DNA methylation. Yet, little is known about DNA methylation variation over evolutionary time. Here we present an investigation of the variation of DNA methylation and transposable element (TE) content in species of the filamentous ascomycetes Neurospora. We generated genome-wide DNA methylation data at single-base resolution, together with genomic TE content and gene expression data, of 10 individuals representing five closely related Neurospora species. We found that the methylation levels were low (ranging from 1.3% to 2.5%) and varied among the genomes in a species-specific way. Furthermore, we found that the TEs over 400 bp long were targeted by DNA methylation, and in all genomes, high methylation correlated with low GC, confirming a conserved link between DNA methylation and Repeat Induced Point (RIP) mutations in this group of fungi. Both TE content and DNA methylation pattern showed phylogenetic signal, and the species with the highest TE load (N. crassa) also exhibited the highest methylation level per TE. Our results suggest that DNA methylation is an evolvable trait and indicate that the genomes of Neurospora are shaped by an evolutionary arms race between TEs and host defence.

Genome Report: Whole Genome Sequence and Annotation of the Parasitoid Jewel Wasp Nasonia giraulti Laboratory Strain RV2X[u]

Jewel wasps in the genus of Nasonia are parasitoids with haplodiploidy sex determination, rapid development and are easy to culture in the laboratory. They are excellent models for insect genetics, genomics, epigenetics, development, and evolution. Nasonia vitripennis (Nv) and N. giraulti (Ng) are closely-related species that can be intercrossed, particularly after removal of the intracellular bacterium Wolbachia, which serve as a powerful tool to map and positionally clone morphological, behavioral, expression and methylation phenotypes. The Nv reference genome was assembled using Sanger, PacBio and Nanopore approaches and annotated with extensive RNA-seq data. In contrast, Ng genome is only available through low coverage resequencing. Therefore, de novo Ng assembly is in urgent need to advance this system. In this study, we report a high-quality Ng assembly using 10X Genomics linked-reads with 670X sequencing depth. The current assembly has a genome size of 259,040,977 bp in 3,160 scaffolds with 38.05% G-C and a 98.6% BUSCO completeness score. 97% of the RNA reads are perfectly aligned to the genome, indicating high quality in contiguity and completeness. A total of 14,777 genes are annotated in the Ng genome, and 72% of the annotated genes have a one-to-one ortholog in the Nv genome. We reported 5 million Ng-Nv SNPs which will facility mapping and population genomic studies in Nasonia In addition, 42 Ng-specific genes were identified by comparing with Nv genome and annotation. This is the first de novo assembly for this important species in the Nasonia model system, providing a useful new genomic toolkit.

Reprograming of epigenetic mechanisms controlling host insect immunity and development in response to egg-laying by a parasitoid wasp

Parasitoids are insects that use other insects as hosts. They sabotage host cellular and humoral defences to promote the survival of their offspring by injecting viruses and venoms along with their eggs. Many pathogens and parasites disrupt host epigenetic mechanisms to overcome immune system defences, and we hypothesized that parasitoids may use the same strategy. We used the ichneumon wasp Pimpla turionellae as a model idiobiont parasitoid to test this hypothesis, with pupae of the greater wax moth Galleria mellonella as the host. We found that parasitoid infestation involves the suppression of host immunity-related effector genes and the modulation of host genes involved in developmental hormone signalling. The transcriptional reprogramming of host genes following the injection of parasitoid eggs was associated with changes in host epigenetic mechanisms. The introduction of parasitoids resulted in a transient decrease in host global DNA methylation and the modulation of acetylation ratios for specific histones. Genes encoding regulators of histone acetylation and deacetylation were mostly downregulated in the parasitized pupae, suggesting that parasitoids can suppress host transcription. We also detected a strong parasitoid-specific effect on host microRNAs regulating gene expression at the post-transcriptional level. Our data therefore support the hypothesis that parasitoids may favour the survival of their offspring by interfering with host epigenetic mechanisms to suppress the immune system and disrupt development.

------Widespread conservation and lineage-specific diversification of genome-wide DNA methylation patterns across arthropods

Cytosine methylation is an ancient epigenetic modification yet its function and extent within genomes is highly variable across eukaryotes. In mammals, methylation controls transposable elements and regulates the promoters of genes. In insects, DNA methylation is generally restricted to a small subset of transcribed genes, with both intergenic regions and transposable elements (TEs) depleted of methylation. The evolutionary origin and the function of these methylation patterns are poorly understood. Here we characterise the evolution of DNA methylation across the arthropod phylum. While the common ancestor of the arthropods had low levels of TE methylation and did not methylate promoters, both of these functions have evolved independently in centipedes and mealybugs. In contrast, methylation of the exons of a subset of transcribed genes is ancestral and widely conserved across the phylum, but has been lost in specific lineages. A similar set of genes is methylated in all species that retained exon-enriched methylation. We show that these genes have characteristic patterns of expression correlating to broad transcription initiation sites and well-positioned nucleosomes, providing new insights into potential mechanisms driving methylation patterns over hundreds of millions of years.

Quantitative genetics of wing morphology in the parasitoid wasp Nasonia vitripennis: hosts increase sibling similarity

The central aim of evolutionary biology is to understand patterns of genetic variation between species and within populations. To quantify the genetic variation underlying intraspecific differences, estimating quantitative genetic parameters of traits is essential. In Pterygota, wing morphology is an important trait affecting flight ability. Moreover, gregarious parasitoids such as Nasonia vitripennis oviposit multiple eggs in the same host, and siblings thus share a common environment during their development. Here we estimate the genetic parameters of wing morphology in the outbred HVRx population of N. vitripennis, using a sire-dam model adapted to haplodiploids and disentangled additive genetic and host effects. The results show that the wing-size traits have low heritability (h² ~ 0.1), while most wing-shape traits have roughly twice the heritability compared with wing-size traits. However, the estimates increased to h² ~ 0.6 for wing-size traits when omitting the host effect from the statistical model, while no meaningful increases were observed for wing-shape traits. Overall, host effects contributed to ~50% of the variation in wing-size traits. This indicates that hosts have a large effect on wing-size traits, about fivefold more than genetics. Moreover, bivariate analyses were conducted to derive the genetic relationships among traits. Overall, we demonstrate the evolutionary potential for morphological traits in the N. vitripennis HVRx-outbred population, and report the host effects on wing morphology. Our findings can contribute to a further dissection of the genetics underlying wing morphology in N. vitripennis, with relevance for gregarious parasitoids and possibly other insects as well.

Distinct epigenomic and transcriptomic modifications associated with Wolbachia-mediated asexuality

Wolbachia are maternally transmitted intracellular bacteria that induce a range of pathogenic and fitness-altering effects on insect and nematode hosts. In parasitoid wasps of the genus Trichogramma, Wolbachia infection induces asexual production of females, thus increasing transmission of Wolbachia. It has been hypothesized that Wolbachia infection accompanies a modification of the host epigenome. However, to date, data on genome-wide epigenomic changes associated with Wolbachia are limited, and are often confounded by background genetic differences. Here, we took sexually reproducing Trichogramma free of Wolbachia and introgressed their genome into a Wolbachia-infected cytoplasm, converting them to Wolbachia-mediated asexuality. Wolbachia was then cured from replicates of these introgressed lines, allowing us to examine the genome-wide effects of wasps newly converted to asexual reproduction while controlling for genetic background. We thus identified gene expression and DNA methylation changes associated with Wolbachia-infection. We found no overlaps between differentially expressed genes and differentially methylated genes, indicating that Wolbachia-infection associated DNA methylation change does not directly modulate levels of gene expression. Furthermore, genes affected by these mechanisms exhibit distinct evolutionary histories. Genes differentially methylated due to the infection tended to be evolutionarily conserved. In contrast, differentially expressed genes were significantly more likely to be unique to the Trichogramma lineage, suggesting host-specific transcriptomic responses to infection. Nevertheless, we identified several novel aspects of Wolbachia-associated DNA methylation changes. Differentially methylated genes included those involved in oocyte development and chromosome segregation. Interestingly, Wolbachia-infection was associated with higher levels of DNA methylation. Additionally, Wolbachia infection reduced overall variability in gene expression, even after accounting for the effect of DNA methylation. We also identified specific cases where alternative exon usage was associated with DNA methylation changes due to Wolbachia infection. These results begin to reveal distinct genes and molecular pathways subject to Wolbachia induced epigenetic modification and/or host responses to Wolbachia-infection.

Wolbachia is an extremely common endosymbiotic infection of arthropods and nematodes. One of the reasons why Wolbachia can so successfully infect diverse species is the bacterium’s ability to profoundly alter the reproductive behavior of its host. It has been proposed that Wolbachia may modify host’s epigenetic programs to alter its reproductive behavior. However, it has been difficult to study how epigenetic programs change with Wolbachia infection, due to the confounding effects of genetic backgrounds. Here, we studied host transcriptome and epigenome changes associated with Wolbachia infection in a homogenous genetic background, by carrying out an innovative introgression scheme. By doing so, we show, for the first time, high-resolution molecular consequences of intracellular infection and offer insights into epigenetic and transcriptomic regulation of invertebrates.

Epigenetic effects of (nano)materials in environmental species - Cu case study in Enchytraeus crypticus

Chemical stressors can induce epigenomic changes, i.e., changes that are transferred to the next generation, even when the stressor is removed. Literature on chemical induced epigenetic effects in environmental species is scarce. We here provide the first results on epigenetic effects caused by nanomaterials with an environmental OECD standard soil model species Enchytraeus crypticus species. We assessed the epigenetic potential in terms of global DNA methylation, gene-specific methylation via bisulfite sequencing and MS-HRM (Methylation Sensitive - High Resolution Melting), and gene expression qPCR for genes involved in DNA methylation, histone modifications, non-coding RNA and stress response mechanisms). We have exposed E. crypticus in a multigenerational (MG) test design to Cu (copper oxide nanomaterials (CuO NMs) and copper salt (CuCl₂)). To link possible epigenetic effects to population changes, we used exposure concentrations (ECx) that caused a 10% and 50% reduction in the reproductive output (10% and 50% are the standards for regulatory Risk Assessment), the organisms were exposed for five consecutive generations (F1-F5) plus two generations after transferring to clean media (F5-F7), 7 generations in a total of 224 days. Results showed that MG exposure to Cu increased global DNA methylation and corresponded with phenotypic effects (reproduction). Gene expression analyses showed changes in the epigenetic, stress and detoxification gene targets, depending on the generation and Cu form, also occurring in post-exposure generations, hence indicative of transgenerational effects. There were in general clear differences between organisms exposed to different Cu-forms, hence indicate nanoparticulate-specific effects.

The Roles of DNA Methyltransferases 1 (DNMT1) in Regulating Sexual Dimorphism in the Cotton Mealybug, Phenacoccus solenopsis

The cotton mealybug, Phenacoccus solenopsis, is an invasive pest that can cause massive damage to many host plants of agricultural importance. P. solenopsis is highly polyphagous, and shows extreme sexual dimorphism between males and females. The functions of DNA methyltransferase (DNMT) enzymes in the cotton mealybug have not been well studied. Here, we carried out an investigation of DNMTs in cotton mealybug to study their roles in sexual dimorphism. We found that the cotton mealybug has two copies of PsDnmt1, but Dnmt3 is absent. We then amplified the full-length cDNAs of PsDnmt1A (2,225 bp) and PsDnmt1B (2,862 bp) using rapid amplification cDNA ends (RACE). Quantitative reverse transcriptase PCR shows that both PsDnmt1A and PsDnmt1B are highly expressed in adult males, while the expression of PsDnmt1B is 30-fold higher in gravid females than in virgin females. We knocked down PsDnmt1A and PsDnmt1B with small interfering RNAs (siRNAs), and both genes were successfully down-regulated after 24 h or 72 h in adult females and pupa (t-test, p < 0.05). Down-regulating the expression of these two DNMT genes led to offspring lethality and abnormal body color in adult females. Furthermore, the silencing of PsDnmt1B induced abnormal wing development in emerged adult males. Our results provide evidence that PsDnmt1 plays a crucial role in regulating sexual dimorphism in the cotton mealybug.

DNA Methylation Is Correlated with Gene Expression during Diapause Termination of Early Embryonic Development in the Silkworm (Bombyx mori)

DNA modification is a naturally occurring DNA modification in prokaryotic and eukaryotic organisms and is involved in several biological processes. Although genome-wide methylation has been studied in many insects, the understanding of global and genomic DNA methylation during insect early embryonic development, is lacking especially for insect diapause. In this study, we analyzed the relationship between DNA methylomes and transcriptomes in diapause-destined eggs compared to diapause-terminated eggs in the silkworm, Bombyx mori (B. mori). The results revealed that methylation was sparse in this species, as previously reported. Moreover, methylation levels in diapause-terminated eggs (HCl-treated) were 0.05% higher than in non-treated eggs, mainly due to the contribution of CG methylation sites. Methylation tends to occur in the coding sequences and promoter regions, especially at transcription initiation sites and short interspersed elements. Additionally, 364 methylome- and transcriptome-associated genes were identified, which showed significant differences in methylation and expression levels in diapause-destined eggs when compared with diapause-terminated eggs, and 74% of methylome and transcriptome associated genes showed both hypermethylation and elevated expression. Most importantly, Kyoto Encyclopaedia of Genes and Genomes (KEGG) analyses showed that methylation may be positively associated with Bombyx mori embryonic development, by regulating cell differentiation, metabolism, apoptosis pathways and phosphorylation. Through analyzing the G2/M phase-specific E3 ubiquitin-protein ligase (G2E3), we speculate that methylation may affect embryo diapause by regulating the cell cycle in Bombyx mori. These findings will help unravel potential linkages between DNA methylation and gene expression during early insect embryonic development and insect diapause.

Sex biased expression and co-expression networks in development, using the hymenopteran Nasonia vitripennis

Sexual dimorphism requires regulation of gene expression in developing organisms. These developmental differences are caused by differential expression of genes and isoforms. The effect of expressing a gene is also influenced by which other genes are simultaneously expressed (functional interactions). However, few studies have described how these processes change across development. We compare the dynamics of differential expression, isoform switching and functional interactions in the sexual development of the model parasitoid wasp Nasonia vitripennis, a system that permits genome wide analysis of sex bias from early embryos to adults. We find relatively little sex-bias in embryos and larvae at the gene level, but several sub-networks show sex-biased functional interactions in early developmental stages. These networks provide new candidates for hymenopteran sex determination, including histone modification. In contrast, sex-bias in pupae and adults is driven by the differential expression of genes. We observe sex-biased isoform switching consistently across development, but mostly in genes that are already differentially expressed. Finally, we discover that sex-biased networks are enriched by genes specific to the Nasonia clade, and that those genes possess the topological properties of key regulators. These findings suggest that regulators in sex-biased networks evolve more rapidly than regulators of other developmental networks.

Evolution of DNA Methylome Diversity in Eukaryotes

Cytosine DNA methylation (5mC) is a widespread base modification in eukaryotic genomes with critical roles in transcriptional regulation. In recent years, our understanding of 5mC has changed because of advances in 5mC detection techniques that allow mapping of this mark on the whole genome scale. Profiling DNA methylomes from organisms across the eukaryotic tree of life has reshaped our views on the evolution of 5mC. In this review, we explore the macroevolution of 5mC in major eukaryotic groups, and then focus on recent advances made in animals. Genomic 5mC patterns as well as the mechanisms of 5mC deposition tend to be evolutionary labile across large phylogenetic distances; however, some common patterns are starting to emerge. Within the animal kingdom, 5mC diversity has proven to be much greater than anticipated. For example, a previously held common view that genome hypermethylation is a trait exclusive to vertebrates has recently been challenged. Also, data from genome-wide studies are starting to yield insights into the potential roles of 5mC in invertebrate cis regulation. Here we provide an evolutionary perspective of both the well-known and enigmatic roles of 5mC across the eukaryotic tree of life.

Analysis of γ-Tocopherol methyl transferase3 promoter activity and study of methylation patterns of the promoter and its gene body

Soybeans are known for its good source of protein (40%), oil (20%) and also serve as a source of nutraceutical compounds including tocopherols (toc). To know the molecular basis of differential α-toc accumulation in two contrasting soybean genotypes: DS74 (low α-toc - 1.36 μg/g and total-toc -29.72 μg/g) and Bragg (high α-toc - 10.48 μg/g and total-toc 178.91 μg/g), the analysis of γ-TMT3 promoter activity and its methylation patterns were carried out. The sequencing results revealed nucleotide variation between Bragg:γ-TMT3-P and DS74:γ-TMT3-P, however none of the variations were found in core-promoter region or in cis-elements. The histochemical GUS assay revealed higher promoter activity of Bragg:γ-TMT3-P than that of DS74:γ-TMT3-P and correlated with significantly higher and lower (P < 0.05) expression of γ-TMT3 gene respectively. To know the molecular basis of differential accumulation of α-toc in these contrasting soybean genotypes, the DNA methylation pattern of γ-TMT3 gene body and its promoter was studied in both varieties. The results showed higher percentage (62.5%) of methylation in DS74:γ-TMT3-P than in Bragg:γ-TMT3-P (50%). Out of all the methylation sites in the promoter region, one of methylation site was found at CAAT box (-190 bp) of DS74:γ-TMT3-P. Further gene body methylation patterns revealed lowest % (40%) of CG methylation in DS74:γ-TMT3 gene as compared to Bragg:γ-TMT3 (64.2%). Thus our study revealed that, expression of γ-TMT3 gene was influenced by its promoter activity and methylation patterns in cis-elements of γ-TMT3 promoter and gene body. This study will help us to understand the possible role of methylation and promoter activity in determining the α-toc content in soybean seeds.

DNA methylation is maintained with high fidelity in the honey bee germline and exhibits global non-functional fluctuations during somatic development

BACKGROUND

DNA methylation of active genes, also known as gene body methylation, is found in many animal and plant genomes. Despite this, the transcriptional and developmental role of such methylation remains poorly understood. Here, we explore the dynamic range of DNA methylation in honey bee, a model organism for gene body methylation.

RESULTS

Our data show that CG methylation in gene bodies globally fluctuates during honey bee development. However, these changes cause no gene expression alterations. Intriguingly, despite the global alterations, tissue-specific CG methylation patterns of complete genes or exons are rare, implying robust maintenance of genic methylation during development. Additionally, we show that CG methylation maintenance fluctuates in somatic cells, while reaching maximum fidelity in sperm cells. Finally, unlike universally present CG methylation, we discovered non-CG methylation specifically in bee heads that resembles such methylation in mammalian brain tissue.

CONCLUSIONS

Based on these results, we propose that gene body CG methylation can oscillate during development if it is kept to a level adequate to preserve function. Additionally, our data suggest that heightened non-CG methylation is a conserved regulator of animal nervous systems.

Methylation and gene expression differences between reproductive and sterile bumblebee workers

Phenotypic plasticity is the production of multiple phenotypes from a single genome and is notably observed in social insects. Multiple epigenetic mechanisms have been associated with social insect plasticity, with DNA methylation being explored to the greatest extent. DNA methylation is thought to play a role in caste determination in Apis mellifera, and other social insects, but there is limited knowledge on its role in other bee species. In this study, we analyzed whole genome bisulfite sequencing and RNA-seq data sets from head tissue of reproductive and sterile castes of the eusocial bumblebee Bombus terrestris. We found that genome-wide methylation in B. terrestris is similar to other holometabolous insects and does not differ between reproductive castes. We did, however, find differentially methylated genes between castes, which are enriched for multiple biological processes including reproduction. However, we found no relationship between differential methylation and differential gene expression or differential exon usage between castes. Our results also indicate high intercolony variation in methylation. These findings suggest that methylation is associated with caste differences but may serve an alternate function, other than direct caste determination in this species. This study provides the first insights into the nature of a bumblebee caste-specific methylome as well as its interaction with gene expression and caste-specific alternative splicing, providing greater understanding of the role of methylation in phenotypic plasticity within social bee species. Future experimental work is needed to determine the function of methylation and other epigenetic mechanisms in insects.

Implication for DNA methylation involved in the host transfer of diamondback moth, Plutella xylostella (L.)

DNA methylation exerts extensive impacts on gene expression of various living organisms exposed to environmental variation. However, little is known whether DNA methylation is involved in the host transfer of diamondback moth, Plutella xylostella (L.), a worldwide destructive pest of crucifers. In this study, we found that P. xylostella genome exhibited a relatively low level of DNA methylation on the basis of the CpG O/E prediction and experimental validation. A significant positive linear correlation was observed between the stage-specific expressions of PxDNMT1 and DNA methylation levels (5mC content). Particularly, high levels of DNA methylation and gene expression of PxDNMT1 were observed in eggs and mature females of P. xylostella. After host transfer of P. xylostella from Raphanus sativus to Arabidopsis thaliana, we identified some potential genomic loci that might have changed methylation levels. Using the method of fluorescence-labeled methylation-sensitive amplified polymorphism (F-MSAP), we also found the corresponding genes primarily involved in neural system and signaling. The expressions of six candidate genes were verified by qRT-PCR. One of the genes, Px009600, might be regulated by a DNA methylation-mediated mechanism in response to host transfer. Our study provides evidence for a functional system of DNA methylation in P. xylostella and its possible role in adaptation during host transfer. Further studies should examine methylation as responsive factors to different host plants and environmental cues in insect pests.

Convergent evolution of a vertebrate-like methylome in a marine sponge

Vertebrates have highly methylated genomes at CpG positions whereas invertebrates have sparsely methylated genomes. This increase in methylation content is considered a major regulatory innovation of vertebrate genomes. However, here we report that a marine sponge, proposed as the sister group to the rest of animals, has a highly methylated genome. Despite major differences in genome size and architecture, we find similarities between the independent acquisitions of the hypermethylated state. Both lineages show genome wide CpG depletion, conserved strong transcription factor methyl-sensitivity, and developmental methylation dynamics at 5-hydroxymethylcytosine enriched regions. Together, our findings trace back patterns associated with DNA methylation in vertebrates to the early steps of animal evolution. Thus, the sponge methylome challenges prior hypotheses concerning the uniqueness of vertebrate genome hypermethylation and its implications for regulatory complexity.

Sex-specific changes in the aphid DNA methylation landscape

Aphids present an ideal system to study epigenetics as they can produce diverse, but genetically identical, morphs in response to environmental stimuli. Here, using whole genome bisulphite sequencing and transcriptome sequencing of the green peach aphid (Myzus persicae), we present the first detailed analysis of cytosine methylation in an aphid and investigate differences in the methylation and transcriptional landscapes of male and asexual female morphs. We found that methylation primarily occurs in a CG dinucleotide (CpG) context and that exons are highly enriched for methylated CpGs, particularly at the 3' end of genes. Methylation is positively associated with gene expression, and methylated genes are more stably expressed than unmethylated genes. Male and asexual female morphs have distinct methylation profiles. Strikingly, these profiles are divergent between the sex chromosome and the autosomes; autosomal genes are hypomethylated in males compared to asexual females, whereas genes belonging to the sex chromosome, which is haploid in males, are hypermethylated. Overall, we found correlated changes in methylation and gene expression between males and asexual females, and this correlation was particularly strong for genes located on the sex chromosome. Our results suggest that differential methylation of sex-biased genes plays a role in aphid sexual differentiation.

DNA Methylation: Shared and Divergent Features across Eukaryotes

Chemical modification of nucleotide bases in DNA provides one mechanism for conveying information in addition to the genetic code. 5-methylcytosine (5mC) represents the most common chemically modified base in eukaryotic genomes. Sometimes referred to simply as DNA methylation, in eukaryotes 5mC is most prevalent at CpG dinucleotides and is frequently associated with transcriptional repression of transposable elements. However, 5mC levels and distributions are variable across phylogenies, and emerging evidence suggests that the functions of DNA methylation may be more diverse and complex than was previously appreciated. We summarize the current understanding of DNA methylation profiles and functions in different eukaryotic lineages.

Detecting differential DNA methylation from sequencing of bisulfite converted DNA of diverse species

DNA methylation is one of the most extensively studied epigenetic modifications of genomic DNA. In recent years, sequencing of bisulfite-converted DNA, particularly via next-generation sequencing technologies, has become a widely popular method to study DNA methylation. This method can be readily applied to a variety of species, dramatically expanding the scope of DNA methylation studies beyond the traditionally studied human and mouse systems. In parallel to the increasing wealth of genomic methylation profiles, many statistical tools have been developed to detect differentially methylated loci (DMLs) or differentially methylated regions (DMRs) between biological conditions. We discuss and summarize several key properties of currently available tools to detect DMLs and DMRs from sequencing of bisulfite-converted DNA. However, the majority of the statistical tools developed for DML/DMR analyses have been validated using only mammalian data sets, and less priority has been placed on the analyses of invertebrate or plant DNA methylation data. We demonstrate that genomic methylation profiles of non-mammalian species are often highly distinct from those of mammalian species using examples of honey bees and humans. We then discuss how such differences in data properties may affect statistical analyses. Based on these differences, we provide three specific recommendations to improve the power and accuracy of DML and DMR analyses of invertebrate data when using currently available statistical tools. These considerations should facilitate systematic and robust analyses of DNA methylation from diverse species, thus advancing our understanding of DNA methylation.

Validating the Demethylating Effects of 5-aza-2'-deoxycytidine in Insects Requires a Whole-Genome Approach

We previously demonstrated that treatment with the demethylating agent 5-aza-2'-deoxycytidine (5-aza-dC) alters the offspring sex ratios produced by females of the parasitoid wasp Nasonia vitripennis. Females allocate offspring sex ratio in line with local mate competition theory, producing more or less female-biased sex ratios as the number of other females laying eggs on a patch varies, thereby reducing competition among their sons for mates. Interestingly, treatment with 5-aza-dC did not ablate the facultative sex allocation response. Instead, sex ratios became less female biased, a shift in the direction of the optimum sex ratio for paternally inherited alleles according to genomic conflict theory. This was the first (albeit indirect) experimental evidence for genomic conflict over sex allocation. In their comment, Ellers and colleagues assayed the effects of 5-aza-dC on DNA methylation in 10 Nasonia genes, finding no evidence of demethylation in these 10 genes, from which they conclude that 5-aza-dC has no demethylating capability in N. vitripennis. Quantifying the efficacy of 5-aza-dC in terms of demethylation is indeed crucial to in-depth interpretation of studies using 5-aza-dC to link phenotypes to epigenetic regulation. Here we outline the mode of action of 5-aza-dC and demonstrate that determining the efficacy of 5-aza-dC in insect systems requires a whole-genome approach.

The Importance of Validating the Demethylating Effect of 5-aza-2'-deoxycytidine in Model Species (A Comment on Cook et al., "DNA Methylation and Sex Allocation in the Parasitoid Wasp Nasonia vitripennis")

The use of DNA demethylating agents has been popular in epigenetic studies. Recently, Cook and colleagues, in a 2015 American Naturalist article, claimed an effect of 5-aza-2'-deoxycytidine (5-aza-dC) on the sex ratio of a parasitoid wasp without verifying its effect on DNA methylation. We repeated the 5-aza-dC feeding treatment to test its effectiveness. We used bisulfite amplicon sequencing of 10 genes that either were heavily methylated, previously showed a response to 5-aza-dC, or were suggested to regulate fatty acid synthesis epigenetically, and we demonstrate that wasps fed 5-aza-dC did not show reduced DNA methylation at these loci. Therefore, the conclusion that demethylation shifts sex ratios upward needs reconsideration.

The effects of the neonicotinoid imidacloprid on gene expression and DNA methylation in the buff-tailed bumblebee Bombus terrestris

Neonicotinoids are effective insecticides used on many important arable and horticultural crops. They are nicotinic acetylcholine receptor agonists which disrupt the function of insect neurons and cause paralysis and death. In addition to direct mortality, there are numerous sublethal effects of low doses of neonicotinoids on bees. We hypothesize that some of these large array of effects could be a consequence of epigenetic changes in bees induced by neonicotinoids. We compared whole methylome (BS-seq) and RNA-seq libraries of the brains of buff-tailed bumblebee Bombus terrestris workers exposed to field-realistic doses of the neonicotinoid imidacloprid to libraries from control workers. We found numerous genes which show differential expression between neonicotinoid-treated bees and control bees, but no differentially methylated cytosines in any context. We found CpG methylation to be focused mainly in exons and associated with highly expressed genes. We discuss the implications of our results for future legislation.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.