Cell-type deconvolution from DNA methylation: a review of recent applications

Inflammatory breast cancer microenvironment repertoire based on DNA methylation data deconvolution reveals actionable targets to enhance the treatment efficacy

BACKGROUND

Although the clinical signs of inflammatory breast cancer (IBC) resemble acute inflammation, the role played by infiltrating immune and stromal cells in this aggressive disease is uncharted. The tumor microenvironment (TME) presents molecular alterations, such as epimutations, prior to morphological abnormalities. These changes affect the distribution and the intricate communication between the TME components related to cancer prognosis and therapy response. Herein, we explored the global DNA methylation profile of IBC and surrounding tissues to estimate the microenvironment cellular composition and identify epigenetically dysregulated markers.

METHODS

We used the HiTIMED algorithm to deconvolve the bulk DNA methylation data of 24 IBC and six surrounding non-tumoral tissues (SNT) (GSE238092) and determine their cellular composition. The prognostic relevance of cell types infiltrating IBC and their relationship with clinicopathological variables were investigated. CD34 (endothelial cell marker) and CD68 (macrophage marker) immunofluorescence staining was evaluated in an independent set of 17 IBC and 16 non-IBC samples.

RESULTS

We found lower infiltration of endothelial, stromal, memory B, dendritic, and natural killer cells in IBC than in SNT samples. Higher endothelial cell (EC) and stromal cell content were related to better overall survival. EC proportions positively correlated with memory B and memory CD8+ T infiltration in IBC. Immune and EC markers exhibited distinct DNA methylation profiles between IBC and SNT samples, revealing hypermethylated regions mapped to six genes (CD40, CD34, EMCN, HLA-G, PDPN, and TEK). We identified significantly higher CD34 and CD68 protein expression in IBC compared to non-IBC.

CONCLUSIONS

Our findings underscored cell subsets that distinguished patients with better survival and dysregulated markers potentially actionable through combinations of immunotherapy and epigenetic drugs.

An assessment of compositional methods for the analysis of DNA methylation-based deconvolution estimates

DNA methylation (DNAm)-based deconvolution estimates contain relative data, forming a composition, that standard methods (testing directly on cell proportions) are ill-suited to handle. In this study we examined the performance of an alternative method, analysis of compositions of microbiomes (ANCOM), for the analysis of DNAm-based deconvolution estimates. We performed two different simulation studies comparing ANCOM to a standard approach (two sample t-test performed directly on cell proportions) and analyzed a real-world data from the Women's Health Initiative to evaluate the applicability of ANCOM to DNAm-based deconvolution estimates. Our findings indicate that ANCOM can effectively account for the compositional nature of DNAm-based deconvolution estimates. ANCOM adequately controls the false discovery rate while maintaining statistical power comparable to that of standard methods.

Biological aging of different blood cell types

Biological age (BA) captures detrimental age-related changes. The best-known and most-used BA indicators include DNA methylation-based epigenetic clocks and telomere length (TL). The most common biological sample material for epidemiological aging studies, whole blood, is composed of different cell types. We aimed to compare differences in BAs between blood cell types and assessed the BA indicators' cell type-specific associations with chronological age (CA). An analysis of DNA methylation-based BA indicators, including TL, methylation level at cg16867657 in ELOVL2, as well as the Hannum, Horvath, DNAmPhenoAge, and DunedinPACE epigenetic clocks, was performed on 428 biological samples of 12 blood cell types. BA values were different in the majority of the pairwise comparisons between cell types, as well as in comparison to whole blood (p < 0.05). DNAmPhenoAge showed the largest cell type differences, up to 44.5 years and DNA methylation-based TL showed the lowest differences. T cells generally had the "youngest" BA values, with differences across subsets, whereas monocytes had the "oldest" values. All BA indicators, except DunedinPACE, strongly correlated with CA within a cell type. Some differences such as DNAmPhenoAge-difference between naïve CD4 + T cells and monocytes were constant regardless of the blood donor's CA (range 20-80 years), while for DunedinPACE they were not. In conclusion, DNA methylation-based indicators of BA exhibit cell type-specific characteristics. Our results have implications for understanding the molecular mechanisms underlying epigenetic clocks and underscore the importance of considering cell composition when utilizing them as indicators for the success of aging interventions.

Tissue of origin detection for cancer tumor using low-depth cfDNA samples through combination of tumor-specific methylation atlas and genome-wide methylation density in graph convolutional neural networks

BACKGROUND

Cell free DNA (cfDNA)-based assays hold great potential in detecting early cancer signals yet determining the tissue-of-origin (TOO) for cancer signals remains a challenging task. Here, we investigated the contribution of a methylation atlas to TOO detection in low depth cfDNA samples.

METHODS

We constructed a tumor-specific methylation atlas (TSMA) using whole-genome bisulfite sequencing (WGBS) data from five types of tumor tissues (breast, colorectal, gastric, liver and lung cancer) and paired white blood cells (WBC). TSMA was used with a non-negative least square matrix factorization (NNLS) deconvolution algorithm to identify the abundance of tumor tissue types in a WGBS sample. We showed that TSMA worked well with tumor tissue but struggled with cfDNA samples due to the overwhelming amount of WBC-derived DNA. To construct a model for TOO, we adopted the multi-modal strategy and used as inputs the combination of deconvolution scores from TSMA with other features of cfDNA.

RESULTS

Our final model comprised of a graph convolutional neural network using deconvolution scores and genome-wide methylation density features, which achieved an accuracy of 69% in a held-out validation dataset of 239 low-depth cfDNA samples.

CONCLUSIONS

In conclusion, we have demonstrated that our TSMA in combination with other cfDNA features can improve TOO detection in low-depth cfDNA samples.

GLDADec: marker-gene guided LDA modeling for bulk gene expression deconvolution

Inferring cell type proportions from bulk transcriptome data is crucial in immunology and oncology. Here, we introduce guided LDA deconvolution (GLDADec), a bulk deconvolution method that guides topics using cell type-specific marker gene names to estimate topic distributions for each sample. Through benchmarking using blood-derived datasets, we demonstrate its high estimation performance and robustness. Moreover, we apply GLDADec to heterogeneous tissue bulk data and perform comprehensive cell type analysis in a data-driven manner. We show that GLDADec outperforms existing methods in estimation performance and evaluate its biological interpretability by examining enrichment of biological processes for topics. Finally, we apply GLDADec to The Cancer Genome Atlas tumor samples, enabling subtype stratification and survival analysis based on estimated cell type proportions, thus proving its practical utility in clinical settings. This approach, utilizing marker gene names as partial prior information, can be applied to various scenarios for bulk data deconvolution. GLDADec is available as an open-source Python package at https://github.com/mizuno-group/GLDADec.

Investigating neonatal health risk variables through cell-type specific methylome-wide association studies

Adverse neonatal outcomes are a prevailing risk factor for both short- and long-term mortality and morbidity in infants. Given the importance of these outcomes, refining their assessment is paramount for improving prevention and care. Here we aim to enhance the assessment of these often correlated and multifaceted neonatal outcomes. To achieve this, we employ factor analysis to identify common and unique effects and further confirm these effects using criterion-related validity testing. This validation leverages methylome-wide profiles from neonatal blood. Specifically, we investigate nine neonatal health risk variables, including gestational age, Apgar score, three indicators of body size, jaundice, birth diagnosis, maternal preeclampsia, and maternal age. The methylomic profiles used for this research capture data from nearly all 28 million methylation sites in human blood, derived from the blood spot collected from 333 neonates, within 72 h post-birth. Our factor analysis revealed two common factors, size factor, that captured the shared effects of weight, head size, height, and gestational age and disease factor capturing the orthogonal shared effects of gestational age, combined with jaundice and birth diagnosis. To minimize false positives in the validation studies, validation was limited to variables with significant cumulative association as estimated through an in-sample replication procedure. This screening resulted in that the two common factors and the unique effects for gestational age, jaundice and Apgar were further investigated with full-scale cell-type specific methylome-wide association analyses. Highly significant, cell-type specific, associations were detected for both common effect factors and for Apgar. Gene Ontology analyses revealed multiple significant biologically relevant terms for the five fully investigated neonatal health risk variables. Given the established links between adverse neonatal outcomes and both immediate and long-term health, the distinct factor effects (representing the common and unique effects of the risk variables) and their biological profiles confirmed in our work, suggest their potential role as clinical biomarkers for assessing health risks and enhancing personalized care.

