Tracing human stem cell lineage during development using DNA methylation

Trauma context exerts intergenerational effects on child mental health via DNA methylation

Many people experience traumatic or negative events, but few develop mental health issues as a result. This study investigated whether newborn DNA methylation (DNAm) previously associated with maternal childhood physical abuse by her father affected the child's mental health and physical growth, as well as whether it mediated or moderated developmental outcomes.

METHODS

Study sample (N = 903) and data came from Bristol University's Avon Longitudinal Study of Parents and Children. DNAm was measured in cord blood at birth. DNAm data was preprocessed, normalized, and quality controlled before subsetting to 60 CpG sites of interest from previous research. Linear regression analysis examined newborn DNAm and child development outcome associations. Sobel test examined the mediating relationship between mother's history of childhood abuse by father, newborn targeted gene DNAm of significant CpG sites, and child's mental health and physical growth. Moderation analyses examined the interaction effects between the significant CpG sites and mothers' physical abuse by their fathers on child's mental health and physical growth.

RESULTS

Full cohort analyses showed that newborn DNAm of several different CpG sites associates with separation anxiety, fear, and unhappy/tearful presentations in children aged 6-7 y. Sex-specific associations emerged with boys showing associations with anxiety and fear, and girls showing associations with fear and unhappiness. In boys only, cord blood DNAm mediates the effect of maternal childhood trauma on offspring mental health. No moderation effects emerged.

CONCLUSION

Intergenerational effects of mother's relationship to her abuser present in newborn DNAm associate with 7-year-old child's mental health, show sex-specific effects, and newborn DNAm does mediate maternal childhood trauma effects on offspring mental health in early-life.

Genetic and epigenetic alterations in aging and rejuvenation of human

All the information essential for life is encoded within our genome and epigenome, which orchestrates diverse cellular states spatially and temporally. In particular, the epigenome interacts with internal and external stimuli, encoding and preserving cellular experiences, and it serves as the regulatory base of the transcriptome across diverse cell types. The emergence of single-cell transcriptomic and epigenomic data collection has revealed unique omics signatures in diverse tissues, highlighting cellular heterogeneity. Recent research has documented age-related epigenetic changes at the single-cell level, alongside the validation of cellular rejuvenation through partial reprogramming, which involves simultaneous epigenetic modifications. These dynamic shifts, primarily fueled by stem cell plasticity, have catalyzed significant interest and cross-disciplinary research endeavors. This review explores the genomic and epigenomic alterations with aging, elucidating their reciprocal interactions. Additionally, it seeks to discuss the evolving landscape of rejuvenation research, with a particular emphasis on dissecting stem cell behavior through the lens of single-cell analysis. Moreover, it proposes potential research methodologies for future studies.

Glioma immune microenvironment composition calculator (GIMiCC): a method of estimating the proportions of eighteen cell types from DNA methylation microarray data

A scalable platform for cell typing in the glioma microenvironment can improve tumor subtyping and immune landscape detection as successful immunotherapy strategies continue to be sought and evaluated. DNA methylation (DNAm) biomarkers for molecular classification of tumor subtypes have been developed for clinical use. However, tools that predict the cellular landscape of the tumor are not well-defined or readily available. We developed the Glioma Immune Microenvironment Composition Calculator (GIMiCC), an approach for deconvolution of cell types in gliomas using DNAm data. Using data from 17 isolated cell types, we describe the derivation of the deconvolution libraries in the biological context of selected genomic regions and validate deconvolution results using independent datasets. We utilize GIMiCC to illustrate that DNAm-based estimates of immune composition are clinically relevant and scalable for potential clinical implementation. In addition, we utilize GIMiCC to identify composition-independent DNAm alterations that are associated with high immune infiltration. Our future work aims to optimize GIMiCC and advance the clinical evaluation of glioma.

Deciphering the role of immune cell composition in epigenetic age acceleration: Insights from cell-type deconvolution applied to human blood epigenetic clocks

