Quantitative assessment of protein activity in orphan tissues and single cells using the metaVIPER algorithm

Epigenetic targeting of PGBD5-dependent DNA damage in SMARCB1-deficient sarcomas

Despite the potential of targeted epigenetic therapies, most cancers do not respond to current epigenetic drugs. The Polycomb repressive complex EZH2 inhibitor tazemetostat was recently approved for the treatment of SMARCB1-deficient epithelioid sarcomas, based on the functional antagonism between PRC2 and loss of SMARCB1. Through the analysis of tazemetostat-treated patient tumors, we recently defined key principles of their response and resistance to EZH2 epigenetic therapy. Here, using transcriptomic inference from SMARCB1-deficient tumor cells, we nominate the DNA damage repair kinase ATR as a target for rational combination EZH2 epigenetic therapy. We show that EZH2 inhibition promotes DNA damage in epithelioid and rhabdoid tumor cells, at least in part via its induction of the transposase-derived PGBD5. We leverage this collateral synthetic lethal dependency to target PGBD5-dependent DNA damage by inhibition of ATR but not CHK1 using elimusertib. Consequently, combined EZH2 and ATR inhibition improves therapeutic responses in diverse patient-derived epithelioid and rhabdoid tumors in vivo. This advances a combination epigenetic therapy based on EZH2-PGBD5 synthetic lethal dependency suitable for immediate translation to clinical trials for patients.

Identification and Pharmacological Targeting of Treatment-Resistant, Stem-like Breast Cancer Cells for Combination Therapy

Tumors frequently harbor isogenic yet epigenetically distinct subpopulations of multi-potent cells with high tumor-initiating potential-often called Cancer Stem-Like Cells (CSLCs). These can display preferential resistance to standard-of-care chemotherapy. Single-cell analyses can help elucidate Master Regulator (MR) proteins responsible for governing the transcriptional state of these cells, thus revealing complementary dependencies that may be leveraged via combination therapy. Interrogation of single-cell RNA sequencing profiles from seven metastatic breast cancer patients, using perturbational profiles of clinically relevant drugs, identified drugs predicted to invert the activity of MR proteins governing the transcriptional state of chemoresistant CSLCs, which were then validated by CROP-seq assays. The top drug, the anthelmintic albendazole, depleted this subpopulation in vivo without noticeable cytotoxicity. Moreover, sequential cycles of albendazole and paclitaxel-a commonly used chemotherapeutic -displayed significant synergy in a patient-derived xenograft (PDX) from a TNBC patient, suggesting that network-based approaches can help develop mechanism-based combinatorial therapies targeting complementary subpopulations.

Statement of significance

Network-based approaches, as shown in a study on metastatic breast cancer, can develop effective combinatorial therapies targeting complementary subpopulations. By analyzing scRNA-seq data and using clinically relevant drugs, researchers identified and depleted chemoresistant Cancer Stem-Like Cells, enhancing the efficacy of standard chemotherapies.

Antileukemia efficacy of the dual BCL2/BCL-XL inhibitor AZD0466 in acute lymphoblastic leukemia preclinical models

ABSTRACT

The upregulation of B-cell lymphoma 2 (BCL2) and B-cell lymphoma-extra large (BCL-XL), 2 proteins in the BCL2 family of proteins, leads to a disproportional expression of prodeath and prosurvival proteins in favor of leukemia survival, tumorigenesis, and chemoresistance. In different subsets of acute lymphoblastic leukemia (ALL), the proportion of these 2 proteins varies, and their potential as therapeutic targets needs detailed characterization. Here, we investigated BCL2 and BCL-XL, the genes that encode BCL2 and BCL-XL, and their expression differences between B-cell acute lymphoblastic leukemia (B-ALL) and T-cell ALL (T-ALL). We also evaluated the therapeutic potential of targeting these proteins with AZD0466, a novel drug-dendrimer conjugate of the BCL2/BCL-XL inhibitor AZD4320, and with BCL2 inhibitor venetoclax (ABT-199). Gene expression and activity analyses supported by the protein expression patterns in ALL cell lines and primary samples demonstrated increased levels of BCL2 expression in B-ALL, with high sensitivity to venetoclax or AZD4320. In contrast, strong BCL-XL expression and sensitivity to dual BCL2/BCL-XL inhibition was observed specifically in T-ALL samples. This observation was confirmed by BH3 profiling, demonstrating BCL2/BCL-XL codependence in T-ALL and BCL2 dependence in B-ALL. In a mouse model of T-ALL, AZD0466 but not venetoclax reduced leukemic burden and prolonged survival without significant toxicities. Our findings therefore suggest that the novel dual BCL2/BCL-XL inhibitor AZD0466 outperforms single BCL2 inhibition by venetoclax in T-ALL. These findings facilitate the translation of dual BCL2/BCL-XL inhibitors into ALL clinical trials, either alone or in combination with standard-of-care chemotherapy and immune therapies.

Exon inclusion signatures enable accurate estimation of splicing factor activity

Splicing factors control exon inclusion in messenger RNA, shaping transcriptome and proteome diversity. Their catalytic activity is regulated by multiple layers, making single-omic measurements on their own fall short in identifying which splicing factors underlie a phenotype. Here, we propose splicing factor activity can be estimated by interpreting changes in exon inclusion. We benchmark methods to construct splicing factor→exon networks and calculate activity. Combining RNA-seq perturbation-based networks with VIPER (virtual inference of protein activity by enriched regulon analysis) accurately captures splicing factor activation modulated by different regulatory layers. This approach consolidates splicing factor regulation into a single score derived solely from exon inclusion signatures, allowing functional interpretation of heterogeneous conditions. As a proof of concept, we identify recurrent cancer splicing programs, revealing oncogenic- and tumor suppressor-like splicing factors missed by conventional methods. These programs correlate with patient survival and key cancer hallmarks: initiation, proliferation, and immune evasion. Altogether, we show splicing factor activity can be accurately estimated from exon inclusion changes, enabling comprehensive analyses of splicing regulation with minimal data requirements.

A computational tool to infer enzyme activity using post-translational modification profiling data

Enzymes play a pivotal role in orchestrating complex cellular responses to external stimuli and environmental changes through signal transduction pathways. Despite their crucial roles, measuring enzyme activities is typically indirect and performed on a smaller scale, unlike protein abundance measured by high-throughput proteomics. Moreover, it is challenging to derive the activity of enzymes from proteome-wide post-translational modification (PTM) profiling data. To address this challenge, we introduce enzyme activity inference with structural equation modeling under the JUMP umbrella (JUMPsem), a novel computational tool designed to infer enzyme activity using PTM profiling data. We demonstrate that the JUMPsem program enables estimating kinase activities using phosphoproteome data, ubiquitin E3 ligase activities from the ubiquitinome, and histone acetyltransferase (HAT) activities based on the acetylome. In addition, JUMPsem is capable of establishing novel enzyme-substrate relationships through searching motif sequences. JUMPsem outperforms widely used kinase activity tools, such as IKAP and KSEA, in terms of the number of kinases and the computational speed. The JUMPsem program is scalable and publicly available as an open-source R package and user-friendly web-based R/Shiny app. Collectively, JUMPsem provides an improved tool for inferring protein enzyme activities, potentially facilitating targeted drug development.

Causal genes identification of giant cell arteritis in CD4+ Memory t cells: an integration of multi-omics and expression quantitative trait locus analysis

BACKGROUND

