Hypoxia drives shared and distinct transcriptomic changes in two invasive glioma stem cell lines.

Glioblastoma (GBM) is the most common primary malignant cancer of the central nervous system. Insufficient oxygenation (hypoxia) has been linked to GBM invasion and aggression, leading to poor patient outcomes. Hypoxia induces gene expression for cellular adaptations. However, GBM is characterized by high intertumoral (molecular subtypes) and intratumoral heterogeneity (cell states), and it is not well understood to what extent hypoxia triggers patient-specific gene responses and cellular diversity in GBM. Here, we surveyed eight patient-derived GBM stem cell lines for invasion phenotypes in 3D culture, which identified two GBM lines showing increased invasiveness in response to hypoxia. RNA-seq analysis of the two patient GBM lines revealed a set of shared hypoxia response genes concerning glucose metabolism, angiogenesis, and autophagy, but also a large set of patient-specific hypoxia-induced genes featuring cell migration and anti-inflammation, highlighting intertumoral diversity of hypoxia responses in GBM. We further applied the Shared GBM Hypoxia gene signature to single cell RNA-seq datasets of glioma patients, which showed that hypoxic cells displayed a shift towards mesenchymal-like (MES) and astrocyte-like (AC) states. Interestingly, in response to hypoxia, tumor cells in IDH-mutant gliomas displayed a strong shift to the AC state, whereas tumor cells in IDH-wildtype gliomas mainly shifted to the MES state. This distinct hypoxia response of IDH-mutant gliomas may contribute to its more favorable prognosis. Our transcriptomic studies provide a basis for future approaches to better understand the diversity of hypoxic niches in gliomas.

Plexin-B2 facilitates glioblastoma infiltration by modulating cell biomechanics.

Infiltrative growth is a major cause of high lethality of malignant brain tumors such as glioblastoma (GBM). We show here that GBM cells upregulate guidance receptor Plexin-B2 to gain invasiveness. Deletion of Plexin-B2 in GBM stem cells limited tumor spread and shifted invasion paths from axon fiber tracts to perivascular routes. On a cellular level, Plexin-B2 adjusts cell adhesiveness, migratory responses to different matrix stiffness, and actomyosin dynamics, thus empowering GBM cells to leave stiff tumor bulk and infiltrate softer brain parenchyma. Correspondingly, gene signatures affected by Plexin-B2 were associated with locomotor regulation, matrix interactions, and cellular biomechanics. On a molecular level, the intracellular Ras-GAP domain contributed to Plexin-B2 function, while the signaling relationship with downstream effectors Rap1/2 appeared variable between GBM stem cell lines, reflecting intertumoral heterogeneity. Our studies establish Plexin-B2 as a modulator of cell biomechanics that is usurped by GBM cells to gain invasiveness.

Akaluc bioluminescence offers superior sensitivity to track in vivo glioma expansion.

BACKGROUND: Longitudinal tracking of tumor growth using noninvasive bioluminescence imaging (BLI) is a key approach for studies of in vivo cancer models, with particular relevance for investigations of malignant gliomas in rodent intracranial transplant paradigms. Akaluciferase (Akaluc) is a new BLI system with higher signal strength than standard firefly luciferase (Fluc). Here, we establish Akaluc BLI as a sensitive method for in vivo tracking of glioma expansion. METHODS: We engineered a lentiviral vector for expression of Akaluc in high-grade glioma cell lines, including patient-derived glioma stem cell (GSC) lines. Akaluc-expressing glioma cells were compared to matching cells expressing Fluc in both in vitro and in vivo BLI assays. We also conducted proof-of-principle BLI studies with intracranial transplant cohorts receiving chemoradiation therapy. RESULTS: Akaluc-expressing glioma cells produced more than 10 times higher BLI signals than Fluc-expressing counterparts when examined in vitro, and more than 100-fold higher signals when compared to Fluc-expressing counterparts in intracranial transplant models in vivo. The high sensitivity of Akaluc permitted detection of intracranial glioma transplants starting as early as 4 h after implantation and with as little as 5000 transplanted cells. The sensitivity of the system allowed us to follow engraftment and expansion of intracranial transplants of GSC lines. Akaluc was also robust for sensitive detection of in vivo tumor regression after therapy and subsequent relapse. CONCLUSION: Akaluc BLI offers superior sensitivity for in vivo tracking of glioma in the intracranial transplant paradigm, facilitating sensitive approaches for the study of glioma growth and response to therapy.

