Low mutation and neoantigen burden and fewer effector tumor infiltrating lymphocytes correlate with breast cancer metastasization to lymph nodes

Predictive value of cuproptosis and disulfidptosis-related lncRNA in head and neck squamous cell carcinoma prognosis and treatment

Objective

Head and neck squamous cell carcinoma (HNSCC) is a highly lethal and prevalent malignant tumor with a poor prognosis due to its high recurrence rate, This study aims to develop a prognostic index for HNSCC patients based on Cuproptosis and Disulfidptosis-related long noncoding RNA.

Methods

Gene expression and clinical data for HNSCC were obtained from The Cancer Genome Atlas (TCGA). Using Lasso regression and multivariate Cox regression, we established a risk scoring model. The predictive ability of the nomogram, based on clinical features and risk scores, was verified using receiver operating characteristics and calibration curves. We compared independent prognostic parameters, risk score distribution, and survival between high-risk and low-risk groups, followed by preliminary validity evaluations of the model.

Results

Our systematic evaluation of prognostic risk provides a new direction for improving the survival prognosis of HNSCC patients in clinical practice, The model effectively categorized patients into high- and low-risk groups with distinct outcomes, identifying numerous gene mutations in these groups, A low-risk score was associated with a better prognosis and higher survival rates.

Conclusion

The risk score prognostic prediction system developed in this study shows potential efficacy in predicting the prognosis of HNSCC patients and has practical applications in clinical settings.

Fasting-mimicking diet remodels gut microbiota and suppresses colorectal cancer progression

The progression of colorectal cancer is closely associated with diet. Fasting-mimicking diet (FMD) is a promising type of dietary intervention that have beneficial effects in the prevention and treatment of various cancers. We investigated the therapeutic effect of 4-day FMD against colorectal cancer in mice through immune cell analysis, microbiota composition analysis and anti-PD-1 treatment. These FMD cycles effectively suppressed colorectal cancer growth, reduced cell proliferation and angiogenesis, increased tumor-infiltration lymphocytes especially CD8+T cells. FMD stimulated protective gut microbiota, especially Lactobacillus. Supplementation of Lactobacillus johnsonii induced similar results as FMD intervention, which also suppressed tumor growth and increased CD45+ and CD8+ T cells. Additionally, FMD synthesizing with anti-PD-1 therapy effectively inhibited CRC progression. These findings suggest that Lactobacillus. johnsonii is necessary for the anticancer process of FMD in CRC. FMD through its effects on both gut microbiota and immune system, effectively suppressed colorectal cancer progression in mouse model.

Neoantigens and cancer-testis antigens as promising vaccine candidates for triple-negative breast cancer: Delivery strategies and clinical trials

Immunotherapy is gaining prominence as a promising strategy for treating triple-negative breast cancer (TNBC). Neoantigens (neoAgs) and cancer-testis antigens (CTAs) are tumor-specific targets originating from somatic mutations and epigenetic changes in cancer cells. These antigens hold great promise for personalized cancer vaccines, as supported by preclinical and early clinical evidence in TNBC. This review delves into the potential of neoAgs and CTAs as vaccine candidates, emphasizing diverse strategies and delivery approaches. It also highlights the current status of vaccination modalities undergoing clinical trials in TNBC therapy. A comprehensive understanding of neoAgs, CTAs, vaccination strategies, and innovative delivery methods is crucial for optimizing neoAg-based immunotherapies in clinical practice.

CD93 regulates breast cancer growth and vasculogenic mimicry through the PI3K/AKT/SP2 signaling pathway activated by integrin β1

In women, breast cancer (BC) accounts for 7%-10% of all cancer cases and is one of the most common cancers. To identify a new method for treating BC, the role of CD93 and its underlying mechanism were explored. MDA-MB-231 cells were used in this study and transfected with si-CD93, si-MMRN2, oe-CD93, si-integrin β1, or oe-SP2 lentivirus. After MDA-MB-231 cells were transfected with si-NC or si-CD93, they were injected into nude mice by subcutaneous injection at a dose of 5 × 106/mouse to construct a BC animal model. The expression of genes and proteins and cell migration, invasion and vasculogenic mimicry were detected by RT‒qPCR, western blot, immunohistochemistry, immunofluorescence, Transwell, and angiogenesis assays. In pathological samples and BC cell lines, CD93 was highly expressed. Functionally, CD93 promoted the proliferation, migration, and vasculogenic mimicry of MDA-MB-231 cells. Moreover, CD93 interacts with MMRN2 and integrin β1. Knockdown of CD93 and MMRN2 can inhibit the activation of integrin β1, thereby inhibiting the PI3K/AKT/SP2 signaling pathway and inhibiting BC growth and vasculogenic mimicry. In conclusion, the binding of CD93 to MMRN2 can activate integrin β1, thereby activating the PI3K/AKT/SP2 signaling pathway and subsequently promoting BC growth and vasculogenic mimicry.

