Deep learning-based cell composition analysis from tissue expression profiles

ARTdeConv: adaptive regularized tri-factor non-negative matrix factorization for cell type deconvolution

Accurate deconvolution of cell types from bulk gene expression is crucial for understanding cellular compositions and uncovering cell-type specific differential expression and physiological states of diseased tissues. Existing deconvolution methods have limitations, such as requiring complete cellular gene expression signatures or neglecting partial biological information. Moreover, these methods often overlook varying cell-type messenger RNA amounts, leading to biased proportion estimates. Additionally, they do not effectively utilize valuable reference information from external studies, such as means and ranges of population cell-type proportions. To address these challenges, we introduce an adaptive regularized tri-factor non-negative matrix factorization approach for deconvolution (ARTdeConv). We rigorously establish the numerical convergence of our algorithm. Through benchmark simulations, we demonstrate the superior performance of ARTdeConv compared to state-of-the-art semi-reference-based and reference-free methods as well as its robustness under challenges to its assumptions. In a real-world application to a dataset from a trivalent influenza vaccine study, our method accurately estimates cellular proportions, as evidenced by the nearly perfect Pearson's correlation between ARTdeConv estimates and flow cytometry measurements. Moreover, our analysis of ARTdeConv estimates in COVID-19 patients reveals patterns consistent with important immunological phenomena observed in other studies. The proposed method, ARTdeConv, is implemented as an R package and can be accessed on GitHub for researchers and practitioners.

DeepDeconUQ estimates malignant cell fraction prediction intervals in bulk RNA-seq tissue

Accurate estimation of malignant cell fractions in tissues plays a critical role in cancer diagnosis, prognosis, and subsequent treatment decisions. However, most currently available methods provide only point estimates, neglecting the quantification of uncertainties, which is essential for both clinical and research applications. This study introduces DeepDeconUQ, a deep neural network model developed to estimate prediction intervals for malignant cell fractions based on bulk RNA-seq data. This approach addresses limitations in current malignant cell fraction estimation methods by integrating uncertainty quantification into predictions of cancer cell fractions. DeepDeconUQ leverages single-cell RNA sequencing (scRNA-seq) data in conjunction with conformalized quantile regression to produce reliable prediction intervals. The model trains a quantile regression neural network to establish upper and lower bounds for cancer cell proportions, followed by a calibration step that refines these intervals to ensure both statistical validity (coverage probability) and discrimination (narrow intervals). Benchmark analyses indicate that DeepDeconUQ consistently surpasses existing methods, achieving high coverage accuracy with tight prediction intervals across simulated and real cancer datasets. The robustness of DeepDeconUQ is further demonstrated by its resilience to various gene expression perturbations. The DeepDeconUQ method is publicly accessible at https://github.com/jiaweih14/DeepDeconUQ.

Microdroplet-Engineered Skeletal Muscle Organoids from Primary Tissue Recapitulate Parental Physiology with High Reproducibility

Achieving high maturity and functionality in in vitro skeletal muscle models is essential for advancing our understanding of muscle biology, disease mechanisms, and drug discovery. However, current models struggle to fully recapitulate key features such as sarcomere structure, muscle fiber composition, and contractile function while also ensuring consistency and rapid production. Adult stem cells residing in muscle tissue are known for their powerful regenerative potential, yet tissue-derived skeletal muscle organoids have not been established. In this study, we introduce droplet-engineered skeletal muscle organoids derived from primary tissue using cascade-tubing microfluidics. These droplet-engineered organoids (DEOs) exhibit high maturity, including well-developed striated sarcomeres, spontaneous and stimulated contractions, and recapitulation of parental muscle fiber types. Notably, DEOs are produced in just 8 d without the need for primary cell culture-substantially accelerating the 50- to 60-d process required by classical organoid models. Additionally, the cascade-tubing microfluidics platform enables high-throughput production of hundreds of uniform DEO replicates from a small tissue sample, providing a scalable and reproducible solution for skeletal muscle research and drug screening.

Unveiling the immune microenvironment of complex tissues and tumors in transcriptomics through a deconvolution approach

Accurately resolving the composition of tumor-infiltrating leukocytes is pivotal for advancing cancer immunotherapy strategies. Despite the success of some clinical trials, applying these strategies remains limited due to the challenges in deciphering the immune microenvironment. In this study, we developed a streamlined, two-step workflow to address the complexity of bioinformatics processes involved in analyzing immune cell composition from transcriptomics data. Our dockerized toolkit, DOCexpress_fastqc, integrates the hisat2-stringtie pipeline with customized scripts within Galaxy/Docker environments, facilitating RNA sequencing (RNA-seq) gene expression profiling. The output from DOCexpress_fastqc is seamlessly formatted with mySORT, a web application that employs a deconvolution algorithm to determine the immune content across 21 cell subclasses. We validated mySORT using synthetic pseudo-bulk data derived from single-cell RNA sequencing (scRNA-seq) datasets. Our predictions exhibit strong concordance with the ground-truth immune cell composition, achieving Pearson's correlation coefficients of 0.871 in melanoma patients and 0.775 in head and neck cancer patients. Additionally, mySORT outperforms existing methods like CIBERSORT in accuracy and provides a wide range of data visualization features, such as hierarchical clustering and cell complexity plots. The toolkit and web application are freely available for the research community, providing enhanced resolution for conventional bulk RNA sequencing data and facilitating the analysis of immune microenvironment responses in immunotherapy. The mySORT demo website and Docker image are free at https://mysort.iis.sinica.edu.tw and https://hub.docker.com/r/lsbnb/mysort_2022 .

Optimized RNA sequencing deconvolution illustrates the impact of obesity and weight loss on cell composition of human adipose tissue

OBJECTIVE

Cellular heterogeneity of human adipose tissue is linked to the pathophysiology of obesity and may impact the response to energy restriction and changes in fat mass. Herein, we provide an optimized pipeline to estimate cellular composition in human abdominal subcutaneous adipose tissue (ASAT) bulk RNA sequencing (RNA-seq) datasets using a single-nuclei RNA-seq signature matrix.

METHODS

A deconvolution pipeline for ASAT was optimized by benchmarking publicly available algorithms using a signature matrix derived from ASAT single-nuclei RNA-seq data from 20 adults and then applied to estimate ASAT cell-type proportions in publicly available obesity and weight loss studies.

RESULTS

Individuals with obesity had greater proportions of macrophages and lower proportions of adipocyte subpopulations and vascular cells compared with lean individuals. Two months of diet-induced weight loss increased the estimated proportions of macrophages; however, 2 years of diet-induced weight loss reduced the estimated proportions of macrophages, thereby suggesting a biphasic nature of cellular remodeling of ASAT during weight loss.

CONCLUSIONS

