Genome-wide association between DNA methylation and alternative splicing in an invertebrate

Epigenetic Regulation of Mammalian Cardiomyocyte Development

The heart is the first organ formed during mammalian development and functions to distribute nutrients and oxygen to other parts of the developing embryo. Cardiomyocytes are the major cell types of the heart and provide both structural support and contractile function to the heart. The successful differentiation of cardiomyocytes during early development is under tight regulation by physical and molecular factors. We have reviewed current studies on epigenetic factors critical for cardiomyocyte differentiation, including DNA methylation, histone modifications, chromatin remodelers, and noncoding RNAs. This review also provides comprehensive details on structural and morphological changes associated with the differentiation of fetal and postnatal cardiomyocytes and highlights their differences. A holistic understanding of all aspects of cardiomyocyte development is critical for the successful in vitro differentiation of cardiomyocytes for therapeutic purposes.

Lack of paternal silencing and ecotype-specific expression in head and body lice hybrids

Paternal genome elimination (PGE) is a non-Mendelian inheritance system, described in numerous arthropod species, in which males develop from fertilized eggs, but their paternally inherited chromosomes are eliminated before or during spermatogenesis. Therefore, PGE males only transmit their maternally inherited set of chromosomes to their offspring. In addition to the elimination of paternal chromosomes, diverse PGE species have also repeatedly evolved the transcriptional silencing of the paternal genome, making males effectively haploid. However, it is unclear if this paternal chromosome silencing is mechanistically linked to the chromosome elimination or has evolved at a later stage, and if so, what drives the haploidization of males under PGE. In order to understand these questions, here we study the human louse, Pediculus humanus, which represents an ideal model system, as it appears to be the only instance of PGE where males eliminate, but not silence their paternal chromosomes, although the latter remains to be shown conclusively. In this study, we analyzed parent-of-origin allele-specific expression patterns in male offspring of crosses between head and body lice ecotypes. We show that hybrid adult males of P. humanus display biparental gene expression, which constitutes the first case of a species with PGE in which genetic activity of paternal chromosomes in the soma is not affected by embryonic silencing or (partial or complete) elimination. We did however also identify a small number of maternally biased genes (potentially imprinted genes), which may be involved in the elimination of paternal chromosomes during spermatogenesis. Finally, we have identified genes that show ecotype-specific expression bias. Given the low genetic diversity between ecotypes, this is suggestive for a role of epigenetic processes in ecotype differences.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Ancestors' Gift: Parental Early Exposure to the Environmentally Realistic Pesticide Mixture Drives Offspring Phenotype in a Larger Extent Than Direct Exposure in the Pacific Oyster, Crassostrea gigas

Marine organisms are threatened by the presence of pesticides in coastal waters. Among them, the Pacific oyster is one of the most studied invertebrates in marine ecotoxicology where numerous studies highlighted the multiscale impacts of pesticides. In the past few years, a growing body of literature has reported the epigenetic outcomes of xenobiotics. Because DNA methylation is an epigenetic mark implicated in organism development and is meiotically heritable, it raises the question of the multigenerational implications of xenobiotic-induced epigenetic alterations. Therefore, we performed a multigenerational exposure to an environmentally relevant mixture of 18 pesticides (nominal sum concentration: 2.85 μg·L-1) during embryo-larval stages (0-48 hpf) of a second generation (F1) for which parents where already exposed or not in F0. Gene expression, DNA methylation, and physiological end points were assessed throughout the life cycle of individuals. Overall, the multigenerational effect has a greater influence on the phenotype than the exposure itself. Thus, multigenerational phenotypic effects were observed: individuals descending from exposed parents exhibited lower epinephrine-induced metamorphosis and field survival rates. At the molecular level, RNA-seq and Methyl-seq data analyses performed in gastrula embryos and metamorphosis-competent pediveliger (MCP) larvae revealed a clear F0 treatment-dependent discrimination. Some genes implicated into shell secretion and immunity exhibited F1:F0 treatment interaction patterns (e.g., Calm and Myd88). Those results suggest that low chronic environmental pesticide contamination can alter organisms beyond the individual scale level and have long-term adaptive implications.

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.

The Evolutionary Complexities of DNA Methylation in Animals: From Plasticity to Genetic Evolution

The importance of DNA methylation in plastic responses to environmental change and evolutionary dynamics is increasingly recognized. Here, we provide a Perspective piece on the diverse roles of DNA methylation on broad evolutionary timescales, including (i) short-term transient acclimation, (ii) stable phenotypic evolution, and (iii) genomic evolution. We show that epigenetic responses vary along a continuum, ranging from short-term acclimatory responses in variable environments within a generation to long-term modifications in populations and species. DNA methylation thus unlocks additional potential for organisms to rapidly acclimate to their environment over short timeframes. If these changes affect fitness, they can circumvent the need for adaptive changes at the genome level. However, methylation has a complex reciprocal relationship with genetic variation as it can be genetically controlled, yet it can also induce point mutations and contribute to genomic evolution. When habitats remain constant over many generations, or populations are separated across habitats, initially plastic phenotypes can become hardwired through epigenetically facilitated mutagenesis. It remains unclear under what circumstances plasticity contributes to evolutionary outcomes, and when plastic changes will become permanently encoded into genotype. We highlight how studies investigating the evolution of epigenetic plasticity need to carefully consider how plasticity in methylation state could evolve among different evolutionary scenarios, the possible phenotypic outcomes, its effects on genomic evolution, and the proximate energetic and ultimate fitness costs of methylation. We argue that accumulating evidence suggests that DNA methylation can contribute toward evolution on various timescales, spanning a continuum from acclimatory plasticity to genomic evolution.

The PWWP domain and the evolution of unique DNA methylation toolkits in Hymenoptera

DNMT3 in Hymenoptera has a unique duplication of the essential PWWP domain. Using GST-tagged PWWP fusion proteins and histone arrays we show that these domains have gained new properties and represent the first case of PWWP domains binding to H3K27 chromatin modifications, including H3K27me3, a key modification that is important during development. Phylogenetic analyses of 107 genomes indicate that the duplicated PWWP domains separated into two sister clades, and their distinct binding capacities are supported by 3D modeling. Other features of this unique DNA methylation system include variable copies, losses, and duplications of DNMT1 and DNMT3, and combinatorial generations of DNMT3 isoforms including variants missing the catalytic domain. Some of these losses and duplications of are found only in parasitic wasps. We discuss our findings in the context of the crosstalk between DNA methylation and histone methylation, and the expanded potential of epigenomic modifications in Hymenoptera to drive evolutionary novelties.

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.

