COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

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.

Highly Multiplexed Spatially Resolved Proteomic and Transcriptional Profiling of the Glioblastoma Microenvironment Using Archived Formalin-Fixed Paraffin-Embedded Specimens.

Glioblastoma is a heterogeneous tumor for which effective treatment options are limited and often insufficient. Few studies have examined the intratumoral transcriptional and proteomic heterogeneity of the glioblastoma microenvironment to characterize the spatial distribution of potential molecular and cellular therapeutic immunooncology targets. We applied an integrated multimodal approach comprised of NanoString GeoMx Digital Spatial Profiling, single-cell RNA-seq (scRNA-seq), and expert neuropathologic assessment to characterize archival formalin-fixed paraffin-embedded glioblastoma specimens. Clustering analysis and spatial cluster maps highlighted the intratumoral heterogeneity of each specimen. Mixed cell deconvolution analysis revealed that neoplastic and vascular cells were the prominent cell types throughout each specimen, with macrophages, oligodendrocyte precursors, neurons, astrocytes, and oligodendrocytes present in lower abundance and illustrated the regional distribution of the respective cellular enrichment scores. The spatial resolution of the actionable immunotherapeutic landscape showed that robust B7H3 gene and protein expression was broadly distributed throughout each specimen and identified STING and VISTA as potential targets. Lastly, we uncovered remarkable variability in VEGFA expression and discovered unanticipated associations between VEGFA, endothelial cell markers, hypoxia, and the expression of immunoregulatory genes, indicative of regionally distinct immunosuppressive microdomains. This work provides an early demonstration of the ability of an integrated panel-based spatial biology approach to characterize and quantify the intrinsic molecular heterogeneity of the glioblastoma microenvironment.

Multiplex Detection of Clinically Relevant Mutations in Liquid Biopsies of Cancer Patients Using a Hybridization-Based Platform.

BACKGROUND: With the advent of precision oncology, liquid biopsies are quickly gaining acceptance in the clinical setting. However, in some cases, the amount of DNA isolated is insufficient for Next-Generation Sequencing (NGS) analysis. The nCounter platform could be an alternative, but it has never been explored for detection of clinically relevant alterations in fluids. METHODS: Circulating-free DNA (cfDNA) was purified from blood, cerebrospinal fluid, and ascites of patients with cancer and analyzed with the nCounter 3 D Single Nucleotide Variant (SNV) Solid Tumor Panel, which allows for detection of 97 driver mutations in 24 genes. RESULTS: Validation experiments revealed that the nCounter SNV panel could detect mutations at allelic fractions of 0.02-2% in samples with ≥5 pg mutant DNA/µL. In a retrospective analysis of 70 cfDNAs from patients with cancer, the panel successfully detected EGFR, KRAS, BRAF, PIK3CA, and NRAS mutations when compared with previous genotyping in the same liquid biopsies and paired tumor tissues [Cohen kappa of 0.96 (CI = 0.92-1.00) and 0.90 (CI = 0.74-1.00), respectively]. In a prospective study including 91 liquid biopsies from patients with different malignancies, 90 yielded valid results with the SNV panel and mutations in EGFR, KRAS, BRAF, PIK3CA, TP53, NFE2L2, CTNNB1, ALK, FBXW7, and PTEN were found. Finally, serial liquid biopsies from a patient with NSCLC revealed that the semiquantitative results of the mutation analysis by the SNV panel correlated with the evolution of the disease. CONCLUSIONS: The nCounter platform requires less DNA than NGS and can be employed for routine mutation testing in liquid biopsies of patients with cancer.

Implementation of a Multiplex and Quantitative Proteomics Platform for Assessing Protein Lysates Using DNA-Barcoded Antibodies.

