Hwang, Tae Hyun

Abstract Background: Recent studies support that high levels of tumor-infiltrating lymphocytes (TILs) were associated with better prognosis in HPV(+) oropharyngeal squamous cell carcinoma (OPSCC). While these studies provide a potential role for TILs as a prognostic biomarker, the analyses were relied on manual quantification performed by pathologists, resulting in inter-observer variability. Deep Learning (DL)-based whole slide hematoxylin and eosin (H&E) image analyses may overcome these challenges. However, this analysis only provides the presence of cell types without consideration for the functional roles of each cell within the tumor immune microenvironment (TME). In this work, we develop a computational pipeline to integrate paired H&E and immunohistochemistry (IHC) images to functionally characterize TILs and investigate their prognostic utility. Methods: We analyzed 88 patients: 42 stage I, 39 stage II, and 7 stage III. Our data contain both H&E and IHCs examining FoxP3, CD3, PD-L1, CD20, etc. on serial sections of the tissue. In-house DL-based H&E analysis used to identify TILs, tumors, and stroma in each tumor, then performed registration between adjacent H&E and IHC images from the same tissue. The patients were then classified into three basic immune phenotypes: immune inflamed (IN; high TILs in the tumor region), immune excluded (EX; TILs are mostly localized in stroma), and immune desert (ID; few/no TILs) based on TIL enrichment in the TME. To functionally characterize TILs, we quantified protein expression from the adjacent IHCs. For example, we further classified each patient into different subtypes based on enriched protein expression (e.g., FoxP3 high IN, FoxP3 low IN). We used the Kaplan-Meier method and multivariate Cox analysis to evaluate the prognostic value of different subtypes enriched with different proteins to predict disease-free survival (DFS). Results: The IN group with 43 patients was significantly associated with good prognosis. Interestingly, further stratification of the IN subgroup based on Foxp3 quantification on TIL regions (i.e., high FoxP3 IN and low FoxP3 IN) showed that high protein expression of FoxP3 in TILs in the IN subgroup is significantly associated with a better prognosis compared to other immune subgroups (HR, 0.16; p-value, 0.003). Multivariate analysis, including other clinical covariates showed that the immune subtypes associated with high FoxP3 are independently associated with DFS. These results demonstrate that DL-based integrative IHC and H&E image analysis could be used to identify subgroups with distinct clinical outcomes. Furthermore, our results reveal unknown roles for Foxp3 expression in the TILs in HPV(+) OPSCC as a prognostic biomarker, a finding which should be evaluated in a larger cohort. Citation Format: Sumanth Reddy Nakkireddy, Inyeop Jang, Minji Kim, Linda X. Yin, Michael Rivera, Joaquin J. Garcia MD, Kathleen R. Bartemes, David M. Routman, Eric J. Moore, Daniel J. Ma, Kathryn M. Van Abel, Tae Hyun Hwang. Integrative spatial analysis of paired IHC and H&E images identifies Foxp3 enriched tumor-infiltrating lymphocytes associated with disease-free survival in human papillomavirus (HPV)-related oropharyngeal squamous cell carcinoma. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4638.

464 Background: A larger proportion of Hispanic gastric cancer patients are diagnosed at a young age than patients of other ethnicities. Despite diagnosis at an early age, young Hispanic patients also present more frequently with advanced disease. Whether there are molecular differences between tumors from young and old Hispanic gastric cancer patients is unknown. Methods: We performed RNA-sequencing on pre-treatment tumor samples of gastric cancer patients of self-identified Hispanic ethnicity. We stratified patients into two groups: <50 years and ≥50 years. Differentially expressed genes between groups were detected using the DESeq2 package (RStudio). We performed Gene Set Enrichment Analysis (GSEA) to identify pathways that were uniquely overexpressed in each age group. Results: 57 patients were included in our study. 23 patients were <50 years and 33 patients were ≥50 years. Comparison of the top 10% most differentially expressed genes showed no clustering based on age. However, visual inspection of the top 40 most differentially expressed genes revealed 2 distinct patterns of gene expression between age groups. GSEA revealed upregulation of DNA replication, cell cycle, and DNA repair pathways in the younger patients The tumors in younger patients also had upregulation of multiple immune activation pathways. Cancers from older patients had upregulation of programmed cell death pathways. Conclusions: Early-onset gastric cancer in Hispanic patients may represent a distinct biologic process.