Relationships between phenotypic plasticity and epigenetic variation in two Caribbean Acropora corals

The plastic ability for a range of phenotypes to be exhibited by the same genotype allows organisms to respond to environmental variation and may modulate fitness in novel environments. Differing capacities for phenotypic plasticity within a population, apparent as genotype by environment interactions (GxE), can therefore have both ecological and evolutionary implications. Epigenetic gene regulation alters gene function in response to environmental cues without changes to the underlying genetic sequence and likely mediates phenotypic variation. DNA methylation is currently the most well described epigenetic mechanism and is related to transcriptional homeostasis in invertebrates. However, evidence quantitatively linking variation in DNA methylation with that of phenotype is lacking in some taxa, including reef-building corals. In this study, spatial and seasonal environmental variation in Bonaire, Caribbean Netherlands was utilized to assess relationships between physiology and DNA methylation profiles within genetic clones across different genotypes of Acropora cervicornis and A. palmata corals. The physiology of both species was highly influenced by environmental variation compared to the effect of genotype. GxE effects on phenotype were only apparent in A. cervicornis. DNA methylation in both species differed between genotypes and seasons and epigenetic variation was significantly related to coral physiological metrics. Furthermore, plastic shifts in physiology across seasons were significantly positively correlated with shifts in DNA methylation profiles in both species. These results highlight the dynamic influence of environmental conditions and genetic constraints on the physiology of two important Caribbean coral species. Additionally, this study provides quantitative support for the role of epigenetic DNA methylation in mediating phenotypic plasticity in invertebrates.

Third-Generation Sequencing Reveals the Adaptive Role of the Epigenome in Three Deep-Sea Polychaetes

The roles of DNA methylation in invertebrates are poorly characterized, and critical data are missing for the phylum Annelida. We fill this knowledge gap by conducting the first genome-wide survey of DNA methylation in the deep-sea polychaetes dominant in deep-sea vents and seeps: Paraescarpia echinospica, Ridgeia piscesae, and Paralvinella palmiformis. DNA methylation calls were inferred from Oxford Nanopore sequencing after assembling high-quality genomes of these animals. The genomes of these worms encode all the key enzymes of the DNA methylation metabolism and possess a mosaic methylome similar to that of other invertebrates. Transcriptomic data of these polychaetes support the hypotheses that gene body methylation strengthens the expression of housekeeping genes and that promoter methylation acts as a silencing mechanism but not the hypothesis that DNA methylation suppresses the activity of transposable elements. The conserved epigenetic profiles of genes responsible for maintaining homeostasis under extreme hydrostatic pressure suggest DNA methylation plays an important adaptive role in these worms.

Associations between DNA methylation and gene regulation depend on chromatin accessibility during transgenerational plasticity

BACKGROUND

Epigenetic processes are proposed to be a mechanism regulating gene expression during phenotypic plasticity. However, environmentally induced changes in DNA methylation exhibit little-to-no association with differential gene expression in metazoans at a transcriptome-wide level. It remains unexplored whether associations between environmentally induced differential methylation and expression are contingent upon other epigenomic processes such as chromatin accessibility. We quantified methylation and gene expression in larvae of the purple sea urchin Strongylocentrotus purpuratus exposed to different ecologically relevant conditions during gametogenesis (maternal conditioning) and modeled changes in gene expression and splicing resulting from maternal conditioning as functions of differential methylation, incorporating covariates for genomic features and chromatin accessibility. We detected significant interactions between differential methylation, chromatin accessibility, and genic feature type associated with differential expression and splicing.

RESULTS

Differential gene body methylation had significantly stronger effects on expression among genes with poorly accessible transcriptional start sites while baseline transcript abundance influenced the direction of this effect. Transcriptional responses to maternal conditioning were 4-13 × more likely when accounting for interactions between methylation and chromatin accessibility, demonstrating that the relationship between differential methylation and gene regulation is partially explained by chromatin state.

CONCLUSIONS

DNA methylation likely possesses multiple associations with gene regulation during transgenerational plasticity in S. purpuratus and potentially other metazoans, but its effects are dependent on chromatin accessibility and underlying genic features.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

DNA methylation is associated with codon degeneracy in a species of bumblebee

Social insects display extreme phenotypic differences between sexes and castes even though the underlying genome can be almost identical. Epigenetic processes have been proposed as a possible mechanism for mediating these phenotypic differences. Using whole genome bisulfite sequencing of queens, males, and reproductive female workers we have characterised the sex- and caste-specific methylome of the bumblebee Bombus terrestris. We have identified a potential role for DNA methylation in histone modification processes which may influence sex and caste phenotypic differences. We also find differentially methylated genes generally show low levels of DNA methylation which may suggest a separate function for lowly methylated genes in mediating transcriptional plasticity, unlike highly methylated genes which are usually involved in housekeeping functions. We also examined the relationship between the underlying genome and the methylome using whole genome re-sequencing of the same queens and males. We find DNA methylation is enriched at zero-fold degenerate sites. We suggest DNA methylation may be acting as a targeted mutagen at these sites, providing substrate for selection via non-synonymous changes in the underlying genome. However, we did not see any relationship between DNA methylation and rates of positive selection in our samples. In order to fully assess a possible role for DNA methylation in adaptive processes a specifically designed study using natural population data is needed.

DNA methylation profile in CpG-depleted regions uncovers a high-risk subtype of early-stage colorectal cancer

BACKGROUND

The current risk stratification system defined by clinicopathological features does not identify the risk of recurrence in early-stage (stage I-II) colorectal cancer (CRC) with sufficient accuracy. We aimed to investigate whether DNA methylation could serve as a novel biomarker for predicting prognosis in early-stage CRC patients.

METHODS

We analyzed the genome-wide methylation status of CpG loci using Infinium MethylationEPIC array run on primary tumor tissues and normal mucosa of early-stage CRC patients to identify potential methylation markers for prognosis. The machine-learning approach was applied to construct a DNA methylation-based prognostic classifier for early-stage CRC (MePEC) using the 4 gene methylation markers FAT3, KAZN, TLE4, and DUSP3. The prognostic value of the classifier was evaluated in 2 independent cohorts (n = 438 and 359, respectively).

