Spatially resolved whole transcriptome profiling in human and mouse tissue using Digital Spatial Profiling.

Emerging spatial profiling technology has enabled high-plex molecular profiling in biological tissues, preserving the spatial and morphological context of gene expression. Here, we describe expanding the chemistry for the Digital Spatial Profiling platform to quantify whole transcriptomes in human and mouse tissues using a wide range of spatial profiling strategies and sample types. We designed multiplexed in situ hybridization probes targeting the protein-coding genes of the human and mouse transcriptomes, referred to as the human or mouse Whole Transcriptome Atlas (WTA). Human and mouse WTAs were validated in cell lines for concordance with orthogonal gene expression profiling methods in regions ranging from ∼10-500 cells. By benchmarking against bulk RNA-seq and fluorescence in situ hybridization, we show robust transcript detection down to ∼100 transcripts per region. To assess the performance of WTA across tissue and sample types, we applied WTA to biological questions in cancer, molecular pathology, and developmental biology. Spatial profiling with WTA detected expected gene expression differences between tumor and tumor microenvironment, identified disease-specific gene expression heterogeneity in histological structures of the human kidney, and comprehensively mapped transcriptional programs in anatomical substructures of nine organs in the developing mouse embryo. Digital Spatial Profiling technology with the WTA assays provides a flexible method for spatial whole transcriptome profiling applicable to diverse tissue types and biological contexts.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Biomarkers Associated with Beneficial PD-1 Checkpoint Blockade in Non-Small Cell Lung Cancer (NSCLC) Identified Using High-Plex Digital Spatial Profiling.

PURPOSE: Only a minority of patients with advanced non-small cell lung cancer (NSCLC) truly benefits from single-agent PD-1 checkpoint blockade, and more robust predictive biomarkers are needed. EXPERIMENTAL DESIGN: We assessed tumor samples from 67 immunotherapy-treated NSCLC cases represented in a tissue microarray, 53 of whom had pretreatment samples and received monotherapy. Using GeoMx Digital Spatial Profiling System (NanoString Technologies), we quantified 39 immune parameters simultaneously in four tissue compartments defined by fluorescence colocalization [tumor (panCK+), leucocytes (CD45+), macrophages (CD68+), and nonimmune stroma]. RESULTS: A total of 156 protein variables were generated per case. In the univariate unadjusted analysis, we found 18 markers associated with outcome in spatial context, five of which remained significant after multiplicity adjustment. In the multivariate analysis, high levels of CD56 and CD4 measured in the CD45 compartment were the only markers that were predictive for all clinical outcomes, including progression-free survival (PFS, HR: 0.24, P = 0.006; and HR: 0.31, P = 0.011, respectively), and overall survival (OS, HR: 0.26, P = 0.014; and HR: 0.23, P = 0.007, respectively). Then, using an orthogonal method based on multiplex immunofluorescence and cell counting (inForm), we validated that high CD56+ immune cell counts in the stroma were associated with PFS and OS in the same cohort. CONCLUSIONS: This pilot scale discovery study shows the potential of the digital spatial profiling technology in the identification of spatially informed biomarkers of response to PD-1 checkpoint blockade in NSCLC. We identified a number of relevant candidate immune predictors in spatial context that deserve validation in larger independent cohorts.

Digital quantitative assessment of PD-L1 using digital spatial profiling.

The assessment of programmed death 1 ligand 1 (PD-L1) expression by Immunohistochemistry (IHC) is the US Food and Drug Administration (FDA)-approved predictive marker to select responders to checkpoint blockade anti-PD-1/PD-L1 axis immunotherapies. Different PD-L1 immunohistochemistry (IHC) assays use different antibodies and different scoring methods in tumor cells and immune cells. Multiple studies have compared the performance of these assays with variable results. Here, we investigate an alternative method for assessment of PD-L1 using a new technology known as digital spatial profiling. We use a previously described standardization tissue microarray (TMA) to assess the accuracy of the method and compare digital spatial profiler (DSP) to each FDA-approved PD-L1 assays, one LDT assay and three quantitative fluorescence assays. The standardized cell line Index tissue microarray contains 10 isogenic cells lines in triplicates expressing various ranges of PD-L1. The dynamic range of PD-L1 digital counts was measured in the ten cell lines on the Index TMA using the GeoMx DSP assay and read on the nCounter platform. The digital method shows very high correlation with immunohistochemistry scored with quantitative software and with quantitative fluorescence. High correlation of PD-L1 digital DSP counts were seen between rows on the same Index TMA. Finally, experiments from two Index TMAs showed reproducibility of DSP counts were independent of variable slide storage time over a three-week period after antibody labeling but before collection of cleaved tags. In summary, DSP appears to have quantitative potential comparable to quantitative immunohistochemistry. It is possible that this technology could be used as a PD-L1 protein measurement system for companion diagnostic testing for immune therapy.

