DNA Methylation Variation Research Articles

Variation in DNA methylation and response to short-term herbivory in Thlaspi arvense

Plant metabolic pathways and gene networks involved in the response to herbivory are well-established, but the impact of epigenetic factors as modulators of those responses is less understood. Here, we use the demethylating agent 5-azacytidine to uncover the role of DNA cytosine methylation on phenotypic responses after short-term herbivory in Thlaspi arvense plants that came from two European populations with contrasting flowering phenotypes expected to differ in the response to experimental demethylation. The experimental design followed a 2 × 3 factorial design, that was replicated for each flowering-type. First, half the seeds were submerged in a water solution of 5-azacytidine and the other half only in water, as controls. Then, we assigned control and demethylated plants to three herbivory categories (i) insect herbivory, (ii) artificial herbivory, and (iii) undamaged plants. The effects of the demethylation and herbivory treatments were assessed by quantifying genome-wide global DNA cytosine methylation, concentration of leaf glucosinolates, final stem biomass, fruit and seed production, and seed size. For most of the plant traits analysed, individuals from the two flowering-types responded differently. In late-flowering plants, global DNA methylation did not differ between control and demethylated plants but it was significantly reduced by herbivory. Conversely, in early-flowering plants, demethylation at seed stage was still evident in leaf DNA of reproductive individuals whereas herbivory did not affect their global DNA methylation. In late-flowering plants, artificial herbivory imposed a stronger reduction than insect herbivory in global DNA methylation and final stem biomass, and induced higher concentration of aliphatic glucosinolates. In early-flowering plants, the effects of herbivory were non-significant for the same traits. Finally, the effect of herbivory on reproductive parameters varied with the level of demethylation and the plant flowering-type. Although further investigations with more populations and families are required to confirm our results, they suggest that the genetic background of experimental plants and timing of damage can affect the response to herbivory and point towards multifaceted genetic-epigenetic interactions in determining herbivory-induced phenotypic plasticity.

Association of childhood BMI trajectory with post-adolescent and adult lung function is mediated by pre-adolescent DNA methylation

BackgroundBody mass index (BMI) has been shown to be associated with lung function. Recent findings showed that DNA methylation (DNAm) variation is likely to be a consequence of changes in BMI. However, whether DNAm mediates the association of BMI with lung function is unknown. We examined the mediating role of DNAm on the association of pre-adolescent BMI trajectories with post-adolescent and adulthood lung function (forced expiratory volume (FEV1), forced vital capacity (FVC), and FEV1/FVC).MethodsAnalyses were undertaken in the Isle of Wight birth cohort (IOWBC). Group-based trajectory modelling was applied to infer latent BMI trajectories from age 1 to 10 years. An R package, ttscreening, was applied to identify CpGs at 10 years potentially associated with BMI trajectories for each sex. Linear regressions were implemented to further screen CpGs for their association with lung function at 18 years. Path analysis, stratified by sex, was applied to each screened CpG to assess its role of mediation. Internal validation was applied to further examine the mediation consistency of the detected CpGs based on lung function at 26 years. Mendelian randomization (MR-base) was used to test possible causal effects of the identified CpGs.ResultsTwo BMI trajectories (high vs. low) were identified. Of the 442,475 CpG sites, 18 CpGs in males and 33 in females passed screening. Eight CpGs in males and 16 CpGs in females (none overlapping) were identified as mediators. For subjects with high BMI trajectory, high DNAm at all CpGs in males were associated with decreased lung function, while 8 CpGs in females were associated with increased lung function at 18 years. At 26 years, 6 CpGs in males and 14 CpGs in females showed the same direction of indirect effects as those at 18 years. DNAm at CpGs cg19088553 (GRIK2) and cg00612625 (HPSE2) showed a potential causal effect on FEV1.ConclusionsThe effects of BMI trajectory in early childhood on post-adolescence lung function were likely to be mediated by pre-adolescence DNAm in both males and females, but such mediation effects were likely to diminish over time.

Temporal variation in DNA methylation during gonadal development in a reptile with temperature-dependent sex determination†.