RESULTS

The comprehensive analysis identified an epigenetic subtype with high risk of recurrence based on a group of CpG loci in the CpG-depleted region. In multivariable analysis, the MePEC classifier was independently and statistically significantly associated with time to recurrence in validation cohort 1 (hazard ratio = 2.35, 95% confidence interval = 1.47 to 3.76, P < .001) and cohort 2 (hazard ratio = 3.20, 95% confidence interval = 1.92 to 5.33, P < .001). All results were further confirmed after each cohort was stratified by clinicopathological variables and molecular subtypes.

CONCLUSIONS

We demonstrated the prognostic statistical significance of a DNA methylation profile in the CpG-depleted region, which may serve as a valuable source for tumor biomarkers. MePEC could identify an epigenetic subtype with high risk of recurrence and improve the prognostic accuracy of current clinical variables in early-stage CRC.

Epigenetic regulations as drivers of insecticide resistance and resilience to climate change in arthropod pests

Arthropod pests are remarkably capable of rapidly adapting to novel forms of environmental stress, including insecticides and climate change. The dynamic interplay between epigenetics and genetics explains the largely unexplored reality underlying rapid climatic adaptation and the development of insecticide resistance in insects. Epigenetic regulation modulates gene expression by methylating DNA and acetylating histones that play an essential role in governing insecticide resistance and adaptation to climate change. This review summarises and discusses the significance of recent advances in epigenetic regulation that facilitate phenotypic plasticity in insects and their symbiotic microbes to cope with selection pressure implied by extensive insecticide applications and climate change. We also discuss how epigenetic changes are passed on to multiple generations through sexual recombination, which remains enigmatic. Finally, we explain how these epigenetic signatures can be utilized to manage insecticide resistance and pest resilience to climate change in Anthropocene.

Adverse maternal environment affects hippocampal HTR2c variant expression and epigenetic characteristics in mouse offspring

BACKGROUND

An adverse maternal environment (AME) predisposes progeny towards cognitive impairment in humans and mice. Cognitive impairment associates with hippocampal dysfunction. An important regulator of hippocampal function is the hippocampal serotonergic system. Dysregulation of hippocampal serotonin receptor 2c (HTR2c) expression is linked with cognitive impairment. HTR2c contains multiple mRNA variants and isoforms that are epigenetically regulated including DNA methylation, histone modifications, and small nucleolar RNA MBII-52. We tested the hypotheses that AME increases HTR2c variant expression and alters epigenetic modifications along the HTR2c gene locus.

METHODS

We create an AME through maternal Western diet and prenatal environmental stress in the mouse. We analyzed hippocampal HTR2c and variants' expression, DNA methylation and histone modifications along the gene locus, and MBII-52 levels in postnatal day 21 offspring.

RESULTS

AME significantly increased the expressions of total HTR2c and full-length variants (V201 and V202) concurrently with an altered epigenetic profile along the HTR2c gene locus in male offspring hippocampi. Moreover, increased full-length variants' expression in AME males was in line with increased MBII-52 levels.

CONCLUSIONS

AME affects male offspring hippocampal expression of HTR2c and full-length variants via epigenetic mechanisms. Altered hippocampal HTR2c expression may contribute to cognitive impairment seen in adult males in this model.

IMPACT

The key message of our article is that an adverse maternal environment increases expression of total HTR2c mRNA and protein, alters proportions of HTR2c mRNA variants, and impacts HTR2c epigenetic modifications in male offspring hippocampi relative to controls. Our findings add to the literature by providing the first report of altered HTR2c mRNA variant expression in association with altered epigenetic modifications in the hippocampus of offspring mice exposed to an adverse maternal environment. Our findings suggest that an adverse maternal environment affects the expression of genes previously determined to regulate cognitive function through an epigenetic mechanism in a sex-specific manner.

Adapting to climate with limited genetic diversity: Nucleotide, DNA methylation and microbiome variation among populations of the social spider Stegodyphus dumicola

Understanding the role of genetic and nongenetic variants in modulating phenotypes is central to our knowledge of adaptive responses to local conditions and environmental change, particularly in species with such low population genetic diversity that it is likely to limit their evolutionary potential. A first step towards uncovering the molecular mechanisms underlying population-specific responses to the environment is to carry out environmental association studies. We associated climatic variation with genetic, epigenetic and microbiome variation in populations of a social spider with extremely low standing genetic diversity. We identified genetic variants that are associated strongly with environmental variation, particularly with average temperature, a pattern consistent with local adaptation. Variation in DNA methylation in many genes was strongly correlated with a wide set of climate parameters, thereby revealing a different pattern of associations than that of genetic variants, which show strong correlations to a more restricted range of climate parameters. DNA methylation levels were largely independent of cis-genetic variation and of overall genetic population structure, suggesting that DNA methylation can work as an independent mechanism. Microbiome composition also correlated with environmental variation, but most strong associations were with precipitation-related climatic factors. Our results suggest a role for both genetic and nongenetic mechanisms in shaping phenotypic responses to local environments.

Comparative Methylome Analysis Reveals Epigenetic Signatures Associated with Growth and Shell Color in the Pacific Oyster, Crassostrea gigas

Fast growth is one of the most important breeding goals for all economic species such as the Pacific oyster (Crassostrea gigas), an aquaculture mollusk with top global production. Although the genetic basis and molecular mechanisms of growth-related traits have been widely investigated in the oyster, the role of DNA methylation involved in growth regulation remains largely unexplored. In this study, we performed a comparative DNA methylome analysis of two selectively bred C. gigas strains with contrasted phenotypes in growth and shell color based on whole-genome bisulfite sequencing (WGBS). Genome-wide profiling of DNA methylation at the single-base resolution revealed that DNA methylations were widely spread across the genome with obvious hotspots, coinciding with the distribution of genes and repetitive elements. Higher methylation levels were observed within genic regions compared with intergenic and promoter regions. Comparative analysis of DNA methylation allowed the identification of 339,604 differentially methylated CpG sites (DMCs) clustering in 27,600 differentially methylated regions (DMRs). Functional annotation analysis identified 11,033 genes from DMRs which were enriched in biological processes including cytoskeleton system, cell cycle, signal transduction, and protein biosynthesis. Integrative analysis of methylome and transcriptome profiles revealed a positive correlation between gene expression and DNA methylation within gene-body regions. Protein-protein interaction (PPI) analysis of differentially expressed and methylated genes allowed for the identification of integrin beta-6 (homolog of human ITGB3) as a hub modulator of the PI3K/Akt signaling pathway that was involved in various growth-related processes. This work provided insights into epigenetic regulation of growth in oysters and will be valuable resources for studying DNA methylation in invertebrates.

