Digital Spatial Profiling Research Articles

CD163+ macrophages in mantle cell lymphoma induce activation of prosurvival pathways and immune suppression

AbstractMantle cell lymphoma (MCL) is dependent on a supportive tumor immune microenvironment (TIME) in which infiltration of CD163+ macrophages has a negative prognostic impact. This study explores how abundance and spatial localization of CD163+ cells are associated with the biology of MCL, using spatial multiomic investigations of tumor and infiltrating CD163+ and CD3+ cells. A total of 63 proteins were measured using GeoMx digital spatial profiling in tissue microarrays from 100 diagnostic MCL tissues. Regions of interest were selected in tumor-rich and tumor-sparse tissue regions. Molecular profiling of CD163+ macrophages, CD20+ MCL cells, and CD3+ T-cells was performed. To validate protein profiles, 1811 messenger RNAs were measured in CD20+ cells and 2 subsets of T cells. Image analysis was used to extract the phenotype and position of each targeted cell, thereby allowing the exploration of cell frequencies and cellular neighborhoods. Proteomic investigations revealed that CD163+ cells modulate their immune profile depending on their localization and that the immune inhibitory molecules, V-domain immunoglobulin suppressor of T-cell activation and B7 homolog 3, have higher expression in tumor-sparse than in tumor-rich tissue regions and that targeting should be explored. We showed that MCL tissues with more abundant infiltration of CD163+ cells have a higher proteomic and transcriptional expression of key components of the MAPK pathway. Thus, the MAPK pathway may be a feasible therapeutic target in patients with MCL with CD163+ cell infiltration. We further showed the independent and combined prognostic values of CD11c and CD163 beyond established risk factors.

A hot and cold tumor‑related prognostic signature for stage II colorectal cancer.

Globally, colorectal cancer (CRC) is one of the most lethal and prevalent malignancies. Based on the presence of immune cell infiltration in the tumor microenvironment, CRC can be divided into immunologically 'hot' or 'cold' tumors, which in turn leads to the differential efficacy of immunotherapy. However, the immune characteristics of hot and cold CRC tumors remain largely elusive, prompting further investigation of their properties regarding the tumor microenvironment. In the present study, a predictive model was developed based on the differential expression of proteins between cold and hot CRC tumors. First, the differentially expressed proteins (DEPs) were identified using digital spatial profiling and mass spectrometry-based proteomics analysis, and the pathway features of the DEPs were analyzed using functional enrichment analysis. A novel eight-gene signature prognostic risk model was developed (IDO1, MAT1A, NPEPL1, NT5C, PTGR2, RPL29, TMEM126A and TUBB4B), which was validated using data obtained from The Cancer Genome Atlas. The results revealed that the risk score of the eight-gene signature acted as an independent prognostic indicator in patients with stage II CRC (T3-4N0M0). It was also found that a high-risk score in the eight-gene signature was associated with high immune cell infiltration in patients with CRC. Taken together, these findings revealed some of the differential immune characteristics of hot and cold CRC tumors, and an eight-gene signature prognostic risk model was developed, which may serve as an independent prognostic indicator for patients with stage II CRC (T3-4N0M0).

Molecular analysis of human tick-bitten skin yields signatures associated with distinct spatial and temporal trajectories - A proof-of-concept study

Tick-associated diseases present challenges due to tridirectional interactions among host-specific responses, tick toxins and salivary proteins as well as microbes. We aimed to uncover molecular mechanisms in tick-bitten skin samples (cases) and contralateral skin samples (controls) collected simultaneously from the same participants, using spatial transcriptomics. Cases and controls analysed using NanoString GeoMx Digital Spatial Profiler identified 274 upregulated and 840 downregulated differentially expressed genes (DEGs), revealing perturbations in keratinization and immune system regulation. Samples of skin biopsies taken within 72 h post tick-bite DEGs had changes in protein metabolism and viral infection pathways as compared to samples taken 3 months post tick-bite, which instead displayed significant perturbations in several epigenetic regulatory pathways, highlighting the temporal nature of the host response following tick bites. Within-individual signatures distinguished tick-bitten samples from controls and identified between-individual signatures, offering promise for future biomarker discovery to guide prognosis and therapy.

Long-term follow-up of levonorgestrel intrauterine device for atypical hyperplasia and early endometrial cancer reveals relapse characterized by immune exhaustion.