Giant cell arteritis (GCA) is a prevalent artery and is strongly correlated with age. The role of CD4+ Memory T cells in giant cell arteritis has not been elucidated.

METHOD

Through single-cell analysis, we focused on the CD4+ Memory T cells in giant cell arteritis. eQTL analysis and mendelian randomization analysis identified the significant genes which have a causal effect on giant cell arteritis risk. CD4+ Memory T cells were subsequently divided into gene-positive and gene-negative groups, then further single-cell analysis was conducted. Mendelian randomization of plasma proteins, blood-urine biomarkers and metabolites were also performed. Eventually, the PMA induced Jurkat cell lines were used for biological experiments to explore the specific functions of significant causal genes in CD4+ Memory T cells.

RESULTS

Similarity of CD4+ Memory T cells in GCA and old samples were explored. DDIT4 and ARHGAP15 were identified as significant risk genes via mendelian randomization. The CD4+ Memory T cells were then divided into DDIT4 ± or ARHGAP15 ± groups, and further single-cell analysis indicated the differences in aspects involving intercellular communication, functional pathways, protein activity, metabolism and drug sensitivity between positive and negative groups. In vitro experiments, including overexpression and knockdown, demonstrated that DDIT4 leading to a chronic, low-intensity inflammatory state in CD4+ Memory T cells, eventually promoting the development of GCA.

CONCLUSION

DDIT4 and ARHGAP15 have significant causal effects on giant cell arteritis risk. Specifically, DDIT4 exhibit pro-inflammatory effects on GCA via promotes chronic, low-intensity inflammatory in CD4+ Memory T cell.

Identification and Targeting of Regulators of SARS-CoV-2-Host Interactions in the Airway Epithelium

Although the impact of SARS-CoV-2 in the lung has been extensively studied, the molecular regulators and targets of the host-cell programs hijacked by the virus in distinct human airway epithelial cell populations remain poorly understood. This is in part ascribed to the use of nonprimary cell systems, overreliance on single-cell gene expression profiling that does not ultimately reflect protein activity, and bias toward the downstream effects rather than their mechanistic determinants. Here we address these issues by network-based analysis of single cell transcriptomic profiles of pathophysiologically relevant human adult basal, ciliated and secretory cells to identify master regulator (MR) protein modules controlling their SARS-CoV-2-mediated reprogramming. This uncovered chromatin remodeling, endosomal sorting, ubiquitin pathways, as well as proviral factors identified by CRISPR analyses as components of the host response collectively or selectively activated in these cells. Large-scale perturbation assays, using a clinically relevant drug library, identified 11 drugs able to invert the entire MR signature activated by SARS-CoV-2 in these cell types. Leveraging MR analysis and perturbational profiles of human primary cells represents a novel mechanism-based approach and resource that can be directly generalized to interrogate signatures of other airway conditions for drug prioritization.

Targeting the immune privilege of tumor-initiating cells to enhance cancer immunotherapy

Tumor-initiating cells (TICs) possess the ability to evade anti-tumor immunity, potentially explaining many failures of cancer immunotherapy. Here, we identify CD49f as a prominent marker for discerning TICs in hepatocellular carcinoma (HCC), outperforming other commonly used TIC markers. CD49f-high TICs specifically recruit tumor-promoting neutrophils via the CXCL2-CXCR2 axis and create an immunosuppressive milieu in the tumor microenvironment (TME). Reciprocally, the neutrophils reprogram nearby tumor cells toward a TIC phenotype via secreting CCL4. These cells can evade CD8+ T cell-mediated killing through CCL4/STAT3-induced and CD49f-stabilized CD155 expression. Notably, while aberrant CD155 expression contributes to immune suppression, it also represents a TIC-specific vulnerability. We demonstrate that either CD155 deletion or antibody blockade significantly enhances sensitivity to anti-PD-1 therapy in preclinical HCC models. Our findings reveal a new mechanism of tumor immune evasion and provide a rationale for combining CD155 blockade with anti-PD-1/PD-L1 therapy in HCC.

NOVEL ACTIVE PROTEINS FOR SEPSIS PROGNOSIS REVEALED THROUGH ScRNA-seq AND QUANTITATIVE PROTEOMICS

ABSTRACT

Objective: To uncover critical active proteins influencing sepsis outcomes through multiomics analysis. Methods: This study collected peripheral blood from sepsis patients (NS = 26, SV = 27) and controls (Con = 16). Cellular heterogeneity was assessed using scRNA-seq. Cellular populations were identified through clustering and annotation. Gene set variation analysis was employed to detect pathway alterations in sepsis, while the Viper algorithm estimated protein activity at the single-cell level. Signaling networks were investigated via cell-cell communication analysis. Differentially expressed proteins were identified by DIA proteomics and confirmed through integrated analysis. Prognostic value was evaluated via meta and survival analyses. Results: scRNA-seq of 22,673 features within 34,228 cells identified five cellular clusters and 253 active proteins via Viper, validated by DIA (FC > 2, P < 0.05). Four proteins (SPI1, MEF2A, CBX3, UBTF) with prognostic significance were discovered and mapped onto the cellular landscape. Gene set variation analysis enrichment analysis revealed that the NS group exhibited significant alterations in pathways related to cellular apoptosis and inflammatory responses, while the SV group displayed increased activity in DNA repair and cellular survival pathways. Conclusion: The study's findings advance the understanding of sepsis pathophysiology by linking differentially active proteins to patient prognosis, paving the way for targeted therapeutic strategies.

A forward genetic screen identifies Sirtuin1 as a driver of neuroendocrine prostate cancer

Although localized prostate cancer is relatively indolent, advanced prostate cancer manifests with aggressive and often lethal variants, including neuroendocrine prostate cancer (NEPC). To identify drivers of aggressive prostate cancer, we leveraged Sleeping Beauty (SB) transposon mutagenesis in a mouse model based on prostate-specific loss-of-function of Pten and Tp53 . Compared with control mice, SB mice developed more aggressive prostate tumors, with increased incidence of metastasis. Notably, a significant percentage of the SB prostate tumors display NEPC phenotypes, and the transcriptomic features of these SB mouse tumors recapitulated those of human NEPC. We identified common SB transposon insertion sites (CIS) and prioritized associated CIS-genes differentially expressed in NEPC versus non-NEPC SB tumors. Integrated analysis of CIS-genes encoding for proteins representing upstream, post-translational modulators of master regulators controlling the transcriptional state of SB -mouse and human NEPC tumors identified sirtuin 1 ( Sirt1 ) as a candidate mechanistic determinant of NEPC. Gain-of-function studies in human prostate cancer cell lines confirmed that SIRT1 promotes NEPC, while its loss-of-function or pharmacological inhibition abrogates NEPC. This integrative analysis is generalizable and can be used to identify novel cancer drivers for other malignancies.

Summary

Using an unbiased forward mutagenesis screen in an autochthonous mouse model, we have investigated mechanistic determinants of aggressive prostate cancer. SIRT1 emerged as a key regulator of neuroendocrine prostate cancer differentiation and a potential target for therapeutic intervention.

Etiology-independent activation of the LTβ-LTβR-RELB axis drives aggressiveness and predicts poor prognosis in HCC

BACKGROUND AND AIMS

HCC is the most common primary liver tumor, with an increasing incidence worldwide. HCC is a heterogeneous malignancy and usually develops in a chronically injured liver. The NF-κB signaling network consists of a canonical and a noncanonical branch. Activation of canonical NF-κB in HCC is documented. However, a functional and clinically relevant role of noncanonical NF-κB and its downstream effectors is not established.