Pollution induces epigenetic effects that are stably transmitted across multiple generations

It has been hypothesized that the effects of pollutants on phenotypes can be passed to subsequent generations through epigenetic inheritance, affecting populations long after the removal of a pollutant. But there is still little evidence that pollutants can induce persistent epigenetic effects in animals. Here, we show that low doses of commonly used pollutants induce genome‐wide differences in cytosine methylation in the freshwater crustacean Daphnia pulex. Uniclonal populations were either continually exposed to pollutants or switched to clean water, and methylation was compared to control populations that did not experience pollutant exposure. Although some direct changes to methylation were only present in the continually exposed populations, others were present in both the continually exposed and switched to clean water treatments, suggesting that these modifications had persisted for 7 months (>15 generations). We also identified modifications that were only present in the populations that had switched to clean water, indicating a long‐term legacy of pollutant exposure distinct from the persistent effects. Pollutant‐induced differential methylation tended to occur at sites that were highly methylated in controls. Modifications that were observed in both continually and switched treatments were highly methylated in controls and showed reduced methylation in the treatments. On the other hand, modifications found just in the switched treatment tended to have lower levels of methylation in the controls and showed increase methylation in the switched treatment. In a second experiment, we confirmed that sublethal doses of the same pollutants generate effects on life histories for at least three generations following the removal of the pollutant. Our results demonstrate that even low doses of pollutants can induce transgenerational epigenetic effects that are stably transmitted over many generations. Persistent effects are likely to influence phenotypic development, which could contribute to the rapid adaptation, or extinction, of populations confronted by anthropogenic stressors.

Phenotypic Plasticity: What Has DNA Methylation Got to Do with It?

How does one genome give rise to multiple, often markedly different, phenotypes in response to an environmental cue? This phenomenon, known as phenotypic plasticity, is common amongst plants and animals, but arguably the most striking examples are seen in insects. Well-known insect examples include seasonal morphs of butterfly wing patterns, sexual and asexual reproduction in aphids, and queen and worker castes of eusocial insects. Ultimately, we need to understand how phenotypic plasticity works at a mechanistic level; how do environmental signals alter gene expression, and how are changes in gene expression translated into novel morphology, physiology and behaviour? Understanding how plasticity works is of major interest in evolutionary-developmental biology and may have implications for understanding how insects respond to global change. It has been proposed that epigenetic mechanisms, specifically DNA methylation, are the key link between environmental cues and changes in gene expression. Here, we review the available evidence on the function of DNA methylation of insects, the possible role(s) for DNA methylation in phenotypic plasticity and also highlight key outstanding questions in this field as well as new experimental approaches to address these questions.

OUP accepted manuscript

DNA methylation, an important component of eukaryotic epigenetics, varies in pattern and function across Metazoa. Notably, bilaterian vertebrates and invertebrates differ dramatically in gene body methylation (GbM). Using the frequency of cytosine-phospho-guanines (CpGs), which are lost through mutation when methylated, we report the first broad survey of DNA methylation in Cnidaria, the ancient sister group to Bilateria. We find that: 1) GbM differentially relates to expression categories as it does in most bilaterian invertebrates, but distributions of GbM are less discretely bimodal. 2) Cnidarians generally have lower CpG frequencies on gene bodies than bilaterian invertebrates potentially suggesting a compensatory mechanism to replace CpG lost to mutation in Bilateria that is lacking in Cnidaria. 3) GbM patterns show some consistency within taxonomic groups such as the Scleractinian corals; however, GbM patterns variation across a range of taxonomic ranks in Cnidaria suggests active evolutionary change in GbM within Cnidaria. 4) Some but not all GbM variation is associated with life history change and genome expansion, whereas GbM loss is evident in endoparasitic cnidarians. 5) Cnidarian repetitive elements are less methylated than gene bodies, and methylation of both correlate with genome repeat content. 6) These observations reinforce claims that GbM evolved in stem Metazoa. Thus, this work supports overlap between DNA methylation processes in Cnidaria and Bilateria, provides a framework to compare methylation within and between Cnidaria and Bilateria, and demonstrates the previously unknown rapid evolution of cnidarian methylation.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Tools and applications for integrative analysis of DNA methylation in social insects

DNA methylation is a common epigenetic signalling tool and an important biological process which is widely studied in a large array of species. The presence, level and function of DNA methylation vary greatly across species. In some insects, DNA methylation systems are minimal, and overall methylation rates tend to be low in all studied insect species. Low methylation levels probed by whole-genome bisulphite sequencing require great care with respect to data quality control and interpretation. Here, we introduce BWASP/R, a complete workflow that allows efficient, scalable and entirely reproducible analyses of raw DNA methylation sequencing data. Consistent application of quality control filters and analysis parameters provides fair comparisons among different studies and an integrated view of all experiments on one species. We describe the capabilities of the BWASP/R workflow by re-analysing several publicly available social insect WGBS data sets, comprising 70 samples and cumulatively 147 replicates from four different species. We show that the CpG methylome comprises only about 1.5% of CpG sites in the honeybee genome and that the cumulative data are consistent with genetic signatures of site accessibility and physiological control of methylation levels.

Invertebrate methylomes provide insight into mechanisms of environmental tolerance and reveal methodological biases

