Dynamic DNA methylation: In the right place at the right time

The synovial and blood monocyte DNA methylomes mirror prognosis, evolution and treatment in early arthritis

Identifying predictive biomarkers at early stages of inflammatory arthritis is crucial for starting appropriate therapies to avoid poor outcomes. Monocytes (MOs) and macrophages, largely associated with arthritis, are contributors and sensors of inflammation through epigenetic modifications. In this study, we investigated associations between clinical features and DNA methylation in blood and synovial fluid (SF) MOs in a prospective cohort of patients with early inflammatory arthritis. DNA methylation profiles of undifferentiated arthritis (UA) blood MOs exhibited marked alterations in comparison with those from healthy donors. We identified additional differences both in blood and SF MOs after comparing patients with UA grouped by their future outcomes, i.e., good versus poor. Patient profiles in subsequent visits revealed a reversion toward a healthy level in both groups, those requiring disease-modifying antirheumatic drugs and those who remitted spontaneously. Changes in disease activity between visits also affected DNA methylation, which was partially concomitant in the SF of UA and in blood MOs of patients with rheumatoid arthritis. Epigenetic similarities between arthritis types allow a common prediction of disease activity. Our results constitute a resource of DNA methylation–based biomarkers of poor prognosis, disease activity, and treatment efficacy for the personalized clinical management of early inflammatory arthritis.

Gene-body methylation in plants: mechanisms, functions and important implications for understanding evolutionary processes

Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants.

Integrating whole genome sequencing, methylation, gene expression, topological associated domain information in regulatory mutation prediction: a study of follicular lymphoma

A major challenge in human genetics is of the analysis of the interplay between genetic and epigenetic factors in a multifactorial disease like cancer. Here, a novel methodology is proposed to investigate genome-wide regulatory mechanisms in cancer, as studied with the example of follicular Lymphoma (FL). In a first phase, a new machine-learning method is designed to identify Differentially Methylated Regions (DMRs) by computing six attributes. In a second phase, an integrative data analysis method is developed to study regulatory mutations in FL, by considering differential methylation information together with DNA sequence variation, differential gene expression, 3D organization of genome (e.g., topologically associated domains), and enriched biological pathways. Resulting mutation block-gene pairs are further ranked to find out the significant ones. By this approach, BCL2 and BCL6 were identified as top-ranking FL-related genes with several mutation blocks and DMRs acting on their regulatory regions. Two additional genes, CDCA4 and CTSO, were also found in top rank with significant DNA sequence variation and differential methylation in neighboring areas, pointing towards their potential use as biomarkers for FL. This work combines both genomic and epigenomic information to investigate genome-wide gene regulatory mechanisms in cancer and contribute to devising novel treatment strategies.

Multi-omic characterization of genome-wide abnormal DNA methylation reveals diagnostic and prognostic markers for esophageal squamous-cell carcinoma

This study investigates aberrant DNA methylations as potential diagnosis and prognosis markers for esophageal squamous-cell carcinoma (ESCC), which if diagnosed at advanced stages has <30% five-year survival rate. Comparing genome-wide methylation sites of 91 ESCC and matched adjacent normal tissues, we identified 35,577 differentially methylated CpG sites (DMCs) and characterized their distribution patterns. Integrating whole-genome DNA and RNA-sequencing data of the same samples, we found multiple dysregulated transcription factors and ESCC-specific genomic correlates of identified DMCs. Using featured DMCs, we developed a 12-marker diagnostic panel with high accuracy in our dataset and the TCGA ESCC dataset, and a 4-marker prognostic panel distinguishing high-risk patients. In-vitro experiments validated the functions of 4 marker host genes. Together these results provide additional evidence for the important roles of aberrant DNA methylations in ESCC development and progression. Our DMC-based diagnostic and prognostic panels have potential values for clinical care of ESCC, laying foundations for developing targeted methylation assays for future non-invasive cancer detection methods.

Higher Expression of DNA (de)methylation-Related Genes Reduces Adipogenicity in Dental Pulp Stem Cells

Obesity is a significant health concern that has reached alarming proportions worldwide. The overconsumption of high-energy foods may cause metabolic dysfunction and promote the generation of new adipocytes by contributing to several obesity-related diseases. Such concerns demand a deeper understanding of the origin of adipocytes if we want to develop new therapeutic approaches. Recent findings indicate that adipocyte development is facilitated by tight epigenetic reprogramming, which is required to activate the gene program to change the fate of mesenchymal stem cells (MSCs) into mature adipocytes. Like adipose tissue, different tissues are also potential sources of adipocyte-generating MSCs, so it is interesting to explore whether the epigenetic mechanisms of adipogenic differentiation vary from one depot to another. To investigate how DNA methylation (an epigenetic mark that plays an essential role in controlling transcription and cellular differentiation) contributes to adipogenic potential, dental pulp stem cells (DPSCs) and periodontal ligament stem cells (PLSCs) were analyzed during adipogenic differentiation in vitro. Here, we show that the capacity to differentiate from DPSCs or PLSCs to adipocytes may be associated with the expression pattern of DNA methylation-related genes acquired during the induction of the adipogenic program. Our study provides insights into the details of DNA methylation during the adipogenic determination of dental stem cells, which can be a starting point to identify the factors that affect the differentiation of these cells and provide new strategies to regulate differentiation and adipocyte expansion.

scMelody: An Enhanced Consensus-Based Clustering Model for Single-Cell Methylation Data by Reconstructing Cell-to-Cell Similarity

Single-cell DNA methylation sequencing technology has brought new perspectives to investigate epigenetic heterogeneity, supporting a need for computational methods to cluster cells based on single-cell methylation profiles. Although several methods have been developed, most of them cluster cells based on single (dis)similarity measures, failing to capture complete cell heterogeneity and resulting in locally optimal solutions. Here, we present scMelody, which utilizes an enhanced consensus-based clustering model to reconstruct cell-to-cell methylation similarity patterns and identifies cell subpopulations with the leveraged information from multiple basic similarity measures. Besides, benefitted from the reconstructed cell-to-cell similarity measure, scMelody could conveniently leverage the clustering validation criteria to determine the optimal number of clusters. Assessments on distinct real datasets showed that scMelody accurately recapitulated methylation subpopulations and outperformed existing methods in terms of both cluster partitions and the number of clusters. Moreover, when benchmarking the clustering stability of scMelody on a variety of synthetic datasets, it achieved significant clustering performance gains over existing methods and robustly maintained its clustering accuracy over a wide range of number of cells, number of clusters and CpG dropout proportions. Finally, the real case studies demonstrated the capability of scMelody to assess known cell types and uncover novel cell clusters.

