Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

The MYCN-HMGA2-CDKN2A pathway in non-small cell lung carcinoma—differences in histological subtypes

BMC Cancer volume 16, Article number: 71 (2016) Cite this article

Abstract

Background

Extensive research has increased our understanding of the molecular alterations needed for non-small cell lung cancer (NSCLC) development. Deregulation of a pathway including MYCN, HMGA2 and CDKN2A, with the participation of DICER1, is of importance in several solid tumours, and may also be of significance in the pathogenesis of NSCLC.

Methods

Gene expression of MYCN, HMGA2, CDKN2A and DICER1 were investigated with RT-qPCR in surgically resected NSCLC tumour tissue from 175 patients. Expression of the let-7 microRNA family was performed in 78 adenocarcinomas and 16 matching normal lung tissue samples using microarrays. The protein levels of HMGA2 were determined by immunohistochemistry in 156 tumour samples and the protein expression was correlated with gene expression. Associations between clinical data, including time to recurrence, and expression of mRNA, protein and microRNAs were analysed.

Results

Compared to adenocarcinomas, squamous cell carcinomas had a median 5-fold increase in mRNA expression of HMGA2 (p = 0.003). A positive correlation (r = 0.513, p < 0.010) between HMGA2 mRNA expression and HMGA2 protein expression was seen. At the protein level, 90 % of the squamous cell carcinomas expressed high levels of the HMGA2 protein compared to 47 % of the adenocarcinomas (p < 0.0001). MYCN was positively correlated with HMGA2 (p < 0.010) and DICER1 mRNA expression (p < 0.010), and the expression of the let-7 microRNAs seemed to be correlated with the genes studied. MYCN expression was associated with time to recurrence in multivariate survival analyses (p = 0.020).

Conclusions

A significant difference in HMGA2 mRNA expression between the histological subtypes of NSCLC was seen with a higher expression in the squamous cell carcinomas. This was also found at the protein level, and we found a good correlation between the mRNA and the protein expression of HMGA2. Moreover, the expression of MYCN, HMGA2, and DICER1 seems to be correlated to each other and the expression of the let7-genes impacted by their expression. MYCN gene expression seems to be of importance in time to recurrence in this patient cohort with resected NSCLC.

Peer Review reports

Background

The prognosis of lung cancer is dismal. For all stages, 5 year relative survival are 19 % for women and 14 % for men in Norway [1]. Worldwide, lung cancer is the most common cause of cancer related death and it is estimated that 1.6 million people die annually due to lung cancer [2].

Non-small cell lung cancer (NSCLC) constitutes over 80 % of all lung cancers and can be divided into histological subtypes. Investigation to further unveil the NSCLC biology includes molecular characterization. Epidermal growth factor receptor (EGFR) mutational testing has been performed routinely in Norway since 2010. Targeted therapies such as EGFR inhibitors and anaplastic lymphoma kinase (ALK) inhibitors are currently in clinical use. Increased understanding of different levels of tumour development, of genetic, epigenetic, protein alterations and their functional influence are of clinical importance. For a disease where a majority of patients are diagnosed in late stages, there is a need to discover biomarkers as well as new targets for therapeutic interventions.

Overexpression of the HMGA2 gene is linked to the development of cancer [3]. The HMGA2 gene encodes a non-histone chromatin modifying protein that binds to AT-rich regions in the DNA; thus leading to a change in DNA structure and interaction with transcription factors influencing cell growth, proliferation, differentiation and cell death [4]. HMGA proteins are abundantly expressed during fetal development, but scarcely present or even absent in normal adult tissue [5]. In NSCLC, however, re-expression of HMGA2 is proposed as a common event and considered as a molecular marker [6]. Indeed, HMGA2 protein expression is in some studies shown related to lung cancer development and progression, and inversely associated with lung cancer survival [7, 8].

The MYC oncogenes are well characterized participants in cell proliferation, growth, differentiation and apoptosis [9]. One of the MYC family members, MYCN, functions as a positive regulator of LIN28B, a known repressor of the let-7 family of microRNAs [10, 11]. MicroRNAs are small non-coding RNAs that regulate gene expression, making them important players in cancer development and progression [12]. LIN28 protein binds to pre-let-7 s and prevents further processing into mature let-7 microRNAs by DICER1 [13]. HMGA2 is a thoroughly described target of let-7 s and a loss of let-7 microRNAs can lead to HMGA2 overexpression, shown both in cell lines and in solid tumours [14, 15]. Furthermore, through downregulation of the tumour suppressor gene CDKN2A, HMGA2 is shown to be involved in stem cell renewal [16]. CDKN2A level is previously demonstrated as an independent prognostic factor in non-small cell lung carcinoma [17].

