Age Prediction Using DNA Methylation Heterogeneity Metrics.

Dynamic changes in genomic DNA methylation patterns govern the epigenetic developmental programs and accompany the organism's aging. Epigenetic clock (eAge) algorithms utilize DNA methylation to estimate the age and risk factors for diseases as well as analyze the impact of various interventions. High-throughput bisulfite sequencing methods, such as reduced-representation bisulfite sequencing (RRBS) or whole genome bisulfite sequencing (WGBS), provide an opportunity to identify the genomic regions of disordered or heterogeneous DNA methylation, which might be associated with cell-type heterogeneity, DNA methylation erosion, and allele-specific methylation. We systematically evaluated the applicability of five scores assessing the variability of methylation patterns by evaluating within-sample heterogeneity (WSH) to construct human blood epigenetic clock models using RRBS data. The best performance was demonstrated by the model based on a metric designed to assess DNA methylation erosion with an MAE of 3.686 years. We also trained a prediction model that uses the average methylation level over genomic regions. Although this region-based model was relatively more efficient than the WSH-based model, the latter required the analysis of just a few short genomic regions and, therefore, could be a useful tool to design a reduced epigenetic clock that is analyzed by targeted next-generation sequencing.

Kernel-based testing for single-cell differential analysis.

Single-cell technologies offer insights into molecular feature distributions, but comparing them poses challenges. We propose a kernel-testing framework for non-linear cell-wise distribution comparison, analyzing gene expression and epigenomic modifications. Our method allows feature-wise and global transcriptome/epigenome comparisons, revealing cell population heterogeneities. Using a classifier based on embedding variability, we identify transitions in cell states, overcoming limitations of traditional single-cell analysis. Applied to single-cell ChIP-Seq data, our approach identifies untreated breast cancer cells with an epigenomic profile resembling persister cells. This demonstrates the effectiveness of kernel testing in uncovering subtle population variations that might be missed by other methods.

Learning Consistency and Specificity of Cells From Single-Cell Multi-Omic Data.

Advancements in single-cell technologies concomitantly develop the epigenomic and transcriptomic profiles at the cell levels, providing opportunities to explore the potential biological mechanisms. Even though significant efforts have been dedicated to them, it remains challenging for the integration analysis of multi-omic data of single-cell because of the heterogeneity, complicated coupling and interpretability of data. To handle these issues, we propose a novel self-representation Learning-based Multi-omics data Integrative Clustering algorithm (sLMIC) for the integration of single-cell epigenomic profiles (DNA methylation or scATAC-seq) and transcriptomic (scRNA-seq), which the consistent and specific features of cells are explicitly extracted facilitating the cell clustering. Specifically, sLMIC constructs a graph for each type of single-cell data, thereby transforming omics data into multi-layer networks, which effectively removes heterogeneity of omic data. Then, sLMIC employs the low-rank and exclusivity constraints to separate the self-representation of cells into two parts, i.e., the shared and specific features, which explicitly characterize the consistency and diversity of omic data, providing an effective strategy to model the structure of cell types. Feature extraction and cell clustering are jointly formulated as an overall objective function, where latent features of data are obtained under the guidance of cell clustering. The extensive experimental results on 13 multi-omics datasets of single-cell from diverse organisms and tissues indicate that sLMIC observably exceeds the advanced algorithms regarding various measurements.

methscore: a comprehensive R function for DNA methylation-based health predictors.