Teaching an Old Drug New Tricks: Dexamethasone as an In Vivo Inhibitor of Glioblastoma Dispersal.

Identifying drugs that can mitigate dispersal of glioblastoma cells, particularly after patients undergo radiotherapy and concomitant chemotherapy, may increase the length of time to recurrence and improve overall survival. Previous studies have shown that dexamethasone (Dex), a drug currently used to treat brain tumor-related edema, which is tapered immediately after the edema has resolved, induces fibronectin matrix assembly (FNMA) and reduces dispersal of primary human glioblastoma multiforme (GBM) cells in vitro and ex vivo. Here, we utilized an in vivo mouse retina dispersal assay to demonstrate that Dex also inhibits dispersal in vivo. We show that 1) Dex significantly reduces z-axis penetration of glioblastoma cells into mouse retina; 2) treatment alters the morphology of dispersal; 3) without Dex, the presence of fibronectin increases dispersal; 4) treatment activates in vivo FNMA by glioblastoma cells, leading to the containment of the tumor mass; and 5) Dex-mediated activation of FNMA is fibronectin dose-dependent. Dispersal inhibition could be achieved at human equivalent doses as low as 1 mg/day, a dose significantly lower than currently used to reduce edema. This is the first step towards future studies in which patients can be potentially maintained on low-dose dexamethasone therapy with the aim of increasing the time between initial resection and recurrence.

Gene signatures of quiescent glioblastoma cells reveal mesenchymal shift and interactions with niche microenvironment.

BACKGROUND: Glioblastoma (GBM), a highly malignant brain tumor, invariably recurs after therapy. Quiescent GBM cells represent a potential source of tumor recurrence, but little is known about their molecular underpinnings. METHODS: Patient-derived GBM cells were engineered by CRISPR/Cas9-assisted knock-in of an inducible histone2B-GFP (iH2B-GFP) reporter to track cell division history. We utilized an in vitro 3D GBM organoid approach to isolate live quiescent GBM (qGBM) cells and their proliferative counterparts (pGBM) to compare stem cell properties and therapy resistance. Gene expression programs of qGBM and pGBM cells were analyzed by RNA-Seq and NanoString platforms. FINDINGS: H2B-GFP-retaining qGBM cells exhibited comparable self-renewal capacity but higher therapy resistance relative to pGBM. Quiescent GBM cells expressed distinct gene programs that affect cell cycle control, metabolic adaptation, and extracellular matrix (ECM) interactions. Transcriptome analysis also revealed a mesenchymal shift in qGBM cells of both proneural and mesenchymal GBM subtypes. Bioinformatic analyses and functional assays in GBM organoids established hypoxia and TGFß signaling as potential niche factors that promote quiescence in GBM. Finally, network co-expression analysis of TCGA glioma patient data identified gene modules that are enriched for qGBM signatures and also associated with survival rate. INTERPRETATION: Our in vitro study in 3D GBM organoids supports the presence of a quiescent cell population that displays self-renewal capacity, high therapy resistance, and mesenchymal gene signatures. It also sheds light on how GBM cells may acquire and maintain quiescence through ECM organization and interaction with niche factors such as TGFß and hypoxia. Our findings provide a starting point for developing strategies to tackle the quiescent population of GBM. FUND: National Institutes of Health (NIH) and Deutsche Forschungsgemeinschaft (DFG).