Analysis of 5-Methylcytosine Regulators and DNA Methylation-Driven Genes in Colon Cancer

Background: Epigenetic-driven events are important molecular mechanisms of carcinogenesis. The 5-methylcytosine (5mC) regulators play important roles in the methylation-driven gene expression. However, the effect of the 5mC regulators on the oncogenic pathways in colon cancer (CC) remains unclear. Also, the clinical value of such epigenetic-driven events needs further research.

Methods: The transcriptome and matching epigenetic data were obtained from The Cancer Genome Atlas dataset. The gene set variation analysis identified the oncogenic pathways adjusted by 5mC regulators. The “edgeR” and “methylmix” package identified the differential expression genes of DNA methylation-driven genes. The correlation between 5mC regulators or transcription factors and shortlisted genes was investigated by calculating the Spearman's rank correlation coefficient. Among them, the genes related to diagnosis were screened out based on differential gene expression in extracellular vesicles (EVs) by the “limma” package and histology by immunohistochemistry. Then, a risk signature was constructed by fitting the generalized linear model and validated by the receiver operating characteristic curve.

Results: MYC targets pathway and phosphatidylinositol-3-kinase–AKT–mammalian target of rapamycin signaling pathway were identified as the hallmark-related pathways associated with 5mC regulators. Also, the P53 pathway was subject to the influence of regulators' expression. A five methylation-driven gene signature (FIRRE, MYBL2, TGFBI, AXIN2, and SLC35D3) was developed as the biomarker for CC diagnosis. Meanwhile, those genes positively related to 5mC regulators and interacted with their relevant or transcription factors.

Conclusion: In general, 5mC regulators are positively related to each other and DNA methylation-driven genes, with the relationship of multiple active and inhibitory pathways related to cancer. Meanwhile, the signature (FIRRE, MYBL2, TGFBI, AXIN2, and SLC35D3) can prefigure prospective diagnosis in CC.

Potential Role of Domains Rearranged Methyltransferase7 in Starch and Chlorophyll Metabolism to Regulate Leaf Senescence in Tomato

Deoxyribonucleic acid (DNA) methylation is an important epigenetic mark involved in diverse biological processes. Here, we report the critical function of tomato (Solanum lycopersicum) Domains Rearranged Methyltransferase7 (SlDRM7) in plant growth and development, especially in leaf interveinal chlorosis and senescence. Using a hairpin RNA-mediated RNA interference (RNAi), we generated SlDRM7-RNAi lines and observed pleiotropic developmental defects including small and interveinal chlorosis leaves. Combined analyses of whole genome bisulfite sequence (WGBS) and RNA-seq revealed that silencing of SlDRM7 caused alterations in both methylation levels and transcript levels of 289 genes, which are involved in chlorophyll synthesis, photosynthesis, and starch degradation. Furthermore, the photosynthetic capacity decreased in SlDRM7-RNAi lines, consistent with the reduced chlorophyll content and repression of genes involved in chlorophyll biosynthesis, photosystem, and photosynthesis. In contrast, starch granules were highly accumulated in chloroplasts of SlDRM7-RNAi lines and associated with lowered expression of genes in the starch degradation pathway. In addition, SlDRM7 was activated by aging- and dark-induced senescence. Collectively, these results demonstrate that SlDRM7 acts as an epi-regulator to modulate the expression of genes related to starch and chlorophyll metabolism, thereby affecting leaf chlorosis and senescence in tomatoes.

Antifibrotic factor KLF4 is repressed by the miR-10/TFAP2A/TBX5 axis in dermal fibroblasts: insights from twins discordant for systemic sclerosis

OBJECTIVES

Systemic sclerosis (SSc) is a complex disease of unknown aetiology in which inflammation and fibrosis lead to multiple organ damage. There is currently no effective therapy that can halt the progression of fibrosis or reverse it, thus studies that provide novel insights into disease pathogenesis and identify novel potential therapeutic targets are critically needed.

METHODS

We used global gene expression and genome-wide DNA methylation analyses of dermal fibroblasts (dFBs) from a unique cohort of twins discordant for SSc to identify molecular features of this pathology. We validated the findings using in vitro, ex vivo and in vivo models.

RESULTS

Our results revealed distinct differentially expressed and methylated genes, including several transcription factors involved in stem cell differentiation and developmental programmes (KLF4, TBX5, TFAP2A and homeobox genes) and the microRNAs miR-10a and miR-10b which target several of these deregulated genes. We show that KLF4 expression is reduced in SSc dFBs and its expression is repressed by TBX5 and TFAP2A. We also show that KLF4 is antifibrotic, and its conditional knockout in fibroblasts promotes a fibrotic phenotype.

CONCLUSIONS

Our data support a role for epigenetic dysregulation in mediating SSc susceptibility in dFBs, illustrating the intricate interplay between CpG methylation, miRNAs and transcription factors in SSc pathogenesis, and highlighting the potential for future use of epigenetic modifiers as therapies.

The biology of time: dynamic responses of cell types to developmental, circadian and environmental cues

As sessile organisms, plants are finely tuned to respond dynamically to developmental, circadian and environmental cues. Genome-wide studies investigating these types of cues have uncovered the intrinsically different ways they can impact gene expression over time. Recent advances in single-cell sequencing and time-based bioinformatic algorithms are now beginning to reveal the dynamics of these time-based responses within individual cells and plant tissues. Here, we review what these techniques have revealed about the spatiotemporal nature of gene regulation, paying particular attention to the three distinct ways in which plant tissues are time sensitive. (i) First, we discuss how studying plant cell identity can reveal developmental trajectories hidden in pseudotime. (ii) Next, we present evidence that indicates that plant cell types keep their own local time through tissue-specific regulation of the circadian clock. (iii) Finally, we review what determines the speed of environmental signaling responses, and how they can be contingent on developmental and circadian time. By these means, this review sheds light on how these different scales of time-based responses can act with tissue and cell-type specificity to elicit changes in whole plant systems.