APPROACH AND RESULTS

Four human HCC cohorts (total n = 1462) and 4 mouse HCC models were assessed for expression and localization of NF-κB signaling components and activating ligands. In vitro , NF-κB signaling, proliferation, and cell death were measured, proving a pro-proliferative role of v-rel avian reticuloendotheliosis viral oncogene homolog B (RELB) activated by means of NF-κB-inducing kinase. In vivo , lymphotoxin beta was identified as the predominant inducer of RELB activation. Importantly, hepatocyte-specific RELB knockout in a murine HCC model led to a lower incidence compared to controls and lower maximal tumor diameters. In silico , RELB activity and RELB-directed transcriptomics were validated on the The Cancer Genome Atlas HCC cohort using inferred protein activity and Gene Set Enrichment Analysis. In RELB-active HCC, pathways mediating proliferation were significantly activated. In contrast to v-rel avian reticuloendotheliosis viral oncogene homolog A, nuclear enrichment of noncanonical RELB expression identified patients with a poor prognosis in an etiology-independent manner. Moreover, RELB activation was associated with malignant features metastasis and recurrence.

CONCLUSIONS

This study demonstrates a prognostically relevant, etiology-independent, and cross-species consistent activation of a lymphotoxin beta/LTβR/RELB axis in hepatocarcinogenesis. These observations may harbor broad implications for HCC, including possible clinical exploitation.

Elucidating Compound Mechanism of Action and Polypharmacology with a Large-scale Perturbational Profile Compendium

The Mechanism of Action (MoA) of a drug is generally represented as a small, non-tissue-specific repertoire of high-affinity binding targets. Yet, drug activity and polypharmacology are increasingly associated with a broad range of off-target and tissue-specific effector proteins. To address this challenge, we have leveraged a microfluidics-based Plate-Seq technology to survey drug perturbational profiles representing >700 FDA-approved and experimental oncology drugs, in cell lines selected as high-fidelity models of 23 aggressive tumor subtypes. Built on this dataset, we implemented an efficient computational framework to define a tissue-specific protein activity landscape of these drugs and reported almost 50 million differential protein activities derived from drug perturbations vs. vehicle controls. These analyses revealed thousands of highly reproducible and novel, drug-mediated modulation of tissue-specific targets, leading to generation of a proteome-wide drug functional network, characterization of MoA-related drug clusters and off-target effects, dramatical expansion of druggable human proteome, and identification and experimental validation of novel, tissue-specific inhibitors of undruggable oncoproteins, most never reported before. The drug perturbation profile resource described here represents the first, large-scale, whole-genome-wide, RNA-Seq based dataset assembled to date, with the proposed framework, which is easily extended to elucidating the MoA of novel small-molecule libraries, facilitates mechanistic exploration of drug functions, supports systematic and quantitative approaches to precision oncology, and serves as a rich data foundation for drug discovery.

Isthmus progenitor cells contribute to homeostatic cellular turnover and support regeneration following intestinal injury

The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.

Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Elucidation and Pharmacologic Targeting of Master Regulator Dependencies in Coexisting Diffuse Midline Glioma Subpopulations

Diffuse Midline Gliomas (DMGs) are universally fatal, primarily pediatric malignancies affecting the midline structures of the central nervous system. Despite decades of clinical trials, treatment remains limited to palliative radiation therapy. A major challenge is the coexistence of molecularly distinct malignant cell states with potentially orthogonal drug sensitivities. To address this challenge, we leveraged established network-based methodologies to elucidate Master Regulator (MR) proteins representing mechanistic, non-oncogene dependencies of seven coexisting subpopulations identified by single-cell analysis-whose enrichment in essential genes was validated by pooled CRISPR/Cas9 screens. Perturbational profiles of 372 clinically relevant drugs helped identify those able to invert the activity of subpopulation-specific MRs for follow-up in vivo validation. While individual drugs predicted to target individual subpopulations-including avapritinib, larotrectinib, and ruxolitinib-produced only modest tumor growth reduction in orthotopic models, systemic co-administration induced significant survival extension, making this approach a valuable contribution to the rational design of combination therapy.

β-cell Jagged1 is sufficient but not necessary for islet Notch activity and insulin secretory defects in obese mice

OBJECTIVE

Notch signaling, re-activated in β cells from obese mice and causal to β cell dysfunction, is determined in part by transmembrane ligand availability in a neighboring cell. We hypothesized that β cell expression of Jagged1 determines the maladaptive Notch response and resultant insulin secretory defects in obese mice.

METHODS

We assessed expression of Notch pathway components in high-fat diet-fed (HFD) or leptin receptor-deficient (db/db) mice, and performed single-cell RNA sequencing (scRNA-Seq) in islets from patients with and without type 2 diabetes (T2D). We generated and performed glucose tolerance testing in inducible, β cell-specific Jagged1 gain-of- and loss-of-function mice. We also tested effects of monoclonal neutralizing antibodies to Jagged1 in glucose-stimulated insulin secretion (GSIS) assays in isolated islets.

RESULTS

Jag1 was the only Notch ligand that tracked with increased Notch activity in HFD-fed and db/db mice, as well as in metabolically-inflexible β cells enriched in patients with T2D. Neutralizing antibodies to block Jagged1 in islets isolated from HFD-fed and db/db mice potentiated GSIS ex vivo. To demonstrate if β cell Jagged1 is sufficient to cause glucose tolerance in vivo, we generated inducible β cell-specific Jag1 transgenic (β-Jag1TG) and loss-of-function (iβ-Jag1KO) mice. While forced Jagged1 impaired glucose intolerance due to reduced GSIS, loss of β cell Jagged1 did not protect against HFD-induced insulin secretory defects.

CONCLUSIONS

Jagged1 is increased in islets from obese mice and in patients with T2D, and neutralizing Jagged1 antibodies lead to improved GSIS, suggesting that inhibition of Jagged1-Notch signaling may have therapeutic benefit. However, genetic loss-of-function experiments suggest that β cells are not a likely source of the Jagged1 signal.

Tumor Explants Elucidate a Cascade of Paracrine SHH, WNT, and VEGF Signals Driving Pancreatic Cancer Angiosuppression

The sparse vascularity of pancreatic ductal adenocarcinoma (PDAC) presents a mystery: What prevents this aggressive malignancy from undergoing neoangiogenesis to counteract hypoxia and better support growth? An incidental finding from prior work on paracrine communication between malignant PDAC cells and fibroblasts revealed that inhibition of the Hedgehog (HH) pathway partially relieved angiosuppression, increasing tumor vascularity through unknown mechanisms. Initial efforts to study this phenotype were hindered by difficulties replicating the complex interactions of multiple cell types in vitro. Here we identify a cascade of paracrine signals between multiple cell types that act sequentially to suppress angiogenesis in PDAC. Malignant epithelial cells promote HH signaling in fibroblasts, leading to inhibition of noncanonical WNT signaling in fibroblasts and epithelial cells, thereby limiting VEGFR2-dependent activation of endothelial hypersprouting. This cascade was elucidated using human and murine PDAC explant models, which effectively retain the complex cellular interactions of native tumor tissues.

SIGNIFICANCE

We present a key mechanism of tumor angiosuppression, a process that sculpts the physiologic, cellular, and metabolic environment of PDAC. We further present a computational and experimental framework for the dissection of complex signaling cascades that propagate among multiple cell types in the tissue environment. This article is featured in Selected Articles from This Issue, p. 201.