Our optimized high-throughput pipeline facilitates the assessment of composition changes of highly characterized cell types in large numbers of ASAT samples using low-cost bulk RNA-seq. Our data reveal novel changes in cellular heterogeneity and its association with cardiometabolic health in humans with obesity and following weight loss.

Advances and applications in single-cell and spatial genomics

The applications of single-cell and spatial technologies in recent times have revolutionized the present understanding of cellular states and the cellular heterogeneity inherent in complex biological systems. These advancements offer unprecedented resolution in the examination of the functional genomics of individual cells and their spatial context within tissues. In this review, we have comprehensively discussed the historical development and recent progress in the field of single-cell and spatial genomics. We have reviewed the breakthroughs in single-cell multi-omics technologies, spatial genomics methods, and the computational strategies employed toward the analyses of single-cell atlas data. Furthermore, we have highlighted the advances made in constructing cellular atlases and their clinical applications, particularly in the context of disease. Finally, we have discussed the emerging trends, challenges, and opportunities in this rapidly evolving field.

Comprehensive evaluation of methods for identifying tissues or cell types of origin of the plasma cell-free transcriptome

Plasma cell-free RNA (cfRNA) is derived from cells in various tissues and organs throughout the body and reflects the physiological and pathological conditions. Identifying the origins of cfRNA is essential for comprehending its variations. Only a few tools are designed for cfRNA deconvolution, and most studies have relied on traditional bulk RNA methods. In this study, we employed human tissue and cell transcriptomic data as reference sets and evaluated the performance of seven deconvolution methods on cfRNA. We compared the analysis results of cell types and tissues of origin of plasma cfRNA and chose to use single-cell RNA sequencing (scRNA-seq) data as reference to conduct further evaluation of deconvolution methods. Subsequently, we assessed the accuracy and robustness of the methods by utilizing simulated cfRNA data generated from scRNA-seq. We also evaluated the methods' accuracy on real plasma cfRNA data by analyzing the correlation between the predicted cell proportions and the corresponding clinical indicators. Moreover, we compared the methods' effectiveness in revealing the impacts of diseases on cells and evaluated the performance of cancer classification models based on the cell origin data they provided. In summary, our study provides valuable insights into cfRNA origin analysis, enhancing its potential in biomedical research.

Astrocytic RNA editing regulates the host immune response to alpha-synuclein

RNA editing is a posttranscriptional mechanism that targets changes in RNA transcripts to modulate innate immune responses. We report the role of astrocyte-specific, ADAR1-mediated RNA editing in neuroinflammation in Parkinson's disease (PD). We generated human induced pluripotent stem cell-derived astrocytes, neurons and cocultures and exposed them to small soluble alpha-synuclein aggregates. Oligomeric alpha-synuclein triggered an inflammatory glial state associated with Toll-like receptor activation, viral responses, and cytokine secretion. This reactive state resulted in loss of neurosupportive functions and the induction of neuronal toxicity. Notably, interferon response pathways were activated leading to up-regulation and isoform switching of the RNA deaminase enzyme, ADAR1. ADAR1 mediates A-to-I RNA editing, and increases in RNA editing were observed in inflammatory pathways in cells, as well as in postmortem human PD brain. Aberrant, or dysregulated, ADAR1 responses and RNA editing may lead to sustained inflammatory reactive states in astrocytes triggered by alpha-synuclein aggregation, and this may drive the neuroinflammatory cascade in Parkinson's.

Crohn's Disease and Ulcerative Colitis Share 2 Molecular Subtypes With Different Mechanisms and Drug Responses

BACKGROUND AND AIMS

Several therapies have been approved to treat Crohn's disease (CD) and ulcerative colitis (UC), indicating that both diseases may share the same molecular subtypes. The aim of this study is to identify shared patient subtypes with common molecular drivers of disease.

METHODS

Five public datasets with 406 CD and 421 UC samples were integrated to identify molecular subtypes. Then, the patient labels from 6 independent datasets and 8 treatment datasets were predicted for validating subtypes and identifying the relationship with response status of corticosteroids, infliximab, vedolizumab, and ustekinumab.

RESULTS

Two molecular subtypes were identified from the training datasets, in which CD and UC patients were relatively evenly represented in each subtype. We found 6 S1-specific gene modules related to innate/adaptive immune responses and tissue remodeling and 9 S1-specific cell types (cycling T cells, Tregs, CD8+ lamina propria, follicular B cells, cycling B cells, plasma cells, inflammatory monocytes, inflammatory fibroblasts, and postcapillary venules). Subtype S2 was associated with 3 modules related to metabolism functions and 4 cell types (immature enterocytes, transit amplifying cells, immature goblet cells, and WNT5B+ cells). The subtypes can be replicated in 6 independent datasets based on a 20-gene classifier. Furthermore, response rates to 4 treatments in subtype S2 were significantly higher than those in subtype S1.

CONCLUSIONS

This study discovered and validated a robust transcriptome-based molecular classification shared by CD and UC and built a 20-gene classifier. Because 2 subtypes have different molecular mechanisms and drug response, our classification may aid interpretation of heterogeneous molecular and clinical information in inflammatory bowel disease patients.

GBMPurity: A Machine Learning Tool for Estimating Glioblastoma Tumour Purity from Bulk RNA-seq Data

BACKGROUND

Glioblastoma (GBM) presents a significant clinical challenge due to its aggressive nature and extensive heterogeneity. Tumour purity, the proportion of malignant cells within a tumour, is an important covariate for understanding the disease, having direct clinical relevance or obscuring signal of the malignant portion in molecular analyses of bulk samples. However, current methods for estimating tumour purity are non-specific and technically demanding. Therefore, we aimed to build a reliable and accessible purity estimator for GBM.

METHODS

We developed GBMPurity, a deep-learning model specifically designed to estimate the purity of IDH-wildtype primary GBM from bulk RNA-seq data. The model was trained using simulated pseudobulk tumours of known purity from labelled single-cell data acquired from the GBmap resource. The performance of GBMPurity was evaluated and compared to several existing tools using independent datasets.

RESULTS

GBMPurity outperformed existing tools, achieving a mean absolute error of 0.15 and a concordance correlation coefficient of 0.88 on validation datasets. We demonstrate the utility of GBMPurity through inference on bulk RNA-seq samples and observe reduced purity of the Proneural molecular subtype relative to the Classical, attributed to the increased presence of healthy brain cells.

CONCLUSIONS

GBMPurity provides a reliable and accessible tool for estimating tumour purity from bulk RNA-seq data, enhancing the interpretation of bulk RNA-seq data and offering valuable insights into GBM biology. To facilitate the use of this model by the wider research community, GBMPurity is available as a web-based tool at: https://gbmdeconvoluter.leeds.ac.uk/.

A deconvolution framework that uses single-cell sequencing plus a small benchmark data set for accurate analysis of cell type ratios in complex tissue samples