There is a growing focus on the role of DNA methylation in the ability of marine invertebrates to rapidly respond to changing environmental factors and anthropogenic impacts. However, genome-wide DNA methylation studies in nonmodel organisms are currently hampered by a limited understanding of methodological biases. Here, we compare three methods for quantifying DNA methylation at single base-pair resolution-whole genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS), and methyl-CpG binding domain bisulfite sequencing (MBDBS)-using multiple individuals from two reef-building coral species with contrasting environmental sensitivity. All methods reveal substantially greater methylation in Montipora capitata (11.4%) than the more sensitive Pocillopora acuta (2.9%). The majority of CpG methylation in both species occurs in gene bodies and flanking regions. In both species, MBDBS has the greatest capacity for detecting CpGs in coding regions at our sequencing depth, but MBDBS may be influenced by intrasample methylation heterogeneity. RRBS yields robust information for specific loci albeit without enrichment of any particular genome feature and with significantly reduced genome coverage. Relative genome size strongly influences the number and location of CpGs detected by each method when sequencing depth is limited, illuminating nuances in cross-species comparisons. As genome-wide methylation differences, supported by data across bisulfite sequencing methods, may contribute to environmental sensitivity phenotypes in critical marine invertebrate taxa, these data provide a genomic resource for investigating the functional role of DNA methylation in environmental tolerance.

Symbiont shuffling induces differential DNA methylation responses to thermal stress in the coral Montastraea cavernosa

Algal symbiont shuffling in favour of more thermotolerant species has been shown to enhance coral resistance to heat-stress. Yet, the mechanistic underpinnings and long-term implications of these changes are poorly understood. This work studied the modifications in coral DNA methylation, an epigenetic mechanism involved in coral acclimatization, in response to symbiont manipulation and subsequent heat stress exposure. Symbiont composition was manipulated in the great star coral Montastraea cavernosa through controlled thermal bleaching and recovery, producing paired ramets of three genets dominated by either their native symbionts (genus Cladocopium) or the thermotolerant species (Durusdinium trenchi). Single-base genome-wide analyses showed significant modifications in DNA methylation concentrated in intergenic regions, introns and transposable elements. Remarkably, DNA methylation changes in response to heat stress were dependent on the dominant symbiont, with twice as many differentially methylated regions found in heat-stressed corals hosting different symbionts (Cladocopium vs. D. trenchii) compared to all other comparisons. Interestingly, while differential gene body methylation was not correlated with gene expression, an enrichment in differentially methylated regions was evident in repetitive genome regions. Overall, these results suggest that changes in algal symbionts favouring heat tolerant associations are accompanied by changes in DNA methylation in the coral host. The implications of these results for coral adaptation, along with future avenues of research based on current knowledge gaps, are discussed in the present work.

Screening bioactive food compounds in honey bees suggests curcumin blocks alcohol-induced damage to longevity and DNA methylation

Various bioactive food compounds may confer health and longevity benefits, possibly through altering or preserving the epigenome. While bioactive food compounds are widely being marketed for human consumption as 'improving health and longevity' by counteracting harmful effects of poor nutrition and lifestyle, claimed effects are often not adequately documented. Using the honey bee (Apis mellifera) as a model species, we here employed a multi-step screening approach to investigate seven compounds for effects on lifespan and DNA methylation using ELISA and whole genome bisulfite sequencing (WGBS). A positive longevity effect was detected for valproic acid, isovaleric acid, and cyanocobalamin. For curcumin, we found that lifespan shortening caused by ethanol intake, was restored when curcumin and ethanol were co-administered. Furthermore, we identified region specific DNA methylation changes as a result of ethanol intake. Ethanol specific changes in DNA methylation were fully or partially blocked in honey bees receiving ethanol and curcumin together. Ethanol-affected and curcumin-blocked differentially methylated regions covered genes involved in fertility, temperature regulation and tubulin transport. Our results demonstrate fundamental negative effects of low dose ethanol consumption on lifespan and associated DNA methylation changes and present a proof-of-principle on how longevity and DNA methylation changes can be negated by the bioactive food component curcumin. Our findings provide a fundament for further studies of curcumin in invertebrates.

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.

Knockdown of DNA methyltransferase 1 reduces DNA methylation and alters expression patterns of cardiac genes in embryonic cardiomyocytes

We previously found that DNA methyltransferase 3a (DNMT3a) plays an important role in regulating embryonic cardiomyocyte gene expression, morphology, and function. In this study, we investigated the role of the most abundant DNMT in mammalian cells, DNMT1, in these processes. It is known that DNMT1 is essential for embryonic development, during which it is involved in regulating cardiomyocyte DNA methylation and gene expression. We used siRNA to knock down DNMT1 expression in primary cultures of mouse embryonic cardiomyocytes. Immunofluorescence staining and multielectrode array were, respectively, utilized to evaluate cardiomyocyte growth and electrophysiology. RNA sequencing (RNA‐Seq) and multiplex bisulfite sequencing were, respectively, performed to examine gene expression and promoter methylation. At 72 h post‐transfection, reduction of DNMT1 expression decreased the number and increased the size of embryonic cardiomyocytes. Beat frequency and the amplitude of field action potentials were decreased by DNMT1 siRNA. RNA‐Seq analysis identified 801 up‐regulated genes and 494 down‐regulated genes in the DNMT1 knockdown cells when compared to controls. Pathway analysis of the differentially expressed genes revealed pathways that were associated with cell death and survival, cell morphology, cardiac function, and cardiac disease. Alternative splicing analysis identified 929 differentially expressed exons, including 583 up‐regulated exons and 308 down‐regulated exons. Moreover, decreased methylation levels were found in the promoters of cardiac genes Myh6, Myh7, Myh7b, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb, mef2c, mef2d, Camta2, Cdkn1A, and Cdkn1C. Of these 13 genes, 6 (Myh6, Tnnc1, Tnni3, Tnnt2, Nppa, Nppb) and 1 (Cdkn1C) had increased or decreased gene expression, respectively. Altogether, these data show that DNMT1 is important in embryonic cardiomyocytes by regulating DNA methylation, gene expression, gene splicing, and cell function.

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Methionine as a methyl donor regulates caste differentiation in the European honey bee (Apis mellifera)

Nutrition contributes to honey bee caste differentiation, but the role of individual nutrients is still unclear. Most essential amino acid contents, except that of methionine (Met), are greater in royal jelly than worker jelly. After ∼3.5 d, the Met content in the latter was slightly greater than in the former. Met is the major raw material used in the synthesis of S-adenosyl-L-methionine, an active methyl donor for DNA methylation, which is an epigenetic driver of caste differentiation. Here, we tested whether Met regulates caste differentiation in honey bees by determining its effects on the caste development of bees receiving four diets: the basic, basic + 0.2% Met, basic + 0.2% Met + 20 mg/kg 5-azacytidine, and basic + 20 mg/kg 5-azacytidine. The presence of Met decreased the adult bee body length and the numbers of ovarioles, indicating that Met may direct the development of female larvae toward worker bees. The upregulated expression of SAMS, Dnmt1, and Dnmt3 caused by Met exposure in 4-d-old larvae indicated that the worker-inductive effects of Met may occur through the promotion of DNA methylation. We investigated the co-effects of Met and glucose on bee development, and found that the effects of an increased glucose level on the number of ovarioles and body length did not strengthen the worker-inductive effects caused by Met. Our results contribute to caste development theory and suggest that Met-as a methyl donor-plays a regulatory, but not decisive, role in caste differentiation.