Nonsurgical treatment options are increasingly needed for endometrial atypical hyperplasia (AH) and endometrioid endometrial cancer (EEC). Despite promising initial response rates, prospective long-term data and determinants for relapse are limited. Follow-up data from patients in our prospective phase II trial of LIUD for AH/G1EEC were collected from medical records. Spatial transcriptomics (Nanostring GeoMX digital spatial profiling) with in silico cell type deconvolution and pathway analyses were employed on longitudinal biopsy samples from five patients across pre-treatment, on-treatment, and relapse. Of 43 participants exhibiting initial response to LIUD, 41 had follow-up data. Sixteen (39%) experienced relapse. Clinical factors associated with shorter response duration included younger age, initial diagnosis of G1EEC, lack of response at six months, premenopausal status, and Hispanic ethnicity (p<0.05), but only six-month response status remained a significant predictor in a multivariate model (p=0.023). LIUD increased abundance of NK cells (DMCP-counter score=46.13, FDR=0.004) and cytotoxic lymphocytes (DMCP-counter score=277.67, FDR=0.004), as well as lymphocyte cytotoxicity markers PRF1 (log2FC=1.62, FDR=0.025) and GZMA (log2FC=2.47, FDR=0.008). NK cells were reduced at relapse (DMCP-counter score=-55.96, FDR=0.02). Immune-related pathways (IFNα-response and TGFβ-signaling) were enriched at relapse (FDR<0.05). IDO1 expression, reflecting immune exhaustion, was upregulated at relapse (FDR<0.05). Upfront resistance and relapse after initial response to LIUD for AH/G1EEC impacts nearly half of patients, remaining a major hurdle for non-surgical treatment of AH/G1EEC. Molecular studies evaluating longitudinal biopsies from a small cohort implicate immune mechanisms at relapse, including reversal of progestin-related immunomodulation and increased immune exhaustion.

Mapping the Spatial Dynamics of the CD4+ T Cell Spectrum in Classical Hodgkin Lymphoma

As around 25% to 30% of classical Hodgkin lymphoma (cHL) patients with advanced stages do not respond to standard therapies, the tumor microenvironment of cHL is one avenue that may be explored with the aim of improving risk stratification. CD4+ T cells are thought to be one of the main cell types in the tumor microenvironment. However, few immune signatures have been studied, and many of these lack related spatial data. Thus, our aim is to spatially resolve the CD4+ T cell subtypes that influence cHL outcome, depicting new immune signatures or transcriptional patterns that are in crosstalk with the tumor cells. This study was conducted using the NanoString GeoMx digital spatial profiling technology, based on the selection of distinct functional areas of patients’ tissues followed by gene-expression profiling. The goals were to assess the differences in CD4+ T cell populations between tumor-rich and immune-predominant areas defined by different CD30 and PD-L1 expression levels and seek correlations with clinical metadata. Our results depict a complex map of CD4+ T cells with different functions and differentiation states that are enriched at distinct locations, the flux of cytokines and chemokines that could be related to these, and the specific relationships with the clinical outcome.

Spatial Transcriptome-Wide Profiling of Small Cell Lung Cancer Reveals Intra-Tumoral Molecular and Subtype Heterogeneity.

Small cell lung cancer (SCLC) is a highly aggressive malignancy characterized by rapid growth and early metastasis and is susceptible to treatment resistance and recurrence. Understanding the intra-tumoral spatial heterogeneity in SCLC is crucial for improving patient outcomes and clinically relevant subtyping. In this study, a spatial whole transcriptome-wide analysis of 25 SCLC patients at sub-histological resolution using GeoMx Digital Spatial Profiling technology is performed. This analysis deciphered intra-tumoral multi-regional heterogeneity, characterized by distinct molecular profiles, biological functions, immune features, and molecular subtypes within spatially localized histological regions. Connections between different transcript-defined intra-tumoral phenotypes and their impact on patient survival and therapeutic response are also established. Finally, a gene signature, termed ITHtyper, based on the prevalence of intra-tumoral heterogeneity levels, which enables patient risk stratification from bulk RNA-seq profiles is identified. The prognostic value of ITHtyper is rigorously validated in independent multicenter patient cohorts. This study introduces a preliminary tumor-centric, regionally targeted spatial transcriptome resource that sheds light on previously unexplored intra-tumoral spatial heterogeneity in SCLC. These findings hold promise to improve tumor reclassification and facilitate the development of personalized treatments for SCLC patients.

Digital spatial profiling of the microenvironment of muscle invasive bladder cancer

