Using singscore to predict mutation status in acute myeloid leukemia from transcriptomic signatures

Integration of single-cell and bulk RNA sequencing data using machine learning identifies oxidative stress-related genes LUM and PCOLCE2 as potential biomarkers for heart failure

Oxidative stress (OS) is a pivotal mechanism driving the progression of cardiovascular diseases, particularly heart failure (HF). However, the comprehensive characterisation of OS-related genes in HF remains largely unexplored. In the present study, we analysed single-cell RNA sequencing datasets from the Gene Expression Omnibus and OS gene sets from GeneCards. We identified 167 OS-related genes potentially linked to HF by applying algorithms, such as AUCell, UCell, singscore, ssgsea, and AddModuleScore, combined with correlation analysis. Subsequently, we used feature selection algorithms, including least absolute shrinkage and selection operator, XGBoost, Boruta, random forest, gradient boosting machines, decision trees, and support vector machine recursive feature elimination, to identify lumican (LUM) and procollagen C-endopeptidase enhancer 2 (PCOLCE2) as key biomarker candidates with significant diagnostic potential. Bulk RNA-sequencing confirmed their elevated expression in patients with HF, highlighting their predictive utility. Single-cell analysis further revealed their upregulation primarily in fibroblasts, emphasising their cell-specific role in HF. To validate these findings, we developed a transverse aortic constriction-induced HF mouse model that showed enhanced cardiac OS activity and significant PCOLCE2 upregulation in the HF group. These results provide strong evidence of the involvement of OS-related mechanisms in HF. Herein, we propose a diagnostic strategy that provides novel insights into the molecular mechanisms underlying HF. However, further studies are required to validate its clinical utility and ensure its application in the diagnosis of HF.

Biopsy Proteome Score Performs Well as an Effect Measure in a Gluten Challenge Trial of Celiac Disease

BACKGROUND AND AIMS

Development of novel treatments for celiac disease is dependent on precise tools to monitor changes in gluten-induced mucosal damage. Current histology measures are subjective and difficult to standardize. Biopsy proteome scoring is an objective alternative to histology which is based on robust changes in biological pathways that directly reflect gluten-induced mucosal damage. In this study, we aimed to evaluate biopsy proteome scoring as an effect measure in a clinical trial setting by measuring intestinal remodeling in response to oral gluten challenge.

METHODS

We analyzed biopsies from a 14-day gluten challenge trial of treated celiac disease patients that consumed 3 g (n = 6) or 10 g (n = 7) gluten per day. Sections from individually embedded biopsies collected before and after challenge were processed for proteome scoring (n = 109) and measurement of villus height-to-crypt depth ratio (n = 58). Proteome scores were compared with histology, intraepithelial lymphocyte frequency and plasma interleukin-2.

RESULTS

Change in proteome scores were significant for the group of patients who consumed 10 g gluten, but not for the group who consumed 3 g gluten. Altogether, 8 of 13 patients had changes in delta proteome scores above the cutoff. Proteome scores correlated with villus height-to-crypt depth ratios both at run-in and at day 15. Proteome scores at day 15 correlated with intraepithelial lymphocyte frequency and with plasma interleukin-2 levels measured 4 hours post-gluten intake.

CONCLUSION

Biopsy proteome scoring is a simple and reliable measure of gluten-induced mucosal remodeling in response to 14-day oral gluten challenge.

CLINICALTRIALS

gov, Number: NCT03409796.

Machine learning algorithms integrate bulk and single-cell RNA data to unveil oxidative stress following intracerebral hemorrhage

BACKGROUND

Increased oxidative stress (OS) activity following intracerebral hemorrhage (ICH) had significantly impacting patient prognosis. Identifying optimal genes associated with OS could enhance the understanding of OS after ICH.

METHODS

We employed single-cell RNA sequencing (scRNA-seq) to investigate the heterogeneity of OS across various cellular tiers following ICH, aiming to acquire biological insights into ICH. We utilized AUCell, Ucell, singscore, ssgsea, and AddModuleScore algorithms, along with correlation analysis, to identify hub genes influencing high OS post-ICH. Furthermore, we employed four machine learning algorithms, eXtreme Gradient Boosting, Boruta, Random Forest, and Least Absolute Shrinkage and Selection Operator, to identify the optimal feature genes. To validate the accuracy of our analysis, we conducted validation in ICH animal experiments.

RESULTS