DNA methylation plays a significant role in transducing external environmental signals to a cellular response in reptiles; however, whether the methylation patterns are conserved across species remains unclear. Here, we examined the genome-wide DNA methylation differentiation between male and female hatchling gonads of the temperature-dependent sex determination (TSD) Mauremys mutica (M. mutica) using methylation-dependent restriction-site associated DNA sequencing (MethylRAD-seq) to test differentially methylated genes underlying sexual development. Several categories, including heat-shock genes (HSP90A, HSP30C), histone- (KDM8) and ubiquitin-related genes (TRIM39), kinases (WNK3), and gonad differentiation or gonadal-development-related genes (HSD17B8, HSD17B12), were identified as candidates for future study. Additionally, we identified several regulatory pathways potentially mediating TSD thermosensitivity such as the GnRH signaling pathway and calcium signaling pathway. These findings provide evidence that sexually dimorphic DNA methylation may be associated with sex determination or sex differentiation in TSD M. mutica.

Evaluating the Microsatellite Instability of Colorectal Cancer Based on Multimodal Deep Learning Integrating Histopathological and Molecular Data.

Microsatellite instability (MSI), an important biomarker for immunotherapy and the diagnosis of Lynch syndrome, refers to the change of microsatellite (MS) sequence length caused by insertion or deletion during DNA replication. However, traditional wet-lab experiment-based MSI detection is time-consuming and relies on experimental conditions. In addition, a comprehensive study on the associations between MSI status and various molecules like mRNA and miRNA has not been performed. In this study, we first studied the association between MSI status and several molecules including mRNA, miRNA, lncRNA, DNA methylation, and copy number variation (CNV) using colorectal cancer data from The Cancer Genome Atlas (TCGA). Then, we developed a novel deep learning framework to predict MSI status based solely on hematoxylin and eosin (H&E) staining images, and combined the H&E image with the above-mentioned molecules by multimodal compact bilinear pooling. Our results showed that there were significant differences in mRNA, miRNA, and lncRNA between the high microsatellite instability (MSI-H) patient group and the low microsatellite instability or microsatellite stability (MSI-L/MSS) patient group. By using the H&E image alone, one can predict MSI status with an acceptable prediction area under the curve (AUC) of 0.809 in 5-fold cross-validation. The fusion models integrating H&E image with a single type of molecule have higher prediction accuracies than that using H&E image alone, with the highest AUC of 0.952 achieved when combining H&E image with DNA methylation data. However, prediction accuracy will decrease when combining H&E image with all types of molecular data. In conclusion, combining H&E image with deep learning can predict the MSI status of colorectal cancer, the accuracy of which can further be improved by integrating appropriate molecular data. This study may have clinical significance in practice.

DNA Methylome and LncRNAome Analysis Provide Insights Into Mechanisms of Genome-Dosage Effects in Autotetraploid Cassava.

Whole genome duplication (WGD) increases the dosage of all coding and non-coding genes, yet the molecular implications of genome-dosage effects remain elusive. In this study, we generated integrated maps of the methylomes and lncRNAomes for diploid and artificially generated autotetraploid cassava (Manihot esculenta Crantz). We found that transposable elements (TEs) suppressed adjacent protein coding gene (PCG)-expression levels, while TEs activated the expression of nearby long non-coding RNAs (lncRNAs) in the cassava genome. The hypermethylation of DNA transposons in mCG and mCHH sites may be an effective way to suppress the expression of nearby PCGs in autotetraploid cassava, resulting in similar expression levels for most of PCGs between autotetraploid and diploid cassava. In the autotetraploid, decreased methylation levels of retrotransposons at mCHG and mCHH sites contributed to reduced methylation of Gypsy-neighboring long intergenic non-coding RNAs, potentially preserving diploid-like expression patterns in the major of lncRNAs. Collectively, our study highlighted that WGD-induced DNA methylation variation in DNA transposons and retrotransposons may be as direct adaptive responses to dosage of all coding-genes and lncRNAs, respectively.

Examining the epigenetic mechanisms of childhood adversity and sensitive periods: A gene set-based approach