ADGRE5-centered Tsurv model in T cells recognizes responders to neoadjuvant cancer immunotherapy

Background

Neoadjuvant immunotherapy with anti-programmed death-1 (neo-antiPD1) has revolutionized perioperative methods for improvement of overall survival (OS), while approaches for major pathologic response patients' (MPR) recognition along with methods for overcoming non-MPR resistance are still in urgent need.

Methods

We utilized and integrated publicly-available immune checkpoint inhibitors regimens (ICIs) single-cell (sc) data as the discovery datasets, and innovatively developed a cell-communication analysis pipeline, along with a VIPER-based-SCENIC process, to thoroughly dissect MPR-responding subsets. Besides, we further employed our own non-small cell lung cancer (NSCLC) ICIs cohort's sc data for validation in-silico. Afterward, we resorted to ICIs-resistant murine models developed by us with multimodal investigation, including bulk-RNA-sequencing, Chip-sequencing and high-dimensional cytometry by time of flight (CYTOF) to consolidate our findings in-vivo. To comprehensively explore mechanisms, we adopted 3D ex-vivo hydrogel models for analysis. Furthermore, we constructed an ADGRE5-centered Tsurv model from our discovery dataset by machine learning (ML) algorithms for a wide range of tumor types (NSCLC, melanoma, urothelial cancer, etc.) and verified it in peripheral blood mononuclear cells (PBMCs) sc datasets.

Results

Through a meta-analysis of multimodal sequential sc sequencing data from pre-ICIs and post-ICIs, we identified an MPR-expanding T cells meta-cluster (MPR-E) in the tumor microenvironment (TME), characterized by a stem-like CD8+ T cluster (survT) with STAT5-ADGRE5 axis enhancement compared to non-MPR or pre-ICIs TME. Through multi-omics analysis of murine TME, we further confirmed the existence of survT with silenced function and immune checkpoints (ICs) in MPR-E. After verification of the STAT5-ADGRE5 axis of survT in independent ICIs cohorts, an ADGRE5-centered Tsurv model was then developed through ML for identification of MPR patients pre-ICIs and post-ICIs, both in TME and PBMCs, which was further verified in pan-cancer immunotherapy cohorts. Mechanistically, we unveiled ICIs stimulated ADGRE5 upregulation in a STAT5-IL32 dependent manner in a 3D ex-vivo system (3D-HYGTIC) developed by us previously, which marked Tsurv with better survival flexibility, enhanced stemness and potential cytotoxicity within TME.

Conclusion

Our research provides insights into mechanisms underlying MPR in neo-antiPD1 and a well-performed model for the identification of non-MPR.

Deciphering driver regulators of cell fate decisions from single-cell transcriptomics data with CEFCON

Single-cell technologies enable the dynamic analyses of cell fate mapping. However, capturing the gene regulatory relationships and identifying the driver factors that control cell fate decisions are still challenging. We present CEFCON, a network-based framework that first uses a graph neural network with attention mechanism to infer a cell-lineage-specific gene regulatory network (GRN) from single-cell RNA-sequencing data, and then models cell fate dynamics through network control theory to identify driver regulators and the associated gene modules, revealing their critical biological processes related to cell states. Extensive benchmarking tests consistently demonstrated the superiority of CEFCON in GRN construction, driver regulator identification, and gene module identification over baseline methods. When applied to the mouse hematopoietic stem cell differentiation data, CEFCON successfully identified driver regulators for three developmental lineages, which offered useful insights into their differentiation from a network control perspective. Overall, CEFCON provides a valuable tool for studying the underlying mechanisms of cell fate decisions from single-cell RNA-seq data.

Anti-PD-1 immunotherapy with androgen deprivation therapy induces robust immune infiltration in metastatic castration-sensitive prostate cancer

When compared to other malignancies, the tumor microenvironment (TME) of primary and castration-resistant prostate cancer (CRPC) is relatively devoid of immune infiltrates. While androgen deprivation therapy (ADT) induces a complex immune infiltrate in localized prostate cancer, the composition of the TME in metastatic castration-sensitive prostate cancer (mCSPC), and the effects of ADT and other treatments in this context are poorly understood. Here, we perform a comprehensive single-cell RNA sequencing (scRNA-seq) profiling of metastatic sites from patients participating in a phase 2 clinical trial (NCT03951831) that evaluated standard-of-care chemo-hormonal therapy combined with anti-PD-1 immunotherapy. We perform a longitudinal, protein activity-based analysis of TME subpopulations, revealing immune subpopulations conserved across multiple metastatic sites. We also observe dynamic changes in these immune subpopulations in response to treatment and a correlation with clinical outcomes. Our study uncovers a therapy-resistant, transcriptionally distinct tumor subpopulation that expands in cell number in treatment-refractory patients.

Master regulator activity QTL protocol to implicate regulatory pathways potentially mediating GWAS signals using eQTL data

Here, we present a protocol to identify transcriptional regulators potentially mediating downstream biological effects of germline variants associated with complex traits of interest, which enables functional hypothesis generation independent of colocalizing expression quantitative trait loci (eQTLs). We describe steps for tissue-/cell-type-specific co-expression network modeling, expression regulator activity inference, and identification of representative phenotypic master regulators. Finally, we detail activity QTL and eQTL analyses. This protocol requires genotype, expression, and relevant covariables and phenotype data from existing eQTL datasets. For complete details on the use and execution of this protocol, please refer to Hoskins et al.1.

Reversing pancreatic β-cell dedifferentiation in the treatment of type 2 diabetes

The maintenance of glucose homeostasis is fundamental for survival and health. Diabetes develops when glucose homeostasis fails. Type 2 diabetes (T2D) is characterized by insulin resistance and pancreatic β-cell failure. The failure of β-cells to compensate for insulin resistance results in hyperglycemia, which in turn drives altered lipid metabolism and β-cell failure. Thus, insulin secretion by pancreatic β-cells is a primary component of glucose homeostasis. Impaired β-cell function and reduced β-cell mass are found in diabetes. Both features stem from a failure to maintain β-cell identity, which causes β-cells to dedifferentiate into nonfunctional endocrine progenitor-like cells or to trans-differentiate into other endocrine cell types. In this regard, one of the key issues in achieving disease modification is how to reestablish β-cell identity. In this review, we focus on the causes and implications of β-cell failure, as well as its potential reversibility as a T2D treatment.

A Transcriptome-Based Precision Oncology Platform for Patient-Therapy Alignment in a Diverse Set of Treatment-Resistant Malignancies

Predicting in vivo response to antineoplastics remains an elusive challenge. We performed a first-of-kind evaluation of two transcriptome-based precision cancer medicine methodologies to predict tumor sensitivity to a comprehensive repertoire of clinically relevant oncology drugs, whose mechanism of action we experimentally assessed in cognate cell lines. We enrolled patients with histologically distinct, poor-prognosis malignancies who had progressed on multiple therapies, and developed low-passage, patient-derived xenograft models that were used to validate 35 patient-specific drug predictions. Both OncoTarget, which identifies high-affinity inhibitors of individual master regulator (MR) proteins, and OncoTreat, which identifies drugs that invert the transcriptional activity of hyperconnected MR modules, produced highly significant 30-day disease control rates (68% and 91%, respectively). Moreover, of 18 OncoTreat-predicted drugs, 15 induced the predicted MR-module activity inversion in vivo. Predicted drugs significantly outperformed antineoplastic drugs selected as unpredicted controls, suggesting these methods may substantively complement existing precision cancer medicine approaches, as also illustrated by a case study.