Gene expression profiling distinguishes prefibrotic from overtly fibrotic myeloproliferative neoplasms and identifies disease subsets with distinct inflammatory signatures.

The Philadelphia chromosome-negative myeloproliferative neoplasms (MPN) share similar molecular characteristics in that they frequently harbor hotspot mutations in JAK2, CALR or MPL, leading to activated JAK/STAT signaling. However, these MPN have distinct symptoms, morphology, and natural histories, including different tendencies to progress to fibrosis. Although the role of inflammation in tissue fibrosis is well recognized, inflammatory gene expression in bone marrows involved by MPN has been understudied. We analyzed the expression of inflammatory genes by directly measuring RNA transcript abundance in bone marrow biopsies of 108 MPN patients. Unsupervised analyses identified gene expression patterns that distinguish prefibrotic (grade 1-2) MPN from overtly fibrotic (grade 2-3) MPN with high sensitivity and specificity in two independent cohorts. Furthermore, prefibrotic and overtly fibrotic MPN are separable into subsets with different activities in biological pathways linked to inflammation, including cytokines, chemokines, interferon response, and toll-like receptor signaling. In conclusion, this study demonstrates that MPN with overt fibrosis is associated with significant inflammatory gene upregulation in the bone marrow, revealing potential roles for multiple pro-inflammatory signaling networks in the development of myelofibrosis and suggesting potential therapeutic mechanisms to alleviate this process in the bone marrow microenvironment.

Seven-Year Follow-Up Analysis of Adjuvant Paclitaxel and Trastuzumab Trial for Node-Negative, Human Epidermal Growth Factor Receptor 2-Positive Breast Cancer.

PURPOSE: The Adjuvant Paclitaxel and Trastuzumab trial was designed to address treatment of patients with small human epidermal growth factor receptor 2 (HER2)-positive breast cancer. The primary analysis of the Adjuvant Paclitaxel and Trastuzumab trial demonstrated a 3-year disease-free survival (DFS) of 98.7%. In this planned secondary analysis, we report longer-term outcomes and exploratory results to characterize the biology of small HER2-positive tumors and genetic factors that may predispose to paclitaxel-induced peripheral neuropathy (TIPN). PATIENTS AND METHODS: In this phase II study, patients with HER2-positive breast cancer with tumors 3 cm or smaller and negative nodes received paclitaxel (80 mg/m2) with trastuzumab for 12 weeks, followed by trastuzumab for 9 months. The primary end point was DFS. Recurrence-free interval (RFI), breast cancer-specific survival, and overall survival (OS) were also analyzed. In an exploratory analysis, intrinsic subtyping by PAM50 (Prosigna) and calculation of the risk of recurrence score were performed on the nCounter analysis system on archival tissue. Genotyping was performed to investigate TIPN. RESULTS: A total of 410 patients were enrolled from October 2007 to September 2010. After a median follow-up of 6.5 years, there were 23 DFS events. The 7-year DFS was 93% (95% CI, 90.4 to 96.2) with four (1.0%) distant recurrences, 7-year OS was 95% (95% CI, 92.4 to 97.7), and 7-year RFI was 97.5% (95% CI, 95.9 to 99.1). PAM50 analyses (n = 278) showed that most tumors were HER2-enriched (66%), followed by luminal B (14%), luminal A (13%), and basal-like (8%). Genotyping (n = 230) identified one single-nucleotide polymorphism, rs3012437, associated with an increased risk of TIPN in patients with grade 2 or greater TIPN (10.4%). CONCLUSION: With longer follow-up, adjuvant paclitaxel and trastuzumab is associated with excellent long-term outcomes. Distribution of PAM50 intrinsic subtypes in small HER2-positive tumors is similar to that previously reported for larger tumors.