BackgroundSensitive periods are developmental stages of heightened plasticity when life experiences, including exposure to childhood adversity, have the potential to exert more lasting impacts. Epigenetic mechanisms, including DNA methylation (DNAm), may provide a pathway through which adversity induces long-term biological changes. DNAm shifts may be more likely to occur during sensitive periods, especially within genes that regulate the timing of sensitive periods. Here, we investigated the possibility that childhood adversity during specific life stages is associated with DNAm changes in genes known to regulate the timing and duration of sensitive periods. MethodsGenome-wide DNAm profiles came from the Avon Longitudinal Study of Parents and Children (n = 785). We first used principal component analysis (PCA) to summarize DNAm variation across 530 CpG sites mapped to the promoters of 58 genes previously-identified as regulating sensitive periods. Gene-level DNAm summaries were calculated for genes regulating sensitive period opening (ngenes = 15), closing (ngenes = 36), and expression (ngenes = 8). We then performed linear discriminant analysis (LDA) to test associations between seven types of parent-reported, time-varying measures of exposure to childhood adversity and DNAm principal components. To our knowledge, this is the first time LDA has been applied to analyze functionally grouped DNAm data to characterize associations between an environmental exposure and epigenetic differences. ResultsSuggestive evidence emerged for associations between sexual or physical abuse as well as financial hardship during middle childhood, and DNAm of genetic pathways regulating sensitive period opening and expression. However, no statistically significant associations were identified after multiple testing correction. ConclusionsOur gene set-based method combining PCA and LDA complements epigenome-wide approaches. Although our results were largely null, these findings provide a proof-of-concept for studying time-varying exposures and gene- or pathway-level epigenetic modifications.

Are Alterations in DNA Methylation Related to CKD Development?

The modifications in genomic DNA methylation are involved in the regulation of normal and pathological cellular processes. The epigenetic regulation stimulates biological plasticity as an adaptive response to variations in environmental factors. The role of epigenetic changes is vital for the development of some diseases, including atherogenesis, cancers, and chronic kidney disease (CKD). The results of studies presented in this review have suggested that altered DNA methylation can modulate the expression of pro-inflammatory and pro-fibrotic genes, as well those essential for kidney development and function, thus stimulating renal disease progression. Abnormally increased homocysteine, hypoxia, and inflammation have been suggested to alter epigenetic regulation of gene expression in CKD. Studies of renal samples have demonstrated the relationship between variations in DNA methylation and fibrosis and variations in estimated glomerular filtration rate (eGFR) in human CKD. The unravelling of the genetic–epigenetic profile would enhance our understanding of processes underlying the development of CKD. The understanding of multifaceted relationship between DNA methylation, genes expression, and disease development and progression could improve the ability to identify individuals at risk of CKD and enable the choice of appropriate disease management.

DNA Methylation and RNA-Sequencing Analysis to Identify Genes Related to Spontaneous Leaf Spots in a Wheat Variety ‘Zhongkenuomai No.1’

Greenish leaf variation has been reported widely as a trait of great interest in wheat for improving photosynthesis. Zhongkenuomai No.1 (ZKNM1), a mutant with recoverable leaf spots, was regarded as a suitable material for studying chlorophyll synthesis-related mechanisms. In this study, transcriptome and DNA methylation analyses were conducted in ZKNM1 leaves to determine the transcriptional regulatory mechanism of leaf spot development. Ultimately, 890 differentially expressed genes (DEGs) were discovered, with chlorophyll biosynthesis pathway genes downregulated and chlorophyll degradation pathway genes upregulated, possibly acting as a double block to chlorophyll accumulation. Among them, HEMA1s (Glutamyl-tRNA reductase family proteins) and PORAs (protochlorophyllide oxidoreductase A) were the most important controlled genes. Furthermore, a genome-wide methylation analysis indicates that a hypermethylated region is present 1690 bp upstream of the transcriptional start sites in spot tissues (SPs), and 131 DNA methylation-mediated DEGs were identified, one of which encoded a putative resistance gene (TraesCS1A02G009500) and was a hub gene in interaction network modules. In the sample groups with leaf spots (SPs), this gene may be involved in the photosynthetic processes. The findings indicated that dynamic variations in DNA methylation play key roles in gene regulation to govern leaf spot development.

Association between DNA methylation variability and self-reported exposure to heavy metals

