Comparative epigenomics: a powerful tool to understand the evolution of DNA methylation

Epigenetics Research in Evolutionary Biology: Perspectives on Timescales and Mechanisms

Epigenetics research in evolutionary biology encompasses a variety of research areas, from regulation of gene expression to inheritance of environmentally mediated phenotypes. Such divergent research foci can occasionally render the umbrella term "epigenetics" ambiguous. Here I discuss several areas of contemporary epigenetics research in the context of evolutionary biology, aiming to provide balanced views across timescales and molecular mechanisms. The importance of epigenetics in development is now being assessed in many nonmodel species. These studies not only confirm the importance of epigenetic marks in developmental processes, but also highlight the significant diversity in epigenetic regulatory mechanisms across taxa. Further, these comparative epigenomic studies have begun to show promise toward enhancing our understanding of how regulatory programs evolve. A key property of epigenetic marks is that they can be inherited along mitotic cell lineages, and epigenetic differences that occur during early development can have lasting consequences on the organismal phenotypes. Thus, epigenetic marks may play roles in short-term (within an organism's lifetime or to the next generation) adaptation and phenotypic plasticity. However, the extent to which observed epigenetic variation occurs independently of genetic influences remains uncertain, due to the widespread impact of genetics on epigenetic variation and the limited availability of comprehensive (epi)genomic resources from most species. While epigenetic marks can be inherited independently of genetic sequences in some species, there is little evidence that such "transgenerational inheritance" is a general phenomenon. Rather, molecular mechanisms of epigenetic inheritance are highly variable between species.

Novel function of the PsDMAP1 protein in regulating the growth and pathogenicity of Phytophthora sojae

The DNA methyltransferase 1-associated protein (DMAP1) was initially identified as an activator of DNA methyltransferase 1 (DNMT1), a conserved eukaryotic enzyme involved in diverse molecular processes, including histone acetylation and chromatin remodeling. However, the roles and regulatory mechanisms of DMAP1 in filamentous pathogens are still largely unknown. Here, employing bioinformatic analysis, we identified PsDMAP1 in P. sojae, which features a canonical histone tail-binding domain, as the ortholog of the human DMAP1. A phylogenetic analysis of DMAP1 protein sequences across diverse eukaryotic organisms revealed the remarkable conservation and distinctiveness of oomycete DMAP1 orthologs. Homozygous knockout of PsDMAP1 resulted in the mortality of P. sojae. Furthermore, silencing of PsDMAP1 caused a pronounced reduction in mycelial growth, production of sporangia and zoospore, cystospore germination, and virulence. PsDMAP1 also played a crucial role in the response of P. sojae to reactive oxygen species (ROS) and osmotic stresses. Moreover, PsDMAP1 interacted with DNA N6-methyladenine (6 mA) methyltransferase PsDAMT1, thereby enhancing its catalytic activity and effectively regulating 6 mA abundance in P. sojae. Our findings reveal the functional importance of PsDAMP1 in the development and infection of P. sojae, and this marks the initial exploration of the novel 6 mA regulator PsDMAP1 in plant pathogens.

Comparative epigenome analysis using Infinium DNA methylation BeadChips

The arrival of the Infinium DNA methylation BeadChips for mice and other nonhuman mammalian species has outpaced the development of the informatics that supports their use for epigenetics study in model organisms. Here, we present informatics infrastructure and methods to allow easy DNA methylation analysis on multiple species, including domesticated animals and inbred laboratory mice (in SeSAMe version 1.16.0+). First, we developed a data-driven analysis pipeline covering species inference, genome-specific data preprocessing and regression modeling. We targeted genomes of 310 species and 37 inbred mouse strains and showed that genome-specific preprocessing prevents artifacts and yields more accurate measurements than generic pipelines. Second, we uncovered the dynamics of the epigenome evolution in different genomic territories and tissue types through comparative analysis. We identified a catalog of inbred mouse strain-specific methylation differences, some of which are linked to the strains' immune, metabolic and neurological phenotypes. By streamlining DNA methylation array analysis for undesigned genomes, our methods extend epigenome research to broad species contexts.

DNA methylation-free Arabidopsis reveals crucial roles of DNA methylation in regulating gene expression and development

A contribution of DNA methylation to defense against invading nucleic acids and maintenance of genome integrity is uncontested; however, our understanding of the extent of involvement of this epigenetic mark in genome-wide gene regulation and plant developmental control is incomplete. Here, we knock out all five known DNA methyltransferases in Arabidopsis, generating DNA methylation-free plants. This quintuple mutant exhibits a suite of developmental defects, unequivocally demonstrating that DNA methylation is essential for multiple aspects of plant development. We show that CG methylation and non-CG methylation are required for a plethora of biological processes, including pavement cell shape, endoreduplication, cell death, flowering, trichome morphology, vasculature and meristem development, and root cell fate determination. Moreover, we find that DNA methylation has a strong dose-dependent effect on gene expression and repression of transposable elements. Taken together, our results demonstrate that DNA methylation is dispensable for Arabidopsis survival but essential for the proper regulation of multiple biological processes.

Our understanding of the extent of involvement of DNA methylation in genome-wide gene regulation and plant developmental control is incomplete. Here, the authors knock out all five known DNA methyltransferases and show the developmental and gene expression changes in the DNA methylation-free Arabidopsis plants.

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.

Transcriptional memory and response to adverse temperatures in plants

Temperature is one of the major environmental signals controlling plant development, geographical distribution, and seasonal behavior. Plants perceive adverse temperatures, such as high, low, and freezing temperatures, as stressful signals that can cause physiological defects and even death. As sessile organisms, plants have evolved sophisticated mechanisms to adapt to recurring stressful environments through changing gene expression or transcriptional reprogramming. Transcriptional memory refers to the ability of primed plants to remember previously experienced stress and acquire enhanced tolerance to similar or different stresses. Epigenetic modifications mediate transcriptional memory and play a key role in adapting to adverse temperatures. Understanding the mechanisms of the formation, maintenance, and resetting of stress-induced transcriptional memory will not only enable us to understand why there is a trade-off between plant defense and growth, but also provide a theoretical basis for generating stress-tolerant crops optimized for future climate change. In this review, we summarize recent advances in dissecting the mechanisms of plant transcriptional memory in response to adverse temperatures, based mainly on studies of the model plant Arabidopsis thaliana. We also discuss remaining questions that are important for further understanding the mechanisms of transcriptional memory during the adverse temperature response.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Epigenetic Regulation and Environmental Sex Determination in Cichlid Fishes

Studying environmental sex determination (ESD) in cichlids provides a phylogenetic and comparative approach to understand the evolution of the underlying mechanisms, their impact on the evolution of the overlying systems, and the neuroethology of life history strategies. Natural selection normally favors parents who invest equally in the development of male and female offspring, but evolution may favor deviations from this 50:50 ratio when environmental conditions produce an advantage for doing so. Many species of cichlids demonstrate ESD in response to water chemistry (temperature, pH, and oxygen concentration). The relative strengths of and the exact interactions between these factors vary between congeners, demonstrating genetic variation in sensitivity. The presence of sizable proportions of the less common sex towards the environmental extremes in most species strongly suggests the presence of some genetic sex-determining loci acting in parallel with the ESD factors. Sex determination and differentiation in these species does not seem to result in the organization of a final and irreversible sexual fate, so much as a life-long ongoing battle between competing male- and female-determining genetic and hormonal networks governed by epigenetic factors. We discuss what is and is not known about the epigenetic mechanism behind the differentiation of both gonads and sex differences in the brain. Beyond the well-studied tilapia species, the 2 best-studied dwarf cichlid systems showing ESD are the South American genus Apistogramma and the West African genus Pelvicachromis. Both species demonstrate male morphs with alternative reproductive tactics. We discuss the further neuroethology opportunities such systems provide to the study of epigenetics of alternative life history strategies and other behavioral variation.

Comparative Epigenomics Reveals Host Diversity of the Trichinella Epigenomes and Their Effects on Differential Parasitism

Comparative epigenomics provides new insights on evolutionary biology in relation with complex interactions between species and their environments. In the present study, we focus on deciphering the conservation and divergence of DNA methylomes during Trichinella evolution. Whole-genome bisulfite sequencing and RNA-seq were performed on the two clades of Trichinella species, in addition to whole-genome sequencing. We demonstrate that methylation patterns of sing-copy orthologous genes (SCOs) of the 12 Trichinella species are host-related and can mirror known phylogenetic relationships. Among these SCOs, we identify a panel of genes exhibiting hyper-/hypo-methylated features in gene-bodies or respective promoters that play pivotal roles in transcriptome regulation. These hyper-/hypo-methylated SCOs are also of functional significance across developmental stages, as they are highly enriched species-specific and stage-specific expressed genes both in Ad and ML stages. We further identify a set of parasitism-related functional genes that exhibit host-related differential methylation and expression among those SCOs, including p53-like transcription factor and Cdc37 that are of functional significance for elucidating differential parasitology between the two clades of Trichinella. This comparative epigenome study can help to decipher the environmental effects on differential adaptation and parasitism of the genus Trichinella.

The epigenome and beyond: How does non-genetic inheritance change our view of evolution?

Evidence from across the tree of life suggests that epigenetic inheritance is more common than previously thought. If epigenetic inheritance is indeed as common as the data suggest, this finding has potentially important implications for evolutionary theory and our understanding of how evolution and adaptation progress. However, we currently lack an understanding of how common various epigenetic inheritance types are, and how they impact phenotypes. In this perspective, we review the open questions that need to be addressed to fully integrate epigenetic inheritance into evolutionary theory and to develop reliable predictive models for phenotypic evolution. We posit that addressing these challenges will require the collaboration of biologists from different disciplines and a focus on the exploration of data and phenomena without preconceived limits on potential mechanisms or outcomes.

F-box protein CFK1 interacts with and degrades de novo DNA methyltransferase in Arabidopsis

DNA methylation plays crucial roles in cellular development and stress responses through gene regulation and genome stability control. Precise regulation of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), the de novo Arabidopsis DNA methyltransferase, is crucial to maintain DNA methylation homeostasis to ensure genome integrity. Compared with the extensive studies on DRM2 targeting mechanisms, little information is known regarding the quality control of DRM2 itself. Here, we conducted yeast two-hybrid screen assay and identified an E3 ligase, COP9 INTERACTING F-BOX KELCH 1 (CFK1), as a novel DRM2-interacting partner and targets DRM2 for degradation via the ubiquitin-26S proteasome pathway in Arabidopsis thaliana. We also performed whole genome bisulfite sequencing (BS-seq) to determine the biological significance of CFK1-mediated DRM2 degradation. Loss-of-function CFK1 leads to increased DRM2 protein abundance and overexpression of CFK1 showed reduced DRM2 protein levels. Consistently, CFK1 overexpression induces genome-wide CHH hypomethylation and transcriptional de-repression at specific DRM2 target loci. This study uncovered a distinct mechanism regulating de novo DNA methyltransferase by CFK1 to control DNA methylation level.