Muscle invasive bladder cancer (MIBC) is a molecularly diverse disease with varied clinical outcomes. Molecular studies typically employ bulk sequencing analysis, giving a transcriptomic snapshot of a section of the tumour. However, tumour tissues are not homogeneous, but are composed of distinct compartments such as the tumour and stroma. To investigate the molecular profiles of bladder cancer, whilst also maintaining the spatial complexity of the tumours, we employed whole transcriptome Digital Spatial Profiling (DSP). With this method we generated a dataset of transcriptomic profiles of tumour epithelium, stroma, and immune infiltrate. With these data we investigate the spatial relationship of molecular subtype signatures and ligand signalling events. We find that Basal/Squamous and Classical subtypes are mostly restricted to tumour regions, while the stroma-rich subtype signatures are abundant within the stroma itself. Additionally, we identify ligand signalling events occurring between tumour, stroma, and immune infiltrate regions, such as immune infiltrate derived GPNMB, which was highly correlated with VEGFA expression within the tumour. These findings give us new insights into the diversity of MIBC at a molecular level and provide a dataset with detailed spatial information that was not available before in bladder cancer research.

Hepatic LRP-1 plays an important role in amyloidosis in Alzheimer's disease mice: Potential role in chronic heavy alcohol feeding

BackgroundHepatic lipoprotein receptor-related protein 1 (LRP-1) plays a central role in peripheral amyloid beta (Aβ) clearance, but its importance in Alzheimer's disease (AD) pathology is understudied. Our previous work showed that intragastric alcohol feeding to C57BL/6 J mice reduced hepatic LRP-1 expression which correlated with significant AD-relevant brain changes. Herein, we examined the role of hepatic LRP-1 in AD pathogenesis in APP/PS1 AD mice using two approaches to modulate hepatic LRP-1, intragastric alcohol feeding to model chronic heavy drinking shown by us to reduce hepatic LRP-1, and hepato-specific LRP-1 silencing. MethodsEight-month-old male APP/PS1 mice were fed ethanol or control diet intragastrically for 5 weeks (n = 7–11/group). Brain and liver Aβ were assessed using immunoassays. Three important mechanisms of brain amyloidosis were investigated: hepatic LRP-1 (major peripheral Aβ regulator), blood-brain barrier (BBB) function (vascular Aβ regulator), and microglia (major brain Aβ regulator) using immunoassays. Spatial LRP-1 gene expression in the periportal versus pericentral hepatic regions was confirmed using NanoString GeoMx Digital Spatial Profiler. Further, hepatic LRP-1 was silenced by injecting LRP-1 microRNA delivered by the adeno-associated virus 8 (AAV8) and the hepato-specific thyroxine-binding globulin (TBG) promoter to 4-month-old male APP/PS1 mice (n = 6). Control male APP/PS1 mice received control AAV8 (n = 6). Spatial memory and locomotion were assessed 12 weeks after LRP-1 silencing using Y-maze and open-field test, respectively, and brain and liver Aβ were measured. ResultsAlcohol feeding reduced plaque-associated microglia in APP/PS1 mice brains and increased aggregated Aβ (p < 0.05) by ELISA and 6E10-positive Aβ load by immunostaining (p < 0.05). Increased brain Aβ corresponded with a significant downregulation of hepatic LRP-1 (p < 0.01) at the protein and transcript level, primarily in pericentral hepatocytes (zone 3) where alcohol-induced injury occurs. Hepato-specific LRP-1 silencing significantly increased brain Aβ and locomotion hyperactivity (p < 0.05) in APP/PS1 mice. ConclusionChronic heavy alcohol intake reduced hepatic LRP-1 expression and increased brain Aβ. The hepato-specific LRP-1 silencing similarly increased brain Aβ which was associated with behavioral deficits in APP/PS1 mice. Collectively, our results suggest that hepatic LRP-1 is a key regulator of brain amyloidosis in alcohol-dependent AD.

Spatial multi-omics of human skin reveals KRAS and inflammatory responses to spaceflight

Spaceflight can change metabolic, immunological, and biological homeostasis and cause skin rashes and irritation, yet the molecular basis remains unclear. To investigate the impact of short-duration spaceflight on the skin, we conducted skin biopsies on the Inspiration4 crew members before (L-44) and after (R + 1) flight. Leveraging multi-omics assays including GeoMx™ Digital Spatial Profiler, single-cell RNA/ATAC-seq, and metagenomics/metatranscriptomics, we assessed spatial gene expressions and associated microbial and immune changes across 95 skin regions in four compartments: outer epidermis, inner epidermis, outer dermis, and vasculature. Post-flight samples showed significant up-regulation of genes related to inflammation and KRAS signaling across all skin regions. These spaceflight-associated changes mapped to specific cellular responses, including altered interferon responses, DNA damage, epithelial barrier disruptions, T-cell migration, and hindered regeneration were located primarily in outer tissue compartments. We also linked epithelial disruption to microbial shifts in skin swab and immune cell activity to PBMC single-cell data from the same crew and timepoints. Our findings present the inaugural collection and examination of astronaut skin, offering insights for future space missions and response countermeasures.

