Cell-type-specific resolution epigenetics without the need for cell sorting or single-cell biology

Methods in DNA methylation array dataset analysis: A review

Understanding the intricate relationships between gene expression levels and epigenetic modifications in a genome is crucial to comprehending the pathogenic mechanisms of many diseases. With the advancement of DNA Methylome Profiling techniques, the emphasis on identifying Differentially Methylated Regions (DMRs/DMGs) has become crucial for biomarker discovery, offering new insights into the etiology of illnesses. This review surveys the current state of computational tools/algorithms for the analysis of microarray-based DNA methylation profiling datasets, focusing on key concepts underlying the diagnostic/prognostic CpG site extraction. It addresses methodological frameworks, algorithms, and pipelines employed by various authors, serving as a roadmap to address challenges and understand changing trends in the methodologies for analyzing array-based DNA methylation profiling datasets derived from diseased genomes. Additionally, it highlights the importance of integrating gene expression and methylation datasets for accurate biomarker identification, explores prognostic prediction models, and discusses molecular subtyping for disease classification. The review also emphasizes the contributions of machine learning, neural networks, and data mining to enhance diagnostic workflow development, thereby improving accuracy, precision, and robustness.

cypress: an R/Bioconductor package for cell-type-specific differential expression analysis power assessment

SUMMARY

Recent methodology advances in computational signal deconvolution have enabled bulk transcriptome data analysis at a finer cell-type level. Through deconvolution, identifying cell-type-specific differentially expressed (csDE) genes is drawing increasing attention in clinical applications. However, researchers still face a number of difficulties in adopting csDE genes detection methods in practice, especially in their experimental design. Here we present cypress, the first experimental design and statistical power analysis tool in csDE genes identification. This tool can reliably model purified cell-type-specific (CTS) profiles, cell-type compositions, biological and technical variations, offering a high-fidelity simulator for bulk RNA-seq convolution and deconvolution. cypress conducts simulation and evaluates the impact of multiple influencing factors, by various statistical metrics, to help researchers optimize experimental design and conduct power analysis.

AVAILABILITY AND IMPLEMENTATION

cypress is an open-source R/Bioconductor package at https://bioconductor.org/packages/cypress/.

Inhibitory mechanisms in the prefrontal-cortex differentially mediate Putamen activity during valence-based learning

Learning from appetitive and aversive stimuli is important for survival. It involves interactions between the prefrontal cortex and subcortical structures, with inhibition playing a crucial role. However, direct evidence for this in humans is limited. Here, we overcome the difficulty of measuring inhibition in the human brain and find that GABA, the main inhibitory neurotransmitter, affects how the dACC interacts with subcortical structures during appetitive and aversive learning differently. We used 7T magnetic resonance spectroscopy (MRS) to track GABA levels in the dACC alongside whole-brain fMRI scans while participants engaged in appetitive and aversive learning tasks. During appetitive learning, dACC GABA levels were negatively correlated with learning performance and BOLD activity measured from the dACC and the Putamen. While under aversive learning, dACC GABA concentration negatively correlated with the functional connectivity between the dACC and the Putamen. Our results show that inhibition in the dACC mediates appetitive and aversive learning in humans through distinct mechanisms.

Chromatin- and nucleosome-associated features in liquid biopsy: implications for cancer biomarker discovery

Cell-free DNA (cfDNA) from the bloodstream has been studied for cancer biomarker discovery, and chromatin-derived epigenetic features have come into the spotlight for their potential to expand clinical applications. Methylation, fragmentation, and nucleosome positioning patterns of cfDNA have previously been shown to reveal epigenomic and inferred transcriptomic information. More recently, histone modifications have emerged as a tool to further identify tumor-specific chromatin variants in plasma. A number of sequencing methods have been developed to analyze these epigenetic markers, offering new insights into tumor biology. Features within cfDNA allow for cancer detection, subtype and tissue of origin classification, and inference of gene expression. These methods provide a window into the complexity of cancer and the dynamic nature of its progression. In this review, we highlight the array of epigenetic features in cfDNA that can be extracted from chromatin- and nucleosome-associated organization and outline potential use cases in cancer management.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods estimate cell-type proportions and gene expressions in bulk tissue samples, yet their performance and biological applications remain unexplored, particularly in human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk tissue RNA sequencing (RNA-seq), single-cell/nuclei (sc/sn) RNA-seq, and immunohistochemistry. A total of 1,130,767 nuclei per cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expressions. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk tissue or single-cell eQTLs did alone. Differential gene expressions associated with Alzheimer's disease, schizophrenia, and brain development were also examined using the deconvoluted data. Our findings, which were replicated in bulk tissue and single-cell data, provided insights into the biological applications of deconvoluted data in multiple brain disorders.

Cell-type specific inference from bulk RNA-sequencing data by integrating single cell reference profiles via EPIC-unmix

Cell type specific (CTS) analysis is essential to reveal biological insights obscured in bulk tissue data. However, single-cell (sc) or single-nuclei (sn) resolution data are still cost-prohibitive for large-scale samples. Thus, computational methods to perform deconvolution from bulk tissue data are highly valuable. We here present EPIC-unmix, a novel two-step empirical Bayesian method integrating reference sc/sn RNA-seq data and bulk RNA-seq data from target samples to enhance the accuracy of CTS inference. We demonstrate through comprehensive simulations across three tissues that EPIC-unmix achieved 4.6% - 109.8% higher accuracy compared to alternative methods. By applying EPIC-unmix to human bulk brain RNA-seq data from the ROSMAP and MSBB cohorts, we identified multiple genes differentially expressed between Alzheimer's disease (AD) cases versus controls in a CTS manner, including 57.4% novel genes not identified using similar sample size sc/snRNA-seq data, indicating the power of our in-silico approach. Among the 6-69% overlapping, 83%-100% are in consistent direction with those from sc/snRNA-seq data, supporting the reliability of our findings. EPIC-unmix inferred CTS expression profiles similarly empowers CTS eQTL analysis. Among the novel eQTLs, we highlight a microglia eQTL for AD risk gene AP3B2, obscured in bulk and missed by sc/snRNA-seq based eQTL analysis. The variant resides in a microglia-specific cCRE, forming chromatin loop with AP3B2 promoter region in microglia. Taken together, we believe EPIC-unmix will be a valuable tool to enable more powerful CTS analysis.