Analysis of chromatin accessibility uncovers TEAD1 as a regulator of migration in human glioblastoma.

The intrinsic drivers of migration in glioblastoma (GBM) are poorly understood. To better capture the native molecular imprint of GBM and its developmental context, here we isolate human stem cell populations from GBM (GSC) and germinal matrix tissues and map their chromatin accessibility via ATAC-seq. We uncover two distinct regulatory GSC signatures, a developmentally shared/proliferative and a tumor-specific/migratory one in which TEAD1/4 motifs are uniquely overrepresented. Using ChIP-PCR, we validate TEAD1 trans occupancy at accessibility sites within AQP4, EGFR, and CDH4. To further characterize TEAD's functional role in GBM, we knockout TEAD1 or TEAD4 in patient-derived GBM lines using CRISPR-Cas9. TEAD1 ablation robustly diminishes migration, both in vitro and in vivo, and alters migratory and EMT transcriptome signatures with consistent downregulation of its target AQP4. TEAD1 overexpression restores AQP4 expression, and both TEAD1 and AQP4 overexpression rescue migratory deficits in TEAD1-knockout cells, implicating a direct regulatory role for TEAD1-AQP4 in GBM migration.

Prospective Isolation and Comparison of Human Germinal Matrix and Glioblastoma EGFR+ Populations with Stem Cell Properties.

Characterization of non-neoplastic and malignant human stem cell populations in their native state can provide new insights into gliomagenesis. Here we developed a purification strategy to directly isolate EGFR+/- populations from human germinal matrix (GM) and adult subventricular zone autopsy tissues, and from de novo glioblastoma (GBM) resections, enriching for cells capable of binding EGF ligand (LBEGFR+), and uniquely compared their functional and molecular properties. LBEGFR+ populations in both GM and GBM encompassed all sphere-forming cells and displayed proliferative stem cell properties in vitro. In xenografts, LBEGFR+ GBM cells showed robust tumor initiation and progression to high-grade, infiltrative gliomas. Whole-transcriptome sequencing analysis confirmed enrichment of proliferative pathways in both developing and neoplastic freshly isolated EGFR+ populations, and identified both unique and shared sets of genes. The ability to prospectively isolate stem cell populations using native ligand-binding capacity opens new doors onto understanding both normal human development and tumor cell biology.

The significance of PIWI family expression in human lung embryogenesis and non-small cell lung cancer.

The expression of Piwi-interacting RNAs, small RNAs that bind to PIWI proteins, was until recently believed to be limited to germinal stem cells. We have studied the expression of PIWI genes during human lung embryogenesis and in paired tumor and normal tissue prospectively collected from 71 resected non-small-cell lung cancer patients. The mRNA expression analysis showed that PIWIL1 was highly expressed in 7-week embryos and downregulated during the subsequent weeks of development. PIWIL1 was expressed in 11 of the tumor samples but in none of the normal tissue samples. These results were validated by immunohistochemistry, showing faint cytoplasmic reactivity in the PIWIL1-positive samples. Interestingly, the patients expressing PIWIL1 had a shorter time to relapse (TTR) (p = 0.006) and overall survival (OS) (p = 0.0076) than those without PIWIL1 expression. PIWIL2 and 4 were downregulated in tumor tissue in comparison to the normal tissue (p < 0.001) and the patients with lower levels of PIWIL4 had shorter TTR (p = 0.048) and OS (p = 0.033). In the multivariate analysis, PIWIL1 expression emerged as an independent prognostic marker. Using 5-Aza-dC treatment and bisulfite sequencing, we observed that PIWIL1 expression could be regulated in part by methylation. Finally, an in silico study identified a stem-cell expression signature associated with PIWIL1 expression.

A serum microRNA signature associated with complete remission and progression after autologous stem-cell transplantation in patients with multiple myeloma.

