Computational deconvolution of transcriptomics data from mixed cell populations

Heterogeneous pseudobulk simulation enables realistic benchmarking of cell-type deconvolution methods

BACKGROUND

Computational cell type deconvolution enables the estimation of cell type abundance from bulk tissues and is important for understanding tissue microenviroment, especially in tumor tissues. With rapid development of deconvolution methods, many benchmarking studies have been published aiming for a comprehensive evaluation for these methods. Benchmarking studies rely on cell-type resolved single-cell RNA-seq data to create simulated pseudobulk datasets by adding individual cells-types in controlled proportions.

RESULTS

In our work, we show that the standard application of this approach, which uses randomly selected single cells, regardless of the intrinsic difference between them, generates synthetic bulk expression values that lack appropriate biological variance. We demonstrate why and how the current bulk simulation pipeline with random cells is unrealistic and propose a heterogeneous simulation strategy as a solution. The heterogeneously simulated bulk samples match up with the variance observed in real bulk datasets and therefore provide concrete benefits for benchmarking in several ways. We demonstrate that conceptual classes of deconvolution methods differ dramatically in their robustness to heterogeneity with reference-free methods performing particularly poorly. For regression-based methods, the heterogeneous simulation provides an explicit framework to disentangle the contributions of reference construction and regression methods to performance. Finally, we perform an extensive benchmark of diverse methods across eight different datasets and find BayesPrism and a hybrid MuSiC/CIBERSORTx approach to be the top performers.

CONCLUSIONS

Our heterogeneous bulk simulation method and the entire benchmarking framework is implemented in a user friendly package https://github.com/humengying0907/deconvBenchmarking and https://doi.org/10.5281/zenodo.8206516 , enabling further developments in deconvolution methods.

Assessing transcriptomic heterogeneity of single-cell RNASeq data by bulk-level gene expression data

BACKGROUND

Single-cell RNA sequencing (sc-RNASeq) data illuminate transcriptomic heterogeneity but also possess a high level of noise, abundant missing entries and sometimes inadequate or no cell type annotations at all. Bulk-level gene expression data lack direct information of cell population composition but are more robust and complete and often better annotated. We propose a modeling framework to integrate bulk-level and single-cell RNASeq data to address the deficiencies and leverage the mutual strengths of each type of data and enable a more comprehensive inference of their transcriptomic heterogeneity. Contrary to the standard approaches of factorizing the bulk-level data with one algorithm and (for some methods) treating single-cell RNASeq data as references to decompose bulk-level data, we employed multiple deconvolution algorithms to factorize the bulk-level data, constructed the probabilistic graphical models of cell-level gene expressions from the decomposition outcomes, and compared the log-likelihood scores of these models in single-cell data. We term this framework backward deconvolution as inference operates from coarse-grained bulk-level data to fine-grained single-cell data. As the abundant missing entries in sc-RNASeq data have a significant effect on log-likelihood scores, we also developed a criterion for inclusion or exclusion of zero entries in log-likelihood score computation.

RESULTS

We selected nine deconvolution algorithms and validated backward deconvolution in five datasets. In the in-silico mixtures of mouse sc-RNASeq data, the log-likelihood scores of the deconvolution algorithms were strongly anticorrelated with their errors of mixture coefficients and cell type specific gene expression signatures. In the true bulk-level mouse data, the sample mixture coefficients were unknown but the log-likelihood scores were strongly correlated with accuracy rates of inferred cell types. In the data of autism spectrum disorder (ASD) and normal controls, we found that ASD brains possessed higher fractions of astrocytes and lower fractions of NRGN-expressing neurons than normal controls. In datasets of breast cancer and low-grade gliomas (LGG), we compared the log-likelihood scores of three simple hypotheses about the gene expression patterns of the cell types underlying the tumor subtypes. The model that tumors of each subtype were dominated by one cell type persistently outperformed an alternative model that each cell type had elevated expression in one gene group and tumors were mixtures of those cell types. Superiority of the former model is also supported by comparing the real breast cancer sc-RNASeq clusters with those generated by simulated sc-RNASeq data.

CONCLUSIONS

The results indicate that backward deconvolution serves as a sensible model selection tool for deconvolution algorithms and facilitates discerning hypotheses about cell type compositions underlying heterogeneous specimens such as tumors.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

DeepDecon accurately estimates cancer cell fractions in bulk RNA-seq data

Understanding the cellular composition of a disease-related tissue is important in disease diagnosis, prognosis, and downstream treatment. Recent advances in single-cell RNA-sequencing (scRNA-seq) technique have allowed the measurement of gene expression profiles for individual cells. However, scRNA-seq is still too expensive to be used for large-scale population studies, and bulk RNA-seq is still widely used in such situations. An essential challenge is to deconvolve cellular composition for bulk RNA-seq data based on scRNA-seq data. Here, we present DeepDecon, a deep neural network model that leverages single-cell gene expression information to accurately predict the fraction of cancer cells in bulk tissues. It provides a refining strategy in which the cancer cell fraction is iteratively estimated by a set of trained models. When applied to simulated and real cancer data, DeepDecon exhibits superior performance compared to existing decomposition methods in terms of accuracy.

Dataset including whole blood gene expression profiles and matched leukocyte counts with utility for benchmarking cellular deconvolution pipelines

OBJECTIVES

Cellular deconvolution is a valuable computational process that can infer the cellular composition of heterogeneous tissue samples from bulk RNA-sequencing data. Benchmark testing is a crucial step in the development and evaluation of new cellular deconvolution algorithms, and also plays a key role in the process of building and optimizing deconvolution pipelines for specific experimental applications. However, few in vivo benchmarking datasets exist, particularly for whole blood, which is the single most profiled human tissue. Here, we describe a unique dataset containing whole blood gene expression profiles and matched circulating leukocyte counts from a large cohort of human donors with utility for benchmarking cellular deconvolution pipelines.

DATA DESCRIPTION

To produce this dataset, venous whole blood was sampled from 138 total donors recruited at an academic medical center. Genome-wide expression profiling was subsequently performed via next-generation RNA sequencing, and white blood cell differentials were collected in parallel using flow cytometry. The resultant final dataset contains donor-level expression data for over 45,000 protein coding and non-protein coding genes, as well as matched neutrophil, lymphocyte, monocyte, and eosinophil counts.

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/ .

Predicting drug response through tumor deconvolution by cancer cell lines