A potential screening method for epigenetic drugs: uncovering stress-induced gene silencing in Chlamydomonas

Histone modifications and DNA methylation together govern promoter availability, thereby influencing gene expression. This study queries the unicellular chlorophyte, Chlamydomonas reinhardtii using a three step "epigenetic assay" design to phenotypically track the variegation of a randomly integrated Paromomycin resistance transgene(s) (Pm^R). Based on its position of integration, the Pm^R gene expression hinged on two epigenetic hallmarks: the spreading of heterochromatin, and the transmissible memory of epigenetic states across generations. Using a spot-dilution analysis, the loss of antibiotic resistance phenotype was scored from 0 to 4, four being maximally silenced. Appropriate construct designs were used to demonstrate that the cis-spread of heterochromatin could be interfered with a stronger euchromatic barrier (TUB2 promoter). When assayed for metal ion stress, a combination of Mn deficiency with excess Cu or Zn stress was shown to induce gene silencing in Chlamydomonas. Cu stress resulted in the accumulation of intracellular ROS, while Zn stress elevated the sensitivity to ROS. As proof of functional conservation, mammalian epigenetic drugs demonstrably interfered with stress-induced gene silencing. Finally, a selected group of transgenic clones responsive to HDACi sodium butyrate, when tested in a gradient plate format showed similarity in phenotype to the plant-derived compound cinnamic acid. This indicated a possible commonality in their mode of action, unlike curcumin which might have a different mechanism. Thus, using binned libraries, based on a common set of responses to known drugs, a cost-effective high-throughput screening strategy for epigenetically active compounds from plants or other sources is described.