After analyzing the scRNA-seq dataset using various algorithms, we found that OS activity exhibited heterogeneity across different cellular layers following ICH, with particularly heightened activity observed in monocytes. Further integration of bulk data and machine learning algorithms revealed that ANXA2 and COTL1 were closely associated with high OS after ICH. Our animal experiments demonstrated an increase in OS expression post-ICH. Additionally, the protein expression of ANXA2 and COTL1 was significantly elevated and co-localized with microglia. Pearson correlation coefficient analysis revealed a significant correlation between ANXA2 and OS, indicating strong consistency (r = 0.84, p < 0.05). Similar results were observed for COTL1 and OS (r = 0.69, p < 0.05).

CONCLUSIONS

Following ICH, ANXA2 and COTL1 might penetrate the brain via monocytes, localize within microglia, and enhance OS activity. This might help us better understand OS after ICH.

Identifying a gene signature of metastatic potential by linking pre-metastatic state to ultimate metastatic fate

Identifying the key molecular pathways that enable metastasis by analyzing the eventual metastatic tumor is challenging because the state of the founder subclone likely changes following metastatic colonization. To address this challenge, we labeled primary mouse pancreatic ductal adenocarcinoma (PDAC) subclones with DNA barcodes to characterize their pre-metastatic state using ATAC-seq and RNA-seq and determine their relative in vivo metastatic potential prospectively. We identified a gene signature separating metastasis-high and metastasis-low subclones orthogonal to the normal-to-PDAC and classical-to-basal axes. The metastasis-high subclones feature activation of IL-1 pathway genes and high NF-κB and Zeb/Snail family activity and the metastasis-low subclones feature activation of neuroendocrine, motility, and Wnt pathway genes and high CDX2 and HOXA13 activity. In a functional screen, we validated novel mediators of PDAC metastasis in the IL-1 pathway, including the NF-κB targets Fos and Il23a, and beyond the IL-1 pathway including Myo1b and Tmem40. We scored human PDAC tumors for our signature of metastatic potential from mouse and found that metastases have higher scores than primary tumors. Moreover, primary tumors with higher scores are associated with worse prognosis. We also found that our metastatic potential signature is enriched in other human carcinomas, suggesting that it is conserved across epithelial malignancies. This work establishes a strategy for linking cancer cell state to future behavior, reveals novel functional regulators of PDAC metastasis, and establishes a method for scoring human carcinomas based on metastatic potential.

Biopsy Proteome Scoring to Determine Mucosal Remodeling in Celiac Disease

BACKGROUND & AIMS

Histologic evaluation of gut biopsies is a cornerstone for diagnosis and management of celiac disease (CeD). Despite its wide use, the method depends on proper biopsy orientation, and it suffers from interobserver variability. Biopsy proteome measurement reporting on the tissue state can be obtained by mass spectrometry analysis of formalin-fixed paraffin-embedded tissue. Here we aimed to transform biopsy proteome data into numerical scores that give observer-independent measures of mucosal remodeling in CeD.

METHODS

A pipeline using glass-mounted formalin-fixed paraffin-embedded sections for mass spectrometry-based proteome analysis was established. Proteome data were converted to numerical scores using 2 complementary approaches: a rank-based enrichment score and a score based on machine learning using logistic regression. The 2 scoring approaches were compared with each other and with histology analyzing 18 patients with CeD with biopsies collected before and after treatment with a gluten-free diet as well as biopsies from patients with CeD with varying degree of remission (n = 22). Biopsies from individuals without CeD (n = 32) were also analyzed.

RESULTS

The method yielded reliable proteome scoring of both unstained and H&E-stained glass-mounted sections. The scores of the 2 approaches were highly correlated, reflecting that both approaches pick up proteome changes in the same biological pathways. The proteome scores correlated with villus height-to-crypt depth ratio. Thus, the method is able to score biopsies with poor orientation.

CONCLUSIONS

Biopsy proteome scores give reliable observer and orientation-independent measures of mucosal remodeling in CeD. The proteomic method can readily be implemented by nonexpert laboratories in parallel to histology assessment and easily scaled for clinical trial settings.

Metabolic heterogeneity in early-stage lung adenocarcinoma revealed by RNA-seq and scRNA-seq

PURPOSE

Cancer cells maintain cell growth, division, and survival through altered energy metabolism. However, research on metabolic reprogramming in lung adenocarcinoma (LUAD) is limited METHODS: We downloaded TCGA and GEO sequencing data. Consistent clustering with the ConsensusClusterPlus package was employed to detect the scores for four metabolism-related pathways. The LUAD samples in the TCGA dataset were clustered with ConsensusClusterPlus, and the optimal number of clusters was determined according to the cumulative distribution function (CDF). The cell score for each sample in the TCGA dataset was calculated using the MCPcounter estimate function of the MCPcounter package.