Identification of cell-specific epigenetic patterns associated with chondroitin sulfate treatment response in an endemic arthritis, Kashin-Beck disease

Aims

To assess the alterations in cell-specific DNA methylation associated with chondroitin sulphate response using peripheral blood collected from Kashin-Beck disease (KBD) patients before initiation of chondroitin sulphate treatment.

Methods

Peripheral blood samples were collected from KBD patients at baseline of chondroitin sulphate treatment. Methylation profiles were generated using reduced representation bisulphite sequencing (RRBS) from peripheral blood. Differentially methylated regions (DMRs) were identified using MethylKit, while DMR-related genes were defined as those annotated to the gene body or 2.2-kilobase upstream regions of DMRs. Selected DMR-related genes were further validated by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) to assess expression levels. Tensor composition analysis was performed to identify cell-specific differential DNA methylation from bulk tissue.

Results

This study revealed 21,060 hypermethylated and 44,472 hypomethylated DMRs, and 13,194 hypermethylated and 22,448 hypomethylated CpG islands for differential global methylation for chondroitin sulphate treatment response. A total of 12,666 DMR-related genes containing DMRs were identified in their promoter regions, such as CHL1 (false discovery rate (FDR) = 2.11 × 10-11), RIC8A (FDR = 7.05 × 10-4), and SOX12 (FDR = 1.43 × 10-3). Additionally, RIC8A and CHL1 were hypermethylated in responders, while SOX12 was hypomethylated in responders, all showing decreased gene expression. The patterns of cell-specific differential global methylation associated with chondroitin sulphate response were observed. Specifically, we found that DMRs located in TESPA1 and ATP11A exhibited differential DNA methylation between responders and non-responders in granulocytes, monocytes, and B cells.

Conclusion

Our study identified cell-specific changes in DNA methylation associated with chondroitin sulphate response in KBD patients.

regionalpcs: improved discovery of DNA methylation associations with complex traits

We have developed the regional principal components (rPCs) method, a novel approach for summarizing gene-level methylation. rPCs address the challenge of deciphering complex epigenetic mechanisms in diseases like Alzheimer's disease (AD). In contrast to traditional averaging, rPCs leverage principal components analysis to capture complex methylation patterns across gene regions. Our method demonstrated a 54% improvement in sensitivity over averaging in simulations, offering a robust framework for identifying subtle epigenetic variations. Applying rPCs to the AD brain methylation data in ROSMAP, combined with cell type deconvolution, we uncovered 838 differentially methylated genes associated with neuritic plaque burden-significantly outperforming conventional methods. Integrating methylation quantitative trait loci (meQTL) with genome-wide association studies (GWAS) identified 17 genes with potential causal roles in AD, including MS4A4A and PICALM. Our approach is available in the Bioconductor package regionalpcs, opening avenues for research and facilitating a deeper understanding of the epigenetic landscape in complex diseases.

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

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

imply: improving cell-type deconvolution accuracy using personalized reference profiles

Using computational tools, bulk transcriptomics can be deconvoluted to estimate the abundance of constituent cell types. However, existing deconvolution methods are conditioned on the assumption that the whole study population is served by a single reference panel, ignoring person-to-person heterogeneity. Here, we present imply, a novel algorithm to deconvolute cell type proportions using personalized reference panels. Simulation studies demonstrate reduced bias compared with existing methods. Real data analyses on longitudinal consortia show disparities in cell type proportions are associated with several disease phenotypes in Type 1 diabetes and Parkinson's disease. imply is available through the R/Bioconductor package ISLET at https://bioconductor.org/packages/ISLET/ .

Benchmark of cellular deconvolution methods using a multi-assay reference dataset from postmortem human prefrontal cortex

Background

Cellular deconvolution of bulk RNA-sequencing (RNA-seq) data using single cell or nuclei RNA-seq (sc/snRNA-seq) reference data is an important strategy for estimating cell type composition in heterogeneous tissues, such as human brain. Computational methods for deconvolution have been developed and benchmarked against simulated data, pseudobulked sc/snRNA-seq data, or immunohistochemistry reference data. A major limitation in developing improved deconvolution algorithms has been the lack of integrated datasets with orthogonal measurements of gene expression and estimates of cell type proportions on the same tissue sample. Deconvolution algorithm performance has not yet been evaluated across different RNA extraction methods (cytosolic, nuclear, or whole cell RNA), different library preparation types (mRNA enrichment vs. ribosomal RNA depletion), or with matched single cell reference datasets.

Results

A rich multi-assay dataset was generated in postmortem human dorsolateral prefrontal cortex (DLPFC) from 22 tissue blocks. Assays included spatially-resolved transcriptomics, snRNA-seq, bulk RNA-seq (across six library/extraction RNA-seq combinations), and RNAScope/Immunofluorescence (RNAScope/IF) for six broad cell types. The Mean Ratio method, implemented in the DeconvoBuddies R package, was developed for selecting cell type marker genes. Six computational deconvolution algorithms were evaluated in DLPFC and predicted cell type proportions were compared to orthogonal RNAScope/IF measurements.

Conclusions

Bisque and hspe were the most accurate methods, were robust to differences in RNA library types and extractions. This multi-assay dataset showed that cell size differences, marker genes differentially quantified across RNA libraries, and cell composition variability in reference snRNA-seq impact the accuracy of current deconvolution methods.

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.

DESP demixes cell-state profiles from dynamic bulk molecular measurements

Understanding the dynamic expression of proteins and other key molecules driving phenotypic remodeling in development and pathobiology has garnered widespread interest, yet the exploration of these systems at the foundational resolution of the underlying cell states has been significantly limited by technical constraints. Here, we present DESP, an algorithm designed to leverage independent estimates of cell-state proportions, such as from single-cell RNA sequencing, to resolve the relative contributions of cell states to bulk molecular measurements, most notably quantitative proteomics, recorded in parallel. We applied DESP to an in vitro model of the epithelial-to-mesenchymal transition and demonstrated its ability to accurately reconstruct cell-state signatures from bulk-level measurements of both the proteome and transcriptome, providing insights into transient regulatory mechanisms. DESP provides a generalizable computational framework for modeling the relationship between bulk and single-cell molecular measurements, enabling the study of proteomes and other molecular profiles at the cell-state level using established bulk-level workflows.

CAM3.0: determining cell type composition and expression from bulk tissues with fully unsupervised deconvolution

MOTIVATION

