- Open Access
Comprehensive analysis of angiogenesis subtype of squamous cell carcinoma
World Journal of Surgical Oncology volume 19, Article number: 275 (2021)
Squamous cell carcinoma (SCC) is a disease with distinct management complexities as it displays a remarkably heterogeneous molecular subtype. However, the landscape of angiogenesis for SCC is not fully investigated.
Method and materials
The angiogenesis-related subtypes of SCC were established by using the ConsensusClusterPlus package based on angiogenesis-related genes and TCGA data. We analyzed the alteration of genes and miRNAs as well as pathways associated with angiogenesis subtypes. Next, the regulation network, the correlation with genomic characteristics, immune microenvironment, and clinical features of the angiogenesis subtypes were further investigated. Finally, the prognostic impact of the angiogenesis-related subtypes for SCC was also analyzed.
A total of 1368 SCC samples were included in this study. Two angiogenesis subtypes were then identified based on the one hundred and sixty-three angiogenesis-related genes with subtype1 (angiogenesis subtype) of 951 SCC patients and subtype2 (non-angiogenesis subtype) of 417 SCC. GSEA revealed that angiogenesis and epithelial-mesenchymal transition, inflammatory response, and hypoxia were enriched in the angiogenesis subtype. Eight of the 15 immune checkpoints (ADORA2A, BTLA, CD276, CYBB, HAVCR2, SIGLEC7, SIGLEC9, and VTCN1) were significantly upregulated while C10orf54 were significantly downregulated in the angiogenesis subtype. The survival analysis revealed that the patients in the angiogenesis subtype have poorer survival outcomes than those in the non-angiogenesis subtype (P = 0.017 for disease-free interval and P = 0.00013 for overall survival).
Our analysis revealed a novel angiogenesis subtype classification in SCC and provides new insights into a hallmark of SCC progression.
Squamous cell carcinoma (SCC) represents the most common human solid tumor and is a major cause of cancer mortality . The occurrence of SCC is closely related to the key changes in genome disturbance, gene mutation, and molecular expression at different stages. Fundamental changes in mesenchymal cells also play an important role in the development of these tumors, which might even be a major determinant of promoting immune escape and chemotherapy drug resistance . PAM50, a classification of breast cancer widely used in gene expression profiling, can divide different clinical outcomes into five subtypes . On the other hand, changes in certain cancer-related genes epidermal growth factor receptor (EGFR) and TP53 pathway present in different SCC types [4, 5]. Also, it is demonstrated that the smoking status is an important factor for the response of immune therapy for SCC like non-small cell lung cancer (NSCLC) . The novel roles of the protein-coding genes were found in SCC regarding angiogenesis and aerobic glycolysis pathway [7,8,9,10]. The recent studies also focus on miRNA  and lncRNAs  involved in the tumorigenesis and progression of SCC. SCCs share common histologic features and have similar molecular patterns, which are different from other cancer types .
Angiogenesis is a change in the balance between pro-angiogenesis and anti-angiogenesis factors . Increased angiogenesis is associated with tumor progression, metastasis, and unfavorable outcome . Numerous studies have shown that solid tumors are “angiogenesis-dependent” . For example, in oral squamous cell carcinoma (OSCC), keratinocytes, and inflammatory cells directly produce a variety of molecules that can induce angiogenesis. Besides, in many different tumors including head and neck squamous cell carcinoma (HNSCC), the increased expression of vascular endothelial growth factor (VEGF) protein contributes to the induction of angiogenesis in tumors [17, 18]. The relationship between angiogenesis and tumor is a recognized factor, while the full biological mechanisms on angiogenesis and SCC are not fully investigated. Although in recent years it has made many pan-cancer classifications based on gene expression, a consensus on angiogenesis molecular for SCC was not yet established. Therefore, it is essential to discover the new molecular subtypes of SCC research.
To determine the association between SCC subtypes and the angiogenesis process, this study established a novel SCC classification based on 474 angiogenesis genes and The Cancer Genome Atlas (TCGA) database. The regulation network, genome mutations, immune characteristics, and prognostic value of subtypes were explored to unveil the potential associations between angiogenesis and SCC.
Materials and methods
Data download and preprocessing
TCGA pan-cancer data were downloaded from the UCSC Genome Browser (https://genome.ucsc.edu), including batch-effect normalized transcription data, clinical data, single-sample gene set enrichment analysis (ssGSEA) score data, drug target data, homologous recombination deficiency (HRD) score and genome-wide DNA damage data, immune signature scores data, and RNA-based stemness scores data. All data processing is described on the official website. The pan-cancer study combined with clinical data included a total of 1368 SCC samples, including 252 cervical squamous cell carcinomas (CESCs), 95 esophageal squamous cell carcinomas (ESCAs), 520 HNSCs, and 501 lung squamous cell carcinomas (LUSCs).
First of all, we obtained 507 angiogenesis-related genes from the AmiGO2 website (http://amigo.geneontology.org/amigo). Combined with gene expression data, 474 angiogenesis genes were finally obtained for analysis. Univariable Cox analysis was used to filter the angiogenesis genes which had the prognostic value for SCC patients (P < 0.05). Based on the prognostic angiogenesis-related genes, ConsensusClusterPlus R-package was used to identify subtypes in SCC tumor samples using 1000 iterations, 80% sample resampling from 2 to 7 clusters (k2 to k7) using kmdist with average linkage algorithm and correlation as the similarity metric.
Gene set enrichment analysis and pathway
To study the changes in gene sets, Gene Set Enrichment Analysis (GSEA) was performed on all genes . We analyzed the correlations between angiogenesis subtypes with cancer hallmark pathways and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway in each tumor sample. GSEA can highlight genes associated with the subtypes through pathway analysis. Gene Ontology (GO) analysis and KEGG pathway enrichment analysis were also performed based on differentially expressed genes (DEGs).
Genomic correlations with angiogenesis subtypes
Aneuploidy and LOH scores and ABSOLUTE purity/ploidy files were obtained from the research by Thorsson et al. . All purity, ploidy, LOH, and copy number variant (CNV) invocation used to create the DNA damage scores were derived by the TCGA Aneuploidy AWG using ABSOLUTE . Moreover, HRD and HRD-loss of heterozygosity (HRD-LOH) scores were obtained from the UCSC genome browser. The copy number burden fraction change and the number of segments represent the base fraction deviating from the baseline multiplicity and the total number of segments in each sample’s copy number profile, respectively. Each fragment was designated as amplification, deletion, or neutral based on its number of copies relative to the circular ploidy of the sample. In addition, we calculated Oncoplot, mutation landscape, and OncogenicPathways based on TCGAmutation and maftools R packages.
Differentially expressed genes and regulations associated with angiogenesis
Mann-Whitney U test was performed to derive the DEGs and miRNAs between subtype1 and subtype2 (FDR < 0.05, absolute logFC > 1). Abnormal vascular networks due to tumor cells can secrete a large number of pro-angiogenesis factors, which are characterized by vascular disease, immaturity, and permeability . To clarify the regulation of angiogenesis subtypes, we performed a computational analysis to identify two “master regulators”: transcription factors (TF) and miRNA. We first downloaded 318 TFs from the Cistrome Bowser (http://cistrome.org/). The correlation was analyzed between 163 angiogenesis genes obtained from survival analysis with TF (person correlation: R2 = 0.4, P < 0.05). The Mann-Whitney U test was utilized to calculate the differential miRNA (FDR = 0.05, log FC > 1). And the target gene of five miRNAs was predicted by TargetScan (http://www.targetscan.org/vert_72/). Each of the predicted miRNA target genes is the result of the correlation test which intersected angiogenesis genes. Based on TransmiR v2.0 database (http://www.cuilab.cn/transmir), the predicted miRNA was then intersected with the TF. Finally, the TF-miRNA-target gene regulatory network was constructed using Cytoscape software (https://cytoscape.org/).
The microenvironment in the angiogenesis subtypes
Tumor-infiltrating leukocytes (TILs) have been shown to be associated with tumor prognosis and treatment response and are a crucial component of the tumor microenvironment . CIBERSORT is a common computational method for the quantification of cellular components by gene expression profiling (GEP) from bulk tissues . Therefore, we upload the gene expression data to the CIBERSORT website, which can be obtained free of charge (https://cibersort.stanford.edu/). We used a leukocyte gene signature matrix, termed LM22 and 1000 permutations. LM22 contains 547 genes and distinguishes 22 kinds of human hematopoietic cell subtypes, including 7 types of T cells, naive and memory B cells, plasma cells, NK cells, and myeloid subpopulations . We rank the gene expression values of a given sample by ssGSEA and then calculated the enrichment score of angiogenesis genes in each sample. In the ssGSEA analysis, we used only genes associated with angiogenesis for our calculations by GSVA, limma, and GSEABase R packages. The Estimating Stromal Cells and Immune Cells in Malignant Tumor Tissues Using Expression Data (ESTIMATE) algorithm was used to calculate stromal cells, immune cells, and estimated scores .
The survival outcome of these patients with different subtypes was calculated by Log-rank test. To investigate the correlation between clinicopathological characteristics and SCC, we further studied the relationship between angiogenesis subtypes and gender, clinical stage (I~IV), tumor status, and histological grade by Chi-square test. Pearson correlation test was performed to evaluate the immune signatures level for each sample. The enrichment levels of 68 immune signatures and angiogenesis were quantified by the heatmap R package. We compared the differences in drug target types between subtype1 and subtype2 and visualized the differences using the “pheatmap” R software package. A two-sided P value less than 0.05 was set as a statistical significance threshold. All the analysis were performed based on R version 3.4.2 (https://www.r-project.org/).
A total of 252 CESCs, 95 ESCAs 520 HNSCs, and 501 LUSCs were included in this study. A total of 90 grade I, 461 grade II, 248 grade III, and 8 grade IV SCC were analyzed (Table 1). Based on univariate Cox analysis for survival analysis, we got one hundred and sixty-three genes with the prognostic impact on SCC (P < 0.05). When K = 2, the module boundaries are clearest and free of crossings. We divided the patients into subtype1 (951 SCC patients) and subtype2 (417 SCC patients) based on the one hundred and sixty-three genes of survival analysis (Fig. 1A).
Genomic correlations with angiogenesis subtypes
As result, the cancer DNA fraction, subclonal genome fraction, and LOH fraction altered were significantly higher in subtype1 (Fig. 1B). Additionally, the tumor purity was significantly lower and the ploidy was significantly higher in subtype1 (Fig. 1C). However, the HRD-LOH presented no significant difference between the two subtypes. HRD in subtype1 was significantly higher than that in subtype2 (Fig. 1D). The more prime LOH was found in subtype1 (Fig. 1E). Moreover, the first 20 highly mutated genes were shown in Oncoplot (Fig. 2A, B). A landscape of numbers of mutations and mutation types was showed in Fig. 2C and D. Ten significantly enriched oncogenic pathways had the greatest impact in subtype1 including RTK-RAS, WNT, Hippo, PI3K, Cell Cycle, MYC, TGF-Beta, TP53, and NRF2 (Fig. 2E). The TGF-Beta and NRF2 pathways had the greatest effect in subtype2 (Fig. 2F).
Differentially expressed genes and network regulation of angiogenesis subtypes
We identified 163 DEGs including 124 upregulated genes and 18 with downregulated genes in the subtype1 (P < 0.05, Fig. 3A). Next, we analyzed the differently expressed miRNAs between the two subtypes. As result, a total of 34 differently expressed miRNAs were found between two subtypes, including 11 downregulated miRNAs and 23 upregulated miRNAs (Fig. 3B). To reveal whether the regulatory relationship between TF and miRNAs has an impact on SCC, we established the regulation network for TF-miRNA-target (Fig. 3C). Finally, most of the drug targets were shown to be also altered by angiogenesis (Fig. 3D).
Enrichment pathways for angiogenesis subtype
According to DEGs, twenty pathways were significantly enriched (P < 0.05) by performing KEGG pathway analysis, which included the PI3K/Akt signaling pathway and ECM receptor interactions (Fig. 4A). Notably, twenty GO terms with DEGs were also found to be enriched, extracellular matrix structural constituent with the most significant P-value, with the maximum difference (Fig. 4B). Figure 4C shows the correlation between different pathways and the angiogenesis subtype. GSEA analysis shows these pathways were strongly associated with angiogenesis in tumor progression (Fig. 4D). Thus, we defined subtype1 as the angiogenesis subtype and subtype2 as the non-angiogenesis subtype. As for the angiogenesis subtype, ten pathways were found to be upregulated including focal adhesion, pathways in cancer, adherens junction, renal cell carcinoma, regulation of actin cytoskeleton, neurotrophin signaling pathway, ECM receptor interaction, TGF-beta signaling pathway, ERBB signaling pathway, and GAP junction. On the other hand, another ten KEGG pathways were upregulated in the non-angiogenesis subtype such as spliceosome, base excision repair DNA replication. Notably, ten hallmarks’ pathways were also found to be enriched in the angiogenesis subtype, including epithelial-mesenchymal transition, inflammatory response, hypoxia, and angiogenesis. As for subtype2, oxidative phosphorylation DNA repair, E2F targets, MCY targets v1, G2M checkpoint, and MYC targets v2 were found to be the enriched hallmarks pathways.
The immune microenvironment in the angiogenesis subtypes
Next, we aimed to explore the microenvironment in these two angiogenesis subtypes. By applying CIBERSORT, we discovered complex associations between 22 different leukocyte subsets and angiogenesis subtypes (Fig. 5A). We also explored the distribution of the immune cells between the angiogenesis subtype and non-angiogenesis subtype (Fig. 5B). Additionally, we also evaluated the immune score by using the ESTIMATE algorithm. As result, the ESTIMATE score immune score and stromal score were significantly higher in the angiogenesis subtype and the tumor purity was significantly lower in the non-angiogenesis subtype (Fig. 5C). Further, we analyzed the expression of 15 immune checkpoints between two angiogenesis subtypes. The results showed that 8 of the 15 immune checkpoints (ADORA2A, BTLA, 276, CYBB, HAVCR2, SIGLEC7, SIGLEC9, and VTCN1) were significantly upregulated while C10orf54 were significantly downregulated in the angiogenesis subtype (Fig. 5D).
The clinical implication of angiogenesis subtypes
To evaluate the prognostic significance of angiogenesis subtype for SCC patients. The KM curves revealed that the patients in the angiogenesis subtype have a lower DFI than that in the non-angiogenesis subtype (P = 0.017, Fig. 6A). Additionally, the patients of the angiogenesis subtype were revealed to have a poor OS outcome than that of the non-angiogenesis subtype (P = 0.00013, Fig. 6B). However, there was no statistical significance was found for the DSS and PFI (Fig. 6C, D).
In this study, we represent the novel insights into the angiogenic features in 1368 samples from 4 SCC types, in which two novel angiogenesis-related subtypes of SCC were established based on TCGA database. The angiogenesis subtype displayed distinct angiogenesis characteristics and had a poorer prognosis, providing new evidence for the progression of SCC.
So far, most studies investigating the subtypes and SCC were mainly focused on the immune microenvironment. A previous pan-cancer study classified the SCCs fell into wound healing and IFN g dominant subtypes based on the immune signature . Another TCGA study has identified six immune subtypes of SCCs including angiogenesis, while the molecular characteristics of angiogenesis were not fully uncovered . Recently, Benjamin et al. developed an endothelial index to estimate the vascular density of 31 solid tumor types and classified the human tumors into 6 vascular microenvironment signatures based on the 24 vascular “hub” genes . While the different histology types of tumors present widely different patterns, our study only performed the analysis for SCCs and aimed to find out similar features for these types of tumors, providing the evidence for the personalized therapy.
Angiogenesis is a complex process that plays a vital role in organ and tissue regeneration, growth and development, and many pathological conditions . Angiogenesis contributes to the development of malignant subtype traits. It is believed that the transition to the angiogenesis subtype is caused by the change in the balance of positive and negative regulators of angiogenesis . Tumors require a blood supply to grow and may be generated by the expression of pro-angiogenesis growth factors. Our study revealed that the angiogenesis subtype was enriched in the epithelial-mesenchymal transition (EMT) pathway. EMT has been closely linked to “stemness” in the development of cancer . The study by Zhang et al. showed that endothelial cells secrete factors which may be attracted to the epithelial tumor cells by blood vessels, allowing it to pass through the EMT connective tissue transfer, and to enhance its potential by imparting tumorigenic tumor cells with a stem cell-like phenotype . Ke et al. found that VEGF-A and Notch signaling pathways were activated by LncRNA NEAT1, leading to promote EMT and repress apoptosis in oral squamous cell carcinoma . Moreover, the results and also the gene set in inflammatory were enriched in the angiogenesis subtype. Inflammatory responses increase the risk of developing many types of cancer. The hallmark of cancer-related inflammation includes inflammatory cells and tissue expression of inflammatory mediators (e.g., chemokines, cytokines, and prostaglandins), and tissue remodeling, and angiogenesis, chronic inflammatory diseases, and similar reactions and tissue repair . Furthermore, hypoxia is a key factor in the tumor angiogenesis process. Studies have shown that tissue is upregulated by HIF-1a promoter steady-state reaction, to adapt to hypoxia, which may further facilitate tumor growth and tumor angiogenesis . In summary, the angiogenesis subtype is linked to the characteristics of angiogenesis closely and promotes the development of cancer.
Also, eight of the 15 immune checkpoint genes were upregulated, whereas C10orf54 was significantly downregulated in the angiogenesis subtype. Cancer cells have various mechanisms to evade the local immune attack, including upregulation of immune checkpoint proteins . Immune checkpoint proteins are activated by ligand-receptor, resulting in a dynamic balance between stimulatory and non-stimulatory, and inhibitory signals that regulate the immune response . This phase of equilibrium is ensured primarily by the PD-1/PD-L pathway, which inhibits the activation and proliferation of T-cell as well as cytokine production. And the cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) pathway, which induces cell cycle arrest and apoptosis in Tregs and T cells activation. However, once tumor antigens expressed by tumor cells can bypass the immune system, the balance is altered [34, 36] and the final result showed the progression of cancerous tumors in clinical practice . The appropriate subtype classification may provide evidence for the selection of appropriate therapy.
In summary, we established a novel subtype classification related to angiogenesis for SCC with distinct molecular characteristics and clinical implications. This work provides new sights into the progression of SCC, contributing to the future design of personalized therapy.
Availability of data and materials
Data and materials of this work are available from the corresponding author on reasonable request.
Squamous cell carcinoma
Oral squamous cell carcinoma
Head and neck squamous cell carcinoma
Single-sample gene set enrichment analysis
Cervical squamous cell carcinoma
Esophageal squamous cell carcinomas
Lung squamous cell carcinoma
Homologous recombination deficiency
Gene set enrichment analysis
Kyoto encyclopedia of genes and genomes
HRD-loss of heterozygosity
False discovery rate
The Cancer Genome Atlas
Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data
Endothelial cells secrete factor
Vascular endothelial growth factor
Transforming growth factor beta
T-lymphocyte-associated antigen 4
Gene expression profiling
Dotto GP, Rustgi AK. Squamous cell cancers: a unified perspective on biology and genetics. Cancer Cell. 2016;29(5):622–37. https://doi.org/10.1016/j.ccell.2016.04.004.
Dotto GP. Multifocal epithelial tumors and field cancerization: stroma as a primary determinant. J Clin Invest. 2014;124(4):1446–53. https://doi.org/10.1172/JCI72589.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7. https://doi.org/10.1200/JCO.2008.18.1370.
Martínez E, Yoshihara K, Kim H, Mills GM, Treviño V, Verhaak RG. Comparison of gene expression patterns across 12 tumor types identifies a cancer supercluster characterized by TP53 mutations and cell cycle defects. Oncogene. 2015;34(21):2732–40. https://doi.org/10.1038/onc.2014.216.
Pandey M, Kannepali KK, Dixit R, Kumar M. Effect of neoadjuvant chemotherapy and its correlation with HPV status, EGFR, Her-2-neu, and GADD45 expression in oral squamous cell carcinoma. World J Surg Oncol. 2018;16(1):20. https://doi.org/10.1186/s12957-018-1308-7.
Mo J, Hu X, Gu L, Chen B, Khadaroo PA, Shen Z, et al. Smokers or non-smokers: who benefits more from immune checkpoint inhibitors in treatment of malignancies? An up-to-date meta-analysis. World J Surg Oncol. 2020;18(1):15. https://doi.org/10.1186/s12957-020-1792-4.
Zhang R, Zhang TT, Zhai GQ, Guo XY, Qin Y, Gan TQ, et al. Evaluation of the HOXA11 level in patients with lung squamous cancer and insights into potential molecular pathways via bioinformatics analysis. World J Surg Oncol. 2018;16(1):109. https://doi.org/10.1186/s12957-018-1375-9.
Mao L, Wu X, Gong Z, Yu M, Huang Z. PDIA6 contributes to aerobic glycolysis and cancer progression in oral squamous cell carcinoma. World J Surg Oncol. 2021;19(1):88. https://doi.org/10.1186/s12957-021-02190-w.
Zheng Z, Tian R, Wang P. Roles of KAI1 and nm23 in lymphangiogenesis and lymph metastasis of laryngeal squamous cell carcinoma. World J Surg Oncol. 2017;15(1):211. https://doi.org/10.1186/s12957-017-1279-0.
Tang J, Yang Z, Wang Z, Li Z, Li H, Yin J, et al. Foxp3 is correlated with VEGF-C expression and lymphangiogenesis in cervical cancer. World J Surg Oncol. 2017;15(1):173. https://doi.org/10.1186/s12957-017-1221-5.
Luo J, Shi K, Yin SY, Tang RX, Chen WJ, Huang LZ, et al. Clinical value of miR-182-5p in lung squamous cell carcinoma: a study combining data from TCGA, GEO, and RT-qPCR validation. World J Surg Oncol. 2018;16(1):76. https://doi.org/10.1186/s12957-018-1378-6.
He K, Zhu ZB, Shu R, Hong A. LncRNA NEAT1 mediates progression of oral squamous cell carcinoma via VEGF-A and Notch signaling pathway. World J Surg Oncol. 2020;18(1):261. https://doi.org/10.1186/s12957-020-02028-x.
Hoadley KA, Yau C, Hinoue T, Wolf DM, Lazar AJ, Drill E, et al. Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer. Cell. 2018;173(2):291–304.e296.
Marla V, Hegde V, Shrestha A. Relationship of angiogenesis and oral squamous cell carcinoma. Kathmandu Univ Med J (KUMJ). 2015;13(50):178–85. https://doi.org/10.3126/kumj.v13i2.16796.
Kabiraj A, Jaiswal R, Singh A, Gupta J, Singh A, Samadi FM. Immunohistochemical evaluation of tumor angiogenesis and the role of mast cells in oral squamous cell carcinoma. J Cancer Res Ther. 2018;14(3):495–502. https://doi.org/10.4103/0973-1482.163693.
Folkman J. Fundamental concepts of the angiogenic process. Curr Mol Med. 2003;3(7):643–51. https://doi.org/10.2174/1566524033479465.
Ascani G, Balercia P, Messi M, Lupi L, Goteri G, Filosa A, et al. Angiogenesis in oral squamous cell carcinoma. Acta Otorhinolaryngol Ital. 2005;25(1):13–7.
Lingen MW. Angiogenesis in the development of head and neck cancer and its inhibition by chemopreventive agents. Crit Rev Oral Biol Med. 1999;10(2):153–64. https://doi.org/10.1177/10454411990100020301.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of cancer. Immunity. 2018;48(4):812–830.e814.
Carter SL, Cibulskis K, Helman E, McKenna A, Shen H, Zack T, et al. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012;30(5):413–21. https://doi.org/10.1038/nbt.2203.
Viallard C, Larrivée B. Tumor angiogenesis and vascular normalization: alternative therapeutic targets. Angiogenesis. 2017;20(4):409–26. https://doi.org/10.1007/s10456-017-9562-9.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7. https://doi.org/10.1038/nmeth.3337.
Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol. 2018;1711:243–59.
Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4(1):2612. https://doi.org/10.1038/ncomms3612.
Li B, Cui Y, Nambiar DK, Sunwoo JB, Li R. The immune subtypes and landscape of squamous cell carcinoma. Clin Cancer Res. 2019;25(12):3528–37. https://doi.org/10.1158/1078-0432.CCR-18-4085.
Kahn BM, Lucas A, Alur RG, Wengyn MD, Schwartz GW, Li J, Sun K, Maurer HC, Olive KP, Faryabi RB, Stanger BZ. The vascular landscape of human cancer. J Clin Invest. 2021:131.
Nowak-Sliwinska P, Alitalo K, Allen E, Anisimov A, Aplin AC, Auerbach R, et al. Consensus guidelines for the use and interpretation of angiogenesis assays. Angiogenesis. 2018;21(3):425–532. https://doi.org/10.1007/s10456-018-9613-x.
Hawighorst T. Angiogenesis, lymphangiogenesis, and tumor progression. Zentralbl Gynakol. 2002;124(11):497–505. https://doi.org/10.1055/s-2002-39653.
Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2014;15(3):178–96. https://doi.org/10.1038/nrm3758.
Zhang Z, Dong Z, Lauxen IS, Filho MS, Nör JE. Endothelial cell-secreted EGF induces epithelial to mesenchymal transition and endows head and neck cancer cells with stem-like phenotype. Cancer Res. 2014;74(10):2869–81. https://doi.org/10.1158/0008-5472.CAN-13-2032.
Mantovani A, Allavena P, Sica A, Balkwill F. Cancer-related inflammation. Nature. 2008;454(7203):436–44. https://doi.org/10.1038/nature07205.
An X, Xu G, Yang L, Wang Y, Li Y, McHepange UO, et al. Expression of hypoxia-inducible factor-1α, vascular endothelial growth factor and prolyl hydroxylase domain protein 2 in cutaneous squamous cell carcinoma and precursor lesions and their relationship with histological stages and clinical features. J Dermatol. 2014;41(1):76–83. https://doi.org/10.1111/1346-8138.12314.
De Felice F, Marchetti C, Palaia I, Ostuni R, Muzii L, Tombolini V, et al. Immune check-point in cervical cancer. Crit Rev Oncol Hematol. 2018;129:40–3. https://doi.org/10.1016/j.critrevonc.2018.06.006.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64. https://doi.org/10.1038/nrc3239.
De Felice F, Marchetti C, Palaia I, Musio D, Muzii L, Tombolini V, et al. Immunotherapy of ovarian cancer: the role of checkpoint inhibitors. J Immunol Res. 2015;2015:191832.
La-Beck NM, Jean GW, Huynh C, Alzghari SK, Lowe DB. Immune Checkpoint Inhibitors: New Insights and Current Place in Cancer Therapy. Pharmacotherapy. 2015;35:963–76.
This study was supported by the National Natural Science Foundation of China (NSFC81660488) and the Guangxi Natural Science Foundation under grant 2017GXNSFAA198123.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Qin, F., Lin, S., Deng, K. et al. Comprehensive analysis of angiogenesis subtype of squamous cell carcinoma. World J Surg Onc 19, 275 (2021). https://doi.org/10.1186/s12957-021-02367-3
- Squamous cell carcinoma
- Comprehensive analysis