SIGNIFICANCE

Complementary precision cancer medicine paradigms are needed to broaden the clinical benefit realized through genetic profiling and immunotherapy. In this first-in-class application, we introduce two transcriptome-based tumor-agnostic systems biology tools to predict drug response in vivo. OncoTarget and OncoTreat are scalable for the design of basket and umbrella clinical trials. This article is highlighted in the In This Issue feature, p. 1275.

NaRnEA: An Information Theoretic Framework for Gene Set Analysis

Gene sets are being increasingly leveraged to make high-level biological inferences from transcriptomic data; however, existing gene set analysis methods rely on overly conservative, heuristic approaches for quantifying the statistical significance of gene set enrichment. We created Nonparametric analytical-Rank-based Enrichment Analysis (NaRnEA) to facilitate accurate and robust gene set analysis with an optimal null model derived using the information theoretic Principle of Maximum Entropy. By measuring the differential activity of ~2500 transcriptional regulatory proteins based on the differential expression of each protein's transcriptional targets between primary tumors and normal tissue samples in three cohorts from The Cancer Genome Atlas (TCGA), we demonstrate that NaRnEA critically improves in two widely used gene set analysis methods: Gene Set Enrichment Analysis (GSEA) and analytical-Rank-based Enrichment Analysis (aREA). We show that the NaRnEA-inferred differential protein activity is significantly correlated with differential protein abundance inferred from independent, phenotype-matched mass spectrometry data in the Clinical Proteomic Tumor Analysis Consortium (CPTAC), confirming the statistical and biological accuracy of our approach. Additionally, our analysis crucially demonstrates that the sample-shuffling empirical null models leveraged by GSEA and aREA for gene set analysis are overly conservative, a shortcoming that is avoided by the newly developed Maximum Entropy analytical null model employed by NaRnEA.

Network-based elucidation of colon cancer drug resistance by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. By leveraging progress in proteomic technologies and network-based methodologies, over the past decade, we developed VESPA-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and used it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogation of tumor-specific enzyme/substrate interactions accurately inferred kinase and phosphatase activity, based on their inferred substrate phosphorylation state, effectively accounting for signal cross-talk and sparse phosphoproteome coverage. The analysis elucidated time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring that was experimentally confirmed by CRISPRko assays, suggesting broad applicability to cancer and other diseases.

Integrated regulatory and metabolic networks of the tumor microenvironment for therapeutic target prioritization

Translation of genomic discovery, such as single-cell sequencing data, to clinical decisions remains a longstanding bottleneck in the field. Meanwhile, computational systems biological models, such as cellular metabolism models and cell signaling pathways, have emerged as powerful approaches to provide efficient predictions in metabolites and gene expression levels, respectively. However, there has been limited research on the integration between these two models. This work develops a methodology for integrating computational models of probabilistic gene regulatory networks with a constraint-based metabolism model. By using probabilistic reasoning with Bayesian Networks, we aim to predict cell-specific changes under different interventions, which are embedded into the constraint-based models of metabolism. Applications to single-cell sequencing data of glioblastoma brain tumors generate predictions about the effects of pharmaceutical interventions on the regulatory network and downstream metabolisms in different cell types from the tumor microenvironment. The model presents possible insights into treatments that could potentially suppress anaerobic metabolism in malignant cells with minimal impact on other cell types' metabolism. The proposed integrated model can guide therapeutic target prioritization, the formulation of combination therapies, and future drug discovery. This model integration framework is also generalizable to other applications, such as different cell types, organisms, and diseases.

Single-cell protein activity analysis reveals a novel subpopulation of chondrocytes and the corresponding key master regulator proteins associated with anti-senescence and OA progression

Background

Osteoarthritis (OA) is a prevalent senescence-related disease with substantial joint pain, loss of joint function, and cartilage degeneration. Because of the paucity of single-cell studies of OA and the gene dropout problem of single-cell RNA sequencing, it is difficult to acquire an in-depth understanding of the molecular characteristics of various chondrocyte clusters.

Methods

Here, we aimed to provide new insights into chondrocyte senescence and a rationale for the development of effective intervention strategies for OA by using published single-cell RNA-sequencing data sets and the metaVIPER algorithm (Virtual Inference of Protein activity by Enriched Regulon). This algorithm was employed to present a proteome catalog of 62,449 chondrocytes from the cartilage of healthy individuals and OA patients at single-cell resolution. Furthermore, histopathologic analysis was carried out in cartilage samples from clinical patients and experimental mouse models of OA to validate above results.

Results

We identified 16 protein-activity-based chondrocyte clusters as well as the underlying master regulators in each cluster. By assessing the enrichment score of each cluster in bulk RNA-sequencing data, followed by gene-set variation analysis, we preliminarily identified a novel subpopulation of chondrocytes (cluster 3). This clinically relevant cluster was predicted to be the main chondrocyte cluster responsible for maintaining cellular homeostasis and anti-senescence. Specifically, we uncovered a set of the key leading-edge proteins of cluster 3 by validating the robustness of the above results using another human chondrocyte single-cell RNA-sequencing data set, consisting of 24,675 chondrocytes. Furthermore, cartilage samples from clinical patients and experimental mouse models of OA were used to evaluate the expression patterns of these leading-edge proteins, and the results indicated that NDRG2, TSPYL2, JMJD6 and HMGB2 are closely associated with OA pathogenesis and might play critical roles in modulating cellular homeostasis and anti-senescence in chondrocytes.

Conclusion

Our study revealed a novel subpopulation of chondrocytes that are critical for anti-progression of OA and the corresponding master regulator proteins, which might serve as therapeutic targets in OA.

What Makes a Functional Gene Regulatory Network? A Circuit Motif Analysis

One of the key questions in systems biology is to understand the roles of gene regulatory circuits in determining cellular states and their functions. In previous studies, some researchers have inferred large gene networks from genome wide genomics/transcriptomics data using the top-down approach, while others have modeled core gene circuits of small sizes using the bottom-up approach. Despite many existing systems biology studies, there is still no general rule on what sizes of gene networks and what types of circuit motifs a system would need to achieve robust biological functions. Here, we adopt a gene circuit motif analysis to discover four-node circuits responsible for multiplicity (rich in dynamical behavior), flexibility (versatile to alter gene expression), or both. We identify the most reoccurring two-node circuit motifs and the co-occurring motif pairs. Furthermore, we investigate the contributing factors of multiplicity and flexibility for large gene networks of different types and sizes. We find that gene networks of intermediate sizes tend to have combined high levels of multiplicity and flexibility. Our study will contribute to a better understanding of the dynamical mechanisms of gene regulatory circuits and provide insights into rational designs of robust gene circuits in synthetic and systems biology.

Elevating microglia TREM2 reduces amyloid seeding and suppresses disease-associated microglia