DNA Epigenetics in Addiction Susceptibility

Addiction is a chronically relapsing neuropsychiatric disease that occurs in some, but not all, individuals who use substances of abuse. Relatively little is known about the mechanisms which contribute to individual differences in susceptibility to addiction. Neural gene expression regulation underlies the pathogenesis of addiction, which is mediated by epigenetic mechanisms, such as DNA modifications. A growing body of work has demonstrated distinct DNA epigenetic signatures in brain reward regions that may be associated with addiction susceptibility. Furthermore, factors that influence addiction susceptibility are also known to have a DNA epigenetic basis. In the present review, we discuss the notion that addiction susceptibility has an underlying DNA epigenetic basis. We focus on major phenotypes of addiction susceptibility and review evidence of cell type-specific, time dependent, and sex biased effects of drug use. We highlight the role of DNA epigenetics in these diverse processes and propose its contribution to addiction susceptibility differences. Given the prevalence and lack of effective treatments for addiction, elucidating the DNA epigenetic mechanism of addiction vulnerability may represent an expeditious approach to relieving the addiction disease burden.

DNA Methylation-Dependent Gene Expression Regulation of Glutamate Transporters in Cultured Radial Glial Cells

Exposure to xenobiotics has a significant impact in brain physiology that could be liked to an excitotoxic process induced by a massive release of the main excitatory neurotransmitter, L-glutamate. Overstimulation of extra-synaptic glutamate receptors, mainly of the N-methyl-D-aspartate subtype leads to a disturbance of intracellular calcium homeostasis that is critically involved in neuronal death. Hence, glutamate extracellular levels are tightly regulated through its uptake by glial glutamate transporters. It has been observed that glutamate regulates its own removal, both in the short-time frame via a transporter-mediated decrease in the uptake, and in the long-term through the transcriptional control of its gene expression, a process mediated by glutamate receptors that involves the Ca²⁺/diacylglycerol-dependent protein kinase and the transcription factor Ying Yang 1. Taking into consideration that this transcription factor is a member of the Polycomb complex and thus, part of repressive and activating chromatin remodeling factors, it might direct the interaction of DNA methyltransferases or dioxygenases of methylated cytosines to their target sequences. Here we explored the role of dynamic DNA methylation in the expression and function of glial glutamate transporters. To this end, we used the well-characterized models of primary cultures of chick cerebellar Bergmann glia cells and a human retina-derived Müller glia cell line. A time and dose-dependent increase in global DNA methylation was evident upon glutamate exposure. Under hypomethylation conditions, the glial glutamate transporter protein levels and uptake activity were increased. These results favor the notion that a dynamic DNA methylation program triggered by glutamate in glial cells modulates one of its major functions: glutamate removal.

The CLASSY family controls tissue-specific DNA methylation patterns in Arabidopsis

DNA methylation shapes the epigenetic landscape of the genome, plays critical roles in regulating gene expression, and ensures transposon silencing. As is evidenced by the numerous defects associated with aberrant DNA methylation landscapes, establishing proper tissue-specific methylation patterns is critical. Yet, how such differences arise remains a largely open question in both plants and animals. Here we demonstrate that CLASSY1-4 (CLSY1-4), four locus-specific regulators of DNA methylation, also control tissue-specific methylation patterns, with the most striking pattern observed in ovules where CLSY3 and CLSY4 control DNA methylation at loci with a highly conserved DNA motif. On a more global scale, we demonstrate that specific clsy mutants are sufficient to shift the epigenetic landscape between tissues. Together, these findings reveal substantial epigenetic diversity between tissues and assign these changes to specific CLSY proteins, elucidating how locus-specific targeting combined with tissue-specific expression enables the CLSYs to generate epigenetic diversity during plant development.

CLASSY (CLSY) proteins regulate DNA methylation at specific loci in the Arabidopsis genome. Here the authors show that the CLSYs also control tissue-specific DNA methylation, including at siren loci in ovules, and that the lack of an individual CLSYs can shift the epigenetic landscape between tissues.

Developmental and Injury-induced Changes in DNA Methylation in Regenerative versus Non-regenerative Regions of the Vertebrate Central Nervous System

BACKGROUND

Because some of its CNS neurons (e.g., retinal ganglion cells after optic nerve crush (ONC)) regenerate axons throughout life, whereas others (e.g., hindbrain neurons after spinal cord injury (SCI)) lose this capacity as tadpoles metamorphose into frogs, the South African claw-toed frog, Xenopus laevis, offers unique opportunities for exploring differences between regenerative and non-regenerative responses to CNS injury within the same organism. An earlier, three-way RNA-seq study (frog ONC eye, tadpole SCI hindbrain, frog SCI hindbrain) identified genes that regulate chromatin accessibility among those that were differentially expressed in regenerative vs non-regenerative CNS [11]. The current study used whole genome bisulfite sequencing (WGBS) of DNA collected from these same animals at the peak period of axon regeneration to study the extent to which DNA methylation could potentially underlie differences in chromatin accessibility between regenerative and non-regenerative CNS.

RESULTS

Consistent with the hypothesis that DNA of regenerative CNS is more accessible than that of non-regenerative CNS, DNA from both the regenerative tadpole hindbrain and frog eye was less methylated than that of the non-regenerative frog hindbrain. Also, consistent with observations of CNS injury in mammals, DNA methylation in non-regenerative frog hindbrain decreased after SCI. However, contrary to expectations that the level of DNA methylation would decrease even further with axotomy in regenerative CNS, DNA methylation in these regions instead increased with injury. Injury-induced differences in CpG methylation in regenerative CNS became especially enriched in gene promoter regions, whereas non-CpG methylation differences were more evenly distributed across promoter regions, intergenic, and intragenic regions. In non-regenerative CNS, tissue-related (i.e., regenerative vs. non-regenerative CNS) and injury-induced decreases in promoter region CpG methylation were significantly correlated with increased RNA expression, but the injury-induced, increased CpG methylation seen in regenerative CNS across promoter regions was not, suggesting it was associated with increased rather than decreased chromatin accessibility. This hypothesis received support from observations that in regenerative CNS, many genes exhibiting increased, injury-induced, promoter-associated CpG-methylation also exhibited increased RNA expression and association with histone markers for active promoters and enhancers. DNA immunoprecipitation for 5hmC in optic nerve regeneration found that the promoter-associated increases seen in CpG methylation were distinct from those exhibiting changes in 5hmC.