The role of a deregulated pathway including MYCN, HMGA2 and CDKN2A with the participation of DICER1 (Fig. 1) has not previously been investigated in NSCLC. In this study we explore the significance of MYCN, HMGA2, CDKN2A and DICER1 in NSCLC tumours by gene expression analysis. HMGA2 protein expression data are evaluated as well as microarray information on let-7 microRNAs, to elucidate possible correlations. Differences in histological or clinical subgroups are assessed. Finally, associations with time to recurrence are examined.

Fig. 1
figure 1

The MYCN, HMGA2, CDKN2A pathway. MYCN functions as a positive regulator of LIN28B, a known repressor of the let-7 family of microRNAs via binding to DICER1. Loss of let-7 microRNAs can lead to HMGA2 overexpression that may result in a downregulation of CDKN2A

Full size image

Methods

Patient cohort

Eligibility criteria for the present study were patients diagnosed with lung cancer stage I-IIIa, treated with curatively intended surgery at Oslo University Hospital – Rikshospitalet, from 2006 to 2010.

Clinical data were collected from the hospital’s medical records. Follow up data procured via questionnaires administered to the patients, through copies of medical records from local hospitals and information from general practitioners. All data were registered in a project database. The patients were followed until death or 5 years past surgery. Tumours were staged according to the Union for International Cancer Control (UICC), TNM 7, and histopathological parameters were retrieved from the pathology reports. Tumours were grouped in three based on histology; adenocarcinomas, squamous cell carcinomas and others (including large cell carcinomas, undifferentiated carcinomas, adeno-squamous carcinomas, small cell carcinomas and carcinoids). Neoadjuvant and adjuvant treatment were given according to guidelines based on TNM stage and the patient’s age.

Tumour samples

The tumour specimens were collected by thoracic surgeons during lung cancer surgery. They were snap-frozen in liquid nitrogen and stored at − 80 °C until further processing.

RNA isolation

Total RNA was extracted from the tumour specimens using standard TRIZOL methods (Invitrogen, Carlsbad, CA). RNA quantity and quality were assessed by NanoDrop ND-1000 spectrometer (NanoDrop technologies) and by RNA integrity numbers (RIN) measured using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) for gene expression

Gene expression of the four selected genes; MYCN, HMGA2, CDKN2A and DICER1 was measured on tumour samples by RT-qPCR using the Applied Biosystems 7900HT Fast Real-Time PCR system (7900HT Fast System, Applied Biosystems, Foster City, CA). Single stranded cDNA was synthesized from totalRNA in samples using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Human ACTB (beta actin) was used as endogenous control and a commercial Ambion breast control (Life Technologies) used as a calibrator in the qPCR. ACTB where chosen as a single endogenous control when the RT-qPCR where performed in 2011–2012 based on earlier satisfactory experience at our institution, and the current perception of endogenous controls in RT-qPCR experiments in general. All reactions were performed as triplicates for each sample. Outliers were omitted from further analysis. Negative and positive controls were included in every run.

The mRNA expression values was quantified as ΔΔCt values, where Ct = threshold cycle, ΔCt = (Ct target mRNA - Ct B-actin) and ΔΔCt = (ΔCt target mRNA - ΔCt calibrator). Relative quantification (fold change) of mRNA expression was calculated by the 2 –ΔΔCt method [18] using the RQ manager software, version 1.2 (Applied Biosystems).

In a total of 175 patients included, gene expression data was generated from 170 tumour samples for MYCN, 132 tumour samples for HMGA2, 171 tumour samples for CDKN2A and 164 tumour samples for DICER1, respectively. A REMARK diagram of the patients is presented in Additional file 1: Figure S1.

Additional RT-qPCR analyses did not pass the quality control or were not possible to perform due to lack of sufficient tumour material.

Immunohistochemistry

HMGA2 protein expression levels in the tumour tissue were determined by immunohistochemistry. We used rabbit anti-HMGA2 antibodies (www.biocheckinc.com) and Dako EnVision Flex + System (K8012) on tissue micro arrays (TMAs) from 156 tumour samples. All samples were represented on the TMAs in duplicates. Protein expression was scored by counting nuclear staining positivity; either as negative (0), less than 10 % tumour cells with strong nuclear staining (1+), 10–50 % tumour cells with strong nuclear staining (2+) or more than 50 % tumour cells with strong nuclear staining (3+). Samples with score 2+ and 3+ were considered overexpressing the HMGA2 protein.

MicroRNA expression

For microRNA profiling, microRNA microarrays from Agilent Technologies (Agilent human microRNA microarray kit release 16.0, 8 x 60 K) was used in 78 tumour tissue samples. For 16 of the cases matching normal lung tissue was also available for analysis. The normal lung tissue was collected from the lung or lobe removed during the operation, at least 10 cm apart from the macroscopic tumour. The microRNA kit encodes 1205 human microRNAs and 144 human viral microRNAs listed in the Sanger miRBase (release 16.0). Arrays were scanned with Agilent Microarray Scanner (Agilent Technologies, Santa Clara, CA) and raw data pre-processed with Agilent Feature Extraction Software (v. 10.7.3.1), with default parameters employed. MicroRNAs detected in less than 10 % of the samples were filtered out, 570 microRNAs, including 10 microRNAs of the let-7family, remained for further analysis. Microarray data were processed by log2 transformation, and normalized between arrays by the 90th percentile method using the Genespring GX analysis Software v.12.1 (Agilent Technology). Some of the results on the microRNA profiling have previously been published by Bjaanæs et al. [19] and characteristics of the patients are presented in Additional file 1: Table S1.