Aging is a significant risk factor for various human disorders, and DNA methylation clocks have emerged as powerful tools for estimating biological age and predicting health-related outcomes. Methylation data from blood DNA has been a focus of more recently developed DNA methylation clocks. However, the impact of immune cell composition on epigenetic age acceleration (EAA) remains unclear as only some clocks incorporate partial cell type composition information when analyzing EAA. We investigated associations of 12 immune cell types measured by cell-type deconvolution with EAA predicted by six widely-used DNA methylation clocks in data from >10,000 blood samples. We observed significant associations of immune cell composition with EAA for all six clocks tested. Across the clocks, nine or more of the 12 cell types tested exhibited significant associations with EAA. Higher memory lymphocyte subtype proportions were associated with increased EAA, and naïve lymphocyte subtypes were associated with decreased EAA. To demonstrate the potential confounding of EAA by immune cell composition, we applied EAA in rheumatoid arthritis. Our research maps immune cell type contributions to EAA in human blood and offers opportunities to adjust for immune cell composition in EAA studies to a significantly more granular level. Understanding associations of EAA with immune profiles has implications for the interpretation of epigenetic age and its relevance in aging and disease research. Our detailed map of immune cell type contributions serves as a resource for studies utilizing epigenetic clocks across diverse research fields, including aging-related diseases, precision medicine, and therapeutic interventions.

Matched analysis of detailed peripheral blood and tumor immune microenvironment profiles in bladder cancer

Background: Bladder cancer and therapy responses hinge on immune profiles in the tumor microenvironment (TME) and blood, yet studies linking tumor-infiltrating immune cells to peripheral immune profiles are limited. Methods: DNA methylation cytometry quantified TME and matched peripheral blood immune cell proportions. With tumor immune profile data as the input, subjects were grouped by immune infiltration status and consensus clustering. Results: Immune hot and cold groups had different immune compositions in the TME but not in circulating blood. Two clusters of patients identified with consensus clustering had different immune compositions not only in the TME but also in blood. Conclusion: Detailed immune profiling via methylation cytometry reveals the significance of understanding tumor and systemic immune relationships in cancer patients.

Single-cell DNA methylation and 3D genome architecture in the human brain

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Wide field-of-view fluorescence imaging for organ-level lineage tracing of rare intestinal stem cell populations

Significance

Lineage tracing using fluorescent reporters is a common tool for monitoring the expression of genes and transcription factors in stem cell populations and their progeny. The zinc-binding protein 89 (ZBP-89/Zfp148 mouse gene) is a transcription factor that plays a role in gastrointestinal (GI) stem cell maintenance and cellular differentiation and has been linked to the progression of colon cancer. While lineage tracing is a useful tool, it is commonly performed with high-magnification microscopy on a small field of view within tissue sections, thereby limiting the ability to resolve reporter expression at the organ level. Furthermore, this technique requires extensive tissue processing, which is time consuming and requires euthanizing the animal. Further knowledge could be elucidated by measuring the expression of fluorescent reporters across entire organs with minimal tissue processing.

Aim

We present the application of wide-field fluorescence imaging for whole-organ lineage tracing of an inducible Zfp148-tdTomato-expressing transgenic mouse line to assess the expression of ZBP-89/Zfp148 in the GI tract.

Approach

We measured tdTomato fluorescence in ex vivo organs at time points between 24 h and 6 months post-induction. Fluctuations in tdTomato expression were validated by fluorescence microscopy of tissue sections.

Results

Quantification of the wide field-of-view images showed a statistically significant increase in fluorescent signal across the GI tract between transgenic mice and littermate controls. The results also showed a gradient of decreasing reporter expression from proximal to distal intestine, suggesting a higher abundance of ZBP-89 expressing stem cells, or higher expression of ZBP-89 within the stem cells, in the proximal intestine.

Conclusions

We demonstrate that wide-field fluorescence imaging is a valuable tool for monitoring whole-organ expression of fluorescent reporters. This technique could potentially be applied in vivo for longitudinal assessment of a single animal, further enhancing our ability to resolve rare stem cell lineages spatially and temporally.

Hierarchical deconvolution for extensive cell type resolution in the human brain using DNA methylation

Introduction

The human brain comprises heterogeneous cell types whose composition can be altered with physiological and pathological conditions. New approaches to discern the diversity and distribution of brain cells associated with neurological conditions would significantly advance the study of brain-related pathophysiology and neuroscience. Unlike single-nuclei approaches, DNA methylation-based deconvolution does not require special sample handling or processing, is cost-effective, and easily scales to large study designs. Existing DNA methylation-based methods for brain cell deconvolution are limited in the number of cell types deconvolved.