Individuals encounter varying environmental exposures throughout their lifetimes. Some exposures such as smoking are readily observed and have high personal recall; others are more indirect or sporadic and might only be inferred from long occupational histories or lifestyles. We evaluated the utility of using lifetime-long self-reported exposures for identifying differential methylation in an amyotrophic lateral sclerosis cases-control cohort of 855 individuals. Individuals submitted paper-based surveys on exposure and occupational histories as well as whole blood samples. Genome-wide DNA methylation levels were quantified using the Illumina Infinium Human Methylation450 array. We analyzed 15 environmental exposures using the OSCA software linear and MOA models, where we regressed exposures individually by methylation adjusted for batch effects and disease status as well as predicted scores for age, sex, cell count, and smoking status. We also regressed on the first principal components on clustered environmental exposures to detect DNA methylation changes associated with a more generalised definition of environmental exposure. Five DNA methylation probes across three environmental exposures (cadmium, mercury and metalwork) were significantly associated using the MOA models and seven through the linear models, with one additionally across a principal component representing chemical exposures. Methylome-wide significance for four of these markers was driven by extreme hyper/hypo-methylation in small numbers of individuals. The results indicate the potential for using self-reported exposure histories in detecting DNA methylation changes in response to the environment, but also highlight the confounded nature of environmental exposure in cohort studies.

Identification of COVID-19-Associated DNA Methylation Variations by Integrating Methylation Array and scRNA-Seq Data at Cell-Type Resolution.

Single-cell transcriptome studies have revealed immune dysfunction in COVID-19 patients, including lymphopenia, T cell exhaustion, and increased levels of pro-inflammatory cytokines, while DNA methylation plays an important role in the regulation of immune response and inflammatory response. The specific cell types of immune responses regulated by DNA methylation in COVID-19 patients will be better understood by exploring the COVID-19 DNA methylation variation at the cell-type level. Here, we developed an analytical pipeline to explore single-cell DNA methylation variations in COVID-19 patients by transferring bulk-tissue-level knowledge to the single-cell level. We discovered that the methylation variations in the whole blood of COVID-19 patients showed significant cell-type specificity with remarkable enrichment in gamma-delta T cells and presented a phenomenon of hypermethylation and low expression. Furthermore, we identified five genes whose methylation variations were associated with several cell types. Among them, S100A9, AHNAK, and CX3CR1 have been reported as potential COVID-19 biomarkers previously, and the others (TRAF3IP3 and LFNG) are closely associated with the immune and virus-related signaling pathways. We propose that they might serve as potential epigenetic biomarkers for COVID-19 and could play roles in important biological processes such as the immune response and antiviral activity.

Role of Epigenetics in Modulating Phenotypic Plasticity against Abiotic Stresses in Plants

Plants being sessile are always exposed to various environmental stresses, and to overcome these stresses, modifications at the epigenetic level can prove vital for their long-term survival. Epigenomics refers to the large-scale study of epigenetic marks on the genome, which include covalent modifications of histone tails (acetylation, methylation, phosphorylation, ubiquitination, and the small RNA machinery). Studies based on epigenetics have evolved over the years especially in understanding the mechanisms at transcriptional and posttranscriptional levels in plants against various environmental stimuli. Epigenomic changes in plants through induced methylation of specific genes that lead to changes in their expression can help to overcome various stress conditions. Recent studies suggested that epigenomics has a significant potential for crop improvement in plants. By the induction and modulation of various cellular processes like DNA methylation, histone modification, and biogenesis of noncoding RNAs, the plant genome can be activated which can help in achieving a quicker response against various plant stresses. Epigenetic modifications in plants allow them to adjust under varied environmental stresses by modulating their phenotypic plasticity and at the same time ensure the quality and yield of crops. The plasticity of the epigenome helps to adapt the plants during pre- and postdevelopmental processes. The variation in DNA methylation in different organisms exhibits variable phenotypic responses. The epigenetic changes also occur sequentially in the genome. Various studies indicated that environmentally stimulated epimutations produce variable responses especially in differentially methylated regions (DMR) that play a major role in the management of stress conditions in plants. Besides, it has been observed that environmental stresses cause specific changes in the epigenome that are closely associated with phenotypic modifications. However, the relationship between epigenetic modifications and phenotypic plasticity is still debatable. In this review, we will be discussing the role of various factors that allow epigenetic changes to modulate phenotypic plasticity against various abiotic stress in plants.

Integrated Analysis of ECT2 and COL17A1 as Potential Biomarkers for Pancreatic Cancer.