We have examined serum microRNA expression in multiple myeloma (MM) patients at diagnosis and at complete response (CR) after autologous stem-cell transplantation (ASCT), in patients with stable monoclonal gammopathy of undetermined significance, and in healthy controls. MicroRNAs were first profiled using TaqMan Human MicroRNA Arrays. Differentially expressed microRNAs were then validated by individual TaqMan MicroRNA assays and correlated with CR and progression-free survival (PFS) after ASCT. Supervised analysis identified a differentially expressed 14-microRNA signature. The differential expression of miR-16 (P = 0.028), miR-17 (P = 0.016), miR-19b (P = 0.009), miR-20a (P = 0.017) and miR-660 (P = 0.048) at diagnosis and CR was then confirmed by individual assays. In addition, high levels of miR-25 were related to the presence of oligoclonal bands (P = 0.002). Longer PFS after ASCT was observed in patients with high levels of miR-19b (6 vs. 1.8 years; P < 0.001) or miR-331 (8.6 vs. 2.9 years; P = 0.001). Low expression of both miR-19b and miR-331 in combination was a marker of shorter PFS (HR 5.3; P = 0.033). We have identified a serum microRNA signature with potential as a diagnostic and prognostic tool in MM.

Differential MIR-21 expression in plasma from mesenteric versus peripheral veins: an observational study of disease-free survival in surgically resected colon cancer patients.

Findings on the role of plasma miR-21 expression in colorectal cancer are contradictory. Before reaching a peripheral vein (PV), microRNAs released by the tumor are dispersed throughout the body. We hypothesized that blood drawn from the mesenteric vein (MV) near the site of the primary tumor could provide more homogeneous information than blood drawn from the PV.We have analyzed miR-21 expression in matched samples of tumor tissue, normal tissue, MV plasma, and PV plasma in 57 surgically resected patients with colon cancer and correlated our findings with clinical characteristics and disease-free survival (DFS).miR-21 expression was higher in MV than PV plasma (P = 0.014) and in tumor than in normal tissue (P < 0.001). Patients with high levels of miR-21 in MV plasma had shorter DFS (P = 0.05) than those with low levels, and those with high levels in both MV and PV plasma had shorter DFS than all other patients (P = 0.01).Our findings suggest that the primary tumor in colon cancer releases high concentrations of miR-21 in the MV but that these concentrations are later diluted in the circulatory system. MV expression of miR-21 may be a stronger prognostic marker than PV expression.

miR-141 and miR-200c as markers of overall survival in early stage non-small cell lung cancer adenocarcinoma.

BACKGROUND: Several treatments in non-small cell lung cancer (NSCLC) are histology-dependent, and the need for histology-related markers is increasing. MicroRNAs (miRNAs) are promising molecular markers in multiple cancers and show differences in expression depending on histological subtype. The miRNA family miR-200 has been associated with the regulation of epithelial-mesenchymal (EMT)/mesenchymal-epithelial transition (MET). EMT involves profound phenotypic changes that include the loss of cell-cell adhesion, the loss of cell polarity, and the acquisition of migratory and invasive properties that facilitates metastasis. A dual role for the miR-200 family in the prognosis of several tumors has been related to tumor cell origin. However, the prognostic role and function of miR-200 family in early-stage NSCLC adenocarcinoma and squamous cell carcinoma (SCC) have not been well established. METHODS: miRNA expression was determined using TaqMan assays in 155 tumors from resected NSCLC patients. Functional studies were conducted in three NSCLC cell lines: H23, A-549 and HCC-44. RESULTS: High miR-200c expression was associated with shorter overall survival (OS) in the entire cohort (p = 0.024). High miR-200c (p = 0.0004) and miR-141 (p = 0.009) expression correlated with shorter OS in adenocarcinoma - but not in SCC. In the multivariate analysis, a risk score based on miR-141 and miR-200c expression emerged as an independent prognostic factor for OS in the entire cohort (OR, 2.787; p = 0.033) and in adenocarcinoma patients (OR, 10.649; p = 0.002). Functional analyses showed that miR-200c, was related to mesenchymal-epithelial transition (MET) and affected cell migration and E-cadherin levels, while overexpression of miR-141 reduced KLF6 protein levels and produced an increase of secretion of VEGFA in vitro (H23, p = 0.04; A-549, p = 0.03; HCC-44, p = 0.02) and was associated with higher blood microvessel density in patient tumor samples (p<0.001). CONCLUSION: High miR-141 and miR-200c expression are associated with shorter OS in NSCLC patients with adenocarcinoma through MET and angiogenesis.