Insects Provide Unique Systems to Investigate How Early-Life Experience Alters the Brain and Behavior

Early-life experiences have strong and long-lasting consequences for behavior in a surprising diversity of animals. Determining which environmental inputs cause behavioral change, how this information becomes neurobiologically encoded, and the functional consequences of these changes remain fundamental puzzles relevant to diverse fields from evolutionary biology to the health sciences. Here we explore how insects provide unique opportunities for comparative study of developmental behavioral plasticity. Insects have sophisticated behavior and cognitive abilities, and they are frequently studied in their natural environments, which provides an ecological and adaptive perspective that is often more limited in lab-based vertebrate models. A range of cues, from relatively simple cues like temperature to complex social information, influence insect behavior. This variety provides experimentally tractable opportunities to study diverse neural plasticity mechanisms. Insects also have a wide range of neurodevelopmental trajectories while sharing many developmental plasticity mechanisms with vertebrates. In addition, some insects retain only subsets of their juvenile neuronal population in adulthood, narrowing the targets for detailed study of cellular plasticity mechanisms. Insects and vertebrates share many of the same knowledge gaps pertaining to developmental behavioral plasticity. Combined with the extensive study of insect behavior under natural conditions and their experimental tractability, insect systems may be uniquely qualified to address some of the biggest unanswered questions in this field.

High-throughput analyses and Bayesian network modeling highlight novel epigenetic Adverse Outcome Pathway networks of DNA methyltransferase inhibitor mediated transgenerational effects

A number of epigenetic modulating chemicals are known to affect multiple generations of a population from a single ancestral exposure, thus posing transgenerational hazards. The present study aimed to establish a high-throughput (HT) analytical workflow for cost-efficient concentration-response analysis of epigenetic and phenotypic effects, and to support the development of novel Adverse Outcome Pathway (AOP) networks for DNA methyltransferase (DNMT) inhibitor-mediated transgenerational effects on aquatic organisms. The model DNMT inhibitor 5-azacytidine (5AC) and the model freshwater crustacean Daphnia magna were used to generate new experimental data and served as prototypes to construct AOPs for aquatic organisms. Targeted HT bioassays (DNMT ELISA, MS-HRM and qPCR) in combination with multigenerational ecotoxicity tests revealed concentration-dependent transgenerational (F0-F3) effects of 5AC on total DNMT activity, DNA promoter methylation, gene body methylation, gene transcription and reproduction. Top sensitive toxicity pathways related to 5AC exposure, such as apoptosis and DNA damage responses were identified in both F0 and F3 using Gaussian Bayesian network modeling. Two novel epigenetic AOP networks on DNMT inhibitor mediated one-generational and transgenerational effects were developed for aquatic organisms and assessed for the weight of evidence. The new HT analytical workflow and AOPs can facilitate future ecological hazard assessment of epigenetic modulating chemicals.

Insecticide exposure affects intergenerational patterns of DNA methylation in the Colorado potato beetle, Leptinotarsa decemlineata

Insecticide use is pervasive as a selective force in modern agroecosystems. Insect herbivores exposed to these insecticides have been able to rapidly evolve resistance to them, but how they are able to do so is poorly understood. One possible but largely unexplored explanation is that exposure to sublethal doses of insecticides may alter epigenetic patterns that are heritable. For instance, epigenetic mechanisms, such as DNA methylation that modifies gene expression without changing the underlying genetic code, may facilitate the emergence of resistant phenotypes in complex ways. We assessed the effects of sublethal insecticide exposure, with the neonicotinoid imidacloprid, on DNA methylation in the Colorado potato beetle, Leptinotarsa decemlineata, examining both global changes in DNA methylation and specific changes found within genes and transposable elements. We found that exposure to insecticide led to decreases in global DNA methylation for parent and F2 generations and that many of the sites of changes in methylation are found within genes associated with insecticide resistance, such as cytochrome P450s, or within transposable elements. Exposure to sublethal doses of insecticide caused heritable changes in DNA methylation in an agricultural insect herbivore. Therefore, epigenetics may play a role in insecticide resistance, highlighting a fundamental mechanism of evolution while informing how we might better coexist with insect species in agroecosystems.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

DNA methylation changes in response to ocean acidification at the time of larval metamorphosis in the edible oyster, Crassostrea hongkongensis

Unprecedented rate of increased CO₂ level in the ocean and the subsequent changes in carbonate system including decreased pH, known as ocean acidification (OA), is predicted to disrupt not only the calcification process but also several other physiological and developmental processes in a variety of marine organisms, including edible oysters. Nonetheless, not all species are vulnerable to those OA threats, e.g. some species may be able to cope with OA stress using environmentally induced modifications on gene and protein expressions. For example, external environmental stressors including OA can influence the addition and removal of methyl groups through epigenetic modification (e.g. DNA methylation) process to turn gene expression "on or off" as part of a rapid adaptive mechanism to cope with OA. In this study, we tested the above hypothesis through testing the effect of OA, using decreased pH 7.4 as proxy, on DNA methylation pattern of an endemic and a commercially important estuary oyster species, Crassostrea hongkongensis at the time of larval habitat selection and metamorphosis. Larval growth rate did not differ between control pH 8.1 and treatment pH 7.4. The metamorphosis rate of the pediveliger larvae was higher at pH 7.4 than those in control pH 8.1, however over one-third of the larvae raised at pH 7.4 failed to attach on optimal substrate as defined by biofilm presence. During larval development, a total of 130 genes were differentially methylated across the two treatments. The differential methylation in the larval genes may have partially accounted for the higher metamorphosis success rate under decreased pH 7.4 but with poor substratum selection ability. Differentially methylated loci were concentrated in the exon regions and appear to be associated with cytoskeletal and signal transduction, oxidative stress, metabolic processes, and larval metamorphosis, which implies the high potential of C. hongkongensis larvae to acclimate and adapt through non-genetic ways to OA threats within a single generation.