Immunomethylomic profiles of long-term head and neck squamous cell carcinoma survivors on immune checkpoint inhibitors

Aim: This study addresses the challenge of predicting the response of head and neck squamous cell carcinoma (HNSCC) patients to immunotherapy. Methods: Using DNA methylation cytometry, we analyzed the immune profiles of six HNSCC patients who showed a positive response to immunotherapy over a year without disease progression. Results: There was an initial increase in CD8 T memory cells and natural killer cells during the first four cycles of immunotherapy, which then returned to baseline levels after a year. Baseline CD8 T cell levels were lower in HNSCC immunotherapy responders but became similar to those in healthy subjects after immunotherapy. Conclusion: These findings suggest that monitoring fluctuations in immune profiles could potentially identify biomarkers for immunotherapy response in HNSCC patients.

Benchmarking of methods for DNA methylome deconvolution

Defining the number and abundance of different cell types in tissues is important for understanding disease mechanisms as well as for diagnostic and prognostic purposes. Typically, this is achieved by immunohistological analyses, cell sorting, or single-cell RNA-sequencing. Alternatively, cell-specific DNA methylome information can be leveraged to deconvolve cell fractions from a bulk DNA mixture. However, comprehensive benchmarking of deconvolution methods and modalities was not yet performed. Here we evaluate 16 deconvolution algorithms, developed either specifically for DNA methylome data or more generically. We assess the performance of these algorithms, and the effect of normalization methods, while modeling variables that impact deconvolution performance, including cell abundance, cell type similarity, reference panel size, method for methylome profiling (array or sequencing), and technical variation. We observe differences in algorithm performance depending on each these variables, emphasizing the need for tailoring deconvolution analyses. The complexity of the reference, marker selection method, number of marker loci and, for sequencing-based assays, sequencing depth have a marked influence on performance. By developing handles to select the optimal analysis configuration, we provide a valuable source of information for studies aiming to deconvolve array- or sequencing-based methylation data.

DNA Methylation-Derived Immune Cell Proportions and Cancer Risk, Including Lung Cancer, in Black Participants

Prior cohort studies assessing cancer risk based on immune cell subtype profiles have predominantly focused on White populations. This limitation obscures vital insights into how cancer risk varies across race. Immune cell subtype proportions were estimated using deconvolution based on leukocyte DNA methylation markers from blood samples collected at baseline on participants without cancer in the Atherosclerosis Risk in Communities (ARIC) Study. Over a mean of 17.5 years of follow-up, 668 incident cancers were diagnosed in 2,467 Black participants. Cox proportional hazards regression was used to examine immune cell subtype proportions and overall cancer incidence and site-specific incidence (lung, breast, and prostate cancers). Higher T regulatory cell proportions were associated with statistically significantly higher lung cancer risk (hazard ratio = 1.22, 95% confidence interval = 1.06-1.41 per percent increase). Increased memory B cell proportions were associated with significantly higher risk of prostate cancer (1.17, 1.04-1.33) and all cancers (1.13, 1.05-1.22). Increased CD8+ naïve cell proportions were associated with significantly lower risk of all cancers in participants ≥55 years (0.91, 0.83-0.98). Other immune cell subtypes did not display statistically significant associations with cancer risk. These results in Black participants align closely with prior findings in largely White populations. Findings from this study could help identify those at high cancer risk and outline risk stratifying to target patients for cancer screening, prevention, and other interventions. Further studies should assess these relationships in other cancer types, better elucidate the interplay of B cells in cancer risk, and identify biomarkers for personalized risk stratification.

Computational deconvolution of DNA methylation data from mixed DNA samples

In this review, we provide a comprehensive overview of the different computational tools that have been published for the deconvolution of bulk DNA methylation (DNAm) data. Here, deconvolution refers to the estimation of cell-type proportions that constitute a mixed sample. The paper reviews and compares 25 deconvolution methods (supervised, unsupervised or hybrid) developed between 2012 and 2023 and compares the strengths and limitations of each approach. Moreover, in this study, we describe the impact of the platform used for the generation of methylation data (including microarrays and sequencing), the applied data pre-processing steps and the used reference dataset on the deconvolution performance. Next to reference-based methods, we also examine methods that require only partial reference datasets or require no reference set at all. In this review, we provide guidelines for the use of specific methods dependent on the DNA methylation data type and data availability.

Cell type signatures in cell-free DNA fragmentation profiles reveal disease biology

Circulating cell-free DNA (cfDNA) fragments have characteristics that are specific to the cell types that release them. Current methods for cfDNA deconvolution typically use disease tailored marker selection in a limited number of bulk tissues or cell lines. Here, we utilize single cell transcriptome data as a comprehensive cellular reference set for disease-agnostic cfDNA cell-of-origin analysis. We correlate cfDNA-inferred nucleosome spacing with gene expression to rank the relative contribution of over 490 cell types to plasma cfDNA. In 744 healthy individuals and patients, we uncover cell type signatures in support of emerging disease paradigms in oncology and prenatal care. We train predictive models that can differentiate patients with colorectal cancer (84.7%), early-stage breast cancer (90.1%), multiple myeloma (AUC 95.0%), and preeclampsia (88.3%) from matched controls. Importantly, our approach performs well in ultra-low coverage cfDNA datasets and can be readily transferred to diverse clinical settings for the expansion of liquid biopsy.

Computationally inferred cell-type specific epigenome-wide DNA methylation analysis unveils distinct methylation patterns among immune cells for HIV infection in three cohorts

BACKGROUND

Epigenome-wide association studies (EWAS) have identified CpG sites associated with HIV infection in blood cells in bulk, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. In this study, we aim to identify differentially methylated CpG sites for HIV infection in immune cell types: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes.

METHODS

Applying a computational deconvolution method, we performed a cell-type based EWAS for HIV infection in three independent cohorts (Ntotal = 1,382). DNA methylation in blood or in peripheral blood mononuclear cells (PBMCs) was profiled by an array-based method and then deconvoluted by Tensor Composition Analysis (TCA). The TCA-computed CpG methylation in each cell type was first benchmarked by bisulfite DNA methylation capture sequencing in a subset of the samples. Cell-type EWAS of HIV infection was performed in each cohort separately and a meta-EWAS was conducted followed by gene set enrichment analysis.

RESULTS

The meta-analysis unveiled a total of 2,021 cell-type unique significant CpG sites for five inferred cell types. Among these inferred cell-type unique CpG sites, the concordance rate in the three cohorts ranged from 96% to 100% in each cell type. Cell-type level meta-EWAS unveiled distinct patterns of HIV-associated differential CpG methylation, where 74% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of unique HIV-associated CpG sites (N = 1,624) compared to any other cell type. Genes harboring significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CD4+ T-cells: NLRC5, CX3CR1, B cells: IFI44L, NK cells: IL12R, monocytes: IRF7), and in oncogenesis (e.g. CD4+ T-cells: BCL family, PRDM16, monocytes: PRDM16, PDCD1LG2). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis that were enriched among interferon-α and -γ, TNF-α, inflammatory response, and apoptotic pathways.

CONCLUSION

Our findings uncovered computationally inferred cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding HIV pathogenesis.

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.

Circulating Leukocyte Subsets Before and After a Breast Cancer Diagnosis and Therapy

Importance

