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.

DNA Methylation Signatures of Multiple Sclerosis Occur Independently of Known Genetic Risk and Are Primarily Attributed to B Cells and Monocytes

Epigenetic mechanisms can regulate how DNA is expressed independently of sequence and are known to be associated with various diseases. Among those epigenetic mechanisms, DNA methylation (DNAm) is influenced by genotype and the environment, making it an important molecular interface for studying disease etiology and progression. In this study, we examined the whole blood DNA methylation profiles of a large group of people with (pw) multiple sclerosis (MS) compared to those of controls. We reveal that methylation differences in pwMS occur independently of known genetic risk loci and show that they more strongly differentiate disease (AUC = 0.85, 95% CI 0.82-0.89, p = 1.22 × 10-29) than known genetic risk loci (AUC = 0.72, 95% CI: 0.66-0.76, p = 9.07 × 10-17). We also show that methylation differences in MS occur predominantly in B cells and monocytes and indicate the involvement of cell-specific biological pathways. Overall, this study comprehensively characterizes the immune cell-specific epigenetic architecture of MS.

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.

Mature neutrophils and a NF-κB-to-IFN transition determine the unifying disease recovery dynamics in COVID-19

Disease recovery dynamics are often difficult to assess, as patients display heterogeneous recovery courses. To model recovery dynamics, exemplified by severe COVID-19, we apply a computational scheme on longitudinally sampled blood transcriptomes, generating recovery states, which we then link to cellular and molecular mechanisms, presenting a framework for studying the kinetics of recovery compared with non-recovery over time and long-term effects of the disease. Specifically, a decrease in mature neutrophils is the strongest cellular effect during recovery, with direct implications on disease outcome. Furthermore, we present strong indications for global regulatory changes in gene programs, decoupled from cell compositional changes, including an early rise in T cell activation and differentiation, resulting in immune rebalancing between interferon and NF-κB activity and restoration of cell homeostasis. Overall, we present a clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

Refined expression quantitative trait locus analysis on adenocarcinoma at the gastroesophageal junction reveals susceptibility and prognostic markers

Objectives: This study aimed to explore cell type level expression quantitative trait loci (eQTL) in adenocarcinoma at the gastroesophageal junction (ACGEJ) and identify susceptibility and prognosis markers. Methods: Whole-genome sequencing (WGS) was performed on 120 paired samples from Chinese ACGEJ patients. Germline mutations were detected by GATK tools. RNA sequencing (RNA-seq) data on ACGEJ samples were taken from our previous studies. Public single-cell RNA sequencing (scRNA-seq) data were used to produce the proportion of epithelial cells. Matrix eQTL and a linear mixed model were used to identify condition-specific cis-eQTLs. The R package coloc was used to perform co-localization analysis with the public data of genome-wide association studies (GWASs). Log-rank and Cox regression tests were used to identify survival-associated eQTL and genes. Functions of candidate risk loci were explored by experimental validation. Results: Refined eQTL analyses of paired ACGEJ samples were performed and 2,036 potential ACGEJ-specific eQTLs with East Asian specificity were identified in total. ACGEJ-gain eQTLs were enriched at promoter regions more than ACGEJ-loss eQTLs. rs658524 was identified as the top eQTL close to the transcription start site of its paired gene (CTSW). rs2240191-RASAL1, rs4236599-FOXP2, rs4947311-PSORS1C1, rs13134812-LOC391674, and rs17508585-CDK13-DT were identified as ACGEJ-specific susceptibility eQTLs. rs309483-LINC01355 was associated with the overall survival of ACGEJ patients. We explored functions of candidate eQTLs such as rs658524, rs309483, rs2240191, and rs4947311 by experimental validation. Conclusion: This study provides new risk loci for ACGEJ susceptibility and effective disease prognosis biomarkers.

Decosus: An R Framework for Universal Integration of Cell Proportion Estimation Methods