Validation

For validation of the survival analyses, two data sets were used. The first set contained data from mainly adenocarcinomas from our own research group, while the other was composed of data from squamous cell carcinomas obtained via the The Cancer Genome Atlas (TCGA) data portal [20].

We used the SurePrint G3 Human Gene Expression 8x60K Microarray Kit (Agilent Technologies) for mRNA profiling in 187 lung cancer samples, including 184 lung adenocarcinomas. The arrays were scanned with Agilent C scanner (Agilent’s Scan control software, version A.8.4.1) and the dataset extracted with the Agilent Feature Extraction Software. In GeneSpring GX 12 (Agilent Technologies) the data was normalized by using the 75 percentile method. Median follow up for patients still alive in this cohort who had not developed metastasis or a local recurrence was 49 months (range 26–60 months). Forty percent of the patients in this cohort had a recurrence of lung cancer disease during the follow-up time.

A dataset of lung squamous cell carcinoma was obtained through the TCGA data portal where mRNA expression data were available in a total of 508 tumour specimens defined as lung squamous cell carcinoma. Clinical data and details of progression and survival were, however, only available for a total of 280 of the patients on the website, and these patients were thus included in the survival analysis. Median follow up time in the TCGA dataset for patients still alive with no metastasis or local recurrence was 20 months, ranging from 2 to 60 months and 27 % of the patients experienced local recurrence or distant metastasis during follow-up.

Clinical details for both datasets are available in Additional file 1: Table S2 and Additional file 1: Table S3.

Statistical analyses

Data are reported using descriptive statistics with percentages, means, standard deviations, medians and ranges.

Relative quantification of mRNA expression (fold change) of the selected genes; MYCN, HMGA2, CDKN2A and DICER1 in the tumour samples was analysed in different clinical subsets such as histology, stage, sex and smoking history. Differences in median values between groups regarding continuous mRNA expression data were analysed using non-parametric tests, Mann–Whitney U test and/or Kruskal Wallis test where appropriate. Fold change values are not normally distributed and analyses of correlations were done with Spearman’s Rank Order correlation. Expression of HMGA2 protein in different histological entities, were analyzed by Chi-square test. To evaluate differences in mean let-7 microRNA expression in tumour and normal lung tissue, and to see whether mRNA expression were associated with mean microRNA let-7 data, independent t-tests were performed.

For the survival analyses, follow up data existed for all included patients. Time to recurrence of lung cancer was calculated as the interval between the date of operation and date of local relapse, distant metastases, or death due to lung cancer. Patients that died of other causes without local relapse or metastases were censored in the survival analyses (n = 38). One of the participants died 13 days after the operation in cardiac arrest, and was excluded from the survival analyses. Of the remaining 37 patients that died of alternative causes it was verified through medical records that 15 patients died from cardiovascular disease. 15 patients died from other verified diseases such as COPD, pneumonia and septicaemia. For 8 patients no definitive cause of death were uncovered, however neither of these patients had a verified relapse of NSCLC prior to death, nor were there any indication of a lung cancer related death in the medical records procured.

Survival curves and estimation of statistical significance between groups were performed using the Kaplan-Meier method and log rank test, respectively; here mRNA expression (fold change) was categorized into high and low values based on the median gene expression value. Factors associated with lung cancer progression were analyzed using the Cox proportional hazard regression model also using dichotomous fold change data for gene expression. Although not significant in the log rank tests performed, age at surgery, sex, stage and histology, were still included in the multivariate Cox regression analysis, due to high clinical relevance.

A p-value ≤ 0.050 was considered statistically significant. All statistical analyses were performed using SPSS version 21.0 (SPSS Inc., Chicago, MO, USA).

Ethical considerations

This project was approved by the regional ethics committee (Regional comittees for medical and health research ethics - South East) (Approval no:S-06402b) and the institutional review board (Protokollutvalget, Radiumhospitalet). The patients were given oral and written information prior to inclusion. A written consent was obtained from all the participants.

Results

Patient characteristics

One hundred and seventy-five patients were included in the study. The median age at surgery was 66 years (range 34 to 82 years) and 54 % were male. The tumour samples included 57 % adenocarcinomas, 25 % squamous cell carcinomas and 18 % others. According to the current TNM classification, 55 % of the patients were in pTNM stage I, 31 % in stage II and 13 % in stage IIIa. Fifty four of the patients (31 %) received adjuvant chemotherapy, 3 patients received neoadjuvant chemotherapy and 9 patients received radiotherapy (adjuvant or neoadjuvant). During follow up, 60 participants (34 %) presented with relapse of disease; either as a local recurrence or metastatic disease. At the end of follow up 86 patients (49 %) had died. For those whose disease recurred, the median time from initial surgery to diagnosis of relapse was 18 months (range 2 to 60 months). Clinical and pathological features and outcome parameters are summarized in Table 1.