Changes in leukocyte composition often precede chronic disease onset. Patients with a history of breast cancer (hereinafter referred to as breast cancer survivors) are at increased risk for subsequent chronic diseases, but the long-term changes in peripheral leukocyte composition following a breast cancer diagnosis and treatment remain unknown.

Objective

To examine longitudinal changes in peripheral leukocyte composition in women who did and did not develop breast cancer and identify whether differences in breast cancer survivors were associated with specific treatments.

Design, Setting, and Participants

In this prospective cohort study, paired blood samples were collected from 2315 women enrolled in The Sister Study, a US-nationwide prospective cohort study of 50 884 women, at baseline (July 2003 to March 2009) and follow-up (October 2013 to March 2015) home visits, with a mean (SD) follow-up interval of 7.6 (1.4) years. By design, approximately half of the included women had been diagnosed and treated for breast cancer after enrollment and before the second blood draw. A total of 410 women were included in the present study, including 185 breast cancer survivors and 225 who remained free of breast cancer over a comparable follow-up period. Data were analyzed from April 21 to September 9, 2022.

Exposures

Breast cancer status and, among breast cancer survivors, cancer treatment type (chemotherapy, radiotherapy, endocrine therapy, or surgery).

Main Outcomes and Measures

Blood DNA methylation data were generated in 2019 using a genome-wide methylation screening tool and deconvolved to estimate percentages of 12 circulating leukocyte subsets.

Results

Of the 410 women included in the analysis, the mean (SD) age at enrollment was 56 (9) years. Compared with breast cancer-free women, breast cancer survivors had decreased percentages of circulating eosinophils (-0.45% [95% CI, -0.87% to -0.03%]; P = .03), total CD4+ helper T cells (-1.50% [95% CI, -2.56% to -0.44%]; P = .01), and memory B cells (-0.22% [95% CI, -0.34% to -0.09%]; P = .001) and increased percentages of circulating naive B cells (0.46% [95% CI, 0.17%-0.75%]; P = .002). In breast cancer survivor-only analyses, radiotherapy was associated with decreases in total CD4+ T cell levels, whereas chemotherapy was associated with increases in naive B cell levels. Surgery and endocrine therapy were not meaningfully associated with leukocyte changes.

Conclusions and Relevance

In this cohort study of 410 women, breast cancer survivors experienced lasting changes in peripheral leukocyte composition compared with women who remained free of breast cancer. These changes may be related to treatment with chemotherapy or radiotherapy and could influence future chronic disease risk.

DNA Methylation Signatures in Paired Placenta and Umbilical Cord Samples: Relationship with Maternal Pregestational Body Mass Index and Offspring Metabolic Outcomes

An epigenomic approach was used to study the impact of maternal pregestational body mass index (BMI) on the placenta and umbilical cord methylomes and their potential effect on the offspring's metabolic phenotype. DNA methylome was assessed in 24 paired placenta and umbilical cord samples. The differentially methylated CpGs associated with maternal pregestational BMI were identified and the metabolic pathways and the potentially related diseases affected by their annotated genes were determined. Two top differentially methylated CpGs were studied in 90 additional samples and the relationship with the offspring's metabolic phenotype was determined. The results showed that maternal pregestational BMI is associated with the methylation of genes involved in endocrine and developmental pathways with potential effects on type 2 diabetes and obesity. The methylation and expression of HADHA and SLC2A8 genes in placenta and umbilical cord were related to several metabolic parameters in the offspring at 6 years (weight SDS, height SDS, BMI SDS, Δ BW-BMI SDS, FM SDS, waist, SBP, TG, HOMA-IR, perirenal fat; all p < 0.05). Our data suggest that epigenetic analysis in placenta and umbilical cord may be useful for identifying individual vulnerability to later metabolic diseases.

Immunological shifts during early-stage Parkinson's disease identified with DNA methylation data on longitudinally collected blood samples

Parkinson's disease (PD) is the second most common neurodegenerative disease in the United States. Decades before motor symptoms manifest, non-motor symptoms such as hyposmia and rapid eye movement (REM) sleep behavior disorder are highly predictive of PD. Previous immune profiling studies have identified alterations to the proportions of immune cells in the blood of clinically defined PD patients. However, it remains unclear if these phenotypes manifest before the clinical diagnosis of PD. We utilized longitudinal DNA methylation (DNAm) microarray data from the Parkinson's Progression Marker's Initiative (PPMI) to perform immune profiling in clinically defined PD and prodromal PD patients (Prod). We identified previously reported changes in neutrophil, monocyte, and T cell numbers in PD patients. Additionally, we noted previously unrecognized decreases in the naive B cell compartment in the defined PD and Prod patient group. Over time, we observed the proportion of innate immune cells in PD blood increased, but the proportion of adaptive immune cells decreased. We identified decreases in T and B cell subsets associated with REM sleep disturbances and early cognitive decline. Lastly, we identified increases in B memory cells associated with both genetic (LRRK2 genotype) and infectious (cytomegalovirus seropositivity) risk factors of PD. Our analysis shows that the peripheral immune system is dynamic as the disease progresses. The study provides a platform to understand how and when peripheral immune alterations occur in PD and whether intervention at particular stages may be therapeutically advantageous.

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.

A noninvasive method for whole-genome skin methylome profiling

BACKGROUND

Ageing, disease and malignant transformation of the skin are associated with changes in DNA methylation. So far, mostly invasive methodologies such as biopsies have been applied in collecting DNA methylation signatures. Tape stripping offers a noninvasive option for skin diagnostics. It enables the easy but robust capture of biologic material in large numbers of participants without the need for specialized medical personnel.

OBJECTIVES

To design and validate a methodology for noninvasive skin sample collection using tape stripping for subsequent DNA -methylation analysis.

METHODS

A total of 175 participants were recruited and provided tape-stripping samples from a sun-exposed area; 92 provided matched tape-stripping samples from a sun-protected area, and an additional 5 provided matched skin-shave biopsies from the same area. Using -enzymatic conversion and whole-genome Illumina sequencing, we generated genome-wide DNA methylation profiles that were used to evaluate the feasibility of noninvasive data acquisition, to compare with established sampling approaches and to investigate biomarker identification for age and ultraviolet (UV) exposure.

RESULTS

We found that tape-stripping samples showed strong concordance in their global DNA methylation landscapes to those of conventional invasive biopsies. Moreover, we showed sample reproducibility and consistent global methylation profiles in skin tape-stripping samples collected from different areas of the body. Using matched samples from sun-protected and sun-exposed areas of the body we were able to validate the capacity of our method to capture the effects of environmental changes and ageing in a cohort covering various ages, ethnicities and skin types. We found DNA methylation changes on the skin resulting from UV exposure and identified significant age-related hypermethylation of CpG islands, with a pronounced peak effect at 50-55 years of age, including methylation changes in well-described markers of ageing.

CONCLUSIONS

These data demonstrate the feasibility of using tape stripping combined with whole-genome sequencing as a noninvasive approach to measuring DNA methylation changes in the skin. In addition, they outline a viable experimental framework for the use of skin tape stripping, particularly when it is performed in large cohorts of patients to identify biomarkers of skin ageing, UV damage and, possibly, to track treatment response to therapeutic interventions.

The application of epiphenotyping approaches to DNA methylation array studies of the human placenta

BACKGROUND

Genome-wide DNA methylation (DNAme) profiling of the placenta with Illumina Infinium Methylation bead arrays is often used to explore the connections between in utero exposures, placental pathology, and fetal development. However, many technical and biological factors can lead to signals of DNAme variation between samples and between cohorts, and understanding and accounting for these factors is essential to ensure meaningful and replicable data analysis. Recently, "epiphenotyping" approaches have been developed whereby DNAme data can be used to impute information about phenotypic variables such as gestational age, sex, cell composition, and ancestry. These epiphenotypes offer avenues to compare phenotypic data across cohorts, and to understand how phenotypic variables relate to DNAme variability. However, the relationships between placental epiphenotyping variables and other technical and biological variables, and their application to downstream epigenome analyses, have not been well studied.