CONCLUSIONS

Although seemingly paradoxical, the increased injury-associated DNA methylation seen in regenerative CNS has many parallels in stem cells and cancer. Thus, these axotomy-induced changes in DNA methylation in regenerative CNS provide evidence for a novel epigenetic state favoring successful over unsuccessful CNS axon regeneration. The datasets described in this study should help lay the foundations for future studies of the molecular and cellular mechanisms involved. The insights gained should, in turn, help point the way to novel therapeutic approaches for treating CNS injury in mammals.

DNA methylation is regulated by both the stability and topology of G-quadruplex

The methylation reaction was regulated by not only the stability of G4 but also the topology of G4.

Bisulfite-free and single-nucleotide resolution sequencing of DNA epigenetic modification of 5-hydroxymethylcytosine by engineered deaminase

The discovery of 5-hydroxymethylcytosine (5hmC) in mammalian genomes is a landmark in epigenomics study. Similar to 5-methylcytosine (5mC), 5hmC is viewed a critical epigenetic modification. Deciphering the functions of 5hmC...

Homologues of epigenetic pyrimidines: 5-alkyl-, 5-hydroxyalkyl and 5-acyluracil and -cytosine nucleotides: synthesis, enzymatic incorporation into DNA and effect on transcription with bacterial RNA polymerase

Homologues of natural epigenetic pyrimidine nucleosides and nucleotides were designed and synthesized. They included 5-ethyl-, 5-propyl-, 5-(1-hydroxyethyl)-, 5-(1-hydroxypropyl)- and 5-acetyl- and 5-propionylcytosine and -uracil 2′-deoxyribonucleosides and their corresponding 5′-O-triphosphates (dN^XTPs). The epimers of 5-(1-hydroxyethyl)- and 5-(1-hydroxypropyl)pyrimidine nucleosides were separated and their absolute configuration was determined by a combination of X-ray and NMR analysis. The modified dN^XTPs were used as substrates for PCR synthesis of modified DNA templates used for the study of transcription with bacterial RNA polymerase. Fundamental differences in transcription efficiency were observed, depending on the various modifications. The most notable effects included pronounced stimulation of transcription from 5-ethyluracil-bearing templates (200% transcription yield compared to natural thymine) and an enhancing effect of 5-acetylcytosine versus inhibiting effect of 5-acetyluracil. In summary, these results reveal that RNA polymerase copes with dramatically altered DNA structure and suggest that these nucleobases could potentially play roles as artificial epigenetic DNA nucleobases.

Nucleotides derived from homologues of epigenetic pyrimidine bases were prepared and used for polymerase synthesis of modified DNA templates. Interesting effects of the substituents on PCR and transcription have been observed.

Genome-wide mapping of N4-methylcytosine at single-base resolution by APOBEC3A-mediated deamination sequencing

N⁴-methylcytosine (4mC) is a natural DNA modification occurring in thermophiles and plays important roles in restriction-modification (R-M) systems in bacterial genomes. However, the precise location and sequence context of 4mC in the whole genome are limited. In this study, we developed an APOBEC3A-mediated deamination sequencing (4mC-AMD-seq) method for genome-wide mapping of 4mC at single-base resolution. In the 4mC-AMD-seq method, cytosine and 5-methylcytosine (5mC) are deaminated by APOBEC3A (A3A) protein to generate uracil and thymine, both of which are read as thymine in sequencing, while 4mC is resistant to deamination and therefore read as cytosine. Thus, the readouts of cytosines from sequencing could manifest the original 4mC sites in genomes. With the 4mC-AMD-seq method, we achieved the genome-wide mapping of 4mC in Deinococcus radiodurans (D. radiodurans). In addition, we confirmed that 4mC, but not 5mC, was the major modification in the D. radiodurans genome. We identified 1586 4mC sites in the genome of D. radiodurans, among which 564 sites were located in the CCGCGG motif. The average methylation levels in the CCGCGG motif and non-CCGCGG sequence were 70.0% and 22.8%, respectively. We envision that the 4mC-AMD-seq method will facilitate the investigation of 4mC functions, including the 4mC-involved R-M systems, in uncharacterized but potentially useful strains.

Genome-wide mapping of N⁴-methylcytosine (4mC) at single-base resolution with APOBEC3A-mediated deamination sequencing (4mC-AMD-seq).

Radiation-induced cardiovascular disease: an overlooked role for DNA methylation?

Radiotherapy in cancer treatment involves the use of ionizing radiation for cancer cell killing. Although radiotherapy has shown significant improvements on cancer recurrence and mortality, several radiation-induced adverse effects have been documented. Of these adverse effects, radiation-induced cardiovascular disease (CVD) is particularly prominent among patients receiving mediastinal radiotherapy, such as breast cancer and Hodgkin's lymphoma patients. A number of mechanisms of radiation-induced CVD pathogenesis have been proposed such as endothelial inflammatory activation, premature endothelial senescence, increased ROS and mitochondrial dysfunction. However, current research seems to point to a so-far unexamined and potentially novel involvement of epigenetics in radiation-induced CVD pathogenesis. Firstly, epigenetic mechanisms have been implicated in CVD pathophysiology. In addition, several studies have shown that ionizing radiation can cause epigenetic modifications, especially DNA methylation alterations. As a result, this review aims to provide a summary of the current literature linking DNA methylation to radiation-induced CVD and thereby explore DNA methylation as a possible contributor to radiation-induced CVD pathogenesis.

Dysmaturational Longitudinal Epigenetic Aging During Transition to Psychosis

