|Ahead of print publication
Is hypermethylation of SOX1 gene an independent prognostic marker in surgically resected non-small cell lung cancer?
Milica Kontic1, Dragana Jovanovic1, Izidor Kern2, Heather H Nelson3, Svetlana Bojic4, Miodrag Ognjanovic5, SImona Ognjanovic6
1 Clinic for Pulmonology, Clinical Center of Serbia; Department of Internal Medicine – Pulmonology, School of Medicine, University of Belgrade, Belgrade, Serbia
2 Department for Pathology, The University Clinic of Pulmonary and Allergic Diseases Golnik, Golnik, Slovenia
3 Division of Epidemiology and Community Health, Masonic Cancer Center, University of Minnesota, Minneapolis, Minnesota
4 Institute for Medical Statistics and Informatics, School of Medicine, University of Belgrade, Belgrade, Serbia
5 International Organization for Cancer Prevention and Research, Belgrade, Serbia
6 Mayo Clinic, Mayo Graduate School, Rochester, New York, USA
|Date of Submission||30-Jan-2020|
|Date of Decision||28-May-2020|
|Date of Acceptance||10-Aug-2020|
|Date of Web Publication||20-Aug-2021|
Clinic for Pulmonology, Clinical Center of Serbia, Belgrade, Serbia
Source of Support: None, Conflict of Interest: None
Background: Promoter hypermethylation of tumor suppressor genes presents promising markers for lung cancer diagnosis and prognosis. The purpose of this study was to determine the association between the promoter hypermethylation of multiple genes and 5-year survival rate in patients with Non-small cell lung cancer (NSCLC).
Materials and Methods: Primary tumor samples (n = 65), corresponding nonmalignant lung tissues (n = 65), and circulating blood were obtained from NSCLC patients who underwent curative surgery. Promoter methylation status in seven genes (RASSF1A, CDH13, MGMT, ESR1, DAPK, SOX1, and HOXA9) was quantified by using bisulfite pyrosequencing. Five-year survival data were obtained by a clinician. Cox proportional hazards models were used to analyze the associations between gene methylation status and overall patient survival.
Results: The 5-year survival of the patients with SOX1 aberrant tumor methylation was found to be statistically significantly shorter than for those patients without aberrant tumor methylation (P = 0.01). This effect was independent of TNM stage. No significant survival differences were associated with aberrant methylation in other genes tested in either of the two tissue types.
Conclusion: Our study shows that SOX1 promoter hypermethylation in NSCLC tumors is significantly associated with inferior survival, showing promise as a useful prognostic biomarker in patients with NSCLC.
Keywords: 5-year survival rate, methylation, non-small cell lung cancer, SOX1
|How to cite this URL:|
Kontic M, Jovanovic D, Kern I, Nelson HH, Bojic S, Ognjanovic M, Ognjanovic S. Is hypermethylation of SOX1 gene an independent prognostic marker in surgically resected non-small cell lung cancer?. J Can Res Ther [Epub ahead of print] [cited 2022 Jul 3]. Available from: https://www.cancerjournal.net/preprintarticle.asp?id=324182
| > Introduction|| |
For decades, lung cancer remains one of the most common causes of cancer-related death worldwide, associated with over 1.3 million deaths per year. While the treatment options for non-small cell lung cancer (NSCLC) continues to evolve, the prognosis is still very poor.
The poor prognosis is associated with late disease diagnosis and the small number of effective drugs. Good prognostic biomarkers for lung cancer patients are still not available-even patients with Stage I NSCLC who undergo surgical resection have high risk of dying from recurrent disease, with a 5-year relapse rate of 35%–50%.
It is well established that epigenetic control of gene expression plays an important role in carcinogenesis. Aberrant methylation of CpG dinucleotides is a commonly observed epigenetic modification in human cancer; Hypermethylation in the promoters of tumor suppressor genes is the mechanism of transcriptional silencing equivalent to mutations. Epigenetic changes appear to be even more frequently found in human cancer than genetic mutations. These methylation changes have been extensively studied in lung cancer, with some of them have been pointed as potential biomarkers with clinical relevance for diagnosis, prognosis, and response to therapy.
The establishment of new prognostic markers holds promise for making more informed choices for cancer treatment and more precise prediction of therapeutic response,, even though our study design described below did not support evaluation of prediction of therapeutic response. To achieve these aims, it will be necessary to demonstrate the specificity of epigenetic changes to lung cancerous cells, and to evaluate the sensitivity of detecting aberrations in easily accessible tissue. This approach requires investigation of specific gene methylation patterns that might be helpful for identification of NSCLC prognostic markers. While a number of genes have been reported as aberrantly methylated in lung tumors,, we have focused on the seven genes frequently hypermethylated in lung cancer (RASSF1A, CDH13, MGMT, ESR1, HOXA9, SOX1, and DAPK), for which we have previously shown that their methylation is associated with clinicopathological features of NSCLC tumors. Once the 5-year survival data collection has been completed, it opened the opportunity for us to investigate whether there is a correlation between promoter hypermethylation of these genes and 5-year survival rate in NSCLC patients after a surgical resection.
| > Materials And Methods|| |
In 2009, patients with resectable non-small cell lung cancer (NSCLC) were approached for participation in a study at the Institute for lung diseases, Clinical Center of Serbia, University of Belgrade. The study received the approval of the center's ethics committee, and all the participants provided informed consent. Sixty-five patients undergoing surgical resection for non-small cell lung cancer consented to tissue collection and a small amount of tumor specimen, normal lung tissue distant to the tumor, and blood samples were collected during surgery and immediately snap frozen for research purposes. None of the patients have received neoadjuvant chemotherapy or radiotherapy.
Demographic information was abstracted from medical records and additional information, was obtained through patient interviews. Patient follow-up was conducted by a research nurse or dedicated physician in order to obtain the data on 5-year survival. Biospecimens collection included primary tumor samples (n = 65), corresponding nonmalignant lung tissues (n = 65), and matching blood samples (n = 51). Five milliliters of blood was collected using EDTA vacutainers and stored at − 20°C until processing. DNA was extracted from whole-blood samples by using a well-established method. Similarly, DNA extraction from fresh frozen tumors using standard methods, which we previously summarized. The quality of DNA was analyzed by electrophoresis on agarose gel, while Nanodrop spectrophotometer (ThermoScientific Inc., Wilmington, DE, USA) was used for DNA quantification. The methylation analyses were carried out at the University of Minnesota Masonic Cancer Center using a well-established method of bisulfide sequencing. Briefly, genomic DNA is initially treated with sodium bisulfite, which converts unmethylated cytosines to uracil, while methylated cytosines stay unaltered. This bisulfite modification step was performed according to the manufacturer's protocol (Zymo Research, Irvine, CA, USA). Upon completion of bisulfite modification of DNA, a strand-specific polymerase chain reaction (PCR) product is generated to provide a suitable DNA template for pyrosequencing. Primer sequences and PCR conditions have been described elsewhere,, and show that multiple products were sequence per each gene analyzed. Amplicons were resolved by agarose electrophoresis to confirm proper amplification and quality of product. Percent methylation for each CpG as well as average methylation across CpG's was calculated for each promoter using PyroMark software (Qiagen, Germantown, MD 20874, Maryland, USA).
Statistical analyses included calculations for each CpG site, as well as average over multiple CpG sites for each gene. Methylation values were not normally distributed, thus Wilcoxon signed-rank test on average methylation over all CpG sites within a gene were used to assess differences in methylation between tumor and normal tissue. When the analysis was repeated using individual CpG sites instead of average for each gene, the obtained results were similar. Compared to the normal tissue, tumor was considered hypermethylated when the average methylation across different CpG sites for a specific gene was more than 3 standard deviations higher than the corresponding average methylation observed in the normal tissue, as previously described. Methylation in blood samples was compared between patients with hypermethylated tumors and those without by using Wilcoxon rank-sum test.
Cox proportional hazard regression was used to evaluate clinicopathologic factors (age, gender, smoking, histology type, histology grade, and TNM stage) as predictors of 5-year survival. Multivariate Cox proportional hazard regression was used to evaluate methylation in specific genes as predictors of 5-year survival. TNM was the covariate used in these models. Kaplan–Meier method (log rank test) was used to evaluate differences in overall survival for the groups, defined as having values above and below the cutoff, for methylation levels of each gene. The optimal cutoff values were previously estimated from the data, as described by Dardis. R package version 0.5.5, which included Miscellaneous Functions for Survival Data package (survMisc) was used for these analyses.
| > Results|| |
Patient characteristics are shown in [Table 1]. The majority were male (76.1%), with long-term smoking history (>40 pack years). At the time of the surgery, 81.3% of the patients were current smokers, while the rest were former smokers. All patients were older than 55 years, with almost half being over 65 years old at the time of the surgery. The predominant histology type was squamous cell carcinoma (56%), followed by adenocarcinoma. The remaining histologies included giant-cell carcinoma (n = 3), atypical carcinoids (n = 1), mucinous carcinoma (n = 1), and carcinoma sarcomatoid (n = 1), while large cell carcinoma was not included in this study.
Average methylation across all CpG sites within a gene was significantly higher in tumor tissue compared to normal adjacent tissue for all genes evaluated in this study [Table 2]. The only exception was ERS1, where the use of the second primer did not reach statistical significance (P = 0.07). When blood samples were compared between the patients with tumor hypermethylation for a specific gene and those without hypermethylation, no statistically significant difference was observed for any of the genes examined. Therefore, these results show that elevated methylation levels observed in genes SOX1, RASSF1A, HOXA9, CDH13, and DAPK in NSCLC were cancer specific, while the lack of reflection of these methylation changes in patients blood indicate their poorly suitability for a screening test.
|Table 2: Mean CpG methylation and standard deviation at all analyzed genes|
Click here to view
Five-year survival was not significantly associated with age, gender, smoking history, tumor histology or grade. The only statistically significant predictor of survival was TNM stage (P = 0.011), as shown in [Table 3]. Thus, TNM stage was the only covariate used to adjust Cox-proportional hazard models which evaluated the effect of tumor methylation level for specific genes on 5-year survival. For SOX1 gene, the model identified a cutoff value for tumor methylation at 29.11, which was higher that the average tumor methylation value of 26.44 [shown in [Table 2]. This cutoff value that clearly distinguished two groups of patients in terms of 5-year survival: those with tumor methylation above the cutoff value [marked as Poz on [Figure 1]a showing survival curves] had poorer survival compared to those marked Neg, which had tumor methylation values below the cutoff. Conversely, when two groups were compared in Cox proportional hazard models adjusted for TNM stage, hazard ratio = 2.371 was observed with 95% confidence interval of 1.001–5.647, which was statistically significant (P = 0.05).
|Table 3: Association between five-year survival and clinical-pathological patients chatacteristics|
Click here to view
|Figure 1: Kaplan–Meier survival curve, according to tumor methylation status for SOX1 (a) and RASSF1A (b) gene. Positive (Poz) and negative (Neg) survival curves depict subjects with tumors with methylation above (Poz) and those below (Neg) respective cutoff values for SOX1 and RASSF1A|
Click here to view
By comparison, the cutoff value of 9.01 identified for RASSF1A gene (second primer), was lower than the average tumor methylation observed for this gene [10.44, shown in [Table 2]. The Kaplan–Meier survival curves [Figure 1]b, did not show any favorable survival of the group with less methylated tumors (those below the 9.01 cutoff).
| > Discussion|| |
We used bisufite pyrosequencing method including 7-gene panel to compare methylation in patient-matched tumor and nontumor tissue. We observed hypermethylation in tumor tissue for all the genes examined, while the methylation in normal tissue (adjacent to tumor) was similar to the one observed in matched patient's blood, rending blood not a suitable indicator of methylation in tumor tissue. We also evaluated whether the observed tumor hypermethylation was associated with 5-year survival of NSCLC in this patient population. Patients with hypermethylation in SOX1 gene in tumor tissue, which was above the average value of methylation across all tumors, had a significantly poorer survival.
SOX1 hypermethylation in NSCLC has been previously reported. These authors have used methylation array method and identified SOX1 as the top five most hypermethylated genes, with over 5-fold larger methylation in tumor compared to normal tissue and confirmed these findings in two independent cohorts. Similarly, Zhao et al. reported SOX1, as one of the several hypermethylated genes in stage I NSCLC. These represent rare reports of SOX1 methylation in NSCLC and none of these examined the relationship between SOX1 hypermethlation and survival. SOX1 antibodies have been reported as a serological marker of a specific small cell carcinoma type, but no relationship between this marker and survival has been observed. In addition to lung cancer, SOX1 hypermethylation has been reported in several different malignancies, including ovarian, esophageal, prostate cancer, and bladder cancer.
SOX1 (sex-determining region Y [SRY] related high-mobility group bo × 1) encodes a transcription factor belongs to a SRY box gene family of genes, which plays an important role in embryonal development. SOX1 is evolutionary conserved displaying similar protein structure in many species and has been implicated as one of the key regulators of neurodevelopment. Its role in cancer has only recently begun to emerge through the reports of genetic silencing in a number of tumor tissues, as discussed above. The plausible mechanism underlying the role of SOX1 in carcinogenesis was proposed by Li and Li, by suggesting that epigenetic silencing of SOX1 in NSCLC was associated with cell migration. Namely, the authors showed that restoration of SOX1 inhibited cell migration by regulating actin cytoskeletal remodeling in NSCLC.
Several limitations of our study need to be noted. All participants in this study were either former or current smokers, thus tissue sample adjacent to tumor, which was considered normal tissue sample, may have been affected by long-term exposure to tobacco carcinogens. To limit the effect of cancer field, the normal tissue sample was attempted to be obtained distant from the resected tumor mass. However, in smokers, all of lung tissue would be exposed to tobacco carcinogens, thus we acknowledge that this approach may not have been impactful. The methylation observed in normal lung tissue was very similar to that in the matched blood sample, while being significantly lower compared to that observed in the tumor tissue, thus it is unlikely that this limitation impacted our results. On the other hand, the approach of comparing methylation in tumor to normal lung tissue has the potential to reflect the changes that occur along the specter of neoplastic transformation and thus increase the utility of such biomarkers for screening. Small sample size and the use of samples from a single center warrant that our results need to be confirmed in our, preferably larger studies.
To our knowledge, this is the first report of association between SOX1 methylation and NSCLC survival, while there are only a handful of publications associating SOX1 tumor methylation and poorer survival in other malignancies, as noted above. In conclusion, hypermethylation of SOX1 gene in NSCLC correlates with poor prognosis of patients. It may serve as novel epigenetic-based diagnostic biomarker with further clinical impact for risk stratification of NSCLC patients. Caution is warranted due to small sample size and further studies are needed to confirm the specificity of this marker to NSCLC and therefore whether this could be utilized as a prognostic biomarker of NSCLC survival.
All authors contributed to this manuscript.
We thank Ms. Jacquelyn Kuriger for performing bisulfide sequencing.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| > References|| |
Jemal A, Murray T, Samuels A, Ghafoor A, Ward E, Thun MJ. Cancer statistics, 2003. CA Cancer J Clin 2003;53:5-26.
Siegel R, Naishadham D, Jemal A. Cancer statistics, 2012. CA Cancer J Clin 2012;62:10-29.
Kim JS, Han J, Shim YM, Park J, Kim DH. Aberrant methylation of H-cadherin (CDH13) promoter is associated with tumor progression in primary nonsmall cell lung carcinoma. Cancer 2005;104:1825-33.
Sandoval J, Mendez-Gonzalez J, Nadal E, Chen G, Carmona FJ, Sayols S, Moran S, et al
. A prognostic DNA methylation signature for stage i non-small-cell lung cancer. J Clin Oncol 2013;31:4140-47.
Bowman RV, Wright CM, Davidson MR, Francis SM, Yang IA, Fong KM. Epigenomic targets for the treatment of respiratory disease. Expert Opin Ther Targets 2009;13:625-40.
Kontic M, Stojsic J, Jovanovic D, Bunjevacki V, Ognjanovic S, Kuriger J, et al
. Aberrant promoter methylation of CDH13 and MGMT genes is associated with clinicopathologic characteristics of primary non-small-cell lung carcinoma. Clin Lung Cancer 2012;13:297-303.
Higuchi R, Bowman B, Freiberger M, Ryder OA, Wilson AC. DNA sequences from the Quagga, an extinct member of the horse family. Nature 1984;312:282-4.
Saito K, Kawakami K, Matsumoto I, Oda M, Watanabe G, Minamoto T. Long interspersed nuclear element 1 hypomethylation is a marker of poor prognosis in stage IA non-small cell lung cancer. Clin Cancer Res 2010;16:2418-26.
Toyooka KO, Toyooka S, Virmani AK, Sathyanarayana UG, Euhus DM, Gilcrease M, et al
. Loss of expression and aberrant methylation of the CDH13 (H-cadherin) gene in breast and lung carcinomas. Cancer Res 2001;61:4556-60.
van Engeland M, Roemen GM, Brink M, Pachen MM, Weijenberg MP, de Bruïne AP, et al
. K-ras mutations and RASSF1A promoter methylation in colorectal cancer. Oncogene 2002;21:3792-5.
Buckingham L, Penfield Faber L, Kim A, Liptay M, Barger C, Basu S, et al
. PTEN, RASSF1 and DAPK site-specific hypermethylation and outcome in surgically treated stage I and II nonsmall cell lung cancer patients. Int J Cancer 2010;126:1630-9.
Nelson HH, Marsit CJ, Christensen BC, Houseman EA, Kontic M, Wiemels JL, et al
. Key epigenetic changes associated with lung cancer development: Results from dense methylation array profiling. Epigenetics 2012;7:559-66.
Zhao Y, Zhou H, Ma K, Sun J, Feng X, Geng J, et al
. Abnormal methylation of seven genes and their associations with clinical characteristics in early stage non-small cell lung cancer. Oncol Lett 2013;5:1211-8.
Titulaer MJ, Klooster R, Potman M, Sabater L, Graus F, Hegeman IM, et al
. SOX antibodies in small-cell lung cancer and Lambert-Eaton myasthenic syndrome: Frequency and relation with survival. J Clin Oncol 2009;27:4260-7.
Su HY, Lai HC, Lin YW, Chou YC, Liu CY, Yu MH. An epigenetic marker panel for screening and prognostic prediction of ovarian cancer. Int J Cancer 2009;124:387-93.
Kuo IY, Chang JM, Jiang SS, Chen CH, Chang IS, Sheu BS, et al
. Prognostic CpG methylation biomarkers identified by methylation array in esophageal squamous cell carcinoma patients. Int J Med Sci 2014;11:779-87.
Mathews LA, Hurt EM, Zhang X, Farrar WL. Epigenetic regulation of CpG promoter methylation in invasive prostate cancer cells. Mol Cancer 2010;9:267-84.
Lopez JI, Angulo JC, Martin A, Sanchez-Chapado M, Gonzalez-Corpas A, et al
. A DNA hypermethylation profile reveals new potential biomarkers for the evaluation of prognosis in urothelial bladder cancer. APMIS 2017;125:787-96.
Bowles J, Schepers G, Koopman P. Phylogeny of the SOX family of developmental transcription factors based on sequence and structural indicators. Dev Biol 2000;15:239-55.
Sarkar A, Hochedlinger K. The sox family of transcription factors: Versatile regulators of stem and progenitor cell fate. Cell Stem Cell 2013;12:15-30.
Li N, Li S. Epigenetic inactivation of SOX1 promotes cell migration in lung cancer. Tumour Biol 2015;36:4603-10.
[Table 1], [Table 2], [Table 3]