Methods

Using DNA methylation profiles of the top cell-type-specific differentially methylated CpGs, we employed a hierarchical modeling approach to deconvolve GABAergic neurons, glutamatergic neurons, astrocytes, microglial cells, oligodendrocytes, endothelial cells, and stromal cells.

Results

We demonstrate the utility of our method by applying it to data on normal tissues from various brain regions and in aging and diseased tissues, including Alzheimer's disease, autism, Huntington's disease, epilepsy, and schizophrenia.

Discussion

We expect that the ability to determine the cellular composition in the brain using only DNA from bulk samples will accelerate understanding brain cell type composition and cell-type-specific epigenetic states in normal and diseased brain tissues.

Down Syndrome Altered Cell Composition in Blood, Brain, and Buccal Swab Samples Profiled by DNA-Methylation-Based Cell-Type Deconvolution

Down syndrome (DS) is a genetic disorder caused by an extra copy of chromosome 21 that presents developmental dysfunction and intellectual disability. To better understand the cellular changes associated with DS, we investigated the cell composition in blood, brain, and buccal swab samples from DS patients and controls using DNA methylation-based cell-type deconvolution. We used genome-scale DNA methylation data from Illumina HumanMethylation450k and HumanMethylationEPIC arrays to profile cell composition and trace fetal lineage cells in blood samples (DS N = 46; control N = 1469), brain samples from various regions (DS N = 71; control N = 101), and buccal swab samples (DS N = 10; control N = 10). In early development, the number of cells from the fetal lineage in the blood is drastically lower in DS patients (Δ = 17.5%), indicating an epigenetically dysregulated maturation process for DS patients. Across sample types, we observed significant alterations in relative cell-type proportions for DS subjects compared with the controls. Cell-type proportion alterations were present in samples from early development and adulthood. Our findings provide insight into DS cellular biology and suggest potential cellular interventional targets for DS.

Evaluation of cross-platform compatibility of a DNA methylation-based glucocorticoid response biomarker

Background

Identifying blood-based DNA methylation patterns is a minimally invasive way to detect biomarkers in predicting age, characteristics of certain diseases and conditions, as well as responses to immunotherapies. As microarray platforms continue to evolve and increase the scope of CpGs measured, new discoveries based on the most recent platform version and how they compare to available data from the previous versions of the platform are unknown. The neutrophil dexamethasone methylation index (NDMI 850) is a blood-based DNA methylation biomarker built on the Illumina MethylationEPIC (850K) array that measures epigenetic responses to dexamethasone (DEX), a synthetic glucocorticoid often administered for inflammation. Here, we compare the NDMI 850 to one we built using data from the Illumina Methylation 450K (NDMI 450).

Results

The NDMI 450 consisted of 22 loci, 15 of which were present on the NDMI 850. In adult whole blood samples, the linear composite scores from NDMI 450 and NDMI 850 were highly correlated and had equivalent predictive accuracy for detecting DEX exposure among adult glioma patients and non-glioma adult controls. However, the NDMI 450 scores of newborn cord blood were significantly lower than NDMI 850 in samples measured with both assays.

Conclusions

We developed an algorithm that reproduces the DNA methylation glucocorticoid response score using 450K data, increasing the accessibility for researchers to assess this biomarker in archived or publicly available datasets that use the 450K version of the Illumina BeadChip array. However, the NDMI850 and NDMI450 do not give similar results in cord blood, and due to data availability limitations, results from sample types of newborn cord blood should be interpreted with care.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01352-1.

Lineage tracing in human tissues

The dynamical process of cell division that underpins homeostasis in the human body cannot be directly observed in vivo, but instead is measurable from the pattern of somatic genetic or epigenetic mutations that accrue in tissues over an individual's lifetime. Because somatic mutations are heritable, they serve as natural lineage tracing markers that delineate clonal expansions. Mathematical analysis of the distribution of somatic clone sizes gives a quantitative readout of the rates of cell birth, death, and replacement. In this review we explore the broad range of somatic mutation types that have been used for lineage tracing in human tissues, introduce the mathematical concepts used to infer dynamical information from these clone size data, and discuss the insights of this lineage tracing approach for our understanding of homeostasis and cancer development. We use the human colon as a particularly instructive exemplar tissue. There is a rich history of human somatic cell dynamics surreptitiously written into the cell genomes that is being uncovered by advances in sequencing and careful mathematical analysis lineage of tracing data. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