Genome-Wide Differential DNA Methylation in Reproductive, Morphological, and Visual System Differences Between Queen Bee and Worker Bee (Apis mellifera)

There are many differences in external morphology and internal physiology between the Apis mellifera queen bee and worker bee, some of which are relevant to beekeeping production. These include reproductive traits, body size, royal jelly secreting properties, and visual system development, among others. The identification of candidate genes that control the differentiation of these traits is critical for selective honeybee breeding programs. In this study, we compared the genomic methylation of queen bee and worker bee larvae at 3, 4, and 5 days of age by whole-genome bisulfite sequencing, and found that the basic characteristics of genomic methylation in queen and worker larvae were the same. There were approximately 49 million cytosines in the Apis larvae genome, of which about 90,000 were methylated. Methylated CpG sites accounted for 99% of the methylated cytosines, and methylation mainly occurred in exons. However, methylation levels of queen and worker larvae showed different trends with age: the methylation level of queen larvae varied with age in an inverted parabola, while the corresponding trend for worker larvae with resembled an exponential curve with a platform. The methylation level of queen larvae was higher than that of worker larvae at 3 days of age, lower than that of worker larvae at 4 days of age, and similar to that of worker larvae at 5 days old. The top 10 differentially methylated genes (DMGs) and 13 caste-specific methylated genes were listed, and correlations with caste determination were speculated. We additionally screened 38 DMGs between queen larvae and worker larvae involved in specific organ differentiation as well as reproduction, morphology, and vision differentiation during caste determination. These genes are potential molecular markers for selective breeding of A. mellifera to improve fecundity, royal jelly production, body size, and foraging, and represent candidate genes for investigating specialized functional segregation during the process of caste differentiation.

The prognostic value of a Methylome-based Malignancy Density Scoring System to predict recurrence risk in early-stage Lung Adenocarcinoma

Current NCCN guidelines do not recommend the use of adjuvant chemotherapy for stage IA lung adenocarcinoma patients with R0 surgery. However, 25% to 40% of patients with stage IA disease experience recurrence. Stratifying patients according to the recurrence risk may tailor adjuvant therapy and surveillance imaging for those with a higher risk. However, prognostic markers are often identified by comparing high-risk and low-risk cases which might introduce bias due to the widespread interpatient heterogeneity. Here, we developed a scoring system quantifying the degree of field cancerization in adjacent normal tissues and revealed its association with disease-free survival (DFS).

Methods: We recruited a cohort of 44 patients with resected stage IA lung adenocarcinoma who did not receive adjuvant therapy. Both tumor and adjacent normal tissues were obtained from each patient and subjected to capture-based targeted genomic and epigenomic profiling. A novel methylome-based scoring system namely malignancy density ratio (MD ratio) was developed based on 39 patients by comparing tumor and corresponding adjacent normal tissues of each patient. A MD score was then obtained by Wald statistics. The correlations of MD ratio, MD score, and genomic features with clinical outcome were investigated.

Results: Patients with a high-risk MD ratio showed a significantly shorter postsurgical DFS compared with those with a low-risk MD ratio (HR=4.47, P=0.01). The MD ratio was not associated with T stage (P=1), tumor cell fraction (P=0.748) nor inflammatory status (p=0.548). Patients with a high-risk MD score also demonstrated an inferior DFS (HR=4.69, P=0.039). In addition, multivariate analysis revealed EGFR 19 del (HR=5.39, P=0.012) and MD score (HR= 7.90, P=0.01) were independent prognostic markers.

Conclusion: The novel methylome-based scoring system, developed by comparing the signatures between tumor and corresponding adjacent normal tissues of individual patients, largely minimizes the bias of interpatient heterogeneity and reveals a robust prognostic value in patients with resected lung adenocarcinoma.

Characterization of Global DNA Methylation in Different Gene Regions Reveals Candidate Biomarkers in Pigs with High and Low Levels of Boar Taint

DNA methylation of different gene components, including different exons and introns, or different lengths of exons and introns is associated with differences in gene expression. To investigate the methylation of porcine gene components associated with the boar taint (BT) trait, this study used reduced representation bisulfite sequencing (RRBS) data from nine porcine testis samples in three BT groups (low, medium and high BT). The results showed that the methylation levels of the first exons and first introns were lower than those of the other exons and introns. The first exons/introns of CpG island regions had even lower levels of methylation. A total of 123 differentially methylated promoters (DMPs), 194 differentially methylated exons (DMEs) and 402 differentially methylated introns (DMIs) were identified, of which 80 DMPs (DMP-CpGis), 112 DMEs (DME-CpGis) and 166 DMIs (DMI-CpGis) were discovered in CpG islands. Importantly, GPX1 contained one each of DMP, DME, DMI, DMP-CpGi, DME-CpGi and DMI-CpGi. Gene-GO term relationships and pathways analysis showed DMP-CpGi-related genes are mainly involved in methylation-related biological functions. In addition, gene-gene interaction networks consisted of nodes that were hypo-methylated GPX1, hypo-methylated APP, hypo-methylated ATOX1, hyper-methylated ADRB2, hyper-methylated RPS6KA1 and hyper-methylated PNMT. They could be used as candidate biomarkers for reducing boar taint in pigs, after further validation in large cohorts.

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.

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.

Myocardial tissue-specific Dnmt1 knockout in rats protects against pathological injury induced by Adriamycin

Novel molecular mechanisms of the pathophysiology of heart failure (HF) are continuously being discovered, including epigenetic regulation. Among epigenetic marks, the role of DNA hypomethylation in shaping heart morphology and function in vivo and the pathogenesis of cardiomyopathy and/or HF, especially in adults, has not been clearly established. Here we show that the strong expression of DNA methyltransferase 1 (Dnmt1) is obviously downregulated in the WT adult rat heart with age. By contrast, the expression of Dnmt1 is upregulated suddenly in heart tissues from pressure overload-induced HF mice and adriamycin-induced cardiac injury and HF mice, consistent with the increased expression of Dnmt1 observed in familial hypertrophic cardiomyopathy (FHCM) patients. To further assess the role of Dnmt1, we generated myocardium-specific Dnmt1 knockout (Dnmt1 KO) rats using CRISPR-Cas9 technology. Echocardiographic and histopathological examinations demonstrated that Dnmt1 deficiency is associated with resistance to cardiac pathological changes and protection at the global and organization levels in response to pathological stress. Furthermore, Dnmt1 deficiency in the myocardium restricts the expressional reprogramming of genes and activates pathways involved in myocardial protection and anti-apoptosis in response to pathological stress. Transcriptome and genome-wide DNA methylation analyses revealed that these changes in regulation are linked to alterations in the methylation status of genes due to Dnmt1 knockout. The present study is the first to investigate in vivo the impact of genome-wide cardiac DNA methyltransferase deficiency on physiological development and the pathological processes of heart tissues in response to stress. The exploration of the role of epigenetics in the development, modification, and prevention of cardiomyopathy and HF is in a very preliminary stage but has an infinite future.