Evidence for divergence of DNA methylation maintenance and a conserved inhibitory mechanism from DNA demethylation in chickens and mammals

BACKGROUND

DNA methylation is a significant epigenetic modification that is evolutionarily conserved in various species and often serves as a repressive mark for transcription. DNA methylation levels and patterns are regulated by a balance of opposing enzyme functions, DNA methyltransferases, DNMT1/3A/3B and methylcytosine dioxygenases, TET1/2/3. In mice, the TET enzyme converts DNA cytosine methylation (5mC) to 5-hydroxymethylcytosine (5hmC) at the beginning of fertilisation and gastrulation and initiates a global loss of 5mC, while the 5mC level is increased on the onset of cell differentiation during early embryonic development.

OBJECTIVE

Global loss and gain of DNA methylation may be differently regulated in diverged species.

METHODS

Chicken B-cell lymphoma DT40 cells were used as an avian model to compare differences in the overall regulation of DNA modification with mammals.

RESULTS

We found that DNA methylation is maintained at high levels in DT40 cells through compact chromatin formation, which inhibits TET-mediated demethylation. Human and mouse chromosomes introduced into DT40 cells by cell fusion lost the majority of 5mC, except for human subtelomeric repeats.

CONCLUSION

Our attempt to elucidate the differences in the epigenetic regulatory mechanisms between birds and mammals explored the evidence that they share a common chromatin-based regulation of TET-DNA access, while chicken DNMT1 is involved in different target sequence recognition systems, suggesting that factors inducing DNMT-DNA association have already diverged.

Taking the Wheel - de novo DNA Methylation as a Driving Force of Plant Embryonic Development

During plant embryogenesis, regardless of whether it begins with a fertilized egg cell (zygotic embryogenesis) or an induced somatic cell (somatic embryogenesis), significant epigenetic reprogramming occurs with the purpose of parental or vegetative transcript silencing and establishment of a next-generation epigenetic patterning. To ensure genome stability of a developing embryo, large-scale transposon silencing occurs by an RNA-directed DNA methylation (RdDM) pathway, which introduces methylation patterns de novo and as such potentially serves as a global mechanism of transcription control during developmental transitions. RdDM is controlled by a two-armed mechanism based around the activity of two RNA polymerases. While PolIV produces siRNAs accompanied by protein complexes comprising the methylation machinery, PolV produces lncRNA which guides the methylation machinery toward specific genomic locations. Recently, RdDM has been proposed as a dominant methylation mechanism during gamete formation and early embryo development in Arabidopsis thaliana, overshadowing all other methylation mechanisms. Here, we bring an overview of current knowledge about different roles of DNA methylation with emphasis on RdDM during plant zygotic and somatic embryogenesis. Based on published chromatin immunoprecipitation data on PolV binding sites within the A. thaliana genome, we uncover groups of auxin metabolism, reproductive development and embryogenesis-related genes, and discuss possible roles of RdDM at the onset of early embryonic development via targeted methylation at sites involved in different embryogenesis-related developmental mechanisms.

Application of the MSAP Technique to Evaluate Epigenetic Changes in Plant Conservation

Epigenetic variation, and particularly DNA methylation, is involved in plasticity and responses to changes in the environment. Conservation biology studies have focused on the measurement of this variation to establish demographic parameters, diversity levels and population structure to design the appropriate conservation strategies. However, in ex situ conservation approaches, the main objective is to guarantee the characteristics of the conserved material (phenotype and epi-genetic). We review the use of the Methylation Sensitive Amplified Polymorphism (MSAP) technique to detect changes in the DNA methylation patterns of plant material conserved by the main ex situ plant conservation methods: seed banks, in vitro slow growth and cryopreservation. Comparison of DNA methylation patterns before and after conservation is a useful tool to check the fidelity of the regenerated plants, and, at the same time, may be related with other genetic variations that might appear during the conservation process (i.e., somaclonal variation). Analyses of MSAP profiles can be useful in the management of ex situ plant conservation but differs in the approach used in the in situ conservation. Likewise, an easy-to-use methodology is necessary for a rapid interpretation of data, in order to be readily implemented by conservation managers.

Comparative epigenetics in animal physiology: An emerging frontier

The unprecedented access to annotated genomes now facilitates the investigation of the molecular basis of epigenetic phenomena in phenotypically diverse animals. In this critical review, we describe the roles of molecular epigenetic mechanisms in regulating mitotically and meiotically stable spatiotemporal gene expression, phenomena that provide the molecular foundation for the intra-, inter-, and trans-generational emergence of physiological phenotypes. By focusing principally on emerging comparative epigenetic roles of DNA-level and transcriptome-level epigenetic mark dynamics in the emergence of phenotypes, we highlight the relationship between evolutionary conservation and innovation of specific epigenetic pathways, and their interplay as a priority for future study. This comparative approach is expected to significantly advance our understanding of epigenetic phenomena, as animals show a diverse array of strategies to epigenetically modify physiological responses. Additionally, we review recent technological advances in the field of molecular epigenetics (single-cell epigenomics and transcriptomics and editing of epigenetic marks) in order to (1) investigate environmental and endogenous factor dependent epigenetic mark dynamics in an integrative manner; (2) functionally test the contribution of specific epigenetic marks for animal phenotypes via genome and transcript-editing tools. Finally, we describe advantages and limitations of emerging animal models, which under the Krogh principle, may be particularly useful in the advancement of comparative epigenomics and its potential translational applications in animal science, ecotoxicology, ecophysiology, climate change science and wild-life conservation, as well as organismal health.