Psychosis frequently occurs during adolescence and young adulthood, possibly as a result of gene-environment interactions, mediated by epigenetic mechanisms such as DNA methylation. Methylation patterns can be leveraged to predict epigenetic age in order to identify anomalies in aging trajectories that may be associated with the emergence of psychosis. Thus, epigenetic age may provide a measurable surrogate of psychotic risk or psychosis' emergence, and shed light on the neurodevelopmental model of psychosis. In this study, we present the first longitudinal analysis of epigenetic age trajectory during conversion to psychosis in a population at ultra-high-risk, with available genome-wide methylation DNA at two time points, at baseline and after one year of follow-up (N = 38 × 2). After predicting epigenetic age, we computed epigenetic age gap as the cross-sectional difference between real age and predicted age, and (longitudinal) epigenetic age acceleration as the derivative of predicted age with respect to time. At baseline, future converters were 2.7 years younger than nonconverters and this difference disappeared at follow-up, when some converted to psychosis. This is because during conversion to psychosis, the epigenetic age of converters accelerated by 2.8 years/year compared to nonconverters. This acceleration was robust with a strictly positive 95% confidence interval, and held its significance after adjustment for age, sex, and cannabis intake. The methylation sites most associated with aging were on genes also linked with schizophrenia and neurodevelopmental disorders. This accelerated age trajectory, following a previous deceleration, may therefore reflect dysmaturational processes.

Correlation between small-scale methylation changes and gene expression during the development of myopia

Visually induced changes in the expression of early growth response-1 (EGR1), FBJ osteosarcoma oncogene (FOS), and NGFI-A binding protein-2 (NAB2) appear to form a part of a retinal network fundamental to ocular growth regulation, and thus, the development of myopia (short-sightedness). However, it is unclear how environmental (visual) cues are translated into these molecular changes. One possibility is through epigenetic modifications such as DNA methylation, a known regulator of such processes. By sequencing bisulfite-converted DNA amplicons, this study examined whether changes in DNA methylation occur within specific regulatory and promoter regions of EGR1, FOS, and NAB2 during the periods of increased and decreased ocular growth in chicks. Visually induced changes in ocular growth rates were associated with single-point, but not large-scale, shifts in methylation levels within the investigated regions. Analysis of methylation pattern variability (entropy) demonstrated that the observed methylation changes are occurring within small subpopulations of retinal cells. This concurs with previous observations that EGR1 and FOS are differentially regulated at the peptide level within specific retinal cell types. Together, the findings of this study support a potential role for DNA methylation in the translation of external visual cues into molecular changes critical for ocular growth regulation and myopia development.

Epigenetic Studies of Chinese Herbal Medicine: Pleiotropic Role of DNA Methylation

Accumulating knowledge has been achieved on DNA methylation participating in numerous cellular processes and multiple human diseases; however, few studies have addressed the pleiotropic role of DNA methylation in Chinese herbal medicine (CHM). CHM has been used worldwide for the prevention and treatment of multiple diseases. Newly developed epigenetic techniques have brought great opportunities for the development of CHM. In this review, we summarize the DNA methylation studies and portray the pleiotropic role of DNA methylation in CHM. DNA methylation serves as a mediator participating in plant responses to environmental factors, and thus affecting CHM medicinal plants growth and bioactive compound biosynthesis which are vital for therapeutic effects. Furthermore, DNA methylation helps to uncover the pharmaceutical mechanisms of CHM formulae, herbs, and herbal-derived compounds. It also provides scientific validation for constitution theory and other essential issues of CHM. This newly developed field of DNA methylation is up-and-coming to address many complicated scientific questions of CHM; it thus not only promotes disease treatment but also facilitates health maintenance.

Epigenetic Age and the Risk of Incident Atrial Fibrillation

BACKGROUND

The most prominent risk factor for atrial fibrillation (AF) is chronological age; however, underlying mechanisms are unexplained. Algorithms using epigenetic modifications to the human genome effectively predict chronological age. Chronological and epigenetic predicted ages may diverge in a phenomenon referred to as epigenetic age acceleration (EAA), which may reflect accelerated biological aging. We sought to evaluate for associations between epigenetic age measures and incident AF.

METHODS

Measures for 4 epigenetic clocks (Horvath, Hannum, DNA methylation [DNAm] PhenoAge, and DNAm GrimAge) and an epigenetic predictor of PAI-1 (plasminogen activator inhibitor-1) levels (ie, DNAm PAI-1) were determined for study participants from 3 population-based cohort studies. Cox models evaluated for associations with incident AF and results were combined via random-effects meta-analyses. Two-sample summary-level Mendelian randomization analyses evaluated for associations between genetic instruments of the EAA measures and AF.

RESULTS

Among 5600 participants (mean age, 65.5 years; female, 60.1%; Black, 50.7%), there were 905 incident AF cases during a mean follow-up of 12.9 years. Unadjusted analyses revealed all 4 epigenetic clocks and the DNAm PAI-1 predictor were associated with statistically significant higher hazards of incident AF, though the magnitudes of their point estimates were smaller relative to the associations observed for chronological age. The pooled EAA estimates for each epigenetic measure, with the exception of Horvath EAA, were associated with incident AF in models adjusted for chronological age, race, sex, and smoking variables. After multivariable adjustment for additional known AF risk factors that could also potentially function as mediators, pooled EAA measures for 2 clocks remained statistically significant. Five-year increases in EAA measures for DNAm GrimAge and DNAm PhenoAge were associated with 19% (adjusted hazard ratio [HR], 1.19 [95% CI, 1.09-1.31]; P<0.01) and 15% (adjusted HR, 1.15 [95% CI, 1.05-1.25]; P<0.01) higher hazards of incident AF, respectively. Mendelian randomization analyses for the 5 EAA measures did not reveal statistically significant associations with AF.

CONCLUSIONS

Our study identified adjusted associations between EAA measures and incident AF, suggesting that biological aging plays an important role independent of chronological age, though a potential underlying causal relationship remains unclear. These aging processes may be modifiable and not constrained by the immutable factor of time.

Coordinated glucocorticoid receptor and MAFB action induces tolerogenesis and epigenome remodeling in dendritic cells