Clinicopathological features and cancer transcriptomic profiling of poorly cohesive gastric carcinoma subtypes.

Gastric poorly cohesive carcinoma (PCC) manifests with a diffuse pattern and diverse tumor cell morphologies, often indicating a more unfavorable prognosis. Recent consensus has reclassified PCC based on the proportion of signet-ring cells (SRCs) in tumors for research purposes. The two most distinct subtypes, poorly cohesive carcinoma not otherwise specified (PCC-NOS) and signet-ring cell carcinoma (SRCC), are characterized by less than 10% and more than 90% SRCs, respectively. However, research comparing the clinicopathological and transcriptomic differences between these subtypes remains limited. In this study, we conducted a comparative analysis of clinicopathological features in 55 advanced-stage PCCs, consisting of 43 PCC-NOS and 12 SRCC cases. Subsequently, 12 PCC-NOS and 5 SRCC cases were randomly selected for initial cancer-related gene expression profiling and pathway enrichment analysis using the GeoMx digital spatial profiler, followed by validation in a separate validation group comprising 16 PCC-NOS and 6 SRCC cases. These transcriptomic findings were then correlated with tumor morphology and clinicopathological data. PCC-NOS cases exhibited larger tumor size, a higher prevalence of pathological N3 disease, and a worse 1-year progression-free survival rate compared to SRCC cases. Clustering of PCC-NOS and SRCC was successfully achieved using the GeoMx Cancer Transcriptome Atlas. Among all studied genes, only MMP7 showed differential expression, with its overexpression significantly associated with the PCC-NOS subtype, increased perineural invasion, and earlier disease progression. Pathway analysis revealed significantly enriched pathways in PCC-NOS related to vesicle-mediated transport, adaptive immune systems, oncogenic signaling, and extracellular matrix organization, while SRCC displayed significant enrichment in pathways associated with respiratory electron transport and the cell cycle. In conclusion, this study compares and correlates clinicopathological features and transcriptomic data between PCC-NOS and SRCC at advanced stages, employing the latest consensus classification and a novel platform for analysis.

Identification of Early Events in Serrated Pathway Colorectal Tumorigenesis by Using Digital Spatial Profiling.

The colorectal serrated pathway involves precursor lesions known as sessile serrated lesions (SSL) and traditional serrated adenomas (TSA). Mutations in BRAF or KRAS are crucial early events in this pathway. Additional genetic and epigenetic changes contribute to the progression of these lesions into high-grade lesions and, eventually, invasive carcinoma. We employed digital spatial profiling to investigate the transcriptional changes associated with SSL and TSA. The genes identified are confirmed by immunohistochemical (IHC) staining. Colorectal cancer (CRC) cell lines with CEACAM6 overexpression and knockdown were established to study the roles of CEACAM6 on tumorigenesis of CRC. Ten genes were upregulated in SSL and TSA, and seven were upregulated in both types of lesions. IHC staining confirmed overexpression of CEACAM6, LCN2, KRT19, and lysozyme in SSL and TSA. CEACAM6 expression is an early event in the serrated pathway but a late event in the conventional pathway. Using cell line models, we confirmed that CEACAM6 promotes CRC cells' proliferation, migration, and invasion abilities. These results highlight that the transcriptional changes in the early stages of tumorigenesis exhibit relative uniformity. Identifying these early events may hold significant promise in elucidating the mechanisms behind tumor initiation.

Multiomic analysis in ovarian clear cell carcinoma to gain insight into biology and therapeutic opportunities.

e17540 Background: Ovarian clear cell carcinoma (OCCC) is a rare subtype with distinct histologic and molecular features. The prognosis of advanced/recurrent OCCC treated with standard of care surgery and platinum-based chemotherapy remains poor due to inherent chemoresistance. Therefore, there is an urgent need to identify novel biological insights and personalized treatment approaches to improve outcomes. Methods: Participants were referred to the SMMART Program at the OHSU Knight Cancer Institute. Archival or fresh tissue was obtained from metastatic sites and processed through a multiomics platform that included immunohistochemistry (Ki67, HER2, CD4, CD8, PD-L1), fluorescent in-situ hybridization, targeted next-generation sequencing panel (225 genes), whole transcriptomic sequencing, chromosomal microarray, and multiplexed analysis of 67 key cancer proteins and phosphoproteins on the Nanostring GeoMx Digital Spatial Profiler. Results: 5 participants were enrolled with an average age of 62 years (range 52-72). Therapeutically relevant findings included detectable HER2 expression by IHC (2+, N=2; 3+, N=2), HER2 amplification by FISH (N=2), and high HER2 and phospho-HER2 protein expression compared to a breast cancer cohort (phospho-HER2 supports active HER2 signaling in the tumor cells). Additional findings include TMB &lt;10 muts/Mb (N=5), microsatellite stable (N=5), Ki67 40-60% (N=4), mild to moderate CD4 and/or CD8 infiltration (N=4), PD-L1 &gt;1% (N=3), low homologous recombination deficiency score (N=4), chromosomal losses including ATM, SMAD4, RB1, and TP53, and single nucleotide variations in ARID1A (N=4). Two participants with HER2 expression were treated with off-label trastuzumab deruxtecan for 6 or 7 cycles yielding significant clinical benefit. Participant A (HER2 = 3+, FISH negative) had an exceptional response while Participant B (HER2 = 2+, FISH amplified) had stable disease. The toxicity profile was consistent with literature. Conclusions: Our OCCC cohort revealed elevated HER2 expression with or without genomic amplification in 4/5 cases. Two participants experienced clinical benefit when treated with trastuzumab deruxtecan, supporting further exploration of HER2 antibody drug conjugates in this aggressive disease. More study is needed to determine the frequency of HER2 overexpression and amplification in OCCC, the importance of HER2 on OCCC pathophysiology, and the therapeutic benefit of targeting HER2 in this population.