Table 1 Characteristics of the patient cohort
Full size table

Gene expression

We discovered a significant difference in the median values of HMGA2 mRNA expression according to tumour histology (p = 0.003) with a considerably higher median fold change value in the squamous cell carcinoma group (70.6 fold) compared to adenocarcinomas (13.1 fold) and others (5.8 fold) (Table 2). Regarding CDKN2A expression, a significantly higher median expression was found in women compared to men (p = 0.013). Boxplots illustrating the distribution of gene expression is presented in Additional file 1: Figure S2.

Table 2 Relative quantification of mRNA expression (fold change) in tumour sample subsets stratified by histology, stage, sex and smoking history
Full size table

HMGA2 protein expression

In total, 86 (55 %) of the tumour samples were scored with a high expression level of the HMGA2 protein. In histological subsets, 90 % of the squamous cell carcinomas expressed high levels of HMGA2, while 47 % of the adenocarcinomas showed high expression of the protein; this difference was highly significant (p < 0.0001). It is worth noting that none of the 8 carcinoid tumours expressed HMGA2 protein. Protein expression levels according to histological subtypes are presented in Table 3. Microscopy images of HMGA2 protein in NSCLC are shown in Fig. 2.

Table 3 Levels of HMGA2 protein in histological subsets of non-small cell lung cancer
Full size table
Fig. 2
figure 2

Expression of HMGA2 protein in NSCLC. Photomicrographs of NSCLC specimens stained with anti-HMGA2 antibodies. a represents a squamous cell carcinoma sample demonstrating strong immunoreactivity whereas negative staining is shown in b, an adenocarcinoma sample

Full size image

Correlations, mRNA expression and protein

A positive correlation coefficient of 0.513 was identified between HMGA2 mRNA expression and HMGA2 protein expression (p < 0.010). MYCN and HMGA2 mRNA expression showed a positive, albeit weaker correlation (correlation coefficient of 0.232, p < 0.010). DICER1 mRNA expression was positively correlated to MYCN, HMGA2 and CDKN2A mRNA expression with correlation coefficients of 0.463 (p < 0.010), 0.174 (p < 0.050) and 0.167 (p < 0.050) respectively (Additional file 1: Figure S3).

A significant correlation between HMGA2 and CDKN2A mRNA was not seen. The correlation coefficient of CDKN2A mRNA expression and HMGA2 protein was weakly negative, however not significant.

Let-7 microRNAs

All of the measured let-7 s had a mean expression level significantly lower in tumour tissue compared to normal lung tissue. Compared to tumours with low HMGA2 mRNA expression, mean expression of let-7a, let-7c, let-7d and let-7f was significantly lower in tumours with high expression of HMGA2 whereas let-7d was higher. Mean expression of let-7a and let-7d was significantly different depending on mRNA expression in all of the four genes in question, while mean let-7f was significantly different in MYCN, HMGA2 and CKDN2A (Table 4). Details on the independent t-tests performed are included in supplementary tables (Additional file 1: Table S4, Additional file 1: Table S5, Additional file 1: Table S6, Additional file 1: Table S7, Additional file 1: Table S8).

Table 4 Mean let-7 microRNA expression associated with mRNA expression in tumour samples
Full size table

Survival analyses

In survival analyses patients with a MYCN mRNA expression level above median had a significantly worse prognosis compared to patients with MYCN mRNA expression level below median, both in the cohort as a whole (log rank test, p = 0.029) and in patients with squamous cell carcinomas (log rank test, p = 0.044) (Fig. 3). In the multivariate Cox regression analysis adjusting for age, gender, stage and histology, the association between MYCN and time to recurrence was confirmed (p = 0.020) (Table 5). These findings were not, however, validated in a squamous cell carcinoma cohort with clinical follow-up data from the TCGA project or in the cohort with mRNA data from adenocarcinomas.

Fig. 3
figure 3

Associations between gene expression and patient outcome. mRNA fold change values dichotomized to high (green) or low (blue) expression based on the median expression level. High expression of MYCN mRNA had a significantly worse prognosis compared to lower expression levels in all tumour samples (a, p = 0.029) and in a subset of squamous cell carcinomas (b, p = 0.044). Lower expression of both MYCN and HMGA2 (co-expressed) had a significantly better prognosis compared to the other alternatives (c, p = 0.040)

Full size image
Table 5 Predictors for Progression Free Survival (n = 174)
Full size table

In tumours with apparent co-regulation of expression of HMGA2 and MYCN, there was a significantly better outcome when both genes were downregulated compared to the alternatives (log rank test, p = 0.040) (Fig. 3).