High levels of global DNA methylation are an independent adverse prognostic factor in a series of 90 patients with de novo myelodysplastic syndrome.

The prognostic impact of global DNA methylation and hydroxymethylation was assessed in 90 patients with de novo myelodysplastic syndrome (MDS). DNA was isolated from bone marrow samples obtained at diagnosis and global methylation and hydroxymethylation were determined by ELISA. Patients with a percentage of methylated DNA above 2.73% had a shorter overall survival than those with lower levels (P=0.018) and presented a negative trend in terms of leukemia-free survival (P=0.084), that was statistically significant after censoring 9 patients that received disease-modifying treatments both in univariate and multivariate analyses. Similarly, the low-risk MDS patients defined by the IPSS, WPSS and IPSS-R with 5-mC percentage in total DNA above 2.73% had a shorter overall survival (P=0.032; P=0.023; P=0.031). No cut-off value for the 5-hydroxymethylcytosine percentage with statistical significance for overall or leukemia-free survival was obtained. This study suggests that global DNA methylation predicts overall survival in myelodysplastic syndromes.

Role of miR-200 family members in survival of colorectal cancer patients treated with fluoropyrimidines.

BACKGROUND AND OBJECTIVES: Surgery is the standard treatment for colorectal cancer (CRC), and adjuvant chemotherapy has been shown to be effective in stage III but less so in stage II. We have analyzed the expression of the miR-200 family in tissue samples from resected CRC patients and correlated our findings with survival to adjuvant treatment with fluoropyrimidines. METHODS: Tumor tissue samples were obtained from 127 surgically resected patients with stage I-III CRC. miRNA detection was performed using TaqMan MicroRNA assays. RESULTS: High levels of miR-200a and miR-200c were associated with longer overall survival, while high levels of miR-429 correlated with longer overall and disease-free survival (DFS). In the subgroup of 56 patients treated with fluoropyrimidines and in the smaller subgroup of 32 stage II patients treated with fluoropyrimidines, those with high levels of miR-200a, miR-200c, miR-141, or miR-429 had significantly longer overall and DFS. Low miR-429 levels were identified as an independent prognostic marker. High levels of miR-429 combined with 5-fluorouracil inhibited cell invasion in LOVO cells. CONCLUSIONS: miR-200a, miR-200c, miR-141, and miR-429 expression levels may identify CRC patients, including those with stage II disease, who are most likely to benefit from adjuvant chemotherapy.

B cell activation through CD40 and IL4R ligation modulates the response of chronic lymphocytic leukaemia cells to BAFF and APRIL.

The two tumour necrosis factor family proteins BAFF (TNFSF13B) and APRIL (TNFSF13) and their receptors [BAFF-R (TNFRSF13C), TACI (TNFRSF13B), BCMA (TNFRSF17)] play a critical role in the survival of normal B cells. The sensitivity of normal B cells to BAFF and APRIL can be modulated by signals regulated by their receptors. This modulation, however, has not been extensively investigated in chronic lymphocytic leukaemia (CLL) cells. We evaluated the expression, regulation and signalling of BAFF and APRIL receptors in normal and in CLL cells upon stimulation through CD40+IL4R and BCR. We further analysed the prognostic value of BAFF and APRIL receptors expression in patients with CLL. BCMA expression was significantly higher on CLL cells than on normal B cells. BCR and CD40+IL4R stimulation promoted an increase in TACI and BCMA expression, cell viability and activation in normal B cells. A similar effect was observed in CLL cells after CD40+IL4R but not BCR stimulation. BCMA expression correlated with unmutated IGHV genes, poor-risk cytogenetics, and short progression-free survival. These findings further characterize the link between CD40+IL4R regulatory signals, BAFF, APRIL and their receptors and the survival of leukaemic cells and clinical features of CLL.