Spatially resolved multi-region proteomics for prediction of clinical outcome in deficient mismatch repair metastatic colorectal cancer treated with PD-1 blockade.

3520 Background: In patients with metastatic colorectal cancer (CRC) and deficient DNA mismatch repair (dMMR), we previously found transcriptomic signatures reflecting stromal compartment composition and proliferative capacity that can predict the efficacy of immune checkpoint blockade (ICB). To identify proteomic markers of response and/or resistance to ICB, we performed spatially resolved proteomic analysis of the tumor microenvironment. Methods: Patients with surgically resected, primary dMMR metastatic CRCs received pembrolizumab (Mayo Clinic, N=30). In FFPE tissues, Digital Spatial Profiling of 71 proteins was performed using panels covering immune cells, immune activation, cell death, PI3/AKT signaling, MAPK signaling, and pan-tumor markers (GeoMx nCounter). In ten regions of interest/slide (~600 µm2) at the tumor invasive margin, protein expression was determined in four cellular compartments: epithelium (panCK+), immune cells (CD45+; leucocyte common antigen), stroma (panCK-), and nonimmune stroma (panCK-/CD45-). After normalization, protein abundance per compartment was analyzed by response (RECIST 1.1) [OmicVerse Python]. Differentially expressed proteins (binary) were examined by progression-free survival (PFS) using multivariable Cox regression with covariate adjustment (Lifelines Python). Results: Pre-treatment panels including 71 proteins were analyzed by treatment response and PFS after PD-1 inhibition. Proteins found to be significantly associated with response (p&lt; 0.01) and PFS are described. Across all compartments, responders showed upregulation of NK cells and T cell subsets (CD3+, CD4+, CD8+), and ICOS; strongest association was found for CD4+: HRadj: 0.16; 95%CI: 0.08-0.31; p&lt; 0.001. Nonresponders showed higher phospho-JNK (HRadj: 2.09; 95%CI: 1.28-3.42; p=0.003) and SMA. In the epithelial compartment, upregulation of the Treg marker CD25 and reduced caspase-9 were observed in responders. Among responders, the stromal compartment was highly enriched in HLA-DR, dendritic cells (CD11c), PD-L1, immune activation markers (CD40, CD44, CD80, CD27), CD95/Fas, and GZMB. In the stromal compartment, downregulation of BCL6 and phospho-GSK3B was observed in responders. The immune compartment of responders showed increased beta2-microglobulin and Ki-67, while the non-immune stromal compartment showed increase in the macrophage marker CD68. Importantly, proteins found to be upregulated in responders were significantly and independently associated with longer PFS whereas the opposite was true for nonresponders. Conclusions: In metastatic dMMR CRCs, spatially defined protein markers were significantly associated with anti-PD-1 response, and independently associated with patient PFS (strongest predictive utility for CD4 in immune stroma and phospho-JNK in epithelium).

Secukinumab and Dead Sea Climatotherapy Impact Resolved Psoriasis Skin Differently Potentially Affecting Disease Memory.

