High-resolution Nanopore methylome-maps reveal random hyper-methylation at CpG-poor regions as driver of chemoresistance in leukemias.

Aberrant DNA methylation at CpG dinucleotides is a cancer hallmark that is associated with the emergence of resistance to anti cancer treatment, though molecular mechanisms and biological significance remain elusive. Genome scale methylation maps by currently used methods are based on chemical modification of DNA and are best suited for analyses of methylation at CpG rich regions (CpG islands). We report the first high coverage whole-genome map in cancer using the long read nanopore technology, which allows simultaneous DNA-sequence and -methylation analyses on native DNA. We analyzed clonal epigenomic/genomic evolution in Acute Myeloid Leukemias (AMLs) at diagnosis and relapse, after chemotherapy. Long read sequencing coupled to a novel computational method allowed definition of differential methylation at unprecedented resolution, and showed that the relapse methylome is characterized by hypermethylation at both CpG islands and sparse CpGs regions. Most differentially methylated genes, however, were not differentially expressed nor enriched for chemoresistance genes. A small fraction of under-expressed and hyper-methylated genes at sparse CpGs, in the gene body, was significantly enriched in transcription factors (TFs). Remarkably, these few TFs supported large gene-regulatory networks including 50% of all differentially expressed genes in the relapsed AMLs and highly-enriched in chemoresistance genes. Notably, hypermethylated regions at sparse CpGs were poorly conserved in the relapsed AMLs, under-represented at their genomic positions and showed higher methylation entropy, as compared to CpG islands. Analyses of available datasets confirmed TF binding to their target genes and conservation of the same gene-regulatory networks in large patient cohorts. Relapsed AMLs carried few patient specific structural variants and DNA mutations, apparently not involved in drug resistance. Thus, drug resistance in AMLs can be mainly ascribed to the selection of random epigenetic alterations at sparse CpGs of a few transcription factors, which then induce reprogramming of the relapsing phenotype, independently of clonal genomic evolution.

Mitochondrial Uncoupling Induces Epigenome Remodeling and Promotes Differentiation in Neuroblastoma.

The Warburg effect is the major metabolic hallmark of cancer. According to Warburg himself, the consequence of the Warburg effect is cell dedifferentiation. Therefore, reversing the Warburg effect might be an approach to restore cell differentiation in cancer. In this study, we used a mitochondrial uncoupler, niclosamide ethanolamine (NEN), to activate mitochondrial respiration, which induced neural differentiation in neuroblastoma cells. NEN treatment increased the NAD+/NADH and pyruvate/lactate ratios and also the α-ketoglutarate/2-hydroxyglutarate (2-HG) ratio. Consequently, NEN treatment induced promoter CpG island demethylation and epigenetic landscape remodeling, activating the neural differentiation program. In addition, NEN treatment upregulated p53 but downregulated N-Myc and ß-catenin signaling in neuroblastoma cells. Importantly, even under hypoxia, NEN treatment remained effective in inhibiting 2-HG generation, promoting DNA demethylation, and suppressing hypoxia-inducible factor signaling. Dietary NEN intervention reduced tumor growth rate, 2-HG levels, and expression of N-Myc and ß-catenin in tumors in an orthotopic neuroblastoma mouse model. Integrative analysis indicated that NEN treatment upregulated favorable prognosis genes and downregulated unfavorable prognosis genes, which were defined using multiple neuroblastoma patient datasets. Altogether, these results suggest that mitochondrial uncoupling is an effective metabolic and epigenetic therapy for reversing the Warburg effect and inducing differentiation in neuroblastoma. SIGNIFICANCE: Targeting cancer metabolism using the mitochondrial uncoupler niclosamide ethanolamine leads to methylome reprogramming and differentiation in neuroblastoma, providing a therapeutic opportunity to reverse the Warburg effect and suppress tumor growth. See related commentary by Byrne and Bell, p.167.

Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism.

Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation; however, this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males and 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci; however, these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16% of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p value = 4.6e-13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

The Transcriptome and Epigenome Reveal Novel Changes in Transcription Regulation During Pancreatic Rat Islet Maturation.