TREM2 is exclusively expressed by microglia in the brain and is strongly linked to the risk for Alzheimer's disease (AD). As microglial responses modulated by TREM2 are central to AD pathogenesis, enhancing TREM2 signaling has been explored as an AD therapeutic strategy. However, the effective therapeutic window targeting TREM2 is unclear. Here, by using microglia-specific inducible mouse models overexpressing human wild-type TREM2 (TREM2-WT) or R47H risk variant (TREM2-R47H), we show that TREM2-WT expression reduces amyloid deposition and neuritic dystrophy only during the early amyloid seeding stage, whereas TREM2-R47H exacerbates amyloid burden during the middle amyloid rapid growth stage. Single-cell RNA sequencing reveals suppressed disease-associated microglia (DAM) signature and reduced DAM population upon TREM2-WT expression in the early stage, whereas upregulated antigen presentation pathway is detected with TREM2-R47H expression in the middle stage. Together, our findings highlight the dynamic effects of TREM2 in modulating AD pathogenesis and emphasize the beneficial effect of enhancing TREM2 function in the early stage of AD development.

Validation of a non-oncogene encoded vulnerability to exportin 1 inhibition in pediatric renal tumors

BACKGROUND

Malignant rhabdoid tumors (MRTs) and Wilms' tumors (WTs) are rare and aggressive renal tumors of infants and young children comprising ∼5% of all pediatric cancers. MRTs are among the most genomically stable cancers, and although WTs are genomically heterogeneous, both generally lack therapeutically targetable genetic mutations.

METHODS

Comparative protein activity analysis of MRTs (n = 68) and WTs (n = 132) across TCGA and TARGET cohorts, using metaVIPER, revealed elevated exportin 1 (XPO1) inferred activity. In vitro studies were performed on a panel of MRT and WT cell lines to evaluate effects on proliferation and cell-cycle progression following treatment with the selective XPO1 inhibitor selinexor. In vivo anti-tumor activity was assessed in patient-derived xenograft (PDX) models of MRTs and WTs.

FINDINGS

metaVIPER analysis identified markedly aberrant activation of XPO1 in MRTs and WTs compared with other tumor types. All MRT and most WT cell lines demonstrated baseline, aberrant XPO1 activity with in vitro sensitivity to selinexor via cell-cycle arrest and induction of apoptosis. In vivo, XPO1 inhibitors significantly abrogated tumor growth in PDX models, inducing effective disease control with sustained treatment. Corroborating human relevance, we present a case report of a child with multiply relapsed WTs with prolonged disease control on selinexor.

CONCLUSIONS

We report on a novel systems-biology-based comparative framework to identify non-genetically encoded vulnerabilities in genomically quiescent pediatric cancers. These results have provided preclinical rationale for investigation of XPO1 inhibitors in an upcoming investigator-initiated clinical trial of selinexor in children with MRTs and WTs and offer opportunities for exploration of inferred XPO1 activity as a potential predictive biomarker for response.

FUNDING

This work was funded by CureSearch for Children's Cancer, Alan B. Slifka Foundation, NIH (U01 CA217858, S10 OD012351, and S10 OD021764), Michael's Miracle Cure, Hyundai Hope on Wheels, Cannonball Kids Cancer, Conquer Cancer the ASCO Foundation, Cycle for Survival, Paulie Strong Foundation, and the Grayson Fund.

Cell Layers: uncovering clustering structure in unsupervised single-cell transcriptomic analysis

Motivation

Unsupervised clustering of single-cell transcriptomics is a powerful method for identifying cell populations. Static visualization techniques for single-cell clustering only display results for a single resolution parameter. Analysts will often evaluate more than one resolution parameter but then only report one.

Results

We developed Cell Layers, an interactive Sankey tool for the quantitative investigation of gene expression, co-expression, biological processes and cluster integrity across clustering resolutions. Cell Layers enhances the interpretability of single-cell clustering by linking molecular data and cluster evaluation metrics, providing novel insight into cell populations.

Availability and implementation

https://github.com/apblair/CellLayers.

A model for network-based identification and pharmacological targeting of aberrant, replication-permissive transcriptional programs induced by viral infection

SARS-CoV-2 hijacks the host cell transcriptional machinery to induce a phenotypic state amenable to its replication. Here we show that analysis of Master Regulator proteins representing mechanistic determinants of the gene expression signature induced by SARS-CoV-2 in infected cells revealed coordinated inactivation of Master Regulators enriched in physical interactions with SARS-CoV-2 proteins, suggesting their mechanistic role in maintaining a host cell state refractory to virus replication. To test their functional relevance, we measured SARS-CoV-2 replication in epithelial cells treated with drugs predicted to activate the entire repertoire of repressed Master Regulators, based on their experimentally elucidated, context-specific mechanism of action. Overall, 15 of the 18 drugs predicted to be effective by this methodology induced significant reduction of SARS-CoV-2 replication, without affecting cell viability. This model for host-directed pharmacological therapy is fully generalizable and can be deployed to identify drugs targeting host cell-based Master Regulator signatures induced by virtually any pathogen.

Clinical assessment and molecular mechanism of the upregulation of Toll-like receptor 2 (TLR2) in myocardial infarction

Objective

The prevalence and mortality of cardiovascular diseases remain ranked first worldwide. Myocardial infarction (MI) is the central cause of death from cardiovascular diseases, seriously endangering human health. The clinical implication of toll-like receptor 2 (TLR2) remains contradictory, and its mechanism is still unknown. Hence, the objective of this study was to elucidate the clinical value and molecular mechanism of TLR2 in MI.

Methods

All high-throughput datasets and eligible literature were screened, and the expression levels of TLR2 were collected from the MI. The integrated expression level of TLR2 was displayed by calculating the standardized mean difference (SMD) and the area under the curve (AUC) of the summary receiver operating characteristic curve (sROC). The related TLR2 genes were sent for pathway analyses by gene ontology (GO), Kyoto encyclopedia of genes and genome (KEGG), and disease ontology (DO). Single-cell RNA-seq was applied to ascertain the molecular mechanism of TLR2 in MI.

Results

Nine microarrays and four reported data were available to calculate the comprehensive expression level of TLR2 in MI, including 325 cases of MI and 306 cases of controls. The SMD was 2.55 (95% CI = 1.35–3.75), and the AUC was 0.76 (95% CI = 0.72–0.79), indicating the upregulation of TLR2 in MI. The related TLR2 genes were primarily enriched in the pathways of atherosclerosis, arteriosclerotic cardiovascular disease, and arteriosclerosis, suggesting the clinical role of TLR2 in the progression of MI. Afterward, TLR2 was upregulated in myeloid cells in MI.

Conclusions

TLR2 may have a crucial role in progressing from coronary atherosclerosis to MI. The upregulation of TLR2 may have a favorable screening value for MI.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12872-022-02754-y.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Adult enteric Dclk1-positive glial and neuronal cells reveal distinct responses to acute intestinal injury

Intestinal ganglionic cells in the adult enteric nervous system (ENS) are continually exposed to stimuli from the surrounding microenvironment and need at times to respond to disturbed homeostasis following acute intestinal injury. The kinase DCLK1 and intestinal Dclk1-positive cells have been reported to contribute to intestinal regeneration. Although Dclk1-positive cells are present in adult enteric ganglia, their cellular identity and response to acute injury have not been investigated in detail. Here, we reveal the presence of distinct Dclk1-tdTom+/CD49b+ glial-like and Dclk1-tdTom+/CD49b- neuronal cell types in adult myenteric ganglia. These ganglionic cells demonstrate distinct patterns of tracing over time yet show a similar expansion in response to elevated serotonergic signaling. Interestingly, Dclk1-tdTom+ glial-like and neuronal cell types appear resistant to acute irradiation injury-mediated cell death. Moreover, Dclk1-tdTom+/CD49b+ glial-like cells show prominent changes in gene expression profiles induced by injury, in contrast to Dclk1-tdTom+/CD49b- neuronal cell types. Finally, subsets of Dclk1-tdTom+/CD49b+ glial-like cells demonstrate prominent overlap with Nestin and p75NTR and strong responses to elevated serotonergic signaling or acute injury. These findings, together with their role in early development and their neural crest-like gene expression signature, suggest the presence of reserve progenitor cells in the adult Dclk1 glial cell lineage.NEW & NOTEWORTHY The kinase DCLK1 identifies glial-like and neuronal cell types in adult murine enteric ganglia, which resist acute injury-mediated cell death yet differ in their cellular response to injury. Interestingly, Dclk1-labeled glial-like cells show prominent transcriptional changes in response to injury and harbor features reminiscent of previously described enteric neural precursor cells. Our data thus add to recently emerging evidence of reserve cellular plasticity in the adult enteric nervous system.