Genetic alterations in the neuronal development genes are associated with changes of the tumor immune microenvironment in pancreatic cancer

Background

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis and is highly metastatic. Our prior studies have demonstrated the critical role of axon guidance pathway genes in PDAC and the connection between neuronal development and the tumor microenvironment. A recent study newly identified 20 neuronal development genes [disks large homolog 2 (DLG2), neuron-glial-related cell adhesion molecule (NRCAM), neurexin3 (NRXN3), mitogen-activated protein kinase 10 (MAPK10), platelet-derived growth factor D (PDGFD), protein kinase C epsilon (PRKCE), potassium calcium-activated channel subfamily M alpha 1 (KCNMA1), polycystic kidney and hepatic disease 1 (PKHD1), neural cell adhesion molecule 1 (NCAM1), neuregulin-1 (NRG1), zinc finger protein 667 (ZNF667), cystic fibrosis transmembrane conductance regulator (CFTR), acyl-CoA medium-chain synthetase-3 (ACSM3), complement 6 (C6), protein tyrosine phosphatase receptor type M (PTPRM), hypoxia-inducible factor 1 alpha (HIF1A), adenylyl cyclase 5 (ADCY5), adherens junctions-associated protein 1 (AJAP1), neurobeachin (NBEA), sodium voltage-gated channel alpha subunit 9 (SCN9A)] that are associated with perineural invasion and poor prognosis of PDAC. The relationship between genetic alterations in these 20 genes and tumor immune microenvironment (TME) has not previously been investigated.

Methods

We hence applied the sequential multiplex immunohistochemistry results of biopsy specimens from 63 PDAC patients to investigate this relationship.

Results

We found that, except for PTPRM and NBEA, genetic alterations involving these 20 genes are associated with significant changes in the densities of major immune cell subtypes. Except for AJAP1, the copy number loss involving this panel of neuronal development genes is significantly associated with changes in immune cell infiltrates. In contrast, the copy number gain in fewer genes, including NRXN3, ZNF667, ACSM3, C6, ADCY5, SCN9A, and PRKCE, is significantly associated with changes in immune cell infiltrates.

Conclusions

Our study suggested that neuronal development genes play a role in modulating TME in a pancreatic cancer setting.

Loss of TTC17 promotes breast cancer metastasis through RAP1/CDC42 signaling and sensitizes it to rapamycin and paclitaxel

BACKGROUND

Breast cancer (BC) metastasis is the leading cause of poor prognosis and therapeutic failure. However, the mechanisms underlying cancer metastasis are far from clear.

METHODS

We screened candidate genes related to metastasis through genome-wide CRISPR screening and high-throughput sequencing of patients with metastatic BC, followed by a panel of metastatic model assays. The effects of tetratricopeptide repeat domain 17 (TTC17) on migration, invasion, and colony formation ability together with the responses to anticancer drugs were investigated in vitro and in vivo. The mechanism mediated by TTC17 was determined by RNA sequencing, Western blotting, immunohistochemistry, and immunofluorescence. The clinical significance of TTC17 was evaluated using BC tissue samples combined with clinicopathological data.

RESULTS