RESULTS

We identified two subtypes by scoring the samples based on the 4 metabolism-related pathways and cluster dimensionality reduction. The prognosis of cluster B was obviously poorer than that of cluster A in patients with LUAD. The analysis of single-nucleotide variation (SNV) data showed that the top 15 genes in the four metabolic pathways with the most mutations were TKTL2, PGK2, HK3, EHHADH, GLUD2, PKLR, TKTL1, HADHB, CPT1C, HK1, HK2, PFKL, SLC2A3, PFKFB1, and CPT1A. The IFNγ score of cluster B was significantly higher than that of cluster A. The immune T-cell lytic activity score of cluster B was significantly higher than that of cluster A. We further identified 5 immune cell subsets from single-cell sequencing data. The top 5 marker genes of B cells were IGHM, JCHAIN, IGLC3, IGHA1, and IGKC. The C0 subgroup of monocytes had a higher pentose phosphate pathway (PPP) score than the C6 subgroup.

CONCLUSIONS

Metabolism-related subtypes could be potential biomarkers in the prognosis prediction and treatment of LUAD.

Computational Screening of Anti-Cancer Drugs Identifies a New BRCA Independent Gene Expression Signature to Predict Breast Cancer Sensitivity to Cisplatin

Simple Summary

Using a collection of publicly available drug screening resources, we identified different partners of genes associated with either sensitivity or resistance to 90 anti-cancer therapies. When subsequently applying these signatures to multiple datasets, we found that these predictive models could predict a large range of drug responses in patient samples. In particular, we discovered a new gene signature to identify breast cancer tumors that are likely to respond to cisplatin in the absence of BRCA1 mutations. This work constitutes an important advance to accelerate the application of platinum-based therapies in patient groups that are not routinely treated with these drugs. In the future, this approach may help to guide the choice of drugs based on the molecular profile of the tumors.

Abstract

The development of therapies that target specific disease subtypes has dramatically improved outcomes for patients with breast cancer. However, survival gains have not been uniform across patients, even within a given molecular subtype. Large collections of publicly available drug screening data matched with transcriptomic measurements have facilitated the development of computational models that predict response to therapy. Here, we generated a series of predictive gene signatures to estimate the sensitivity of breast cancer samples to 90 drugs, comprising FDA-approved drugs or compounds in early development. To achieve this, we used a cell line-based drug screen with matched transcriptomic data to derive in silico models that we validated in large independent datasets obtained from cell lines and patient-derived xenograft (PDX) models. Robust computational signatures were obtained for 28 drugs and used to predict drug efficacy in a set of PDX models. We found that our signature for cisplatin can be used to identify tumors that are likely to respond to this drug, even in absence of the BRCA-1 mutation routinely used to select patients for platinum-based therapies. This clinically relevant observation was confirmed in multiple PDXs. Our study foreshadows an effective delivery approach for precision medicine.

Gene expression signatures as candidate biomarkers of response to PD-1 blockade in non-small cell lung cancers

Although anti-PD-1/PD-L1 monotherapy has achieved clinical success in non-small cell lung cancer (NSCLC), definitive predictive biomarkers remain to be elucidated. In this study, we performed whole-transcriptome sequencing of pretreatment tumor tissue samples and pretreatment and on-treatment whole blood samples (WB) samples obtained from a clinically annotated cohort of NSCLC patients (n = 40) treated with nivolumab (anti-PD-1) monotherapy. Using a single-sample gene set enrichment scoring method, we found that the tumors of responders with lung adenocarcinoma (LUAD, n = 20) are inherently immunogenic to promote antitumor immunity, whereas those with lung squamous cell carcinoma (LUSC, n = 18) have a less immunosuppressive tumor microenvironment. These findings suggested that nivolumab may function as a molecular targeted agent in LUAD and as an immunomodulating agent in LUSC. In addition, our study explains why the reliability of PD-L1 expression on tumor cells as a predictive biomarker for the response to nivolumab monotherapy is quite different between LUAD and LUSC.

Oxford Nanopore MinION Direct RNA-Seq for Systems Biology