Molecular analysis of tumors forms the basis for personalized cancer medicine and increasingly guides patient selection for targeted therapy. Future opportunities for personalized medicine are highlighted by the measurement of protein expression levels via immunohistochemistry, protein arrays, and other approaches; however, sample type, sample quantity, batch effects, and "time to result" are limiting factors for clinical application. Here, we present a development pipeline for a novel multiplexed DNA-labeled antibody platform which digitally quantifies protein expression from lysate samples. We implemented a rigorous validation process for each antibody and show that the platform is amenable to multiple protocols covering nitrocellulose and plate-based methods. Results are highly reproducible across technical and biological replicates, and there are no observed "batch effects" which are common for most multiplex molecular assays. Tests from basal and perturbed cancer cell lines indicate that this platform is comparable to orthogonal proteomic assays such as Reverse-Phase Protein Array, and applicable to measuring the pharmacodynamic effects of clinically-relevant cancer therapeutics. Furthermore, we demonstrate the potential clinical utility of the platform with protein profiling from breast cancer patient samples to identify molecular subtypes. Together, these findings highlight the potential of this platform for enhancing our understanding of cancer biology in a clinical translation setting.

Continuously tunable nucleic acid hybridization probes.

In silico-designed nucleic acid probes and primers often do not achieve favorable specificity and sensitivity tradeoffs on the first try, and iterative empirical sequence-based optimization is needed, particularly in multiplexed assays. We present a novel, on-the-fly method of tuning probe affinity and selectivity by adjusting the stoichiometry of auxiliary species, which allows for independent and decoupled adjustment of the hybridization yield for different probes in multiplexed assays. Using this method, we achieved near-continuous tuning of probe effective free energy. To demonstrate our approach, we enforced uniform capture efficiency of 31 DNA molecules (GC content, 0-100%), maximized the signal difference for 11 pairs of single-nucleotide variants and performed tunable hybrid capture of mRNA from total RNA. Using the Nanostring nCounter platform, we applied stoichiometric tuning to simultaneously adjust yields for a 24-plex assay, and we show multiplexed quantitation of RNA sequences and variants from formalin-fixed, paraffin-embedded samples.

Modeling bi-modality improves characterization of cell cycle on gene expression in single cells.

Advances in high-throughput, single cell gene expression are allowing interrogation of cell heterogeneity. However, there is concern that the cell cycle phase of a cell might bias characterizations of gene expression at the single-cell level. We assess the effect of cell cycle phase on gene expression in single cells by measuring 333 genes in 930 cells across three phases and three cell lines. We determine each cell's phase non-invasively without chemical arrest and use it as a covariate in tests of differential expression. We observe bi-modal gene expression, a previously-described phenomenon, wherein the expression of otherwise abundant genes is either strongly positive, or undetectable within individual cells. This bi-modality is likely both biologically and technically driven. Irrespective of its source, we show that it should be modeled to draw accurate inferences from single cell expression experiments. To this end, we propose a semi-continuous modeling framework based on the generalized linear model, and use it to characterize genes with consistent cell cycle effects across three cell lines. Our new computational framework improves the detection of previously characterized cell-cycle genes compared to approaches that do not account for the bi-modality of single-cell data. We use our semi-continuous modelling framework to estimate single cell gene co-expression networks. These networks suggest that in addition to having phase-dependent shifts in expression (when averaged over many cells), some, but not all, canonical cell cycle genes tend to be co-expressed in groups in single cells. We estimate the amount of single cell expression variability attributable to the cell cycle. We find that the cell cycle explains only 5%-17% of expression variability, suggesting that the cell cycle will not tend to be a large nuisance factor in analysis of the single cell transcriptome.

Spatially resolved analysis of pancreatic cancer identifies therapy-associated remodeling of the tumor microenvironment.

In combination with cell-intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell-cell interactions in human pancreatic cancer associated with neoadjuvant chemotherapy and radiotherapy. We developed spatially constrained optimal transport interaction analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand-receptor gene expression. Our results uncovered a marked change in ligand-receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid coculture system. We identified enrichment in interleukin-6 family signaling that functionally confers resistance to chemotherapy. Overall, this study demonstrates that characterization of the tumor microenvironment using single-cell spatial transcriptomics allows for the identification of molecular interactions that may play a role in the emergence of therapeutic resistance and offers a spatially based analysis framework that can be broadly applied to other contexts.

An information theoretic approach to detecting spatially varying genes.