Islet function is critical for normal glucose homeostasis. Unlike adult ß cells, fetal and neonatal islets are more proliferative and have decreased insulin secretion in response to stimuli. However, the underlying mechanisms governing functional maturity of islets have not been completely elucidated. Pancreatic islets comprise different cell types. The microenvironment of islets and interactions between these cell types are critical for ß-cell development and maturation. Thus, the study of intact islets is optimal to identify novel molecular mechanisms controlling islet functional development. Transcriptomes and genome-wide histone landscapes of H3K4me3, H3K27me3, and H3K27Ac from intact islets isolated from 2- and 10-week-old Sprague-Dawley rats were integrated to elucidate genes and pathways modulating islet development, as well as the contribution of epigenetic regulation. A total of 4489 differentially expressed genes were identified; 2289 and 2200 of them were up- and down-regulated in 10-week islets, respectively. Ingenuity Pathway Analysis revealed critical pathways regulating functional maturation of islets, including nutrient sensing, neuronal function, immune function, cell replication, and extracellular matrix. Furthermore, we identified significant changes in enrichment of H3K4me3, H3K27me3, and H3K27Ac marks, which correlated with expression changes of genes critical for islet function. These histone marks were enriched at critical transcription factor-binding motifs, such as Hoxa9, C/EBP-ß, Gata1, Foxo1, E2f1, E2f3, and Mafb. In addition, our chromatin immunoprecipitation sequencing data revealed multiple potential bivalent genes whose poised states changed with maturation. Collectively, our current study identified critical novel pathways for mature islet function and suggested a role for histone modifications in regulating islet development and maturation.

Identification of a methylomics-associated nomogram for predicting overall survival of stage I-II lung adenocarcinoma.

The aim of this paper was to identify DNA methylation based biomarkers for predicting overall survival (OS) of stage I-II lung adenocarcinoma (LUAD) patients. Methylation profile data of patients with stage I-II LUAD from The Cancer Genome Atlas (TCGA) database was used to determine methylation sites-based hallmark for stage I-II LUAD patients' OS. The patients were separated into training and validation datasets by using median risk score as cutoff. Univariate Cox, least absolute shrinkage and selection operator (LASSO) and multivariate Cox analyses were employed to develop a DNA methylation signature for OS of patients with stage I-II LUAD. As a result, an 11-DNA methylation signature was determined to be critically associated with the OS of patients with stage I-II LUAD. Analysis of receiver operating characteristics (ROC) suggested a high prognostic effectiveness of the 11-DNA methylation signature in patients with stage I-II LUAD (AUC at 1, 3, 5 years in training set were (0.849, 0.879, 0.831, respectively), validation set (0.742, 0.807, 0.904, respectively), entire TCGA dataset (0.747, 0.818, 0.870, respectively). Kaplan-Meier survival analyses exhibited that survival was significantly longer in the low-risk cohort compared to the high-risk cohort in the training dataset (P = 7e - 07), in the validation dataset (P = 1e - 08), and in the all-cohort dataset (P = 6e - 14). In addition, a nomogram was developed based on molecular factor (methylation risk score) as well as clinical factors (age and cancer status) (AUC at 1, 3, 5 years entire TCGA dataset were 0.770, 0.849, 0.979, respectively). The result verified that our methylomics-associated nomogram had a strong robustness for predicting stage I-II LUAD patients' OS. Furthermore, the nomogram combined clinical and molecular factors to determine an individualized probability of recurrence for patients with stage I-II LUAD, which stood for a major advance in the field of personalized medicine for pulmonary oncology. Collectively, we successfully identified a DNA methylation biomarker and a DNA methylation-based nomogram to predict the OS of patients with stage I-II LUAD.

Comprehensive Analysis of Differential m6A RNA Methylomes in the Hippocampus of Cocaine-Conditioned Mice.