We identified the loss of TTC17 as a metastasis driver in BC, and its expression was negatively correlated with malignancy and positively correlated with patient prognosis. TTC17 loss in BC cells promoted their migration, invasion, and colony formation capacity in vitro and lung metastasis in vivo. Conversely, overexpressing TTC17 suppressed these aggressive phenotypes. Mechanistically, TTC17 knockdown in BC cells resulted in the activation of the RAP1/CDC42 pathway along with a disordered cytoskeleton in BC cells, and pharmacological blockade of CDC42 abolished the potentiation of motility and invasiveness caused by TTC17 silencing. Research on BC specimens demonstrated reduced TTC17 and increased CDC42 in metastatic tumors and lymph nodes, and low TTC17 expression was linked to more aggressive clinicopathologic characteristics. Through screening the anticancer drug library, the CDC42 inhibitor rapamycin and the microtubule-stabilizing drug paclitaxel showed stronger inhibition of TTC17-silenced BC cells, which was confirmed by more favorable efficacy in BC patients and tumor-bearing mice receiving rapamycin or paclitaxel in the TTC17^Low arm.

CONCLUSIONS

TTC17 loss is a novel factor promoting BC metastasis, that enhances migration and invasion by activating RAP1/CDC42 signaling and sensitizes BC to rapamycin and paclitaxel, which may improve stratified treatment strategies under the concept of molecular phenotyping-based precision therapy of BC.

Newly identified form of phenotypic plasticity of cancer: immunogenic mimicry

Cancer plasticity is now a recognized new hallmark of cancer which is due to disturbances of cell differentiation programs. It is manifested not only in various forms like the best-known epithelial-mesenchymal transition (EMT) but also in vasculogenic and megakaryocytic mimicries regulated by EMT-specific or less-specific transcription factors such as HIF1a or STAT1/2. Studies in the past decades provided ample data that cancer plasticity can be manifested also in the expression of a vast array of immune cell genes; best-known examples are PDL1/CD274, CD47, or IDO, and we termed it immunogenic mimicry (IGM). However, unlike other types of plasticities which are epigenetically regulated, expression of IGM genes are frequently due to gene amplifications. It is important that the majority of the IGM genes are regulated by interferons (IFNs) suggesting that their protein expressions are regulated by the immune microenvironment. Most of the IGM genes have been shown to be involved in immune escape of cancers broadening the repertoire of these mechanisms and offering novel targets for immunotherapeutics.

Derivation and validation of an epigenetic frailty risk score in population-based cohorts of older adults

DNA methylation (DNAm) patterns in peripheral blood have been shown to be associated with aging related health outcomes. We perform an epigenome-wide screening to identify CpGs related to frailty, defined by a frailty index (FI), in a large population-based cohort of older adults from Germany, the ESTHER study. Sixty-five CpGs are identified as frailty related methylation loci. Using LASSO regression, 20 CpGs are selected to derive a DNAm based algorithm for predicting frailty, the epigenetic frailty risk score (eFRS). The eFRS exhibits strong associations with frailty at baseline and after up to five-years of follow-up independently of established frailty risk factors. These associations are confirmed in another independent population-based cohort study, the KORA-Age study, conducted in older adults. In conclusion, we identify 65 CpGs as frailty-related loci, of which 20 CpGs are used to calculate the eFRS with predictive performance for frailty over long-term follow-up.

Multi-Omics Profiling Suggesting Intratumoral Mast Cells as Predictive Index of Breast Cancer Lung Metastasis

Breast cancer lung metastasis has a high mortality rate and lacks effective treatments, for the factors that determine breast cancer lung metastasis are not yet well understood. In this study, data from 1067 primary tumors in four public datasets revealed the distinct microenvironments and immune composition among patients with or without lung metastasis. We used multi-omics data of the TCGA cohort to emphasize the following characteristics that may lead to lung metastasis: more aggressive tumor malignant behaviors, severer genomic instability, higher immunogenicity but showed generalized inhibition of effector functions of immune cells. Furthermore, we found that mast cell fraction can be used as an index for individual lung metastasis status prediction and verified in the 20 human breast cancer samples. The lower mast cell infiltrations correlated with tumors that were more malignant and prone to have lung metastasis. This study is the first comprehensive analysis of the molecular and cellular characteristics and mutation profiles of breast cancer lung metastasis, which may be applicable for prognostic prediction and aid in choosing appropriate medical examinations and therapeutic regimens.

Resistance to Immunotherapy: Mechanisms and Means for Overcoming