RESULTS

Using DNAme data from 204 placentas across three cohorts, we applied the PlaNET R package to estimate epiphenotypes gestational age, ancestry, and cell composition in these samples. PlaNET ancestry estimates were highly correlated with independent polymorphic ancestry-informative markers, and epigenetic gestational age, on average, was estimated within 4 days of reported gestational age, underscoring the accuracy of these tools. Cell composition estimates varied both within and between cohorts, as well as over very long placental processing times. Interestingly, the ratio of cytotrophoblast to syncytiotrophoblast proportion decreased with increasing gestational age, and differed slightly by both maternal ethnicity (lower in white vs. non-white) and genetic ancestry (lower in higher probability European ancestry). The cohort of origin and cytotrophoblast proportion were the largest drivers of DNAme variation in this dataset, based on their associations with the first principal component.

CONCLUSIONS

This work confirms that cohort, array (technical) batch, cell type proportion, self-reported ethnicity, genetic ancestry, and biological sex are important variables to consider in any analyses of Illumina DNAme data. We further demonstrate the specific utility of epiphenotyping tools developed for use with placental DNAme data, and show that these variables (i) provide an independent check of clinically obtained data and (ii) provide a robust approach to compare variables across different datasets. Finally, we present a general framework for the processing and analysis of placental DNAme data, integrating the epiphenotype variables discussed here.

Deconvolution of immune cell composition and biological age of hepatocellular carcinoma using DNA methylation

Hepatocellular carcinoma (HCC) has been an approved indication for the administration of immunotherapy since 2017, but biomarkers that predict therapeutic response have remained limited. Understanding and characterizing the tumor immune microenvironment enables better classification of these tumors and may reveal biomarkers that predict immunotherapeutic efficacy. In this paper, we applied a cell-type deconvolution algorithm using DNA methylation array data to investigate the composition of the tumor microenvironment in HCC. Using publicly available and in-house datasets with a total cohort size of 57 patients, each with tumor and matched normal tissue samples, we identified key differences in immune cell composition. We found that NK cell abundance was significantly decreased in HCC tumors compared to adjacent normal tissue. We also applied DNA methylation "clocks" which estimate phenotypic aging and compared these findings to expression-based determinations of cellular senescence. Senescence and epigenetic aging were significantly increased in HCC tumors, and the degree of age acceleration and senescence was strongly associated with decreased NK cell abundance. In summary, we found that NK cell infiltration in the tumor microenvironment is significantly diminished, and that this loss of NK abundance is strongly associated with increased senescence and age-related phenotype. These findings point to key interactions between NK cells and the senescent tumor microenvironment and offer insights into the pathogenesis of HCC as well as potential biomarkers of therapeutic efficacy.

The role of DNA methylation in personalized medicine for immune-related diseases

Epigenetics functions as a bridge between host genetic & environmental factors, aiding in human health and diseases. Many immune-related diseases, including infectious and allergic diseases, have been linked to epigenetic mechanisms, particularly DNA methylation. In this review, we summarized an updated overview of DNA methylation and its importance in personalized medicine, and demonstrated that DNA methylation has excellent potential for disease prevention, diagnosis, and treatment in a personalized manner. The future implications and limitations of the DNA methylation study have also been well-discussed.

Epigenome-wide DNA methylation profiling in comparison between pathological and physiological hypertrophy of human cardiomyocytes

Background: Physiological and pathological stimuli result in distinct forms of cardiac hypertrophy, but the molecular regulation comparing the two, especially at the DNA methylation level, is not well understood. Methods: We conducted an in vitro study using human cardiomyocytes exposed to angiotensin II (AngII) and insulin-like growth factor 1 (IGF-1) to mimic pathologically and physiologically hypertrophic heart models, respectively. Whole genome DNA methylation patterns were profiled by the Infinium human MethylationEPIC platform with >850 K DNA methylation loci. Two external datasets were used for comparisons and qRT-PCR was performed for examining expression of associated genes of those identified DNA methylation loci. Results: We detected 194 loci that are significantly differentially methylated after AngII treatment, and 206 significant loci after IGF-1 treatment. Mapping the significant loci to genes, we identified 158 genes corresponding to AngII treatment and 175 genes to IGF-1 treatment. Using the gene-set enrichment analysis, the PI3K-Akt signaling pathway was identified to be significantly enriched for both AngII and IGF-1 treatment. The Hippo signaling pathway was enriched after IGF-1 treatment, but not for AngII treatment. CDK6 and RPTOR are components of the PI3K-Akt pathway but have different DNA methylation patterns in response to AngII and IGF-1. qRT-PCR confirmed the different gene expressions of CDK6 and PRTOR. Conclusion: Our study is pioneering in profiling epigenome DNA methylation changes in adult human cardiomyocytes under distinct stress conditions: pathological (AngII) and physiological (IGF-1). The identified DNA methylation loci, genes, and pathways might have the potential to distinguish between pathological and physiological cardiac hypertrophy.

Genome-Wide DNA Methylation Profiles in Whole-Blood and Buccal Samples-Cross-Sectional, Longitudinal, and across Platforms

The field of DNA methylation research is rapidly evolving, focusing on disease and phenotype changes over time using methylation measurements from diverse tissue sources and multiple array platforms. Consequently, identifying the extent of longitudinal, inter-tissue, and inter-platform variation in DNA methylation is crucial for future advancement. DNA methylation was measured in 375 individuals, with 197 of those having 2 blood sample measurements ~10 years apart. Whole-blood samples were measured on Illumina Infinium 450K and EPIC methylation arrays, and buccal samples from a subset of 58 participants were measured on EPIC array. The data were analyzed with the aims to examine the correlation between methylation levels in longitudinal blood samples in 197 individuals, examine the correlation between methylation levels in the blood and buccal samples in 58 individuals, and examine the correlation between blood methylation profiles assessed on the EPIC and 450K arrays in 83 individuals. We identified 136,833, 7674, and 96,891 CpGs significantly and strongly correlated (>0.50) longitudinally, across blood and buccal samples as well as array platforms, respectively. A total of 3674 of these CpGs were shared across all three sets. Analysis of these shared CpGs identified previously found associations with aging, ancestry, and 7016 mQTLs as well.

The mechanisms and factors that induce trained immunity in arthropods and mollusks

Besides dividing the organism's immune system into adaptive and innate immunity, it has long been thought that only adaptive immunity can establish immune memory. However, many studies have shown that innate immunity can also build immunological memory through epigenetic reprogramming and modifications to resist pathogens' reinfection, known as trained immunity. This paper reviews the role of mitochondrial metabolism and epigenetic modifications and describes the molecular foundation in the trained immunity of arthropods and mollusks. Mitochondrial metabolism and epigenetic modifications complement each other and play a key role in trained immunity.

DeCOOC Deconvoluted Hi-C Map Characterizes the Chromatin Architecture of Cells in Physiologically Distinctive Tissues

Deciphering variations in chromosome conformations based on bulk three-dimensional (3D) genomic data from heterogenous tissues is a key to understanding cell-type specific genome architecture and dynamics. Surprisingly, computational deconvolution methods for high-throughput chromosome conformation capture (Hi-C) data remain very rare in the literature. Here, a deep convolutional neural network (CNN), deconvolve bulk Hi-C data (deCOOC) that remarkably outperformed all the state-of-the-art tools in the deconvolution task is developed. Interestingly, it is noticed that the chromatin accessibility or the Hi-C contact frequency alone is insufficient to explain the power of deCOOC, suggesting the existence of a latent embedded layer of information pertaining to the cell type specific 3D genome architecture. By applying deCOOC to in-house-generated bulk Hi-C data from visceral and subcutaneous adipose tissues, it is found that the characteristic chromatin features of M2 cells in the two anatomical loci are distinctively bound to different physiological functionalities. Taken together, deCOOC is both a reliable Hi-C data deconvolution method and a powerful tool for functional extraction of 3D genome architecture.