N6-methyladenosine (m6A) is the most prevalent internal modification found in mRNAs and lncRNA and plays a vital role in posttranscriptional regulation in mammals. m6A is abundant in the nervous system, where it modulates neuronal development and hippocampus-dependent learning and memory. However, the roles of RNAs m6A modification and its related enzymes in cocaine reward are still not fully understood. In this study, we found that the fat mass and obesity-associated gene (FTO) demethylase, but not methyltransferase-like 3 (METTL3) and 14 (METTL14), was downregulated in the hippocampus following cocaine-induced conditioned place preference (CPP), and the level of m6A is notably higher in the hippocampus of cocaine CPP training mice. Using methylated m6A RNA immunoprecipitation sequencing (MeRIP-m6A-seq), we identified a total of 6516 m6A peaks within 4460 mRNAs, and 3083 m6A peaks within 850 lncRNAs were significantly dysregulated. Intriguingly, the altered m6A peaks within mRNAs and lncRNAs were enriched in synapse maturation and localization processes. Our study uncovers a critical role for an m6A epitranscriptomic dysregulation and downregulation of FTO expression in the hippocampus following cocaine-induced CPP.

Epigenome-wide study of brain DNA methylation following acute opioid intoxication.

BACKGROUND: Opioid abuse poses significant risk to individuals in the United States and epigenetic changes are a leading potential biomarker of opioid abuse. Current evidence, however, is mostly limited to candidate gene analysis in whole blood. To clarify the association between opioid abuse and DNA methylation, we conducted an epigenome-wide analysis of DNA methylation in brain samples of individuals who died from acute opioid intoxication and group-matched controls. METHODS: Tissue samples were extracted from the dorsolateral prefrontal cortex of 153 deceased individuals (Mage = 35.42; 62 % male; 77 % European ancestry). The study included 72 opioid samples, 53 psychiatric controls, and 28 normal controls. The epigenome-wide analysis was implemented using the Illumina MethylationEPIC BeadChip; analyses adjusted for sociodemographic characteristics, negative control principal components, ancestry principal components, cellular composition, and surrogate variables. Horvath's epigenetic age and Levine's PhenoAge were calculated, and gene set enrichment analyses were performed. RESULTS: Although no CpG sites survived false-discovery rate correction for multiple testing, 13 sites surpassed a relaxed significance threshold (p < 1.0 × 10-5). One of these sites was located within Netrin-1, a gene implicated in kappa opioid receptor activity. There was an association between opioid use and accelerated PhenoAge (b = 2.24, se = 1.11, p = .045). Gene set enrichment analyses revealed enrichment of differential methylation in GO and KEGG pathways broadly related to substance use. CONCLUSIONS: Netrin-1 may be associated with opioid overdose, and future research with larger samples across stages of opioid use will elucidate the complex genomics of opioid abuse.

Single-Cell Sequencing of Brain Cell Transcriptomes and Epigenomes.

Single-cell sequencing technologies, including transcriptomic and epigenomic assays, are transforming our understanding of the cellular building blocks of neural circuits. By directly measuring multiple molecular signatures in thousands to millions of individual cells, single-cell sequencing methods can comprehensively characterize the diversity of brain cell types. These measurements uncover gene regulatory mechanisms that shape cellular identity and provide insight into developmental and evolutionary relationships between brain cell populations. Single-cell sequencing data can aid the design of tools for targeted functional studies of brain circuit components, linking molecular signatures with anatomy, connectivity, morphology, and physiology. Here, we discuss the fundamental principles of single-cell transcriptome and epigenome sequencing, integrative computational analysis of the data, and key applications in neuroscience.

Epigenome-Wide Association Study of Thyroid Function Traits Identifies Novel Associations of fT3 With KLF9 and DOT1L.