Targeted bisulfite sequencing for biomarker discovery

Cytosine methylation is one of the best studied epigenetic modifications. In mammals, DNA methylation patterns vary among cells and is mainly found in the CpG context. DNA methylation is involved in important processes during development and differentiation and its dysregulation can lead to or is associated with diseases, such as cancer, loss-of-imprinting syndromes and neurological disorders. It has been also shown that DNA methylation at the cellular, tissue and organism level varies with age. To overcome the costs of Whole-Genome Bisulfite Sequencing, the gold standard method to detect 5-methylcytosines at a single base resolution, DNA methylation arrays have been developed and extensively used. This method allows one to assess the status of a fraction of the CpG sites present in the genome of an organism. In order to combine the relatively low cost of Methylation Arrays and digital signals of bisulfite sequencing, we developed a Targeted Bisulfite Sequencing method that can be applied to biomarker discovery for virtually any phenotype. Here we describe a comprehensive step-by-step protocol to build a DNA methylation-based epigenetic clock.

Remodeling the epigenome and (epi)cytoskeleton: a new paradigm for co-regulation by methylation

The epigenome determines heritable patterns of gene expression in the absence of changes in DNA sequence. The result is programming of different cellular-, tissue- and organ-specific phenotypes from a single organismic genome. Epigenetic marks that comprise the epigenome (e.g. methylation) are placed upon or removed from chromatin (histones and DNA) to direct the activity of effectors that regulate gene expression and chromatin structure. Recently, the cytoskeleton has been identified as a second target for the cell's epigenetic machinery. Several epigenetic 'readers, writers and erasers' that remodel chromatin have been discovered to also remodel the cytoskeleton, regulating structure and function of microtubules and actin filaments. This points to an emerging paradigm for dual-function remodelers with 'chromatocytoskeletal' activity that can integrate cytoplasmic and nuclear functions. For example, the SET domain-containing 2 methyltransferase (SETD2) has chromatocytoskeletal activity, methylating both histones and microtubules. The SETD2 methyl mark on chromatin is required for efficient DNA repair, and its microtubule methyl mark is required for proper chromosome segregation during mitosis. This unexpected convergence of SETD2 activity on histones and microtubules to maintain genomic stability suggests the intriguing possibility of an expanded role in the cell for chromatocytoskeletal proteins that read, write and erase methyl marks on the cytoskeleton as well as chromatin. Coordinated use of methyl marks to remodel both the epigenome and the (epi)cytoskeleton opens the possibility for integrated regulation (which we refer to as 'epiregulation') of other higher-level functions, such as muscle contraction or learning and memory, and could even have evolutionary implications.

Epigenomic Regulatory Mechanism in Vegetative Phase Transition of Malus hupehensis

In woody plants, phase transitions substantially affect growth and development. Although there has been considerable interest in the regulatory mechanisms underlying phase changes, the associated epigenetic modifications remain relatively uncharacterized. We examined the DNA methylation changes and the transcriptional responses in adult and juvenile Malus hupehensis leaves. The DNA methylations were 66.61% and 68.3% in the CG context, 49.12% and 52.44% in the CHG context, and 7.02% and 8.22% in the CHH context for the adult and juvenile leaves, respectively. The number of differentially methylated regions in all contexts distributed in the genic regions varied. Additionally, inhibited DNA methylation in adult leaves activated the transcription of indole-3-acetic acid related genes in the signaling, response, and transport pathways. Moreover, the opposite methylation and expression patterns were observed for the SPL and AP2 family genes between the adult and juvenile leaves. Both gene families contribute to the M. hupehensis vegetative phase transition. Furthermore, the hyper-/hypomethylation of the gene body or promoter of transcription factor genes may lead to up-/downregulated gene expression. The methylation levels of the WRKY (22), NAC (21), ERF (8), WOX (2), KNAT (6), EIN3 (2), SCL (7), ZAT (7), and HSF (4) genes were higher in the adult leaves than in the juvenile leaves, whereas the opposite pattern was observed for the TCP (2), MADS-box (11), and DOF (3) genes. An analysis of the correlation between methylation and transcription indicated the methylation of the gene body in all contexts and the methylation of the promoter in the CG and CHG contexts are negatively correlated with gene expression. However, the methylation of the promoter in the CHH context is positively correlated with gene expression. These findings reflect the diversity in the epigenetic regulation of gene expression and may be useful for elucidating the epigenetic regulatory mechanism underlying the M. hupehensis vegetative phase transition.

Prospects for incorporation of epigenetic biomarkers in human health and environmental risk assessment of chemicals

Epigenetic mechanisms have gained relevance in human health and environmental studies, due to their pivotal role in disease, gene × environment interactions and adaptation to environmental change and/or contamination. Epigenetic mechanisms are highly responsive to external stimuli and a wide range of chemicals has been shown to determine specific epigenetic patterns in several organisms. Furthermore, the mitotic/meiotic inheritance of such epigenetic marks as well as the resulting changes in gene expression and cell/organismal phenotypes has now been demonstrated. Therefore, epigenetic signatures are interesting candidates for linking environmental exposures to disease as well as informing on past exposures to stressors. Accordingly, epigenetic biomarkers could be useful tools in both prospective and retrospective risk assessment but epigenetic endpoints are currently not yet incorporated into risk assessments. Achieving a better understanding on this apparent impasse, as well as identifying routes to promote the application of epigenetic biomarkers within environmental risk assessment frameworks are the objectives of this review. We first compile evidence from human health studies supporting the use of epigenetic exposure-associated changes as reliable biomarkers of exposure. Then, specifically focusing on environmental science, we examine the potential and challenges of developing epigenetic biomarkers for environmental fields, and discuss useful organisms and appropriate sequencing techniques to foster their development in this context. Finally, we discuss the practical incorporation of epigenetic biomarkers in the environmental risk assessment of chemicals, highlighting critical data gaps and making key recommendations for future research within a regulatory context.