Secukinumab and Dead Sea treatment result in clear skin for many psoriasis patients, through distinct mechanisms. However, recurrence in the same areas after treatments suggests the existence of a molecular scar. We aimed to compare the molecular and genetic differences in psoriasis patients who achieved complete response from secukinumab and Dead Sea climatotherapy treatments. We performed quantitative immunohistochemical and transcriptomic analysis, in addition to digital spatial profiling of skin punch biopsies. Histologically, both treatments resulted in a normalization of the lesional skin to a level resembling nonlesional skin. Interestingly, the transcriptome was not normalized by either treatments. We revealed 479 differentially expressed genes between secukinumab and Dead Sea climatotherapy at the end of treatment, with a psoriasis panel identifying SERPINB4, SERPINB13, IL36G, IL36RN, and AKR1B10 as upregulated in Dead Sea climatotherapy compared with secukinumab. Using digital spatial profiling, pan-RAS was observed to be differentially expressed in the microenvironment surrounding CD103+ cells, and IDO1 was differentially expressed in the dermis when comparing the two treatments. The differences observed between secukinumab and Dead Sea climatotherapy suggest the presence of a molecular scar, which may stem from mechanistically different pathways and potentially contribute to disease recurrence. This may be important for determining treatment response duration and disease memory.

Age-related Differences in Spatially Resolved Transcriptomic Profiles of Patients With Hormone Receptor-positive Breast Carcinoma.

Patients' age may influence the response to chemotherapy and the clinical course of breast carcinoma. This study aimed to compare the spatial transcriptomic profiles between younger (≤50 years) and older (>50 years) patients with hormone receptor (HR)-positive breast carcinoma. Seven cases of breast carcinoma were included. We performed digital spatial profiling and bioinformatic analysis to investigate the spatial transcriptomes of the epithelial and stromal compartments. In the epithelial compartment of three young-age breast carcinoma (YABC) cases, we found 21 up-regulated and 7 down-regulated genes. The top two most up-regulated genes were serpin peptidase inhibitor clade A member 1 and serine protease. The gene ontology enrichment analysis revealed a significant up-regulation of genes defining ribosomal structures and functions in YABCs. The gene set enrichment analysis revealed that gene sets defining early and late responses to estrogen, response to interferon-α, and tumor necrosis factor-α signaling were significantly enriched in YABCs. We described for the first time the age-related differences in spatially resolved transcriptomic profiles and up-regulated transcriptional pathways of HR-positive breast carcinoma. Our observations highlight the critical need for age-specific treatment strategies for breast carcinoma management.

#163 Multiome-seq and spatial transcriptomics of refractory PLA2RAb(+) membranous nephropathy reveals importance of immunomodulatory function of B cells

Abstract Background and Aims PLA2RAb-associated membranous nephropathy (MN) is a glomerulonephritis that may require immunosuppression to achieve proteinuria remission. However, the disease course is heterogeneous, and some patients exhibit refractory or relapsing nephrotic-range proteinuria. Additional biomarker studies are necessary to dissect the molecular profiles that reflect substantial differences in treatment response and clinical courses. Multiome-seq includes the investigation of gene expression profiles and ATAC-sequencing for open-chromatin regions, which may reveal pathophysiologic gene expression biomarkers and pathways. Method We performed multiome-sequencing of peripheral blood mononuclear cells (PBMCs) from 6 refractory or relapsing primary PLA2RAb-associated MN cases, 7 PLA2RAb-associated MN cases that reached complete remission, and disease controls including 6 minimal change disease and 6 diabetic kidney disease cases. We also collected PBMCs from 7 healthy controls. Cell type clustering, comparison of cell type proportions, and differential gene expression analysis were performed. Additionally, we conducted ATAC-sequencing analysis and identified open-chromatin profiles specific to the refractory or relapsing MN group. Lastly, we investigated the glomerular spatial transcriptome using GeoMx digital spatial profiling to determine whether transcriptomic profiles differ based on the disease courses of MN in the initial kidney tissue. Results We successfully identified diverse PBMCs in the multiome-seq data. In the DEG analysis, we identified a number of cell-type-specific DEGs in the refractory or relapsing MN cases. Of these, CD83 was the top gene that was lowly expressed in the B cells of refractory/relapsing MN compared to the remission cases. In the ATAC data, among the number of differentially expressed peaks, NFKB2 was the region that consistently showed a low open chromatin profile in the B cells of CD83. When comparing the cell subtype proportions, refractory or relapsing MN cases showed low regulatory T cells in the peripheral blood. Lastly, in the kidney tissue spatial transcriptomic profiling, there were no significant DEGs between the refractory/relapsing MN and MN with remission in the glomerulus. Conclusion The clinical course or risk of remission of PLA2RAb-associated MN was more likely a result of differences in systemic immunity rather than differences in the glomerular transcriptome at the initial diagnosis. Both CD83 and NFκB2 are well-known biomarkers of activated B cells that may have immunomodulatory roles leading to the enrichment of regulatory T cells. Such a deficient immunomodulatory status in refractory/relapsing PLA2RAb-associated MN may explain the disease course.

#1169 Class-specific gene expression profiling of glomerulus of lupus nephritis by spatial transcriptomics

