Identified IGSF9 association with prognosis and hypoxia in nasopharyngeal carcinoma by bioinformatics analysis
Cancer Cell International volume 20, Article number: 498 (2020)
Despite improvements in nasopharyngeal carcinoma (NPC) treatment, patients with recurrence and metastasis still have a poor prognosis. Thus, the identification of novel biomarkers is urgently needed to predict outcomes and tailor treatment for NPC.
Four data sets were downloaded from Gene Expression Omnibus, and one data set GSE68799 of which was applied to filtrate key modules and hub genes by construction of a co-expression network. Other data sets (GSE12452 and GSE53819) were used to verify hub genes. The data set GSE102349 was devoted to identify prognostic hub genes by survival analysis. To explored whether prognostic hub genes are related to hypoxia signatures in NPC, correlation analysis was carried out, and followed by functional verification experiments of those genes in vitro.
By co-expression network analysis, blue module was regarded as a key module in the benign and malignant group, and IGSF9 of the blue module was identified as a prognostic hub gene. Moreover, IGSF9 is expected to be a innovative hypoxia-related gene in NPC based on the strong associativity between expression of IGSF9 and hypoxia scores of three signatures (99-gene, 26-gene and 15-gene). Further functional studies verified that down-regulated expression of IGSF9 could reduce the proliferation, migration and invasion ability of NPC cells, and hypoxia could induce the expression of IGSF9.
IGSF9 was identified to be relevant to prognosis and involved in hypoxia in NPC. IGSF9 might serve as one novel prognostic indicator of NPC in the future.
Nasopharyngeal carcinoma (NPC) can be found most frequently in southern China and southeast Asian countries. According to the global cancer statistics of 2018, it was estimated that there were 129,079 new NPC cases and 72,987 deaths that year . Despite improved treatment for NPC patients, especially concurrent chemoradiation therapy combined with adjuvant chemotherapy, those patients with recurrence and metastasis still have worse overall survival . Tumor-node-metastasis staging system can provide an effective prognostic prediction for NPC patients . However, this factor fails to be accurate enough for evaluating the survival outcomes.
It is well documented that the number of biomarkers might predict the prognosis of NPC, such as the EBV DNA levels, lactate dehydrogenase (LDH), SOX2, NGX6, FMNL3 and ANXA2 [4,5,6,7,8,9]. However, the molecular mechanism of NPC has not been completely established. Therefore, the identification of novel biomarkers that reflect tumors heterogeneity is urgently needed to predict outcomes and tailor treatment for NPC patients.
It is noted that bioinformatics analysis is commonly used to identify candidate biomarkers from high-throughput data. Weighted gene co-expression network analysis (WGCNA) , which is mainly applied to filtrate hub genes from large-scale data sets, has been successfully applied in several human tumors [11,12,13,14]. WGCNA had been used to analyze the NPC data sets and identify some candidate diagnostic genes [15, 16], which have limited application in the predication of prognosis in NPC patient.
In this study, we performed WGCNA to identify hub genes and validate these hub genes using three data sets (GSE68799, GSE12452 and GSE53819). Another one data set GSE102349 that contains follow-up information was used for univariate Cox regression analysis to screen hub genes with prognostic values. Numerous reports have testified that tumor hypoxia leads to poor prognosis [17,18,19]. We hypothesized that one prognostic hub gene (IGSF9) might be involved in hypoxia in NPC. Hence, Pearson’s correlation analysis was conducted and the result showed that expression of IGSF9 was related to the scores of three hypoxia signatures. Furthermore, functional studies confirmed that the elevated expression of IGSF9 induced by hypoxia increased NPC cells proliferation, invasion and migration in vitro.
Materials and methods
Data download and pre-processing
Four NPC micro-array data sets (Additional file 1: Table S1) were downloaded from the well-known gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). Since the numerical distribution of raw data value of four data sets is too wide to be used for further analysis, we formulated each raw data value in log2 (X + 1), where X is the raw data value, followed by quantiles of them. After normalization, expression values of each gene of data set GSE68799 were ranked from large to small by variance analysis report (VAR) value. The expression profile of data set GSE68799 which contains the top 6000 genes was used for WGCNA.
Co-expression network construction and identification of key modules and hub genes
The processed matrix of the data set GSE68799 was exploited to construct an unsigned gene co-expression network by WGCNA. Genes were clustered by using the hierarchical clustering method according to the topological overlap measure (TOM). In accordance with the criterion of hybrid dynamic tree, each gene module should contain no less than 50 genes. In line with the dissimilarity of module eigengenes (MEs), the height 0.25 was set as a cut off line to combine some modules. In the process of generating the network, 1000 genes were randomly selected every time. Then, modules with the strongest correlation between MEs in each module and clinical traits were identified as key modules. In a similar way, genes with higher module membership (MM) and higher gene significance (GS) were considered as hub genes, filtered criteria was set as MM > 0.8 and GS >|0.7| .
Validation of hub genes
Differentially expressed genes (DEGs) between NPC group and non-NPC group of three data sets (GSE68799, GSE12452 and GSE53819) were screened by the “limma” R package. Cut-off value was determined by adjusted P < 0.05. Validated hub genes were obtained by overlapping hub genes and common DEGs. Comparison analysis for expression of validated hub genes was conducted in benign and malignant groups of three data sets (GSE68799, GSE12452 and GSE53819), and the same analysis for expression of validated hub genes was performed in different clinical stages of the data set GSE102349.
Identification and functional analysis of prognostic hub genes
The profile contains 23,706 genes from the data set GSE102349 was exploited for univariate Cox regression analysis by the “survival” R package. Cut-off value was set as hazard ratio (HR) >|1|, and adjusted P < 0.05. If one validated hub gene was significantly associated with progression-free survival (PFS), it could be identified as a prognostic hub gene. Functional analysis of them was proceeded with gene set enrichment analysis (GSEA) (https://www.gsea-msigdb.org/gsea/index.jsp) . This analysis was performed on the version of MSigDB.v7.0, and each analysis would conduct on 1000 times gene set substitution. Significant pathways were identified by follow criterion: absolute value of normalized enrichment score (NES) > 1, P < 0.05.
Correlation of hypoxia signatures and prognostic hub genes
To explored whether prognostic hub genes are related to hypoxia signatures in NPC, we focused on three previously reported hypoxia signatures (99-gene, 26-gene and 15-gene) [21,22,23], which have been regarded as effective prognostic tools in head and neck cancer. The data set GSE102349 was used to summarize the scores of hypoxia signatures. For each gene of hypoxia signatures, each sample in the data set was given a primary score of either + 1 or − 1 as the following criteria : If the expression of that gene of one sample is greater than the median value of that gene of all samples in the data set, this sample would be marked as + 1 score. If not, the sample would be marked as − 1 score. For each sample, the summed value was calculated from the + 1′s and − 1′s of total genes of hypoxia signatures, which was the scores of hypoxia signatures for that sample in data set GSE102349. After that, the scores of hypoxia signatures were exploited for survival analysis and correlation analysis between the scores and the expression of prognostic hub genes and parallel analysis was performed in P4HA1 which was one of the common genes from three signatures.
Cell culture and siRNA transfection
CNE-2 and 5-8F cells were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum (FBS, Gibco BRL, NY, USA) and incubated in a humidified atmosphere at 37 °C with 5% CO2. For hypoxia exposure, they were incubated in a anoxic chamber at 37 °C with 5% CO2 and 1% O2 for 24 h. Two siRNA oligonucleotides (siIGSF9#1 siIGSF9#2) and negative control were purchased from Sangon Biotech (Shanghai, China) and LipofectamineTM 2000 transfection reagent were from Invitrogen (Camarillo, USA).
Quantitative real time polymerase chain reaction (qRT-PCR)
Total RNA was extracted from cells using TRIzol reagent (Invitrogen) based on the product specification. qRT-PCR was executed using SYBR green qPCR Master Mix (Takara Bio, Otsu, Japan). Comparative CT method was used for quantitative analysis and GAPDH was served as the internal reference. Primer sequences as follow: IGSF9 sense 5′-TCTGTGGACGAGAACTATGAGTG-3′, antisense 5′-AGCCGAGCCCTAGTAGCAT-3′; P4HA1 sense 5′-GGGACAAGCCTCGTATT-3′, antisense 5′-TTGGGTTTGAAATGGTG-3′; GAPDH sense 5′-GACTCATGACCACAGTCCATGC-3′, antisense 5′-AGAGGCAGGGATGATGTTCTG-3′.
Proteins from cells were lysed using RIPA lysate (Beyotime, Shanghai, China), and quantified using Pierce BCA Protein Assay Kit (Thermo Fisher Scientific, MA, USA). Protein lysates were separated by 10% sodium dodecyl sulfate containing polyacrylamide gels and transferred onto a polyvinylidene difluoride membrane (Millipore, MA, USA) by electrophoresis. GAPDH antibody (1:5000, Abcam, MA, USA), IGSF9 antibody (1:1000, Ybscience, Shanghai, China) and P4HA1 antibody (1:1000, Invitrogen) were incubated at 4℃ overnight. Horseradish Both peroxidase-labeled anti-rabbit (1:8000, Santa Cruz, CA, USA) and anti-mouse (1:8000, Santa Cruz) were used at 4 °C for 2 h. The bands on the membrane were visualized by enhanced chemiluminescence (Thermo Fisher Scientific).
Cell proliferation, migration and invasion assays
Cell proliferation was estimate using cell counting kit-8 (CCK8) assay. 1000 cells were seeded in per well of 96-well plates. The 450 nm absorbance of each well was detected using CCK8 reagent (Dojindo, Kumamoto, Japan) incubation for 0, 24, 48, 72 and 96 h. Cell migration was evaluated using wound-healing assays. An artificial scratch was conducted on fused monolayer of adherent cells. 48 h later, photographs were got from an inverted microscope (Olympus,Tokyo, Japan). For cells invasion assay, cells were seeded on the matrigel (Corning, NY, USA) of the upper layer of transwell chamber without serum medium. Then, FBS was added into the lower layer of the transwell chamber. After 48 h, non-invasive cells and matrigel were softly removed, while invasive cells that pass through matrigel into the lower layer were colored with crystal violet, counted and imaged.
Data were expressed as mean ± standard deviation (SD) and the difference between two groups was estimated using two-tailed Students’ t-test. Survival difference was evaluated using Log-rank test. Pearson’s correlation analysis was used to calculate the correlation between two continuous variables. P < 0.05 was defined to be statistically significant.
Identification of two key modules and 79 hub genes
Expression value of all genes of four NPC data sets were normalized by log2 (X + 1) and quantiles (Additional file 1: Fig. S1). The data profile of the top 6,000 genes collated by the VAR values from the data set GSE68799 was used for WGCNA, and 46 samples of this data set could be divided into two clusters (Fig. 1a). Four was recognized as a soft-thresholding power value to create unsigned co-expression network (Fig. 1b, c), and the network contains 22 independent modules (Fig. 1d, f), which were split into two clusters (Fig. 1e). The blue module contains 869 genes, while the dark red module contains 56 genes (Additional file 1: Table S2). These two modules were identified as key modules (Fig. 2a-c). 74 hub genes and five hub genes were found in the blue and dark red module respectively (Additional file 1: Fig. S2).
Validation of 79 hub genes
To verify the 79 hub genes, we screened DEGs between NPC samples and non-NPC samples of three data sets (GSE68799, GSE12452 and GSE53819). 220 up-regulated common DEGs (Fig. 3a) and 250 down- regulated common DEGs (Fig. 3c) were found. Seven of 74 hub genes of the blue module were significantly up-regulated (Fig. 3b) and (Additional file 1: Fig. S3a) in the NPC patients of three data sets, while none of hub genes of the dark red module was significantly down-regulated (Fig. 3d). In addition, expression of four genes (IGSF9, AGRN, STAP2 and PKP4) was significantly up-regulated in advanced NPC patients of the data set GSE102349 (Additional file 1: Fig. S3b), but expression of other three genes (TRIP6, EPCAM and IRF6) was not.
Identification and functional analysis of a prognostic hub gene (IGSF9)
Univariate Cox regression analysis was carried out in the data set GSE102349 using survival package in R. IGSF9 included in seven validated hub genes was identified as a prognostic hub gene (Fig. 4a, b), and an up regulation of IGSF9 could lead to poor PFS (Fig. 4c). Moreover, functional analysis results showed that IGSF9 may promote metastasis of NPC cells through Akt signaling pathway, including XU AKT1 TARGETS 6HR, LEE METASTASIS AND RNA PROCESSING UP and MTOR UP. N4.V1 UP (Fig. 4d, e).
Correlation of the scores of hypoxia signatures and the expression of IGSF9
The scores of three hypoxia signatures (99-gene, 26-gene and 15-gene) were calculated for survival and correlation analysis in the data set GSE102349. Patients with high scores of these hypoxia signatures had unfavorable PFS (Fig. 5a), and a parallel result was found in P4HA1 included in four common genes (Fig. 5d, e); in other words, an over expression of P4HA1 was associated with a poor PFS. Moreover, patients with high scores of these hypoxia signatures owned an over expression of IGSF9 (Fig. 5b, c), which could also be seen in patients with high expression of P4HA1 (Fig. 5f, g).
In vitro functional analysis of IGSF9
To investigate the effect of IGSF9 on the aggressive phenotype, we transfected two different IGSF9 siRNAs (siIGSF9#1 and siIGSF9#2) into CNE-2 and 5-8F cells to down-regulate IGSF9 (Fig. 6a) and performed CCK8, wound-healing and invasion assays to evaluate whether regulation of IGSF9 expression affects the proliferation, invasion and migration of NPC cells. The results showed a significant inhibition of proliferation, invasion and migration in the cells transfected with IGSF9 siRNAs compared with the control cells (Fig. 6b–d). These data indicated that IGSF9 could promote NPC cell malignant biological behaviors that might be associated with poor prognosis. Next, we tried to verify the interaction between the expression of IGSF9 and hypoxia. qRT-PCR and Western blotting were used to measure the expression of IGSF9 and P4HA1 in CNE-2 and 5-8F cells exposed to hypoxia. We observed that expression levels of IGSF9 and P4HA1 were significantly increased in NPC cells exposed to anoxic conditions (Fig. 6e, f). There was an evidence that expression of IGSF9 might be induced by hypoxia, which was consistent with the result from our bioinformatics analysis.
To identify prognostic hub genes of high throughput data sets, WGCNA and univariate Cox regression analysis are very powerful tools. There are two other studies [15, 16] aiming at identifying candidate genes by WGCNA. Besides this way, our research added univariate Cox regression analysis via the mRNA-seq data set GSE102349 that contains follow-up information. By these two tools, IGSF9 was identified as a prognostic hub gene in NPC.
IGSF9 was highly expressed in NPC patients, and it was an appropriate prognostic gene. IGSF9 belongs to the immunoglobulin superfamily , which plays a key role in inhibiting synaptic development by regulating calmodulin-like activity . Calmodulin is involved in tumor metastasis, and our GSEA result also shows that IGSF9 may cause metastasis of NPC cells through Akt signaling pathway. Meanwhile, we verified that IGSF9 promoted the proliferation, migration and invasion of NPC cells in vitro. To date, there is one published report  on over-expression of IGSF9 associated with poor PFS in endometrial cancer tissue. Similar result was found in NPC tissue in our research. Thus, IGSF9 is likely to be a prognostic gene which would promote the invasion and metastasis of NPC cells.
Our study revealed that three hypoxia signatures (99-gene, 26-gene and 15-gene) have prognostic value in NPC patients using the data set GSE102349, though previous evidences have proved that those three hypoxia signatures can successfully discriminate the prognosis of patients with head and neck cancer [21,22,23, 28]. Interestingly, one research suggested that those three hypoxia signatures might be replaced by P4HA1  and another study also confirmed that P4HA1 was a significantly prognostic factor for oral squamous cell carcinoma patients . We also got a similar result in NPC patients. It is clear that our research provides valuable support for consideration of three hypoxia signatures and one single-gene (P4HA1) hypoxia signature as new prognostic tools like IGSF9 in NPC patients. Moreover, we found the expression of IGSFP was positively associated with both the scores of three hypoxia signatures and the expression of P4HA1 utilizing data set GSE102349, and the expression of IGSFP was significantly increased in NPC cells exposed to anoxic microenvironment. These data revealed that the expression of IGSF9 in NPC cells might be affected by hypoxia.
Taken together, we identified IGSF9 as a prognostic hub gene by WGCNA and univariate Cox regression analysis and found the expression of IGSF9 might be induced by hypoxia. We then performed cytology experiments and verified the role of IGSF9 in NPC cells. However, the investigation of deeper mechanism and in vivo experiment of IGSF9 needs to be supplemented further.
Our study suggests that the elevated expression of IGSF9 will be result in a poor PFS and involved in hypoxia in NPC, implying that IGSF9 might be expected to be a new prognostic indicator of NPC in the future.
Availability of data and materials
Four data sets (GSE68799, GSE12452, GSE53819 and GSE102349) used in this article were derived from in GEO database.
Weighted gene co-expression network analysis
Gene expression omnibus
Variance analysis report
Topological overlap measure
Differentially expressed genes
Gene set enrichment analysis
Normalized enrichment score
Fetal bovine serum
Quantitative real time polymerase chain reaction
Cell counting kit 8
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.
Chua MLK, Wee JTS, Hui EP, Chan ATC. Nasopharyngeal carcinoma. Lancet. 2016;387(10022):1012–24. https://doi.org/10.1016/S0140-6736(15)00055-0.
Tang LL, Chen YP, Mao YP, Wang ZX, Guo R, Chen L, et al. Validation of the 8th Edition of the UICC/AJCC staging system for nasopharyngeal carcinoma from endemic areas in the intensity-modulated radiotherapy era. J Natl Compr Canc Netw. 2017;15(7):913–9. https://doi.org/10.6004/jnccn.2017.0121.
Qu H, Huang Y, Zhao S, Zhou Y, Lv W. Prognostic value of Epstein-Barr virus DNA level for nasopharyngeal carcinoma: a meta-analysis of 8128 cases. Eur Arch Otorhinolaryngol. 2020;277(1):9–18. https://doi.org/10.1007/s00405-019-05699-9.
Zhang M, Wei S, Su L, Lv W, Hong J. Prognostic significance of pretreated serum lactate dehydrogenase level in nasopharyngeal carcinoma among Chinese population: a meta-analysis. Medicine (Baltimore). 2016;95(35):e4494. https://doi.org/10.1097/MD.0000000000004494.
Tang J, Zhong G, Wu J, Chen H, Jia Y. SOX2 recruits KLF4 to regulate nasopharyngeal carcinoma proliferation via PI3K/AKT signaling. Oncogenesis. 2018;7(8):61. https://doi.org/10.1038/s41389-018-0074-2.
Wang YT, Sun XY, Wang H, Huang TT, Wang Y, Yang MX. Effects of NGX6 expression on proliferation and invasion of nasopharyngeal carcinoma cells and survival of patients. Eur Rev Med Pharmacol Sci. 2017;21(23):5378–85. https://doi.org/10.26355/eurrev_201712_13923.
Wu Y, Shen Z, Wang K, Ha Y, Lei H, Jia Y, et al. High FMNL3 expression promotes nasopharyngeal carcinoma cell metastasis: role in TGF-beta1-induced epithelia-to-mesenchymal transition. Sci Rep. 2017;7:42507. https://doi.org/10.1038/srep42507.
Chen CY, Lin YS, Chen CH, Chen YJ. Annexin A2-mediated cancer progression and therapeutic resistance in nasopharyngeal carcinoma. J Biomed Sci. 2018;25(1):30. https://doi.org/10.1186/s12929-018-0430-8.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 2008;9:559. https://doi.org/10.1186/1471-2105-9-559.
Li B, Pu K, Wu X. Identifying novel biomarkers in hepatocellular carcinoma by weighted gene co-expression network analysis. J Cell Biochem. 2019. https://doi.org/10.1002/jcb.28420.
Qiu J, Du Z, Wang Y, Zhou Y, Zhang Y, Xie Y, et al. Weighted gene co-expression network analysis reveals modules and hub genes associated with the development of breast cancer. Medicine. 2019;98(6):e14345. https://doi.org/10.1097/MD.0000000000014345.
Zhang X, Feng H, Li Z, Li D, Liu S, Huang H, et al. Application of weighted gene co-expression network analysis to identify key modules and hub genes in oral squamous cell carcinoma tumorigenesis. Onco Targets Ther. 2018;11:6001–21. https://doi.org/10.2147/OTT.S171791.
Tang J, Kong D, Cui Q, Wang K, Zhang D, Gong Y, et al. Prognostic genes of breast cancer identified by gene co-expression network analysis. Front Oncol. 2018;8:374. https://doi.org/10.3389/fonc.2018.00374.
Xue K, Cao J, Wang Y, Zhao X, Yu D, Jin C, et al. Identification of potential therapeutic gene markers in nasopharyngeal carcinoma based on bioinformatics analysis. Clin Transl Sci. 2019. https://doi.org/10.1111/cts.12690.
Ge Y, He Z, Xiang Y, Wang D, Yang Y, Qiu J, et al. The identification of key genes in nasopharyngeal carcinoma by bioinformatics analysis of high-throughput data. Mol Biol Rep. 2019;46(3):2829–40. https://doi.org/10.1007/s11033-019-04729-3.
Harris AL. Hypoxia–a key regulatory factor in tumour growth. Nat Rev Cancer. 2002;2(1):38–47. https://doi.org/10.1038/nrc704.
Liu Y, Zuo D, Huang C, Zhao M, Hou Y, Chang J, et al. Relationship between hypoxia and carcinoembryonic antigen and upregulated carcinoembryonic antigen is associated with poor prognosis in breast cancer patients. Clin Lab. 2019;65:12. https://doi.org/10.7754/Clin.Lab.2019.190337.
Nordsmark M, Overgaard J. Tumor hypoxia is independent of hemoglobin and prognostic for loco-regional tumor control after primary radiotherapy in advanced head and neck cancer. Acta Oncol. 2004;43(4):396–403. https://doi.org/10.1080/02841860410026189.
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 USA. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Winter SC, Buffa FM, Silva P, Miller C, Valentine HR, Turley H, et al. Relation of a hypoxia metagene derived from head and neck cancer to prognosis of multiple cancers. Cancer Res. 2007;67(7):3441–9. https://doi.org/10.1158/0008-5472.CAN-06-3322.
Eustace A, Mani N, Span PN, Irlam JJ, Taylor J, Betts GN, et al. A 26-gene hypoxia signature predicts benefit from hypoxia-modifying therapy in laryngeal cancer but not bladder cancer. Clin Cancer Res. 2013;19(17):4879–88. https://doi.org/10.1158/1078-0432.CCR-13-0542.
Toustrup K, Sorensen BS, Nordsmark M, Busk M, Wiuf C, Alsner J, et al. Development of a hypoxia gene expression classifier with predictive impact for hypoxic modification of radiotherapy in head and neck cancer. Cancer Res. 2011;71(17):5923–31. https://doi.org/10.1158/0008-5472.CAN-11-1182.
van den Beucken T, Koch E, Chu K, Rupaimoole R, Prickaerts P, Adriaens M, et al. Hypoxia promotes stem cell phenotypes and poor prognosis through epigenetic regulation of DICER. Nat Commun. 2014;5:5203. https://doi.org/10.1038/ncomms6203.
Doudney K, Murdoch JN, Braybrook C, Paternotte C, Bentley L, Copp AJ, et al. Cloning and characterization of Igsf9 in mouse and human: a new member of the immunoglobulin superfamily expressed in the developing nervous system. Genomics. 2002;79(5):663–70. https://doi.org/10.1006/geno.2002.6757.
Mishra A, Traut MH, Becker L, Klopstock T, Stein V, Klein R. Genetic evidence for the adhesion protein IgSF9/Dasm1 to regulate inhibitory synapse development independent of its intracellular domain. J Neurosci. 2014;34(12):4187–99. https://doi.org/10.1523/JNEUROSCI.3671-13.2014.
Shi Z, Li C, Tarwater L, Li J, Li Y, Kaliney W, et al. RNA-seq reveals the overexpression of IGSF9 in endometrial cancer. J Oncol. 2018;2018:2439527. https://doi.org/10.1155/2018/2439527.
Tawk B, Schwager C, Deffaa O, Dyckhoff G, Warta R, Linge A, et al. Comparative analysis of transcriptomics based hypoxia signatures in head- and neck squamous cell carcinoma. Radiother Oncol. 2016;118(2):350–8. https://doi.org/10.1016/j.radonc.2015.11.027.
Kappler M, Kotrba J, Kaune T, Bache M, Rot S, Bethmann D, et al. P4HA1: a single-gene surrogate of hypoxia signatures in oral squamous cell carcinoma patients. Clin Transl Radiat Oncol. 2017;5:6–11. https://doi.org/10.1016/j.ctro.2017.05.002.
This study was funded by National Natural Science Foundation of China (81672452); National Natural Science Foundation of Guangdong Province (2018B0303110015); Science and Technology Program of Guangzhou (201707010354).
Ethics approval and consent to participate
Consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Data distributions with or without normalization of the NPC data sets. Fig. S2 Hub genes of two key modules. Fig. S3 Comparison of expression of the seven validated hub genes. Table S1 Details for the NPC data sets from GEO. Table S2 Number of genes of 22 modules.
About this article
Cite this article
Huang, D., Liu, Q., Zhang, W. et al. Identified IGSF9 association with prognosis and hypoxia in nasopharyngeal carcinoma by bioinformatics analysis. Cancer Cell Int 20, 498 (2020). https://doi.org/10.1186/s12935-020-01587-z