Immune checkpoint blockade transformed cancer therapy during the last decade. However, durable responses remain uncommon, early and late relapses occur over the course of treatment, and many patients with PD-L1-expressing tumors do not respond to PD-(L)1 blockade. In addition, while some malignancies exhibit inherent resistance to treatment, others develop adaptations that allow them to evade antitumor immunity after a period of response. It is crucial to understand the pathophysiology of the tumor-immune system interplay and the mechanisms of immune escape in order to circumvent primary and acquired resistance. Here we provide an outline of the most well-defined mechanisms of resistance and shed light on ongoing efforts to reinvigorate immunoreactivity.

Potential Novel Modules and Hub Genes as Prognostic Candidates of Thyroid Cancer by Weighted Gene Co-Expression Network Analysis

Background

Although thyroid cancer (THCA) is one of the most common type of endocrine malignancy, its highly complex molecular mechanisms of carcinogenesis are not completely known.

Materials and Methods

In this study, weighted gene co-expression network analysis (WGCNA) was utilized to construct gene co-expression networks and evaluate the relations between modules and clinical traits to identify potential prognostic biomarkers for THCA patients. RNA-seq data and clinical data were downloaded from The Cancer Genome Atlas (TCGA). Other independent datasets from the Gene Expression Omnibus (GEO) database and the Human Protein Atlas database were performed to validate findings.

Results

Finally, 11 co-expression modules were constructed and four hub genes, CCDC146, SLC4A4, TDRD9 and MUM1L1, were identified and validated statistically, which were considerably interrelated to worse survival of THCA patients.

Conclusion

This research study revealed four hub genes may be considered candidate prognostic biomarkers and potential therapeutic targets for THCA patients in the future.

MAPK10 Expression as a Prognostic Marker of the Immunosuppressive Tumor Microenvironment in Human Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) remains a devastating malignancy worldwide due to lack of effective therapy. The immune-rich contexture of HCC tumor microenvironment (TME) makes this tumor an appealing target for immune-based therapies; however, the immunosuppressive TME is still a major challenge for more efficient immunotherapy in HCC. Using bioinformatics analysis based on the TCGA database, here we found that MAPK10 is frequently down-regulated in HCC tumors and significantly correlates with poor survival of HCC patients. HCC patients with low MAPK10 expression have lower expression scores of tumor infiltration lymphocytes (TILs) and stromal cells in the TME and increased scores of tumor cells than those with high MAPK10 expression. Further transcriptomic analyses revealed that the immune activity in the TME of HCC was markedly reduced in the low-MAPK10 group of HCC patients compared to the high-MAPK10 group. Additionally, we identified 495 differentially expressed immune-associated genes (DIGs), with 482 genes down-regulated and 13 genes up-regulated in parallel with the decrease of MAPK10 expression. GO enrichment and KEGG pathway analyses indicated that the biological functions of these DIGs included cell chemotaxis, leukocyte migration and positive regulation of the response to cytokine–cytokine receptor interaction, T cell receptor activation and MAPK signaling pathway. Protein–protein interaction (PPI) analyses of the 495 DIGs revealed five potential downstream hub genes of MAPK10, including SYK, CBL, VAV1, LCK, and CD3G. Several hub genes such as SYK, LCK, and VAV1 could respond to the immunological costimulatory signaling mediated by the transmembrane protein ICAM1, which was identified as a down-regulated DIG associated with low-MAPK10 expression. Moreover, ectopic overexpression or knock-down of MAPK10 could up-regulate or down-regulate ICAM1 expression via phosphorylation of c-jun at Ser63 in HCC cell lines, respectively. Collectively, our results demonstrated that MAPK10 down-regulation likely contributes to the immunosuppressive TME of HCC, and this gene might serve as a potential immunotherapeutic target and a prognostic factor for HCC patients.

The Multifaceted Effects of Breast Cancer on Tumor-Draining Lymph Nodes

Breast cancer (BC) accounts for significant morbidity and mortality among women worldwide. About one in three patients with breast cancer present with lymph node (LN) metastasis and LN status is one of the most important prognostic predictors in patients with BC. In addition to their prognostic value, LNs initiate adaptive immunity against BC. Yet, BC cells often avoid immune-mediated destruction in LNs. This review provides an overview of the ways by which BC cells modulate LN stromal and hematopoietic cells to promote metastasis and immune evasion.

ADRB1 was identified as a potential biomarker for breast cancer by the co-analysis of tumor mutational burden and immune infiltration