Functions and mechanisms of plant histone deacetylases

Lysine acetylation, one of the major types of post-translational modifications, plays critical roles in regulating gene expression and protein function. Histone deacetylases (HDACs) are responsible for removing acetyl groups from lysines of both histone and non-histone proteins. While tremendous progress has been made in understanding the function and mechanism of HDACs in animals in the past two decades, nearly half of the HDAC studies in plants were reported within the past five years. In this review, we summarize the major findings on plant HDACs, with a focus on the model plant Arabidopsis thaliana, and highlight the components, regulatory mechanisms, and biological functions of HDAC complexes.

Establishment, maintenance, and biological roles of non-CG methylation in plants

Cytosine DNA methylation is prevalent throughout eukaryotes and prokaryotes. While most commonly thought of as being localized to dinucleotide CpG sites, non-CG sites can also be modified. Such non-CG methylation is widespread in plants, occurring at trinucleotide CHG and CHH (H = A, T, or C) sequence contexts. The prevalence of non-CG methylation in plants is due to the plant-specific CHROMOMETHYLASE (CMT) and RNA-directed DNA Methylation (RdDM) pathways. These pathways have evolved through multiple rounds of gene duplication and gene loss, generating epigenomic variation both within and between species. They regulate both transposable elements and genes, ensure genome integrity, and ultimately influence development and environmental responses. In these capacities, non-CG methylation influence and shape plant genomes.

DNA methylation-regulated QPCT promotes sunitinib resistance by increasing HRAS stability in renal cell carcinoma

Rationale: Although sunitinib has been shown to improve the survival rate of advanced renal cell carcinoma (RCC) patients, poor drug response is a major challenge that reduces patient benefit. It is important to elucidate the underlying mechanism so that the therapeutic response to sunitinib can be restored. Methods: We used an Illumina HumanMethylation 850K microarray to find methylation-differentiated CpG sites between sunitinib-nonresponsive and -responsive RCC tissues and Sequenom MassARRAY methylation analysis to verify the methylation chip results. We verified glutaminyl peptide cyclotransferase (QPCT) expression in sunitinib-nonresponsive and -responsive RCC tissues via qRT-PCR, western blot and immunohistochemical assays. Then, cell counting kit 8 (CCK-8), plate colony formation and flow cytometric assays were used to verify the function of QPCT in RCC sunitinib resistance after QPCT intervention or overexpression. Chromatin immunoprecipitation (ChIP) was performed to clarify the upstream regulatory mechanism of QPCT. A human proteome microarray assay was used to identify downstream proteins that interact with QPCT, and co-immunoprecipitation (co-IP) and confocal laser microscopy were used to verify the protein chip results. Results: We found that the degree of methylation in the QPCT promoter region was significantly different between sunitinib-nonresponsive and -responsive RCC tissues. In the sunitinib-nonresponsive tissues, the degree of methylation in the QPCT promoter region was significantly reduced, and the expression of QPCT was upregulated, which correlated with a clinically poor response to sunitinib. A knockdown of QPCT conferred sunitinib sensitivity traits to RCC cells, whereas an overexpression of QPCT restored sunitinib resistance in RCC cells. Mechanistically, reducing the methylation degree of the QPCT promoter region by 5-aza-2'-deoxycytidine (decitabine) in RCC cells could increase the expression of QPCT and NF-κB (p65) bound to the QPCT promoter region, positively regulating its expression, while the hypermethylation in the QPCT promoter region could inhibit the binding of NF-κB (p65). QPCT could bind to HRAS and attenuate the ubiquitination of HRAS, thus increasing its stability and leading to the activation of the ERK pathway in RCC cells. Conclusion: QPCT may be a novel predictor of the response to sunitinib therapy in RCC patients and a potential therapeutic target.

Comparative RNA-Sequencing and DNA Methylation Analyses of Apple (Malus domestica Borkh.) Buds with Diverse Flowering Capabilities Reveal Novel Insights into the Regulatory Mechanisms of Flower Bud Formation

In plants, DNA methylation (i.e. chromatin modification) is important for various biological processes, including growth, development and flowering. Because 'Fuji' apple trees are alternate bearing and have a long ripening period and poor-quality flower buds, we used bud types with diverse flowering capabilities to investigate the epigenetic regulatory mechanisms influencing flower bud formation. We examined the DNA methylation changes and the transcriptional responses in the selected apple bud types. We observed that in the apple genome, approximately 79.5%, 67.4% and 23.7% of the CG, CHG and CHH sequences are methylated, respectively. For each sequence context, differentially methylated regions exhibited distinct methylation patterns among the analyzed apple bud types. Global methylation and transcriptional analyses revealed that nonexpressed genes or genes expressed at low levels were highly methylated in the gene-body regions, suggesting that gene-body methylation is negatively correlated with gene expression. Moreover, genes with methylated promoters were more highly expressed than genes with unmethylated promoters, implying promoter methylation and gene expression are positively correlated. Additionally, flowering-related genes (e.g. SOC1, AP1 and SPLs) and some transcription factor genes (e.g. GATA, bHLH, bZIP and WOX) were highly expressed in spur buds (highest flowering rate), but were associated with low methylation levels in the gene-body regions. Our findings indicate a potential correlation between DNA methylation and gene expression in apple buds with diverse flowering capabilities, suggesting an epigenetic regulatory mechanism influences apple flower bud formation.