Complex tissues are dynamic ecosystems consisting of molecularly distinct yet interacting cell types. Computational deconvolution aims to dissect bulk tissue data into cell type compositions and cell-specific expressions. With few exceptions, most existing deconvolution tools exploit supervised approaches requiring various types of references that may be unreliable or even unavailable for specific tissue microenvironments.

RESULTS

We previously developed a fully unsupervised deconvolution method-Convex Analysis of Mixtures (CAM), that enables estimation of cell type composition and expression from bulk tissues. We now introduce CAM3.0 tool that improves this framework with three new and highly efficient algorithms, namely, radius-fixed clustering to identify reliable markers, linear programming to detect an initial scatter simplex, and a smart floating search for the optimum latent variable model. The comparative experimental results obtained from both realistic simulations and case studies show that the CAM3.0 tool can help biologists more accurately identify known or novel cell markers, determine cell proportions, and estimate cell-specific expressions, complementing the existing tools particularly when study- or datatype-specific references are unreliable or unavailable.

AVAILABILITY AND IMPLEMENTATION

The open-source R Scripts of CAM3.0 is freely available at https://github.com/ChiungTingWu/CAM3/(https://github.com/Bioconductor/Contributions/issues/3205). A user's guide and a vignette are provided.

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.

Aspirin Challenge-Induced Genome-Wide DNA Methylation Profile of Peripheral Blood Lymphocytes in Aspirin-Exacerbated Respiratory Disease

Genetic variation and epigenetic factors are thought to contribute to the development of hypersensitivity to aspirin. DNA methylation fluctuates dynamically throughout the day. To discover new CpG methylation in lymphocytes associated with aspirin-exacerbated respiratory disease (AERD), we evaluated changes in global CpG methylation profiles from before to after an oral aspirin challenge in patients with AERD and aspirin-tolerant asthma (ATA). Whole-genome CpG methylation levels of peripheral blood mononuclear cells were quantified with an Illumina 860K Infinium Methylation EPIC BeadChip array and then adjusted for inferred lymphocyte fraction (ILF) with GLINT and Tensor Composition Analysis. Among the 866,091 CpGs in the array, differentially methylated CpGs (DMCs) were found in 6 CpGs in samples from all 12 patients with asthma included in the study (AERD, n = 6; ATA, n = 6). DMCs were found in 3 CpGs in the 6 ATA samples and in 615 CpGs in the 6 AERD samples. A total of 663 DMCs in 415 genes and 214 intergenic regions differed significantly in the AERD compared with the ATA. In promoters, 126 CpG loci were predicted to bind to 38 transcription factors (TFs), many of which were factors already known to be involved in the pathogenesis of asthma and immune responses. In conclusion, we identified 615 new CpGs methylated in peripheral blood lymphocytes by oral aspirin challenge in AERD but not in ATA. These findings indicate that oral aspirin challenge induces epigenetic changes in ILFs, specifically in AERD patients, possibly via changes in TF binding, which may have epigenetic effects on the development of AERD.

HLA-DRB1*15:01 and the MERTK Gene Interact to Selectively Influence the Profile of MERTK-Expressing Monocytes in Both Health and MS

BACKGROUND AND OBJECTIVES

HLA-DRB1*15:01 (DR15) and MERTK are 2 risk genes for multiple sclerosis (MS). The variant rs7422195 is an expression quantitative trait locus for MERTK in CD14+ monocytes; cells with phagocytic and immunomodulatory potential. We aimed to understand how drivers of disease risk and pathogenesis vary with HLA and MERTK genotype and disease activity.

METHODS

We investigated how proportions of monocytes vary with HLA and MERTK genotype and disease activity in MS. CD14+ monocytes were isolated from patients with MS at relapse (n = 40) and 3 months later (n = 23). Healthy controls (HCs) underwent 2 blood collections 3 months apart. Immunophenotypic profiling of monocytes was performed by flow cytometry. Methylation of 35 CpG sites within and near the MERTK gene was assessed in whole blood samples of individuals experiencing their first episode of clinical CNS demyelination (n = 204) and matched HCs (n = 345) using an Illumina EPIC array.

RESULTS

DR15-positive patients had lower proportions of CD14+ MERTK+ monocytes than DR15-negative patients, independent of genotype at the MERTK SNP rs7422195. Proportions of CD14+ MERTK+ monocytes were further reduced during relapse in DR15-positive but not DR15-negative patients. Patients homozygous for the major G allele at rs7422195 exhibited higher proportions of CD14+ MERTK+ monocytes at both relapse and remission compared with controls. We observed that increased methylation of the MERTK gene was significantly associated with the presence of DR15.

DISCUSSION

DR15 and MERTK genotype independently influence proportions of CD14+ MERTK+ monocytes in MS. We confirmed previous observations that the MERTK risk SNP rs7422195 is associated with altered MERTK expression in monocytes. We identified that expression of MERTK is stratified by disease in people homozygous for the major G allele of rs7422195. The finding that the proportion of CD14+ MERTK+ monocytes is reduced in DR15-positive individuals supports prior data identifying genetic links between these 2 loci in influencing MS risk. DR15 genotype-dependent alterations in methylation of the MERTK gene provides a molecular link between these loci and identifies a potential mechanism by which MERTK expression is influenced by DR15. This links DR15 haplotype to MS susceptibility beyond direct influence on antigen presentation and suggests the need for HLA-based stratification of approaches to MERTK as a therapeutic target.

Quantifying the proportion of different cell types in the human cortex using DNA methylation profiles

BACKGROUND

Due to interindividual variation in the cellular composition of the human cortex, it is essential that covariates that capture these differences are included in epigenome-wide association studies using bulk tissue. As experimentally derived cell counts are often unavailable, computational solutions have been adopted to estimate the proportion of different cell types using DNA methylation data. Here, we validate and profile the use of an expanded reference DNA methylation dataset incorporating two neuronal and three glial cell subtypes for quantifying the cellular composition of the human cortex.

RESULTS

We tested eight reference panels containing different combinations of neuronal- and glial cell types and characterised their performance in deconvoluting cell proportions from computationally reconstructed or empirically derived human cortex DNA methylation data. Our analyses demonstrate that while these novel brain deconvolution models produce accurate estimates of cellular proportions from profiles generated on postnatal human cortex samples, they are not appropriate for the use in prenatal cortex or cerebellum tissue samples. Applying our models to an extensive collection of empirical datasets, we show that glial cells are twice as abundant as neuronal cells in the human cortex and identify significant associations between increased Alzheimer's disease neuropathology and the proportion of specific cell types including a decrease in NeuNNeg/SOX10Neg nuclei and an increase of NeuNNeg/SOX10Pos nuclei.