Abstract Neuropathological assessment at autopsy is the gold standard for diagnosing neurodegenerative disorders. We aimed to develop a pipeline for diagnosing Alzheimer's disease and other tauopathies, including corticobasal degeneration, globular glial tauopathy, Pick’s disease, and progressive supranuclear palsy. We used deep learning (DL)-based approach called clustering-constrained-attention multiple instance learning (CLAM) on whole slide images (WSIs) of tau immunohistochemistry in three brain regions from 120 patients. We also augmented gradient-weighted class activation mapping (Grad-CAM) to the model for visualizing cellular-level evidence in the model’s decisions. The model using the sections of cingulate and superior frontal gyri achieved the highest area under the curve (0.970±0.037) and diagnostic accuracy (0.873±0.087). Grad-CAM showed the highest attention in known pathognomonic tau lesions for each disease (e.g., Pick bodies for Pick’s disease). Our findings supported the feasibility of the DL-based approach for the classification task on WSIs, which encouraged further investigation, especially focusing on clinicopathological correlation studies.

Abstract Neuropathological assessment at autopsy is the gold standard for diagnosing neurodegenerative disorders. We aimed to develop a pipeline for diagnosing Alzheimer's disease and other tauopathies, including corticobasal degeneration, globular glial tauopathy, Pick’s disease, and progressive supranuclear palsy. We used deep learning (DL)-based approach called clustering-constrained-attention multiple instance learning (CLAM) on whole slide images (WSIs) of tau immunohistochemistry in three brain regions from 120 patients. We also augmented gradient-weighted class activation mapping (Grad-CAM) to the model for visualizing cellular-level evidence in the model’s decisions. The model using the sections of cingulate and superior frontal gyri achieved the highest area under the curve (0.970±0.037) and diagnostic accuracy (0.873±0.087). Grad-CAM showed the highest attention in known pathognomonic tau lesions for each disease (e.g., Pick bodies for Pick’s disease). Our findings supported the feasibility of the DL-based approach for the classification task on WSIs, which encouraged further investigation, especially focusing on clinicopathological correlation studies.

Abstract Purpose: We sought to identify biomarkers that predict overall survival (OS) and response to immune checkpoint inhibitors (ICI) for patients with gastric cancer. Experimental Design: This was a retrospective study of multiple independent cohorts of patients with gastric cancer. The association between tumor ACTA2 expression and OS and ICI response were determined in patients whose tumors were analyzed with bulk mRNA sequencing. Single-cell RNA sequencing (scRNA-seq) and digital spatial profiling data were used to compare tumors from patients with gastric cancer who did and did not respond to ICI. Results: Increasing tumor ACTA2 expression was independently associated with worse OS in a 567-patient discovery cohort [HR, 1.28 per unit increase; 95% confidence interval (CI), 1.02–1.62]. This finding was validated in three independent cohorts (n = 974; HR, 1.52 per unit increase; 95% CI, 1.34–1.73). Of the 108 patients treated with ICI, 56% of patients with low ACTA2 expression responded to ICI versus 25% of patients with high ACTA2 expression (P = 0.004). In an analysis of a publicly available scRNA-seq dataset of 5 microsatellite instability-high patients treated with ICI, the patient who responded to ICI had lower tumor stromal ACTA2 expression than the 4 nonresponders. Digital spatial profiling of tumor samples from 4 ICI responders and 5 ICI nonresponders revealed that responders may have lower ACTA2 expression in α-SMA–positive cancer-associated fibroblasts (CAF) than nonresponders (median: 5.00 vs. 5.50). Conclusions: ACTA2 expression is associated with survival and ICI response in patients with gastric cancer. ACTA2 expression in CAFs, but not in other cellular compartments, appears to be associated with ICI response.