Bulk deconvolution with single-cell/nucleus RNA-seq data is critical for understanding heterogeneity in complex biological samples, yet the technological discrepancy across sequencing platforms limits deconvolution accuracy. To address this, we utilize an experimental design to match inter-platform biological signals, hence revealing the technological discrepancy, and then develop a deconvolution framework called DeMixSC using this well-matched, that is, benchmark, data. Built upon a novel weighted nonnegative least-squares framework, DeMixSC identifies and adjusts genes with high technological discrepancy and aligns the benchmark data with large patient cohorts of matched-tissue-type for large-scale deconvolution. Our results using two benchmark data sets of healthy retinas and ovarian cancer tissues suggest much-improved deconvolution accuracy. Leveraging tissue-specific benchmark data sets, we applied DeMixSC to a large cohort of 453 age-related macular degeneration patients and a cohort of 30 ovarian cancer patients with various responses to neoadjuvant chemotherapy. Only DeMixSC successfully unveiled biologically meaningful differences across patient groups, demonstrating its broad applicability in diverse real-world clinical scenarios. Our findings reveal the impact of technological discrepancy on deconvolution performance and underscore the importance of a well-matched data set to resolve this challenge. The developed DeMixSC framework is generally applicable for accurately deconvolving large cohorts of disease tissues, including cancers, when a well-matched benchmark data set is available.

Senescence-related genes are associated with the immunopathology signature of American tegumentary leishmaniasis lesions and may predict progression to mucosal leishmaniasis

The American tegumentary leishmaniasis (ATL) is caused by protozoans of the genus Leishmania and varies from mild localized cutaneous leishmaniasis (LCL) form to more severe manifestations such as the diffuse cutaneous leishmaniasis (DCL) form and the mucosal leishmaniasis (ML) form. Previously, we demonstrated the accumulation of senescent cells in skin lesions of patients with LCL. Moreover, lesional transcriptomic analyses revealed a robust co-induction of senescence and pro-inflammatory gene signatures, highlighting the critical role of senescent T cells in orchestrating pathology. In this work we hypothesized that senescent cells might operate differently among the ATL spectrum, potentially influencing immunopathological mechanisms and clinical outcome. We analysed previously published RNA-Seq datasets of skin biopsies of healthy subjects and lesional skin from DCL patients, LCL patients, and LCL patients that, after treatment, progressed to mucosal leishmaniasis (MLP). Our findings demonstrate a robust presence of a CD8 T-cell signature associated with both LCL and MLP lesions. Moreover, both inflammatory and cytotoxic signatures were significantly upregulated, showing a strong increase in MLP and LCL groups, but not DCL. The senescence signature was elevated between LCL and MLP groups, representing the only distinguishable signature of immunopathology between them. Interestingly, our analyses further revealed the senescence signature's capacity to predict progression from LCL to mucosal forms, which was not observed with other signatures. Both the senescence-signature score and specific senescence-associated genes demonstrated an increased capacity to predict mucosal progression, with correct predictions exceeding 97% of cases. Collectively, our findings contribute to a comprehensive understanding of immunosenescence in ATL and suggest that senescence may represent the latest and most important signature of the immunopathogenisis. This highlights its potential value in predicting disease severity.

BuDDI: Bulk Deconvolution with Domain Invariance to predict cell-type-specific perturbations from bulk

While single-cell experiments provide deep cellular resolution within a single sample, some single-cell experiments are inherently more challenging than bulk experiments due to dissociation difficulties, cost, or limited tissue availability. This creates a situation where we have deep cellular profiles of one sample or condition, and bulk profiles across multiple samples and conditions. To bridge this gap, we propose BuDDI (BUlk Deconvolution with Domain Invariance). BuDDI utilizes domain adaptation techniques to effectively integrate available corpora of case-control bulk and reference scRNA-seq observations to infer cell-type-specific perturbation effects. BuDDI achieves this by learning independent latent spaces within a single variational autoencoder (VAE) encompassing at least four sources of variability: 1) cell type proportion, 2) perturbation effect, 3) structured experimental variability, and 4) remaining variability. Since each latent space is encouraged to be independent, we simulate perturbation responses by independently composing each latent space to simulate cell-type-specific perturbation responses. We evaluated BuDDI's performance on simulated and real data with experimental designs of increasing complexity. We first validated that BuDDI could learn domain invariant latent spaces on data with matched samples across each source of variability. Then we validated that BuDDI could accurately predict cell-type-specific perturbation response when no single-cell perturbed profiles were used during training; instead, only bulk samples had both perturbed and non-perturbed observations. Finally, we validated BuDDI on predicting sex-specific differences, an experimental design where it is not possible to have matched samples. In each experiment, BuDDI outperformed all other comparative methods and baselines. As more reference atlases are completed, BuDDI provides a path to combine these resources with bulk-profiled treatment or disease signatures to study perturbations, sex differences, or other factors at single-cell resolution.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

SurfR: Riding the wave of RNA-seq data with a comprehensive bioconductor package to identify surface protein-coding genes

Motivation

Proteins at the cell surface connect signaling networks and largely determine a cell's capacity to communicate and interact with its environment. In particular, variations in transcriptomic profiles are often observed between healthy and diseased cells, leading to distinct sets of cell-surface proteins. For these reasons, cell-surface proteins may act as biomarkers for the detection of cells of interest in tissues or body fluids, are often the target of pharmaceutical agents, and hold significant promise in the clinical practice for diagnosis, prognosis, treatment development, and evaluation of therapy response. Therefore, implementing robust methods to identify condition-specific cell-surface proteins is of pivotal importance to advance biomedical research.

Results

We developed SurfR, an R/Bioconductor package providing a streamlined end-to-end workflow for computationally identifying surface protein-coding genes from expression data. Our user-friendly, comprehensive workflow performs systematic expression data retrieval from public databases, differential gene expression across conditions, integration of datasets, enrichment analysis, identification of targetable proteins on a condition of interest, and data visualization.

Availability and implementation