CONTEXT: Circulating concentrations of free triiodothyronine (fT3), free thyroxine (fT4), and thyrotropin (TSH) are partly heritable traits. Recent studies have advanced knowledge of their genetic architecture. Epigenetic modifications, such as DNA methylation (DNAm), may be important in pituitary-thyroid axis regulation and action, but data are limited. OBJECTIVE: To identify novel associations between fT3, fT4, and TSH and differentially methylated positions (DMPs) in the genome in subjects from 2 Australian cohorts. METHOD: We performed an epigenome-wide association study (EWAS) of thyroid function parameters and DNAm using participants from: Brisbane Systems Genetics Study (median age 14.2 years, n = 563) and the Raine Study (median age 17.0 years, n = 863). Plasma fT3, fT4, and TSH were measured by immunoassay. DNAm levels in blood were assessed using Illumina HumanMethylation450 BeadChip arrays. Analyses employed generalized linear mixed models to test association between DNAm and thyroid function parameters. Data from the 2 cohorts were meta-analyzed. RESULTS: We identified 2 DMPs with epigenome-wide significant (P < 2.4E-7) associations with TSH and 6 with fT3, including cg00049440 in KLF9 (P = 2.88E-10) and cg04173586 in DOT1L (P = 2.09E-16), both genes known to be induced by fT3. All DMPs had a positive association between DNAm and TSH and a negative association between DNAm and fT3. There were no DMPs significantly associated with fT4. We identified 23 differentially methylated regions associated with fT3, fT4, or TSH. CONCLUSIONS: This study has demonstrated associations between blood-based DNAm and both fT3 and TSH. This may provide insight into mechanisms underlying thyroid hormone action and/or pituitary-thyroid axis function.

Methylome-wide association study of first-episode schizophrenia reveals a hypermethylated CpG site in the promoter region of the TNIK susceptibility gene.

Accumulating evidence suggests that epigenetics plays an important role in the etiology of schizophrenia. Here, we performed a methylome-wide association study (MWAS) of first-onset schizophrenia patients and controls from the Han Chinese population using microarray technology. The DNA methylation profiles revealed 4494 differentially methylated CpG sites. Gene ontology (GO) analysis showed that the functions of differentially methylated genes were primarily involved in enzymatic activity, cytoskeleton organization and cell adhesion, and the TNIK (encoding TRAF2- and NCK-interacting kinase) gene was enriched in most of these terms. By combining the MWAS results with those of previous genome-wide association studies (GWASs), we identified 72 candidate genes located in 49 human genome loci. Among the overlapping genes, the most significantly methylated CpG sites were in the transcriptional start site (TSS) 200 region (cg21413905, Punadjusted = 3.20 × 10-5) of TNIK. TNIK was listed in the top 50 differentially methylated loci. The results of pyrosequencing and TNIK mRNA expression were consistent with those of the microarray study. Bioinformatics analyses, dual-luciferase reporter assays and chromatin immunoprecipitation (ChIP) studies showed that TNIK interacted with genes associated with schizophrenia and NRF1 was identified as a novel transcription factor (TF) that binds to TNIK in its TSS200 region. Thus, the regulatory function of NRF1 may be influenced by the status of the methylated CpG site in this region. In summary, our study provides new insights into the epigenetic mechanisms that regulate schizophrenia. Studies of the functions of TNIK methylation should be performed in vitro and in vivo to provide a better understanding of the pathophysiology of schizophrenia.

Associations of Alcohol Consumption With Epigenome-Wide DNA Methylation and Epigenetic Age Acceleration: Individual-Level and Co-twin Comparison Analyses.

BACKGROUND: DNA methylation may play a role in the progression from normative to problematic drinking and underlie adverse health outcomes associated with alcohol misuse. We examined the association between alcohol consumption and DNA methylation patterns using 3 approaches: a conventional epigenome-wide association study (EWAS); a co-twin comparison design, which controls for genetic and environmental influences that twins share; and a regression of age acceleration, defined as a discrepancy between chronological age and DNA methylation age, on alcohol consumption. METHODS: Participants came from the Finnish Twin Cohorts (FinnTwin12/FinnTwin16; N = 1,004; 55% female; average age = 23 years). Individuals reported the number of alcoholic beverages consumed in the past week, and epigenome-wide DNA methylation was assessed in whole blood using the Infinium HumanMethylation450 BeadChip. RESULTS: In the EWAS, alcohol consumption was significantly related to methylation at 24 CpG sites. When evaluating whether differences between twin siblings (185 monozygotic pairs) in alcohol consumption predicted differences in DNA methylation, co-twin comparisons replicated 4 CpG sites from the EWAS and identified 23 additional sites. However, when we examined qualitative differences in drinking patterns between twins (heavy drinker vs. light drinker/abstainer or moderate drinker vs. abstainer; 44 pairs), methylation patterns did not significantly differ within twin pairs. Finally, individuals who reported higher alcohol consumption also exhibited greater age acceleration, though results were no longer significant after controlling for genetic and environmental influences shared by co-twins. CONCLUSIONS: Our analyses offer insight into the associations between epigenetic variation and levels of alcohol consumption in young adulthood.