Tracing the origin and evolution history of methylation-related genes in plants

BACKGROUND

DNA methylation is a crucial epigenetic modification, which is involved in many biological processes, including gene expression regulation, embryonic development, cell differentiation and genomic imprinting etc. And it also involves many key regulatory genes in eukaryotes. By tracing the evolutionary history of methylation-related genes, we can understand the origin and expansion time of these genes, which helps to understand the evolutionary history of plants, and we can also understand the changes of DNA methylation patterns in different species. However, most studies on the evolution of methylation-related genes failed to be carried out for the whole DNA methylation pathway.

RESULTS

In this study, we conducted a comprehensive identification of 33 methylation-related genes in 77 species, and investigated gene origin and evolution throughout the plant kingdom. We found that the origin of genes responsible for methylation maintenance and demethylation evolved early, while most de novo methylation-related genes appeared late. The methylation-related genes were expanded by whole genome duplication and tandem replication, but were also accompanied by a large number of gene absence events in different species. The gene length and intron length varied a lot in different species, but exon structure and functional domains were relatively conserved. The phylogenetic relationships of methylation-related genes were traced to reveal the evolution history of DNA methylation in different species. The expression patterns of methylation-related genes have changed during the evolution of species, and the expression patterns of these genes in different species can be clustered into four categories.

CONCLUSIONS

The study describes a global characterization of DNA methylation-related genes in the plant kingdom. The similarities and differences in origin time, gene structure and phylogenetic relationship of these genes lead us to understand the evolutionary conservation and dynamics of DNA methylation in plants.

Retroelement-Linked Transcription Factor Binding Patterns Point to Quickly Developing Molecular Pathways in Human Evolution

Background: Retroelements (REs) are transposable elements occupying ~40% of the human genome that can regulate genes by providing transcription factor binding sites (TFBS). RE-linked TFBS profile can serve as a marker of gene transcriptional regulation evolution. This approach allows for interrogating the regulatory evolution of organisms with RE-rich genomes. We aimed to characterize the evolution of transcriptional regulation for human genes and molecular pathways using RE-linked TFBS accumulation as a metric. Methods: We characterized human genes and molecular pathways either enriched or deficient in RE-linked TFBS regulation. We used ENCODE database with mapped TFBS for 563 transcription factors in 13 human cell lines. For 24,389 genes and 3124 molecular pathways, we calculated the score of RE-linked TFBS regulation reflecting the regulatory evolution rate at the level of individual genes and molecular pathways. Results: The major groups enriched by RE regulation deal with gene regulation by microRNAs, olfaction, color vision, fertilization, cellular immune response, and amino acids and fatty acids metabolism and detoxication. The deficient groups were involved in translation, RNA transcription and processing, chromatin organization, and molecular signaling. Conclusion: We identified genes and molecular processes that have characteristics of especially high or low evolutionary rates at the level of RE-linked TFBS regulation in human lineage.

Phytophthora methylomes are modulated by 6mA methyltransferases and associated with adaptive genome regions

BACKGROUND

Filamentous plant pathogen genomes often display a bipartite architecture with gene-sparse, repeat-rich compartments serving as a cradle for adaptive evolution. The extent to which this two-speed genome architecture is associated with genome-wide DNA modifications is unknown.

RESULTS

We show that the oomycetes Phytophthora infestans and Phytophthora sojae possess functional adenine N6-methylation (6mA) methyltransferases that modulate patterns of 6mA marks across the genome. In contrast, 5-methylcytosine could not be detected in these species. Methylated DNA IP sequencing (MeDIP-seq) of each species reveals 6mA is depleted around the transcription start sites (TSSs) and is associated with lowly expressed genes, particularly transposable elements. Genes occupying the gene-sparse regions have higher levels of 6mA in both genomes, possibly implicating the methylome in adaptive evolution. All six putative adenine methyltransferases from P. infestans and P. sojae, except PsDAMT2, display robust enzymatic activities. Surprisingly, single knockouts in P. sojae significantly reduce in vivo 6mA levels, indicating that the three enzymes are not fully redundant. MeDIP-seq of the psdamt3 mutant reveals uneven 6mA methylation reduction across genes, suggesting that PsDAMT3 may have a preference for gene body methylation after the TSS. Furthermore, transposable elements such as DNA elements are more active in the psdamt3 mutant. A large number of genes, particularly those from the adaptive genomic compartment, are differentially expressed.

CONCLUSIONS

Our findings provide evidence that 6mA modification is potentially an epigenetic mark in Phytophthora genomes, and complex patterns of 6mA methylation may be associated with adaptive evolution in these important plant pathogens.

Transcriptomic analysis of the salivary gland of medicinal leech Hirudo nipponia

Hirudo nipponia (known as Shui Zhi in Chinese) is a well-known Chinese medicine with numerous active ingredients in its body, especially in its saliva. This native Chinese blood-sucking leech has been used for therapeutic purposes since before 100 AD. Modern Chinese physicians use it for a wide range of diseases. Genomic data and molecular information about the pharmacologically active substances produced by this medicinal leech are presently unavailable despite this organism’s medicinal importance. In this study, we performed transcriptome profiling of the salivary glands of medicinal leech H. nipponia using the Illumina platform. In total, 84,657,362 clean reads were assembled into 50,535 unigenes. The obtained unigenes were compared to public databases. Furthermore, a unigene sequence similarity search and comparisons with the whole transcriptome of medical leech were performed to identify potential proteins. Finally, more than 21 genes were predicted to be involved in anticoagulatory, antithrombotic, antibacterial, anti-inflammatory and antitumor processes, which might play important roles in the treatment of various diseases. This study is the first analysis of a sialotranscriptome in H. nipponia. The transcriptome profile will shed light on its genetic background and provide a useful tool to deepen our understanding of the medical value of H. nipponia.