MOTIVATION: DNA methylation-based predictors of various biological metrics have been widely published and are becoming valuable tools in epidemiologic studies of epigenetics and personalized medicine. However, generating these predictors from original source software and web servers is complex and time consuming. Furthermore, different predictors were often derived based on data from different types of arrays, where array differences and batch effects can make predictors difficult to compare across studies. RESULTS: We integrate these published methods into a single R function to produce 158 previously published predictors for chronological age, biological age, exposures, lifestyle traits and serum protein levels using both classical and principal component-based methods. To mitigate batch and array differences, we also provide a modified RCP method (ref-RCP) that normalize input DNA methylation data to reference data prior to estimation. Evaluations in real datasets show that this approach improves estimate precision and comparability across studies. AVAILABILITY AND IMPLEMENTATION: The function was included in software package ENmix, and is freely available from Bioconductor website (https://www.bioconductor.org/packages/release/bioc/html/ENmix.html).

scMNMF: a novel method for single-cell multi-omics clustering based on matrix factorization.

MOTIVATION: The technology for analyzing single-cell multi-omics data has advanced rapidly and has provided comprehensive and accurate cellular information by exploring cell heterogeneity in genomics, transcriptomics, epigenomics, metabolomics and proteomics data. However, because of the high-dimensional and sparse characteristics of single-cell multi-omics data, as well as the limitations of various analysis algorithms, the clustering performance is generally poor. Matrix factorization is an unsupervised, dimensionality reduction-based method that can cluster individuals and discover related omics variables from different blocks. Here, we present a novel algorithm that performs joint dimensionality reduction learning and cell clustering analysis on single-cell multi-omics data using non-negative matrix factorization that we named scMNMF. We formulate the objective function of joint learning as a constrained optimization problem and derive the corresponding iterative formulas through alternating iterative algorithms. The major advantage of the scMNMF algorithm remains its capability to explore hidden related features among omics data. Additionally, the feature selection for dimensionality reduction and cell clustering mutually influence each other iteratively, leading to a more effective discovery of cell types. We validated the performance of the scMNMF algorithm using two simulated and five real datasets. The results show that scMNMF outperformed seven other state-of-the-art algorithms in various measurements. AVAILABILITY AND IMPLEMENTATION: scMNMF code can be found at https://github.com/yushanqiu/scMNMF.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality.

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Disease related changes in ATAC-seq of iPSC-derived motor neuron lines from ALS patients and controls.

Amyotrophic Lateral Sclerosis (ALS), like many other neurodegenerative diseases, is highly heritable, but with only a small fraction of cases explained by monogenic disease alleles. To better understand sporadic ALS, we report epigenomic profiles, as measured by ATAC-seq, of motor neuron cultures derived from a diverse group of 380 ALS patients and 80 healthy controls. We find that chromatin accessibility is heavily influenced by sex, the iPSC cell type of origin, ancestry, and the inherent variance arising from sequencing. Once these covariates are corrected for, we are able to identify ALS-specific signals in the data. Additionally, we find that the ATAC-seq data is able to predict ALS disease progression rates with similar accuracy to methods based on biomarkers and clinical status. These results suggest that iPSC-derived motor neurons recapitulate important disease-relevant epigenomic changes.

Epigenomic mechanism regulating the quality and ripeness of apple fruit with differing harvest maturity.

Harvest maturity significantly affects the quality of apple fruit in post-harvest storage process. Although the regulatory mechanisms underlying fruit ripening have been studied, the associated epigenetic modifications remain unclear. Thus, we compared the DNA methylation changes and the transcriptional responses of mature fruit (MF) and immature fruit (NF). There were significant correlations between DNA methylation and gene expression. Moreover, the sugar contents (sucrose, glucose, and fructose) were higher in MF than in NF, whereas the opposite pattern was detected for the starch content. The expression-level differences were due to DNA methylations and ultimately resulted in diverse fruit textures and ripeness. Furthermore, the higher ethylene, auxin, and abscisic acid levels in MF than in NF, which influenced the fruit texture and ripening, were associated with multiple differentially expressed genes in hormone synthesis, signaling, and response pathways (ACS, ACO, ZEP, NCED, and ABA2) that were regulated by DNA methylations. Multiple transcription factor genes involved in regulating fruit ripening and quality via changes in DNA methylation were identified, including MIKCC-type MADS-box genes and fruit ripening-related genes (NAP, SPL, WRKY, and NAC genes). These findings reflect the diversity in the epigenetic regulation of gene expression and may be relevant for elucidating the epigenetic regulatory mechanism underlying the ripening and quality of apple fruit with differing harvest maturity.

DNA Methylation in Ctenophores.

Epigenomic regulation and dynamic DNA methylation, in particular, are widespread mechanisms orchestrating the genome operation across time and species. Whole-genome bisulfite sequencing (WGBS) is currently the only method for unbiasedly capturing the presence of 5-methylcytosine (5-mC) DNA methylation patterns across an entire genome with single-nucleotide resolution. Bisulfite treatment converts unmethylated cytosines to uracils but leaves methylated cytosines intact, thereby creating a map of all methylated cytosines across a genome also known as a methylome. These epigenomic patterns of DNA methylation have been found to regulate gene expression and influence gene evolution rates between species. While protocols have been optimized for vertebrate methylome production, little adaptation has been done for invertebrates. Creating a methylome reference allows comparisons to be made between rates of transcription and epigenomic patterning in animals. Here we present a method of library construction for bisulfite sequencing optimized for non-bilateral metazoans such as the ctenophore, Mnemiopsis leidyi. We have improved upon our previously published method by including spike-in genomic DNA controls to measure methylation conversion rates. By pooling two bisulfite conversion reactions from the same individual, we also produced sequencing libraries that yielded a higher percentage of sequenced reads uniquely mapping to the reference genome. We successfully detected 5-mC in whole-animal methylomes at CpG, CHG, and CHH sites and visualized datasets using circos diagrams. The proof-of-concept tests were performed both under control conditions and following injury tests with changes in methylation patterns of genes encoding innexins, toxins and neuropeptides. Our approach can be easily adapted to produce epigenomes from other fragile marine animals.

Moving toward the Inclusion of Epigenomics in Bacterial Genome Evolution: Perspectives and Challenges.

The universality of DNA methylation as an epigenetic regulatory mechanism belongs to all biological kingdoms. However, while eukaryotic systems have been the primary focus of DNA methylation studies, the molecular mechanisms in prokaryotes are less known. Nevertheless, DNA methylation in prokaryotes plays a pivotal role in many cellular processes such as defense systems against exogenous DNA, cell cycle dynamics, and gene expression, including virulence. Thanks to single-molecule DNA sequencing technologies, genome-wide identification of methylated DNA is becoming feasible on a large scale, providing the possibility to investigate more deeply the presence, variability, and roles of DNA methylation. Here, we present an overview of the multifaceted roles of DNA methylation in prokaryotes and suggest research directions and tools which can enable us to better understand the contribution of DNA methylation to prokaryotic genome evolution and adaptation. In particular, we emphasize the need to understand the presence and role of transgenerational inheritance, as well as the impact of epigenomic signatures on adaptation and genome evolution. Research directions and the importance of novel computational tools are underlined.

Building Minimized Epigenetic Clock by iPlex MassARRAY Platform.

Epigenetic clocks are valuable tools for estimating both chronological and biological age by assessing DNA methylation levels at specific CpG dinucleotides. While conventional epigenetic clocks rely on genome-wide methylation data, targeted approaches offer a more efficient alternative. In this study, we explored the feasibility of constructing a minimized epigenetic clock utilizing data acquired through the iPlex MassARRAY technology. The study enrolled a cohort of relatively healthy individuals, and their methylation levels of eight specific CpG dinucleotides in genes SLC12A5, LDB2, FIGN, ACSS3, FHL2, and EPHX3 were evaluated using the iPlex MassARRAY system and the Illumina EPIC array. The methylation level of five studied CpG sites demonstrated significant correlations with chronological age and an acceptable convergence of data obtained by the iPlex MassARRAY and Illumina EPIC array. At the same time, the methylation level of three CpG sites showed a weak relationship with age and exhibited a low concordance between the data obtained from the two technologies. The construction of the epigenetic clock involved the utilization of different machine-learning models, including linear models, deep neural networks (DNN), and gradient-boosted decision trees (GBDT). The results obtained from these models were compared with each other and with the outcomes generated by other well-established epigenetic clocks. In our study, the TabNet architecture (deep tabular data learning architecture) exhibited the best performance (best MAE = 5.99). Although our minimized epigenetic clock yielded slightly higher age prediction errors compared to other epigenetic clocks, it still represents a viable alternative to the genome-wide epigenotyping array.

MammalMethylClock R package: software for DNA methylation-based epigenetic clocks in mammals.

MOTIVATION: Epigenetic clocks are prediction methods based on DNA methylation levels in a given species or set of species. Defined as multivariate regression models, these DNA methylation-based biomarkers of age or mortality risk are useful in species conservation efforts and in preclinical studies. RESULTS: We present an R package called MammalMethylClock for the construction, assessment, and application of epigenetic clocks in different mammalian species. The R package includes the utility for implementing pre-existing mammalian clocks from the Mammalian Methylation Consortium. AVAILABILITY AND IMPLEMENTATION: The source code and documentation manual for MammalMethylClock, and clock coefficient .csv files that are included within this software package, can be found on Zenodo at https://doi.org/10.5281/zenodo.10971037.

Adding the epigenomic signature to the prognostic jigsaw of myelodysplastic neoplasm?

Defining mechanisms of resistance to hypomethylating agents (HMAs) and biomarkers predictive of treatment response remains challenging in myelodysplastic neoplasm (MDS). Currently available prognostic tools that predict overall survival and transformation to acute myeloid leukaemia have not been powered to predict responses to HMAs. Noguera-Castells et al. comprehensively characterized the epigenomic profile in patients with MDS treated with azacitidine and described a methylation signature-based prognostic tool in predicting responses to azacitidine. Commentary on: Noguera-Castells et al. DNA methylation profiling of myelodysplastic syndromes and clinical response to azacitidine: a multicentre retrospective study. Br J Haematol 2024;204:1838-1843.

Impact of ambient air pollution on colorectal cancer risk and survival: insights from a prospective cohort and epigenetic Mendelian randomization study.

BACKGROUND: This study investigates the associations between air pollution and colorectal cancer (CRC) risk and survival from an epigenomic perspective. METHODS: Using a newly developed Air Pollutants Exposure Score (APES), we utilized a prospective cohort study (UK Biobank) to investigate the associations of individual and combined air pollution exposures with CRC incidence and survival, followed by an up-to-date systematic review with meta-analysis to verify the associations. In epigenetic two-sample Mendelian randomization analyses, we examine the associations between genetically predicted DNA methylation related to air pollution and CRC risk. Further genetic colocalization and gene-environment interaction analyses provided different insights to disentangle pathogenic effects of air pollution via epigenetic modification. FINDINGS: During a median 12.97-year follow-up, 5767 incident CRC cases among 428,632 participants free of baseline CRC and 533 deaths in 2401 patients with CRC were documented in the UK Biobank. A higher APES score was associated with an increased CRC risk (HR, 1.03, 95% CI = 1.01-1.06; P = 0.016) and poorer survival (HR, 1.13, 95% CI = 1.03-1.23; P = 0.010), particularly among participants with insufficient physical activity and ever smokers (Pinteraction > 0.05). A subsequent meta-analysis of seven observational studies, including UK Biobank data, corroborated the association between PM2.5 exposure (per 10 µg/m3 increment) and elevated CRC risk (RR,1.42, 95% CI = 1.12-1.79; P = 0.004; I2 = 90.8%). Genetically predicted methylation at PM2.5-related CpG site cg13835894 near TMBIM1/PNKD and cg16235962 near CXCR5, and NO2-related cg16947394 near TMEM110 were associated with an increased CRC risk. Gene-environment interaction analysis confirmed the epigenetic modification of aforementioned CpG sites with CRC risk and survival. INTERPRETATION: Our study suggests the association between air pollution and CRC incidence and survival, underscoring the possible modifying roles of epigenomic factors. Methylation may partly mediate pathogenic effects of air pollution on CRC, with annotation to epigenetic alterations in protein-coding genes TMBIM1/PNKD, CXCR5 and TMEM110. FUNDING: Xue Li is supported by the Natural Science Fund for Distinguished Young Scholars of Zhejiang Province (LR22H260001), the National Nature Science Foundation of China (No. 82204019) and Healthy Zhejiang One Million People Cohort (K-20230085). ET is supported by a Cancer Research UK Career Development Fellowship (C31250/A22804). MGD is supported by the MRC Human Genetics Unit Centre Grant (U127527198).

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors.

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Dissecting the Immune System through Gene Regulation.

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

Genomic and Epigenomic Changes in the Progeny of Cold-Stressed Arabidopsis thaliana Plants.

Plants are continuously exposed to various environmental stresses. Because they can not escape stress, they have to develop mechanisms of remembering stress exposures somatically and passing it to the progeny. We studied the Arabidopsis thaliana ecotype Columbia plants exposed to cold stress for 25 continuous generations. Our study revealed that multigenerational exposure to cold stress resulted in the changes in the genome and epigenome (DNA methylation) across generations. Main changes in the progeny were due to the high frequency of genetic mutations rather than epigenetic changes; the difference was primarily in single nucleotide substitutions and deletions. The progeny of cold-stressed plants exhibited the higher rate of missense non-synonymous mutations as compared to the progeny of control plants. At the same time, epigenetic changes were more common in the CHG (C = cytosine, H = cytosine, adenine or thymine, G = guanine) and CHH contexts and favored hypomethylation. There was an increase in the frequency of C to T (thymine) transitions at the CHH positions in the progeny of cold stressed plants; because this type of mutations is often due to the deamination of the methylated cytosines, it can be hypothesized that environment-induced changes in methylation contribute to mutagenesis and may be to microevolution processes and that RNA-dependent DNA methylation plays a crucial role. Our work supports the existence of heritable stress response in plants and demonstrates that genetic changes prevail.

Epigenomic tomography for probing spatially defined chromatin state in the brain.

Spatially resolved epigenomic profiling is critical for understanding biology in the mammalian brain. Single-cell spatial epigenomic assays were developed recently for this purpose, but they remain costly and labor intensive for examining brain tissues across substantial dimensions and surveying a collection of brain samples. Here, we demonstrate an approach, epigenomic tomography, that maps spatial epigenomes of mouse brain at the scale of centimeters. We individually profiled neuronal and glial fractions of mouse neocortex slices with 0.5 mm thickness. Tri-methylation of histone 3 at lysine 27 (H3K27me3) or acetylation of histone 3 at lysine 27 (H3K27ac) features across these slices were grouped into clusters based on their spatial variation patterns to form epigenomic brain maps. As a proof of principle, our approach reveals striking dynamics in the frontal cortex due to kainic-acid-induced seizure, linked with transmembrane ion transporters, exocytosis of synaptic vesicles, and secretion of neurotransmitters. Epigenomic tomography provides a powerful and cost-effective tool for characterizing brain disorders based on the spatial epigenome.

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection.

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

EpiVar Browser: advanced exploration of epigenomics data under controlled access.

MOTIVATION: Human epigenomic data has been generated by large consortia for thousands of cell types to be used as a reference map of normal and disease chromatin states. Since epigenetic data contains potentially identifiable information, similarly to genetic data, most raw files generated by these consortia are stored in controlled-access databases. It is important to protect identifiable information, but this should not hinder secure sharing of these valuable datasets. RESULTS: Guided by the Framework for responsible sharing of genomic and health-related data from the Global Alliance for Genomics and Health (GA4GH), we have developed an approach and a tool to facilitate the exploration of epigenomics datasets' aggregate results, while filtering out identifiable information. Specifically, the EpiVar Browser allows a user to navigate an epigenetic dataset from a cohort of individuals and enables direct exploration of genotype-chromatin phenotype relationships. Because individual genotypes and epigenetic signal tracks are not directly accessible, and rather aggregated in the portal output, no identifiable data is released, yet the interface allows for dynamic genotype-epigenome interrogation. This approach has the potential to accelerate analyses that would otherwise require a lengthy multi-step approval process and provides a generalizable strategy to facilitate responsible access to sensitive epigenomics data. AVAILABILITY AND IMPLEMENTATION: Online portal: https://computationalgenomics.ca/tools/epivar; EpiVar Browser source code: https://github.com/c3g/epivar-browser; bw-merge-window tool source code: https://github.com/c3g/bw-merge-window.