Immunostimulatory Cancer-Associated Fibroblast Subpopulations Can Predict Immunotherapy Response in Head and Neck Cancer

PURPOSE

Cancer-associated fibroblasts (CAF) have been implicated as potential mediators of checkpoint immunotherapy response. However, the extensive heterogeneity of these cells has precluded rigorous understanding of their immunoregulatory role in the tumor microenvironment.

EXPERIMENTAL DESIGN

We performed high-dimensional single-cell RNA sequencing (scRNA-seq) on four patient tumors pretreatment and posttreatment from a neoadjuvant trial of patients with advanced-stage head and neck squamous cell carcinoma that were treated with the αPD-1 therapy, nivolumab. The head and neck CAF (HNCAF) protein activity profiles, derived from this cohort of paired scRNA-seq, were used to perform protein activity enrichment analysis on the 28-patient parental cohort of clinically annotated bulk transcriptomic profiles. Ex vivo coculture assays were used to test functional relevance of HNCAF subtypes.

RESULTS

Fourteen distinct cell types were identified with the fibroblast population showing significant changes in abundance following nivolumab treatment. Among the fibroblast subtypes, HNCAF-0/3 emerged as predictive of nivolumab response, while HNCAF-1 was associated with immunosuppression. Functionally, HNCAF-0/3 were found to reduce TGFβ-dependent PD-1+TIM-3+ exhaustion of CD8 T cells, increase CD103+NKG2A+ resident memory phenotypes, and enhance the overall cytolytic profile of T cells.

CONCLUSIONS

Our findings demonstrate the functional importance of distinct HNCAF subsets in modulating the immunoregulatory milieu of human HNSCC. In addition, we have identified clinically actionable HNCAF subtypes that can be used as a biomarker of response and resistance in future clinical trials.

Single-cell analysis reveals Comma-1D as a unique cell model for mammary gland development and breast cancer

The mammary epithelial tree contains two distinct populations, luminal and basal. The investigation of how this heterogeneity is developed and how it influences tumorigenesis has been hampered by the need to perform these studies using animal models. Comma-1D is an immortalized mouse mammary epithelial cell line that has unique morphogenetic properties. By performing single-cell RNA-seq studies we found that Comma-1D cultures consist of two main populations with luminal and basal features and a smaller population with mixed lineage and bipotent characteristics. We demonstrated that multiple transcription factors associated with the differentiation of the mammary epithelium in vivo also modulate this process in Comma-1D cultures. Additionally, we found that only cells with luminal features were able to acquire transformed characteristics after an oncogenic HER2 mutant was introduced in their genomes. Overall, our studies characterize at a single-cell level the heterogeneity of the Comma-1D cell line and illustrate how Comma-1D cells can be used as an experimental model to study both the differentiation and the transformation processes in vitro.

Notch-mediated Ephrin signaling disrupts islet architecture and β cell function

Altered islet architecture is associated with β cell dysfunction and type 2 diabetes (T2D) progression, but molecular effectors of islet spatial organization remain mostly unknown. Although Notch signaling is known to regulate pancreatic development, we observed “reactivated” β cell Notch activity in obese mouse models. To test the repercussions and reversibility of Notch effects, we generated doxycycline-dependent, β cell–specific Notch gain-of-function mice. As predicted, we found that Notch activation in postnatal β cells impaired glucose-stimulated insulin secretion and glucose intolerance, but we observed a surprising remnant glucose intolerance after doxycycline withdrawal and cessation of Notch activity, associated with a marked disruption of normal islet architecture. Transcriptomic screening of Notch-active islets revealed increased Ephrin signaling. Commensurately, exposure to Ephrin ligands increased β cell repulsion and impaired murine and human pseudoislet formation. Consistent with our mouse data, Notch and Ephrin signaling were increased in metabolically inflexible β cells in patients with T2D. These studies suggest that β cell Notch/Ephrin signaling can permanently alter islet architecture during a morphogenetic window in early life.

CITEMOXMBD: A flexible single-cell multimodal omics analysis framework to reveal the heterogeneity of immune cells

Simultaneous measurement of multiple modalities in single-cell analysis, represented by CITE-seq, is a promising approach to link transcriptional changes to cellular phenotype and function, requiring new computational methods to define cellular subtypes and states based on multiple data types. Here, we design a flexible single-cell multimodal analysis framework, called CITEMO, to integrate the transcriptome and antibody-derived tags (ADT) data to capture cell heterogeneity from the multi omics perspective. CITEMO uses Principal Component Analysis (PCA) to obtain a low-dimensional representation of the transcriptome and ADT, respectively, and then employs PCA again to integrate these low-dimensional multimodal data for downstream analysis. To investigate the effectiveness of the CITEMO framework, we apply CITEMO to analyse the cell subtypes of Cord Blood Mononuclear Cells (CBMC) samples. Results show that the CITEMO framework can comprehensively analyse single-cell multimodal samples and accurately identify cell subtypes. Besides, we find some specific immune cells that co-express multiple ADT markers. To better describe the co-expression phenomenon, we introduce the co-expression entropy to measure the heterogeneous distribution of the ADT combinations. To further validate the robustness of the CITEMO framework, we analyse Human Bone Marrow Cell (HBMC) samples and identify different states of the same cell type. CITEMO has an excellent performance in identifying cell subtypes and states for multimodal omics data. We suggest that the flexible design idea of CITEMO can be an inspiration for other single-cell multimodal tasks. The complete source code and dataset of the CITEMO framework can be obtained from https://github.com/studentiz/CITEMO.

Transforming Growth Factor-β Activated Kinase 1 (Tak1) Is Activated in Hepatocellular Carcinoma, Mediates Tumor Progression, and Predicts Unfavorable Outcome

Simple Summary

Chronic inflammation is known to drive cancer initiation and progression in the liver and other organs. In different genetic mouse models, the role of the pro-inflammatory kinase Tak1 in liver cancer development has been controversial so far. To clarify the role of Tak1 in human hepatocellular carcinoma (HCC), we investigated the expression of Tak1 in a large and clinicopathologically well-characterized patient cohort with HCC. In human livers and HCCs, Tak1 is predominantly present in its isoform Tak1A localizing to the cell nucleus. Tak1 is upregulated in HCCs of the diethylnitrosamine mouse model as well as in human HCCs, independent of etiology, and is further induced in distant metastases. Overexpression of the isoform Tak1A in the HCC cell line Huh7 resulted in increased tumor cell migration, whereas overexpression of full-length Tak1 had no significant effect. In human HCCs, high nuclear Tak1 expression is associated with vascular invasion and short overall survival.

Abstract