A Review of the Epigenetic Clock: Emerging Biomarkers for Asthma and Allergic Disease

DNA methylation (DNAm) is a dynamic, age-dependent epigenetic modification that can be used to study interactions between genetic and environmental factors. Environmental exposures during critical periods of growth and development may alter DNAm patterns, leading to increased susceptibility to diseases such as asthma and allergies. One method to study the role of DNAm is the epigenetic clock-an algorithm that uses DNAm levels at select age-informative Cytosine-phosphate-Guanine (CpG) dinucleotides to predict epigenetic age (EA). The difference between EA and calendar age (CA) is termed epigenetic age acceleration (EAA) and reveals information about the biological capacity of an individual. Associations between EAA and disease susceptibility have been demonstrated for a variety of age-related conditions and, more recently, phenotypes such as asthma and allergic diseases, which often begin in childhood and progress throughout the lifespan. In this review, we explore different epigenetic clocks and how they have been applied, particularly as related to childhood asthma. We delve into how in utero and early life exposures (e.g., smoking, air pollution, maternal BMI) result in methylation changes. Furthermore, we explore the potential for EAA to be used as a biomarker for asthma and allergic diseases and identify areas for further study.

Implication of DNA methylation during lifestyle mediated weight loss

Over the past 50 years, the number of overweight/obese people increased significantly, making obesity a global public health challenge. Apart from rare monogenic forms, obesity is a multifactorial disease, most likely resulting from a concerted interaction of genetic, epigenetic and environmental factors. Although recent studies opened new avenues in elucidating the complex genetics behind obesity, the biological mechanisms contributing to individual's risk to become obese are not yet fully understood. Non-genetic factors such as eating behaviour or physical activity are strong contributing factors for the onset of obesity. These factors may interact with genetic predispositions most likely via epigenetic mechanisms. Epigenome-wide association studies or methylome-wide association studies are measuring DNA methylation at single CpGs across thousands of genes and capture associations to obesity phenotypes such as BMI. However, they only represent a snapshot in the complex biological network and cannot distinguish between causes and consequences. Intervention studies are therefore a suitable method to control for confounding factors and to avoid possible sources of bias. In particular, intervention studies documenting changes in obesity-associated epigenetic markers during lifestyle driven weight loss, make an important contribution to a better understanding of epigenetic reprogramming in obesity. To investigate the impact of lifestyle in obesity state specific DNA methylation, especially concerning the development of new strategies for prevention and individual therapy, we reviewed 19 most recent human intervention studies. In summary, this review highlights the huge potential of targeted interventions to alter disease-associated epigenetic patterns. However, there is an urgent need for further robust and larger studies to identify the specific DNA methylation biomarkers which influence obesity.

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.

Differences in white blood cell proportions between schizophrenia cases and controls are influenced by medication and variations in time of day

Cases with schizophrenia (SCZ) and healthy controls show differences in white blood cell (WBC) counts and blood inflammation markers. Here, we investigate whether time of blood draw and treatment with psychiatric medications are related to differences in estimated WBC proportions between SCZ cases and controls. DNA methylation data from whole blood was used to estimate proportions of six subtypes of WBCs in SCZ patients (n = 333) and healthy controls (n = 396). We tested the association of case-control status with estimated cell-type proportions and the neutrophil-to-lymphocyte ratio (NLR) in 4 models: with/without adjusting for time of blood draw, and then compared results from blood samples drawn during a 12-h (07:00-19:00) or 7-h (07:00-14:00) period. We also investigated WBC proportions in a subgroup of medication-free patients (n = 51). Neutrophil proportions were significantly higher in SCZ cases (mean=54.1%) vs. controls (mean=51.1%; p = <0.001), and CD8⁺T lymphocyte proportions were lower in SCZ cases (mean=12.1%) vs. controls (mean=13.2%; p = 0.001). The effect sizes in the 12-h sample (07:00-19:00) showed a significant difference between SCZ vs. controls for neutrophils, CD4⁺T, CD8⁺T, and B-cells, which remained significant after adjusting for time of blood draw. In the samples matched for time of blood draw during 07.00-14.00, we also observed an association with neutrophils, CD4⁺T, CD8⁺T, and B-cells that was unaffected by further adjustment for time of blood draw. In the medication-free patients, we observed differences that remained significant in neutrophils (p = 0.01) and CD4⁺T (p = 0.01) after adjusting for time of day. The association of SCZ with NLR was significant in all models (range: p < 0.001 to p = 0.03) in both medicated and unmedicated patients. In conclusion, controlling for pharmacological treatment and circadian cycling of WBC is necessary for unbiased estimates in case-control studies. Nevertheless, the association of WBC with SCZ remains, even after adjusting for the time of day.

Cell-free DNA-based liquid biopsies in neurology

This article reviews recent developments in the application of cell-free DNA-based liquid biopsies to neurological diseases. Over the past few decades, an explosion of interest in the use of accessible biofluids to identify and track molecular disease has revolutionized the fields of oncology, prenatal medicine and others. More recently, technological advances in signal detection have allowed for informative analysis of biofluids that are typically sparse in cells and other circulating components, such as CSF. In parallel, advancements in epigenetic profiling have allowed for novel applications of liquid biopsies to diseases without characteristic mutational profiles, including many degenerative, autoimmune, inflammatory, ischaemic and infectious disorders. These events have paved the way for a wide array of neurological conditions to benefit from enhanced diagnostic, prognostic, and treatment abilities through the use of liquid biomarkers: a 'liquid biopsy' approach. This review includes an overview of types of liquid biopsy targets with a focus on circulating cell-free DNA, methods used to identify and probe potential liquid biomarkers, and recent applications of such biomarkers to a variety of complex neurological conditions including CNS tumours, stroke, traumatic brain injury, Alzheimer's disease, epilepsy, multiple sclerosis and neuroinfectious disease. Finally, the challenges of translating liquid biopsies to use in clinical neurology settings-and the opportunities for improvement in disease management that such translation may provide-are discussed.

Genes implicated by a methylome-wide schizophrenia study in neonatal blood show differential expression in adult brain samples

Schizophrenia is a disabling disorder involving genetic predisposition in combination with environmental influences that likely act via dynamic alterations of the epigenome and the transcriptome but its detailed pathophysiology is largely unknown. We performed cell-type specific methylome-wide association study of neonatal blood (N = 333) from individuals who later in life developed schizophrenia and controls. Suggestively significant associations (P < 1.0 × 10-6) were detected in all cell-types and in whole blood with methylome-wide significant associations in monocytes (P = 2.85 × 10-9-4.87 × 10-9), natural killer cells (P = 1.72 × 10-9-7.82 × 10-9) and B cells (P = 3.8 × 10-9). Validation of methylation findings in post-mortem brains (N = 596) from independent schizophrenia cases and controls showed significant enrichment of transcriptional differences (enrichment ratio = 1.98-3.23, P = 2.3 × 10-3-1.0 × 10-5), with specific highly significant differential expression for, for example, BDNF (t = -6.11, P = 1.90 × 10-9). In addition, expression difference in brain significantly predicted schizophrenia (multiple correlation = 0.15-0.22, P = 3.6 × 10-4-4.5 × 10-8). In summary, using a unique design combining pre-disease onset (neonatal) blood methylomic data and post-disease onset (post-mortem) brain transcriptional data, we have identified genes of likely functional relevance that are associated with schizophrenia susceptibility, rather than confounding disease associated artifacts. The identified loci may be of clinical value as a methylation-based biomarker for early detection of increased schizophrenia susceptibility.

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.

Calculating detection limits and uncertainty of reference-based deconvolution of whole-blood DNA methylation data