MiR-SNPs as markers of toxicity and clinical outcome in Hodgkin lymphoma patients.

BACKGROUND: In recent years, microRNA (miRNA) pathways have emerged as a crucial system for the regulation of tumorogenesis. miR-SNPs are a novel class of single nucleotide polymorphisms that can affect miRNA pathways. DESIGN AND METHODS: We analyzed eight miR-SNPs by allelic discrimination in 141 patients with Hodgkin lymphoma and correlated the results with treatment-related toxicity, response, disease-free survival (DFS) and overall survival (OS). RESULTS: The KRT81 (rs3660) GG genotype was associated with an increased risk of neurological toxicity (P = 0.016), while patients with XPO5 (rs11077) AA or CC genotypes had a higher rate of bleomycin-associated pulmonary toxicity (P = 0.048). Both miR-SNPs emerged as independent factors in the multivariate analysis. The XPO5 AA and CC genotypes were also associated with a lower response rate (P = 0.036). XPO5 (P = 0.039) and TRBP (rs784567) (P = 0.022) genotypes emerged as prognostic markers for DFS, and XPO5 was also associated with OS (P = 0.033). In the multivariate analysis, only XPO5 emerged as an independent prognostic factor for DFS (HR: 2.622; 95%CI 1.039-6.620; P = 0.041). Given the influence of XPO5 and TRBP as individual markers, we then investigated the combined effect of these miR-SNPs. Patients with both the XPO5 AA/CC and TRBP TT/TC genotypes had the shortest DFS (P = 0.008) and OS (P = 0.008). CONCLUSION: miR-SNPs can add useful prognostic information on treatment-related toxicity and clinical outcome in Hodgkin lymphoma and can be used to identify patients likely to be chemoresistant or to relapse.

Low miR-145 and high miR-367 are associated with unfavourable prognosis in resected nonsmall cell lung cancer.

The transcription factors SRY-related HMG box (SOX)2 and octamer-binding transcription factor (OCT)4 regulate the expression of the miR-302-367 cluster. miR-145 regulates SOX2 and OCT4 translation and p53 regulates miR-145 expression. We analysed the expression of the miR-302-367 cluster and miR-145 and the mutational status of p53 in resected nonsmall cell lung cancer (NSCLC) patients and correlated results with time to relapse (TTR). Tumour and paired normal tissue samples were obtained from 70 NSCLC patients. MicroRNA expression was assessed with TaqMan MicroRNA Assays. p53 exons 5 to 8 were sequenced. miR-145 was downregulated (p<0.0001) and miR-367 was upregulated (p<0.0001) in tumour compared with normal tissue. Mean TTR was 18.4 months for patients with low miR-145 levels and 28.2 months for those with high levels (p=0.015). Mean TTR was 29.1 months for patients with low miR-367 levels and 23.4 months for those with high levels (p=0.048). TTR was shorter for patients with both unfavourable variables (p=0.009). Low miR-145 expression (p=0.049), the combination of unfavourable microRNA levels (p=0.02) and the combination of low miR-145 with p53 mutations (p=0.011) were independent markers of shorter TTR. In conclusion, miR-145 and miR-367 expression could be novel markers for relapse in surgically treated NSCLC. p53 may play a role in modulating miR-145 expression in NSCLC.

Impact of MiRSNPs on survival and progression in patients with multiple myeloma undergoing autologous stem cell transplantation.