Current evidence for a role of epigenetic mechanisms in response to ionizing radiation in an ecotoxicological context

The issue of potential long-term or hereditary effects for both humans and wildlife exposed to low doses (or dose rates) of ionising radiation is a major concern. Chronic exposure to ionising radiation, defined as an exposure over a large fraction of the organism's lifespan or even over several generations, can possibly have consequences in the progeny. Recent work has begun to show that epigenetics plays an important role in adaptation of organisms challenged to environmental stimulae. Changes to so-called epigenetic marks such as histone modifications, DNA methylation and non-coding RNAs result in altered transcriptomes and proteomes, without directly changing the DNA sequence. Moreover, some of these environmentally-induced epigenetic changes tend to persist over generations, and thus, epigenetic modifications are regarded as the conduits for environmental influence on the genome. Here, we review the current knowledge of possible involvement of epigenetics in the cascade of responses resulting from environmental exposure to ionising radiation. In addition, from a comparison of lab and field obtained data, we investigate evidence on radiation-induced changes in the epigenome and in particular the total or locus specific levels of DNA methylation. The challenges for future research and possible use of changes as an early warning (biomarker) of radiosensitivity and individual exposure is discussed. Such a biomarker could be used to detect and better understand the mechanisms of toxic action and inter/intra-species susceptibility to radiation within an environmental risk assessment and management context.

Programmable targeted epigenetic editing using CRISPR system in Bombyx mori

DNA methylation has been proven to play roles in regulating gene expression, cell fate, disease determination, and chromatin architecture organization in mammals and plants, and is a significant component of epigenetic modification. Compared to mammals or plants, the status and function of DNA methylation are poorly understood in insects, which is partially due to the lack of efficient manipulation tools. In this study, we show that fusion protein of catalytically inactive Cas9 (dCas9) with TET1 can efficiently demethylate genomic DNA of silkworm Bombyx mori, in a programmable target region specific manner. We first developed an all-in-one vector to maximize the targeting efficiency of dCas9-TET1. Then we selected 3 endogenous genes that were previously found to harbor methylated DNA, and designed gRNAs within the methylated region. Co-transfection of dCas9-TET1 and gRNA successfully erased methylation marks near the targeting region, with efficiencies from about 17.50% to 40.00%. Furthermore, targeted demethylation on gene body resulted in increased mRNA transcription level. Unlike the previously widely used decitabine, a methylation inhibitor, dCas9-TET1 is more effective and specific, and has no unwanted impact on whole-genome methylation. DCas9-TET1 provides a powerful tool for investigating the functional significance of DNA methylation in a locus-specific manner, and for exploring the unknown links between methylation and development in insects.

Epigenetic, transcriptional and phenotypic responses in two generations of Daphnia magna exposed to the DNA methylation inhibitor 5-azacytidine

The water flea Daphnia magna is a keystone species in freshwater ecosystems and has been widely used as a model organism in environmental ecotoxicology. This aquatic crustacean is sensitive to environmental stressors and displays considerable plasticity in adapting to changing environmental conditions. Part of this plasticity may be due to epigenetic regulation of gene expression, including changes to DNA methylation and histone modifications. Because of the generally hypomethylated genome of this species, we hypothesized that the histone code may have an essential role in the epigenetic control and that histone modifications might be an early marker for stress. This study aims to characterize the epigenetic, transcriptional and phenotypic responses and their causal linkages in directly exposed adult (F0) Daphnia and peritoneal exposed neonates (F1) after a chronic (7-day) exposure to a sublethal concentration (10 mg/l) of 5-azacytidine, a well-studied vertebrate DNA methylation inhibitor. Exposure of the F0 generation significantly reduced the cumulative fecundity, accompanied with differential expression of genes in the one-carbon-cycle metabolic pathway. In the epigenome of the F0 generation, a decrease in global DNA methylation, but no significant changes on H3K4me3 or H3K27me3, were observed. In the F1 offspring generation, changes in gene expression, a significant reduction in global DNA methylation and changes in histone modifications were identified. The results indicate that exposure during adulthood may result in more pronounced effects on early development in the offspring generation, though interpretation of the data should be carefully done since both the exposure regime and developmental period is different in the two generations examined. The obtained results improve our understanding of crustacean epigenetics and the tools developed may promote use of epigenetic markers in hazard assessment of environmental stressors.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii)

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

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.

Genome-wide analysis reveals that exon methylation facilitates its selective usage in the human transcriptome

DNA methylation, especially in promoter regions, is a well-characterized epigenetic marker related to gene expression regulation in eukaryotes. However, the role of intragenic DNA methylation in the usage of corresponding exons still remains elusive. In this study, we described the DNA methylome across 10 human tissues. The human genome showed both conserved and varied methylation levels among these tissues. We found that the methylation densities in promoters and first exons were negatively correlated with the corresponding gene expression level. Nevertheless, the methylation densities within introns, internal exons and down 1 kb regions showed weak correlation with gene expression levels. Importantly, we observed a remarkably positive relationship between methylation density and exon expression level of intragenic exons. Notably, skip-in exons were much more methylated than skip-out exons. We also identified 260 exons that showed both differential methylation levels and differential expression levels in lung cancer. Genes harboring these differentially regulated exons were significantly enriched in the cancer hallmark-related biological process. Moreover, a 10-exon signature was identified as a promising prognostic predictor for lung cancer. Our study illuminates the DNA methylome, describes its relationship with gene expression across human tissues and provides new insights into intragenic DNA methylation and exon usage during the transcriptional/alternative splicing process and in cancer.