CancerNet: a unified deep learning network for pan-cancer diagnostics

Background

Despite remarkable advances in cancer research, cancer remains one of the leading causes of death worldwide. Early detection of cancer and localization of the tissue of its origin are key to effective treatment. Here, we leverage technological advances in machine learning or artificial intelligence to design a novel framework for cancer diagnostics. Our proposed framework detects cancers and their tissues of origin using a unified model of cancers encompassing 33 cancers represented in The Cancer Genome Atlas (TCGA). Our model exploits the learned features of different cancers reflected in the respective dysregulated epigenomes, which arise early in carcinogenesis and differ remarkably between different cancer types or subtypes, thus holding a great promise in early cancer detection.

Results

Our comprehensive assessment of the proposed model on the 33 different tissues of origin demonstrates its ability to detect and classify cancers to a high accuracy (> 99% overall F-measure). Furthermore, our model distinguishes cancers from pre-cancerous lesions to metastatic tumors and discriminates between hypomethylation changes due to age related epigenetic drift and true cancer.

Conclusions

Beyond detection of primary cancers, our proposed computational model also robustly detects tissues of origin of secondary cancers, including metastatic cancers, second primary cancers, and cancers of unknown primaries. Our assessment revealed the ability of this model to characterize pre-cancer samples, a significant step forward in early cancer detection. Deployed broadly this model can deliver accurate diagnosis for a greatly expanded target patient population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04783-y.

Interactions of Age and Blood Immune Factors and Noninvasive Prediction of Glioma Survival

BACKGROUND

Tumor-based classification of human glioma portends patient prognosis, but considerable unexplained survival variability remains. Host factors (eg, age) also strongly influence survival times, partly reflecting a compromised immune system. How blood epigenetic measures of immune characteristics and age augment molecular classifications in glioma survival has not been investigated. We assess the prognostic impact of immune cell fractions and epigenetic age in archived blood across glioma molecular subtypes for the first time.

METHODS

We evaluated immune cell fractions and epigenetic age in archived blood from the University of California San Francisco Adult Glioma Study, which included a training set of 197 patients with IDH-wild type, 1p19q intact, TERT wild type (IDH/1p19q/TERT-WT) glioma, an evaluation set of 350 patients with other subtypes of glioma, and 454 patients without glioma.

RESULTS

IDH/1p19q/TERT-WT patients had lower lymphocyte fractions (CD4+ T, CD8+ T, natural killer, and B cells) and higher neutrophil fractions than people without glioma. Recursive partitioning analysis delineated 4 statistically significantly different survival groups for patients with IDH/1p19q/TERT-WT based on an interaction between chronological age and 2 blood immune factors, CD4+ T cells, and neutrophils. Median overall survival ranged from 0.76 years (95% confidence interval = 0.55-0.99) for the worst survival group (n = 28) to 9.72 years (95% confidence interval = 6.18 to not available) for the best (n = 33). The recursive partitioning analysis also statistically significantly delineated 4 risk groups in patients with other glioma subtypes.

CONCLUSIONS

The delineation of different survival groups in the training and evaluation sets based on an interaction between chronological age and blood immune characteristics suggests that common host immune factors among different glioma types may affect survival. The ability of DNA methylation-based markers of immune status to capture diverse, clinically relevant information may facilitate noninvasive, personalized patient evaluation in the neuro-oncology clinic.

Experimental development of the epigenomic library construction method to elucidate the epigenetic diversity and causal relationship between epigenome and transcriptome at a single-cell level