PURPOSE: A distinctive new group of polymorphisms is constituted by single-nucleotide polymorphism (SNP) in miRNA processing machinery in miRNA precursor molecules and in miRNA-binding sites, known as miRSNPs. The aim of this study was to ascertain the prognostic impact of six miRSNPs in patients with multiple myeloma and analyze the functional consequences. EXPERIMENTAL DESIGN: One hundred and thirty-seven patients with chemosensitive multiple myeloma (73M/64F) intensified with autologous stem cell transplantation (ASCT) were studied. The median follow-up was 4 years. The genes and SNPs evaluated in genomic DNA by allelic discrimination were KRT81 (rs3660), AFF1 (rs17703261), FAM179b (rs1053667), and MIR196A2 (rs11614913) for miRNA target genes and TRBP (rs784567) and XPO5 (rs11077) for miRNA biogenesis pathway. RESULTS: Overall survival (OS) was significantly longer in patients with KRT81 rs3660 C/C variant (P = 0.037). Functional analysis showed that the presence of C variant in KRT81 3' untranslated region (UTR) is related with a reduction of the protein levels. Moreover, the reduction of KRT81 protein levels by siRNA in multiple myeloma cell lines is related to a decreased proliferation. On the other hand, OS was significantly longer in patients with C/C or A/C variant in XPO5 rs11077 (P = 0.012). There was also a significantly longer progression-free survival (PFS) for this SNP (P = 0.013). This SNP retained its prognostic impact on PFS and OS in a multivariate regression analysis (P = 0.028 and P = 0.014, respectively). CONCLUSION: This is the first report that relates miRSNPs with prognosis in multiple myeloma either in a keratin gene (KRT81), target of diverse miRNA multiple myeloma clusters, or in the miRNA biogenesis pathway-related protein exportin-5.

A dual role for KRT81: a miR-SNP associated with recurrence in non-small-cell lung cancer and a novel marker of squamous cell lung carcinoma.

MicroRNAs (miRNAs) play an important role in carcinogenesis through the regulation of their target genes. miRNA-related single nucleotide polymorphisms (miR-SNPs) can affect miRNA biogenesis and target sites and can alter microRNA expression and functions. We examined 11 miR-SNPs, including 5 in microRNA genes, 3 in microRNA binding sites and 3 in microRNA-processing machinery components, and evaluated time to recurrence (TTR) according to miR-SNP genotypes in 175 surgically resected non-small-cell lung cancer (NSCLC) patients. Significant differences in TTR were found according to KRT81 rs3660 (median TTR: 20.3 months for the CC genotype versus 86.8 months for the CG or GG genotype; P = 0.003) and XPO5 rs11077 (median TTR: 24.7 months for the AA genotype versus 73.1 months for the AC or CC genotypes; P = 0.029). Moreover, when patients were divided according to stage, these differences were maintained for stage I patients (P = 0.002 for KRT81 rs3660; P<0.001 for XPO5 rs11077). When patients were divided into sub-groups according to histology, the effect of the KRT81 rs3660 genotype on TTR was significant in patients with squamous cell carcinoma (P = 0.004) but not in those with adenocarcinoma. In the multivariate analyses, the KRT81 rs3660 CC genotype (OR = 1.8; P = 0.023) and the XPO5 rs11077 AA genotype (OR = 1.77; P = 0.026) emerged as independent variables influencing TTR. Immunohistochemical analyses in 80 lung specimens showed that 95% of squamous cell carcinomas were positive for KRT81, compared to only 19% of adenocarcinomas (P<0.0001). In conclusion, miR-SNPs are a novel class of SNPs that can add useful prognostic information on the clinical outcome of resected NSCLC patients and may be a potential key tool for selecting high-risk stage I patients. Moreover, KRT81 has emerged as a promising immunohistochemical marker for the identification of squamous cell lung carcinoma.