Abstract Pediatric Acute Myeloid Leukemia (AML) with a FLT3 internal tandem duplication (FLT3-ITD) is a challenging disease due to poor outcomes in many patients. The 4-year progression-free survival is still only 31%. Current biomarkers are insufficient to predict why certain patients with FLT3-ITD AML relapse and others do not. The development of prognostic biomarkers in FLT3-ITD pediatric AML may help improve the outcomes and management of these patients. We acquired a panel of 37 diagnostic samples with 18 samples from FLT3-ITD pediatric AML patients with poor outcomes (relapse within 3 years) and 19 from favorable outcome patients that did not exhibit relapse for at least 5 years. Next, we performed FACS to isolate CD34+CD38dim and CD34+CD38+ cells. These two populations were then sequenced using 10x Genomics 3’ end single-cell RNA sequencing (scRNAseq). Using this scRNAseq dataset comprising over 250k single cells, we first investigated if the frequency of specific RNA clusters of AML cells may predict patient outcome. We found 5 clusters had significantly different frequencies between patients with favorable and poor outcomes (p &lt; 0.05). Using these clusters, we performed CIBERSORTx deconvolution on publicly available FLT3-ITD AML bulk RNAseq samples from TCGA (n = 36), BeatAML (n = 47), and the TARGET project (ntrain = 72; ntest = 26) to identify the inferred cluster composition present in each AML patient from these datasets. Next we built LASSO prediction models using the inferred scRNAseq clusters by CIBERSORTx and incorporating clinical features including leukemic blast percentage, bone marrow leukemic blast percentage, and WBC at diagnosis. We found these prediction models to be highly specific as a prognostic biomarker for pediatric AML (p &lt;&lt; 0.01). We were unable to train a model directly using genes from the publicly available bulk RNAseq data for pediatric FLT3-ITD AML, highlighting the relevance of using our deconvolution based biomarker. Between the high (&gt; median) and low (&lt; median) LASSO score groups, there was ~50% difference in the number of patients relapsing after 4 years. When we included this LASSO score in a multivariate cox regression (likelihood ratio p value = 0.007), the LASSO score (HR = 6.8, p = 0.026, CI = 1.26 - 37.28) showed additional predictive ability when including FLT3-ITD allelic ratio (HR = 1.83, p = 0.38, CI = 0.47 - 7.15) and MRD (HR = 3.95, p = 0.021 CI = 0.051 - 0.79). These results demonstrate that incorporating rare AML cell populations into a prognostic biomarker could add substantial clinical information to existing biomarkers. This indicates single-cell derived prognostic models can and should be explored further for pediatric FLT3-ITD AML and demonstrates a platform that could be used for other malignancies, including non-hematologic tumors, to create cell-type biomarkers. Citation Format: Robert Schauner, Zachary Jackson, Nethrie Idippily, Grace Lee, Sheela Karunanithi, Shivaprasaad Manjappa, Tae Hyun Hwang, David Wald. New platform for novel biomarker discovery using single-cell RNA sequencing on pediatric acute myeloid leukemia samples [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 4080.

Abstract Histological distribution and density of tumor infiltrating lymphocytes (TILs) are known to be informative to cancer patients’ outcome. However, due to the heterogeneity and huge size of digitized pathology slides, it has been a challenging problem to quantify TILs-related features in whole slide images (WSI). This study aimed to explore the impact of spatial distribution of TILs quantified by deep learning (DL) approaches based on digitalized WSI stained with hematoxylin and eosin in patients with colorectal cancer (CRC). The cancer and TILs regions are first separately identified by deep learning models in the WSI. The TILs densities at tumor invasive margins are then quantified by image analysis algorithms and used for prognosis of progression-free survival (PFS) in CRC patients. The prognostic impact of spatial distributions of TILs in CRC patients was explored in the Yonsei cohort (n=180) and validated in the TCGA cohort (n=268). Two experienced pathologists manually measured TILs at the most invasive margin as 0-3 by the Klintrup-Mäkinen (KM) grading method and compared to DL approaches. On multivariate analysis, TILs densities within 200µm of the invasive margin (f_im200) was remained as the most significant prognostic factor for PFS (HR 0.004 [95% CI, 0.0001-0.1502], p=.002) in the Yonsei cohort. On multivariate analysis using the TCGA dataset, f_im200 retained prognostic significance for PFS (HR 0.031, [95% CI 0.001-0.645], p=.024). Inter-observer agreement of manual KM grading was insignificant in both cohorts. The survival analysis by pathologists’ manual KM grading showed statistically significant different PFS from the TCGA cohort, but not the Yonsei cohort. In summary, automated quantification of TILs at the invasive margin showed a prognostic utility to predict PFS, and could provide reproducible TILs density measurement in CRC patients. Citation Format: Hongming Xu, Yoon Jin Cha, Jean Clemenceau, Jinhwan Choi, Sung Hak Lee, Jeonghyun Kang, Tae Hyun Hwang. Spatial analysis of tumor infiltrating lymphocytes in histological images using deep learning predicts progression-free survival in colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 642.

Abstract Background: Microsatellite Instability (MSI) is a prognostic and predictive biomarker which can guide treatments including immunotherapy for patients with gastrointestinal cancer. Universal MSI testing would benefit to patients to receive a better therapy, but many patients remain untested. To deliver broadly accessible MSI testing, deep learning models have demonstrated their feasibility for MSI prediction using H&E-stained whole-slide images (WSIs). However, to be deployed into clinical routine care, prediction models should be validated in datasets from multicenter and multiethnic groups. In addition, these models should provide uncertainty of prediction to help clinicians to make informed decisions. Method: We develop a MSI prediction model using WSIs based on Deep Gaussian process (DGP), a Bayesian model being able to model the uncertainty in prediction. We implement a DGP model in the transfer learning, where a WPI is decomposed into multiple non-overlap image patches which are converted to feature vectors by a pretrained convolutional neural network (CNN). Then, the DGP model makes a prediction for an image by averaging score function values at all the feature vectors in the image in weakly supervised learning (WSL). To test our method, we collected H&E stained colorectal/stomach cancer slides (n=1,619) from multi-institutions and multiethnic groups, including TCGA Colorectal (CRC; n=351) and Stomach Adenocarcinoma (STAD; n=174) and six datasets from tertiary hospitals in Korea (i.e., Yonsei-1 (n=174), Yonsei-1-remade (n=146), Yonsei-2 (n=95), St. Mary-1 (n=48), St. Mary-2 (n=50), and Yonsei-CLASSIC trial (n=581)). We also compare our method to multiple deep learning methods, including Densenet, Resnet, Shufflenet, Googlenet, Squeezenet and two WSL models, MIL (Campanell et al., 2019) and CLAM (Lu et al., 2021). Result: In each cancer type, we train each model using TCGA dataset as discovery cohort and tested on the rest of datasets as validation cohort. We use the area under of roc curve (AUC) for the evaluation metric. For the CRC cohorts, our DGP method achieves overall best performance with 0.81 (Yonsei-1), 0.82 (Yonsei-1-remade), 0.89 (Yonsei-2), 0.85 (St. Mary-1) 0.75 (St. Mary-2) AUCs compared to the state-of-the art methods. For the STAD cohorts, our method also achieves overall best performance with 0.74 AUC (Yonsei-classic). In addition, incorporating the uncertainty in prediction measured by our method improves the model’s performance to predict MSI. Specifically, removing cases with prediction results with high uncertainty, which could lead false positive prediction, significantly improves MSI prediction performance. Taken together, these results demonstrate the robustness and generalisability of the DGP model for MSI prediction across multicenter and multiethnic datasets. Citation Format: Sunho Park, Hongming Xu, Sung Hak Lee, Jeonghyun Kang, Tae Hyun Hwang. Deep Gaussian process with uncertainty estimation Improves microsatellite instability prediction based on whole slide image: A retrospective multicenter and multiethnic cohort study [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5012.

Abstract Background: Gastric carcinogenesis has been found to be associated with intragastric microbiome, gastric mucosal gene expression, and immune cells comprising the tumor microenvironment separately in independent studies. However, to date, there has been no multi-omic investigations analyzing these biological traits simultaneously. Thus, there is a paucity in our understanding of the intricate relationship between intragastric microbes and the host in gastric disease. To address this gap, we aimed to comprehensively investigate the microbiome, host gene expression, and immune cells in the stomach across a range of different states along the gastric carcinogenesis pathway. Methods: Gastric biopsy samples were obtained from healthy individuals or patients with gastritis or gastric cancer (biopsy cohort, n=30). Surgical samples (of cancer and adjacent severe gastritis tissues) were obtained from gastric cancer patients (surgery cohort, n=40). Using 16S rRNA gene sequencing, RNA-seq, and cell-type enrichment analysis of each participant’s biospecimens, we identified associations between the microbiome and host transcriptome or inferred immune cell types. Results: Microbiome composition and gene expression patterns showed significant differences between disease states (healthy, gastritis, or cancer). Expression of ‘Cell cycle’ pathway genes was enriched in gastric cancer in both biopsy and surgery cohorts. In the biopsy cohort, Helicobacteraceae, which was abundant in gastritis, was highly correlated with FAM3D, which was previously found to be involved in gastrointestinal inflammation. In the surgery cohort, Lachnospiraceae, which was abundant in gastric cancer, was highly correlated with UBD, which is known to regulate mitosis and reduce cell cycle time. Cell-type enrichment analysis revealed lower B cell infiltration in gastric cancer than in gastritis (both cohorts). Tissue-infiltrating B cells were associated with Helicobacteraceae. Conclusions: Associations between the intragastric microbiome, human gene expression, and tissue-infiltrating immune cells are unique to gastric disease state. The findings herein may aid the discovery of novel biomarkers or therapeutic targets of gastric cancer. Citation Format: Chan Hyuk Park, Changjin Hong, A-reum Lee, Jaeyun Sung, Tae Hyun Hwang. Integrated analysis of intragastric microbiome and human gene expression uncovers genetic, microbial, and immunological associations in gastritis and gastric cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 5013.