DNA methylation (DNAm)-based cell mixture deconvolution (CMD) has become a quintessential part of epigenome-wide association studies where DNAm is profiled in heterogeneous tissue types. Despite being introduced over a decade ago, detection limits, which represent the smallest fraction of a cell type in a mixed biospecimen that can be reliably detected, have yet to be determined in the context of DNAm-based CMD. Moreover, there has been little attention given to approaches for quantifying the uncertainty associated with DNAm-based CMD. Here, analytical frameworks for determining both cell-specific limits of detection and quantification of uncertainty associated with DNAm-based CMD are described. This work may contribute to improved rigor, reproducibility and replicability of epigenome-wide association studies involving CMD.

FOXP3 TSDR Measurement Could Assist Variant Classification and Diagnosis of IPEX Syndrome

Pathogenic FOXP3 variants cause immune dysregulation polyendocrinopathy enteropathy X-linked (IPEX) syndrome, a progressive autoimmune disease resulting from disruption of the regulatory T cell (Treg) compartment. Assigning pathogenicity to novel variants in FOXP3 is challenging due to the heterogeneous phenotype and variable immunological abnormalities. The number of cells with demethylation at the Treg cell-specific demethylated region (TSDR) is an independent biomarker of IPEX. We aimed to investigate if diagnosing IPEX at presentation with isolated diabetes could allow for effective monitoring of disease progression and assess whether TSDR analysis can aid FOXP3 variant classification and predict disease course. We describe a large genetically diagnosed IPEX cohort (n = 65) and 13 individuals with other monogenic autoimmunity subtypes in whom we quantified the proportion of cells with FOXP3 TSDR demethylation, normalized to the number with CD4 demethylation (%TSDR/CD4) and compare them to 29 unaffected controls. IPEX patients presenting with isolated diabetes (50/65, 77%) often later developed enteropathy (20/50, 40%) with a median interval of 23.5 weeks. %TSDR/CD4 was a good discriminator of IPEX vs. unaffected controls (ROC-AUC 0.81, median 13.6% vs. 8.5%, p < 0.0001) with higher levels of demethylation associated with more severe disease. Patients with other monogenic autoimmunity had a similar %TSDR/CD4 to controls (median 8.7%, p = 1.0). Identifying increased %TSDR/CD4 in patients with novel FOXP3 mutations presenting with isolated diabetes facilitates diagnosis and could offer an opportunity to monitor patients and begin immune modulatory treatment before onset of severe enteropathy.

Brain-derived neurotrophic factor (BDNF) epigenomic modifications and brain-related phenotypes in humans: A systematic review

Epigenomic modifications of the brain-derived neurotrophic factor (BDNF) gene have been postulated to underlie the pathogenesis of neurodevelopmental, psychiatric, and neurological conditions. This systematic review summarizes current evidence investigating the association of BDNF epigenomic modifications (DNA methylation, non-coding RNA, histone modifications) with brain-related phenotypes in humans. A novel contribution is our creation of an open access web-based application, the BDNF DNA Methylation Map, to interactively visualize specific positions of CpG sites investigated across all studies for which relevant data were available. Our literature search of four databases through September 27, 2021 returned 1701 articles, of which 153 met inclusion criteria. Our review revealed exceptional heterogeneity in methodological approaches, hindering the identification of clear patterns of robust and/or replicated results. We summarize key findings and provide recommendations for future epigenomic research. The existing literature appears to remain in its infancy and requires additional rigorous research to fulfill its potential to explain BDNF-linked risk for brain-related conditions and improve our understanding of the molecular mechanisms underlying their pathogenesis.

Cell-type specific EWAS identifies genes involved in HIV pathogenesis and oncogenesis among people with HIV infection

Epigenome-wide association studies (EWAS) of heterogenous blood cells have identified CpG sites associated with chronic HIV infection, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. Applying a computational deconvolution method validated by capture bisulfite DNA methylation sequencing, we conducted a cell type-based EWAS and identified differentially methylated CpG sites specific for chronic HIV infection among five immune cell types in blood: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes in two independent cohorts (N total =1,134). Differentially methylated CpG sites for HIV-infection were highly concordant between the two cohorts. Cell-type level meta-EWAS revealed distinct patterns of HIV-associated differential CpG methylation, where 67% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of HIV-associated CpG sites (N=1,472) compared to any other cell type. Genes harboring statistically significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CX3CR1 in CD4+ T-cells, CCR7 in B cells, IL12R in NK cells, LCK in monocytes). More importantly, HIV-associated CpG sites were overrepresented for hallmark genes involved in cancer pathology ( FDR <0.05) (e.g. BCL family, PRDM16, PDCD1LGD, ESR1, DNMT3A, NOTCH2 ). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis such as Kras-signaling, interferon-α and -γ, TNF-α, inflammatory, and apoptotic pathways. Our findings are novel, uncovering cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding pathogen-induced epigenetic oncogenicity, specifically on HIV and its comorbidity with cancers.

Mitochondrial DNA copy number dynamics and associations with the prenatal environment from birth through adolescence in a population of Dominican and African American children

Mitochondrial DNA copy number (mtDNAcn) dynamics throughout childhood are poorly understood. We profiled mtDNAcn from birth through adolescence and evaluated how the prenatal environment influences mtDNAcn across childhood. Data were collected from children from New York City followed through 18 years. Using duplexed qRT-PCR, we quantified mtDNAcn relative to nuclear DNA in blood collected from the umbilical cord (n = 450), children aged 5-7 (n = 510), and adolescents aged 15-18 (n = 278). We examined mtDNAcn across childhood with linear mixed-effects models (LMM). Relative mtDNAcn was lowest at birth (mean ± SD: 0.67 ± 0.35) and increased in childhood (1.24 ± 0.50) then slightly declined in adolescence (1.13 ± 0.44). We observed no differences in mtDNAcn by sex or race/ethnicity. mtDNAcn was positively associated with prenatal environmental tobacco smoke exposure (0.077 [ 0.01, 0.14] change in relative mtDNAcn) but negatively associated with maternal completion of high school (-0.066 [-0.13, 0.00]), with the receipt of public assistance at birth (-0.074 [-0.14, -0.01]), and when mother born outside the U.S (-0.061 [-0.13, 0.003]). Infant birth outcomes were not associated with mtDNAcn. MtDNAcn levels were dynamic through childhood and associated with some prenatal factors, underscoring the need for the investigation of longitudinal mtDNAcn for human health research.

Nucleated red blood cells explain most of the association between DNA methylation and gestational age

Determining if specific cell type(s) are responsible for an association between DNA methylation (DNAm) and a given phenotype is important for understanding the biological mechanisms underlying the association. Our EWAS of gestational age (GA) in 953 newborns from the Norwegian MoBa study identified 13,660 CpGs significantly associated with GA (p_Bonferroni<0.05) after adjustment for cell type composition. When the CellDMC algorithm was applied to explore cell-type specific effects, 2,330 CpGs were significantly associated with GA, mostly in nucleated red blood cells [nRBCs; n = 2,030 (87%)]. Similar patterns were found in another dataset based on a different array and when applying an alternative algorithm to CellDMC called Tensor Composition Analysis (TCA). Our findings point to nRBCs as the main cell type driving the DNAm-GA association, implicating an epigenetic signature of erythropoiesis as a likely mechanism. They also explain the poor correlation observed between epigenetic age clocks for newborns and those for adults.

Epigenetic clocks derived from western samples differentially reflect Taiwanese health outcomes