CONCLUSIONS

Our novel deconvolution models produce accurate estimates for cell proportions in the human cortex. These models are available as a resource to the community enabling the control of cellular heterogeneity in epigenetic studies of brain disorders performed on bulk cortex tissue.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

Engineering Heterogeneous Tumor Models for Biomedical Applications

Tumor tissue engineering holds great promise for replicating the physiological and behavioral characteristics of tumors in vitro. Advances in this field have led to new opportunities for studying the tumor microenvironment and exploring potential anti-cancer therapeutics. However, the main obstacle to the widespread adoption of tumor models is the poor understanding and insufficient reconstruction of tumor heterogeneity. In this review, the current progress of engineering heterogeneous tumor models is discussed. First, the major components of tumor heterogeneity are summarized, which encompasses various signaling pathways, cell proliferations, and spatial configurations. Then, contemporary approaches are elucidated in tumor engineering that are guided by fundamental principles of tumor biology, and the potential of a bottom-up approach in tumor engineering is highlighted. Additionally, the characterization approaches and biomedical applications of tumor models are discussed, emphasizing the significant role of engineered tumor models in scientific research and clinical trials. Lastly, the challenges of heterogeneous tumor models in promoting oncology research and tumor therapy are described and key directions for future research are provided.

Fluorescent Reporter Systems to Investigate Chromatin Effector Proteins in Living Cells

Epigenetic research faces the challenge of the high complexity and tight regulation in chromatin modification networks. Although many isolated mechanisms of chromatin-mediated gene regulation have been described, solid approaches for the comprehensive analysis of specific processes as parts of the bigger epigenome network are missing. In order to expand the toolbox of methods by a system that will help to capture and describe the complexity of transcriptional regulation, we describe here a robust protocol for the generation of stable reporter systems for transcriptional activity and summarize their applications. The system allows for the induced recruitment of a chromatin regulator to a fluorescent reporter gene, followed by the detection of transcriptional changes using flow cytometry. The reporter gene is integrated into an endogenous chromatin environment, thus enabling the detection of regulatory dependencies of the investigated chromatin regulator on endogenous cofactors. The system allows for an easy and dynamic readout at the single-cell level and the ability to compensate for cell-to-cell variances of transcription. The modular design of the system enables the simple adjustment of the method for the investigation of different chromatin regulators in a broad panel of cell lines. We also summarize applications of this technology to characterize the silencing velocity of different chromatin effectors, removal of activating histone modifications, analysis of stability and reversibility of epigenome modifications, the investigation of the effects of small molecule on chromatin effectors and of functional effector-coregulator relationships. The presented method allows to investigate the complexity of transcriptional regulation by epigenetic effector proteins in living cells.

Uncertainty quantification of reference-based cellular deconvolution algorithms

The majority of epigenetic epidemiology studies to date have generated genome-wide profiles from bulk tissues (e.g., whole blood) however these are vulnerable to confounding from variation in cellular composition. Proxies for cellular composition can be mathematically derived from the bulk tissue profiles using a deconvolution algorithm; however, there is no method to assess the validity of these estimates for a dataset where the true cellular proportions are unknown. In this study, we describe, validate and characterize a sample level accuracy metric for derived cellular heterogeneity variables. The CETYGO score captures the deviation between a sample's DNA methylation profile and its expected profile given the estimated cellular proportions and cell type reference profiles. We demonstrate that the CETYGO score consistently distinguishes inaccurate and incomplete deconvolutions when applied to reconstructed whole blood profiles. By applying our novel metric to >6,300 empirical whole blood profiles, we find that estimating accurate cellular composition is influenced by both technical and biological variation. In particular, we show that when using a common reference panel for whole blood, less accurate estimates are generated for females, neonates, older individuals and smokers. Our results highlight the utility of a metric to assess the accuracy of cellular deconvolution, and describe how it can enhance studies of DNA methylation that are reliant on statistical proxies for cellular heterogeneity. To facilitate incorporating our methodology into existing pipelines, we have made it freely available as an R package (https://github.com/ds420/CETYGO).

Linking Prenatal Environmental Exposures to Lifetime Health with Epigenome-Wide Association Studies: State-of-the-Science Review and Future Recommendations

BACKGROUND

The prenatal environment influences lifetime health; epigenetic mechanisms likely predominate. In 2016, the first international consortium paper on cigarette smoking during pregnancy and offspring DNA methylation identified extensive, reproducible exposure signals. This finding raised expectations for epigenome-wide association studies (EWAS) of other exposures.

OBJECTIVE

We review the current state-of-the-science for DNA methylation associations across prenatal exposures in humans and provide future recommendations.

METHODS

We reviewed 134 prenatal environmental EWAS of DNA methylation in newborns, focusing on 51 epidemiological studies with meta-analysis or replication testing. Exposures spanned cigarette smoking, alcohol consumption, air pollution, dietary factors, psychosocial stress, metals, other chemicals, and other exogenous factors. Of the reproducible DNA methylation signatures, we examined implementation as exposure biomarkers.

RESULTS

Only 19 (14%) of these prenatal EWAS were conducted in cohorts of 1,000 or more individuals, reflecting the still early stage of the field. To date, the largest perinatal EWAS sample size was 6,685 participants. For comparison, the most recent genome-wide association study for birth weight included more than 300,000 individuals. Replication, at some level, was successful with exposures to cigarette smoking, folate, dietary glycemic index, particulate matter with aerodynamic diameter < 10 μ m and < 2.5 μ m , nitrogen dioxide, mercury, cadmium, arsenic, electronic waste, PFAS, and DDT. Reproducible effects of a more limited set of prenatal exposures (smoking, folate) enabled robust methylation biomarker creation.

DISCUSSION

Current evidence demonstrates the scientific premise for reproducible DNA methylation exposure signatures. Better powered EWAS could identify signatures across many exposures and enable comprehensive biomarker development. Whether methylation biomarkers of exposures themselves cause health effects remains unclear. We expect that larger EWAS with enhanced coverage of epigenome and exposome, along with improved single-cell technologies and evolving methods for integrative multi-omics analyses and causal inference, will expand mechanistic understanding of causal links between environmental exposures, the epigenome, and health outcomes throughout the life course. https://doi.org/10.1289/EHP12956.