5-Azacytidine: A Promoter of Epigenetic Changes in the Quest to Improve Plant Somatic Embryogenesis

Somatic embryogenesis (SE) is a widely studied process due to its biotechnological potential to generate large quantities of plants in short time frames and from different sources of explants. The success of SE depends on many factors, such as the nature of the explant, the microenvironment generated by in vitro culture conditions, and the regulation of gene expression, among others. Epigenetics has recently been identified as an important factor influencing SE outcome. DNA methylation is one of the most studied epigenetic mechanisms due to its essential role in gene expression, and its participation in SE is crucial. DNA methylation levels can be modified through the use of drugs such as 5-Azacytidine (5-AzaC), an inhibitor of DNA methylation, which has been used during SE protocols. The balance between hypomethylation and hypermethylation seems to be the key to SE success. Here, we discuss the most prominent recent research on the role of 5-AzaC in the regulation of DNA methylation, highlighting its importance during the SE process. Also, the molecular implications that this inhibitor might have for the increase or decrease in the embryogenic potential of various explants are reviewed.

Statistical modelling of CG interdistance across multiple organisms

Background

Statistical approaches to genetic sequences have revealed helpful to gain deeper insight into biological and structural functionalities, using ideas coming from information theory and stochastic modelling of symbolic sequences. In particular, previous analyses on CG dinucleotide position along the genome allowed to highlight its epigenetic role in DNA methylation, showing a different distribution tail as compared to other dinucleotides.

In this paper we extend the analysis to the whole CG distance distribution over a selected set of higher-order organisms. Then we apply the best fitting probability density function to a large range of organisms (>4400) of different complexity (from bacteria to mammals) and we characterize some emerging global features.

Results

We find that the Gamma distribution is optimal for the selected subset as compared to a group of several distributions, chosen for their physical meaning or because recently used in literature for similar studies. The parameters of this distribution, when applied to our larger set of organisms, allows to highlight some biologically relavant features for the considered organism classes, that can be useful also for classification purposes.

Conclusions

The quantification of statistical properties of CG dinucleotide positioning along the genome is confirmed as a useful tool to characterize broad classes of organisms, spanning the whole range of biological complexity.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2303-2) contains supplementary material, which is available to authorized users.

Pattern of DNA Methylation in Daphnia: Evolutionary Perspective

DNA methylation is an evolutionary ancient epigenetic modification that is phylogenetically widespread. Comparative studies of the methylome across a diverse range of non-conventional and conventional model organisms is expected to help reveal how the landscape of DNA methylation and its functions have evolved. Here, we explore the DNA methylation profile of two species of the crustacean Daphnia using whole genome bisulfite sequencing. We then compare our data with the methylomes of two insects and two mammals to achieve a better understanding of the function of DNA methylation in Daphnia. Using RNA-sequencing data for all six species, we investigate the correlation between DNA methylation and gene expression. DNA methylation in Daphnia is mainly enriched within the coding regions of genes, with the highest methylation levels observed at exons 2–4. In contrast, vertebrate genomes are globally methylated, and increase towards the highest methylation levels observed at exon 2, and maintained across the rest of the gene body. Although DNA methylation patterns differ among all species, their methylation profiles share a bimodal distribution across the genomes. Genes with low levels of CpG methylation and gene expression are mainly enriched for species specific genes. In contrast, genes associated with high methylated CpG sites are highly transcribed and evolutionary conserved across all species. Finally, the positive correlation between internal exons and gene expression potentially points to an evolutionary conserved mechanism, whereas the negative regulation of gene expression via methylation of promoters and exon 1 is potentially a secondary mechanism that has been evolved in vertebrates.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

Recent Insights on Biological and Ecological Aspects of Ectomycorrhizal Fungi and Their Interactions

The roots of most terrestrial plants are colonized by mycorrhizal fungi. They play a key role in terrestrial environments influencing soil structure and ecosystem functionality. Around them a peculiar region, the mycorrhizosphere, develops. This is a very dynamic environment where plants, soil and microorganisms interact. Interest in this fascinating environment has increased over the years. For a long period the knowledge of the microbial populations in the rhizosphere has been limited, because they have always been studied by traditional culture-based techniques. These methods, which only allow the study of cultured microorganisms, do not allow the characterization of most organisms existing in nature. The introduction in the last few years of methodologies that are independent of culture techniques has bypassed this limitation. This together with the development of high-throughput molecular tools has given new insights into the biology, evolution, and biodiversity of mycorrhizal associations, as well as, the molecular dialog between plants and fungi. The genomes of many mycorrhizal fungal species have been sequenced so far allowing to better understanding the lifestyle of these fungi, their sexual reproduction modalities and metabolic functions. The possibility to detect the mycelium and the mycorrhizae of heterothallic fungi has also allowed to follow the spatial and temporal distributional patterns of strains of different mating types. On the other hand, the availability of the genome sequencing from several mycorrhizal fungi with a different lifestyle, or belonging to different groups, allowed to verify the common feature of the mycorrhizal symbiosis as well as the differences on how different mycorrhizal species interact and dialog with the plant. Here, we will consider the aspects described before, mainly focusing on ectomycorrhizal fungi and their interactions with plants and other soil microorganisms.

Single-base methylome analysis reveals dynamic epigenomic differences associated with water deficit in apple

