The complexity of the tumor microenvironment poses significant challenges in cancer therapy. Here, to comprehensively investigate the tumor-normal ecosystems, we perform an integrative analysis of 4.9 million single-cell transcriptomes from 1070 tumor and 493 normal samples in combination with pan-cancer 137 spatial transcriptomics, 8887 TCGA, and 1261 checkpoint inhibitor-treated bulk tumors. We define a myriad of cell states constituting the tumor-normal ecosystems and also identify hallmark gene signatures across different cell types and organs. Our atlas characterizes distinctions between inflammatory fibroblasts marked by AKR1C1 or WNT5A in terms of cellular interactions and spatial co-localization patterns. Co-occurrence analysis reveals interferon-enriched community states including tertiary lymphoid structure (TLS) components, which exhibit differential rewiring between tumor, adjacent normal, and healthy normal tissues. The favorable response of interferon-enriched community states to immunotherapy is validated using immunotherapy-treated cancers (n = 1261) including our lung cancer cohort (n = 497). Deconvolution of spatial transcriptomes discriminates TLS-enriched from non-enriched cell types among immunotherapy-favorable components. Our systematic dissection of tumor-normal ecosystems provides a deeper understanding of inter- and intra-tumoral heterogeneity.

Pancreatic ductal adenocarcinoma (PDAC) is a lethal and clinically heterogeneous disease with a limited benefit from human epidermal growth factor receptor 2 (HER2)-targeted therapy. Recently, some studies have addressed the antitumoral effect of novel anti-HER2 drugs in HER2 low-expressing tumors. However, there have been few studies on the significance of low HER2 expression and genetic heterogeneity in PDAC. Using immunohistochemistry and dual-color silver-enhanced in situ hybridization based on the Trastuzumab for a gastric cancer scoring scheme, we evaluated HER2 protein expression, gene amplification, and genetic heterogeneity in three groups (HER2-neg, HER2-low, HER2-pos) of 55 patients. Among the 55 cases, 41.8% (23/55) showed HER2 expression of any intensity. HER2 amplification independent of HER2 expression was 25.5% (14/55). Patients in both these groups had a shorter overall survival than did patients in the HER2-neg group. HER2 genetic heterogeneity was identified in 37 (70.9%) of the 55 cases, mainly in HER2-neg and HER2-low groups. HER2 genetic heterogeneity significantly correlated with worse survival in the HER2-low and HER2-neg groups of PDAC. These findings support the hypothesis that low-level HER2 expression and heterogeneity have significant clinical implications in PDAC. HER2 heterogeneity might indicate the best strategies of combination therapies to prevent the development of subdominant clones with resistance potential.

Ovarian cancer is one of the leading causes of deaths among patients with gynecological malignancies worldwide. In order to identify prognostic markers for ovarian cancer, we performed RNA-sequencing and analyzed the transcriptome data from 51 patients who received conventional therapies for high-grade serous ovarian carcinoma (HGSC). Patients with early-stage (I or II) HGSC exhibited higher immune gene expression than patients with advanced stage (III or IV) HGSC. In order to predict the prognosis of patients with HGSC, we created machine learning-based models and identified USP19 and RPL23 as candidate prognostic markers. Specifically, patients with lower USP19 mRNA levels and those with higher RPL23 mRNA levels had worse prognoses. This model was then used to analyze the data of patients with HGSC hosted on The Cancer Genome Atlas; this analysis validated the prognostic abilities of these two genes with respect to patient survival. Taken together, the transcriptome profiles of USP19 and RPL23 determined using a machine-learning model could serve as prognostic markers for patients with HGSC receiving conventional therapy.

The Ki-67 labeling index (LI) is an important prognostic factor in breast carcinoma. The Ki-67 LI is traditionally calculated via unaided microscopic estimation; however, inter-observer and intra-observer variability and low reproducibility are problems with this visual assessment (VA) method. For more accurate assessment and better reproducibility with Ki-67 LI, digital image analysis was introduced recently. We used both VA and automated digital image analysis (ADIA) (Ventana Virtuoso image management software) to estimate Ki-67 LI for 997 cases of breast carcinoma, and compared VA and ADIA results. VA and ADIA were highly correlated (intraclass correlation coefficient 0.982, and Spearman’s correlation coefficient 0.966, p<0.05). We retrospectively analyzed cases with a greater than 5% difference between VA and ADIA results. The cause of these differences was: (1) tumor heterogeneity (98 cases, 56.0%), (2) VA interpretation error (32 cases, 18.3%), (3) misidentification of tumor cells (26 cases, 14.9%), (4) poor immunostaining or slide quality (16 cases, 9.1%), and (5) Estimation of non-tumor cells (3 cases, 1.7%). There were more discrepancies between VA and ADIA results in the group with a VA value of 10–20% compared to groups with <10% and ≥20%. Although ADIA is more accurate than VA, there are some limitations. Therefore, ADIA findings require confirmation by a pathologist.

Alteration in expression of miRNAs can cause various malignant changes and the metastatic process. Our aim was to identify the miRNAs involved in cervical squamous cell carcinoma (SqCC) and metastasis, and to test their utility as indicators of metastasis and survival. Using microarray technology, we performed miRNA expression profiling on primary cervical SqCC tissue (n = 6) compared with normal control (NC) tissue and compared SqCC that had (SqC-M; n = 3) and had not (SqC-NM; n = 3) metastasized. Four miRNAs were selected for validation by qRT-PCR on 29 SqC-NM and 27 SqC-M samples, and nine metastatic lesions (ML-SqC), from a total of 56 patients. Correlation of miRNA expression and clinicopathological parameters was analyzed to evaluate the clinical impact of candidate miRNAs. We found 40 miRNAs differentially altered in cervical SqCC tissue: 21 miRNAs were upregulated and 19 were downregulated (≥2-fold, p < 0.05). Eight were differentially altered in SqC-M compared with SqC-NM samples: four were upregulated (miR-494, miR-92a-3p, miR-205-5p, and miR-221-3p), and four were downregulated (miR-574-3p, miR-4769-3p, miR-1281, and miR-1825) (≥1.5-fold, p < 0.05). MiR-22-3p might be a metastamiR, which was gradually further downregulated in SqC-NM > SqC-M > ML-SqC. Downregulation of miR-30e-5p significantly correlated with high stage, lymph node metastasis, and low survival rate, suggesting an independent poor prognostic factor.

Lee GH, et al. Yonsei Med J. 2021 Apr;62(4):288-297. https://doi.org/10.3349/ymj.2021.62.4.288