Large-scale cancer drug sensitivity data have become available for a collection of cancer cell lines, but only limited drug response data from patients are available. Bridging the gap in pharmacogenomics knowledge between in vitro and in vivo datasets remains challenging. In this study, we trained a deep learning model, Scaden-CA, for deconvoluting tumor data into proportions of cancer-type-specific cell lines. Then, we developed a drug response prediction method using the deconvoluted proportions and the drug sensitivity data from cell lines. The Scaden-CA model showed excellent performance in terms of concordance correlation coefficients (>0.9 for model testing) and the correctly deconvoluted rate (>70% across most cancers) for model validation using Cancer Cell Line Encyclopedia (CCLE) bulk RNA data. We applied the model to tumors in The Cancer Genome Atlas (TCGA) dataset and examined associations between predicted cell viability and mutation status or gene expression levels to understand underlying mechanisms of potential value for drug repurposing.

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.

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.

Measuring the composition of the tumor microenvironment with transcriptome analysis: past, present and future

Interactions between tumor cells and immune cells in the tumor microenvironment (TME) play a vital role the mechanisms of immune evasion, by which cancer cells escape immune elimination. Thus, the characterization and quantification of different components in the TME is a hot topic in molecular biology and drug discovery. Since the development of transcriptome sequencing in bulk tissue, single cells and spatial dimensions, there are increasing methods emerging to deconvolute and subtype the TME. This review discusses and compares such computational strategies and downstream subtyping analyses. Integrative analyses of the transcriptome with other data, such as epigenetics and T-cell receptor sequencing, are needed to obtain comprehensive knowledge of the dynamic TME.

SpatialDDLS: an R package to deconvolute spatial transcriptomics data using neural networks

SUMMARY

Spatial transcriptomics has changed our way to study tissue structure and cellular organization. However, there are still limitations in its resolution, and most available platforms do not reach a single cell resolution. To address this issue, we introduce SpatialDDLS, a fast neural network-based algorithm for cell type deconvolution of spatial transcriptomics data. SpatialDDLS leverages single-cell RNA sequencing data to simulate mixed transcriptional profiles with predefined cellular composition, which are subsequently used to train a fully connected neural network to uncover cell type diversity within each spot. By comparing it with two state-of-the-art spatial deconvolution methods, we demonstrate that SpatialDDLS is an accurate and fast alternative to the available state-of-the art tools.

AVAILABILITY AND IMPLEMENTATION

The R package SpatialDDLS is available via CRAN-The Comprehensive R Archive Network: https://CRAN.R-project.org/package=SpatialDDLS. A detailed manual of the main functionalities implemented in the package can be found at https://diegommcc.github.io/SpatialDDLS.

Highly Accurate Estimation of Cell Type Abundance in Bulk Tissues Based on Single-Cell Reference and Domain Adaptive Matching

Accurately identifies the cellular composition of complex tissues, which is critical for understanding disease pathogenesis, early diagnosis, and prevention. However, current methods for deconvoluting bulk RNA sequencing (RNA-seq) typically rely on matched single-cell RNA sequencing (scRNA-seq) as a reference, which can be limiting due to differences in sequencing distribution and the potential for invalid information from single-cell references. Hence, a novel computational method named SCROAM is introduced to address these challenges. SCROAM transforms scRNA-seq and bulk RNA-seq into a shared feature space, effectively eliminating distributional differences in the latent space. Subsequently, cell-type-specific expression matrices are generated from the scRNA-seq data, facilitating the precise identification of cell types within bulk tissues. The performance of SCROAM is assessed through benchmarking against simulated and real datasets, demonstrating its accuracy and robustness. To further validate SCROAM's performance, single-cell and bulk RNA-seq experiments are conducted on mouse spinal cord tissue, with SCROAM applied to identify cell types in bulk tissue. Results indicate that SCROAM is a highly effective tool for identifying similar cell types. An integrated analysis of liver cancer and primary glioblastoma is then performed. Overall, this research offers a novel perspective for delivering precise insights into disease pathogenesis and potential therapeutic strategies.

Adjustment of scRNA-seq data to improve cell-type decomposition of spatial transcriptomics

Most sequencing-based spatial transcriptomics (ST) technologies do not achieve single-cell resolution where each captured location (spot) may contain a mixture of cells from heterogeneous cell types, and several cell-type decomposition methods have been proposed to estimate cell type proportions of each spot by integrating with single-cell RNA sequencing (scRNA-seq) data. However, these existing methods did not fully consider the effect of distribution difference between scRNA-seq and ST data for decomposition, leading to biased cell-type-specific genes derived from scRNA-seq for ST data. To address this issue, we develop an instance-based transfer learning framework to adjust scRNA-seq data by ST data to correctly match cell-type-specific gene expression. We evaluate the effect of raw and adjusted scRNA-seq data on cell-type decomposition by eight leading decomposition methods using both simulated and real datasets. Experimental results show that data adjustment can effectively reduce distribution difference and improve decomposition, thus enabling for a more precise depiction on spatial organization of cell types. We highlight the importance of data adjustment in integrative analysis of scRNA-seq with ST data and provide guidance for improved cell-type decomposition.

Historical perspective and future directions: computational science in immuno-oncology

Immuno-oncology holds promise for transforming patient care having achieved durable clinical response rates across a variety of advanced and metastatic cancers. Despite these achievements, only a minority of patients respond to immunotherapy, underscoring the importance of elucidating molecular mechanisms responsible for response and resistance to inform the development and selection of treatments. Breakthroughs in molecular sequencing technologies have led to the generation of an immense amount of genomic and transcriptomic sequencing data that can be mined to uncover complex tumor-immune interactions using computational tools. In this review, we discuss existing and emerging computational methods that contextualize the composition and functional state of the tumor microenvironment, infer the reactivity and clonal dynamics from reconstructed immune cell receptor repertoires, and predict the antigenic landscape for immune cell recognition. We further describe the advantage of multi-omics analyses for capturing multidimensional relationships and artificial intelligence techniques for integrating omics data with histopathological and radiological images to encapsulate patterns of treatment response and tumor-immune biology. Finally, we discuss key challenges impeding their widespread use and clinical application and conclude with future perspectives. We are hopeful that this review will both serve as a guide for prospective researchers seeking to use existing tools for scientific discoveries and inspire the optimization or development of novel tools to enhance precision, ultimately expediting advancements in immunotherapy that improve patient survival and quality of life.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

Cell-Type Composition Affects Adipose Gene Expression Associations With Cardiometabolic Traits

Understanding differences in adipose gene expression between individuals with different levels of clinical traits may reveal the genes and mechanisms leading to cardiometabolic diseases. However, adipose is a heterogeneous tissue. To account for cell-type heterogeneity, we estimated cell-type proportions in 859 subcutaneous adipose tissue samples with bulk RNA sequencing (RNA-seq) using a reference single-nuclear RNA-seq data set. Cell-type proportions were associated with cardiometabolic traits; for example, higher macrophage and adipocyte proportions were associated with higher and lower BMI, respectively. We evaluated cell-type proportions and BMI as covariates in tests of association between >25,000 gene expression levels and 22 cardiometabolic traits. For >95% of genes, the optimal, or best-fit, models included BMI as a covariate, and for 79% of associations, the optimal models also included cell type. After adjusting for the optimal covariates, we identified 2,664 significant associations (P ≤ 2e-6) for 1,252 genes and 14 traits. Among genes proposed to affect cardiometabolic traits based on colocalized genome-wide association study and adipose expression quantitative trait locus signals, 25 showed a corresponding association between trait and gene expression levels. Overall, these results suggest the importance of modeling cell-type proportion when identifying gene expression associations with cardiometabolic traits.