imply: improving cell-type deconvolution accuracy using personalized reference profiles

Real-world clinical samples are often admixtures of signal mosaics from multiple pure cell types. Using computational tools, bulk transcriptomics can be deconvoluted to solve for the abundance of constituent cell types. However, existing deconvolution methods are conditioned on the assumption that the whole study population is served by a single reference panel, which ignores person-to-person heterogeneity. Here we present imply, a novel algorithm to deconvolute cell type proportions using personalized reference panels. imply can borrow information across repeatedly measured samples for each subject, and obtain precise cell type proportion estimations. Simulation studies demonstrate reduced bias in cell type abundance estimation compared with existing methods. Real data analyses on large longitudinal consortia show more realistic deconvolution results that align with biological facts. Our results suggest that disparities in cell type proportions are associated with several disease phenotypes in type 1 diabetes and Parkinson's disease. Our proposed tool imply is available through the R/Bioconductor package ISLET at https://bioconductor.org/packages/ISLET/.

Approximate estimation of cell-type resolution transcriptome in bulk tissue through matrix completion

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for uncovering cellular heterogeneity. However, the high costs associated with this technique have rendered it impractical for studying large patient cohorts. We introduce ENIGMA (Deconvolution based on Regularized Matrix Completion), a method that addresses this limitation through accurately deconvoluting bulk tissue RNA-seq data into a readout with cell-type resolution by leveraging information from scRNA-seq data. By employing a matrix completion strategy, ENIGMA minimizes the distance between the mixture transcriptome obtained with bulk sequencing and a weighted combination of cell-type-specific expression. This allows the quantification of cell-type proportions and reconstruction of cell-type-specific transcriptomes. To validate its performance, ENIGMA was tested on both simulated and real datasets, including disease-related tissues, demonstrating its ability in uncovering novel biological insights.

ISLET: individual-specific reference panel recovery improves cell-type-specific inference

We propose a statistical framework ISLET to infer individual-specific and cell-type-specific transcriptome reference panels. ISLET models the repeatedly measured bulk gene expression data, to optimize the usage of shared information within each subject. ISLET is the first available method to achieve individual-specific reference estimation in repeated samples. Using simulation studies, we show outstanding performance of ISLET in the reference estimation and downstream cell-type-specific differentially expressed genes testing. We apply ISLET to longitudinal transcriptomes profiled from blood samples in a large observational study of young children and confirm the cell-type-specific gene signatures for pancreatic islet autoantibody. ISLET is available at https://bioconductor.org/packages/ISLET .

EpiMix is an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer and immunological and cardiovascular diseases. Recent technological advances have enabled genome-wide profiling of DNAme in large human cohorts. There is a need for analytical methods that can more sensitively detect differential methylation profiles present in subsets of individuals from these heterogeneous, population-level datasets. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared with existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and long non-coding RNAs (lncRNAs). Using cell-type-specific data from two separate studies, we discover epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven ncRNAs in non-small cell lung cancer.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis

OBJECTIVE

We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment.

METHODS

Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software.

RESULTS

We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples.

CONCLUSION

We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

Methylation-Regulated Long Non-Coding RNA Expression in Ulcerative Colitis

Long non-coding RNAs (lncRNAs) have been shown to play a role in the pathogenesis of ulcerative colitis (UC). Although epigenetic processes such as DNA methylation and lncRNA expression are well studied in UC, the importance of the interplay between the two processes has not yet been fully explored. It is, therefore, believed that interactions between environmental factors and epigenetics contribute to disease development. Mucosal biopsies from 11 treatment-naïve UC patients and 13 normal controls were used in this study. From each individual sample, both whole-genome bisulfite sequencing data (WGBS) and lncRNA expression data were analyzed. Correlation analysis between lncRNA expression and upstream differentially methylated regions (DMRs) was used to identify lncRNAs that might be regulated by DMRs. Furthermore, proximal protein-coding genes associated with DMR-regulated lncRNAs were identified by correlating their expression. The study identified UC-associated lncRNAs such as MIR4435-2HG, ZFAS1, IL6-AS1, and Pvt1, which may be regulated by DMRs. Several genes that are involved in inflammatory immune responses were found downstream of DMR-regulated lncRNAs, including SERPINB1, CCL18, and SLC15A4. The interplay between lncRNA expression regulated by DNA methylation in UC might improve our understanding of UC pathogenesis.

In search of environmental risk factors for obsessive-compulsive disorder: study protocol for the OCDTWIN project

BACKGROUND

The causes of obsessive-compulsive disorder (OCD) remain unknown. Gene-searching efforts are well underway, but the identification of environmental risk factors is at least as important and should be a priority because some of them may be amenable to prevention or early intervention strategies. Genetically informative studies, particularly those employing the discordant monozygotic (MZ) twin design, are ideally suited to study environmental risk factors. This protocol paper describes the study rationale, aims, and methods of OCDTWIN, an open cohort of MZ twin pairs who are discordant for the diagnosis of OCD.

METHODS

OCDTWIN has two broad aims. In Aim 1, we are recruiting MZ twin pairs from across Sweden, conducting thorough clinical assessments, and building a biobank of biological specimens, including blood, saliva, urine, stool, hair, nails, and multimodal brain imaging. A wealth of early life exposures (e.g., perinatal variables, health-related information, psychosocial stressors) are available through linkage with the nationwide registers and the Swedish Twin Registry. Blood spots stored in the Swedish phenylketonuria (PKU) biobank will be available to extract DNA, proteins, and metabolites, providing an invaluable source of biomaterial taken at birth. In Aim 2, we will perform within-pair comparisons of discordant MZ twins, which will allow us to isolate unique environmental risk factors that are in the causal pathway to OCD, while strictly controlling for genetic and early shared environmental influences. To date (May 2023), 43 pairs of twins (21 discordant for OCD) have been recruited.

DISCUSSION

OCDTWIN hopes to generate unique insights into environmental risk factors that are in the causal pathway to OCD, some of which have the potential of being actionable targets.

Dynamics and consequences of nutrition-related microbial dysbiosis in early life: study protocol of the VITERBI GUT project

Introduction

