Impact of microvessel patterns and immune status in NSCLC: a non-angiogenic vasculature is an independent negative prognostic factor in lung adenocarcinoma

Introduction: Non-small cell lung carcinomas (NSCLC) exhibit different microvessel patterns (MVPs). Basal (BA), diffuse (DA) and papillary (PA) patterns show signs of angiogenesis (new blood vessels), while an alveolar pattern indicates that tumors are co-opting existing normal vessels (non-angiogenic alveolar, NAA). NAA tumor growth is known to exist in NSCLC, but little is known about its prognostic impact in different histological subgroups, and about associations between MVPs and immune cell infiltration.

Methods: Detailed patterns of angiogenic and non-angiogenic tumor growth were evaluated by CD34 immunohistochemistry in whole tissue slides from 553 surgically treated patients with NSCLC stage I-IIIB disease. Associations with clinicopathological variables and markers related to tumor immunology-, angiogenesis- and hypoxia/metabolism were explored, and disease-specific survival (DSS) was analyzed according to histological subtypes.

Results: The predominant MVP was angiogenic in 82% of tumors: BA 40%, DA 34%, PA 8%, while a NAA pattern dominated in 18%. A contribution of the NAA pattern >5% (NAA+), i.e., either dominant or minority, was observed in 40.1% of tumors and was associated with poor disease-specific survival (DSS) (p=0.015). When stratified by histology, a significantly decreased DSS for NAA+ was found for adenocarcinomas (LUAD) only (p< 0.003). In multivariate analyses, LUAD NAA+ pattern was a significant independent prognostic factor; HR 2.37 (CI 95%, 1.50-3.73, p< 0.001). The immune cell density (CD3, CD4, CD8, CD45RO, CD204, PD1) added prognostic value in squamous cell carcinoma (LUSC) and LUAD with 0-5% NAA (NAA-), but not in LUAD NAA+. In correlation analyses, there were several significant associations between markers related to tumor metabolism (MCT1, MCT4, GLUT1) and different MVPs.

Conclusion: The NAA+ pattern is an independent poor prognostic factor in LUAD. In NAA+ tumors, several immunological markers add prognostic impact in LUSC but not in LUAD.

Introduction

Angiogenesis is the process by which new blood vessels are formed. The induction of angiogenesis, in order to ensure sufficient supply of oxygen and nutrients in growing tumors, was one of the original Hallmarks of cancer (1). However, in the recent update “Hallmark of cancer: new dimensions” by Hanahan, “angiogenesis” was finally modified to “inducing or accessing vasculature” (1–3). After years of increasing reports and a long-lasting debate in the cancer community, there is now an appreciation that tumors can become sufficiently vascularized either by “switching on” angiogenesis or by co-opting normal tissue vessels (2). We and others have previously reviewed convincing data concluding that many tumors, partly or completely, grow without angiogenesis (4–6). Non-angiogenic tumors may also be a cause of resistance to anti-angiogenic therapy (4–6).

Non-small cell lung cancer is a common and deadly disease (7). An emerging understanding of the tumor immune microenvironment, and the introduction of immunotherapy, has led to improved outcomes for many non-small cell lung cancer (NSCLC) patients (8). Yet, the majority of patients progress due to resistance mechanisms. There is a complex interplay between tumor blood supply and tumor immunity, and potential clinical implications are being explored, including treatments combining angiogenesis inhibitors and immunotherapy (9, 10). Hence, there is an increasing need for knowledge of the biology of this interplay, and of the significance of different subtypes of angiogenic or non-angiogenic vessel patterns.

Along with kidneys, liver and brain, lungs are extensively vascularized and non-angiogenic tumor growth is frequently reported (4, 5). In NSCLC, early studies have demonstrated that as many as 10-30% of tumors grow in a non-angiogenic fashion (4, 5, 11). NSCLC data also indicate that an angiogenic microvessel pattern (MVP) may induce a more extensive immunological response when compared to a non-angiogenic MVP, while tumors with co-opted blood supply were associated efficient mitochondrial metabolism (4, 5, 12).

Lung adenocarcinomas (LUAD) and squamous cell carcinomas (LUSC) may be considered two different entities, distinctive both in biological characteristics and treatment approaches. We wanted to explore the prognostic impact of MVPs stratified by histology, and whether combinations with immunological markers added prognostic impact. We have previously explored multiple markers related to angiogenesis, hypoxia/metabolism and tumor immunology in the same cohort. Herein we wanted to examine associations between MVPs and these markers, with emphasis on immunological markers representing essential innate and adaptive immune cells (13–26).

Material and methods

Study design and patients

Tumor specimens from consecutive patients undergoing radical resection for NSCLC, pathological stage I to IIIB, at the University Hospital of North Norway (UNN, n=295) and the Nordland Hospital (NH, n=258) from 1990 through 2010, were collected retrospectively. In total, 553 patients with complete medical records were included, as previously reported (20). Follow-up was completed on October 1, 2013, and median follow-up of survivors was 86 months (range, 34-267 months).

Tumors were restaged according to the 8th International Union Against Cancer TNM classification, and histologically re-evaluated according to the 2015 WHO classification (27, 28). Patients with missing or inadequate tissue quality on stained slides (n=69) were excluded, leaving 484 patients available for analyses. Each tumor was represented by a single whole slide image, and adjacent hematoxylin and eosin-stained sections were available for morphologic review.