The method of single-cell RNA sequencing has been rapidly developed, and numerous experiments have been conducted over the past decade. Their results allow us to recognize various subpopulations and rare cell states in tissues, tumors, and immune systems that are previously unidentified, and guide us to understand fundamental biological processes that determine cell identity based on single-cell gene expression profiles. However, it is still challenging to understand the principle of comprehensive gene regulation that determines the cell fate only with transcriptome, a consequential output of the gene expression program. To elucidate the mechanisms related to the origin and maintenance of comprehensive single-cell transcriptome, we require a corresponding single-cell epigenome, which is a differentiated information of each cell with an identical genome. This review deals with the current development of single-cell epigenomic library construction methods, including multi-omics tools with crucial factors and additional requirements in the future focusing on DNA methylation, chromatin accessibility, and histone post-translational modifications. The study of cellular differentiation and the disease occurrence at a single-cell level has taken the first step with single-cell transcriptome and is now taking the next step with single-cell epigenome.

Immune profiles and DNA methylation alterations related with non-muscle-invasive bladder cancer outcomes

BACKGROUND

Non-muscle-invasive bladder cancer (NMIBC) patients receive frequent monitoring because ≥ 70% will have recurrent disease. However, screening is invasive, expensive, and associated with significant morbidity making bladder cancer the most expensive cancer to treat per capita. There is an urgent need to expand the understanding of markers related to recurrence and survival outcomes of NMIBC.

METHODS AND RESULTS

We used the Illumina HumanMethylationEPIC array to measure peripheral blood DNA methylation profiles of NMIBC patients (N = 603) enrolled in a population-based cohort study in New Hampshire and applied cell type deconvolution to estimate immune cell-type proportions. Using Cox proportional hazard models, we identified that increasing CD4T and CD8T cell proportions were associated with a statistically significant decreased hazard of tumor recurrence or death (CD4T: HR = 0.98, 95% CI = 0.97-1.00; CD8T: HR = 0.97, 95% CI = 0.95-1.00), whereas increasing monocyte proportion and methylation-derived neutrophil-to-lymphocyte ratio (mdNLR) were associated with the increased hazard of tumor recurrence or death (monocyte: HR = 1.04, 95% CI = 1.00-1.07; mdNLR: HR = 1.12, 95% CI = 1.04-1.20). Then, using an epigenome-wide association study (EWAS) approach adjusting for age, sex, smoking status, BCG treatment status, and immune cell profiles, we identified 2528 CpGs associated with the hazard of tumor recurrence or death (P < 0.005). Among these CpGs, the 1572 were associated with an increased hazard and were significantly enriched in open sea regions; the 956 remaining CpGs were associated with a decreased hazard and were significantly enriched in enhancer regions and DNase hypersensitive sites.

CONCLUSIONS

Our results expand on the knowledge of immune profiles and methylation alteration associated with NMIBC outcomes and represent a first step toward the development of DNA methylation-based biomarkers of tumor recurrence.

A transdisciplinary approach to understand the epigenetic basis of race/ethnicity health disparities

Health disparities correspond to differences in disease burden and mortality among socially defined population groups. Such disparities may emerge according to race/ethnicity, socioeconomic status and a variety of other social contexts, and are documented for a wide range of diseases. Here, we provide a transdisciplinary perspective on the contribution of epigenetics to the understanding of health disparities, with a special emphasis on disparities across socially defined racial/ethnic groups. Scientists in the fields of biological anthropology, bioinformatics and molecular epidemiology provide a summary of theoretical, statistical and practical considerations for conducting epigenetic health disparities research, and provide examples of successful applications from cancer research using this approach.

An Optimized Method for the Construction of a DNA Methylome from Small Quantities of Tissue or Purified DNA from Arabidopsis Embryo

DNA methylation is an important epigenetic mechanism affecting genome structure, gene regulation, and the silencing of transposable elements. Cell- and tissue-specific methylation patterns are critical for differentiation and development in eukaryotes. Dynamic spatiotemporal methylation data in these cells or tissues is, therefore, of great interest. However, the construction of bisulfite sequencing libraries can be challenging if the starting material is limited or the genome size is small, such as in Arabidopsis. Here, we describe detailed methods for the purification of Arabidopsis embryos at all stages, and the construction of comprehensive bisulfite libraries from small quantities of input. We constructed bisulfite libraries by releasing embryos from intact seeds, using a different approach for each developmental stage, and manually picking single-embryo with microcapillaries. From these libraries, reliable Arabidopsis methylome data were collected allowing, on average, 11-fold coverage of the genome using as few as five globular, heart, and torpedo embryos as raw input material without the need for DNA purification step. On the other hand, purified DNA from as few as eight bending torpedo embryos or a single mature embryo is sufficient for library construction when RNase A is treated before DNA extraction. This method can be broadly applied to cells from different tissues or cells from other model organisms. Methylome construction can be achieved using a minimal amount of input material using our method; thereby, it has the potential to increase our understanding of dynamic spatiotemporal methylation patterns in model organisms.