Background Pancreatic cancer (PC) is a malignant tumor of the digestive tract. It presents with atypical clinical symptoms and lacks specific diagnostic indicators. This study is aimed at exploring the potential biomarkers of PC. Methods TCGA database pancreatic cancer dataset was normalized and used to identify differentially expressed genes (DEGs). Survival, independent prognostic, and clinical correlation analyses were performed on DEGs to screen for key genes. DNA methylation, mutation, and copy number variation (CNV) analyses were used to analyze genetic variants in key genes. GSEA was performed to explore the functional enrichment of the key genes. Based on the expression of key genes, construction of a competing endogenous RNA (ceRNA) network, analysis of the tumor microenvironment (TME), and prediction of chemotherapeutic drug sensitivity were performed. Furthermore, the GEO database was used to validate the reliability of key genes. Results Two key genes (ECT2 and COL17A1) were identified, which were highly expressed in PC. The mRNA expression of ECT2 and COL17A1 was associated with DNA methylation and CNV. The cell cycle, proteasome, and pathways in cancer were enriched in the high-COL17A1 and ECT2 groups. The TME results showed that immune scores were decreased in the high-ECT2 group. CeRNA network results showed that there were eleven miRNAs were involved in the regulation of ECT2 and COL17A1. Moreover, pRRophetic analysis showed that 20 chemotherapeutic drugs were associated with ECT2 and COL17A1 expression. Conclusions Collectively, ECT2 and COL17A1 may be potential biomarkers for PC, providing a new direction for clinical diagnosis and treatment.

Genome-wide DNA methylation reveals potential epigenetic mechanism of age-dependent viral susceptibility in grass carp

BackgroundGrass carp are an important farmed fish in China that are infected by many pathogens, especially grass carp reovirus (GCRV). Notably, grass carp showed age-dependent susceptibility to GCRV; that is, grass carp not older than one year were sensitive to GCRV, while those over three years old were resistant to this virus. However, the underlying mechanism remains unclear. Herein, whole genome-wide DNA methylation and gene expression variations between susceptible five-month-old (FMO) and resistant three-year-old (TYO) grass carp were investigated aiming to uncover potential epigenetic mechanisms.ResultsColorimetric quantification revealed that the global methylation level in TYO fish was higher than that in FMO fish. Whole-genome bisulfite sequencing (WGBS) of the two groups revealed 6214 differentially methylated regions (DMRs) and 4052 differentially methylated genes (DMGs), with most DMRs and DMGs showing hypermethylation patterns in TYO fish. Correlation analysis revealed that DNA hypomethylation in promoter regions and DNA hypermethylation in gene body regions were associated with gene expression. Enrichment analysis revealed that promoter hypo-DMGs in TYO fish were significantly enriched in typical immune response pathways, whereas gene body hyper-DMGs in TYO fish were significantly enriched in terms related to RNA transcription, biosynthesis, and energy production. RNA-seq analysis of the corresponding samples indicated that most of the genes in the above terms were upregulated in TYO fish. Moreover, gene function analysis revealed that the two genes involved in energy metabolism displayed antiviral effects.ConclusionsCollectively, these results revealed genome-wide variations in DNA methylation between grass carp of different ages. DNA methylation and gene expression variations in genes involved in immune response, biosynthesis, and energy production may contribute to age-dependent susceptibility to GCRV in grass carp. Our results provide important information for disease-resistant breeding programs for grass carp and may also benefit research on age-dependent diseases in humans.

Social Factors, Epigenomics and Lupus in African American Women (SELA) Study: protocol for an observational mechanistic study examining the interplay of multiple individual and social factors on lupus outcomes in a health disparity population

IntroductionDespite the disproportional impact of SLE on historically marginalised communities, the individual and sociocultural factors underlying these health disparities remain elusive. We report the design and methods for a study...

Genetic and epigenetic variations in BDNF gene involved in Anorexia Nervosa

IntroductionAnorexia nervosa (AN) is a chronic psychiatric disorder resulting from abnormal eating habits with a high prevalence (0.5%). AN involves genetic and epigenetic factors supporting that AN is a metabo-psychiatric disorder. One candidate gene for AN, validated by meta-analyzes, is BDNF which encodes the brain-derived neurotrophic factor. BDNF negatively modulates the central control of food intake and its injection in rodents induces weight loss and anorexia. In humans, we observed an association of its functional variant Val66Met/rs6265 and electrodermal response to images of thinness suggesting an association between rs6265 and a reward effect of weight loss in AN.ObjectivesThis work study the impact of the functional polymorphism at risk rs6265, epigenetic variations in DNA methylation of BDNF gene and consequences on the concentrations of BDNF in AN patients.MethodsDNA was isolated from 24 AN patients and 48 controls. DNA methylation was measured for sites spanning the BDNF gene using Infinium HumanMethylation450 BeadChip technology. The genotyping of rs6265 was performed by Taqman-SNP assay. The BDNF was dosaged by ELISA from plasmas.ResultsWe observe that several sites are significantly hypermethylated in AN patients compared to controls. AN patients show significantly higher BDNF levels than controls. Finally, this BDNF concentration is significantly higher in AN carrying the risk Met66 allele.ConclusionsThis work demonstrates the effects of genetic and epigenetic variations of BDNF, which could constitute relevant diagnostic biomarkers of AN, and their likely consequences in the pathophysiology of AN. This work was supported by the Nestlé Foundation.DisclosureNo significant relationships.