The assessment of the cellular heterogeneity and abundance in bulk tissue samples is essential for characterising cellular and organismal states. Computational approaches to estimate cellular abundance from bulk RNA-Seq datasets have variable performances, often requiring benchmarking matrices to select the best performing methods for individual studies. However, such benchmarking investigations are difficult to perform and assess in typical applications because of the absence of gold standard/ground-truth cellular measurements. Here we describe Decosus, an R package that integrates seven methods and signatures for deconvoluting cell types from gene expression profiles (GEP). Benchmark analysis on a range of datasets with ground-truth measurements revealed that our integrated estimates consistently exhibited stable performances across datasets than individual methods and signatures. We further applied Decosus to characterise the immune compartment of skin samples in different settings, confirming the well-established Th1 and Th2 polarisation in psoriasis and atopic dermatitis, respectively. Secondly, we revealed immune system-related UV-induced changes in sun-exposed skin. Furthermore, a significant motivation in the design of Decosus is flexibility and the ability for the user to include new gene signatures, algorithms, and integration methods at run time.

Cell deconvolution-based integrated time-series network of whole blood transcriptome reveals systemic antiviral activities and cell-specific immunological changes against PRRSV infection

Porcine reproductive and respiratory syndrome (PRRS) causes significant economic losses in the swine industry. However, the molecular mechanisms behind the common and cell type-specific systemic responses during PRRS virus (PRRSV) infection are not well understood. In this study, we collected viremia data, antibody levels, and whole-blood RNA-seq data obtained from eight PRRSV-infected piglets. We utilised a cell deconvolution approach to calculate cell type enrichment, constructed a time-serial gene co-expression network with differentially expressed genes, and conducted functional annotations. Three significant modules were identified within the network. The changes associated with viremia revealed an upregulated expression of genes related to antiviral activity. In the T-cell- and NK-cell-specific modules, infection led to an increased T-cell population and upregulation of genes related to T-cell defence responses. Conversely, in the monocyte- and neutrophil-specific module, genes involved in inflammatory responses were downregulated due to a decrease in monocyte proportion. This study highlights the time-series antiviral activities associated with viremia and the transcriptomic changes associated with immune responses in specific cell types. The findings provide comprehensive insights into host responses to PRRSV infection, including diagnostic biomarkers.

Network analyses of upper and lower airway transcriptomes identify shared mechanisms among children with recurrent wheezing and school-age asthma

Background

Predicting which preschool children with recurrent wheezing (RW) will develop school-age asthma (SA) is difficult, highlighting the critical need to clarify the pathogenesis of RW and the mechanistic relationship between RW and SA. Despite shared environmental exposures and genetic determinants, RW and SA are usually studied in isolation. Based on network analysis of nasal and tracheal transcriptomes, we aimed to identify convergent transcriptomic mechanisms in RW and SA.

Methods

RNA-sequencing data from nasal and tracheal brushing samples were acquired from the Gene Expression Omnibus. Combined with single-cell transcriptome data, cell deconvolution was used to infer the composition of 18 cellular components within the airway. Consensus weighted gene co-expression network analysis was performed to identify consensus modules closely related to both RW and SA. Shared pathways underlying consensus modules between RW and SA were explored by enrichment analysis. Hub genes between RW and SA were identified using machine learning strategies and validated using external datasets and quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Finally, the potential value of hub genes in defining RW subsets was determined using nasal and tracheal transcriptome data.

Results

Co-expression network analysis revealed similarities in the transcriptional networks of RW and SA in the upper and lower airways. Cell deconvolution analysis revealed an increase in mast cell fraction but decrease in club cell fraction in both RW and SA airways compared to controls. Consensus network analysis identified two consensus modules highly associated with both RW and SA. Enrichment analysis of the two consensus modules indicated that fatty acid metabolism-related pathways were shared key signals between RW and SA. Furthermore, machine learning strategies identified five hub genes, i.e., CST1, CST2, CST4, POSTN, and NRTK2, with the up-regulated hub genes in RW and SA validated using three independent external datasets and qRT-PCR. The gene signatures of the five hub genes could potentially be used to determine type 2 (T2)-high and T2-low subsets in preschoolers with RW.