Glucocorticoids (GCs) exert potent anti-inflammatory effects in immune cells through the glucocorticoid receptor (GR). Dendritic cells (DCs), central actors for coordinating immune responses, acquire tolerogenic properties in response to GCs. Tolerogenic DCs (tolDCs) have emerged as a potential treatment for various inflammatory diseases. To date, the underlying cell type-specific regulatory mechanisms orchestrating GC-mediated acquisition of immunosuppressive properties remain poorly understood. In this study, we investigated the transcriptomic and epigenomic remodeling associated with differentiation to DCs in the presence of GCs. Our analysis demonstrates a major role of MAFB in this process, in synergy with GR. GR and MAFB both interact with methylcytosine dioxygenase TET2 and bind to genomic loci that undergo specific demethylation in tolDCs. We also show that the role of MAFB is more extensive, binding to thousands of genomic loci in tolDCs. Finally, MAFB knockdown erases the tolerogenic properties of tolDCs and reverts the specific DNA demethylation and gene upregulation. The preeminent role of MAFB is also demonstrated in vivo for myeloid cells from synovium in rheumatoid arthritis following GC treatment. Our results imply that, once directly activated by GR, MAFB plays a critical role in orchestrating the epigenomic and transcriptomic remodeling that define the tolerogenic phenotype.

Various Uses of PD1/PD-L1 Inhibitor in Oncology: Opportunities and Challenges

The occurrence and development of cancer are closely related to the immune escape of tumor cells and immune tolerance. Unlike previous surgical, chemotherapy, radiotherapy and targeted therapy, tumor immunotherapy is a therapeutic strategy that uses various means to stimulate and enhance the immune function of the body, and ultimately achieves the goal of controlling tumor cells.With the in-depth understanding of tumor immune escape mechanism and tumor microenvironment, and the in-depth study of tumor immunotherapy, immune checkpoint inhibitors represented by Programmed Death 1/Programmed cell Death-Ligand 1(PD-1/PD-L1) inhibitors are becoming increasingly significant in cancer medication treatment. employ a variety of ways to avoid detection by the immune system, a single strategy is not more effective in overcoming tumor immune evasion and metastasis. Combining different immune agents or other drugs can effectively address situations where immunotherapy is not efficacious, thereby increasing the chances of success and alternative access to alternative immunotherapy. Immune combination therapies for cancer have become a hot topic in cancer treatment today. In this paper, several combination therapeutic modalities of PD1/PD-L1 inhibitors are systematically reviewed. Finally, an analysis and outlook are provided in the context of the recent advances in combination therapy with PD1/PD-L1 inhibitors and the pressing issues in this field.

Epigenetic inheritance and reproductive mode in plants and animals

Epigenetic inheritance is another piece of the puzzle of nongenetic inheritance, although the prevalence, sources, persistence, and phenotypic consequences of heritable epigenetic marks across taxa remain unclear. We systematically reviewed over 500 studies from the past 5 years to identify trends in the frequency of epigenetic inheritance due to differences in reproductive mode and germline development. Genetic, intrinsic (e.g., disease), and extrinsic (e.g., environmental) factors were identified as sources of epigenetic inheritance, with impacts on phenotype and adaptation depending on environmental predictability. Our review shows that multigenerational persistence of epigenomic patterns is common in both plants and animals, but also highlights many knowledge gaps that remain to be filled. We provide a framework to guide future studies towards understanding the generational persistence and eco-evolutionary significance of epigenomic patterns.

Pterostilbene leads to DNMT3B-mediated DNA methylation and silencing of OCT1-targeted oncogenes in breast cancer cells

Transcription factor (TF)-mediated regulation of genes is often disrupted during carcinogenesis. The DNA methylation state of TF-binding sites may dictate transcriptional activity of corresponding genes. Stilbenoid polyphenols, such as pterostilbene (PTS), have been shown to exert anticancer action by remodeling DNA methylation and gene expression. However, the mechanisms behind these effects still remain unclear. Here, the dynamics between oncogenic TF OCT1 binding and de novo DNA methyltransferase DNMT3B binding in PTS-treated MCF10CA1a invasive breast cancer cells has been explored. Using chromatin immunoprecipitation (ChIP) followed by next generation sequencing, we determined 47 gene regulatory regions with decreased OCT1 binding and enriched DNMT3B binding in response to PTS. Most of those genes were found to have oncogenic functions. We selected three candidates, PRKCA, TNNT2, and DANT2, for further mechanistic investigation taking into account PRKCA functional and regulatory connection with numerous cancer-driving processes and pathways, and some of the highest increase in DNMT3B occupancy within TNNT2 and DANT2 enhancers. PTS led to DNMT3B recruitment within PRKCA, TNNT2, and DANT2 at loci that also displayed reduced OCT1 binding. Substantial decrease in OCT1 with increased DNMT3B binding was accompanied by PRKCA promoter and TNNT2 and DANT2 enhancer hypermethylation, and gene silencing. Interestingly, DNA hypermethylation of the genes was not detected in response to PTS in DNMT3B-CRISPR knockout MCF10CA1a breast cancer cells. It indicates DNMT3B-dependent methylation of PRKCA, TNNT2, and DANT2 upon PTS. Our findings provide a better understanding of mechanistic players and their gene targets that possibly contribute to the anticancer action of stilbenoid polyphenols.

Effect of DNA modifications on the transition between canonical and non-canonical DNA structures in CpG islands during senescence

Patterns and levels of DNA modifications play important roles in senescence. Two major epigenetic modifications of DNA, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), target CpG sites. Importantly, CpG concentrated regions, known as CpG islands, contain GC-rich sequences, which have the potential to fold into non-canonical DNA structures such as i-motifs and G-quadruplexes. In this study, we investigated the effect of 5mC and 5hmC modifications on the transition between a duplex, and i-motif and G-quadruplexes. To examine the transition, we firstly investigated the stability and structure of the i-motif and G-quadruplexes, considering the molecular environment in senescent cells. Analyses of their stability showed that the modifications did not drastically affect the stability. However, noteworthily, the modification can weaken the (de)stabilisation effect on G-quadruplexes caused by cosolute(s) and cations. Circular dichroism analyses indicated that the surrounding environments, including the molecular crowding and the type of cations such as K⁺ and Na⁺, regulate the topology of G-quadruplexes, while neither 5mC nor 5hmC had a drastic effect. On the other hand, the modifications changed the transition between duplexes and quadruplexes. Unmodified DNA preferred to fold into quadruplexes, whereas DNA with 5mC and 5hmC preferred to fold into duplexes in the absence of PEG200; on the other hand, DNA with or without modifications tended to fold into i-motifs under crowded conditions. Furthermore, an investigation of quadruplexes forming sequences in CpG islands, which are hyper- or hypomethylated during senescence, followed by gene ontology enrichment analysis for each gene group classified by the presence of quadruplexes, showed a difference in function between genes with and without quadruplexes in the CpG region. These results indicate that it is important to consider the effects of patterns and levels of DNA modifications on the transition between canonical and non-canonical DNA structures to understand gene regulation by epigenetic modification during senescence.