Early life under- and overnutrition (jointly termed malnutrition) is increasingly recognized as an important risk factor for adult obesity and metabolic syndrome, a diet-related cluster of conditions including high blood sugar, fat and cholesterol. Nevertheless, the exact factors linking early life malnutrition with metabolic syndrome remain poorly characterized. We hypothesize that the microbiota plays a crucial role in this trajectory and that the pathophysiological mechanisms underlying under- and overnutrition are, to some extent, shared. We further hypothesize that a "dysbiotic seed microbiota" is transmitted to children during the birth process, altering the children's microbiota composition and metabolic health. The overall objective of this project is to understand the precise causes and biological mechanisms linking prenatal or early life under- or overnutrition with the predisposition to develop overnutrition and/or metabolic disease in later life, as well as to investigate the possibility of a dysbiotic seed microbiota inheritance in the context of maternal malnutrition.

Methods/design

VITERBI GUT is a prospective birth cohort allowing to study the link between early life malnutrition, the microbiota and metabolic health. VITERBI GUT will include 100 undernourished, 100 normally nourished and 100 overnourished pregnant women living in Vientiane, Lao People's Democratic Republic (PDR). Women will be recruited during their third trimester of pregnancy and followed with their child until its second birthday. Anthropometric, clinical, metabolic and nutritional data are collected from both the mother and the child. The microbiota composition of maternal and child's fecal and oral samples as well as maternal vaginal and breast milk samples will be determined using amplicon and shotgun metagenomic sequencing. Epigenetic modifications and lipid profiles will be assessed in the child's blood at 2 years of age. We will investigate for possible associations between metabolic health, epigenetics, and microbial changes.

Discussion

We expect the VITERBI GUT project to contribute to the emerging literature linking the early life microbiota, epigenetic changes and growth/metabolic health. We also expect this project to give new (molecular) insights into the mechanisms linking malnutrition-induced early life dysbiosis and metabolic health in later life, opening new avenues for microbiota-engineering using microbiota-targeted interventions.

In search of environmental risk factors for obsessive-compulsive disorder: Study protocol for the OCDTWIN project

Background The causes of obsessive-compulsive disorder (OCD) remain unknown. Gene-searching efforts are well underway, but the identification of environmental risk factors is at least as important and should be a priority because some of them may be amenable to prevention or early intervention strategies. Genetically informative studies, particularly those employing the discordant monozygotic (MZ) twin design, are ideally suited to study environmental risk factors. This protocol paper describes the study rationale, aims, and methods of OCDTWIN, an open cohort of MZ twin pairs who are discordant for the diagnosis of OCD. Methods OCDTWIN has two broad aims. In Aim 1, we are recruiting MZ twin pairs from across Sweden, conducting thorough clinical assessments, and building a biobank of biological specimens, including blood, saliva, urine, stool, hair, nails, and multimodal brain imaging. A wealth of early life exposures (e.g., perinatal variables, health-related information, psychosocial stressors) are available through linkage with the nationwide registers and the Swedish Twin Registry. Blood spots stored in the Swedish phenylketonuria (PKU) biobank will be available to extract DNA, proteins, and metabolites, providing an invaluable source of biomaterial taken at birth. In Aim 2, we will perform within-pair comparisons of discordant MZ twins, which will allow us to isolate unique environmental risk factors that are in the causal pathway to OCD, while strictly controlling for genetic and early shared environmental influences. To date (May 2023), 43 pairs of twins (21 discordant for OCD) have been recruited. Discussion OCDTWIN hopes to generate unique insights into environmental risk factors that are in the causal pathway to OCD, some of which have the potential of being actionable targets.

Asthma Exacerbations: The Genes Behind the Scenes

The clinical and socioeconomic burden of asthma exacerbations (AEs) constitutes a major public health problem. In the last 4 years, there has been an increase in ethnic diversity in candidate-gene and genome-wide association studies of AEs, which in the latter case led to the identification of novel genes and underlying pathobiological processes. Pharmacogenomics, admixture mapping analyses, and the combination of multiple "omics" layers have helped to prioritize genomic regions of interest and/or facilitated our understanding of the functional consequences of genetic variation. Nevertheless, the field still lags behind the genomics of asthma, where a vast compendium of genetic approaches has been used (eg, gene-environment nteractions, next-generation sequencing, and polygenic risk scores). Furthermore, the roles of the DNA methylome and histone modifications in AEs have received little attention, and microRNA findings remain to be validated in independent studies. Likewise, the most recent transcriptomic studies highlight the importance of the host-airway microbiome interaction in the modulation of risk of AEs. Leveraging -omics and deep-phenotyping data from subtypes or homogenous subgroups of patients will be crucial if we are to overcome the inherent heterogeneity of AEs, boost the identification of potential therapeutic targets, and implement precision medicine approaches to AEs in clinical practice.

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.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods have been developed to estimate cell-type proportions and gene expressions in bulk-tissue samples. However, the performance of these methods and their biological applications has not been evaluated, particularly on human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk-tissue RNAseq, single-cell/nuclei (sc/sn) RNAseq, and immunohistochemistry. A total of 1,130,767 nuclei/cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expression. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk-tissue or single-cell eQTLs alone. Differential gene expression associated with multiple phenotypes were also examined using the deconvoluted data. Our findings, which were replicated in bulk-tissue RNAseq and sc/snRNAseq data, provided new insights into the biological applications of deconvoluted data.

Hydroxymethylation alterations in progenitor-like cell types of pediatric central nervous system tumors are associated with cell type-specific transcriptional changes

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

CeDAR: incorporating cell type hierarchy improves cell type-specific differential analyses in bulk omics data

Bulk high-throughput omics data contain signals from a mixture of cell types. Recent developments of deconvolution methods facilitate cell type-specific inferences from bulk data. Our real data exploration suggests that differential expression or methylation status is often correlated among cell types. Based on this observation, we develop a novel statistical method named CeDAR to incorporate the cell type hierarchy in cell type-specific differential analyses of bulk data. Extensive simulation and real data analyses demonstrate that this approach significantly improves the accuracy and power in detecting cell type-specific differential signals compared with existing methods, especially in low-abundance cell types.

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.

An overview of DNA methylation-derived trait score methods and applications

Microarray technology has been used to measure genome-wide DNA methylation in thousands of individuals. These studies typically test the associations between individual DNA methylation sites ("probes") and complex traits or diseases. The results can be used to generate methylation profile scores (MPS) to predict outcomes in independent data sets. Although there are many parallels between MPS and polygenic (risk) scores (PGS), there are key differences. Here, we review motivations, methods, and applications of DNA methylation-based trait prediction, with a focus on common diseases. We contrast MPS with PGS, highlighting where assumptions made in genetic modeling may not hold in epigenetic data.