Abstract Background and Aims Lupus nephritis exhibit various clinical features and subclassifications. Subclassification-specific alteration of transcriptomic profile of lupus nephritis needs additional study which may provide new insights regarding the subclass-specific pathophysiologic mechanism. Method Spatial transcriptomic profiling was conducted using GeoMx Digital Spatial Profiler on formalin-fixed paraffin-embedded kidney biopsy specimens obtained from control (n = 7), LN class I to V (n = 24). We also profiled native kidney primary glomerulonephritis (e.g. anti-PLA2R-associated MN) as additional disease control group. Subclass-specific transcriptomic alteration of glomerulus was compared between the subclasses and primary glomerulonephritis cases. Gene Ontology (GO) terms for DEGs were annotated using the ToppGene suite. Results We first focused on the glomerular transcriptomic alteration in class IV lupus nephritis, which was the most active form with high degree of proteinuria. We found 104 upregulated DEGs and 166 downregulated DEGs which were specifically different in the class IV lupus nephritis cases but not in other subclasses of the disease. The upregulated GO terms included interleukin 10 receptor activity, T cell activation, leukocyte activation, and vesicle membrane. On the other hand, the downregulated DEGs were mostly annotated with GO terms of double-stranded DNA binding domain. In our analysis regarding class V lupus nephritis, total of 128 DEGs were downregulated in class V LN compared to MN. These genes were also downregulated in class V LN when compared to control samples, while no significant expression differences were observed in other LN classes compared to controls. Eleven of 128 downregulated DEGs were annotated with biological process GO terms related to Golgi vesicle transport. Regarding cellular component GO terms, Golgi apparatus (30/128 DEGs), endoplasmic reticulum (ER) membrane (19/128 DEGs), and ER-Golgi intermediate compartment (6/128 DEGs) were among the annotated GO terms associated with the DEGs. DEGs in annotation included GOLGA3, TMED7, BET1L, and RAB10. Conclusion We profiled the class-specific transcriptomic changes in glomerulus of lupus nephritis. Class IV LN was associated upregulation of interleukin 10-mediated inflammatory pathway and downregulation of double stranded-DNA binding domain. The ER-Golgi network and related molecular pathway may be involved in the distinct pathogenesis and clinical features of class V LN. The findings provide novel view of the class-specific pathophysiologic mechanism of lupus nephritis.

HLA-A+ tertiary lymphoid structures with reactivated tumor infiltrating lymphocytes are associated with a positive immunotherapy response in esophageal squamous cell carcinoma.

Immune checkpoint blockade (ICB) therapy provides remarkable clinical benefits for multiple cancer types. However, the overall response rate to ICB therapy remains low in esophageal squamous cell carcinoma (ESCC). This study aimed to identify biomarkers of ICB therapy for ESCC and interrogate its potential clinical relevance. We investigated gene expression in 42 treatment-naïve ESCC tumor tissues and identified differentially expressed genes, tumor-infiltrating lymphocytes and immune-related genes signatures associated with differential immunotherapy responses. We systematically assessed the tumor microenvironment using the NanoString GeoMx digital spatial profiler, single-cell RNA-seq and multiplex immunohistochemistry in ESCC. Finally, we evaluated the associations between HLA-A-positive tertiary lymphoid structures (TLSs) and patients' responses to ICB in 60 ESCC patients. Tumor infiltrating B lymphocytes and several immune-related gene signatures, such as the antigen presenting machinery (APM) signature, are significantly elevated in ICB treatment responders. Multiplex immunohistochemistry identified the presence of HLA-A+ TLSs and showed that TLS-resident cells increasingly express HLA-A as TLSs mature. Most TLS-resident HLA-A+ cells are tumor-infiltrating T (TIL-T) or tumor-infiltrating B (TIL-B) lymphocytes. Digital spatial profiling of spatially distinct TIL-T lymphocytes and single-cell RNA-seq data from 60 ESCC tumor tissues revealed that CXCL13-expressing exhausted TIL-Ts inside TLSs are reactivated with elevated expression of the APM signature as TLSs mature. Finally, we demonstrated that HLA-A+ TLSs and their major cellular components, TIL-Ts and TIL-Bs, are associated with a clinical benefit from ICB treatment for ESCC. HLA-A+ TLSs are present in ESCC tumor tissues. TLS-resident TIL-Ts with elevated expression of the APM signature may be reactivated. HLA-A+ TLSs and their major cellular components, TIL-Ts and TIL-Bs, may serve as biomarkers for ICB-treated ESCC patients.