ARTICLE HIGHLIGHTS

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/.

A fast and globally optimal solution for RNA-seq quantification

Alignment-based RNA-seq quantification methods typically involve a time-consuming alignment process prior to estimating transcript abundances. In contrast, alignment-free RNA-seq quantification methods bypass this step, resulting in significant speed improvements. Existing alignment-free methods rely on the Expectation-Maximization (EM) algorithm for estimating transcript abundances. However, EM algorithms only guarantee locally optimal solutions, leaving room for further accuracy improvement by finding a globally optimal solution. In this study, we present TQSLE, the first alignment-free RNA-seq quantification method that provides a globally optimal solution for transcript abundances estimation. TQSLE adopts a two-step approach: first, it constructs a k-mer frequency matrix A for the reference transcriptome and a k-mer frequency vector b for the RNA-seq reads; then, it directly estimates transcript abundances by solving the linear equation ATAx = ATb. We evaluated the performance of TQSLE using simulated and real RNA-seq data sets and observed that, despite comparable speed to other alignment-free methods, TQSLE outperforms them in terms of accuracy. TQSLE is freely available at https://github.com/yhg926/TQSLE.

Bulk brain tissue cell-type deconvolution with bias correction for single-nuclei RNA sequencing data using DeTREM

BACKGROUND

Quantifying cell-type abundance in bulk tissue RNA-sequencing enables researchers to better understand complex systems. Newer deconvolution methodologies, such as MuSiC, use cell-type signatures derived from single-cell RNA-sequencing (scRNA-seq) data to make these calculations. Single-nuclei RNA-sequencing (snRNA-seq) reference data can be used instead of scRNA-seq data for tissues such as human brain where single-cell data are difficult to obtain, but accuracy suffers due to sequencing differences between the technologies.

RESULTS

We propose a modification to MuSiC entitled 'DeTREM' which compensates for sequencing differences between the cell-type signature and bulk RNA-seq datasets in order to better predict cell-type fractions. We show DeTREM to be more accurate than MuSiC in simulated and real human brain bulk RNA-sequencing datasets with various cell-type abundance estimates. We also compare DeTREM to SCDC and CIBERSORTx, two recent deconvolution methods that use scRNA-seq cell-type signatures. We find that they perform well in simulated data but produce less accurate results than DeTREM when used to deconvolute human brain data.

CONCLUSION

DeTREM improves the deconvolution accuracy of MuSiC and outperforms other deconvolution methods when applied to snRNA-seq data. DeTREM enables accurate cell-type deconvolution in situations where scRNA-seq data are not available. This modification improves characterization cell-type specific effects in brain tissue and identification of cell-type abundance differences under various conditions.

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

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

The aging kidney is characterized by tubuloinflammaging, a phenotype associated with MHC-II gene expression

Introduction

Even during physiologic aging, the kidney experiences a loss of mass and a progressive functional decline. This is clinically relevant as it leads to an increased risk of acute and chronic kidney disease. The kidney tubular system plays an important role in the underlying aging process, but the involved cellular mechanisms remain largely elusive.

Methods

Kidneys of 3-, 12- and 24-month-old male C57BL/6J mice were used for RNA sequencing, histological examination, immunostaining and RNA-in-situ-hybridization. Single cell RNA sequencing data of differentially aged murine and human kidneys was analyzed to identify age-dependent expression patterns in tubular epithelial cells. Senescent and non-senescent primary tubular epithelial cells from mouse kidney were used for in vitro experiments.

Results

During normal kidney aging, tubular cells adopt an inflammatory phenotype, characterized by the expression of MHC class II related genes. In our analysis of bulk and single cell transcriptional data we found that subsets of tubular cells show an age-related expression of Cd74, H2-Eb1 and H2-Ab1 in mice and CD74, HLA-DQB1 and HLADRB1 in humans. Expression of MHC class II related genes was associated with a phenotype of tubular cell senescence, and the selective elimination of senescent cells reversed the phenotype. Exposure to the Cd74 ligand MIF promoted a prosenescent phenotype in tubular cell cultures.

Discussion

Together, these data suggest that during normal renal aging tubular cells activate a program of 'tubuloinflammaging', which might contribute to age-related phenotypical changes and to increased disease susceptibility.

Whole genome deconvolution unveils Alzheimer's resilient epigenetic signature

Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) accurately depicts the chromatin regulatory state and altered mechanisms guiding gene expression in disease. However, bulk sequencing entangles information from different cell types and obscures cellular heterogeneity. To address this, we developed Cellformer, a deep learning method that deconvolutes bulk ATAC-seq into cell type-specific expression across the whole genome. Cellformer enables cost-effective cell type-specific open chromatin profiling in large cohorts. Applied to 191 bulk samples from 3 brain regions, Cellformer identifies cell type-specific gene regulatory mechanisms involved in resilience to Alzheimer's disease, an uncommon group of cognitively healthy individuals that harbor a high pathological load of Alzheimer's disease. Cell type-resolved chromatin profiling unveils cell type-specific pathways and nominates potential epigenetic mediators underlying resilience that may illuminate therapeutic opportunities to limit the cognitive impact of the disease. Cellformer is freely available to facilitate future investigations using high-throughput bulk ATAC-seq data.

BEDwARS: a robust Bayesian approach to bulk gene expression deconvolution with noisy reference signatures

Differential gene expression in bulk transcriptomics data can reflect change of transcript abundance within a cell type and/or change in the proportions of cell types. Expression deconvolution methods can help differentiate these scenarios. BEDwARS is a Bayesian deconvolution method designed to address differences between reference signatures of cell types and corresponding true signatures underlying bulk transcriptomic profiles. BEDwARS is more robust to noisy reference signatures and outperforms leading in-class methods for estimating cell type proportions and signatures. Application of BEDwARS to dihydropyridine dehydrogenase deficiency identified the possible involvement of ciliopathy and impaired translational control in the etiology of the disorder.