Association of DNA methylation and genetic variations of the APOE gene with the risk of diabetic dyslipidemia.

Apolipoprotein E (APOE) is a key regulator of lipoprotein metabolism, and consequently, affects the plasma and tissue lipid contents. The aim of the present study was to investigate the parallel effects of APOE genetic variants and promoter methylation levels of six CpGs on the risk of diabetic dyslipidemia. A total of 204 Palestinian type 2 diabetes (T2D) patients (mean age ± SD: 62.7±10.2) were enrolled in the present study (n=96 with dyslipidemia and n=108 without dyslipidemia). Next generation sequencing was performed to analyze five single nucleotide polymorphisms: Two variants rs7412 and rs429358 that determine APOE ε alleles, and three variants in the promoter region (rs769446, rs449647, and rs405509). For all subjects, the most common genotype was ε3/ε3 (79.4%). No statistical differences were observed in the APOE ε polymorphisms and the three promoter variants among T2D patients with and without dyslipidemia (P>0.05). A comparison of lipid parameters between ε3/ε3 subjects and ε4 carriers in both groups revealed no significant differences in the mean values of LDL-C, HDL-C, TG, and TC levels (P>0.05). Six CpG sites in the APOE promoter on chromosome 19:44905755-44906078 were identified, and differential DNA methylation in these CpGs were observed between the study groups. Logistic regression analysis revealed a significant association of DNA methylation level at the six CpGs with an increased risk of diabetic dyslipidemia (odds ratio, 1.038; 95% confidence interval, 1.012-1.064; P=0.004). In conclusion, the present study revealed that DNA methylation levels in six CpGs in the APOE promoter region was associated with the risk of diabetic dyslipidemia independently of the APOE ε4 variant which could be a potential therapeutic target to reverse the methylation of the APOE promoter.

Abstract IA009: Cancer is a disease of epigenetic stochasticity

Abstract My colleagues and I proposed in 2006 (1) that increased epigenetic stochasticity is a driving force of tumor progression from its origin to metastasis and would allow rapid selection for tumor cell survival at the expense of the host. More recently (2), we have been pursuing an idea that natural selection will favor the emergence of genetic loci for epigenetic variance, not just mean, for loci in which the environment changes unpredictably but often enough, and these epigenetically variable loci are critical to normal embryonic development and injury response. The idea is also relevant to cancer, in that increased epigenetic stochasticity would allow rapid selection for tumor cell survival at the expense of the host. This model puts epigenetic instability at the heart of tumor progression and is the primary target of cancer mutations. Several recent observations from the laboratory point to a genome-scale disruption of the epigenome that involves large blocks of variable DNA methylation and chromatin affecting shores and islands and increasing gene expression variability (3,4). We also find large-scale reprogramming of chromatin and DNA modifications during the natural evolution of distant metastasis, with dependence on the oxidative branch of the pentose phosphate pathway (5). We have also developed mathematically rigorous Gibbs-Boltzmann-style epigenetic landscapes incorporating stochasticity, and we have shown its relationship to entropy in information theory (6). Recent data show that this approach identifies epigenetic and genetic drivers of cancer, using acute lymphoblastic leukemia as a model (7).

Evaluation of nanopore sequencing for epigenetic epidemiology: a comparison with DNA methylation microarrays.