Topographic analysis of pancreatic cancer by TMA and digital spatial profiling reveals biological complexity with potential therapeutic implications

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal human malignancies. Tissue microarrays (TMA) are an established method of high throughput biomarker interrogation in tissues but may not capture histological features of cancer with potential biological relevance. Topographic TMAs (T-TMAs) representing pathophysiological hallmarks of cancer were constructed from representative, retrospective PDAC diagnostic material, including 72 individual core tissue samples. The T-TMA was interrogated with tissue hybridization-based experiments to confirm the accuracy of the topographic sampling, expression of pro-tumourigenic and immune mediators of cancer, totalling more than 750 individual biomarker analyses. A custom designed Next Generation Sequencing (NGS) panel and a spatial distribution-specific transcriptomic evaluation were also employed. The morphological choice of the pathophysiological hallmarks of cancer was confirmed by protein-specific expression. Quantitative analysis identified topography-specific patterns of expression in the IDO/TGF-β axis; with a heterogeneous relationship of inflammation and desmoplasia across hallmark areas and a general but variable protein and gene expression of c-MET. NGS results highlighted underlying genetic heterogeneity within samples, which may have a confounding influence on the expression of a particular biomarker. T-TMAs, integrated with quantitative biomarker digital scoring, are useful tools to identify hallmark specific expression of biomarkers in pancreatic cancer.

Abstract PR004: Spatial proteomics and transcriptomics reveal an altered immune cell landscape in bladder cancer patients unresponsive to BCG treatment

Abstract Introduction: Patients with high-risk non-muscle invasive bladder cancer (NMIBC) are recommended treatment with Bacillus Calmette-Guérin (BCG). The therapeutic effect of BCG is highly dependent on the host immune system and the tumor microenvironment (TME). The cellular composition and the functional status of the TME have been shown to play crucial roles for treatment efficacy. However, much is still unknown regarding the cellular interactions and mechanisms of BCG. Using spatial proteomics and transcriptomics, we investigated the immune cell landscape in a cohort of BCG treated patients. Materials and methods: A total of 183 tumors from 111 patients with NMIBC treated with BCG were included on a tissue microarray (TMA). The TMA included paired pre- and post-BCG tumors. We analyzed the tumor tissue using Imaging Mass Cytometry (IMC; StandardBioTools) for simultaneous detection of 35 immune and tumor-related proteins at single cell resolution. Ilastik was used to generate probability maps and DeepCell Mesmer was used for cell segmentation. Markers were quantified for the mean intensity per cell. GeoMX Digital Spatial Profiling (DSP) was used for spatial proteomics and transcriptomics profiling of tumor and stromal areas, targeting 63 protein targets and whole transcriptome amplification, respectively. Results: We identified a total of 364,504 cells in the tumor tissue using IMC. After computing a neighborhood graph using the log normalised and z-scaled signal and embedding it using UMAP, leiden clustering was performed to detect communities of unique cell types by interrogating the presence of phenotypic markers. We identified major lineages such as macrophages, CD8 T cells, dendritic cells (DCs), NK cells and neutrophils in the TME of analyzed tumors. Specifically, we observed more CD8 T cells after BCG (p=0.015), especially in females (p=0.048). Females also had higher levels of DCs and macrophages after treatment compared to males (p=0.024 and p=0.048). Patients who progressed to MIBC after BCG had higher numbers of CD8 T cells and NK cells prior to BCG (p=0.039 and p=0.021). Prolonged exposure to interferons has been associated with immune evasion and tumor cell proliferation. In concordance, we found high pretreatment protein expression of IFNα and IFNγ in tumor regions from patients with a pronounced signature of CD8 T cell exhaustion after BCG treatment using GeoMX DSP proteomic analysis. Data analysis is currently ongoing, and additional results will be presented at the conference. Conclusion: The composition and functional status of the TME is associated with clinical and biological features such as immune cell abundance, CD8 T cell exhaustion, sex and progression. A greater understanding of the TME may help identify patients unresponsive to BCG earlier and improve the understanding of biological differences in tumor development and aggressiveness. This may ultimately improve patient outcomes. Citation Format: Trine Strandgaard, Tine Ginnerup Andreasen, Tessa Jane Divita, Liina Salminen, Sean Houghton, Kasper Thorsen, Iver Nordentoft, Iyinyeoluwa Okulate, Alexander Schmitz, Søren Riis Paludan, John Sfakianos, Amir Horowitz, Jørgen Bjerggaard Jensen, Lars Dyrskjøt. Spatial proteomics and transcriptomics reveal an altered immune cell landscape in bladder cancer patients unresponsive to BCG treatment [abstract]. In: Proceedings of the AACR Special Conference on Bladder Cancer: Transforming the Field; 2024 May 17-20; Charlotte, NC. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(10_Suppl):Abstract nr PR004.