SurfR is released under GNU-GPL-v3.0 License. Source code, documentation, examples, and tutorials are available through Bioconductor (http://www.bioconductor.org/packages/SurfR). RMD notebooks with the use cases code described in the manuscript can be found on GitHub (https://github.com/auroramaurizio/SurfR_UseCases).

DeSide: A unified deep learning approach for cellular deconvolution of tumor microenvironment

Cellular deconvolution via bulk RNA sequencing (RNA-seq) presents a cost-effective and efficient alternative to experimental methods such as flow cytometry and single-cell RNA-seq (scRNA-seq) for analyzing the complex cellular composition of tumor microenvironments. Despite challenges due to heterogeneity within and among tumors, our innovative deep learning-based approach, DeSide, shows exceptional accuracy in estimating the proportions of 16 distinct cell types and subtypes within solid tumors. DeSide integrates biological pathways and assesses noncancerous cell types first, effectively sidestepping the issue of highly variable gene expression profiles (GEPs) associated with cancer cells. By leveraging scRNA-seq data from six cancer types and 185 cancer cell lines across 22 cancer types as references, our method introduces distinctive sampling and filtering techniques to generate a high-quality training set that closely replicates real tumor GEPs, based on The Cancer Genome Atlas (TCGA) bulk RNA-seq data. With this model and high-quality training set, DeSide outperforms existing methods in estimating tumor purity and the proportions of noncancerous cells within solid tumors. Our model precisely predicts cellular compositions across 19 cancer types from TCGA and proves its effectiveness with multiple additional external datasets. Crucially, DeSide enables the identification and analysis of combinatorial cell type pairs, facilitating the stratification of cancer patients into prognostically significant groups. This approach not only provides deeper insights into the dynamics of tumor biology but also highlights potential therapeutic targets by underscoring the importance of specific cell type or subtype interactions.

Droplet-engineered organoids recapitulate parental tissue transcriptome with inter-organoid homogeneity and inter-tumor cell heterogeneity

Organoids are expected to function as effective human organ models for precision cancer studies and drug development. Currently, primary tissue-derived organoids, termed non-engineered organoids (NEOs), are produced by manual pipetting or liquid handling that compromises organoid-organoid homogeneity and organoid-tissue consistency. Droplet-based microfluidics enables automated organoid production with high organoid-organoid homogeneity, organoid-tissue consistency, and a significantly improved production spectrum. It takes advantage of droplet-encapsulation of defined populations of cells and droplet-rendered microstructures that guide cell self-organization. Herein, we studied the droplet-engineered organoids (DEOs), derived from mouse liver tissues and human liver tumors, by using transcriptional analysis and cellular deconvolution on bulk RNA-seq data. The characteristics of DEOs are compared with the parental liver tissues (or tumors) and NEOs. The DEOs are proven higher reproducibility and consistency with the parental tissues, have a high production spectrum and shortened modeling time, and possess inter-organoid homogeneity and inter-tumor cell heterogeneity.

Recent Advances in Artificial Intelligence to Improve Immunotherapy and the Use of Digital Twins to Identify Prognosis of Patients with Solid Tumors

To date, the public health system has been impacted by the increasing costs of many diagnostic and therapeutic pathways due to limited resources. At the same time, we are constantly seeking to improve these paths through approaches aimed at personalized medicine. To achieve the required levels of diagnostic and therapeutic precision, it is necessary to integrate data from different sources and simulation platforms. Today, artificial intelligence (AI), machine learning (ML), and predictive computer models are more efficient at guiding decisions regarding better therapies and medical procedures. The evolution of these multiparametric and multimodal systems has led to the creation of digital twins (DTs). The goal of our review is to summarize AI applications in discovering new immunotherapies and developing predictive models for more precise immunotherapeutic decision-making. The findings from this literature review highlight that DTs, particularly predictive mathematical models, will be pivotal in advancing healthcare outcomes. Over time, DTs will indeed bring the benefits of diagnostic precision and personalized treatment to a broader spectrum of patients.

Identification of nitric oxide-mediated necroptosis as the predominant death route in Parkinson's disease

Parkinson's disease (PD) involves multiple forms of neuronal cell death, but the dominant pathway involved in disease progression remains unclear. This study employed RNA sequencing (RNA-seq) of brain tissue to explore gene expression patterns across different stages of PD. Using the Scaden deep learning algorithm, we predicted neurocyte subtypes and modelled dynamic interactions for five classic cell death pathways to identify the predominant routes of neuronal death during PD progression. Our cell type-specific analysis revealed an increasing shift towards necroptosis, which was strongly correlated with nitric oxide synthase (NOS) expression across most neuronal subtypes. In vitro experiments confirmed that nitric oxide (NO) is a key mediator of necroptosis, leading to nuclear shrinkage and decreased mitochondrial membrane potential via phosphorylation of the PIP1/PIP3/MLKL signalling cascade. Importantly, specific necroptosis inhibitors significantly mitigated neuronal damage in both in vitro and in vivo PD models. Further analysis revealed that NO-mediated necroptosis is prevalent in early-onset Alzheimer's disease (AD) and amyotrophic lateral sclerosis (ALS) and across multiple brain regions but not in brain tumours. Our findings suggest that NO-mediated necroptosis is a critical pathway in PD and other neurodegenerative disorders, providing potential targets for therapeutic intervention.

Lesional senescent CD4+ T cells mediate bystander cytolysis and contribute to the skin pathology of human cutaneous leishmaniasis

Cytotoxic activity is a hallmark of the immunopathogenesis in human cutaneous leishmaniasis (CL). In this study, we identified accumulation of CD4+ granzyme B producing T cells with increased cytotoxic capacity in CL lesions. These cells showed enhanced expression of activating NK receptors (NKG2D and NKG2C), diminished expression of inhibitory NKG2A, along with the upregulation of the senescence marker CD57. Notably, CD4+ T cells freshly isolated from CL lesions demonstrated remarkable capacity to mediate NL-like bystander cytolysis. Phenotypic analyses revealed that lesional CD4+ T cells are mainly composed of late-differentiated effector (CD27-CD45RA-) and terminally differentiated (senescent) TEMRA (CD27-CD45RA+) subsets. Interestingly, the TEMRA CD4+ T cells exhibited higher expression of granzyme B and CD107a. Collectively, our results provide the first evidence that senescent cytotoxic CD4+ T cells may support the skin pathology of human cutaneous leishmaniasis and, together with our previous findings, support the notion that multiple subsets of cytotoxic senescent cells may be involved in inducing the skin lesions in these patients.

Cross-tissue human fibroblast atlas reveals myofibroblast subtypes with distinct roles in immune modulation

Fibroblasts, known for their functional diversity, play crucial roles in inflammation and cancer. In this study, we conduct comprehensive single-cell RNA sequencing analyses on fibroblast cells from 517 human samples, spanning 11 tissue types and diverse pathological states. We identify distinct fibroblast subpopulations with universal and tissue-specific characteristics. Pathological conditions lead to significant shifts in fibroblast compositions, including the expansion of immune-modulating fibroblasts during inflammation and tissue-remodeling myofibroblasts in cancer. Within the myofibroblast category, we identify four transcriptionally distinct subpopulations originating from different developmental origins, with LRRC15+ myofibroblasts displaying terminally differentiated features. Both LRRC15+ and MMP1+ myofibroblasts demonstrate pro-tumor potential that contribute to the immune-excluded and immune-suppressive tumor microenvironments (TMEs), whereas PI16+ fibroblasts show potential anti-tumor functions in adjacent non-cancerous regions. Fibroblast-subtype compositions define patient subtypes with distinct clinical outcomes. This study advances our understanding of fibroblast biology and suggests potential therapeutic strategies for targeting specific fibroblast subsets in cancer treatment.

Artificial intelligence and machine learning applications for cultured meat

Cultured meat has the potential to provide a complementary meat industry with reduced environmental, ethical, and health impacts. However, major technological challenges remain which require time-and resource-intensive research and development efforts. Machine learning has the potential to accelerate cultured meat technology by streamlining experiments, predicting optimal results, and reducing experimentation time and resources. However, the use of machine learning in cultured meat is in its infancy. This review covers the work available to date on the use of machine learning in cultured meat and explores future possibilities. We address four major areas of cultured meat research and development: establishing cell lines, cell culture media design, microscopy and image analysis, and bioprocessing and food processing optimization. In addition, we have included a survey of datasets relevant to CM research. This review aims to provide the foundation necessary for both cultured meat and machine learning scientists to identify research opportunities at the intersection between cultured meat and machine learning.

Spatial transcriptomics elucidates medulla niche supporting germinal center response in myasthenia gravis-associated thymoma

Myasthenia gravis (MG) is etiologically associated with thymus abnormalities, but its pathology in the thymus remains unclear. In this study, we attempt to narrow down the features associated with MG using spatial transcriptome analysis of thymoma and thymic hyperplasia samples. We find that the majority of thymomas are constituted by the cortical region. However, the small medullary region is enlarged in seropositive thymomas and contains polygenic enrichment and MG-specific germinal center structures. Neuromuscular medullary thymic epithelial cells, previously identified as MG-specific autoantigen-producing cells, are enriched in the cortico-medullary junction. The medulla is characterized by a specific chemokine pattern and immune cell composition, including migratory dendritic cells and effector regulatory T cells. Similar germinal center structures and immune microenvironments are also observed in the thymic hyperplasia medulla. This study shows that the medulla and junction areas are linked to MG pathology and provides insights into future MG research.

Machine learning and related approaches in transcriptomics

Data acquisition for transcriptomic studies used to be the bottleneck in the transcriptomic analytical pipeline. However, recent developments in transcriptome profiling technologies have increased researchers' ability to obtain data, resulting in a shift in focus to data analysis. Incorporating machine learning to traditional analytical methods allows the possibility of handling larger volumes of complex data more efficiently. Many bioinformaticians, especially those unfamiliar with ML in the study of human transcriptomics and complex biological systems, face a significant barrier stemming from their limited awareness of the current landscape of ML utilisation in this field. To address this gap, this review endeavours to introduce those individuals to the general types of ML, followed by a comprehensive range of more specific techniques, demonstrated through examples of their incorporation into analytical pipelines for human transcriptome investigations. Important computational aspects such as data pre-processing, task formulation, results (performance of ML models), and validation methods are encompassed. In hope of better practical relevance, there is a strong focus on studies published within the last five years, almost exclusively examining human transcriptomes, with outcomes compared with standard non-ML tools.

Evaluating cell type deconvolution in FFPE breast tissue: application to benign breast disease

Transcriptome profiling using RNA sequencing (RNA-seq) of bulk formalin-fixed paraffin-embedded (FFPE) tissue blocks is a standard method in biomedical research. However, when used on tissues with diverse cell type compositions, it yields averaged gene expression profiles, complicating biomarker identification due to variations in cell proportions. To address the need for optimized strategies for defining individual cell type compositions from bulk FFPE samples, we constructed single-cell RNA-seq reference data for breast tissue and tested cell type deconvolution methods. Initial simulation experiments showed similar performances across multiple commonly used deconvolution methods. However, the introduction of FFPE artifacts significantly impacted their performances, with a root mean squared error (RMSE) ranging between 0.04 and 0.17. Scaden, a deep learning-based method, consistently outperformed the others, demonstrating robustness against FFPE artifacts. Testing these methods on our 62-sample RNA-seq benign breast disease cohort in which cell type composition was estimated using digital pathology approaches, we found that pre-filtering of the reference data enhanced the accuracy of most methods, realizing up to a 32% reduction in RMSE. To support further research efforts in this domain, we introduce SCdeconR, an R package designed for streamlined cell type deconvolution assessments and downstream analyses.

Pan-cancer integrative analyses dissect the remodeling of endothelial cells in human cancers

Therapeutics targeting tumor endothelial cells (TECs) have been explored for decades, with only suboptimal efficacy achieved, partly due to an insufficient understanding of the TEC heterogeneity across cancer patients. We integrated single-cell RNA-seq data of 575 cancer patients from 19 solid tumor types, comprehensively charting the TEC phenotypic diversities. Our analyses uncovered underappreciated compositional and functional heterogeneity in TECs from a pan-cancer perspective. Two subsets, CXCR4 + tip cells and SELE + veins, represented the prominent angiogenic and proinflammatory phenotypes of TECs, respectively. They exhibited distinct spatial organization patterns, and compared to adjacent non-tumor tissues, tumor tissue showed an increased prevalence of CXCR4 + tip cells, yet with SELE + veins depleted. Such functional and spatial characteristics underlie their differential associations with the response of anti-angiogenic therapies and immunotherapies. Our integrative resources and findings open new avenues to understand and clinically intervene in the tumor vasculature.

Community assessment of methods to deconvolve cellular composition from bulk gene expression

We evaluate deconvolution methods, which infer levels of immune infiltration from bulk expression of tumor samples, through a community-wide DREAM Challenge. We assess six published and 22 community-contributed methods using in vitro and in silico transcriptional profiles of admixed cancer and healthy immune cells. Several published methods predict most cell types well, though they either were not trained to evaluate all functional CD8+ T cell states or do so with low accuracy. Several community-contributed methods address this gap, including a deep learning-based approach, whose strong performance establishes the applicability of this paradigm to deconvolution. Despite being developed largely using immune cells from healthy tissues, deconvolution methods predict levels of tumor-derived immune cells well. Our admixed and purified transcriptional profiles will be a valuable resource for developing deconvolution methods, including in response to common challenges we observe across methods, such as sensitive identification of functional CD4+ T cell states.

sNucConv: A bulk RNA-seq deconvolution method trained on single-nucleus RNA-seq data to estimate cell-type composition of human adipose tissues

Deconvolution algorithms mostly rely on single-cell RNA-sequencing (scRNA-seq) data applied onto bulk RNA-sequencing (bulk RNA-seq) to estimate tissues' cell-type composition, with performance accuracy validated on deposited databases. Adipose tissues' cellular composition is highly variable, and adipocytes can only be captured by single-nucleus RNA-sequencing (snRNA-seq). Here we report the development of sNucConv, a Scaden deep-learning-based deconvolution tool, trained using 5 hSAT and 7 hVAT snRNA-seq-based data corrected by (i) snRNA-seq/bulk RNA-seq highly correlated genes and (ii) individual cell-type regression models. Applying sNucConv on our bulk RNA-seq data resulted in cell-type proportion estimation of 15 and 13 cell types, with accuracy of R = 0.93 (range: 0.76-0.97) and R = 0.95 (range: 0.92-0.98) for hVAT and hSAT, respectively. This performance level was further validated on an independent set of samples (5 hSAT; 5 hVAT). The resulting model was depot specific, reflecting depot differences in gene expression patterns. Jointly, sNucConv provides proof-of-concept for producing validated deconvolution models for tissues un-amenable to scRNA-seq.

scSemiProfiler: Advancing large-scale single-cell studies through semi-profiling with deep generative models and active learning

Single-cell sequencing is a crucial tool for dissecting the cellular intricacies of complex diseases. Its prohibitive cost, however, hampers its application in expansive biomedical studies. Traditional cellular deconvolution approaches can infer cell type proportions from more affordable bulk sequencing data, yet they fall short in providing the detailed resolution required for single-cell-level analyses. To overcome this challenge, we introduce "scSemiProfiler", an innovative computational framework that marries deep generative models with active learning strategies. This method adeptly infers single-cell profiles across large cohorts by fusing bulk sequencing data with targeted single-cell sequencing from a few rigorously chosen representatives. Extensive validation across heterogeneous datasets verifies the precision of our semi-profiling approach, aligning closely with true single-cell profiling data and empowering refined cellular analyses. Originally developed for extensive disease cohorts, "scSemiProfiler" is adaptable for broad applications. It provides a scalable, cost-effective solution for single-cell profiling, facilitating in-depth cellular investigation in various biological domains.

Unbiasedly decoding the tumor microenvironment with single-cell multiomics analysis in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive malignancy with a poor prognosis and limited therapeutic options. Research on the tumor microenvironment (TME) of PDAC has propelled the development of immunotherapeutic and targeted therapeutic strategies with a promising future. The emergence of single-cell sequencing and mass spectrometry technologies, coupled with spatial omics, has collectively revealed the heterogeneity of the TME from a multiomics perspective, outlined the development trajectories of cell lineages, and revealed important functions of previously underrated myeloid cells and tumor stroma cells. Concurrently, these findings necessitated more refined annotations of biological functions at the cell cluster or single-cell level. Precise identification of all cell clusters is urgently needed to determine whether they have been investigated adequately and to identify target cell clusters with antitumor potential, design compatible treatment strategies, and determine treatment resistance. Here, we summarize recent research on the PDAC TME at the single-cell multiomics level, with an unbiased focus on the functions and potential classification bases of every cellular component within the TME, and look forward to the prospects of integrating single-cell multiomics data and retrospectively reusing bulk sequencing data, hoping to provide new insights into the PDAC TME.

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.

Blood-based molecular and cellular biomarkers of early response to neoadjuvant PD-1 blockade in patients with non-small cell lung cancer

BACKGROUND

Despite the improved survival observed in PD-1/PD-L1 blockade therapy, a substantial proportion of cancer patients, including those with non-small cell lung cancer (NSCLC), still lack a response.

METHODS

Transcriptomic profiling was conducted on a discovery cohort comprising 100 whole blood samples, as collected multiple times from 48 healthy controls (including 43 published data) and 31 NSCLC patients that under treatment with a combination of anti-PD-1 Tislelizumab and chemotherapy. Differentially expressed genes (DEGs), simulated immune cell subsets, and germline DNA mutational markers were identified from patients achieved a pathological complete response during the early treatment cycles. The predictive values of mutational markers were further validated in an independent immunotherapy cohort of 1661 subjects, and then confirmed in genetically matched lung cancer cell lines by a co-culturing model.

RESULTS

The gene expression of hundreds of DEGs (FDR p < 0.05, fold change < -2 or > 2) distinguished responders from healthy controls, indicating the potential to stratify patients utilizing early on-treatment features from blood. PD-1-mediated cell abundance changes in memory CD4 + and regulatory T cell subset were more significant or exclusively observed in responders. A panel of top-ranked genetic alterations showed significant associations with improved survival (p < 0.05) and heightened responsiveness to anti-PD-1 treatment in patient cohort and co-cultured cell lines.

CONCLUSION

This study discovered and validated peripheral blood-based biomarkers with evident predictive efficacy for early therapy response and patient stratification before treatment for neoadjuvant PD-1 blockade in NSCLC patients.

Adaptive digital tissue deconvolution

MOTIVATION

The inference of cellular compositions from bulk and spatial transcriptomics data increasingly complements data analyses. Multiple computational approaches were suggested and recently, machine learning techniques were developed to systematically improve estimates. Such approaches allow to infer additional, less abundant cell types. However, they rely on training data which do not capture the full biological diversity encountered in transcriptomics analyses; data can contain cellular contributions not seen in the training data and as such, analyses can be biased or blurred. Thus, computational approaches have to deal with unknown, hidden contributions. Moreover, most methods are based on cellular archetypes which serve as a reference; e.g. a generic T-cell profile is used to infer the proportion of T-cells. It is well known that cells adapt their molecular phenotype to the environment and that pre-specified cell archetypes can distort the inference of cellular compositions.

RESULTS

We propose Adaptive Digital Tissue Deconvolution (ADTD) to estimate cellular proportions of pre-selected cell types together with possibly unknown and hidden background contributions. Moreover, ADTD adapts prototypic reference profiles to the molecular environment of the cells, which further resolves cell-type specific gene regulation from bulk transcriptomics data. We verify this in simulation studies and demonstrate that ADTD improves existing approaches in estimating cellular compositions. In an application to bulk transcriptomics data from breast cancer patients, we demonstrate that ADTD provides insights into cell-type specific molecular differences between breast cancer subtypes.

AVAILABILITY AND IMPLEMENTATION

A python implementation of ADTD and a tutorial are available at Gitlab and zenodo (doi:10.5281/zenodo.7548362).

NNICE: a deep quantile neural network algorithm for expression deconvolution

The composition of cell-type is a key indicator of health. Advancements in bulk gene expression data curation, single cell RNA-sequencing technologies, and computational deconvolution approaches offer a new perspective to learn about the composition of different cell types in a quick and affordable way. In this study, we developed a quantile regression and deep learning-based method called Neural Network Immune Contexture Estimator (NNICE) to estimate the cell type abundance and its uncertainty by automatically deconvolving bulk RNA-seq data. The proposed NNICE model was able to successfully recover ground-truth cell type fraction values given unseen bulk mixture gene expression profiles from the same dataset it was trained on. Compared with baseline methods, NNICE achieved better performance on deconvolve both pseudo-bulk gene expressions (Pearson correlation R = 0.9) and real bulk gene expression data (Pearson correlation R = 0.9) across all cell types. In conclusion, NNICE combines statistic inference with deep learning to provide accurate and interpretable cell type deconvolution from bulk gene expression.

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.

Multiomic ALS signatures highlight subclusters and sex differences suggesting the MAPK pathway as therapeutic target

Amyotrophic lateral sclerosis (ALS) is a debilitating motor neuron disease and lacks effective disease-modifying treatments. This study utilizes a comprehensive multiomic approach to investigate the early and sex-specific molecular mechanisms underlying ALS. By analyzing the prefrontal cortex of 51 patients with sporadic ALS and 50 control subjects, alongside four transgenic mouse models (C9orf72-, SOD1-, TDP-43-, and FUS-ALS), we have uncovered significant molecular alterations associated with the disease. Here, we show that males exhibit more pronounced changes in molecular pathways compared to females. Our integrated analysis of transcriptomes, (phospho)proteomes, and miRNAomes also identified distinct ALS subclusters in humans, characterized by variations in immune response, extracellular matrix composition, mitochondrial function, and RNA processing. The molecular signatures of human subclusters were reflected in specific mouse models. Our study highlighted the mitogen-activated protein kinase (MAPK) pathway as an early disease mechanism. We further demonstrate that trametinib, a MAPK inhibitor, has potential therapeutic benefits in vitro and in vivo, particularly in females, suggesting a direction for developing targeted ALS treatments.

Neural-net-based cell deconvolution from DNA methylation reveals tumor microenvironment associated with cancer prognosis

DNA methylation is a pivotal epigenetic modification that defines cellular identity. While cell deconvolution utilizing this information is considered useful for clinical practice, current methods for deconvolution are limited in their accuracy and resolution. In this study, we collected DNA methylation data from 945 human samples derived from various tissues and tumor-infiltrating immune cells and trained a neural network model with them. The model, termed MEnet, predicted abundance of cell population together with the detailed immune cell status from bulk DNA methylation data, and showed consistency to those of flow cytometry and histochemistry. MEnet was superior to the existing methods in the accuracy, speed, and detectable cell diversity, and could be applicable for peripheral blood, tumors, cell-free DNA, and formalin-fixed paraffin-embedded sections. Furthermore, by applying MEnet to 72 intrahepatic cholangiocarcinoma samples, we identified immune cell profiles associated with cancer prognosis. We believe that cell deconvolution by MEnet has the potential for use in clinical settings.

Fourteen years of cellular deconvolution: methodology, applications, technical evaluation and outstanding challenges

Single-cell RNA sequencing (scRNA-Seq) is a recent technology that allows for the measurement of the expression of all genes in each individual cell contained in a sample. Information at the single-cell level has been shown to be extremely useful in many areas. However, performing single-cell experiments is expensive. Although cellular deconvolution cannot provide the same comprehensive information as single-cell experiments, it can extract cell-type information from bulk RNA data, and therefore it allows researchers to conduct studies at cell-type resolution from existing bulk datasets. For these reasons, a great effort has been made to develop such methods for cellular deconvolution. The large number of methods available, the requirement of coding skills, inadequate documentation, and lack of performance assessment all make it extremely difficult for life scientists to choose a suitable method for their experiment. This paper aims to fill this gap by providing a comprehensive review of 53 deconvolution methods regarding their methodology, applications, performance, and outstanding challenges. More importantly, the article presents a benchmarking of all these 53 methods using 283 cell types from 30 tissues of 63 individuals. We also provide an R package named DeconBenchmark that allows readers to execute and benchmark the reviewed methods (https://github.com/tinnlab/DeconBenchmark).

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.

Multi-Omics Characterization of Colon Mucosa and Submucosa/Wall from Crohn's Disease Patients

Crohn's disease (CD) is a subtype of inflammatory bowel disease (IBD) characterized by transmural disease. The concept of transmural healing (TH) has been proposed as an indicator of deep clinical remission of CD and as a predictor of favorable treatment endpoints. Understanding the pathophysiology involved in transmural disease is critical to achieving these endpoints. However, most studies have focused on the intestinal mucosa, overlooking the contribution of the intestinal wall in Crohn's disease. Multi-omics approaches have provided new avenues for exploring the pathogenesis of Crohn's disease and identifying potential biomarkers. We aimed to use transcriptomic and proteomic technologies to compare immune and mesenchymal cell profiles and pathways in the mucosal and submucosa/wall compartments to better understand chronic refractory disease elements to achieve transmural healing. The results revealed similarities and differences in gene and protein expression profiles, metabolic mechanisms, and immune and non-immune pathways between these two compartments. Additionally, the identification of protein isoforms highlights the complex molecular mechanisms underlying this disease, such as decreased RTN4 isoforms (RTN4B2 and RTN4C) in the submucosa/wall, which may be related to the dysregulation of enteric neural processes. These findings have the potential to inform the development of novel therapeutic strategies to achieve TH.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

DISSECT: deep semi-supervised consistency regularization for accurate cell type fraction and gene expression estimation

Cell deconvolution is the estimation of cell type fractions and cell type-specific gene expression from mixed data. An unmet challenge in cell deconvolution is the scarcity of realistic training data and the domain shift often observed in synthetic training data. Here, we show that two novel deep neural networks with simultaneous consistency regularization of the target and training domains significantly improve deconvolution performance. Our algorithm, DISSECT, outperforms competing algorithms in cell fraction and gene expression estimation by up to 14 percentage points. DISSECT can be easily adapted to other biomedical data types, as exemplified by our proteomic deconvolution experiments.

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

scRNA-seq characterizing the heterogeneity of fibroblasts in breast cancer reveals a novel subtype SFRP4+ CAF that inhibits migration and predicts prognosis

Introduction

Cancer-associated fibroblasts (CAFs) are a diverse group of cells that significantly impact the tumor microenvironment and therapeutic responses in breast cancer (BC). Despite their importance, the comprehensive profile of CAFs in BC remains to be fully elucidated.

Methods

To address this gap, we utilized single-cell RNA sequencing (scRNA-seq) to delineate the CAF landscape within 14 BC normal-tumor paired samples. We further corroborated our findings by analyzing several public datasets, thereby validating the newly identified CAF subtype. Additionally, we conducted coculture experiments with BC cells to assess the functional implications of this CAF subtype.

Results

Our scRNA-seq analysis unveiled eight distinct CAF subtypes across five tumor and six adjacent normal tissue samples. Notably, we discovered a novel subtype, designated as SFRP4+ CAFs, which was predominantly observed in normal tissues. The presence of SFRP4+ CAFs was substantiated by two independent scRNA-seq datasets and a spatial transcriptomics dataset. Functionally, SFRP4+ CAFs were found to impede BC cell migration and the epithelial-mesenchymal transition (EMT) process by secreting SFRP4, thereby modulating the WNT signaling pathway. Furthermore, we established that elevated expression levels of SFRP4+ CAF markers correlate with improved survival outcomes in BC patients, yet paradoxically, they predict a diminished response to neoadjuvant chemotherapy in cases of triple-negative breast cancer.

Conclusion

This investigation sheds light on the heterogeneity of CAFs in BC and introduces a novel SFRP4+ CAF subtype that hinders BC cell migration. This discovery holds promise as a potential biomarker for refined prognostic assessment and therapeutic intervention in BC.

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.

BuDDI: BulkDeconvolution withDomainInvariance to predict cell-type-specific perturbations from bulk

While single-cell experiments provide deep cellular resolution within a single sample, some single-cell experiments are inherently more challenging than bulk experiments due to dissociation difficulties, cost, or limited tissue availability. This creates a situation where we have deep cellular profiles of one sample or condition, and bulk profiles across multiple samples and conditions. To bridge this gap, we propose BuDDI (BUlk Deconvolution with Domain Invariance). BuDDI utilizes domain adaptation techniques to effectively integrate available corpora of case-control bulk and reference scRNA-seq observations to infer cell-type-specific perturbation effects. BuDDI achieves this by learning independent latent spaces within a single variational autoencoder (VAE) encompassing at least four sources of variability: 1) cell type proportion, 2) perturbation effect, 3) structured experimental variability, and 4) remaining variability. Since each latent space is encouraged to be independent, we simulate perturbation responses by independently composing each latent space to simulate cell-type-specific perturbation responses. We evaluated BuDDI's performance on simulated and real data with experimental designs of increasing complexity. We first validated that BuDDI could learn domain invariant latent spaces on data with matched samples across each source of variability. Then we validated that BuDDI could accurately predict cell-type-specific perturbation response when no single-cell perturbed profiles were used during training; instead, only bulk samples had both perturbed and non-perturbed observations. Finally, we validated BuDDI on predicting sex-specific differences, an experimental design where it is not possible to have matched samples. In each experiment, BuDDI outperformed all other comparative methods and baselines. As more reference atlases are completed, BuDDI provides a path to combine these resources with bulk-profiled treatment or disease signatures to study perturbations, sex differences, or other factors at single-cell resolution.