The study was approved by the Norwegian Data Protection Authority and the Regional Committee for Medical and Health Research Ethics (protocol ID, 2011/2503). Clinicopathologic variables, survival data and biomarker expressions are reported in accordance with the REMARK (reporting recommendations for tumor marker prognostic studies) guidelines (29).

Immunohistochemistry CD34

IHC staining was performed using the Discovery Ultra autostainer (Ventana, Roche, Tucson, AZ) on formalin-fixed paraffin-embedded whole tissue sections. The slides were deparaffinized in the instrument (68°C, 3 cycles á 12 min). Antigen retrieval was applied using ULTRA Cell Conditioning-1 (Ventana, Roche) for 32 minutes at 95°C. The sections were incubated for 32 minutes with IVD primary mouse monoclonal antibody anti-CD34 (clone QBEnd/10; cat# 790-2927; Ventana, Roche). As a secondary antibody, OmniMap anti-Mouse HRP (cat# 760-4310, Ventana, Roche) was loaded for 16 minutes, followed by 4 minutes of HQ HRP amplification (cat# 760-052, Ventana, Roche). The detection chromogen used was the DAB kit (cat#760-159; Ventana, Roche) with 32 minutes of incubation. Counterstaining was performed using hematoxylin II (cat#790-2208, Ventana, Roche) for 12 minutes and then loading bluing reagent for 4 minutes.

Slides were digitized using a Pannoramic 250 Flash III (3DHistech, Budapest, Hungary), slide scanner.

Markers of tumor immunity, angiogenesis and hypoxia/metabolism

We have previously reported the distribution and prognostic impact of >100 markers of angiogenesis, immunology, hypoxia/metabolism and EMT expressed in the NSCLC tumor microenvironment, in a tissue microarray (TMA) material in the same patient cohort. Markers were analyzed by immunohistochemistry or in situ hybridization and scored in a manual, semi-quantitative or digital, quantitative fashion, as previously published (13–18, 20–24, 30–49). In the present study, we included data for the most relevant markers involved in tumor immunity (n=22), angiogenesis (n=30) and hypoxia/metabolism (n=14) (Supplementary Table S1). Markers were evaluated in both stromal (_S) and tumor epithelial (_T) compartments, separately or combined (_TS). Initially markers were analyzed in a TMA material representing a cohort of 335 patients surgically treated between 1990 and 2005. Additional analyzes, including most immunology markers, followed a cohort expansion by 218 patients treated between 2005 and 2011 (553 in total).

Microvessel pattern evaluation

All tissue slides were digitally reviewed by an experienced pathologist (FP) and an oncologist (EP) using QuPath software v.0.1.2 (Queen’s University, Belfast, Northern Ireland) (50). Microvessel patterns (MVPs) were identified by CD34 IHC staining and quantified as a percentage of the tumor tissue represented in the whole slide image. The MVPs were defined as angiogenic, including the basal (BA), diffuse (DA) and papillary (PA) subtypes, or non-angiogenic alveolar (NAA), in accordance with Pezzella et al, 1997 (11). Detailed descriptions of the MVPs are described in Figure 1. The contribution of each of the four MVPs was manually estimated in 5% increments; all MVPs could be represented in each tumor.

FIGURE 1

Figure 1 Description of microvessel patterns (MVPs) and their quantification. Microvessels were identified by CD34 IHC DAB staining (brown) of endothelial cells. (A) Basal angiogenic (BA) microvessel pattern (MVP). Tumor cells are arranged in epithelial nests surrounded by connective tissue. Most vessels are in the connective tissue beneath the neoplastic epithelium. (B) Diffuse angiogenic (DA) MVP. Normal lung architecture is diffusely replaced. New vessels, stroma and neoplastic cells are mixed and without recognizable architectural structure. (C) Papillary angiogenic (PA) MVP. Vessels are present in a stromal stalk covered by neoplastic cells, mostly as a monolayer. Normal lung architecture is recognizable in some areas. Remodeling of the alveolar structure occurs. (D) Non-angiogenic alveolar (NAA) MVP Vessels arise from alveolar septa only, normal lung tissue architecture is preserved (“chicken wire”). Tumor cells grow in a solid fashion filling alveolar spaces. No endothelial or stromal cells present among neoplastic cells. Anthracitic pigment may be trapped in the parenchyma. (E) Normal lung. The contribution of each of the four different MVPs was quantified as a percentage of the tumor tissue represented in the whole slide image, estimated in 5% increments.

Statistical methods

Statistical analyses were performed using IBM SPSS Statistics (Version 29) or RStudio with R version 4.1.2 (utilizing packages survival, surminer, ggplot, dplyr, tidyverse, ggcorrplot, ggstatplot). The χ² test or Fischer’s exact test was used to investigate the associations between clinicopathological variables and MVPs. Associations between MVPs and other markers (n=123) were evaluated by Spearman’s rank correlation (all non-parametric variables). The contribution of each MVP (continuous, interval variable) was tested against the mean value (categorical, ordinal ranked variable) of all markers. To adjust for multiple testing, Bonferroni correction was conducted. Disease-specific survival (DSS) was defined as the time from surgery to lung cancer death. Univariate analysis of survival according to MVPs was evaluated by the Kaplan-Meier method. The log-rank test was applied to assess statistically significant differences between survival curves. For survival analyses, the cutoff for dichotomization of MVP variables was selected by applying an “optimal p-value” approach. Cox proportional hazards ratios were calculated. Multivariate analysis included relevant covariates including, but not limited to, variables found significant in univariate analyses. The backward conditional method was used for model fitting. Probability for stepwise entry and removal was set at 0.05 and 0.10, respectively. For statistical tests, a significance level <0.05 was considered statistically significant.

Results

Patient characteristics

Demographic, clinical, and histopathologic variables and their associations with DSS are presented in Table 1. The median age was 67 years (range, 28-85), 33% were female. Seventy-six patients (14%) received postoperative radiotherapy because of pN2 disease or nonradical surgical margins. Following its introduction into Norwegian national guidelines in 2005, 43 patients (8%) received adjuvant systemic therapy. None of the patients received immune checkpoint inhibitors or antiangiogenic treatment during follow-up. Clinicopathological variables did not differ significantly between patients included in the MVP analysis (scored) and patients where no adequate tissue was available (missing), except for frequencies of pneumonectomy and vascular infiltration (24 vs 36% and 16 vs 29%).

TABLE 1

Table 1 Clinicopathological variables and non-angiogenic (NAA) microvessel pattern (MVP) as predictors of disease-specific survival in resected NSCLC patients and in LUSC and LUAD histological subgroups (univariate analyses; log rank test; unadjusted proportional hazard ratios).

Microvessel pattern distribution

We observed two or more MVPs in 316 of the 484 tumors. In the 170 cases where a single pattern covered the entire tumor, an angiogenic pattern was more frequently found (BA n=79, DA n=65, PA n=11) than a non-angiogenic (NAA n=15). Any presence (5-100% of the tumor area) of the BA, DA, PA, and NAA MVPs was observed in 257, 274, 68 and 260 of the 484 tumors, respectively. The predominant MVP in each tumor represented from 40% (in tumors with multiple MVPs) to 100% of the area. The predominant MVP was of an angiogenic subtype in 82% (BA 40%, DA 34%, PA 8%) of cases, while the NAA pattern was predominant in 18% of tumors. The distribution of MVPs in tumors is presented in detail Supplementary Figure S1.

Correlations analyses

Correlation matrices presented in Supplementary Figure S2 summarize all correlations between MVPs and markers related to immunology, angiogenesis and hypoxia/metabolism. After Bonferroni correction for multiple testing, p-values <.0004 were considered significant. Correlation coefficients >.2 were regarded as relevant in this explorative setting. A selection of significant correlations (p<0.01) between MVPs and included markers are presented in Table 2.

TABLE 2

Table 2 Correlations between MVP and other markers.

Correlations with markers related to tumor immunity

Few immune markers were significantly correlated with MVPs. CD138 (plasma cells) in tumor and stroma was positively correlated with BA (.293-.298). For DA MVP, positive correlations with CSF1R were found (monocyte and myeloid dendritic cells) (.231).

Correlations with markers related to angiogenesis

A previous assessment of microvessel density, also by CD34 staining, showed a positive correlation with NAA (_T:.308, _S:.215). Except for PDGFRβ_S (-.205 in NAA) and FGFR1_T (-.307 in BA and.210 in PA), markers of angiogenesis (VEGF-A, -C, -D, VEGFR-1, -R2, -R3, PDGF-A, -B, -C, -D, PDGFR-α, bFGF, Notch1, -4, Jagged1, DLL4, Ang1, -2, -4, Tie2), were not significantly correlated with MVPs.

Correlations with markers related to hypoxia and metabolism

Markers involved in hypoxia and metabolism were divergently expressed in different MVPs. A significant positive correlation with the BA MVP was found for MCT1_T (.340), and GLUT1_T (.426), and a negative correlation for MCT4_T (-.268). Oppositely, significant negative correlations were found for GLUT1_T in PA (-.282), MCT4_S in NAA (-.247) and for MCT1_ in both PA (-.339) and NAA (-.232).

Associations with clinicopathological variables

Tumors with a predominant BA MVP were significantly associated with male sex and LUSC, and the NAA pattern with female sex and with LUAD histology. No other associations with other clinicopathological variables were found.

Univariate survival analyses

Pathological TNM stage, T- and N-stage, tumor differentiation, vascular infiltration, surgical procedure, sex and ECOG status were significant prognostic indicators for DSS in all patients (all p< 0.001) (Table 1).

The angiogenic MVPs (BA, DA or PA) were not significantly associated with DSS, when analyzed separately as dichotomized variables (none vs any for BA, DA or PA) or as predominant patterns (Supplementary Table S2).

The NAA MVP was associated with poor prognosis in all patients when applying an optimal cutoff of 0-5% (NAA-) vs >5% (NAA+); (HR 1.42, 95%CI 1.07-1.90, p=.016). When stratified by histological subgroup, there was no significant difference in DSS associated with the NAA MVP in patients with LUSC, while for LUAD patients, DSS for the NAA+ subgroup was significantly decreased (HR 1.88, 95%CI 1.24-2.85, p=0.003).

Results of survival analyses for NAA are presented in Table 1 and Figure 2.

FIGURE 2

Figure 2 Disease specific survival curves for NAA MVP in (A) all patients, (B) LUAD and (C) LUSC subgroups.

Similar associations with poor DSS for the NAA pattern in LUAD patients, not LUSC, was found for other exploratory cutoffs, including none vs any NAA (HR 1.52, 95%CI 0.97-2.36, p=0.058) and 0-10% vs >10% NAA (HR 1.70, 95%CI 1.14-2.53, p=0.009) as well as for predominant NAA MVP vs angiogenic MVPs (HR 1.65, 1.1-2.48, p=.017) (Supplementary Table S2).

Univariate analyses of the prognostic impact of immune cell density combined with NAA MVP status on DSS was assessed in all patients and LUSC and LUAD subgroups. For the LUSC subgroup, high density of CD8_TS contributed to a significantly improved DSS for patients with both NAA+ and NAA- tumors. For LUAD tumors, a high immune cell density added no prognostic impact for NAA+ tumors (Figure 3). Similar results are also shown for CD3, CD4, CD45RO, PD1, CD204 (Supplementary Figure S3).

FIGURE 3

Figure 3 Disease specific survival curves for NAA MVP combined with CD8 density in all (A) patients, (B) LUAD and (C) LUSC subgroups.

Multivariable analysis

Results from the multivariable Cox regression analysis are presented in Table 3. NAA+ MVP was identified as an independent negative prognostic factor for DSS (HR 1.67, 95%CI 1.24-2.24, p=.0006). In histological subgroup analysis, the prognostic impact was limited to LUAD patients (HR 2.37 95%CI 1.50-3.73, p=.0002). The prognostic impact was significant for LUAD patients also when applying other exploratory cutoffs for NAA (none vs any NAA, HR 1.84, 95%CI 1.14-2.96, p=0.012; 0-10% vs >10% NAA, HR 1.78, 95%CI 1.21-2.83, p=0.015; predominant NAA vs angiogenic MVPs, HR 2.09, 95%CI 1.34-3.27, p=0.001) (Supplementary Table S3).

TABLE 3

Table 3 Results of Cox Regression analysis summarizing significant independent prognostic factors for DSS in all patients and LUSC and LUAD histological subgroups.

Discussion

In our cohort of resected NSCLC patients, a detailed evaluation of the distribution of angiogenic and non-angiogenic microvessel patterns (MVPs) was performed. Most tumors presented multiple MVPs, and a predominant non-angiogenic alveolar (NAA) MVP was rare (18%). A NAA pattern was more frequently observed in LUAD than LUSC tumors. The NAA pattern was associated with poor prognosis, even when it occupied a limited proportion of the tumor area (>5%, NAA+), in 40% of tumors. In multivariate analysis, the NAA+ pattern was a significant independent negative prognostic factor, and when stratified by histological subgroups, the prognostic effect was limited to patients with LUAD. The immune cell density (CD3, CD4, CD8, CD45RO, PD1, CD204) added prognostic value for all LUSC tumors and in NAA- LUAD tumors. No prognostic effect, however, was found for immune cell density in NAA+ LUAD tumors. In addition, markers associated with metabolism (MCT1, MCT4 and GLUT1) was divergently associated between the four MVPs.

The main weakness of this study is the explorative design. Associations between MVPs and molecular markers were analyzed by multiple testing and without a predefined hypothesis. Adjustment for multiple testing was performed (Bonferroni) in order to avoid type I errors in correlation analyses. Yet, we may describe significant associations but must be careful suggesting causality. Further experimental studies are needed to validate potential causality. Applying an “optimal cut-off” approach for dichotomization of variables in survival analyses introduces a risk of false positive results, but for NAA several cutoffs yielded significant results, in both uni- and multivariate analyses for LUAD. Another potential weakness arises when comparing whole slide scoring (MVPs, TILs) with TMA (other markers) scoring approaches. In our experience, the potential pitfall of tumor heterogeneity is largely avoided when applying mean scores of multiple TMA cores from each tumor (24).

The use of an unselected cohort strengthens the study. Further, an experienced pathologist (FP) in the field of MVPs was in charge of the scoring, and the frequencies of the different MVPs are similar to previous reports (11). Finally, being able to explore associations with a variety of molecular markers previously analyzed in the same tumors is a great asset to the study.

Results from previously published studies reporting the prognostic impact of the NAA pattern in NSCLC have been conflicting. Some studies, however, have indicated that the NAA pattern is associated with a poor prognosis (4, 51, 52). To our knowledge, this is the first study identifying NAA as an independent negative prognostic factor, limited to the LUAD subgroup. As a consequence, one may consider whether patients with LUAD tumors with a NAA MVP are in special need of close follow-up or adjuvant therapy.

This is also the first study to evaluate the combined prognostic impact of MVP and immunological markers. Intriguingly, when a NAA MVP was present (>5%), several immunological markers added prognostic impact in LUSC, but not in LUAD tumors (CD3, CD4, CD8, CD45RO, PD1, CD204). A biological rationale which explains these differences in LUAD vs LUSC can, based on our results, only be speculative. However, a search for immunological biomarkers of prognostic impact for NSCLC patients is ongoing (53, 54). Our group has previously found that immune cells, including CD8, are strong predictors of survival, especially for patients with resected LUSC tumors (16, 18). The present results indicate that some immunological markers, including CD8, may not be of prognostic significance in resected LUAD with NAA+. On the contrary, in NAA- LUAD tumors, a high CD8 has a strong prognostic impact. This may have clinical implications when optimizing the use of CD8 in a resected NSCLC immune score setting.

It has previously been speculated that NAA tumors are less immunogenic compared to DA tumors (4, 5). In line with this hypothesis, we observe several innate and adaptive immune markers with weak associations with a DA MVP. In contrast, CD8, CSF1R and CD138 are weakly inversely correlated to a NAA+ MVP (12). However, after adjustment for multiple testing these results are not significant. Hence, further investigations may improve our understanding of the immunogenic status in NAA MVP. In light of the above-mentioned lack of prognostic impact of immune cells in LUAD NAA+, one may speculate that neither angiogenesis inhibition nor immunotherapy are optimal treatment strategies in NAA+ LUAD patients.

Interestingly, some of the stronger significant associations were between markers related to metabolism (MCT1, MCT4 and GLUT1) and MVPs. Notably, significant positive correlations are found only with the BA pattern. In hypoxic conditions, a metabolic adaption to produce energy by glycolysis occurs in normal cells (42). Cancer cells, on the other hand, seem to prefer glycolysis despite the presence of oxygen (Warburg effect) (42). In this complex process, lactate and proton transmembrane symporters, MCT1 and MCT4, are pivotal. They prevent intracellular acidification and shuffle the lactate as an energy substrate for the respiratory chain in more oxygenated tissue (42). Although we simply observe associations, it is tempting to hypothesize that glycolysis and degree of hypoxia in the tumor microenvironment, or at least the regulation of these metabolic markers, differs significantly between MVPs. The positive correlation between GLUT1 and MCT1 and negative MCT4 correlation with BA pattern points to a very metabolic active tissue as GLUT1 and MCT1 are energy importers. However, as there is a tendency (though not significant when adjusted for multiple testing) towards an association with CAIX (endogen hypoxic marker) the BA tissue may be more hypoxic. This seems intuitive in a morphological sense, as the distance between vessels is often greater than in the other MVPs (as illustrated by an inverse correlation with microvessel density, MVD). From a clinical point of view, hypoxic tumors are known to be less sensitive to radiotherapy. Hence, it would be interesting to explore whether there is a relationship between metabolic markers, MVPs and radiotherapy response.

In conclusion, NSCLC tumors display different MVPs and several patterns often coexist. NAA+ MVP is an independent negative prognostic factor in LUAD, but not LUSC tumors. Several immunological markers add prognostic impact in LUSC and LUAD tumors with NAA- MVP, but not in LUAD with NAA+ MVP. Further investigation is needed to interpret these novel results as we need further knowledge regarding the relationship between tumor characteristics and AI and IO efficiency.

Our next step is to validate the prognostic impact of MVPs in combination with immunological markers in our comprehensive multicenter prospective NSCLC trial (NCT03299478). A predefined, potentially digital, scoring approach is being considered, as well as including gene expression and mutation data in multivariate analyses. Further studies are needed to increase biological understanding and to potentially establish MVPs, in combination with immunological status, as clinically implementable predictive factors in NSCLC treatment.

Data availability statement

Most relevant data supporting the conclusions of this article are included in Supplementary Figures/Tables. The complete datasets will be made available by the authors upon reasonable request. Requests to access the datasets should be directed to Tom Donnem (tom.donnem@uit.no).

Ethics statement

The studies involving human participants were reviewed and approved by Northern Norway Regional Ethical Committee. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

E-EP, FP, and TD conceptualized and designed the study. E-EP, SA, MR, MIP, AL, MP, TK, L-TB, FP and TD were responsible or contributed to acquisition, analysis, or interpretation of data. E-EP and MR performed the statistical analysis. MIP and AL provided technical support. TD obtained funding. E-EP and TD drafted the manuscript. All authors contributed to the article and approved the submitted version.

Funding

The Norwegian Cancer Society and (Project number 197879-2018) and Northern Norway Health Region Authority (Project number HNF 1378-17) have contributed to this study. Funding for open access publication is provided by UiT, The Artic University of Norway.

Acknowledgments

We want to thank the Norwegian Cancer Society and Northern Norway Health Region Authority for financial support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1157461/full#supplementary-material

References

1. Hanahan D, Weinberg RA. The hallmarks of cancer. Cell (2000) 100(1):57–70. doi: 10.1016/S0092-8674(00)81683-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Hanahan D. Hallmarks of cancer: new dimensions. Cancer Discov (2022) 12(1):31–46. doi: 10.1158/2159-8290.CD-21-1059

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell (2011) 144(5):646–74. doi: 10.1016/j.cell.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Donnem T, Hu J, Ferguson M, Adighibe O, Snell C, Harris AL, et al. Vessel co-option in primary human tumors and metastases: an obstacle to effective anti-angiogenic treatment? Cancer Med (2013) 2(4):427–36. doi: 10.1002/cam4.105

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Donnem T, Reynolds AR, Kuczynski EA, Gatter K, Vermeulen PB, Kerbel RS, et al. Non-angiogenic tumours and their influence on cancer biology. Nat Rev Cancer (2018) 18(5):323–36. doi: 10.1038/nrc.2018.14

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Kuczynski EA, Vermeulen PB, Pezzella F, Kerbel RS, Reynolds AR. Vessel co-option in cancer. Nat Rev Clin Oncol (2019) 16(8):469–93. doi: 10.1038/s41571-019-0181-9

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. CA Cancer J Clin (2022) 72(1):7–33. doi: 10.3322/caac.21708

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Otano I, Ucero AC, Zugazagoitia J, Paz-Ares L. At The crossroads of immunotherapy for oncogene-addicted subsets of NSCLC. Nat Rev Clin Oncol (2023) 20(3):143–59. doi: 10.1038/s41571-022-00718-x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Khan KA, Kerbel RS. Improving immunotherapy outcomes with anti-angiogenic treatments and vice versa. Nat Rev Clin Oncol (2018) 15(5):310–24. doi: 10.1038/nrclinonc.2018.9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhao Y, Guo S, Deng J, Shen J, Du F, Wu X, et al. VEGF/VEGFR-targeted therapy and immunotherapy in non-small cell lung cancer: targeting the tumor microenvironment. Int J Biol Sci (2022) 18(9):3845–58. doi: 10.7150/ijbs.70958

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Pezzella F, Pastorino U, Tagliabue E, Andreola S, Sozzi G, Gasparini G, et al. Non-small-cell lung carcinoma tumor growth without morphological evidence of neo-angiogenesis. Am J Pathol (1997) 151(5):1417–23.

PubMed Abstract | Google Scholar

12. Hu J, Bianchi F, Ferguson M, Cesario A, Margaritora S, Granone P, et al. Gene expression signature for angiogenic and nonangiogenic non-small-cell lung cancer. Oncogene (2005) 24(7):1212–9. doi: 10.1038/sj.onc.1208242

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Al-Shibli K, Al-Saad S, Andersen S, Donnem T, Bremnes RM, Busund LT. The prognostic value of intraepithelial and stromal CD3-, CD117- and CD138-positive cells in non-small cell lung carcinoma. APMIS (2010) 118(5):371–82. doi: 10.1111/j.1600-0463.2010.02609.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Al-Shibli K, Al-Saad S, Donnem T, Persson M, Bremnes RM, Busund LT. The prognostic value of intraepithelial and stromal innate immune system cells in non-small cell lung carcinoma. Histopathology (2009) 55(3):301–12. doi: 10.1111/j.1365-2559.2009.03379.x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Al-Shibli KI, Donnem T, Al-Saad S, Persson M, Bremnes RM, Busund LT. Prognostic effect of epithelial and stromal lymphocyte infiltration in non-small cell lung cancer. Clin Cancer Res (2008) 14(16):5220–7. doi: 10.1158/1078-0432.CCR-08-0133

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Donnem T, Hald SM, Paulsen EE, Richardsen E, Al-Saad S, Kilvaer TK, et al. Stromal CD8+ T-cell density-a promising supplement to TNM staging in non-small cell lung cancer. Clin Cancer Res (2015) 21(11):2635–43. doi: 10.1158/1078-0432.CCR-14-1905

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hald SM, Rakaee M, Martinez I, Richardsen E, Al-Saad S, Paulsen EE, et al. LAG-3 in non-small-cell lung cancer: expression in primary tumors and metastatic lymph nodes is associated with improved survival. Clin Lung Cancer (2018) 19(3):249–59.e2. doi: 10.1016/j.cllc.2017.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kilvaer TK, Paulsen EE, Andersen S, Rakaee M, Bremnes RM, Busund LR, et al. Digitally quantified CD8+ cells: the best candidate marker for an immune cell score in non-small cell lung cancer? Carcinogenesis (2020) 41(12):1671–81. doi: 10.1093/carcin/bgaa105

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Kilvaer TK, Paulsen EE, Khanehkenari MR, Al-Saad S, Johansen RM, Al-Shibli K, et al. The presence of intraepithelial CD45RO+ cells in resected lymph nodes with metastases from NSCLC patients is an independent predictor of disease-specific survival. Br J Cancer (2016) 114(10):1145–51. doi: 10.1038/bjc.2016.92

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Paulsen EE, Kilvaer T, Khanehkenari MR, Maurseth RJ, Al-Saad S, Hald SM, et al. CD45RO(+) memory T lymphocytes–a candidate marker for TNM-immunoscore in squamous non-small cell lung cancer. Neoplasia (2015) 17(11):839–48. doi: 10.1016/j.neo.2015.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Paulsen EE, Kilvaer TK, Khanehkenari MR, Al-Saad S, Hald SM, Andersen S, et al. Assessing PDL-1 and PD-1 in non-small cell lung cancer: a novel immunoscore approach. Clin Lung Cancer (2017) 18(2):220–33.e8. doi: 10.1016/j.cllc.2016.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Paulsen EE, Kilvaer TK, Rakaee M, Richardsen E, Hald SM, Andersen S, et al. CTLA-4 expression in the non-small cell lung cancer patient tumor microenvironment: diverging prognostic impact in primary tumors and lymph node metastases. Cancer Immunol Immunother (2017) 66(11):1449–61. doi: 10.1007/s00262-017-2039-2

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rakaee M, Busund LR, Jamaly S, Paulsen EE, Richardsen E, Andersen S, et al. Prognostic value of macrophage phenotypes in resectable non-small cell lung cancer assessed by multiplex immunohistochemistry. Neoplasia (2019) 21(3):282–93. doi: 10.1016/j.neo.2019.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Rakaee M, Busund LT, Paulsen EE, Richardsen E, Al-Saad S, Andersen S, et al. Prognostic effect of intratumoral neutrophils across histological subtypes of non-small cell lung cancer. Oncotarget (2016) 7(44):72184–96. doi: 10.18632/oncotarget.12360

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Rakaee M, Kilvaer TK, Dalen SM, Richardsen E, Paulsen EE, Hald SM, et al. Evaluation of tumor-infiltrating lymphocytes using routine H&E slides predicts patient survival in resected non-small cell lung cancer. Hum Pathol (2018) 79:188–98. doi: 10.1016/j.humpath.2018.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Rakaee M, Kilvaer TK, Jamaly S, Berg T, Paulsen EE, Berglund M, et al. Tertiary lymphoid structure score: a promising approach to refine the TNM staging in resected non-small cell lung cancer. Br J Cancer (2021) 124(10):1680–9. doi: 10.1038/s41416-021-01307-y

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Goldstraw P, Chansky K, Crowley J, Rami-Porta R, Asamura H, Eberhardt WE, et al. The IASLC lung cancer staging project: proposals for revision of the TNM stage groupings in the forthcoming (Eighth) edition of the TNM classification for lung cancer. J Thorac Oncol (2016) 11(1):39–51. doi: 10.1016/j.jtho.2015.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, et al. The 2015 world health organization classification of lung tumors: impact of genetic, clinical and radiologic advances since the 2004 classification. J Thorac Oncol (2015) 10(9):1243–60. doi: 10.1097/JTO.0000000000000630

PubMed Abstract | CrossRef Full Text | Google Scholar

29. McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM, et al. REporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer (2005) 93(4):387–91. doi: 10.1038/sj.bjc.6602678

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Al-Saad S, Al-Shibli K, Donnem T, Andersen S, Bremnes RM, Busund LT. Clinical significance of epidermal growth factor receptors in non-small cell lung cancer and a prognostic role for HER2 gene copy number in female patients. J Thorac Oncol (2010) 5(10):1536–43. doi: 10.1097/JTO.0b013e3181ea510a

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Al-Saad S, Al-Shibli K, Donnem T, Persson M, Bremnes RM, Busund LT. The prognostic impact of NF-kappaB p105, vimentin, e-cadherin and Par6 expression in epithelial and stromal compartment in non-small-cell lung cancer. Br J Cancer (2008) 99(9):1476–83. doi: 10.1038/sj.bjc.6604713

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Al-Saad S, Donnem T, Al-Shibli K, Persson M, Bremnes RM, Busund LT. Diverse prognostic roles of akt isoforms, PTEN and PI3K in tumor epithelial cells and stromal compartment in non-small cell lung cancer. Anticancer Res (2009) 29(10):4175–83.

PubMed Abstract | Google Scholar

33. Andersen S, Donnem T, Al-Shibli K, Al-Saad S, Stenvold H, Busund LT, et al. Prognostic impacts of angiopoietins in NSCLC tumor cells and stroma: VEGF-a impact is strongly associated with ang-2. PloS One (2011) 6(5):e19773. doi: 10.1371/journal.pone.0019773

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Andersen S, Donnem T, Al-Saad S, Al-Shibli K, Stenvold H, Busund LT, et al. Correlation and coexpression of HIFs and NOTCH markers in NSCLC. Anticancer Res (2011) 31(5):1603–6.

PubMed Abstract | Google Scholar

35. Andersen S, Donnem T, Stenvold H, Al-Saad S, Al-Shibli K, Busund LT, et al. Overexpression of the HIF hydroxylases PHD1, PHD2, PHD3 and FIH are individually and collectively unfavorable prognosticators for NSCLC survival. PloS One (2011) 6(8):e23847. doi: 10.1371/journal.pone.0023847

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Donnem T, Al-Shibli K, Al-Saad S, Busund LT, Bremnes RM. Prognostic impact of fibroblast growth factor 2 in non-small cell lung cancer: coexpression with VEGFR-3 and PDGF-b predicts poor survival. J Thorac Oncol (2009) 4(5):578–85. doi: 10.1097/JTO.0b013e31819f2e38

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Donnem T, Al-Saad S, Al-Shibli K, Andersen S, Busund LT, Bremnes RM. Prognostic impact of platelet-derived growth factors in non-small cell lung cancer tumor and stromal cells. J Thorac Oncol (2008) 3(9):963–70. doi: 10.1097/JTO.0b013e3181834f52

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Donnem T, Al-Saad S, Al-Shibli K, Busund LT, Bremnes RM. Co-Expression of PDGF-b and VEGFR-3 strongly correlates with lymph node metastasis and poor survival in non-small-cell lung cancer. Ann Oncol (2010) 21(2):223–31. doi: 10.1093/annonc/mdp296

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Donnem T, Andersen S, Al-Shibli K, Al-Saad S, Busund LT, Bremnes RM. Prognostic impact of notch ligands and receptors in nonsmall cell lung cancer: coexpression of notch-1 and vascular endothelial growth factor-a predicts poor survival. Cancer (2010) 116(24):5676–85. doi: 10.1002/cncr.25551

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Donnem T, Eklo K, Berg T, Sorbye SW, Lonvik K, Al-Saad S, et al. Prognostic impact of MiR-155 in non-small cell lung cancer evaluated by in situ hybridization. J Transl Med (2011) 9:6. doi: 10.1186/1479-5876-9-6

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Donnem T, Lonvik K, Eklo K, Berg T, Sorbye SW, Al-Shibli K, et al. Independent and tissue-specific prognostic impact of miR-126 in nonsmall cell lung cancer: coexpression with vascular endothelial growth factor-a predicts poor survival. Cancer (2011) 117(14):3193–200. doi: 10.1002/cncr.25907

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Eilertsen M, Andersen S, Al-Saad S, Kiselev Y, Donnem T, Stenvold H, et al. Monocarboxylate transporters 1-4 in NSCLC: MCT1 is an independent prognostic marker for survival. PloS One (2014) 9(9):e105038. doi: 10.1371/journal.pone.0105038

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Eilertsen M, Andersen S, Al-Saad S, Richardsen E, Stenvold H, Hald SM, et al. Positive prognostic impact of miR-210 in non-small cell lung cancer. Lung Cancer (2014) 83(2):272–8. doi: 10.1016/j.lungcan.2013.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Hald SM, Kiselev Y, Al-Saad S, Richardsen E, Johannessen C, Eilertsen M, et al. Prognostic impact of CXCL16 and CXCR6 in non-small cell lung cancer: combined high CXCL16 expression in tumor stroma and cancer cells yields improved survival. BMC Cancer (2015) 15:441. doi: 10.1186/s12885-015-1446-z

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Kilvaer TK, Khanehkenari MR, Hellevik T, Al-Saad S, Paulsen EE, Bremnes RM, et al. Cancer associated fibroblasts in stage I-IIIA NSCLC: prognostic impact and their correlations with tumor molecular markers. PloS One (2015) 10(8):e0134965. doi: 10.1371/journal.pone.0134965

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Kilvaer TK, Rakaee M, Hellevik T, Ostman A, Strell C, Bremnes RM, et al. Tissue analyses reveal a potential immune-adjuvant function of FAP-1 positive fibroblasts in non-small cell lung cancer. PloS One (2018) 13(2):e0192157. doi: 10.1371/journal.pone.0192157

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Stenvold H, Donnem T, Andersen S, Al-Saad S, Al-Shibli K, Busund LT, et al. Overexpression of matrix metalloproteinase-7 and -9 in NSCLC tumor and stromal cells: correlation with a favorable clinical outcome. Lung Cancer (2012) 75(2):235–41. doi: 10.1016/j.lungcan.2011.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Stenvold H, Donnem T, Andersen S, Al-Saad S, Busund LT, Bremnes RM. Stage and tissue-specific prognostic impact of miR-182 in NSCLC. BMC Cancer (2014) 14:138. doi: 10.1186/1471-2407-14-138

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Stenvold H, Donnem T, Andersen S, Al-Saad S, Valkov A, Pedersen MI, et al. High tumor cell expression of microRNA-21 in node positive non-small cell lung cancer predicts a favorable clinical outcome. BMC Clin Pathol (2014) 14(1):9. doi: 10.1186/1472-6890-14-9

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: open source software for digital pathology image analysis. Sci Rep (2017) 7(1):16878. doi: 10.1038/s41598-017-17204-5

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Sardari Nia P, Colpaert C, Vermeulen P, Weyler J, Pezzella F, Van Schil P, et al. Different growth patterns of non-small cell lung cancer represent distinct biologic subtypes. Ann Thorac Surg (2008) 85(2):395–405. doi: 10.1016/j.athoracsur.2007.08.054

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Latacz E, Caspani E, Barnhill R, Lugassy C, Verhoef C, Grunhagen D, et al. Pathological features of vessel co-option versus sprouting angiogenesis. Angiogenesis (2020) 23(1):43–54. doi: 10.1007/s10456-019-09690-0

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Donnem T, Kilvaer TK, Andersen S, Richardsen E, Paulsen EE, Hald SM, et al. Strategies for clinical implementation of TNM-immunoscore in resected nonsmall-cell lung cancer. Ann Oncol (2016) 27(2):225–32. doi: 10.1093/annonc/mdv560

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Rakaee M, Adib E, Ricciuti B, Sholl LM, Shi W, Alessi JV, et al. Association of machine learning-based assessment of tumor-infiltrating lymphocytes on standard histologic images with outcomes of immunotherapy in patients with NSCLC. JAMA Oncol (2022) 9(1):51–60. doi: 10.1001/jamaoncol.2022.4933

CrossRef Full Text | Google Scholar