A key step in spatial transcriptomics is identifying genes with spatially varying expression patterns. We adopt an information theoretic perspective to this problem by equating the degree of spatial coherence with the Jensen-Shannon divergence between pairs of nearby cells and pairs of distant cells. To avoid the notoriously difficult problem of estimating information theoretic divergences, we use modern approximation techniques to implement a computationally efficient algorithm designed to scale with in situ spatial transcriptomics technologies. In addition to being highly scalable, we show that our method, which we call maximization of spatial information (Maxspin), improves accuracy across several spatial transcriptomics platforms and a variety of simulations when compared with a variety of state-of-the-art methods. To further demonstrate the method, we generated in situ spatial transcriptomics data in a renal cell carcinoma sample using the CosMx Spatial Molecular Imager and used Maxspin to reveal novel spatial patterns of tumor cell gene expression.

Ultra High-plex Spatial Proteogenomic Investigation of Giant Cell Glioblastoma Multiforme Immune Infiltrates Reveals Distinct Protein and RNA Expression Profiles.

A deeper understanding of complex biological processes, including tumor development and immune response, requires ultra high-plex, spatial interrogation of multiple "omes". Here we present the development and implementation of a novel spatial proteogenomic (SPG) assay on the GeoMx Digital Spatial Profiler platform with next-generation sequencing readout that enables ultra high-plex digital quantitation of proteins (>100-plex) and RNA (whole transcriptome, >18,000-plex) from a single formalin-fixed paraffin-embedded (FFPE) sample. This study highlighted the high concordance, R > 0.85 and <15% change in sensitivity between the SPG assay and the single-analyte assays on various cell lines and tissues from human and mouse. Furthermore, we demonstrate that the SPG assay was reproducible across multiple users. When used in conjunction with advanced cellular neighborhood segmentation, distinct immune or tumor RNA and protein targets were spatially resolved within individual cell subpopulations in human colorectal cancer and non-small cell lung cancer. We used the SPG assay to interrogate 23 different glioblastoma multiforme (GBM) samples across four pathologies. The study revealed distinct clustering of both RNA and protein based on pathology and anatomic location. The in-depth investigation of giant cell glioblastoma multiforme (gcGBM) revealed distinct protein and RNA expression profiles compared with that of the more common GBM. More importantly, the use of spatial proteogenomics allowed simultaneous interrogation of critical protein posttranslational modifications alongside whole transcriptomic profiles within the same distinct cellular neighborhoods. Significance: We describe ultra high-plex spatial proteogenomics; profiling whole transcriptome and high-plex proteomics on a single FFPE tissue section with spatial resolution. Investigation of gcGBM versus GBM revealed distinct protein and RNA expression profiles.

Multiplatform Analysis of Intratumoral PTEN Heterogeneity in Melanoma.

Loss of protein expression of the tumor suppressor PTEN is associated with increased cancer aggressiveness, decreased tumor immune infiltration, and resistance to immune and targeted therapies in melanoma. We assessed a unique cohort of eight melanoma samples with focal loss of PTEN protein expression to understand the features and mechanisms of PTEN loss in this disease. We compared the PTEN-negative (PTEN[-]) areas to their adjacent PTEN-positive (PTEN[+]) areas using DNA sequencing, DNA methylation, RNA expression, digital spatial profiling, and immunohistochemical platforms. Variations or homozygous deletions of PTEN were identified in PTEN(-) areas that were not detected in the adjacent PTEN(+) areas in three cases (37.5%), but no clear genomic or DNA methylation basis for loss was identified in the remaining PTEN(-) samples. RNA expression data from two independent platforms identified a consistent increase in chromosome segregation gene expression in PTEN(-) versus adjacent PTEN(+) areas. Proteomic analysis showed a relative paucity of tumor-infiltrating lymphocytes in PTEN(-) versus adjacent PTEN(+) areas. The findings add to our understanding of potential molecular intratumoral heterogeneity in melanoma and the features associated with the loss of PTEN protein in this disease.

Organ-specific immunity: A tissue analysis framework for investigating local immune responses to SARS-CoV-2.