Among patients with stage I disease, low HMGA2 protein expression was associated with better prognosis compared to the patients with overexpression of the protein (log rank test, p = 0.034) (Additional file 1: Figure S4). Patients with large cell carcinomas and low expression of the HMGA2 protein also had a better outcome (log rank test, p = 0.014). There was no significant relationship between HMGA2 protein expression and progression free survival among all patients included as a whole. Moreover, the mRNA levels of CDKN2A or DICER did not influence the prognosis in the different histological subgroups.

Discussion

The development of NSCLC includes multiple genetic and epigenetic alterations that may lead to activation of pathways promoting tumour growth as well as inhibition of pathways of tumour suppression. NSCLC has among the greatest number of genetic aberrations of all malignant tumours [21]. This is the first study investigating the pathway involving genes MYCN, HMGA2, CDKN2A, DICER1 and the let-7 family of microRNAs in NSCLC, while several findings have indicated a connected pathway involving these players [14–17].

The HMGA2 protein plays an important role in growth during embryonic development. It is mainly expressed in embryos, and mutant mice with HMGA2 deficiency develop a pygmy phenotype [22]. Close to absent in adult tissue, re-expression in tumours has led to investigations of causality. It is previously known that HMGA2 chromosomal rearrangements are implicated in benign tumours of mesenchymal origin [23] and overexpression of the HMGA2 protein is found in a variety of malignant tumour types such as breast cancer [24], pancreatic cancer [25], oral squamous cell carcinomas [26] as well as NSCLC [6–8]. In our study, analysis of fresh-frozen NSCLC tumour tissue by RT-qPCR revealed a significant difference in HMGA2 mRNA expression in histological subsets; with an elevated median expression value in the squamous cell carcinomas of approximately 5-fold over the median value of adenocarcinomas. A similar pattern was seen for HMGA2 protein; with a significant distinction of nuclear staining in the different histological groups. Close to all squamous cell carcinomas expressed high levels of HMGA2 compared to 46 % of the adenocarcinomas, and this difference was highly significant. Moreover, a strong positive correlation was demonstrated between HMGA2 mRNA and HMGA2 protein expression. Previous studies have indicated a role for thyroid transcription factor-1 (TTF-1) in the regulation of HMGA2 expression, as a loss of TTF-1 triggers overexpression of HMGA2 [27]. Indeed, squamous cell carcinomas are usually negative for TTF-1, while the opposite is shown in adenocarcinomas [28].

An increase in HMGA2 expression (RT-qPCR) in tumour matching non-tumour tissue, with a significant higher increase in squamous cell carcinoma compared to adenocarcinoma, was previously reported by Meyer et al. [6]. These results are consistent with our findings. Many studies highlight the role of HMGA2 protein in cancer progression, HMGA2 protein levels in primary lung tumours has been shown to correlate with increasing tumour grade and has furthermore been regarded as a necessity for a transformed phenotype in metastatic NSCLC cell lines [7]. Accordingly, HMGA2 protein expression was entirely devoid in the slow growing carcinoids in our study and HMGA2 protein levels in stage I patients were associated with progression free survival supporting the role of HMGA2 as a marker of aggressiveness.

We identified a positive correlation between MYCN and HMGA2 mRNA. In addition, DICER1 expression was positively correlated to MYCN, HMGA2 and CDKN2A, supporting the notion of the involvement of this pathway. When focusing on the tumours with seemingly co-regulated HMGA2 and MYCN expression, we found a significantly better outcome when both genes were downregulated compared to the remainders (Fig. 3). This also supports the involvement of this pathway in lung carcinomas, and has not been demonstrated previously. However, it is known in general that the genes investigated are also regulated by several other factors and that the contribution of HMGA2, in particular, is complex, involving several other genes and proteins [4].

Let-7 microRNA expression is scarce in embryonic stages, but increases in tissue following mature differentiation and is suggested as putative tumor suppressors [29]. As confirmed by our study, normal lung tissue showed a higher level of let-7 expression compared to adenocarcinoma tumor samples. Moreover, HMGA2 is a well characterized target of the let-7 microRNA family and known to be inversely correlated with HMGA2 expression in NSCLC cells [30, 31]. We found several let-7 microRNAs to be differentially expressed comparing high and low HMGA2 gene expression; let-7a, let-7c and let-7f was seen inversely correlated, consistent with the biological presumption that loss of let-7 inhibition leads to HMGA2 overexpression in cancer. It has also been shown that aberrant expression of let-7 is more common in squamous cell carcinoma compared to adenocarcinomas [32]. Our cohort included only a few squamous cell carcinoma samples, but we still received a similar end result.