The modification of DNA can regulate the transition between a duplex and quadruplexes during senescence responding to surrounding environments.

Comparison of EM-seq and PBAT methylome library methods for low-input DNA

DNA methylation is the most studied epigenetic mark involved in regulation of gene expression. For low input samples, a limited number of methods for quantifying DNA methylation genome-wide has been evaluated. Here, we compared a series of input DNA amounts (1–10ng) from two methylome library preparation protocols, enzymatic methyl-seq (EM-seq) and post-bisulfite adaptor tagging (PBAT) adapted from single-cell PBAT. EM-seq takes advantage of enzymatic activity while PBAT relies on conventional bisulfite conversion for detection of DNA methylation. We found that both methods accurately quantified DNA methylation genome-wide. They produced expected distribution patterns around genomic features, high C-T transition efficiency at non-CpG sites and high correlation between input amounts. However, EM-seq performed better in regard to library and sequencing quality, i.e. EM-seq produced larger insert sizes, higher alignment rates and higher library complexity with lower duplication rate compared to PBAT. Moreover, EM-seq demonstrated higher CpG coverage, better CpG site overlap and higher consistency between input series. In summary, our data suggests that EM-seq overall performed better than PBAT in whole-genome methylation quantification of low input samples.

An enzyme-mediated bioorthogonal labeling method for genome-wide mapping of 5-hydroxymethyluracil

DNA 5-hydroxymethyluracil (5hmU) is a thymine modification existing in the genomes of various organisms. The post-replicative formation of 5hmU occurs via hydroxylation of thymine by ten-eleven translocation (TET) dioxygenases in mammals and J-binding proteins (JBPs) in protozoans, respectively. In addition, 5hmU can also be generated through oxidation of thymine by reactive oxygen species or deamination of 5hmC by cytidine deaminase. While the biological roles of 5hmU have not yet been fully explored, determining its genomic location will highly assist in elucidating its functions. Herein, we report a novel enzyme-mediated bioorthogonal labeling method for selective enrichment of 5hmU in genomes. 5hmU DNA kinase (5hmUDK) was utilized to selectively install an azide (N₃) group or alkynyl group into the hydroxyl moiety of 5hmU followed by incorporation of the biotin linker through click chemistry, which enabled the capture of 5hmU-containing DNA fragments via streptavidin pull-down. The enriched fragments were applied to deep sequencing to determine the genomic distribution of 5hmU. With this established enzyme-mediated bioorthogonal labeling strategy, we achieved the genome-wide mapping of 5hmU in Trypanosoma brucei. The method described here will allow for a better understanding of the functional roles and dynamics of 5hmU in genomes.

Epigenetics of alcohol use disorder-A review of recent advances in DNA methylation profiling

Alcohol use disorder (AUD) is a major contributor to morbidity and mortality worldwide. Although there is a heritable component, the etiology of AUD is complex and can involve environmental exposures like trauma and can be associated with many different patterns of alcohol consumption. Epigenetic modifications, which can mediate the influence of genetic variants and environmental variables on gene expression, have emerged as an important area of AUD research. Over the past decade, the number of studies investigating AUD and DNA methylation, a form of epigenetic modification, has grown rapidly. Yet we are still far from understanding how DNA methylation contributes to or reflects aspects of AUD. In this paper, we reviewed studies of DNA methylation and AUD and discussed how the field has evolved. We found that global DNA and candidate DNA methylation studies did not produce replicable results. To assess whether findings of epigenome-wide association studies (EWAS) were replicated, we aggregated significant findings across studies and identified 184 genes and 15 gene ontological pathways that were differentially methylated in at least two studies and four genes and three gene ontological pathways that were differentially methylated in three studies. These genes and pathways repeatedly found enrichment of immune processes, which is in line with recent developments suggesting that the immune system may be altered in AUD. Finally, we assess the current limitations of studies of DNA methylation and AUD and make recommendations on how to design future studies to resolve outstanding questions.

singlecellVR: Interactive Visualization of Single-Cell Data in Virtual Reality

Single-cell assays have transformed our ability to model heterogeneity within cell populations. As these assays have advanced in their ability to measure various aspects of molecular processes in cells, computational methods to analyze and meaningfully visualize such data have required matched innovation. Independently, Virtual Reality (VR) has recently emerged as a powerful technology to dynamically explore complex data and shows promise for adaptation to challenges in single-cell data visualization. However, adopting VR for single-cell data visualization has thus far been hindered by expensive prerequisite hardware or advanced data preprocessing skills. To address current shortcomings, we present singlecellVR, a user-friendly web application for visualizing single-cell data, designed for cheap and easily available virtual reality hardware (e.g., Google Cardboard, ∼$8). singlecellVR can visualize data from a variety of sequencing-based technologies including transcriptomic, epigenomic, and proteomic data as well as combinations thereof. Analysis modalities supported include approaches to clustering as well as trajectory inference and visualization of dynamical changes discovered through modelling RNA velocity. We provide a companion software package, scvr to streamline data conversion from the most widely-adopted single-cell analysis tools as well as a growing database of pre-analyzed datasets to which users can contribute.

Recent technical advances in the study of nucleic acid modifications