Differential DNA methylation and changing cell-type proportions as fibrotic stage progresses in NAFLD

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is characterized by changes in cell composition that occur throughout disease pathogenesis, which includes the development of fibrosis in a subset of patients. DNA methylation (DNAm) is a plausible mechanism underlying these shifts, considering that DNAm profiles differ across tissues and cell types, and DNAm may play a role in cell-type differentiation. Previous work investigating the relationship between DNAm and fibrosis in NAFLD has been limited by sample size and the number of CpG sites interrogated.

RESULTS

Here, we performed an epigenome-wide analysis using Infinium MethylationEPIC array data from 325 individuals with NAFLD, including 119 with severe fibrosis and 206 with no histological evidence of fibrosis. After adjustment for latent confounders, we identified 7 CpG sites whose DNAm associated with fibrosis (p < 5.96 × 10^-8). Analysis of RNA-seq data collected from a subset of individuals (N = 56) revealed that gene expression at 288 genes associated with DNAm at one or more of the 7 fibrosis-related CpGs. DNAm-based estimates of cell-type proportions showed that estimated proportions of natural killer cells increased, while epithelial cell proportions decreased with disease stage. Finally, we used an elastic net regression model to assess DNAm as a biomarker of fibrotic stage and found that our model predicted fibrosis with a sensitivity of 0.93 and provided information beyond a model based solely on cell-type proportions.

CONCLUSION

These findings are consistent with DNAm as a mechanism underpinning or marking fibrosis-related shifts in cell composition and demonstrate the potential of DNAm as a possible biomarker of NAFLD fibrosis.

Estimands in epigenome-wide association studies

Background

In DNA methylation analyses like epigenome-wide association studies, effects in differentially methylated CpG sites are assessed. Two kinds of outcomes can be used for statistical analysis: Beta-values and M-values. M-values follow a normal distribution and help to detect differentially methylated CpG sites. As biological effect measures, differences of M-values are more or less meaningless. Beta-values are of more interest since they can be interpreted directly as differences in percentage of DNA methylation at a given CpG site, but they have poor statistical properties. Different frameworks are proposed for reporting estimands in DNA methylation analysis, relying on Beta-values, M-values, or both.

Results

We present and discuss four possible approaches of achieving estimands in DNA methylation analysis. In addition, we present the usage of M-values or Beta-values in the context of bioinformatical pipelines, which often demand a predefined outcome. We show the dependencies between the differences in M-values to differences in Beta-values in two data simulations: a analysis with and without confounder effect. Without present confounder effects, M-values can be used for the statistical analysis and Beta-values statistics for the reporting. If confounder effects exist, we demonstrate the deviations and correct the effects by the intercept method. Finally, we demonstrate the theoretical problem on two large human genome-wide DNA methylation datasets to verify the results.

Conclusions

The usage of M-values in the analysis of DNA methylation data will produce effect estimates, which cannot be biologically interpreted. The parallel usage of Beta-value statistics ignores possible confounder effects and can therefore not be recommended. Hence, if the differences in Beta-values are the focus of the study, the intercept method is recommendable. Hyper- or hypomethylated CpG sites must then be carefully evaluated. If an exploratory analysis of possible CpG sites is the aim of the study, M-values can be used for inference.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01083-9.

No Time to Age: Uncoupling Aging from Chronological Time