Although knowledge on inflammatory signaling pathways driving cancer initiation and progression has been increasing, molecular mechanisms in hepatocarcinogenesis are still far from being completely understood. Hepatocyte-specific deletion of the MAPKKK Tak1 in mice recapitulates important steps of hepatocellular carcinoma (HCC) development, including the occurrence of cell death, steatohepatitis, dysplastic nodules, and HCCs. However, overactivation of Tak1 in mice upon deletion of its deubiquitinase Cyld also results in steatohepatitis and HCC development. To investigate Tak1 and Cyld in human HCCs, we created a tissue microarray to analyze their expression by immunohistochemistry in a large and well-characterized cohort of 871 HCCs of 561 patients. In the human liver and HCC, Tak1 is predominantly present as its isoform Tak1A and predominantly localizes to cell nuclei. Tak1 is upregulated in diethylnitrosamine-induced mouse HCCs as well as in human HCCs independent of etiology and is further induced in distant metastases. A high nuclear Tak1 expression is associated with short survival and vascular invasion. When we overexpressed Tak1A in Huh7 cells, we observed increased tumor cell migration, whereas overexpression of full-length Tak1 had no significant effect. A combined score of low Cyld and high Tak1 expression was an independent prognostic marker in a multivariate Cox regression model.

BACH2 inhibition reverses β cell failure in type 2 diabetes models

Type 2 diabetes (T2D) is associated with defective insulin secretion and reduced β cell mass. Available treatments provide a temporary reprieve, but secondary failure rates are high, making insulin supplementation necessary. Reversibility of β cell failure is a key translational question. Here, we reverse engineered and interrogated pancreatic islet-specific regulatory networks to discover T2D-specific subpopulations characterized by metabolic inflexibility and endocrine progenitor/stem cell features. Single-cell gain- and loss-of-function and glucose-induced Ca2+ flux analyses of top candidate master regulatory (MR) proteins in islet cells validated transcription factor BACH2 and associated epigenetic effectors as key drivers of T2D cell states. BACH2 knockout in T2D islets reversed cellular features of the disease, restoring a nondiabetic phenotype. BACH2-immunoreactive islet cells increased approximately 4-fold in diabetic patients, confirming the algorithmic prediction of clinically relevant subpopulations. Treatment with a BACH inhibitor lowered glycemia and increased plasma insulin levels in diabetic mice, and restored insulin secretion in diabetic mice and human islets. The findings suggest that T2D-specific populations of failing β cells can be reversed and indicate pathways for pharmacological intervention, including via BACH2 inhibition.

Deciphering regulatory protein activity in human pancreatic islets via reverse engineering of single-cell sequencing data

The loss of functional β cell mass contributes to development and progression of type 2 diabetes (T2D). However, the molecular mechanisms differentiating islet dysfunction in T2D from nondiabetic states remain elusive. In this issue of the JCI, Son et al. applied reverse engineering to obtain the activity of gene expression regulatory proteins from single-cell RNA sequencing data of nondiabetic and T2D human islets. The authors identify unique patterns of regulatory protein activities associated with T2D. Furthermore, BACH2 emerged as a potential transcription factor that drives activation of T2D-associated regulatory proteins in human islets.

Inferred expression regulator activities suggest genes mediating cardiometabolic genetic signals

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator's target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals.

Interpretation of cancer mutations using a multiscale map of protein systems

A major goal of cancer research is to understand how mutations distributed across diverse genes affect common cellular systems, including multiprotein complexes and assemblies. Two challenges—how to comprehensively map such systems and how to identify which are under mutational selection—have hindered this understanding. Accordingly, we created a comprehensive map of cancer protein systems integrating both new and published multi-omic interaction data at multiple scales of analysis. We then developed a unified statistical model that pinpoints 395 specific systems under mutational selection across 13 cancer types. This map, called NeST (Nested Systems in Tumors), incorporates canonical processes and notable discoveries, including a PIK3CA-actomyosin complex that inhibits phosphatidylinositol 3-kinase signaling and recurrent mutations in collagen complexes that promote tumor proliferation. These systems can be used as clinical biomarkers and implicate a total of 548 genes in cancer evolution and progression. This work shows how disparate tumor mutations converge on protein assemblies at different scales.

VEGA is an interpretable generative model for inferring biological network activity in single-cell transcriptomics

Deep learning architectures such as variational autoencoders have revolutionized the analysis of transcriptomics data. However, the latent space of these variational autoencoders offers little to no interpretability. To provide further biological insights, we introduce a novel sparse Variational Autoencoder architecture, VEGA (VAE Enhanced by Gene Annotations), whose decoder wiring mirrors user-provided gene modules, providing direct interpretability to the latent variables. We demonstrate the performance of VEGA in diverse biological contexts using pathways, gene regulatory networks and cell type identities as the gene modules that define its latent space. VEGA successfully recapitulates the mechanism of cellular-specific response to treatments, the status of master regulators as well as jointly revealing the cell type and cellular state identity in developing cells. We envision the approach could serve as an explanatory biological model for development and drug treatment experiments.

A molecular single-cell lung atlas of lethal COVID-19

Respiratory failure is the leading cause of death in patients with severe SARS-CoV-2 infection^1,2, but the host response at the lung tissue level is poorly understood. Here we performed single-nucleus RNA sequencing of about 116,000 nuclei from the lungs of nineteen individuals who died of COVID-19 and underwent rapid autopsy and seven control individuals. Integrated analyses identified substantial alterations in cellular composition, transcriptional cell states, and cell-to-cell interactions, thereby providing insight into the biology of lethal COVID-19. The lungs from individuals with COVID-19 were highly inflamed, with dense infiltration of aberrantly activated monocyte-derived macrophages and alveolar macrophages, but had impaired T cell responses. Monocyte/macrophage-derived interleukin-1β and epithelial cell-derived interleukin-6 were unique features of SARS-CoV-2 infection compared to other viral and bacterial causes of pneumonia. Alveolar type 2 cells adopted an inflammation-associated transient progenitor cell state and failed to undergo full transition into alveolar type 1 cells, resulting in impaired lung regeneration. Furthermore, we identified expansion of recently described CTHRC1⁺ pathological fibroblasts³ contributing to rapidly ensuing pulmonary fibrosis in COVID-19. Inference of protein activity and ligand-receptor interactions identified putative drug targets to disrupt deleterious circuits. This atlas enables the dissection of lethal COVID-19, may inform our understanding of long-term complications of COVID-19 survivors, and provides an important resource for therapeutic development.

Elucidating the cellular dynamics of the brain with single-cell RNA sequencing

Single-cell RNA-sequencing (scRNA-seq) has emerged in recent years as a breakthrough technology to understand RNA metabolism at cellular resolution. In addition to allowing new cell types and states to be identified, scRNA-seq can permit cell-type specific differential gene expression changes, pre-mRNA processing events, gene regulatory networks and single-cell developmental trajectories to be uncovered. More recently, a new wave of multi-omic adaptations and complementary spatial transcriptomics workflows have been developed that facilitate the collection of even more holistic information from individual cells. These developments have unprecedented potential to provide penetrating new insights into the basic neural cell dynamics and molecular mechanisms relevant to the nervous system in both health and disease. In this review we discuss this maturation of single-cell RNA-sequencing over the past decade, and review the different adaptations of the technology that can now be applied both at different scales and for different purposes. We conclude by highlighting how these methods have already led to many exciting discoveries across neuroscience that have furthered our cellular understanding of the neurological disease.