Breast cancer (BRCA) has traditionally been considered as having poor immunogenicity and is characterized by relatively low tumor mutational burden (TMB). Improving immunogenicity may improve the response to clinical immunotherapy of BRCA. However, the relationship between TMB, immune infiltration, and prognosis in BRCA remains unclear. We aimed to explore their interrelations and potential biomarkers. In this study, based on somatic mutation data of BRCA from The Cancer Genome Atlas (TCGA), patients were categorized into high and low TMB groups utilizing the TMB values. CIBERSOFT algorithm indicated significant infiltration of activated partial immune cells in high TMB group. Besides, ADRB1 had been identified as a prognosis-related immune gene in the mutant genes by the combination of the ImmPort database and the univariate Cox analysis. ADRB1 mutation was associated with lower TMB and manifested a satisfactory clinical prognosis. Various database applications (Gene Set Enrichment Analysis, Tumor IMmune Estimation Resource, Connectivity Map, KnockTF) supported the selection of treatment strategies targeting ADRB1. In conclusion, TMB was not an independent prognostic factor for BRCA and high TMB was more likely to activate a partial immune response. ADRB1 was identified as a potential biomarker and may provide new insights for co-therapy of BRCA.

Low mutation burden and differential tumor-infiltrating immune cells correlate with lymph node metastasis in colorectal cancer

BACKGROUND

Tumor immunotherapy has become an important means of cancer treatment. A response depends on the interaction of tumor cells with immune regulators in the tumor microenvironment, which plays an important role in inhibiting or enhancing the immune response. However, lymph node (LN) metastasis leads to major changes in the tumor microenvironment of patients with colorectal cancer, directly affecting prognosis.

METHODS

Using data downloaded from the Cancer Genome Atlas (TCGA) database, we studied the microenvironmental differences between LN-negative and positive populations by bioinformatic methods.

RESULTS

Patients in the LN-positive group had significantly lower immune scores, cytolytic activity scores, and overall survival than the LN-negative group. In addition, a high mutation burden and a new antigen burden could inhibit lymph node metastasis of CRC. In particular, in the LN positive group, the ratio of monocytes to M1 macrophages was significantly downregulated. After the differentially expressed mRNAs between the LN positive and negative groups were determined, a new CRC model was constructed based on multivariate Cox proportional hazard regression analysis to examine the prognosis of patients. The analyses showed that the model was stable and robust.

CONCLUSIONS

We used multiple scores and details of immune cell infiltration as indicators to assess changes in the tumor microenvironment of CRC patients before and after lymph node metastasis, and quantify and model the immune cells in the microenvironment to predict the overall survival of CRC patients.

Cancer neoantigen: Boosting immunotherapy

Tumor neoantigen has a high degree of immunogenicity. As one of the emerging methods of tumor immunotherapy, the vaccine developed against it has served to clinical trials of various solid tumors, especially in the treatment of melanoma. Currently, a variety of immunotherapy methods have been applied to the treatment of the tumor. However, other therapeutic methods have the disadvantages of low specificity and prominent side effects. Treatments require tumor antigen with higher immunogenicity as the target of immune attack. This review will recommend the identification of neoantigen, the influencing factors of neoantigen, and the application of personalized vaccines for neoantigen in metastatic tumors such as malignant melanoma.

Atypical cyclins in cancer: New kids on the block?

Atypical cyclins have recently emerged as a new subfamily of cyclins characterized by common structural features and interactor pattern. Interestingly, atypical cyclins are phylogenetically close to canonical cyclins, which have well-established roles in cell cycle regulation and cancer. Therefore, although the function of atypical cyclins is still poorly characterized, it seems likely that they are involved in cancer pathogenesis as well. Here, we coupled gene expression and prognostic significance analysis to bibliographic search in order to provide new insights into the role of atypical cyclins in cancer. The information gathered suggests that atypical cyclins intervene in critical processes to sustain cancer growth and have potential to become novel prognostic markers and drug targets in cancer.