Roles of DEMETER in regulating DNA methylation in vegetative tissues and pathogen resistance.

DNA methylation is an epigenetic mark important for genome stability and gene expression. In Arabidopsis thaliana, the 5-methylcytosine DNA glycosylase/demethylase DEMETER (DME) controls active DNA demethylation during the reproductive stage; however, the lethality of loss-of-function dme mutations has made it difficult to assess DME function in vegetative tissues. Here, we edited DME using clustered regularly interspaced short palindromic repeats (CRISPR) /CRISPR-associated protein 9 and created three weak dme mutants that produced a few viable seeds. We also performed central cell-specific complementation in a strong dme mutant and combined this line with mutations in the other three Arabidopsis demethylase genes to generate the dme ros1 dml2 dml3 (drdd) quadruple mutant. A DNA methylome analysis showed that DME is required for DNA demethylation at hundreds of genomic regions in vegetative tissues. A transcriptome analysis of the drdd mutant revealed that DME and the other three demethylases are important for plant responses to biotic and abiotic stresses in vegetative tissues. Despite the limited role of DME in regulating DNA methylation in vegetative tissues, the dme mutants showed increased susceptibility to bacterial and fungal pathogens. Our study highlights the important functions of DME in vegetative tissues and provides valuable genetic tools for future investigations of DNA demethylation in plants.

Epigenetic mechanisms within the sperm epigenome and their diagnostic potential.

The sperm epigenome contains a highly unique and specialized epigenetic landscape. Insightful questions need be asked about these epigenetic signatures and their predictive potential to assess the approximately 1 in 6 couples who experience infertility. Among those couples that do experience infertility, approximately half of the cases involve a male factor. Unfortunately, there is a significant lack of effective diagnostic tools in the male infertility space and thus clinicians are left with little data upon which they can formulate data driven treatment plans. Taking together this information and the striking prevalence of male infertility it's obvious that there is a need for improved diagnostic techniques for male infertility. Many studies have identified what appear to be clinically meaningful epigenetic alterations in sperm that may add utility in the diagnoses of infertility and improvement of pregnancy outcomes. Many researchers believe that continued analysis of these various epigenetic mechanisms may provide powerful predictive insight. In fact, there is promising current data suggesting that the predictive power of DNA methylation, Nuclear Proteins, and miRNA signatures in sperm likely can improve what is currently found with traditional diagnosis of male infertility. The focus of this review is to give a brief understanding to the field of epigenetics and the potential predictive power the sperm epigenome may hold in relation to improving the treatment and diagnosis of male infertility patients.

Maternal dysglycaemia, changes in the infant's epigenome modified with a diet and physical activity intervention in pregnancy: Secondary analysis of a randomised control trial.