Local immune activation at mucosal surfaces, mediated by mucosal lymphoid tissues, is vital for effective immune responses against pathogens. While pathogens like severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can spread to multiple organs, patients with coronavirus disease 2019 (COVID-19) primarily experience inflammation and damage in their lungs. To investigate this apparent organ-specific immune response, we develop an analytical framework that recognizes the significance of mucosal lymphoid tissues. This framework combines histology, immunofluorescence, spatial transcript profiling, and mathematical modeling to identify cellular and gene expression differences between the lymphoid tissues of the lung and the gut and predict the determinants of those differences. Our findings indicate that mucosal lymphoid tissues are pivotal in organ-specific immune response to SARS-CoV-2, mediating local inflammation and tissue damage and contributing to immune dysfunction. The framework developed here has potential utility in the study of long COVID and may streamline biomarker discovery and treatment design for diseases with differential pathologies at the organ level.

Therapy-associated remodeling of pancreatic cancer revealed by single-cell spatial transcriptomics and optimal transport analysis.

In combination with cell intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged high-plex single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell-cell interactions in human pancreatic cancer associated with specific malignant subtypes and neoadjuvant chemotherapy/radiotherapy. We developed Spatially Constrained Optimal Transport Interaction Analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand-receptor gene expression. Our results uncovered a marked change in ligand-receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid co-culture system. Overall, this study demonstrates that characterization of the tumor microenvironment using high-plex single-cell spatial transcriptomics allows for identification of molecular interactions that may play a role in the emergence of chemoresistance and establishes a translational spatial biology paradigm that can be broadly applied to other malignancies, diseases, and treatments.

Advances in mixed cell deconvolution enable quantification of cell types in spatial transcriptomic data.

Mapping cell types across a tissue is a central concern of spatial biology, but cell type abundance is difficult to extract from spatial gene expression data. We introduce SpatialDecon, an algorithm for quantifying cell populations defined by single cell sequencing within the regions of spatial gene expression studies. SpatialDecon incorporates several advancements in gene expression deconvolution. We propose an algorithm harnessing log-normal regression and modelling background, outperforming classical least-squares methods. We compile cell profile matrices for 75 tissue types. We identify genes whose minimal expression by cancer cells makes them suitable for immune deconvolution in tumors. Using lung tumors, we create a dataset for benchmarking deconvolution methods against marker proteins. SpatialDecon is a simple and flexible tool for mapping cell types in spatial gene expression studies. It obtains cell abundance estimates that are spatially resolved, granular, and paired with highly multiplexed gene expression data.

High-Plex Spatial RNA Profiling Reveals Cell Type‒Specific Biomarker Expression during Melanoma Development.

Early diagnosis of melanoma is critical for improved survival. However, the biomarkers of early melanoma evolution and their origin within the tumor and its microenvironment, including the keratinocytes, are poorly defined. To address this, we used spatial transcript profiling that maintains the morphological tumor context to measure the expression of >1,000 RNAs in situ in patient-derived formalin-fixed, paraffin-embedded tissue sections in primary melanoma and melanocytic nevi. We profiled 134 regions of interest (each 200 µm in diameter) enriched in melanocytes, neighboring keratinocytes, or immune cells. This approach captured distinct expression patterns across cell types and tumor types during melanoma development. Unexpectedly, we discovered that S100A8 is expressed by keratinocytes within the tumor microenvironment during melanoma growth. Immunohistochemistry of 252 tumors showed prominent keratinocyte-derived S100A8 expression in melanoma but not in benign tumors and confirmed the same pattern for S100A8's binding partner S100A9, suggesting that injury to the epidermis may be an early and readily detectable indicator of melanoma development. Together, our results establish a framework for high-plex, spatial, and cell type‒specific resolution of gene expression in archival tissue applicable to the development of biomarkers and characterization of tumor microenvironment interactions in tumor evolution.

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.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients.

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

High-plex imaging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging.