Multicellular life evolved from simple unicellular organisms that could replicate indefinitely, being essentially ageless. At this point, life split into two fundamentally different cell types: the immortal germline representing an unbroken lineage of cell division with no intrinsic endpoint and the mortal soma, which ages and dies. In this review, we describe the germline as clock-free and the soma as clock-bound and discuss aging with respect to three DNA-based cellular clocks (telomeric, DNA methylation, and transposable element). The ticking of these clocks corresponds to the stepwise progressive limitation of growth and regeneration of somatic cells that we term somatic restriction. Somatic restriction acts in opposition to strategies that ensure continued germline replication and regeneration. We thus consider the plasticity of aging as a process not fixed to the pace of chronological time but one that can speed up or slow down depending on the rate of intrinsic cellular clocks. We further describe how germline factor reprogramming might be used to slow the rate of aging and potentially reverse it by causing the clocks to tick backward. Therefore, reprogramming may eventually lead to therapeutic strategies to treat degenerative diseases by altering aging itself, the one condition common to us all.

Identification of a foetal epigenetic compartment in adult human kidney

The mammalian kidney has extensive repair capacity; however, identifying adult renal stem cells has proven elusive. We applied an epigenetic marker of foetal cell origin (FCO) in diverse human tissues as a probe for developmental cell persistence, finding a 5.4-fold greater FCO proportion in kidney. Normal kidney FCO proportions averaged 49% with extensive interindividual variation. FCO proportions were significantly negatively correlated with immune-related gene expression and positively correlated with genes expressed in the renal medulla, including those involved in renal organogenesis (e.g., FGF2, PAX8, and HOXB7). FCO associated genes also mapped to medullary nephron segments in mouse and rat, suggesting evolutionary conservation of this cellular compartment. Renal cancer patients whose tumours contained non-zero FCO scores survived longer. The kidney appears unique in possessing substantial foetal epigenetic features. Further study of FCO-related gene methylation may elucidate regenerative regulatory programmes in tissues without apparent discrete stem cell compartments.

DNA methylation in blood-Potential to provide new insights into cell biology

Epigenetics plays a fundamental role in cellular development and differentiation; epigenetic mechanisms, such as DNA methylation, are involved in gene regulation and the exquisite nuance of expression changes seen in the journey from pluripotency to final differentiation. Thus, DNA methylation as a marker of cell identify has the potential to reveal new insights into cell biology. We mined publicly available DNA methylation data with a machine-learning approach to identify differentially methylated loci between different white blood cell types. We then interrogated the DNA methylation and mRNA expression of candidate loci in CD4⁺, CD8⁺, CD14⁺, CD19⁺ and CD56⁺ fractions from 12 additional, independent healthy individuals (6 male, 6 female). ‘Classic’ immune cell markers such as CD8 and CD19 showed expected methylation/expression associations fitting with established dogma that hypermethylation is associated with the repression of gene expression. We also observed large differential methylation at loci which are not established immune cell markers; some of these loci showed inverse correlations between methylation and mRNA expression (such as PARK2, DCP2). Furthermore, we validated these observations further in publicly available DNA methylation and RNA sequencing datasets. Our results highlight the value of mining publicly available data, the utility of DNA methylation as a discriminatory marker and the potential value of DNA methylation to provide additional insights into cell biology and developmental processes.

Cell Lineage Tracing and Cellular Diversity in Humans

Tracing cell lineages is fundamental for understanding the rules governing development in multicellular organisms and delineating complex biological processes involving the differentiation of multiple cell types with distinct lineage hierarchies. In humans, experimental lineage tracing is unethical, and one has to rely on natural-mutation markers that are created within cells as they proliferate and age. Recent studies have demonstrated that it is now possible to trace lineages in normal, noncancerous cells with a variety of data types using natural variations in the nuclear and mitochondrial DNA as well as variations in DNA methylation status. It is also apparent that the scientific community is on the verge of being able to make a comprehensive and detailed cell lineage map of human embryonic and fetal development. In this review, we discuss the advantages and disadvantages of different approaches and markers for lineage tracing. We also describe the general conceptual design for how to derive a lineage map for humans.

MethylNet: an automated and modular deep learning approach for DNA methylation analysis

BACKGROUND

DNA methylation (DNAm) is an epigenetic regulator of gene expression programs that can be altered by environmental exposures, aging, and in pathogenesis. Traditional analyses that associate DNAm alterations with phenotypes suffer from multiple hypothesis testing and multi-collinearity due to the high-dimensional, continuous, interacting and non-linear nature of the data. Deep learning analyses have shown much promise to study disease heterogeneity. DNAm deep learning approaches have not yet been formalized into user-friendly frameworks for execution, training, and interpreting models. Here, we describe MethylNet, a DNAm deep learning method that can construct embeddings, make predictions, generate new data, and uncover unknown heterogeneity with minimal user supervision.