Introduction: Several epigenetic clocks have been developed, with five measures of epigenetic age acceleration (EAA) especially receiving extensive investigations: HannumEAA, IEAA, PhenoEAA, GrimEAA, and DunedinPACE. These epigenetic clocks were mainly developed by individuals of European or Hispanic ancestry. It remains unclear whether they can reflect disease morbidity and physiological conditions in Asian populations. Methods: I here investigated five measures of EAA of 2,474 Taiwan Biobank participants with DNA methylation data. Using logistic regressions, I sequentially regressed various health outcomes on each of the five measures of EAA while adjusting for chronological age, sex, body mass index, the number of smoking pack-years, drinking status, regular exercise, educational attainment, and six cell-type proportions. Results: Except for IEAA, all measures of EAA reflected the obesity of Taiwanese (p < 4.0E-4). Diabetes was reflected by DunedinPACE (p = 5.4E-6) and GrimEAA (p = 5.8E-5). Moreover, DunedinPACE was associated with dyslipidemia, including hypertriglyceridemia (p = 1.1E-5), low high-density lipoprotein cholesterol (HDL-C) (p = 4.0E-5), and high triglyceride to HDL-C ratio (p = 1.6E-7). Discussion: This is one of the first studies to show that epigenetic clocks (developed by individuals of European or Hispanic ancestry) can reflect Taiwanese physiological conditions. DunedinPACE was associated with more Taiwanese health outcomes than the other four measures of EAA.

Epigenome-wide DNA methylation in externalizing behaviours: A review and combined analysis

DNA methylation (DNAm) is one of the most frequently studied epigenetic mechanisms facilitating the interplay of genomic and environmental factors, which can contribute to externalizing behaviours and related psychiatric disorders. Previous epigenome-wide association studies (EWAS) for externalizing behaviours have been limited in sample size, and, therefore, candidate genes and biomarkers with robust evidence are still lacking. We 1) performed a systematic literature review of EWAS of attention-deficit/hyperactivity disorder (ADHD)- and aggression-related behaviours conducted in peripheral tissue and cord blood and 2) combined the most strongly associated DNAm sites observed in individual studies (p < 10^-3) to identify candidate genes and biological systems for ADHD and aggressive behaviours. We observed enrichment for neuronal processes and neuronal cell marker genes for ADHD. Astrocyte and granulocytes cell markers among genes annotated to DNAm sites were relevant for both ADHD and aggression-related behaviours. Only 1 % of the most significant epigenetic findings for ADHD/ADHD symptoms were likely to be directly explained by genetic factors involved in ADHD. Finally, we discuss how the field would greatly benefit from larger sample sizes and harmonization of assessment instruments.

Mediation analysis of multiple mediators with incomplete omics data

There is an increasing interest in using multiple types of omics features (e.g., DNA sequences, RNA expressions, methylation, protein expressions, and metabolic profiles) to study how the relationships between phenotypes and genotypes may be mediated by other omics markers. Genotypes and phenotypes are typically available for all subjects in genetic studies, but typically, some omics data will be missing for some subjects, due to limitations such as cost and sample quality. In this article, we propose a powerful approach for mediation analysis that accommodates missing data among multiple mediators and allows for various interaction effects. We formulate the relationships among genetic variants, other omics measurements, and phenotypes through linear regression models. We derive the joint likelihood for models with two mediators, accounting for arbitrary patterns of missing values. Utilizing computationally efficient and stable algorithms, we conduct maximum likelihood estimation. Our methods produce unbiased and statistically efficient estimators. We demonstrate the usefulness of our methods through simulation studies and an application to the Metabolic Syndrome in Men study.

A DNA methylation atlas of normal human cell types

DNA methylation is a fundamental epigenetic mark that governs gene expression and chromatin organization, thus providing a window into cellular identity and developmental processes1. Current datasets typically include only a fraction of methylation sites and are often based either on cell lines that underwent massive changes in culture or on tissues containing unspecified mixtures of cells2-5. Here we describe a human methylome atlas, based on deep whole-genome bisulfite sequencing, allowing fragment-level analysis across thousands of unique markers for 39 cell types sorted from 205 healthy tissue samples. Replicates of the same cell type are more than 99.5% identical, demonstrating the robustness of cell identity programmes to environmental perturbation. Unsupervised clustering of the atlas recapitulates key elements of tissue ontogeny and identifies methylation patterns retained since embryonic development. Loci uniquely unmethylated in an individual cell type often reside in transcriptional enhancers and contain DNA binding sites for tissue-specific transcriptional regulators. Uniquely hypermethylated loci are rare and are enriched for CpG islands, Polycomb targets and CTCF binding sites, suggesting a new role in shaping cell-type-specific chromatin looping. The atlas provides an essential resource for study of gene regulation and disease-associated genetic variants, and a wealth of potential tissue-specific biomarkers for use in liquid biopsies.

Methylation differences in Alzheimer's disease neuropathologic change in the aged human brain

Alzheimer's disease (AD) is the most common cause of dementia with advancing age as its strongest risk factor. AD neuropathologic change (ADNC) is known to be associated with numerous DNA methylation changes in the human brain, but the oldest old (> 90 years) have so far been underrepresented in epigenetic studies of ADNC. Our study participants were individuals aged over 90 years (n = 47) from The 90+ Study. We analyzed DNA methylation from bulk samples in eight precisely dissected regions of the human brain: middle frontal gyrus, cingulate gyrus, entorhinal cortex, dentate gyrus, CA1, substantia nigra, locus coeruleus and cerebellar cortex. We deconvolved our bulk data into cell-type-specific (CTS) signals using computational methods. CTS methylation differences were analyzed across different levels of ADNC. The highest amount of ADNC related methylation differences was found in the dentate gyrus, a region that has so far been underrepresented in large scale multi-omic studies. In neurons of the dentate gyrus, DNA methylation significantly differed with increased burden of amyloid beta (Aβ) plaques at 5897 promoter regions of protein-coding genes. Amongst these, higher Aβ plaque burden was associated with promoter hypomethylation of the Presenilin enhancer 2 (PEN-2) gene, one of the rate limiting genes in the formation of gamma-secretase, a multicomponent complex that is responsible in part for the endoproteolytic cleavage of amyloid precursor protein into Aβ peptides. In addition to novel ADNC related DNA methylation changes, we present the most detailed array-based methylation survey of the old aged human brain to date. Our open-sourced dataset can serve as a brain region reference panel for future studies and help advance research in aging and neurodegenerative diseases.

A systematic assessment of cell type deconvolution algorithms for DNA methylation data

We performed systematic assessment of computational deconvolution methods that play an important role in the estimation of cell type proportions from bulk methylation data. The proposed framework methylDeConv (available as an R package) integrates several deconvolution methods for methylation profiles (Illumina HumanMethylation450 and MethylationEPIC arrays) and offers different cell-type-specific CpG selection to construct the extended reference library which incorporates the main immune cell subsets, epithelial cells and cell-free DNAs. We compared the performance of different deconvolution algorithms via simulations and benchmark datasets and further investigated the associations of the estimated cell type proportions to cancer therapy in breast cancer and subtypes in melanoma methylation case studies. Our results indicated that the deconvolution based on the extended reference library is critical to obtain accurate estimates of cell proportions in non-blood tissues.

HiTIMED: hierarchical tumor immune microenvironment epigenetic deconvolution for accurate cell type resolution in the tumor microenvironment using tumor-type-specific DNA methylation data

Background

Cellular compositions of solid tumor microenvironments are heterogeneous, varying across patients and tumor types. High-resolution profiling of the tumor microenvironment cell composition is crucial to understanding its biological and clinical implications. Previously, tumor microenvironment gene expression and DNA methylation-based deconvolution approaches have been shown to deconvolve major cell types. However, existing methods lack accuracy and specificity to tumor type and include limited identification of individual cell types.

Results

We employed a novel tumor-type-specific hierarchical model using DNA methylation data to deconvolve the tumor microenvironment with high resolution, accuracy, and specificity. The deconvolution algorithm is named HiTIMED. Seventeen cell types from three major tumor microenvironment components can be profiled (tumor, immune, angiogenic) by HiTIMED, and it provides tumor-type-specific models for twenty carcinoma types. We demonstrate the prognostic significance of cell types that other tumor microenvironment deconvolution methods do not capture.