Cytosine methylation is an essential feature of epigenetic regulation and is involved in various biological processes. Although cytosine methylation has been analysed at the genomic scale for several plant species, there is a general lack of understanding of the dynamics of global and genic DNA methylation in plants growing in environments challenged with biotic and abiotic stresses. In this study, we mapped cytosine methylation at single-base resolution in the genome of commercial apple (Malus x domestica), and analysed changes in methylation patterns associated with water deficit in representative drought-sensitive and drought-tolerant cultivars. We found that the apple genome exhibits ~54%, ~38% and ~8.5% methylation at CG, CHG and CHH sequence contexts, respectively. We additionally documented changes in gene expression associated with water deficit in an attempt to link methylation and gene expression changes. Global methylation and transcription analysis revealed that promoter-unmethylated genes showed higher expression levels than promoter-methylated genes. Gene body methylation appears to be positively correlated with gene expression. Water deficit stress was associated with changes in methylation at a multitude of genes, including those encoding transcription factors (TFs) and transposable elements (TEs). These results present a methylome map of the apple genome and reveal widespread DNA methylation alterations in response to water deficit stress. These data will be helpful for understanding potential linkages between DNA methylation and gene expression in plants growing in natural environments and challenged with abiotic and biotic stresses.

Epigenetics and DOHaD: from basics to birth and beyond

Developmental origins of health and disease (DOHaD) is the study of how the early life environment can impact the risk of chronic diseases from childhood to adulthood and the mechanisms involved. Epigenetic modifications such as DNA methylation, histone modifications and non-coding RNAs are involved in mediating how early life environment impacts later health. This review is a summary of the Epigenetics and DOHaD workshop held at the 2016 DOHaD Society of Australia and New Zealand Conference. Our extensive knowledge of how the early life environment impacts later risk for chronic disease would not have been possible without animal models. In this review we highlight some animal model examples that demonstrate how an adverse early life exposure results in epigenetic and gene expression changes that may contribute to increased risk of chronic disease later in life. Type 2 diabetes and cardiovascular disease are chronic diseases with an increasing incidence due to the increased number of children and adults that are obese. Epigenetic changes such as DNA methylation have been shown to be associated with metabolic health measures and potentially predict future metabolic health status. Although more difficult to elucidate in humans, recent studies suggest that DNA methylation may be one of the epigenetic mechanisms that mediates the effects of early life exposures on later life risk of obesity and obesity related diseases. Finally, we discuss the role of the microbiome and how it is a new player in developmental programming and mediating early life exposures on later risk of chronic disease.

De novo transcriptome assembly, functional annotation and differential gene expression analysis of juvenile and adult E. fetida, a model oligochaete used in ecotoxicological studies

Background

Earthworms are sensitive to toxic chemicals present in the soil and so are useful indicator organisms for soil health. Eisenia fetida are commonly used in ecotoxicological studies; therefore the assembly of a baseline transcriptome is important for subsequent analyses exploring the impact of toxin exposure on genome wide gene expression.

Results

This paper reports on the de novo transcriptome assembly of E. fetida using Trinity, a freely available software tool. Trinotate was used to carry out functional annotation of the Trinity generated transcriptome file and the transdecoder generated peptide sequence file along with BLASTX, BLASTP and HMMER searches and were loaded into a Sqlite3 database. To identify differentially expressed transcripts; each of the original sequence files were aligned to the de novo assembled transcriptome using Bowtie and then RSEM was used to estimate expression values based on the alignment. EdgeR was used to calculate differential expression between the two conditions, with an FDR corrected P value cut off of 0.001, this returned six significantly differentially expressed genes. Initial BLASTX hits of these putative genes included hits with annelid ferritin and lysozyme proteins, as well as fungal NADH cytochrome b5 reductase and senescence associated proteins. At a cut off of P = 0.01 there were a further 26 differentially expressed genes.

Conclusion

These data have been made publicly available, and to our knowledge represent the most comprehensive available transcriptome for E. fetida assembled from RNA sequencing data. This provides important groundwork for subsequent ecotoxicogenomic studies exploring the impact of the environment on global gene expression in E. fetida and other earthworm species.

Electronic supplementary material

The online version of this article (doi:10.1186/s40659-017-0114-y) contains supplementary material, which is available to authorized users.

Putting DNA methylation in context: from genomes to gene expression in plants

Plant DNA methylation is its own language, interpreted by the cell to maintain silencing of transposons, facilitate chromatin structure, and to ensure proper expression of some genes. Just as in any language, context is important. Rather than being a simple "on-off switch", DNA methylation has a range of "meanings" dependent upon the underlying sequence and its location in the genome. Differences in the sequence context of individual sites are established, maintained, and interpreted by differing molecular pathways. Varying patterns of methylation within genes and surrounding sequences are associated with a continuous range of expression differences, from silencing to constitutive expression. These often-subtle differences have been pieced together from years of effort, but have taken off with the advent of methods for assessing methylation across entire genomes. Recognizing these patterns and identifying underlying causes is essential for understanding the function of DNA methylation and its systems-wide contribution to a range of processes in plant genomes. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Epigenetic Inheritance and Its Role in Evolutionary Biology: Re-Evaluation and New Perspectives

Epigenetics increasingly occupies a pivotal position in our understanding of inheritance, natural selection and, perhaps, even evolution. A survey of the PubMed database, however, reveals that the great majority (>93%) of epigenetic papers have an intra-, rather than an inter-generational focus, primarily on mechanisms and disease. Approximately ~1% of epigenetic papers even mention the nexus of epigenetics, natural selection and evolution. Yet, when environments are dynamic (e.g., climate change effects), there may be an “epigenetic advantage” to phenotypic switching by epigenetic inheritance, rather than by gene mutation. An epigenetically-inherited trait can arise simultaneously in many individuals, as opposed to a single individual with a gene mutation. Moreover, a transient epigenetically-modified phenotype can be quickly “sunsetted”, with individuals reverting to the original phenotype. Thus, epigenetic phenotype switching is dynamic and temporary and can help bridge periods of environmental stress. Epigenetic inheritance likely contributes to evolution both directly and indirectly. While there is as yet incomplete evidence of direct permanent incorporation of a complex epigenetic phenotype into the genome, doubtlessly, the presence of epigenetic markers and the phenotypes they create (which may sort quite separately from the genotype within a population) will influence natural selection and, so, drive the collective genotype of a population.