Rare diseases of epigenetic origin: Challenges and opportunities

Rare diseases (RDs), more than 80% of which have a genetic origin, collectively affect approximately 350 million people worldwide. Progress in next-generation sequencing technology has both greatly accelerated the pace of discovery of novel RDs and provided more accurate means for their diagnosis. RDs that are driven by altered epigenetic regulation with an underlying genetic basis are referred to as rare diseases of epigenetic origin (RDEOs). These diseases pose unique challenges in research, as they often show complex genetic and clinical heterogeneity arising from unknown gene-disease mechanisms. Furthermore, multiple other factors, including cell type and developmental time point, can confound attempts to deconvolute the pathophysiology of these disorders. These challenges are further exacerbated by factors that contribute to epigenetic variability and the difficulty of collecting sufficient participant numbers in human studies. However, new molecular and bioinformatics techniques will provide insight into how these disorders manifest over time. This review highlights recent studies addressing these challenges with innovative solutions. Further research will elucidate the mechanisms of action underlying unique RDEOs and facilitate the discovery of treatments and diagnostic biomarkers for screening, thereby improving health trajectories and clinical outcomes of affected patients.

MiXcan: a framework for cell-type-aware transcriptome-wide association studies with an application to breast cancer

Human bulk tissue samples comprise multiple cell types with diverse roles in disease etiology. Conventional transcriptome-wide association study approaches predict genetically regulated gene expression at the tissue level, without considering cell-type heterogeneity, and test associations of predicted tissue-level expression with disease. Here we develop MiXcan, a cell-type-aware transcriptome-wide association study approach that predicts cell-type-level expression, identifies disease-associated genes via combination of cell-type-level association signals for multiple cell types, and provides insight into the disease-critical cell type. As a proof of concept, we conducted cell-type-aware analyses of breast cancer in 58,648 women and identified 12 transcriptome-wide significant genes using MiXcan compared with only eight genes using conventional approaches. Importantly, MiXcan identified genes with distinct associations in mammary epithelial versus stromal cells, including three new breast cancer susceptibility genes. These findings demonstrate that cell-type-aware transcriptome-wide analyses can reveal new insights into the genetic and cellular etiology of breast cancer and other diseases.

A comprehensive assessment of cell type-specific differential expression methods in bulk data

Accounting for cell type compositions has been very successful at analyzing high-throughput data from heterogeneous tissues. Differential gene expression analysis at cell type level is becoming increasingly popular, yielding biomarker discovery in a finer granularity within a particular cell type. Although several computational methods have been developed to identify cell type-specific differentially expressed genes (csDEG) from RNA-seq data, a systematic evaluation is yet to be performed. Here, we thoroughly benchmark six recently published methods: CellDMC, CARseq, TOAST, LRCDE, CeDAR and TCA, together with two classical methods, csSAM and DESeq2, for a comprehensive comparison. We aim to systematically evaluate the performance of popular csDEG detection methods and provide guidance to researchers. In simulation studies, we benchmark available methods under various scenarios of baseline expression levels, sample sizes, cell type compositions, expression level alterations, technical noises and biological dispersions. Real data analyses of three large datasets on inflammatory bowel disease, lung cancer and autism provide evaluation in both the gene level and the pathway level. We find that csDEG calling is strongly affected by effect size, baseline expression level and cell type compositions. Results imply that csDEG discovery is a challenging task itself, with room to improvements on handling low signal-to-noise ratio and low expression genes.

EpiMix: an integrative tool for epigenomic subtyping using DNA methylation

DNA methylation (DNAme) is a major epigenetic factor influencing gene expression with alterations leading to cancer, immunological, and cardiovascular diseases. Recent technological advances enable genome-wide quantification of DNAme in large human cohorts. So far, existing methods have not been evaluated to identify differential DNAme present in large and heterogeneous patient cohorts. We developed an end-to-end analytical framework named "EpiMix" for population-level analysis of DNAme and gene expression. Compared to existing methods, EpiMix showed higher sensitivity in detecting abnormal DNAme that was present in only small patient subsets. We extended the model-based analyses of EpiMix to cis-regulatory elements within protein-coding genes, distal enhancers, and genes encoding microRNAs and lncRNAs. Using cell-type specific data from two separate studies, we discovered novel epigenetic mechanisms underlying childhood food allergy and survival-associated, methylation-driven non-coding RNAs in non-small cell lung cancer.

Benchmarking and integration of methods for deconvoluting spatial transcriptomic data

MOTIVATION

The rapid development of spatial transcriptomics (ST) approaches has provided new insights into understanding tissue architecture and function. However, the gene expressions measured at a spot may contain contributions from multiple cells due to the low-resolution of current ST technologies. Although many computational methods have been developed to disentangle discrete cell types from spatial mixtures, the community lacks a thorough evaluation of the performance of those deconvolution methods.

RESULTS

Here, we present a comprehensive benchmarking of 14 deconvolution methods on four datasets. Furthermore, we investigate the robustness of different methods to sequencing depth, spot size and the choice of normalization. Moreover, we propose a new ensemble learning-based deconvolution method (EnDecon) by integrating multiple individual methods for more accurate deconvolution. The major new findings include: (i) cell2loction, RCTD and spatialDWLS are more accurate than other ST deconvolution methods, based on the evaluation of three metrics: RMSE, PCC and JSD; (ii) cell2location and spatialDWLS are more robust to the variation of sequencing depth than RCTD; (iii) the accuracy of the existing methods tends to decrease as the spot size becomes smaller; (iv) most deconvolution methods perform best when they normalize ST data using the method described in their original papers; and (v) the integrative method, EnDecon, could achieve more accurate ST deconvolution. Our study provides valuable information and guideline for practically applying ST deconvolution tools and developing new and more effective methods.

AVAILABILITY AND IMPLEMENTATION

The benchmarking pipeline is available at https://github.com/SunXQlab/ST-deconvoulution. An R package for EnDecon is available at https://github.com/SunXQlab/EnDecon.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

The role of maternal DNA methylation in pregnancies complicated by gestational diabetes