BACKGROUND: Higher maternal plasma glucose (PG) concentrations, even below gestational diabetes mellitus (GDM) thresholds, are associated with adverse offspring outcomes, with DNA methylation proposed as a mediating mechanism. Here, we examined the relationships between maternal dysglycaemia at 24 to 28 weeks' gestation and DNA methylation in neonates and whether a dietary and physical activity intervention in pregnant women with obesity modified the methylation signatures associated with maternal dysglycaemia. METHODS AND FINDINGS: We investigated 557 women, recruited between 2009 and 2014 from the UK Pregnancies Better Eating and Activity Trial (UPBEAT), a randomised controlled trial (RCT), of a lifestyle intervention (low glycaemic index (GI) diet plus physical activity) in pregnant women with obesity (294 contol, 263 intervention). Between 27 and 28 weeks of pregnancy, participants had an oral glucose (75 g) tolerance test (OGTT), and GDM diagnosis was based on diagnostic criteria recommended by the International Association of Diabetes and Pregnancy Study Groups (IADPSG), with 159 women having a diagnosis of GDM. Cord blood DNA samples from the infants were interrogated for genome-wide DNA methylation levels using the Infinium Human MethylationEPIC BeadChip array. Robust regression was carried out, adjusting for maternal age, smoking, parity, ethnicity, neonate sex, and predicted cell-type composition. Maternal GDM, fasting glucose, 1-h, and 2-h glucose concentrations following an OGTT were associated with 242, 1, 592, and 17 differentially methylated cytosine-phosphate-guanine (dmCpG) sites (false discovery rate (FDR) ≤ 0.05), respectively, in the infant's cord blood DNA. The most significantly GDM-associated CpG was cg03566881 located within the leucine-rich repeat-containing G-protein coupled receptor 6 (LGR6) (FDR = 0.0002). Moreover, we show that the GDM and 1-h glucose-associated methylation signatures in the cord blood of the infant appeared to be attenuated by the dietary and physical activity intervention during pregnancy; in the intervention arm, there were no GDM and two 1-h glucose-associated dmCpGs, whereas in the standard care arm, there were 41 GDM and 160 1-h glucose-associated dmCpGs. A total of 87% of the GDM and 77% of the 1-h glucose-associated dmCpGs had smaller effect sizes in the intervention compared to the standard care arm; the adjusted r2 for the association of LGR6 cg03566881 with GDM was 0.317 (95% confidence interval (CI) 0.012, 0.022) in the standard care and 0.240 (95% CI 0.001, 0.015) in the intervention arm. Limitations included measurement of DNA methylation in cord blood, where the functional significance of such changes are unclear, and because of the strong collinearity between treatment modality and severity of hyperglycaemia, we cannot exclude that treatment-related differences are potential confounders. CONCLUSIONS: Maternal dysglycaemia was associated with significant changes in the epigenome of the infants. Moreover, we found that the epigenetic impact of a dysglycaemic prenatal maternal environment appeared to be modified by a lifestyle intervention in pregnancy. Further research will be needed to investigate possible medical implications of the findings. TRIAL REGISTRATION: ISRCTN89971375.

Rebooting the Epigenomes during Mammalian Early Embryogenesis.

Upon fertilization, terminally differentiated gametes are transformed to a totipotent zygote, which gives rise to an embryo. How parental epigenetic memories are inherited and reprogrammed to accommodate parental-to-zygotic transition remains a fundamental question in developmental biology, epigenetics, and stem cell biology. With the rapid advancement of ultra-sensitive or single-cell epigenome analysis methods, unusual principles of epigenetic reprogramming begin to be unveiled. Emerging data reveal that in many species, the parental epigenome undergoes dramatic reprogramming followed by subsequent re-establishment of the embryo epigenome, leading to epigenetic "rebooting." Here, we discuss recent progress in understanding epigenetic reprogramming and their functions during mammalian early development. We also highlight the conserved and species-specific principles underlying diverse regulation of the epigenome in early embryos during evolution.

Exploration of the sputum methylome and omics deconvolution by quadratic programming in molecular profiling of asthma and COPD: the road to sputum omics 2.0.

BACKGROUND: To date, most studies involving high-throughput analyses of sputum in asthma and COPD have focused on identifying transcriptomic signatures of disease. No whole-genome methylation analysis of sputum cells has been performed yet. In this context, the highly variable cellular composition of sputum has potential to confound the molecular analyses. METHODS: Whole-genome transcription (Agilent Human 4 × 44 k array) and methylation (Illumina 450 k BeadChip) analyses were performed on sputum samples of 9 asthmatics, 10 healthy and 10 COPD subjects. RNA integrity was checked by capillary electrophoresis and used to correct in silico for bias conferred by RNA degradation during biobank sample storage. Estimates of cell type-specific molecular profiles were derived via regression by quadratic programming based on sputum differential cell counts. All analyses were conducted using the open-source R/Bioconductor software framework. RESULTS: A linear regression step was found to perform well in removing RNA degradation-related bias among the main principal components of the gene expression data, increasing the number of genes detectable as differentially expressed in asthma and COPD sputa (compared to controls). We observed a strong influence of the cellular composition on the results of mixed-cell sputum analyses. Exemplarily, upregulated genes derived from mixed-cell data in asthma were dominated by genes predominantly expressed in eosinophils after deconvolution. The deconvolution, however, allowed to perform differential expression and methylation analyses on the level of individual cell types and, though we only analyzed a limited number of biological replicates, was found to provide good estimates compared to previously published data about gene expression in lung eosinophils in asthma. Analysis of the sputum methylome indicated presence of differential methylation in genomic regions of interest, e.g. mapping to a number of human leukocyte antigen (HLA) genes related to both major histocompatibility complex (MHC) class I and II molecules in asthma and COPD macrophages. Furthermore, we found the SMAD3 (SMAD family member 3) gene, among others, to lie within differentially methylated regions which has been previously reported in the context of asthma. CONCLUSIONS: In this methodology-oriented study, we show that methylation profiling can be easily integrated into sputum analysis workflows and exhibits a strong potential to contribute to the profiling and understanding of pulmonary inflammation. Wherever RNA degradation is of concern, in silico correction can be effective in improving both sensitivity and specificity of downstream analyses. We suggest that deconvolution methods should be integrated in sputum omics analysis workflows whenever possible in order to facilitate the unbiased discovery and interpretation of molecular patterns of inflammation.