Conclusions

These findings improve our understanding of the molecular pathogenesis of RW and provide a rationale for future exploration of the mechanistic relationship between RW and SA.

scDeconv: an R package to deconvolve bulk DNA methylation data with scRNA-seq data and paired bulk RNA-DNA methylation data

Many DNA methylation (DNAm) data are from tissues composed of various cell types, and hence cell deconvolution methods are needed to infer their cell compositions accurately. However, a bottleneck for DNAm data is the lack of cell-type-specific DNAm references. On the other hand, scRNA-seq data are being accumulated rapidly with various cell-type transcriptomic signatures characterized, and also, many paired bulk RNA-DNAm data are publicly available currently. Hence, we developed the R package scDeconv to use these resources to solve the reference deficiency problem of DNAm data and deconvolve them from scRNA-seq data in a trans-omics manner. It assumes that paired samples have similar cell compositions. So the cell content information deconvolved from the scRNA-seq and paired RNA data can be transferred to the paired DNAm samples. Then an ensemble model is trained to fit these cell contents with DNAm features and adjust the paired RNA deconvolution in a co-training manner. Finally, the model can be used on other bulk DNAm data to predict their relative cell-type abundances. The effectiveness of this method is proved by its accurate deconvolution on the three testing datasets here, and if given an appropriate paired dataset, scDeconv can also deconvolve other omics, such as ATAC-seq data. Furthermore, the package also contains other functions, such as identifying cell-type-specific inter-group differential features from bulk DNAm data. scDeconv is available at: https://github.com/yuabrahamliu/scDeconv.

Transcriptomics profiling of the non-small cell lung cancer microenvironment across disease stages reveals dual immune cell-type behaviors

Background

Lung cancer is the leading cause of cancer death worldwide, with poor survival despite recent therapeutic advances. A better understanding of the complexity of the tumor microenvironment is needed to improve patients' outcome.

Methods

We applied a computational immunology approach (involving immune cell proportion estimation by deconvolution, transcription factor activity inference, pathways and immune scores estimations) in order to characterize bulk transcriptomics of 62 primary lung adenocarcinoma (LUAD) samples from patients across disease stages. Focusing specifically on early stage samples, we validated our findings using an independent LUAD cohort with 70 bulk RNAseq and 15 scRNAseq datasets and on TCGA datasets.

Results

Through our methodology and feature integration pipeline, we identified groups of immune cells related to disease stage as well as potential immune response or evasion and survival. More specifically, we reported a duality in the behavior of immune cells, notably natural killer (NK) cells, which was shown to be associated with survival and could be relevant for immune response or evasion. These distinct NK cell populations were further characterized using scRNAseq data, showing potential differences in their cytotoxic activity.

Conclusion

The dual profile of several immune cells, most notably T-cell populations, have been discussed in the context of diseases such as cancer. Here, we report the duality of NK cells which should be taken into account in conjunction with other immune cell populations and behaviors in predicting prognosis, immune response or evasion.

Information about immune cell proportions and tumor stage improves the prediction of recurrence in patients with colorectal cancer

Predicting cancer recurrence is essential to improving the clinical outcomes of patients with colorectal cancer (CRC). Although tumor stage information has been used as a guideline to predict CRC recurrence, patients with the same stage show different clinical outcomes. Therefore, there is a need to develop a method to identify additional features for CRC recurrence prediction. Here, we developed a network-integrated multiomics (NIMO) approach to select appropriate transcriptome signatures for better CRC recurrence prediction by comparing the methylation signatures of immune cells. We validated the performance of the CRC recurrence prediction based on two independent retrospective cohorts of 114 and 110 patients. Moreover, to confirm that the prediction was improved, we used both NIMO-based immune cell proportions and TNM (tumor, node, metastasis) stage data. This work demonstrates the importance of (1) using both immune cell composition and TNM stage data and (2) identifying robust immune cell marker genes to improve CRC recurrence prediction.

Characterization of Alternative Splicing in High-Risk Wilms' Tumors

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.