Resolving the spatial distribution of RNA and protein in tissues at subcellular resolution is a challenge in the field of spatial biology. We describe spatial molecular imaging, a system that measures RNAs and proteins in intact biological samples at subcellular resolution by performing multiple cycles of nucleic acid hybridization of fluorescent molecular barcodes. We demonstrate that spatial molecular imaging has high sensitivity (one or two copies per cell) and very low error rate (0.0092 false calls per cell) and background (~0.04 counts per cell). The imaging system generates three-dimensional, super-resolution localization of analytes at ~2 million cells per sample. Cell segmentation is morphology based using antibodies, compatible with formalin-fixed, paraffin-embedded samples. We measured multiomic data (980 RNAs and 108 proteins) at subcellular resolution in formalin-fixed, paraffin-embedded tissues (nonsmall cell lung and breast cancer) and identified >18 distinct cell types, ten unique tumor microenvironments and 100 pairwise ligand-receptor interactions. Data on >800,000 single cells and ~260 million transcripts can be accessed at http://nanostring.com/CosMx-dataset .

Inter- and intra-tumor heterogeneity of metastatic prostate cancer determined by digital spatial gene expression profiling.

Metastatic prostate cancer (mPC) comprises a spectrum of diverse phenotypes. However, the extent of inter- and intra-tumor heterogeneity is not established. Here we use digital spatial profiling (DSP) technology to quantitate transcript and protein abundance in spatially-distinct regions of mPCs. By assessing multiple discrete areas across multiple metastases, we find a high level of intra-patient homogeneity with respect to tumor phenotype. However, there are notable exceptions including tumors comprised of regions with high and low androgen receptor (AR) and neuroendocrine activity. While the vast majority of metastases examined are devoid of significant inflammatory infiltrates and lack PD1, PD-L1 and CTLA4, the B7-H3/CD276 immune checkpoint protein is highly expressed, particularly in mPCs with high AR activity. Our results demonstrate the utility of DSP for accurately classifying tumor phenotype, assessing tumor heterogeneity, and identifying aspects of tumor biology involving the immunological composition of metastases.

CTLA-4 blockade and interferon-α induce proinflammatory transcriptional changes in the tumor immune landscape that correlate with pathologic response in melanoma.

Patients with locally/regionally advanced melanoma were treated with neoadjuvant combination immunotherapy with high-dose interferon α-2b (HDI) and ipilimumab in a phase I clinical trial. Tumor specimens were obtained prior to the initiation of neoadjuvant therapy, at the time of surgery and progression if available. In this study, gene expression profiles of tumor specimens (N = 27) were investigated using the NanoString nCounter® platform to evaluate associations with clinical outcomes (pathologic response, radiologic response, relapse-free survival (RFS), and overall survival (OS)) and define biomarkers associated with tumor response. The Tumor Inflammation Signature (TIS), an 18-gene signature that enriches for response to Programmed cell death protein 1 (PD-1) checkpoint blockade, was also evaluated for association with clinical response and survival. It was observed that neoadjuvant ipilimumab-HDI therapy demonstrated an upregulation of immune-related genes, chemokines, and transcription regulator genes involved in immune cell activation, function, or cell proliferation. Importantly, increased expression of baseline pro-inflammatory genes CCL19, CD3D, CD8A, CD22, LY9, IL12RB1, C1S, C7, AMICA1, TIAM1, TIGIT, THY1 was associated with longer OS (p < 0.05). In addition, multiple genes that encode a component or a regulator of the extracellular matrix such as MMP2 and COL1A2 were identified post-treatment as being associated with longer RFS and OS. In all baseline tissues, high TIS scores were associated with longer OS (p = 0.0166). Also, downregulated expression of cell proliferation-related genes such as CUL1, CCND1 and AAMP at baseline was associated with pathological and radiological response (unadjusted p < 0.01). In conclusion, we identified numerous genes that play roles in multiple biological pathways involved in immune activation, immune suppression and cell proliferation correlating with pathological/radiological responses following neoadjuvant immunotherapy highlighting the complexity of immune responses modulated by immunotherapy. Our observations suggest that TIS may be a useful biomarker for predicting survival outcomes with combination immunotherapy.