Enzyme-mediated chemical modifications of nucleic acids are indispensable regulators of gene expression. Our understanding of the biochemistry and biological significance of these modifications has largely been driven by an ever-evolving landscape of technologies that enable accurate detection, mapping, and manipulation of these marks. Here we provide a summary of recent technical advances in the study of nucleic acid modifications with a focus on techniques that allow accurate detection and mapping of these modifications. For each modification discussed (N6-methyladenosine, 5-methylcytidine, inosine, pseudouridine, and N4-acetylcytidine), we begin by introducing the "gold standard" technique for its mapping and detection, followed by a discussion of techniques developed to address any shortcomings of the gold standard. By highlighting the commonalities and differences of these techniques, we hope to provide a perspective on the current state of the field and to lay out a guideline for development of future technologies.

Insights into the Role of Oxidative Stress in Ovarian Cancer

Oxidative stress (OS) arises when the body is subjected to harmful endogenous or exogenous factors that overwhelm the antioxidant system. There is increasing evidence that OS is involved in a number of diseases, including ovarian cancer (OC). OC is the most lethal gynecological malignancy, and risk factors include genetic factors, age, infertility, nulliparity, microbial infections, obesity, smoking, etc. OS can promote the proliferation, metastasis, and therapy resistance of OC, while high levels of OS have cytotoxic effects and induce apoptosis in OC cells. This review focuses on the relationship between OS and the development of OC from four aspects: genetic alterations, signaling pathways, transcription factors, and the tumor microenvironment. Furthermore, strategies to target aberrant OS in OC are summarized and discussed, with a view to providing new ideas for clinical treatment.

Intragenomic Decarboxylation of 5-Carboxy-2'-deoxycytidine

Cellular DNA is composed of four canonical nucleosides (dA, dC, dG and T), which form two Watson-Crick base pairs. In addition, 5-methylcytosine (mdC) may be present. The methylation of dC to mdC is known to regulate transcriptional activity. Next to these five nucleosides, the genome, particularly of stem cells, contains three additional dC derivatives, which are formed by stepwise oxidation of the methyl group of mdC with the help of Tet enzymes. These are 5-hydroxymethyl-dC (hmdC), 5-formyl-dC (fdC), and 5-carboxy-dC (cadC). It is believed that fdC and cadC are converted back into dC, which establishes an epigenetic control cycle that starts with methylation of dC to mdC, followed by oxidation and removal of fdC and cadC. While fdC was shown to undergo intragenomic deformylation to give dC directly, a similar decarboxylation of cadC was postulated but not yet observed on the genomic level. By using metabolic labelling, we show here that cadC decarboxylates in several cell types, which confirms that both fdC and cadC are nucleosides that are directly converted back to dC within the genome by C-C bond cleavage.

Genes, environment, and developmental timing: New insights from translational approaches to understand early origins of respiratory diseases

Over the past decade, "omics" approaches have advanced our understanding of the molecular programming of the airways in humans. Several studies have identified potential molecular mechanisms that contribute to early life epigenetic reprogramming, including DNA methylation, histone modifications, microRNAs, and the homeostasis of the respiratory mucosa (epithelial function and microbiota). Current evidence supports the notion that early infancy is characterized by heightened susceptibility to airway genetic reprogramming in response to the first exposures in life, some of which can have life-long consequences. Here, we summarize and analyze the latest insights from studies that support a novel epigenetic paradigm centered on human maturational and developmental programs including three cardinal elements: genes, environment, and developmental timing. The combination of these factors is likely responsible for the functional trajectory of the respiratory system at the molecular, functional, and clinical levels.

DNA methylation atlas of the mouse brain at single-cell resolution

Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing1,2 to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data3 enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments4. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.

DNA methylation dynamics of long noncoding RNA during human fetal development

Aim: To determine whether the promoters of long noncoding RNAs (lncRNAs) undergo dynamic changes in DNA methylation during fetal development. Methods: ANOVA and the tissue specificity index were used to identify and validate tissue-specific methylation sites. Age-associated DNA methylation signatures were identified by applying the elastic net method. Results: The lncRNA methylome landscape was characterized in four types of fetal tissue and at three gestational time points, and specific characteristics relative to the tissue of origin and developmental age were identified. Higher levels of lncRNA methylation might be involved in tissue differentiation. LncRNAs harboring age-associated methylation signatures may participate in the fetal developmental process. Conclusion: This study provides novel insights into the role of lncRNA methylomes in fetal tissue specification and development.

Cell-free DNA TAPS provides multimodal information for early cancer detection

Multimodal, genome-wide characterization of epigenetic and genetic information in circulating cell-free DNA (cfDNA) could enable more sensitive early cancer detection, but it is technologically challenging. Recently, we developed TET-assisted pyridine borane sequencing (TAPS), which is a mild, bisulfite-free method for base-resolution direct DNA methylation sequencing. Here, we optimized TAPS for cfDNA (cfTAPS) to provide high-quality and high-depth whole-genome cell-free methylomes. We applied cfTAPS to 85 cfDNA samples from patients with hepatocellular carcinoma (HCC) or pancreatic ductal adenocarcinoma (PDAC) and noncancer controls. From only 10 ng of cfDNA (1 to 3 ml of plasma), we generated the most comprehensive cfDNA methylome to date. We demonstrated that cfTAPS provides multimodal information about cfDNA characteristics, including DNA methylation, tissue of origin, and DNA fragmentation. Integrated analysis of these epigenetic and genetic features enables accurate identification of early HCC and PDAC.

Extended-representation bisulfite sequencing of gene regulatory elements in multiplexed samples and single cells

The biological roles of DNA methylation have been elucidated by profiling methods based on whole-genome or reduced-representation bisulfite sequencing, but these approaches do not efficiently survey the vast numbers of non-coding regulatory elements in mammalian genomes. Here we present an extended-representation bisulfite sequencing (XRBS) method for targeted profiling of DNA methylation. Our design strikes a balance between expanding coverage of regulatory elements and reproducibly enriching informative CpG dinucleotides in promoters, enhancers and CTCF binding sites. Barcoded DNA fragments are pooled before bisulfite conversion, allowing multiplex processing and technical consistency in low-input samples. Application of XRBS to single leukemia cells enabled us to evaluate genetic copy number variations and methylation variability across individual cells. Our analysis highlights heterochromatic H3K9me3 regions as having the highest cell-to-cell variability in their methylation, likely reflecting inherent epigenetic instability of these late-replicating regions, compounded by differences in cell cycle stages among sampled cells.