Although we identified a significant impact of MYCN-expression levels on progression free survival in both the cohort investigated as a whole, and in the squamous cell carcinomas alone, this was not validated in the TCGA data. There may be several explanations for this lack of validation. The TCGA project was initiated for molecular analyses, and the clinical data was missing in many of the available cases. Moreover, a fairly large proportion of patients with existing clinical data had only a short term follow-up. Our cohort in the present study consisted of several histological subgroups, and each subgroup had limited number of cases. We propose that similar survival analyses need to be further investigated in larger cohorts of homogenous subtypes in future studies.

The mRNA level of HMGA2, CDKN2A and DICER1 did not influence survival in our study significantly. Previous studies have not investigated this fully in specific histological lung cancer subgroups, but have shown that the HMGA2 protein is involved in the transformation of lung cancer cells. Its role might therefore very well be in the initiation and establishment of the cancer. In neuroblastoma, MYCN, LIN28B and let-7 seem to be involved and MYCN amplification is correlated with adverse clinical outcome [33, 34]. This is a more benign disease although the location of tumours may unfavourably impact the prognosis.

One of the limitations to this study is the low number of samples in our subgroup analyses. In addition, cohorts with solid clinical data and follow up would ensure further investigations of the prognostic impact activation of this pathway in lung cancers, and in the different histological subtypes.

Conclusion

In this study we have identified a significant difference between the histological subtypes of NSCLC and HMGA2 expression at both the mRNA and protein level. The more benign histology of carcinoids lacked HMGA2 expression, while among squamous cell carcinomas, most samples showed high expression. We have also demonstrated a correlation between mRNA expression of HMGA2 and HMGA2 protein expression. The mRNA expression levels of MYCN and HMGA2 and DICER1 were significantly correlated, and co-regulation of expression of HMGA2 and MYCN, had a significant impact on survival. The impact on survival is probably complex, and needs to be investigated in larger cohorts of specific histological non-small cell lung cancer subgroups.

Abbreviations

ACTB:

Actin-beta

ALK:

Anaplastic lymphoma kinase

CDKN2A:

Cyclin-dependent kinase inhibitor 2A

EGFR:

Epidermal growth factor receptor

HMGA2:

High mobility group AT-hook 2

MYCN:

V-myc avian myelocytomatosis viral oncogene neuroblastoma derived homolog

NSCLC:

Non-small cell lung cancer

RT-qPCR:

Reverse transcription quantitative polymerase chain reaction

TCGA:

The Cancer Genome Atlas

TMA:

Tissue micro array

TNM:

TNM classification of malignant tumours

TTF-1:

Thyroid transcription factor-1

UICC:

Union for International Cancer Control