Diabetes in pregnancy is associated with adverse pregnancy outcomes and poses a serious threat to the health of mother and child. Although the pathophysiological mechanisms that underlie the association between maternal diabetes and pregnancy complications have not yet been elucidated, it has been suggested that the frequency and severity of pregnancy complications are linked to the degree of hyperglycemia. Epigenetic mechanisms reflect gene-environment interactions and have emerged as key players in metabolic adaptation to pregnancy and the development of complications. DNA methylation, the best characterized epigenetic mechanism, has been reported to be dysregulated during various pregnancy complications, including pre-eclampsia, hypertension, diabetes, early pregnancy loss and preterm birth. The identification of altered DNA methylation patterns may serve to elucidate the pathophysiological mechanisms that underlie the different types of maternal diabetes during pregnancy. This review aims to provide a summary of existing knowledge on DNA methylation patterns in pregnancies complicated by pregestational type 1 (T1DM) and type 2 diabetes mellitus (T2DM), and gestational diabetes mellitus (GDM). Four databases, CINAHL, Scopus, PubMed and Google Scholar, were searched for studies on DNA methylation profiling in pregnancies complicated with diabetes. A total of 1985 articles were identified, of which 32 met the inclusion criteria and are included in this review. All studies profiled DNA methylation during GDM or impaired glucose tolerance (IGT), while no studies investigated T1DM or T2DM. We highlight the increased methylation of two genes, Hypoxia‐inducible Factor‐3α (HIF3α) and Peroxisome Proliferator-activated Receptor Gamma-coactivator-Alpha (PGC1-α), and the decreased methylation of one gene, Peroxisome Proliferator Activated Receptor Alpha (PPARα), in women with GDM compared to pregnant women with normoglycemia that were consistently methylated across diverse populations with varying pregnancy durations, and using different diagnostic criteria, methodologies and biological sources. These findings support the candidacy of these three differentially methylated genes as biomarkers for GDM. Furthermore, these genes may provide insight into the pathways that are epigenetically influenced during maternal diabetes and which should be prioritized and replicated in longitudinal studies and in larger populations to ensure their clinical applicability. Finally, we discuss the challenges and limitations of DNA methylation analysis, and the need for DNA methylation profiling to be conducted in different types of maternal diabetes in pregnancy.

Epigenetic signature of exposure to maternal Trypanosoma cruzi infection in cord blood cells from uninfected newborns

Aims: To assess the epigenetic effects of in utero exposure to maternal Trypanosoma cruzi infection. Methods: We performed an epigenome-wide association study to compare the DNA methylation patterns of umbilical cord blood cells from uninfected babies from chagasic and uninfected mothers. DNA methylation was measured using Infinium EPIC arrays. Results: We identified a differential DNA methylation signature of fetal exposure to maternal T. cruzi infection, in the absence of parasite transmission, with 12 differentially methylated sites in B cells and CD4⁺ T cells, including eight protein-coding genes. Conclusion: These genes participate in hematopoietic cell differentiation and the immune response and may be involved in immune disorders. They also have been associated with several developmental disorders and syndromes.

Schizophrenia-associated differential DNA methylation in brain is distributed across the genome and annotated to MAD1L1, a locus at which DNA methylation and transcription phenotypes share genetic variation with schizophrenia risk

DNA methylation (DNAm), the addition of a methyl group to a cytosine in DNA, plays an important role in the regulation of gene expression. Single-nucleotide polymorphisms (SNPs) associated with schizophrenia (SZ) by genome-wide association studies (GWAS) often influence local DNAm levels. Thus, DNAm alterations, acting through effects on gene expression, represent one potential mechanism by which SZ-associated SNPs confer risk. In this study, we investigated genome-wide DNAm in postmortem superior temporal gyrus from 44 subjects with SZ and 44 non-psychiatric comparison subjects using Illumina Infinium MethylationEPIC BeadChip microarrays, and extracted cell-type-specific methylation signals by applying tensor composition analysis. We identified SZ-associated differential methylation at 242 sites, and 44 regions containing two or more sites (FDR cutoff of q = 0.1) and determined a subset of these were cell-type specific. We found mitotic arrest deficient 1-like 1 (MAD1L1), a gene within an established GWAS risk locus, harbored robust SZ-associated differential methylation. We investigated the potential role of MAD1L1 DNAm in conferring SZ risk by assessing for colocalization among quantitative trait loci for methylation and gene transcripts (mQTLs and tQTLs) in brain tissue and GWAS signal at the locus using multiple-trait-colocalization analysis. We found that mQTLs and tQTLs colocalized with the GWAS signal (posterior probability >0.8). Our findings suggest that alterations in MAD1L1 methylation and transcription may mediate risk for SZ at the MAD1L1-containing locus. Future studies to identify how SZ-associated differential methylation affects MAD1L1 biological function are indicated.

Genome-wide DNA methylation profile of peripheral blood lymphocytes from subjects with nonsteroidal anti-inflammatory drug-induced respiratory diseases

BACKGROUND

Significant changes in CpG methylation have been identified in nasal polyps, which are the main targets of nonsteroidal anti-inflammatory drug-exacerbated respiratory disease (NERD); however, these polyps are composed of various cellular components. In the present study, whole-genome CpG methylation in peripheral blood lymphocytes (PBLs) was analyzed to define the epigenetic changes in lymphocytes, which are the primary immune cells involved in NERD.

MATERIALS AND METHODS

Genomic DNA from peripheral blood mononuclear cells from 27 NERD and 24 aspirin-tolerant asthma (ATA) was subjected to bisulfate conversion and a methylation array. Quantitative CpG methylation, the β-values as a quantitative measure of DNA methylation, in lymphocytes were calculated after adjustments for cellular composition.

RESULTS

Fifty-six hypermethylated and three hypomethylated differentially methylated CpGs (DMCs) in PBLs in the NERD compared with ATA. The top 10 CpG loci predicted the methylation risk score, with a positive predictive value of 91.3%, a negative predictive value of 81.5% and an accuracy of 84.3%. As demonstrated in the nasal polyps, 30 DMCs were predicted to bind to the following 10 transcription factors, ranked in descending order: AP-2alphaA, TFII-1, STAT4, FOXP3, GR, c-Est-1, E2F-1, XBP1, ENKTF-1 and NF-1. Gene ontology analysis identified 13 categories such as regulation of T-helper 17 cell differentiation, including SMAD7 and NFKBIZ. PBLs in NERD contained no DMCs in genes associated with the prostaglandin and leukotriene pathways, which were found in ATA.

CONCLUSION

PBLs in NERD form a unique pattern of DNA CpG methylation, and the combined analysis may provide predictive values for NERD.