The role of the hypothalamus and pituitary epigenomes in central activation of the reproductive axis at puberty.

Puberty is programmed through a multifactorial gene network which works to activate the pulsatile secretion of the gonadotropin releasing hormone (GnRH), and subsequently elevate circulating levels of the pituitary gonadotropins that stimulate gonadal activity. Although this developmental transition normally occurs at a limited age-range in individuals of the same genetic background and environment, pubertal onset can occur prematurely or be delayed following changes in ambient conditions, or due to genetic variations or mutations, many of which have remained elusive due to their location in distal regulatory elements. Growing evidence is pointing to a pivotal role for the epigenome in regulating key genes in the reproductive hypothalamus and pituitary at this time, which might mediate some of the plasticity of pubertal timing. This review will address epigenetic mechanisms which have been demonstrated in the KNDy neurons that increase the output of pulsatile GnRH, and those involved in activation of the GnRH gene and its receptor, and describes how GnRH utilizes epigenetic mechanisms to stimulate transcription of the pituitary gonadotropin genes in the context of the chromatin landscape.

Transcriptome and epigenome analyses of vernalization in Arabidopsis thaliana.

Vernalization accelerates flowering after prolonged winter cold. Transcriptional and epigenetic changes are known to be involved in the regulation of the vernalization response. Despite intensive applications of next-generation sequencing in diverse aspects of plant research, genome-wide transcriptome and epigenome profiling during the vernalization response has not been conducted. In this work, to our knowledge, we present the first comprehensive analyses of transcriptomic and epigenomic dynamics during the vernalization process in Arabidopsis thaliana. Six major clusters of genes exhibiting distinctive features were identified. Temporary changes in histone H3K4me3 levels were observed that likely coordinate photosynthesis and prevent oxidative damage during cold exposure. In addition, vernalization induced a stable accumulation of H3K27me3 over genes encoding many development-related transcription factors, which resulted in either inhibition of transcription or a bivalent status of the genes. Lastly, FLC-like and VIN3-like genes were identified that appear to be novel components of the vernalization pathway.

Role of the Epigenome in Heart Failure.

Gene expression is needed for the maintenance of heart function under normal conditions and in response to stress. Each cell type of the heart has a specific program controlling transcription. Different types of stress induce modifications of these programs and, if prolonged, can lead to altered cardiac phenotype and, eventually, to heart failure. The transcriptional status of a gene is regulated by the epigenome, a complex network of DNA and histone modifications. Until a few years ago, our understanding of the role of the epigenome in heart disease was limited to that played by histone deacetylation. But over the last decade, the consequences for the maintenance of homeostasis in the heart and for the development of cardiac hypertrophy of a number of other modifications, including DNA methylation and hydroxymethylation, histone methylation and acetylation, and changes in chromatin architecture, have become better understood. Indeed, it is now clear that many levels of regulation contribute to defining the epigenetic landscape required for correct cardiomyocyte function, and that their perturbation is responsible for cardiac hypertrophy and fibrosis. Here, we review these aspects and draw a picture of what epigenetic modification may imply at the therapeutic level for heart failure.

Distinct epigenomic and transcriptomic modifications associated with Wolbachia-mediated asexuality.

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