Most epigenetic epidemiology to date has utilized microarrays to identify positions in the genome where variation in DNA methylation is associated with environmental exposures or disease. However, these profile less than 3% of DNA methylation sites in the human genome, potentially missing affected loci and preventing the discovery of disrupted biological pathways. Third generation sequencing technologies, including Nanopore sequencing, have the potential to revolutionize the generation of epigenetic data, not only by providing genuine genome-wide coverage but profiling epigenetic modifications direct from native DNA. Here we assess the viability of using Nanopore sequencing for epidemiology by performing a comparison with DNA methylation quantified using the most comprehensive microarray available, the Illumina EPIC array. We implemented a CRISPR-Cas9 targeted sequencing approach in concert with Nanopore sequencing to profile DNA methylation in three genomic regions to attempt to rediscover genomic positions that existing technologies have shown are differentially methylated in tobacco smokers. Using Nanopore sequencing reads, DNA methylation was quantified at 1779 CpGs across three regions, providing a finer resolution of DNA methylation patterns compared to the EPIC array. The correlation of estimated levels of DNA methylation between platforms was high. Furthermore, we identified 12 CpGs where hypomethylation was significantly associated with smoking status, including 10 within the AHRR gene. In summary, Nanopore sequencing is a valid option for identifying genomic loci where large differences in DNAm are associated with a phenotype and has the potential to advance our understanding of the role differential methylation plays in the etiology of complex disease.

Genetic and Epigenetic Differentiation Across Intertidal Gradients in the Foundation Plant Spartina alterniflora

Ecological genomics approaches have informed us about the structure of genetic diversity in natural populations that might underlie patterns in trait variation. However, we still know surprisingly little about the mechanisms that permit organisms to adapt to variable environmental conditions. The salt marsh foundation plant Spartina alterniflora exhibits a dramatic range in phenotype that is associated with a pronounced intertidal environmental gradient across a narrow spatial scale. Both genetic and non-genetic molecular mechanisms might underlie this phenotypic variation. To investigate both, we used epigenotyping-by-sequencing (epiGBS) to evaluate the make-up of natural populations across the intertidal environmental gradient. Based on recent findings, we expected that both DNA sequence and DNA methylation diversity would be explained by source population and habitat within populations. However, we predicted that epigenetic variation might be more strongly associated with habitat since similar epigenetic modifications could be rapidly elicited across different genetic backgrounds by similar environmental conditions. Overall, with PERMANOVA we found that population of origin explained a significant amount of the genetic (8.6%) and epigenetic (3.2%) variance. In addition, we found that a small but significant amount of genetic and epigenetic variance (<1%) was explained by habitat within populations. The interaction of population and habitat explained an additional 2.9% of the genetic variance and 1.4% of the epigenetic variance. By examining genetic and epigenetic variation within the same fragments (variation in close-cis), we found that population explained epigenetic variation in 9.2% of 8,960 tested loci, even after accounting for differences in the DNA sequence of the fragment. Habitat alone explained very little (<0.1%) of the variation in these close-cis comparisons, but the interaction of population and habitat explained 2.1% of the epigenetic variation in these loci. Using multiple matrix regression with randomization (MMRR) we found that phenotypic differences in natural populations were correlated with epigenetic and environmental differences even when accounting for genetic differences. Our results support the contention that sequence variation explains most of the variation in DNA methylation, but we have provided evidence that DNA methylation distinctly contributes to plant responses in natural populations.

Global DNA Methylation and mRNA-miRNA Variations Activated by Heat Shock Boost Early Microspore Embryogenesis in Cabbage (Brassica oleracea).

Microspore culture, a type of haploid breeding, is extensively used in the cultivation of cruciferous crops such as cabbage. Heat shock (HS) treatment is essential to improve the embryo rate during the culture process; however, its molecular role in boosting early microspore embryogenesis (ME) remains unknown. Here we combined DNA methylation levels, miRNAs, and transcriptome profiles in isolated microspores of cabbage ‘01-88’ under HS (32 °C for 24 h) and normal temperature (25 °C for 24 h) to investigate the regulatory roles of DNA methylation and miRNA in early ME. Global methylation levels were significantly different in the two pre-treatments, and 508 differentially methylated regions (DMRs) were identified; 59.92% of DMRs were correlated with transcripts, and 39.43% of miRNA locus were associated with methylation levels. Significantly, the association analysis revealed that 31 differentially expressed genes (DEGs) were targeted by methylation and miRNA and were mainly involved in the reactive oxygen species (ROS) response and abscisic acid (ABA) signaling, indicating that HS induced DNA methylation, and miRNA might affect ME by influencing ROS and ABA. This study revealed that DNA methylation and miRNA interfered with ME by modulating key genes and pathways, which could broaden our understanding of the molecular regulation of ME induced by HS pre-treatment.