References

  1. Cancer in Norway 2013 [http://www.kreftregisteret.no/no/Generelt/Publikasjoner/Cancer-in-Norway/]

  2. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: Sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2014;136:E359–86.

    Article  PubMed  Google Scholar 

  3. Fedele M, Fusco A. HMGA and cancer. Biochim Biophys Acta. 2010;1799:48–54.

    Article  CAS  PubMed  Google Scholar 

  4. Cleynen I, Van den Ven W. The HMGA proteins: a myriad of functions (Review). Int J Oncol. 2008;32:289–305.

    CAS  PubMed  Google Scholar 

  5. Rogalla P, Drechsler K, Frey G, Hennig Y, Helmke B, Bonk U, et al. HMGI-C expression patterns in human tissues. Implications for the genesis of frequent mesenchymal tumors. Am J Pathol. 1996;149:775–9.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. Meyer B, Loeschke S, Schultze A, Weigel T, Sandkamp M. HMGA2 Overexpression in non-small cell lung cancer. Mol Carcinog. 2007;46:503–11.

    Article  CAS  PubMed  Google Scholar 

  7. Di Cello F, Hillion J, Hristov A, Wood LJ, Schuldenfrei A, Kowalski J, et al. HMGA2 Participates in transformation in human lung cancer. Mol Cancer Res. 2008;6:743–50.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Sarhadi VK, Wikman H, Salmenkivi K, Kuosma E, Sioris T, Salo J, et al. Increased expression of high mobility group a proteins in lung cancer. J Pathol. 2006;209:206–12.

    Article  CAS  PubMed  Google Scholar 

  9. Meyer N, Penn LZ. Reflecting on 25 years with MYC. Nat Rev Cancer. 2008;8:976–90.

    Article  CAS  PubMed  Google Scholar 

  10. Cotterman R, Knoepfler PS. N-Myc regulates expression of pluripotency genes in neuroblastoma including lif, klf2, klf4, and lin28b. PLoS One. 2009;4(6):e5799.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Lockhart VL, Shah SP, Tanwar PS, Mermel CH, Beroukhim R, Azam M, et al. Lin28 promotes transformation and is associated with advanced human malignancies. Nat Genet. 2009;41:843–8.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Bartels CL, Tsongalis GJ. Mini-reviews micrornas:novel biomarkers for human cancer. Clin Chem. 2009;55:623–31.

    Article  CAS  PubMed  Google Scholar 

  13. Heo I, Joo C, Cho J, Ha M, Han J, Kim VN. Lin28 mediates the Terminal Uridylation of let-7 Precursor MicroRNA. Mol Cell. 2008;32:276–84.

    Article  CAS  PubMed  Google Scholar 

  14. Mayr C, Hemann MT, Bartel DP. Disrupting the pairing between let-7 and Hmga2 enhances oncogenic transformation. Science. 2007;315:1576–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Johnsen H, Holm R, Kristensen GB, Birrer MJ, Pearson RB, Børresen-Dale A-L, et al. Deregulation of MYCN, LIN28B and LET7 in a molecular subtype of aggressive high-grade serous ovarian cancers. PLoS One. 2011;6:e18064.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Nishino J, Kim I, Chada K, Morrison S. Hmga2 promotes neural stem cell self-renewal in young but not old mice by reducing p16Ink4a and p19Arf expression. Cell. 2008;135:227–39.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Tong J, Sun X, Cheng H, Zhao D, Ma J, Zhen Q, et al. Expression of p16 in non-small cell lung cancer and its prognostic significance: a meta-analysis of published literatures. Lung Cancer. 2011;74:155–63.

    Article  PubMed  Google Scholar 

  18. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

  19. Bjaanaes MM, Halvorsen AR, Solberg S, Jørgensen L, Dragani TA, Galvan A, et al. Unique microRNA-profiles in EGFR-mutated lung adenocarcinomas. Int J Cancer. 2014;135:1812–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. The Cancer Genome Atlas [http://cancergenome.nih.gov]

  21. Welch JS, Walter MJ, Wendl MC, Ley TJ, Wilson RK, Raphael BJ, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502:333–9.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Zhou X, Benson KF, Ashar HR, Chada K. Mutation responsible for the mouse pygmy phenotype in the developmentally regulated factor HMGI-C. Nature. 1995;376(6543):771–4.

    Article  CAS  PubMed  Google Scholar 

  23. Schoenmakers EF, Wanschura S, Mols R, Bullerdiek J, Van den Berghe H, Van de Ven WJ. Recurrent rearrangements in the high mobility group protein gene, HMGI-C, in benign mesenchymal tumours. Nat Genet. 1995;10:436–44.

    Article  CAS  PubMed  Google Scholar 

  24. Rogalla P, Drechsler K, Kazmierczak B, Rippe V, Bonk U, Bullerdiek J. Expression of HMGI-C, a member of the high mobility group protein family, in a subset of breast cancers: relationship to histologic grade. Mol Carcinog. 1997;19:153–6.

    Article  CAS  PubMed  Google Scholar 

  25. Atomi Y, Fusco A, Chiappetta G, Sugiyama M, Abe N, Watanabe T, et al. An increased high-mobility group A2 expression level is associated with malignant phenotype in pancreatic exocrine tissue. Br J Cancer. 2003;89:2104–9.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Miyazawa J, Mitoro A, Kawashiri S, Chada KK, Imai K. Expression of Mesenchyme-Specific gene HMGA2 in squamous cell carcinomas of the oral cavity. Cancer Res. 2004;64:2024–9.

    Article  CAS  PubMed  Google Scholar 

  27. Rice SJ, Lai SC, Wood LW, Helsley KR, Runkle EA, Winslow MM, et al. MicroRNA-33a mediates the regulation of high mobility group AT-hook 2 gene (HMGA2) by thyroid transcription factor 1 (TTF-1/NKX2-1). J Biol Chem. 2013;288:16348–60.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Bresaola E, Valduga F, Lupo C, Viale G, Terzi A, Iannucci A, et al. Immunoreactivity for thyroid transcription factor-1 in stage I non-small cell carcinomas of the lung. Am J Surg Pathol. 2001;25:363–72.

    Article  PubMed  Google Scholar 

  29. Jacks T, Sharp PA, Ebert MS, Chen CY, Pester RE, Kumar MS, et al. Suppression of non-small cell lung tumor development by the let-7 microRNA family. Proc Natl Acad Sci U S A. 2008;105:3903–8.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Lee YS, Dutta A. The tumor suppressor microRNA let-7 represses the HMGA2 oncogene. Genes Dev. 2007;21:1025–30.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Park SM, Shell S, Radjabi AR, Schickel R, Feig C, Boyerinas B, et al. Let-7 prevents early cancer progression by suppressing expression of the embryonic gene HMGA2. Cell Cycle. 2007;6:2585–90.

    Article  CAS  PubMed  Google Scholar 

  32. Capodanno A, Boldrini L, Proietti A, Alì G, Pelliccioni S, Niccoli C, et al. Let-7 g and miR-21 expression in non-small cell lung cancer: Correlation with clinicopathological and molecular features. Int J Oncol. 2013;43:765–74.

    CAS  PubMed  Google Scholar 

  33. Molenaar JJ, Domingo-Fernández R, Ebus ME, Lindner S, Koster J, Drabek K, et al. LIN28B induces neuroblastoma and enhances MYCN levels via let-7 suppression. Nat Genet. 2012;44:1199–206.

    Article  CAS  PubMed  Google Scholar 

  34. Schwab M, Westermann F, Hero B, Berthold F. Neuroblastoma: biology and molecular and chromosomal pathology. Lancet Oncol. 2003;4:472–80.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Ingjerd Solvoll, study nurse, for help in collecting tumour samples and obtaining the clinical information for the patients included, Ellen Hellesylt for doing the immunohistochemical staining on tumour samples and Marius Lund-Iversen for providing the photomicrographs of NSCLC specimens stained with anti-HMGA2 antibodies. Sources of support: This work was supported by the Norwegian Cancer Society (grant # 51120 to ÅH).

Author information

Authors and Affiliations

  1. Department of Cancer Genetics, Institute for Cancer Research, Oslo University Hospital-The Norwegian Radium Hospital, Oslo, Norway

    Hanne A. Eide, Ann Rita Halvorsen, Maria Moksnes Bjaanæs, Odd Terje Brustugun & Åslaug Helland

  2. Department of Oncology, Oslo University Hospital-The Norwegian Radium Hospital, Oslo, Norway

    Hanne A. Eide, Maria Moksnes Bjaanæs, Odd Terje Brustugun, Cecilie Essholt Kiserud & Åslaug Helland

  3. Cellular and Molecular Research Center, Qazvin University of Medical Sciences, Qazvin, Iran

    Hossein Piri

  4. Department of Pathology, Oslo University Hospital-The Norwegian Radium Hospital, Oslo, Norway

    Ruth Holm

  5. Department of Cardiothoracic Surgery, Oslo University Hospital-Rikshospitalet, Oslo, Norway

    Steinar Solberg & Lars Jørgensen

  6. Department of Oncology, National Advisory Unit on Late Effects After Cancer Treatment, Oslo University Hospital-The Norwegian Radium Hospital, Oslo, Norway

    Cecilie Essholt Kiserud

Authors
  1. Hanne A. Eide
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ann Rita Halvorsen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Maria Moksnes Bjaanæs
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Hossein Piri
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ruth Holm
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Steinar Solberg
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lars Jørgensen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Odd Terje Brustugun
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Cecilie Essholt Kiserud
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Åslaug Helland
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Åslaug Helland.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HAE performed the data analysis and wrote the manuscript. ARH and HP performed the RT-qPCR. ARH and MMB performed the microRNA microarrays and provided the adenocarcinoma validation set. RH evaluated the immunostained sections. SS, LJ and OTB provided patient material and patient data. CEK participated in the data analysis and in writing the manuscript. ÅH conceived the study, provided patient material and data, and participated in writing the manuscript. All authors read and approved the final manuscript.

Additional file

Additional file 1: Table S1.

Characteristics of the microRNA sample set. Table S2: Characteristics of the adenocarcinoma validation sample set. Table S3: Characteristics of the squamous cell carcinoma (TCGA) validation sample set. Table S4: Independent samples t-test. Mean let-7 microRNA expression in tumour and normal tissue samples. Table S5: Independent samples t-test. Mean let-7 microRNA expression and association with MYCN mRNA expression in tumour samples. Table S6: Independent samples t-test. Mean let-7 microRNA expression and association with HMGA2 mRNA expression in tumour samples. Table S7: Independent samples t-test. Mean let-7 microRNA expression and association with CDKN2A mRNA expression in tumour samples. Table S8: Independent samples t-test. Mean let-7 microRNA expression and association with DICER1 mRNA expression in tumour samples and additional figures (Figure S1: REMARK diagram detailing sample availability and use of different analytical techniques in the present study. Figure S2: Boxplots illustrating the distribution of gene expression (fold change) from the RT-qPCR . Figure S3: Scatterplots illustrating significant correlations. Figure S4: Associations between HMGA2 protein expression and patient outcome). Figure S1. REMARK diagram detailing sample availability and use of different analytical techniques in the present study. Figure S2. Boxplots illustrating the distribution of gene expression from the RT-qPCR. Figure S3. Scatterplots illustrating significant correlations. Figure S4. Associations between NSCLC stage, HMGA2 protein expression and patient outcome. HMGA2 protein expression values dichotomized to low expression (blue) and overexpression (green) based on immunohistochemistry. Low expression of HMGA2 protein had a significantly better prognosis compared to overexpression in stage I non-small cell lung cancer tumour samples (A, p = 0.034). No significant association in stage II (B, p = 0.492) or stage III (C, p = 0.862) patients was seen. (DOCX 401 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Eide, H.A., Halvorsen, A.R., Bjaanæs, M.M. et al. The MYCN-HMGA2-CDKN2A pathway in non-small cell lung carcinoma—differences in histological subtypes. BMC Cancer 16, 71 (2016). https://doi.org/10.1186/s12885-016-2104-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-016-2104-9

Keywords

BMC Cancer

ISSN: 1471-2407

Contact us