Long-read direct RNA sequencing developed by Oxford Nanopore Technologies (ONT) is quickly gaining popularity for transcriptome studies, while fast turnaround time and low cost make it an attractive instrument for clinical applications. There is a growing interest to utilize transcriptome data to unravel activated biological processes responsible for disease progression and response to therapies. This trend is of particular interest for precision medicine which aims at single-patient analysis. Here we evaluated whether gene abundances measured by MinION direct RNA sequencing are suited to produce robust estimates of pathway activation for single sample scoring methods. We performed multiple RNA-seq analyses for a single sample that originated from the HepG2 cell line, namely five ONT replicates, and three replicates using Illumina NovaSeq. Two pathway scoring methods were employed-ssGSEA and singscore. We estimated the ONT performance in terms of detected protein-coding genes and average pairwise correlation between pathway activation scores using an exhaustive computational scheme for all combinations of replicates. In brief, we found that at least two ONT replicates are required to obtain reproducible pathway scores for both algorithms. We hope that our findings may be of interest to researchers planning their ONT direct RNA-seq experiments.

SLUG-related partial epithelial-to-mesenchymal transition is a transcriptomic prognosticator of head and neck cancer survival

Partial epithelial‐to‐mesenchymal transition (pEMT) contributes to cellular heterogeneity that is associated with nodal metastases and unfavorable clinical parameters in head and neck squamous cell carcinomas (HNSCCs). We developed a single‐cell RNA sequencing signature‐based pEMT quantification through cell type‐dependent deconvolution of bulk RNA sequencing and microarray data combined with single‐sample scoring of molecular phenotypes (Singscoring). Clinical pEMT‐Singscores served as molecular classifiers in multivariable Cox proportional hazard models and high scores prognosticated poor overall survival and reduced response to irradiation as independent parameters in large HNSCC cohorts [The Cancer Genome Atlas (TCGA), MD Anderson Cancer Centre (MDACC), Fred Hutchinson Cancer Research Center (FHCRC)]. Differentially expressed genes confirmed enhanced cell motility and reduced oxidative phosphorylation and epithelial differentiation in pEMT^high patients. In patients and cell lines, the EMT transcription factor SLUG correlated most strongly with pEMT‐Singscores and promoted pEMT, enhanced invasion, and resistance to irradiation in vitro. SLUG protein levels in HNSCC predicted disease‐free survival, and its peripheral expression at the interphase to the tumor microenvironment was significantly increased in relapsing patients. Hence, pEMT‐Singscores represent a novel risk predictor for HNSCC stratification regarding clinical outcome and therapy response that is partly controlled by SLUG.

Stable gene expression for normalisation and single-sample scoring

Gene expression signatures have been critical in defining the molecular phenotypes of cells, tissues, and patient samples. Their most notable and widespread clinical application is stratification of breast cancer patients into molecular (PAM50) subtypes. The cost and relatively large amounts of fresh starting material required for whole-transcriptome sequencing has limited clinical application of thousands of existing gene signatures captured in repositories such as the Molecular Signature Database. We identified genes with stable expression across a range of abundances, and with a preserved relative ordering across thousands of samples, allowing signature scoring and supporting general data normalisation for transcriptomic data. Our new method, stingscore, quantifies and summarises relative expression levels of signature genes from individual samples through the inclusion of these ‘stably-expressed genes’. We show that our list of stable genes has better stability across cancer and normal tissue data than previously proposed gene sets. Additionally, we show that signature scores computed from targeted transcript measurements using stingscore can predict docetaxel response in breast cancer patients. This new approach to gene expression signature analysis will facilitate the development of panel-type tests for gene expression signatures, thus supporting clinical translation of the powerful insights gained from cancer transcriptomic studies.

Differential co-expression-based detection of conditional relationships in transcriptional data: comparative analysis and application to breast cancer

BACKGROUND

Elucidation of regulatory networks, including identification of regulatory mechanisms specific to a given biological context, is a key aim in systems biology. This has motivated the move from co-expression to differential co-expression analysis and numerous methods have been developed subsequently to address this task; however, evaluation of methods and interpretation of the resulting networks has been hindered by the lack of known context-specific regulatory interactions.

RESULTS

In this study, we develop a simulator based on dynamical systems modelling capable of simulating differential co-expression patterns. With the simulator and an evaluation framework, we benchmark and characterise the performance of inference methods. Defining three different levels of "true" networks for each simulation, we show that accurate inference of causation is difficult for all methods, compared to inference of associations. We show that a z-score-based method has the best general performance. Further, analysis of simulation parameters reveals five network and simulation properties that explained the performance of methods. The evaluation framework and inference methods used in this study are available in the dcanr R/Bioconductor package.

CONCLUSIONS

Our analysis of networks inferred from simulated data show that hub nodes are more likely to be differentially regulated targets than transcription factors. Based on this observation, we propose an interpretation of the inferred differential network that can reconstruct a putative causal network.