Conclusion

We developed HiTIMED, a DNA methylation-based algorithm, to estimate cell proportions in the tumor microenvironment with high resolution and accuracy. HiTIMED deconvolution is amenable to archival biospecimens providing high-resolution profiles enabling to study of clinical and biological implications of variation and composition of the tumor microenvironment.

Ecoimmunology in the field: Measuring multiple dimensions of immune function with minimally invasive, field-adapted techniques

OBJECTIVE

Immune function is multifaceted and characterizations based on single biomarkers may be uninformative or misleading, particularly when considered across ecological contexts. However, measuring the many facets of immunity in the field can be challenging, since many measures cannot be obtained on-site, necessitating sample preservation and transport. Here we assess state-of-the-art methods for measuring immunity, focusing on measures that require a minimal blood sample obtained from a finger prick, which can be: (1) dried on filter paper, (2) frozen in liquid nitrogen, or (3) stabilized with chemical reagents.

RESULTS

We review immune measures that can be obtained from point-of-care devices or from immunoassays of dried blood spots (DBSs), field methods for flow cytometry, the use of RNA or DNA sequencing and quantification, and the application of immune activation assays under field conditions.

CONCLUSIONS

Stable protein products, such as immunoglobulins and C-reactive protein are reliably measured in DBSs. Because less stable proteins, such as cytokines, may be problematic to measure even in fresh blood, mRNA from stabilized blood may provide a cleaner measure of cytokine and broader immune-related gene expression. Gene methylation assays or mRNA sequencing also allow for the quantification of many other parameters, including the inference of leukocyte subsets, though with less accuracy than with flow cytometry. Combining these techniques provides an improvement over single-marker studies, allowing for a more nuanced understanding of how social and ecological variables are linked to immune measures and disease risk in diverse populations and settings.

A detection panel of novel methylated DNA markers for malignant pleural effusion

Background

Cytology remains the gold standard for the detection of malignant cells in pleural effusion. However, its sensitivity is limited. The aim of this study was to establish a novel panel of cancer-specific methylated genes for the differential diagnosis of malignant pleural effusion (MPE).

Methods

A cohort of 100 cancer patients (68 lung cancer, 32 other malignant tumors) and 48 patients with benign disease presenting with pleural effusion was prospectively enrolled. Pleural effusion was evaluated by means of cytopathological investigation and DNA methylation of SHOX2, RASSF1A, SEPTIN9 and HOXA9 in the cellular fraction. DNA methylation in bisulfite-converted DNA was determined using quantitative methylation-specific real-time PCR (MS-PCR). Cytopathological and DNA methylation results were evaluated with regard to the final clinical diagnosis.

Results

The LungMe® SHOX2 and RASSF1A Assay (Tellgen Corporation, China) has been reported to be highly sensitive and specific for lung cancer using bronchial aspirates. As expected, LungMe® detected metastases of lung cancer (sensitivity: 76.5%) as well as metastases of other malignant tumors (sensitivity: 68.8%). OncoMe, a novel combination of SHOX2, RASSF1A, SEPTIN9 and HOXA9 methylation, led to an additional 11% increase in the detection rate of MPE, resulting in a sensitivity of 85% and a specificity of 96%. Overall, OncoMe showed a higher positive detection rate in SCLC (100%), LUAC (87%), OC (100%), BC (92.9%), GC (80.0%), and MESO (80%) than in LUSC (50%). Cytopathological analyses only detected 23 positive samples, which were all positively measured by both LungMe® and OncoMe.

Conclusion

OncoMe has potential for use as a biomarker for the detection of MPE, even not limited to lung cancer.

Maternal obesity alters methylation level of cytosine in CpG island for epigenetic inheritance in fetal umbilical cord blood

BACKGROUND

Over the past few decades, global maternal obesity prevalence has rapidly increased. This condition may induce long-lasting pathophysiological effects on either fetal or infant health that could be attributable to unknown unique changes in the umbilical blood composition.

METHODS

A total of 34 overweight/obese and 32 normal-weight pregnant women were recruited. Fifteen umbilical blood samples including 8 overweight/obese subjects and 7 normal weight women were sequenced using Targeted Bisulfite Sequencing technology to detect the average methylation level of cytosine and identify the differentially methylated region (DMR). GO and KEGG analyses were then employed to perform pathway enrichment analysis of DMR-related genes and promoters. Moreover, the mRNA levels of methylation-related genes histone deacetylases (HDACs) and DNA methyltransferases (DNMTs) were characterized in the samples obtained from these two groups.

RESULTS

Average methylated cytosine levels in both the CpG islands (CGI) and promoter significantly decreased in overweight/obese groups. A total of 1669 DMRs exhibited differences in their DNA methylation status between the overweight/obese and control groups. GO and KEGG analyses revealed that DMR-related genes and promoters were enriched in the metabolism, cancer and cardiomyopathy signaling pathways. Furthermore, the HDACs and DNMTs mRNA levels trended to decline in overweight/obese groups.

CONCLUSIONS

Decreased methylated cytosine levels in overweight/obese women induce the gene expression activity at a higher level than in the control group. DMRs between these two groups in the fetal blood may contribute to the changes in gene transcription that underlie the increased risk of metabolic disorders, cancers and cardiomyopathy in their offspring.

Metabolomic predictors of phenotypic traits can replace and complement measured clinical variables in population-scale expression profiling studies

Population-scale expression profiling studies can provide valuable insights into biological and disease-underlying mechanisms. The availability of phenotypic traits is essential for studying clinical effects. Therefore, missing, incomplete, or inaccurate phenotypic information can make analyses challenging and prevent RNA-seq or other omics data to be reused. A possible solution are predictors that infer clinical or behavioral phenotypic traits from molecular data. While such predictors have been developed based on different omics data types and are being applied in various studies, metabolomics-based surrogates are less commonly used than predictors based on DNA methylation profiles.In this study, we inferred 17 traits, including diabetes status and exposure to lipid medication, using previously trained metabolomic predictors. We evaluated whether these metabolomic surrogates can be used as an alternative to reported information for studying the respective phenotypes using expression profiling data of four population cohorts. For the majority of the 17 traits, the metabolomic surrogates performed similarly to the reported phenotypes in terms of effect sizes, number of significant associations, replication rates, and significantly enriched pathways.The application of metabolomics-derived surrogate outcomes opens new possibilities for reuse of multi-omics data sets. In studies where availability of clinical metadata is limited, missing or incomplete information can be complemented by these surrogates, thereby increasing the size of available data sets. Additionally, the availability of such surrogates could be used to correct for potential biological confounding. In the future, it would be interesting to further investigate the use of molecular predictors across different omics types and cohorts.

Tracing the Origin of Cell-Free DNA Molecules through Tissue-Specific Epigenetic Signatures

All cell and tissue types constantly release DNA fragments into human body fluids by various mechanisms including programmed cell death, accidental cell degradation and active extrusion. Particularly, cell-free DNA (cfDNA) in plasma or serum has been utilized for minimally invasive molecular diagnostics. Disease onset or pathological conditions that lead to increased cell death alter the contribution of different tissues to the total pool of cfDNA. Because cfDNA molecules retain cell-type specific epigenetic features, it is possible to infer tissue-of-origin from epigenetic characteristics. Recent research efforts demonstrated that analysis of, e.g., methylation patterns, nucleosome occupancy, and fragmentomics determined the cell- or tissue-of-origin of individual cfDNA molecules. This novel tissue-of origin-analysis enables to estimate the contributions of different tissues to the total cfDNA pool in body fluids and find tissues with increased cell death (pathologic condition), expanding the portfolio of liquid biopsies towards a wide range of pathologies and early diagnosis. In this review, we summarize the currently available tissue-of-origin approaches and point out the next steps towards clinical implementation.