ScType enables fast and accurate cell type identification from spatial transcriptomics data.

SUMMARY: The limited resolution of spatial transcriptomics (ST) assays in the past has led to the development of cell type annotation methods that separate the convolved signal based on available external atlas data. In light of the rapidly increasing resolution of the ST assay technologies, we made available and investigated the performance of a deconvolution-free marker-based cell annotation method called scType. In contrast to existing methods, the spatial application of scType does not require computationally strenuous deconvolution, nor large single-cell reference atlases. We show that scType enables ultra-fast and accurate identification of abundant cell types from ST data, especially when a large enough panel of genes is detected. Examples of such assays are Visium and Slide-seq, which currently offer the best trade-off between high resolution and number of genes detected by the assay for cell type annotation. AVAILABILITY AND IMPLEMENTATION: scType source R and python codes for spatial data are openly available in GitHub (https://github.com/kris-nader/sp-type or https://github.com/kris-nader/sc-type-py). Step-by-step tutorials for R and python spatial data analysis can be found in https://github.com/kris-nader/sp-type and https://github.com/kris-nader/sc-type-py/blob/main/spatial_tutorial.md, respectively.

A community challenge to predict clinical outcomes after immune checkpoint blockade in non-small cell lung cancer.

BACKGROUND: Predictive biomarkers of immune checkpoint inhibitor (ICI) efficacy are currently lacking for non-small cell lung cancer (NSCLC). Here, we describe the results from the Anti-PD-1 Response Prediction DREAM Challenge, a crowdsourced initiative that enabled the assessment of predictive models by using data from two randomized controlled clinical trials (RCTs) of ICIs in first-line metastatic NSCLC. METHODS: Participants developed and trained models using public resources. These were evaluated with data from the CheckMate 026 trial (NCT02041533), according to the model-to-data paradigm to maintain patient confidentiality. The generalizability of the models with the best predictive performance was assessed using data from the CheckMate 227 trial (NCT02477826). Both trials were phase III RCTs with a chemotherapy control arm, which supported the differentiation between predictive and prognostic models. Isolated model containers were evaluated using a bespoke strategy that considered the challenges of handling transcriptome data from clinical trials. RESULTS: A total of 59 teams participated, with 417 models submitted. Multiple predictive models, as opposed to a prognostic model, were generated for predicting overall survival, progression-free survival, and progressive disease status with ICIs. Variables within the models submitted by participants included tumor mutational burden (TMB), programmed death ligand 1 (PD-L1) expression, and gene-expression-based signatures. The best-performing models showed improved predictive power over reference variables, including TMB or PD-L1. CONCLUSIONS: This DREAM Challenge is the first successful attempt to use protected phase III clinical data for a crowdsourced effort towards generating predictive models for ICI clinical outcomes and could serve as a blueprint for similar efforts in other tumor types and disease states, setting a benchmark for future studies aiming to identify biomarkers predictive of ICI efficacy. TRIAL REGISTRATION: CheckMate 026; NCT02041533, registered January 22, 2014. CheckMate 227; NCT02477826, registered June 23, 2015.

SOX9 has distinct roles in the formation and progression of different non-small cell lung cancer histotypes.

The transcription factor SOX9 is a key regulator of multiple developmental processes and is frequently re-expressed in non-small cell lung cancer (NSCLC). Its precise role in the progression of NSCLC histotypes has, however, remained elusive. We show that SOX9 expression relates to poor overall survival and invasive histopathology in human non-mucinous adenocarcinoma and is absent in murine early minimally invasive and low in human in situ adenocarcinoma. Interestingly, despite wide SOX9 expression across advanced NSCLC histotypes, its genetic deletion in the murine KrasG12D ;Lkb1fl/fl model selectively disrupted only the growth of papillary NSCLC, without affecting the initiation of precursor lesions or growth of mucinous or squamous tissue. Spatial tissue phenotyping indicated a requirement of SOX9 expression for the progression of surfactant protein C-expressing progenitor cells, which gave rise to papillary tumours. Intriguingly, while SOX9 expression was dispensable for squamous tissue formation, its loss in fact led to enhanced squamous tumour metastasis, which was associated with altered collagen IV deposition in the basement membrane. Our work therefore demonstrates histopathology-selective roles for SOX9 in NSCLC progression, namely as a promoter for papillary adenocarcinoma progression, but an opposing metastasis-suppressing role in squamous histotype tissue. This attests to a pleiotropic SOX9 function, linked to the cell of origin and microenvironmental tissue contexts. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Spatial aspects of oncogenic signalling determine the response to combination therapy in slice explants from Kras-driven lung tumours.

A key question in precision medicine is how functional heterogeneity in solid tumours informs therapeutic sensitivity. We demonstrate that spatial characteristics of oncogenic signalling and therapy response can be modelled in precision-cut slices from Kras-driven non-small-cell lung cancer with varying histopathologies. Unexpectedly, profiling of in situ tumours demonstrated that signalling stratifies mostly according to histopathology, showing enhanced AKT and SRC activity in adenosquamous carcinoma, and mitogen-activated protein kinase (MAPK) activity in adenocarcinoma. In addition, high intertumour and intratumour variability was detected, particularly of MAPK and mammalian target of rapamycin (mTOR) complex 1 activity. Using short-term treatment of slice explants, we showed that cytotoxic responses to combination MAPK and phosphoinositide 3-kinase-mTOR inhibition correlate with the spatially defined activities of both pathways. Thus, whereas genetic drivers determine histopathology spectra, histopathology-associated and spatially variable signalling activities determine drug sensitivity. Our study is in support of spatial aspects of signalling heterogeneity being considered in clinical diagnostic settings, particularly to guide the selection of drug combinations. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Extendable blocking probe in reverse transcription for analysis of RNA variants with superior selectivity.

Here we provide the first strategy to use a competitive Extendable Blocking Probe (ExBP) for allele-specific priming with superior selectivity at the stage of reverse transcription. In order to analyze highly similar RNA variants, a reverse-transcriptase primer whose sequence matches a specific variant selectively primes only that variant, whereas mismatch priming to the alternative variant is suppressed by virtue of hybridization and subsequent extension of the perfectly matched ExBP on that alternative variant template to form a cDNA-RNA hybrid. This hybrid will render the alternative RNA template unavailable for mismatch priming initiated by the specific primer in a hot-start protocol of reverse transcription when the temperature decreases to a level where such mismatch priming could occur. The ExBP-based reverse transcription assay detected BRAF and KRAS mutations in at least 1000-fold excess of wild-type RNA and detection was linear over a 4-log dynamic range. This novel strategy not only reveals the presence or absence of rare mutations with an exceptionally high selectivity, but also provides a convenient tool for accurate determination of RNA variants in different settings, such as quantification of allele-specific expression.

Temporal dissection of p53 function in vitro and in vivo.

To investigate the functions of the p53 tumor suppressor, we created a new knock-in gene replacement mouse model in which the endogenous Trp53 gene is substituted by one encoding p53ER(TAM), a p53 fusion protein whose function is completely dependent on ectopic provision of 4-hydroxytamoxifen. We show here that both tissues in vivo and cells in vitro derived from such mice can be rapidly toggled between wild-type and p53 knockout states. Using this rapid perturbation model, we define the kinetics, dependence, persistence and reversibility of p53-mediated responses to DNA damage in tissues in vivo and to activation of the Ras oncoprotein and stress in vitro. This is the first example to our knowledge of a new class of genetic model that allows the specific, rapid and reversible perturbation of the function of a single endogenous gene in vivo.

PI3Kß inhibition enhances ALK-inhibitor sensitivity in ALK-rearranged lung cancer.

Treatment with anaplastic lymphoma kinase (ALK) inhibitors significantly improves outcome for non-small-cell lung cancer (NSCLC) patients with ALK-rearranged tumors. However, clinical resistance typically develops over time and, in the majority of cases, resistance mechanisms are ALK-independent. We generated tumor cell cultures from multiple regions of an ALK-rearranged clinical tumor specimen and deployed functional drug screens to identify modulators of ALK-inhibitor response. This identified a role for PI3Kß and EGFR inhibition in sensitizing the response regulating resistance to ALK inhibition. Inhibition of ALK elicited activation of EGFR, and subsequent MAPK and PI3K-AKT pathway reactivation. Sensitivity to ALK targeting was enhanced by inhibition or knockdown of PI3Kß. In ALK-rearranged primary cultures, the combined inhibition of ALK and PI3Kß prevented the EGFR-mediated ALK-inhibitor resistance, and selectively targeted the cancer cells. The combinatorial effect was seen also in the background of TP53 mutations and in epithelial-to-mesenchymal transformed cells. In conclusion, combinatorial ALK- and PI3Kß-inhibitor treatment carries promise as a treatment for ALK-rearranged NSCLC.

Effective intravenous delivery of adenovirus armed with TNFα and IL-2 improves anti-PD-1 checkpoint blockade in non-small cell lung cancer.

Lung cancer remains among the most difficult-to-treat malignancies and is the leading cause of cancer-related deaths worldwide. The introduction of targeted therapies and checkpoint inhibitors has improved treatment outcomes; however, most patients with advanced-stage non-small cell lung cancer (NSCLC) eventually fail these therapies. Therefore, there is a major unmet clinical need for checkpoint refractory/resistant NSCLC. Here, we tested the combination of aPD-1 and adenovirus armed with TNFα and IL-2 (Ad5-CMV-mTNFα/mIL-2) in an immunocompetent murine NSCLC model. Moreover, although local delivery has been standard for virotherapy, treatment was administered intravenously to facilitate clinical translation and putative routine use. We showed that treatment of tumor-bearing animals with aPD-1 in combination with intravenously injected armed adenovirus significantly decreased cancer growth, even in the presence of neutralizing antibodies. We observed an increased frequency of cytotoxic tumor-infiltrating lymphocytes, including tumor-specific cells. Combination treatment led to a decreased percentage of immunosuppressive tumor-associated macrophages and an improvement in dendritic cell maturation. Moreover, we observed expansion of the tumor-specific memory T cell compartment in secondary lymphoid organs in the group that received aPD-1 with the virus. However, although the non-replicative Ad5-CMV-mTNFα/mIL-2 virus allows high transgene expression in the murine model, it does not fully reflect the clinical outcome in humans. Thus, we complemented our findings using NSCLC ex vivo models fully permissive for the TNFα and IL-2- armed oncolytic adenovirus TILT-123. Overall, our data demonstrate the ability of systemically administered adenovirus armed with TNFα and IL-2 to potentiate the anti-tumor efficacy of aPD-1 and warrant further investigation in clinical trials.

Putting transcription repression and protein destruction in pRb's pocket.

The tumor suppressor function of the retinoblastoma protein (pRb) is historically attributed to inhibition of E2F gene transcription. In a recent issue of Nature Cell Biology, Binné and colleagues show that pRB is physically linked to the active anaphase-promoting complex/cyclosome (APC/C) ubiquitin ligase, suggesting that pRB-mediated tumor suppression may also function by directing the local degradation of E2F targets, including Skp2 (Binné et al., 2006).

Prophase I arrest and progression to metaphase I in mouse oocytes are controlled by Emi1-dependent regulation of APC(Cdh1).

Mammalian oocytes are arrested in prophase of the first meiotic division. Progression into the first meiotic division is driven by an increase in the activity of maturation-promoting factor (MPF). In mouse oocytes, we find that early mitotic inhibitor 1 (Emi1), an inhibitor of the anaphase-promoting complex (APC) that is responsible for cyclin B destruction and inactivation of MPF, is present at prophase I and undergoes Skp1-Cul1-F-box/betaTrCP-mediated destruction immediately after germinal vesicle breakdown (GVBD). Exogenous Emi1 or the inhibition of Emi1 destruction in prophase-arrested oocytes leads to a stabilization of cyclin B1-GFP that is sufficient to trigger GVBD. In contrast, the depletion of Emi1 using morpholino oligonucleotides increases cyclin B1-GFP destruction, resulting in an attenuation of MPF activation and a delay of entry into the first meiotic division. Finally, we show that Emi1-dependent effects on meiosis I require the presence of Cdh1. These observations reveal a novel mechanism for the control of entry into the first meiotic division: an Emi1-dependent inhibition of APC(Cdh1).

The oncogenic potential of Kaposi's sarcoma-associated herpesvirus cyclin is exposed by p53 loss in vitro and in vivo.

Expression of the Kaposi's sarcoma-associated herpesvirus (KSHV) cyclin D homolog, K cyclin, is thought to contribute to viral oncogenesis. We show that K cyclin expression in primary cells sensitizes to apoptosis and induces growth arrest, both of which are dependent on p53 but independent of E2F1 or p19(ARF). DNA synthesis, but not cytokinesis, continues in K cyclin-expressing cells, leading to multinucleation and polyploidy. Such polyploid cells exhibit pronounced centrosome amplification and consequent aneuploidy. Our data suggest that K cyclin expression leads to cytokinesis defects and polyploidy, which activates p53. However, in the absence of p53, such cells survive and expand as an aneuploid population. Corroborating these findings, in vivo Emu; K cyclin expression cooperates with p53 loss in the induction of lymphomas.

Protocol to utilize fresh uncultured human lung tumor cells for personalized functional diagnostics.

Drug sensitivity data acquired from solid tumor-derived cultures are often unsuitable for personalized treatment guidance due to the lengthy turnaround time. Here, we present a protocol for determining ex vivo drug sensitivities using fresh uncultured human lung tumor-derived EpCAM+ epithelial cells (FUTCs). We describe steps for drug testing in FUTCs to identify tumor cell-selective single or combination therapy in 72 h of sample processing. The FUTC-based approach can also be used to predict in vivo resistance to known targeted therapies. For complete details on the use and execution of this protocol, please refer to Talwelkar et al. (2021).

Development of an adenosquamous carcinoma histopathology - selective lung metastasis model.

Preclinical tumor models with native tissue microenvironments provide essential tools to understand how heterogeneous tumor phenotypes relate to drug response. Here we present syngeneic graft models of aggressive, metastasis-prone histopathology-specific NSCLC tumor types driven by KRAS mutation and loss of LKB1 (KL): adenosquamous carcinoma (ASC) and adenocarcinoma (AC). We show that subcutaneous injection of primary KL; ASC cells results in squamous cell carcinoma (SCC) tumors with high levels of stromal infiltrates, lacking the source heterogeneous histotype. Despite forming subcutaneous tumors, intravenously injected KL;AC cells were unable to form lung tumors. In contrast, intravenous injection of KL;ASC cells leads to their lung re-colonization and lesions recapitulating the mixed AC and SCC histopathology, tumor immune suppressive microenvironment and oncogenic signaling profile of source tumors, demonstrating histopathology-selective phenotypic dominance over genetic drivers. Pan-ERBB inhibition increased survival, while selective ERBB1/EGFR inhibition did not, suggesting a role of the ERBB network crosstalk in resistance to ERBB1/EGFR. This immunocompetent NSCLC lung colonization model hence phenocopies key properties of the metastasis-prone ASC histopathology, and serves as a preclinical model to dissect therapy responses and metastasis-associated processes.

Inactivation of AMPK Leads to Attenuation of Antigen Presentation and Immune Evasion in Lung Adenocarcinoma.

PURPOSE: Mutations in STK11 (LKB1) occur in 17% of lung adenocarcinoma (LUAD) and drive a suppressive (cold) tumor immune microenvironment (TIME) and resistance to immunotherapy. The mechanisms underpinning the establishment and maintenance of a cold TIME in LKB1-mutant LUAD remain poorly understood. In this study, we investigated the role of the LKB1 substrate AMPK in immune evasion in human non-small cell lung cancer (NSCLC) and mouse models and explored the mechanisms involved. EXPERIMENTAL DESIGN: We addressed the role of AMPK in immune evasion in NSCLC by correlating AMPK phosphorylation and immune-suppressive signatures and by deleting AMPKα1 (Prkaa1) and AMPKα2 (Prkaa2) in a KrasG12D -driven LUAD. Furthermore, we dissected the molecular mechanisms involved in immune evasion by comparing gene-expression signatures, AMPK activity, and immune infiltration in mouse and human LUAD and gain or loss-of-function experiments with LKB1- or AMPK-deficient cell lines. RESULTS: Inactivation of both AMPKα1 and AMPKα2 together with Kras activation accelerated tumorigenesis and led to tumors with reduced infiltration of CD8+/CD4+ T cells and gene signatures associated with a suppressive TIME. These signatures recapitulate those in Lkb1-deleted murine LUAD and in LKB1-deficient human NSCLC. Interestingly, a similar signature is noted in human NSCLC with low AMPK activity. In mechanistic studies, we find that compromised LKB1 and AMPK activity leads to attenuated antigen presentation in both LUAD mouse models and human NSCLC. CONCLUSIONS: The results provide evidence that the immune evasion noted in LKB1-inactivated lung cancer is due to subsequent inactivation of AMPK and attenuation of antigen presentation.

Functional diagnostics using fresh uncultured lung tumor cells to guide personalized treatments.

Functional profiling of a cancer patient's tumor cells holds potential to tailor personalized cancer treatment. Here, we report the utility of fresh uncultured tumor-derived EpCAM+ epithelial cells (FUTCs) for ex vivo drug-response interrogation. Analysis of murine Kras mutant FUTCs demonstrates pharmacological and adaptive signaling profiles comparable to subtype-matched cultured cells. By applying FUTC profiling on non-small-cell lung cancer patient samples, we report robust drug-response data in 19 of 20 cases, with cells exhibiting targeted drug sensitivities corresponding to their oncogenic drivers. In one of these cases, an EGFR mutant lung adenocarcinoma patient refractory to osimertinib, FUTC profiling is used to guide compassionate treatment. FUTC profiling identifies selective sensitivity to disulfiram and the combination of carboplatin plus etoposide, and the patient receives substantial clinical benefit from treatment with these agents. We conclude that FUTC profiling provides a robust, rapid, and actionable assessment of personalized cancer treatment options.

Oncogenic Ras Disrupts Epithelial Integrity by Activating the Transmembrane Serine Protease Hepsin.

Ras proteins play a causal role in human cancer by activating multiple pathways that promote cancer growth and invasion. However, little is known about how Ras induces the first diagnostic features of invasion in solid tumors, including loss of epithelial integrity and breaching of the basement membrane (BM). In this study, we found that oncogenic Ras strongly promotes the activation of hepsin, a member of the hepsin/TMPRSS type II transmembrane serine protease family. Mechanistically, the Ras-dependent hepsin activation was mediated via Raf-MEK-ERK signaling, which controlled hepsin protein stability through the heat shock transcription factor-1 stress pathway. In Ras-transformed three-dimensional mammary epithelial culture, ablation of hepsin restored desmosomal cell-cell junctions, hemidesmosomes, and BM integrity and epithelial cohesion. In tumor xenografts harboring mutant KRas, silencing of hepsin increased local invasion concomitantly with accumulation of collagen IV. These findings suggest that hepsin is a critical protease for Ras-dependent tumorigenesis, executing cell-cell and cell-matrix pathologies important for early tumor dissemination. SIGNIFICANCE: These findings identify the cell-surface serine protease hepsin as a potential therapeutic target for its role in oncogenic Ras-mediated deregulation of epithelial cell-cell and cell-matrix interactions and cohesion of epithelial structure.

Receptor Tyrosine Kinase Signaling Networks Define Sensitivity to ERBB Inhibition and Stratify Kras-Mutant Lung Cancers.

Most non-small cell lung cancers (NSCLC) contain nontargetable mutations, including KRAS, TP53, or STK11/LKB1 alterations. By coupling ex vivo drug sensitivity profiling with in vivo drug response studies, we aimed to identify drug vulnerabilities for these NSCLC subtypes. Primary adenosquamous carcinoma (ASC) or adenocarcinoma (AC) cultures were established from KrasG12D/+;Lkb1fl/fl (KL) tumors or AC cultures from KrasG12D/+;p53fl/fl (KP) tumors. Although p53-null cells readily propagated as conventional cultures, Lkb1-null cells required conditional reprograming for establishment. Drug response profiling revealed short-term response to MEK inhibition, yet long-term clonogenic assays demonstrated resistance, associated with sustained or adaptive activation of receptor tyrosine kinases (RTK): activation of ERBBs in KL cultures, or FGFR in AC cultures. Furthermore, pan-ERBB inhibition reduced the clonogenicity of KL cultures, which was exacerbated by combinatorial MEK inhibition, whereas combinatorial MEK and FGFR inhibition suppressed clonogenicity of AC cultures. Importantly, in vivo studies confirmed KL-selective sensitivity to pan-ERBB inhibition, which correlated with high ERBB ligand expression and activation of ERBB receptors, implying that ERBB network activity may serve as a predictive biomarker of drug response. Interestingly, in human NSCLCs, phosphorylation of EGFR or ERBB3 was frequently detected in ASCs and squamous cell carcinomas. We conclude that analysis of in situ ERBB signaling networks in conjunction with ex vivo drug response profiling and biochemical dissection of adaptive RTK activities may serve as a valid diagnostic approach to identify tumors sensitive to ERBB network inhibition.

Spa-RQ: an Image Analysis Tool to Visualise and Quantify Spatial Phenotypes Applied to Non-Small Cell Lung Cancer.

To facilitate analysis of spatial tissue phenotypes, we created an open-source tool package named 'Spa-RQ' for 'Spatial tissue analysis: image Registration & Quantification'. Spa-RQ contains software for image registration (Spa-R) and quantitative analysis of DAB staining overlap (Spa-Q). It provides an easy-to-implement workflow for serial sectioning and staining as an alternative to multiplexed techniques. To demonstrate Spa-RQ's applicability, we analysed the spatial aspects of oncogenic KRAS-related signalling activities in non-small cell lung cancer (NSCLC). Using Spa-R in conjunction with ImageJ/Fiji, we first performed annotation-guided tumour-by-tumour phenotyping using multiple signalling markers. This analysis showed histopathology-selective activation of PI3K/AKT and MAPK signalling in Kras mutant murine tumours, as well as high p38MAPK stress signalling in p53 null murine NSCLC. Subsequently, Spa-RQ was applied to measure the co-activation of MAPK, AKT, and their mutual effector mTOR pathway in individual tumours. Both murine and clinical NSCLC samples could be stratified into 'MAPK/mTOR', 'AKT/mTOR', and 'Null' signature subclasses, suggesting mutually exclusive MAPK and AKT signalling activities. Spa-RQ thus provides a robust and easy to use tool that can be employed to identify spatially-distributed tissue phenotypes.

Establishment and Analysis of Tumor Slice Explants As a Prerequisite for Diagnostic Testing.

Organotypic primary tissue explant cultures, which include precision-cut slices, represent the three-dimensional (3-D) tissue architecture as well as the multicellular interactions of native tissue. Tissue slices immediately cut from freshly resected tumors preserve spatial aspects of intratumor heterogeneity (ITH), thus making them useful surrogates of in vivo biology. Careful optimization of tissue slice preparation and cultivation conditions is fundamental for the predictive diagnostic potential of tumor slice explants. In this regard, murine models are valuable, as these provide a consistent flow of tumor material to perform replicate and reproducible experiments. This protocol describes the culturing of murine lung tumor-derived tissue slices using a rotating incubation unit, a system that enables intermittent exposure of tissues to oxygen and nutrients. Our previous work showed that the use of rotating incubation units improves the viability of tissue compared to other culture methods, particularly floating slices and stagnant filter supports. Here, we further show that slice thickness influences the viability of cultured slices, suggesting that optimization of slice thickness should be done for different tissue types. Pronounced ITH in relevant oncogenic functions, such as signaling activities, stromal cell infiltration or expression of differentiation markers, necessitates evaluation of adjacent tissue slices for the expression of markers altered by drug treatment or cultivation itself. In summary, this protocol describes the generation of murine lung tumor slices and their culture on a rotating incubation unit and demonstrates how slices should be systematically analyzed for the expression of heterogeneous tissue markers, as a prerequisite prior to drug response studies.