Peripheral blood cellular dynamics of rheumatoid arthritis treatment informs about efficacy of response to disease modifying drugs

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4⁺ cells (T_h1 and T_h2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 T_h cells (T_h1 and T_h17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Molecular mechanisms of environmental exposures and human disease

A substantial proportion of disease risk for common complex disorders is attributable to environmental exposures and pollutants. An appreciation of how environmental pollutants act on our cells to produce deleterious health effects has led to advances in our understanding of the molecular mechanisms underlying the pathogenesis of chronic diseases, including cancer and cardiovascular, neurodegenerative and respiratory diseases. Here, we discuss emerging research on the interplay of environmental pollutants with the human genome and epigenome. We review evidence showing the environmental impact on gene expression through epigenetic modifications, including DNA methylation, histone modification and non-coding RNAs. We also highlight recent studies that evaluate recently discovered molecular processes through which the environment can exert its effects, including extracellular vesicles, the epitranscriptome and the mitochondrial genome. Finally, we discuss current challenges when studying the exposome - the cumulative measure of environmental influences over the lifespan - and its integration into future environmental health research.

Integrating genetics and transcriptomics to study major depressive disorder: a conceptual framework, bioinformatic approaches, and recent findings

Major depressive disorder (MDD) is a complex and heterogeneous psychiatric syndrome with genetic and environmental influences. In addition to neuroanatomical and circuit-level disturbances, dysregulation of the brain transcriptome is a key phenotypic signature of MDD. Postmortem brain gene expression data are uniquely valuable resources for identifying this signature and key genomic drivers in human depression; however, the scarcity of brain tissue limits our capacity to observe the dynamic transcriptional landscape of MDD. It is therefore crucial to explore and integrate depression and stress transcriptomic data from numerous, complementary perspectives to construct a richer understanding of the pathophysiology of depression. In this review, we discuss multiple approaches for exploring the brain transcriptome reflecting dynamic stages of MDD: predisposition, onset, and illness. We next highlight bioinformatic approaches for hypothesis-free, genome-wide analyses of genomic and transcriptomic data and their integration. Last, we summarize the findings of recent genetic and transcriptomic studies within this conceptual framework.

PUREE: accurate pan-cancer tumor purity estimation from gene expression data

Tumors are complex masses composed of malignant and non-malignant cells. Variation in tumor purity (proportion of cancer cells in a sample) can both confound integrative analysis and enable studies of tumor heterogeneity. Here we developed PUREE, which uses a weakly supervised learning approach to infer tumor purity from a tumor gene expression profile. PUREE was trained on gene expression data and genomic consensus purity estimates from 7864 solid tumor samples. PUREE predicted purity with high accuracy across distinct solid tumor types and generalized to tumor samples from unseen tumor types and cohorts. Gene features of PUREE were further validated using single-cell RNA-seq data from distinct tumor types. In a comprehensive benchmark, PUREE outperformed existing transcriptome-based purity estimation approaches. Overall, PUREE is a highly accurate and versatile method for estimating tumor purity and interrogating tumor heterogeneity from bulk tumor gene expression data, which can complement genomics-based approaches or be used in settings where genomic data is unavailable.

Sex-specific gene expression patterns in head and neck squamous cell carcinomas

Objective

The head and neck squamous cell carcinomas (HNSCCs) have higher incidence rates in men, but the reasons are still obscure. This study aimed to investigate the sex-specific gene expression patterns and predict the regulatory mechanisms.

Design

Data including clinical, survival, RNA-seq, miRNA, and methylation information were derived from The Cancer Genome Atlas (TCGA). A total of 131 paired male and female cases were included based on propensity score matching. We concentrated on the prognostic values of the sex-specific pathways enriched by differentially expressed genes (DEGs) and predicted the potential regulatory mechanisms from immune cell infiltration, ceRNA regulatory network, methylation, and differential coexpression analysis.

Results

Compared with females, males exhibited a lower activity of immune-related functions and higher activities of mitochondrial and ubiquitination functions. The pathway activities were associated with the prognosis of males but less relevant to females. We extracted eight pathways with sex-biased survival patterns, of which five were about down-regulated immune functions, and three were up-regulated pathways (GTP biosynthetic, DNA polymerase, and spliceosomal complex assembly). The five immune pathways were moderately or strongly correlated with the proportion of macrophages. We identified six over-expressed lncRNAs that might be involved in the regulation of the three up-regulated pathways. These lncRNAs exhibited a lower methylation density in males, which might account for their over-expression.

Conclusions

For HNSCCs, males were characterized by immunosuppression. It was a sign of unfavorable prognosis and might be associated the proportion of macrophages. LncRNAs and methylation might be involved in the regulation of these pathways.

Spatial RNA sequencing methods show high resolution of single cell in cancer metastasis and the formation of tumor microenvironment

Cancer metastasis often leads to death and therapeutic resistance. This process involves the participation of a variety of cell components, especially cellular and intercellular communications in the tumor microenvironment (TME). Using genetic sequencing technology to comprehensively characterize the tumor and TME is therefore key to understanding metastasis and therapeutic resistance. The use of spatial transcriptome sequencing enables the localization of gene expressions and cell activities in tissue sections. By examining the localization change as well as gene expression of these cells, it is possible to characterize the progress of tumor metastasis and TME formation. With improvements of this technology, spatial transcriptome sequencing technology has been extended from local regions to whole tissues, and from single sequencing technology to multimodal analysis combined with a variety of datasets. This has enabled the detection of every single cell in tissue slides, with high resolution, to provide more accurate predictive information for tumor treatments. In this review, we summarize the results of recent studies dealing with new multimodal methods and spatial transcriptome sequencing methods in tumors to illustrate recent developments in the imaging resolution of micro-tissues.

Clustering-independent estimation of cell abundances in bulk tissues using single-cell RNA-seq data

Single-cell RNA-sequencing has transformed the study of biological tissues by enabling transcriptomic characterizations of their constituent cell states. Computational methods for gene expression deconvolution use this information to infer the cell composition of related tissues profiled at the bulk level. However, current deconvolution methods are restricted to discrete cell types and have limited power to make inferences about continuous cellular processes like cell differentiation or immune cell activation. We present ConDecon, a clustering-independent method for inferring the likelihood for each cell in a single-cell dataset to be present in a bulk tissue. ConDecon represents an improvement in functionality and accuracy with respect to current deconvolution methods. Using ConDecon, we discover the implication of neurodegenerative microglial inflammatory pathways in the mesenchymal transformation of ependymoma, recapitulate spatial patterns of cell differentiation during zebrafish embryogenesis, and make temporal inferences from bulk ATAC-seq data. Overall, ConDecon significantly enhances our understanding of dynamic cellular processes within bulk tissue samples.

Gene Expression Changes in Cultured Reactive Rat Astrocyte Models and Comparison to Device-Associated Effects in the Brain

Implanted microelectrode arrays hold immense therapeutic potential for many neurodegenerative diseases. However, a foreign body response limits long-term device performance. Recent literature supports the role of astrocytes in the response to damage to the central nervous system (CNS) and suggests that reactive astrocytes exist on a spectrum of phenotypes, from beneficial to neurotoxic. The goal of our study was to gain insight into the subtypes of reactive astrocytes responding to electrodes implanted in the brain. In this study, we tested the transcriptomic profile of two reactive astrocyte culture models (cytokine cocktail or lipopolysaccharide, LPS) utilizing RNA sequencing, which we then compared to differential gene expression surrounding devices inserted into rat motor cortex via spatial transcriptomics. We interpreted changes in the genetic expression of the culture models to that of 24 hour, 1 week and 6 week rat tissue samples at multiple distances radiating from the injury site. We found overlapping expression of up to ∼250 genes between in vitro models and in vivo effects, depending on duration of implantation. Cytokine-induced cells shared more genes in common with chronically implanted tissue (≥1 week) in comparison to LPS-exposed cells. We revealed localized expression of a subset of these intersecting genes (e.g., Serping1, Chi3l1, and Cyp7b1) in regions of device-encapsulating, glial fibrillary acidic protein (GFAP)-expressing astrocytes identified with immunohistochemistry. We applied a factorization approach to assess the strength of the relationship between reactivity markers and the spatial distribution of GFAP-expressing astrocytes in vivo . We also provide lists of hundreds of differentially expressed genes between reactive culture models and untreated controls, and we observed 311 shared genes between the cytokine induced model and the LPS-reaction induced control model. Our results show that comparisons of reactive astrocyte culture models with spatial transcriptomics data can reveal new biomarkers of the foreign body response to implantable neurotechnology. These comparisons also provide a strategy to assess the development of in vitro models of the tissue response to implanted electrodes.

Bioinformatics Tools for Bulk Gene Expression Deconvolution in Diabetic Retinopathy

Retinal neovascularization is one of the leading causes of vision loss and a hallmark of proliferative diabetic retinopathy (PDR). The immune system is observed to be involved in the pathogenesis of diabetic retinopathy (DR). The specific immune cell type that contributes to retinal neovascularization can be identified via a bioinformatics analysis of RNA sequencing (RNA-seq) data, known as deconvolution analysis. Previous study has identified the infiltration of macrophages in the retina of rats with hypoxia-induced retinal neovascularization and patients with PDR through a deconvolution algorithm, known as CIBERSORTx. Here, we describe the protocols of using CIBERSORTx to perform the deconvolution analysis and downstream analysis of RNA-seq data.

Deconvolution analysis of cell-type expression from bulk tissues by integrating with single-cell expression reference

To understand phenotypic variations and key factors which affect disease susceptibility of complex traits, it is important to decipher cell-type tissue compositions. To study cellular compositions of bulk tissue samples, one can evaluate cellular abundances and cell-type-specific gene expression patterns from the tissue transcriptome profiles. We develop both fixed and mixed models to reconstruct cellular expression fractions for bulk-profiled samples by using reference single-cell (sc) RNA-sequencing (RNA-seq) reference data. In benchmark evaluations of estimating cellular expression fractions, the mixed-effect models provide similar results as an elegant machine learning algorithm named cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORTx), which is a well-known and reliable procedure to reconstruct cell-type abundances and cell-type-specific gene expression profiles. In real data analysis, the mixed-effect models outperform or perform similarly as CIBERSORTx. The mixed models perform better than the fixed models in both benchmark evaluations and data analysis. In simulation studies, we show that if the heterogeneity exists in scRNA-seq data, it is better to use mixed models with heterogeneous mean and variance-covariance. As a byproduct, the mixed models provide fractions of covariance between subject-specific gene expression and cell types to measure their correlations. The proposed mixed models provide a complementary tool to dissect bulk tissues using scRNA-seq data.

Toward Precision Radiotherapy: A Nonlinear Optimization Framework and an Accelerated Machine Learning Algorithm for the Deconvolution of Tumor-Infiltrating Immune Cells

Tumor-infiltrating immune cells (TIICs) form a critical part of the ecosystem surrounding a cancerous tumor. Recent advances in radiobiology have shown that, in addition to damaging cancerous cells, radiotherapy drives the upregulation of immunosuppressive and immunostimulatory TIICs, which in turn impacts treatment response. Quantifying TIICs in tumor samples could form an important predictive biomarker guiding patient stratification and the design of radiotherapy regimens and combined immune-radiation treatments. As a result of several limitations associated with experimental methods for quantifying TIICs and the availability of extensive gene sequencing data, deconvolution-based computational methods have appeared as a suitable alternative for quantifying TIICs. Accordingly, we introduce and discuss a nonlinear regression approach (remarkably different from the traditional linear modeling approach of current deconvolution-based methods) and a machine learning algorithm for approximating the solution of the resulting constrained optimization problem. This way, the deconvolution problem is treated naturally, given that the gene expression levels of pure and heterogenous samples do not have a strictly linear relationship. When applied across transcriptomics datasets, our approach, which also allows the coupling of different loss functions, yields results that closely match ground-truth values from experimental methods and exhibits superior performance over popular deconvolution-based methods.

Niche deconvolution of the glioblastoma proteome reveals a distinct infiltrative phenotype within the proneural transcriptomic subgroup

Glioblastoma is often subdivided into three transcriptional subtypes (classical, proneural, mesenchymal) based on bulk RNA signatures that correlate with distinct genetic and clinical features. Potential cellular-level differences of these subgroups, such as the relative proportions of glioblastoma’s hallmark histopathologic features (e.g. brain infiltration, microvascular proliferation), may provide insight into their distinct phenotypes but are, however, not well understood. Here we leverage machine learning and reference proteomic profiles derived from micro-dissected samples of these major histomorphologic glioblastoma features to deconvolute and estimate niche proportions in an independent proteogenomically-characterized cohort. This approach revealed a strong association of the proneural transcriptional subtype with a diffusely infiltrating phenotype. Similarly, enrichment of a microvascular proliferation proteomic signature was seen within the mesenchymal subtype. This study is the first to link differences in the cellular pathology signatures and transcriptional profiles of glioblastoma, providing potential new insights into the genetic drivers and poor treatment response of specific subsets of glioblastomas.

Early transcriptional responses of bronchial epithelial cells to whole cigarette smoke mirror those of in-vivo exposed human bronchial mucosa

BACKGROUND

Despite the well-known detrimental effects of cigarette smoke (CS), little is known about the complex gene expression dynamics in the early stages after exposure. This study aims to investigate early transcriptomic responses following CS exposure of airway epithelial cells in culture and compare these to those found in human CS exposure studies.

METHODS

Primary bronchial epithelial cells (PBEC) were differentiated at the air-liquid interface (ALI) and exposed to whole CS. Bulk RNA-sequencing was performed at 1 h, 4 h, and 24 h hereafter, followed by differential gene expression analysis. Results were additionally compared to data retrieved from human CS studies.

RESULTS

ALI-PBEC gene expression in response to CS was most significantly changed at 4 h after exposure. Early transcriptomic changes (1 h, 4 h post CS exposure) were related to oxidative stress, xenobiotic metabolism, higher expression of immediate early genes and pro-inflammatory pathways (i.e., Nrf2, AP-1, AhR). At 24 h, ferroptosis-associated genes were significantly increased, whereas PRKN, involved in removing dysfunctional mitochondria, was downregulated. Importantly, the transcriptome dynamics of the current study mirrored in-vivo human studies of acute CS exposure, chronic smokers, and inversely mirrored smoking cessation.

CONCLUSION

These findings show that early after CS exposure xenobiotic metabolism and pro-inflammatory pathways were activated, followed by activation of the ferroptosis-related cell death pathway. Moreover, significant overlap between these transcriptomic responses in the in-vitro model and human in-vivo studies was found, with an early response of ciliated cells. These results provide validation for the use of ALI-PBEC cultures to study the human lung epithelial response to inhaled toxicants.

GEOMETRIC STRUCTURE GUIDED MODEL AND ALGORITHMS FOR COMPLETE DECONVOLUTION OF GENE EXPRESSION DATA

Complete deconvolution analysis for bulk RNA-seq data is important and helpful to distinguish whether the differences of disease-associated GEPs (gene expression profiles) in tissues of patients and normal controls are due to changes in cellular composition of tissue samples, or due to GEPs changes in specific cells. One of the major techniques to perform complete deconvolution is nonnegative matrix factorization (NMF), which also has a wide-range of applications in the machine learning community. However, the NMF is a well-known strongly ill-posed problem, so a direct application of NMF to RNA-seq data will suffer severe difficulties in the interpretability of solutions. In this paper, we develop an NMF-based mathematical model and corresponding computational algorithms to improve the solution identifiability of deconvoluting bulk RNA-seq data. In our approach, we combine the biological concept of marker genes with the solvability conditions of the NMF theories, and develop a geometric structures guided optimization model. In this strategy, the geometric structure of bulk tissue data is first explored by the spectral clustering technique. Then, the identified information of marker genes is integrated as solvability constraints, while the overall correlation graph is used as manifold regularization. Both synthetic and biological data are used to validate the proposed model and algorithms, from which solution interpretability and accuracy are significantly improved.

Integrative analysis of drug response and clinical outcome in acute myeloid leukemia

Acute myeloid leukemia (AML) is a cancer of myeloid-lineage cells with limited therapeutic options. We previously combined ex vivo drug sensitivity with genomic, transcriptomic, and clinical annotations for a large cohort of AML patients, which facilitated discovery of functional genomic correlates. Here, we present a dataset that has been harmonized with our initial report to yield a cumulative cohort of 805 patients (942 specimens). We show strong cross-cohort concordance and identify features of drug response. Further, deconvoluting transcriptomic data shows that drug sensitivity is governed broadly by AML cell differentiation state, sometimes conditionally affecting other correlates of response. Finally, modeling of clinical outcome reveals a single gene, PEAR1, to be among the strongest predictors of patient survival, especially for young patients. Collectively, this report expands a large functional genomic resource, offers avenues for mechanistic exploration and drug development, and reveals tools for predicting outcome in AML.

Future directions for precision oncology in prostate cancer

Clinical genomic testing is becoming routine in prostate cancer, as biomarker-driven therapies such as poly-ADP ribose polymerase (PARP) inhibitors and anti-PD1 immunotherapy are now approved for select men with castration-resistant prostate cancer harboring alterations in DNA repair genes. Challenges for precision medicine in prostate cancer include an overall low prevalence of actionable genomic alterations and a still limited understanding of the impact of tumor heterogeneity and co-occurring alterations on treatment response and outcomes across diverse patient populations. Expanded tissue-based technologies such as whole-genome sequencing, transcriptome analysis, epigenetic analysis, and single-cell RNA sequencing have not yet entered the clinical realm and could potentially improve upon our understanding of how molecular features of tumors, intratumoral heterogeneity, and the tumor microenvironment impact therapy response and resistance. Blood-based technologies including cell-free DNA, circulating tumor cells (CTCs), and extracellular vesicles (EVs) are less invasive molecular profiling resources that could also help capture intraindividual tumor heterogeneity and track dynamic changes that occur in the context of specific therapies. Furthermore, molecular imaging is an important biomarker tool within the framework of prostate cancer precision medicine with a capability to detect heterogeneity across metastases and potential therapeutic targets less invasively. Here, we review recent technological advances that may help promote the future implementation and value of precision oncology testing for patients with advanced prostate cancer.

Spatial Transcriptomics as a Novel Approach to Redefine Electrical Stimulation Safety

Current standards for safe delivery of electrical stimulation to the central nervous system are based on foundational studies which examined post-mortem tissue for histological signs of damage. This set of observations and the subsequently proposed limits to safe stimulation, termed the “Shannon limits,” allow for a simple calculation (using charge per phase and charge density) to determine the intensity of electrical stimulation that can be delivered safely to brain tissue. In the three decades since the Shannon limits were reported, advances in molecular biology have allowed for more nuanced and detailed approaches to be used to expand current understanding of the physiological effects of stimulation. Here, we demonstrate the use of spatial transcriptomics (ST) in an exploratory investigation to assess the biological response to electrical stimulation in the brain. Electrical stimulation was delivered to the rat visual cortex with either acute or chronic electrode implantation procedures. To explore the influence of device type and stimulation parameters, we used carbon fiber ultramicroelectrode arrays (7 μm diameter) and microwire electrode arrays (50 μm diameter) delivering charge and charge density levels selected above and below reported tissue damage thresholds (range: 2–20 nC, 0.1–1 mC/cm²). Spatial transcriptomics was performed using Visium Spatial Gene Expression Slides (10x Genomics, Pleasanton, CA, United States), which enabled simultaneous immunohistochemistry and ST to directly compare traditional histological metrics to transcriptional profiles within each tissue sample. Our data give a first look at unique spatial patterns of gene expression that are related to cellular processes including inflammation, cell cycle progression, and neuronal plasticity. At the acute timepoint, an increase in inflammatory and plasticity related genes was observed surrounding a stimulating electrode compared to a craniotomy control. At the chronic timepoint, an increase in inflammatory and cell cycle progression related genes was observed both in the stimulating vs. non-stimulating microwire electrode comparison and in the stimulating microwire vs. carbon fiber comparison. Using the spatial aspect of this method as well as the within-sample link to traditional metrics of tissue damage, we demonstrate how these data may be analyzed and used to generate new hypotheses and inform safety standards for stimulation in cortex.

Deconvolution of a Large Cohort of Placental Microarray Data Reveals Clinically Distinct Subtypes of Preeclampsia

It has been well established that the dysfunctional placenta plays an important role in the pathogenesis of preeclampsia (PE), a hypertensive disorder in pregnancy. However, it is not well understood how individual cell types in the placenta are involved in placenta dysfunction because of limited single-cell studies of placenta with PE. Given that a high-resolution single-cell atlas in the placenta is now available, deconvolution of publicly available bulk PE transcriptome data may provide us with the opportunity to investigate the contribution of individual placental cell types to PE. Recent benchmark studies on deconvolution have provided suggestions on the strategy of marker gene selection and the choice of methodologies. In this study, we experimented with these suggestions by using real bulk data with known cell-type proportions and established a deconvolution pipeline using CIBERSORT. Applying the deconvolution pipeline to a large cohort of PE placental microarray data, we found that the proportions of trophoblast cells in the placenta were significantly different between PE and normal controls. We then predicted cell-type-level expression profiles for each sample using CIBERSORTx and found that the activities of several canonical PE-related pathways were significantly altered in specific subtypes of trophoblasts in PE. Finally, we constructed an integrated expression profile for each PE sample by combining the predicted cell-type-level expression profiles of several clinically relevant placental cell types and identified four clusters likely representing four PE subtypes with clinically distinct features. As such, our study showed that deconvolution of a large cohort of placental microarray provided new insights about the molecular mechanism of PE that would not be obtained by analyzing bulk expression profiles.

Distinct cell type-specific protein signatures in GRN and MAPT genetic subtypes of frontotemporal dementia

Frontotemporal dementia is characterized by progressive atrophy of frontal and/or temporal cortices at an early age of onset. The disorder shows considerable clinical, pathological, and genetic heterogeneity. Here we investigated the proteomic signatures of frontal and temporal cortex from brains with frontotemporal dementia due to GRN and MAPT mutations to identify the key cell types and molecular pathways in their pathophysiology. We compared patients with mutations in the GRN gene (n = 9) or with mutations in the MAPT gene (n = 13) with non-demented controls (n = 11). Using quantitative proteomic analysis on laser-dissected tissues we identified brain region-specific protein signatures for both genetic subtypes. Using published single cell RNA expression data resources we deduced the involvement of major brain cell types in driving these different protein signatures. Subsequent gene ontology analysis identified distinct genetic subtype- and cell type-specific biological processes. For the GRN subtype, we observed a distinct role for immune processes related to endothelial cells and for mitochondrial dysregulation in neurons. For the MAPT subtype, we observed distinct involvement of dysregulated RNA processing, oligodendrocyte dysfunction, and axonal impairments. Comparison with an in-house protein signature of Alzheimer’s disease brains indicated that the observed alterations in RNA processing and oligodendrocyte function are distinct for the frontotemporal dementia MAPT subtype. Taken together, our results indicate the involvement of different brain cell types and biological mechanisms in genetic subtypes of frontotemporal dementia. Furthermore, we demonstrate that comparison of proteomic profiles of different disease entities can separate general neurodegenerative processes from disease-specific pathways, which may aid the development of disease subtype-specific treatment strategies.

Predicting Algorithm of Tissue Cell Ratio Based on Deep Learning Using Single-Cell RNA Sequencing

Background: Understanding the proportion of cell types in heterogeneous tissue samples is important in bioinformatics. It is a challenge to infer the proportion of tissues using bulk RNA sequencing data in bioinformatics because most traditional algorithms for predicting tissue cell ratios heavily rely on standardized specific cell-type gene expression profiles, and do not consider tissue heterogeneity. The prediction accuracy of algorithms is limited, and robustness is lacking. This means that new approaches are needed urgently. Methods: In this study, we introduced an algorithm that automatically predicts tissue cell ratios named Autoptcr. The algorithm uses the data simulated by single-cell RNA sequencing (ScRNA-Seq) for model training, using convolutional neural networks (CNNs) to extract intrinsic relationships between genes and predict the cell proportions of tissues. Results: We trained the algorithm using simulated bulk samples and made predictions using real bulk PBMC data. Comparing Autoptcr with existing advanced algorithms, the Pearson correlation coefficient between the actual value of Autoptcr and the predicted value was the highest, reaching 0.903. Tested on a bulk sample, the correlation coefficient of Lin was 41% higher than that of CSx. The algorithm can infer tissue cell proportions directly from tissue gene expression data. Conclusions: The Autoptcr algorithm uses simulated ScRNA-Seq data for training to solve the problem of specific cell-type gene expression profiles. It also has high prediction accuracy and strong noise resistance for the tissue cell ratio. This work is expected to provide new research ideas for the prediction of tissue cell proportions.

A novel integrative computational framework for breast cancer radiogenomic biomarker discovery

•

Bayesian tensor factorization is used to integrate multiomics data for radiogenomics analysis.

•

A regression framework is proposed to handle the unmatched data issue in radiogenomics analysis.

•

Deep learning is used to identify prognostic meaningful radiogenomic biomarkers for cancer.

In precise medicine, it is with great value to develop computational frameworks for identifying prognostic biomarkers which can capture both multi-genomic and phenotypic heterogeneity of breast cancer (BC). Radiogenomics is a field where medical images and genomic measurements are integrated and mined to solve challenging clinical problems. Previous radiogenomic studies suffered from data incompleteness, feature subjectivity and low interpretability. For example, the majority of the radiogenomic studies miss one or two of medical imaging data, genomic data, and clinical outcome data, which results in the data incomplete issue. Feature subjectivity issue comes from the extraction of imaging features with significant human involvement. Thus, there is an urgent need to address above-mentioned limitations so that fully automatic and transparent radiogenomic prognostic biomarkers could be identified for BC.

We proposed a novel framework for BC prognostic radiogenomic biomarker identification. This framework involves an explainable DL model for image feature extraction, a Bayesian tensor factorization (BTF) processing for multi-genomic feature extraction, a leverage strategy to utilize unpaired imaging, genomic, and survival outcome data, and a mediation analysis to provide further interpretation for identified biomarkers. This work provided a new perspective for conducting a comprehensive radiogenomic study when only limited resources are given. Compared with baseline traditional radiogenomic biomarkers, the 23 biomarkers identified by the proposed framework performed better in indicating patients’ survival outcome. And their interpretability is guaranteed by different levels of build-in and follow-up analyses.

Coding and non-coding co-expression network analysis identifies key modules and driver genes associated with precursor lesions of gastric cancer

BACKGROUND

Helicobacter pylori infection is the most important risk factor for gastric cancer (GC). Human gastric adenocarcinoma develops after long-term H. pylori infection via the Correa cascade. This carcinogenic pathway describes the progression from gastritis to atrophy, intestinal metaplasia (IM), dysplasia and GC. Patients with atrophy and intestinal metaplasia are considered to have precancerous lesions of GC (PLGC). H. pylori eradication and endoscopy surveillance are currently the only interventions for preventing GC. Better knowledge of the biology of human PLGC may help find stratification markers and contribute to better understanding of biological mechanisms. One way to achieve this is by using co-expression network analysis. Weighted gene co-expression network analysis (WGCNA) is often used to identify modules from co-expression networks and relate them to clinical traits. It also allows identification of driver genes that may be critical for PLGC.

AIM

The purpose of this study was to identify co-expression modules and differential gene expression in dyspeptic patients at different stages of the Correa pathway.

METHODS

We studied 96 gastric biopsies from 78 patients that were clinically classified as: non-active (n = 10) and chronic-active gastritis (n = 20), atrophy (n = 12), and IM (n = 36). Gene expression of coding RNAs was determined by microarrays and non-coding RNAs by RNA-seq. The WGCNA package was used for network construction, module detection, module preservation and hub and driver gene selection.

RESULTS

WGCNA identified 20 modules for coding RNAs and 4 for each miRNA and small RNA class. Modules were associated with antrum and corpus gastric locations, chronic gastritis and IM. Notably, coding RNA modules correlated with the Correa cascade. One was associated with the presence of H. pylori. In three modules, the module eigengene (ME) gradually increased in the stages toward IM, while in three others the inverse relationship was found. One miRNA module was negatively correlated to IM and was used for a mRNA-miRNA integration analysis. WGCNA also uncovered driver genes. Driver genes show both high connectivity within a module and are significantly associated with clinical traits. Some of those genes have been previously involved in H. pylori carcinogenesis, but others are new. Lastly, using similar external transcriptomic data, we confirmed that the discovered mRNA modules were highly preserved.

CONCLUSION

Our analysis captured co-expression modules that provide valuable information to understand the pathogenesis of the progression of PLGC.

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.

Tumor Microenvironment Profiling Identifies Prognostic Signatures and Suggests Immunotherapeutic Benefits in Neuroblastoma

The tumor microenvironment (TME) influences disease initiation and progression. Cross-talks of cells within TME can affect the efficacy of immunotherapies. However, a precise, concise, and comprehensive TME landscape in neuroblastoma (NB) has not been established. Here, we profiled the TME landscape of 498 NB-related patients on a self-curated gene list and identified three prognostic TMEsubgroups. The differentially expressed genes in these three TMEsubgroups were used to construct a genetic signature of the TME landscape and characterize three GeneSubgroups. The subgroup with the worst overall survival prognosis, the TMEsubgroup/GeneSubgroup3, lacked immune cell infiltration and received the highest scores of MYCN- and ALK-related signatures and lowest scores of immune pathways. Additionally, we found that the GeneSubgroup3 might be benefited from anti-GD2 instead of anti-PD-1 therapy. We further created a 48-gene signature, the TMEscore, to infer prognosis and validated it in three independent NB cohorts and a pan-cancer cohort of 9,460 patients. We did RNA-seq on 16 samples and verified that TMEscore was higher in patients with stage 3/4 than stage 1/2 diseases. The TMEscore could also predict responses for several immunotherapies. After adding clinical features, we found that the nomogram-based score system, the TMEIndex, surpassed the current risk system at predicting survivals. Our analysis explained TME at the transcriptome level and paved the way for immunotherapies in NB.

The Effect of Hypoxia and Hypoxia-Associated Pathways in the Regulation of Antitumor Response: Friends or Foes?

Hypoxia is an environmental stressor that is instigated by low oxygen availability. It fuels the progression of solid tumors by driving tumor plasticity, heterogeneity, stemness and genomic instability. Hypoxia metabolically reprograms the tumor microenvironment (TME), adding insult to injury to the acidic, nutrient deprived and poorly vascularized conditions that act to dampen immune cell function. Through its impact on key cancer hallmarks and by creating a physical barrier conducive to tumor survival, hypoxia modulates tumor cell escape from the mounted immune response. The tumor cell-immune cell crosstalk in the context of a hypoxic TME tips the balance towards a cold and immunosuppressed microenvironment that is resistant to immune checkpoint inhibitors (ICI). Nonetheless, evidence is emerging that could make hypoxia an asset for improving response to ICI. Tackling the tumor immune contexture has taken on an in silico, digitalized approach with an increasing number of studies applying bioinformatics to deconvolute the cellular and non-cellular elements of the TME. Such approaches have additionally been combined with signature-based proxies of hypoxia to further dissect the turbulent hypoxia-immune relationship. In this review we will be highlighting the mechanisms by which hypoxia impacts immune cell functions and how that could translate to predicting response to immunotherapy in an era of machine learning and computational biology.

Comprehensive evaluation of deconvolution methods for human brain gene expression

Transcriptome deconvolution aims to estimate the cellular composition of an RNA sample from its gene expression data, which in turn can be used to correct for composition differences across samples. The human brain is unique in its transcriptomic diversity, and comprises a complex mixture of cell-types, including transcriptionally similar subtypes of neurons. Here, we carry out a comprehensive evaluation of deconvolution methods for human brain transcriptome data, and assess the tissue-specificity of our key observations by comparison with human pancreas and heart. We evaluate eight transcriptome deconvolution approaches and nine cell-type signatures, testing the accuracy of deconvolution using in silico mixtures of single-cell RNA-seq data, RNA mixtures, as well as nearly 2000 human brain samples. Our results identify the main factors that drive deconvolution accuracy for brain data, and highlight the importance of biological factors influencing cell-type signatures, such as brain region and in vitro cell culturing.

Transcriptome deconvolution aims to estimate cellular composition based on gene expression data. Here the authors evaluate deconvolution methods for human brain transcriptome and conclude that partial deconvolution algorithms work best, but that appropriate cell-type signatures are also important.

A protocol to extract cell-type-specific signatures from differentially expressed genes in bulk-tissue RNA-seq

Bulk-tissue RNA-seq is widely used to dissect variation in gene expression levels across tissues and under different experimental conditions. Here, we introduce a protocol that leverages existing single-cell expression data to deconvolve patterns of cell-type-specific gene expression in differentially expressed gene lists from highly heterogeneous tissue. We apply this protocol to interrogate cell-type-specific gene expression and variation in cell type composition between the distinct sublayers of the hippocampal CA1 region of the brain in a rodent model of epilepsy.

For complete details on the use and execution of this protocol, please refer to Cid et al. (2021).

•

A protocol to explore gene signatures from bulk RNA-seq at the cell-type-specific level

•

Deconvolution of complex gene signatures from highly heterogeneous tissues

•

Publicly available single-cell gene expression dataset is retrieved and curated

•

Gene signatures across brain regions and disease states are surveyed in scRNA-seq data