Immunotherapy in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is considered an immunogenic tumor that arises in chronically inflamed livers due to underlying chronic liver disease caused by viral and non-viral pathogenesis. This inflammation leads to tumor development and is associated to higher tumor immunogenicity. For this reason immunotherapeutic approaches may be suitable therapeutic strategies for HCC. Indeed, several preclinical and clinical data support this hypothesis showing that immunotherapy and even more their combination may be a good alternative candidate for the treatment of HCC patients. However, considering that the liver plays a central role in host defense as well as in the maintenance of self-tolerance, it is characterized by a strong intrinsic immune suppressive microenvironment as well as by a high immune evasion, which may represent a major impediment for an effective immune response against tumor. Furthermore, the low expression of tumor antigens on liver cancer cells leads to a lower T-cell activation and tumor infiltration, resulting in a less efficient control of the tumor growth and, consequently, in a worse clinical outcome. For this reason, strategies should be developed to counteract the different factors in the HCC tumor microenvironment playing a major role in reducing the effects of immunotherapy.

Atypical cyclins: the extended family portrait

Regulation of cell division is orchestrated by cyclins, which bind and activate their catalytic workmates, the cyclin-dependent kinases (CDKs). Cyclins have been traditionally defined by an oscillating (cyclic) pattern of expression and by the presence of a characteristic "cyclin box" that determines binding to the CDKs. Noteworthy, the Human Genome Sequence Project unveiled the existence of several other proteins containing the "cyclin box" domain. These potential "cyclins" have been named new, orphan or atypical, creating a conundrum in cyclins nomenclature. Moreover, although many years have passed after their discovery, the scarcity of information regarding these possible members of the family has hampered the establishment of criteria for systematization. Here, we discuss the criteria that define cyclins and we propose a classification and nomenclature update based on structural features, interactors, and phylogenetic information. The application of these criteria allows to systematically define, for the first time, the subfamily of atypical cyclins and enables the use of a common nomenclature for this extended family.

Immunotherapy: A Challenge of Breast Cancer Treatment

Breast cancer is the most commonly diagnosed cancer in women and is a leading cause of cancer death in women worldwide. Despite the significant benefit of the use of conventional chemotherapy and monoclonal antibodies in the prognosis of breast cancer patients and although the recent approval of the anti-PD-L1 antibody atezolizumab in combination with chemotherapy has been a milestone for the treatment of patients with metastatic triple-negative breast cancer, immunologic treatment of breast tumors remains a great challenge. In this review, we summarize current breast cancer classification and standard of care, the main obstacles that hinder the success of immunotherapies in breast cancer patients, as well as different approaches that could be useful to enhance the response of breast tumors to immunotherapies.

Tumor antigens heterogeneity and immune response-targeting neoantigens in breast cancer

Breast cancer is both the most common type of cancer and the most frequent cause of cancer mortality in women, mainly because of its heterogeneity and limited immunogenicity. The aim of specific active cancer immunotherapy is to stimulate the host's immune response against cancer cells directly using a vaccine platform carrying one or more tumor antigens. In particular, the ideal tumor antigen should be able to elicit T cell and B cell responses, be specific for the tumor and be expressed at high levels on cancer cells. Neoantigens are ideal targets for immunotherapy because they are exclusive to individual patient's tumors, are absent in healthy tissues and are not subject to immune tolerance mechanisms. Thus, neoantigens should generate a specific reaction towards tumors since they constitute the largest fraction of targets of tumor-infiltrating T cells. In this review, we describe the technologies used for neoantigen discovery, the heterogeneity of neoantigens in breast cancer and recent studies of breast cancer immunotherapy targeting neoantigens.

Immunologic and immunogenomic aspects of tumor progression

Tumor progression to metastatic disease is characterized by continuous genetic alterations due to instability of the genome. Immune sensitivity was found to be linked to tumor mutational burden (TMB) and the resulting amount of neoantigens. However, APOBEC activity resulting in increase in TMB causes immune evasion. On the other hand, clonal or acquired genetic loss of HLA class I also hampers immune sensitivity of tumors. Rare amplification of the PD-L1 gene in cancers may render them sensitive to immune checkpoint inhibitors but involvement of broader regions of chromosome 9p may ultimately lead again to immune evasion due to inactivation of the IFN-γ signaling pathway. Such genetic changes may occur not only in the primary tumor but at any phase of progression: in lymphatic as well as in visceral metastases. Accordingly, it is rational to monitor these changes continuously during disease progression similar to target therapies. Moreover, beside temporal variability, genomic features of tumors such as mutation profiles, as well as the tumor immune microenvironment also show considerable inter- and intratumoral spatial heterogeneity, suggesting the necessity of multiple sampling in biomarker studies.