RESULTS

The results of our experiments indicate that MethylNet can study cellular differences, grasp higher order information of cancer sub-types, estimate age and capture factors associated with smoking in concordance with known differences.

CONCLUSION

The ability of MethylNet to capture nonlinear interactions presents an opportunity for further study of unknown disease, cellular heterogeneity and aging processes.

Unravelling cellular relationships during development and regeneration using genetic lineage tracing

Tracking the progeny of single cells is necessary for building lineage trees that recapitulate processes such as embryonic development and stem cell differentiation. In classical lineage tracing experiments, cells are fluorescently labelled to allow identification by microscopy of a limited number of cell clones. To track a larger number of clones in complex tissues, fluorescent proteins are now replaced by heritable DNA barcodes that are read using next-generation sequencing. In prospective lineage tracing, unique DNA barcodes are introduced into single cells through genetic manipulation (using, for example, Cre-mediated recombination or CRISPR-Cas9-mediated editing) and tracked over time. Alternatively, in retrospective lineage tracing, naturally occurring somatic mutations can be used as endogenous DNA barcodes. Finally, single-cell mRNA-sequencing datasets that capture different cell states within a developmental or differentiation trajectory can be used to recapitulate lineages. In this Review, we discuss methods for prospective or retrospective lineage tracing and demonstrate how trajectory reconstruction algorithms can be applied to single-cell mRNA-sequencing datasets to infer developmental or differentiation tracks. We discuss how these approaches are used to understand cell fate during embryogenesis, cell differentiation and tissue regeneration.

Absence of an embryonic stem cell DNA methylation signature in human cancer

BACKGROUND

Differentiated cells that arise from stem cells in early development contain DNA methylation features that provide a memory trace of their fetal cell origin (FCO). The FCO signature was developed to estimate the proportion of cells in a mixture of cell types that are of fetal origin and are reminiscent of embryonic stem cell lineage. Here we implemented the FCO signature estimation method to compare the fraction of cells with the FCO signature in tumor tissues and their corresponding nontumor normal tissues.

METHODS

We applied our FCO algorithm to discovery data sets obtained from The Cancer Genome Atlas (TCGA) and replication data sets obtained from the Gene Expression Omnibus (GEO) data repository. Wilcoxon rank sum tests, linear regression models with adjustments for potential confounders and non-parametric randomization-based tests were used to test the association of FCO proportion between tumor tissues and nontumor normal tissues. P-values of < 0.05 were considered statistically significant.

RESULTS

Across 20 different tumor types we observed a consistently lower FCO signature in tumor tissues compared with nontumor normal tissues, with 18 observed to have significantly lower FCO fractions in tumor tissue (total n = 6,795 tumor, n = 922 nontumor, P < 0.05). We replicated our findings in 15 tumor types using data from independent subjects in 15 publicly available data sets (total n = 740 tumor, n = 424 nontumor, P < 0.05).

CONCLUSIONS

The results suggest that cancer development itself is substantially devoid of recapitulation of normal embryologic processes. Our results emphasize the distinction between DNA methylation in normal tightly regulated stem cell driven differentiation and cancer stem cell reprogramming that involves altered methylation in the service of great cell heterogeneity and plasticity.

Recording development with single cell dynamic lineage tracing

Every animal grows from a single fertilized egg into an intricate network of cell types and organ systems. This process is captured in a lineage tree: a diagram of every cell's ancestry back to the founding zygote. Biologists have long sought to trace this cell lineage tree in individual organisms and have developed a variety of technologies to map the progeny of specific cells. However, there are billions to trillions of cells in complex organisms, and conventional approaches can only map a limited number of clonal populations per experiment. A new generation of tools that use molecular recording methods integrated with single cell profiling technologies may provide a solution. Here, we summarize recent breakthroughs in these technologies, outline experimental and computational challenges, and discuss biological questions that can be addressed using single cell dynamic lineage tracing.