IBDTransDB: a manually curated transcriptomic database for inflammatory bowel disease

Inflammatory Bowel Disease (IBD) therapies are ineffective in at least 40% patients, and transcriptomic datasets have been widely used to reveal the pathogenesis and to identify the novel drug targets for these patients. Although public IBD transcriptomic datasets are available from many web-based tools/databases, due to the unstructured metadata and data description of these public datasets, most of these tools/databases do not allow querying datasets based on multiple keywords (e.g. colon and infliximab). Furthermore, few tools/databases can compare and integrate the datasets from the query results. To fill these gaps, we have developed IBDTransDB (https://abbviegrc.shinyapps.io/ibdtransdb/), a manually curated transcriptomic database for IBD. IBDTransDB includes a manually curated database with 34 transcriptomic datasets (2932 samples, 122 differential comparisons) and a query system supporting 35 keywords from 5 attributes (e.g. tissue and treatment). IBDTransDB also provides three modules for data analyses and integration. IBDExplore allows interactive visualization of differential gene list, pathway enrichment, gene signature and cell deconvolution analyses from a single dataset. IBDCompare supports comparisons of selected genes or pathways from multiple datasets across different conditions. IBDIntegrate performs meta-analysis to prioritize a list of genes/pathways based on user-selected datasets and conditions. Using two case studies related to infliximab treatment, we demonstrated that IBDTransDB provides a unique platform for biologists and clinicians to reveal IBD pathogenesis and identify the novel targets by integrating with other omics data. Database URL: https://abbviegrc.shinyapps.io/ibdtransdb/.

CATD: a reproducible pipeline for selecting cell-type deconvolution methods across tissues

Motivation

Cell-type deconvolution methods aim to infer cell composition from bulk transcriptomic data. The proliferation of developed methods coupled with inconsistent results obtained in many cases, highlights the pressing need for guidance in the selection of appropriate methods. Additionally, the growing accessibility of single-cell RNA sequencing datasets, often accompanied by bulk expression from related samples enable the benchmark of existing methods.

Results

In this study, we conduct a comprehensive assessment of 31 methods, utilizing single-cell RNA-sequencing data from diverse human and mouse tissues. Employing various simulation scenarios, we reveal the efficacy of regression-based deconvolution methods, highlighting their sensitivity to reference choices. We investigate the impact of bulk-reference differences, incorporating variables such as sample, study and technology. We provide validation using a gold standard dataset from mononuclear cells and suggest a consensus prediction of proportions when ground truth is not available. We validated the consensus method on data from the stomach and studied its spillover effect. Importantly, we propose the use of the critical assessment of transcriptomic deconvolution (CATD) pipeline which encompasses functionalities for generating references and pseudo-bulks and running implemented deconvolution methods. CATD streamlines simultaneous deconvolution of numerous bulk samples, providing a practical solution for speeding up the evaluation of newly developed methods.

Availability and implementation

https://github.com/Papatheodorou-Group/CATD_snakemake.

CaSSiDI: novel single-cell "Cluster Similarity Scoring and Distinction Index" reveals critical functions for PirB and context-dependent Cebpb repression

PirB is an inhibitory cell surface receptor particularly prominent on myeloid cells. PirB curtails the phenotypes of activated macrophages during inflammation or tumorigenesis, but its functions in macrophage homeostasis are obscure. To elucidate PirB-related functions in macrophages at steady-state, we generated and compared single-cell RNA-sequencing (scRNAseq) datasets obtained from myeloid cell subsets of wild type (WT) and PirB-deficient knockout (PirB KO) mice. To facilitate this analysis, we developed a novel approach to clustering parameter optimization called "Cluster Similarity Scoring and Distinction Index" (CaSSiDI). We demonstrate that CaSSiDI is an adaptable computational framework that facilitates tandem analysis of two scRNAseq datasets by optimizing clustering parameters. We further show that CaSSiDI offers more advantages than a standard Seurat analysis because it allows direct comparison of two or more independently clustered datasets, thereby alleviating the need for batch-correction while identifying the most similar and different clusters. Using CaSSiDI, we found that PirB is a novel regulator of Cebpb expression that controls the generation of Ly6Clo patrolling monocytes and the expansion properties of peritoneal macrophages. PirB's effect on Cebpb is tissue-specific since it was not observed in splenic red pulp macrophages (RPMs). However, CaSSiDI revealed a segregation of the WT RPM population into a CD68loIrf8+ "neuronal-primed" subset and an CD68hiFtl1+ "iron-loaded" subset. Our results establish the utility of CaSSiDI for single-cell assay analyses and the determination of optimal clustering parameters. Our application of CaSSiDI in this study has revealed previously unknown roles for PirB in myeloid cell populations. In particular, we have discovered homeostatic functions for PirB that are related to Cebpb expression in distinct macrophage subsets.