Abstract
MicroRNAs, negative post-transcriptional regulators of gene expression, are involved in cancer.
Their role in early bronchial carcinogenesis was analysed in 60 biopsies obtained by fluorescence bronchoscopy (six per stage: normal tissue of nonsmokers, normal normofluorescent and hypofluorescent bronchial tissue of smokers, hyperplasia, metaplasia, mild, moderate and severe dysplasia, in situ carcinoma and invasive squamous cell carcinoma (SQCC)).
In total, 69 microRNAs were found to be differentially expressed in the course of bronchial carcinogenesis. Among them, some microRNAs showed a linear evolution of their expression level, such as miR-32 and miR-34c, whose expression progressively decreased from normal bronchial tissues of nonsmokers to SQCC. Others behaved differently at successive stages, such as miR-142-3p or miR-9, or are only altered from a specific stage, such as miR-199a or miR-139. MicroRNAs globally followed a two-step evolution, first decreasing (a reverse of their increase during embryogenesis) during the earliest morphological modifications of bronchial epithelium, and thereafter increasing at later stages of lung carcinogenesis. Moreover, microRNA expression was very efficient for the prediction of the histological classification between low- and high-grade lesions and between in situ and invasive carcinoma.
The present data show, for the first time, that microRNAs are involved in bronchial carcinogenesis from the very early steps of this process and, thus, could provide tools for early detection of lung cancer.
Each year, lung cancer kills 1,300,000 persons worldwide and its incidence is steadily increasing. Improving its cure rate is thus a major public health objective. The bad prognosis of this cancer is mainly explained by the fact that the diagnosis is generally made only at advanced stages. Therefore, identification of early biomarkers, especially if they are detectable by noninvasive methods (serum, sputum, etc.) in high-risk patients (mainly smokers), is mandatory. Lung squamous cell carcinoma (SQCC) is the end-point of a whole range of morphological abnormalities 1 that are diffusely displayed in the bronchial airways of smokers (concept of field cancerisation) and that could be used to identify markers of the ongoing cancer process.
Discovery of microRNAs (miRNAs), negative post-transcriptional regulators of gene expression, was a landmark in molecular biology. Recent studies demonstrated that mutations in miRNAs or their aberrant expression are associated with diverse human diseases, including cancers. Indeed, miRNAs might act as oncogenes or tumour suppressor genes. Many miRNAs are located on fragile sites and genomic regions involved in cancer 2. They play important roles in cell differentiation, cell growth and cell death 3. They have been found to be involved in known oncogenic pathways like those of p53 4, Bcl2 5 or K-Ras 6. Finally, miRNAs seem to be very significant prognostic factors in patients with different tumours 7–10 and could be useful for treatment 11.
Interestingly, miRNAs maintain, even in cancer, the specificity of their developmental lineage 7. Expression of miRNAs is tissue and cell-type regulated 12. Moreover, it has also been shown that, even if different tumour types display highly different patterns of miRNA expression, tumours with a common embryonic precursor still express common miRNA patterns 7, 13. Therefore, their pattern of expression very efficiently classifies different tumour types in a more robust way than gene expression microarrays 7.
The current authors studied the evolution of miRNA levels of expression in biopsies at the successive stages of human bronchial squamous carcinogenesis in order to determine their potential role in this pathological process and to identify potential tools for early detection of lung cancer.
MATERIAL AND METHODS
Sample collection and RNA isolation
All biopsies were obtained between 2003 and 2007 after informed consent from the patients and approval by the local ethical committee (Institut Jules Bordet, Brussels, Belgium). Selected patients had a minimum smoking exposure of 30 pack-yrs and/or a history of lung or head and neck cancer. Fluorescence bronchoscopy was performed under local anaesthesia and all hypofluorescent areas were biopsied for histopathological diagnosis. Biopsies were classified, using the 1999 histological World Health Organization/International Association for the Study of Lung Cancer classification of pre-invasive and invasive squamous lesions of the bronchus 14, into microscopically normal bronchial epithelium, hyperplasia, metaplasia, mild, moderate and severe dysplasia, carcinoma in situ (CIS) and invasive SQCC. Normal bronchial epithelium was biopsied both in normo- and hypofluorescent areas in smokers and in six nonsmokers with healthy lungs. A total of 60 biopsies were included, with six samples of each category: normal bronchial epithelium of nonsmokers, normal normofluorescent bronchial epithelial tissue of smokers, histologically normal but hypofluorescent bronchial epithelium of smokers, hyperplasia, metaplasia, mild, moderate and severe dysplasia, CIS and SQCC. In each case, two biopsies were taken with a clean forceps in the same area. One biopsy was used for histopathology and the second biopsy was immediately dropped into Tripure Isolation Reagent (Roche Diagnostics, Indianapolis, IN, USA), lysed by ultraturaxing and frozen in liquid nitrogen. All samples were kept at −80°C. RNA extraction was performed according to a classical phenol/chloroform protocol as described by the Tripure provider. As carrier, 20 µg of glycogen (Roche Diagnostics) was added, and the separation between the organic and the aqueous phase was performed on Phase Lock Gel (Eppendorf, Hamburg, Germany) to optimise the recovery of nucleic acids. Isolated RNAs were assessed for quantity and purity on the NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Rockland, DE, USA) and for quality on the Agilent 2100 bioanalyser with RNA 6000 NanoAssay (Agilent Technologies, Palo Alto, CA, USA). RNAs were stored at −80°C.
MiRNA expression analyses
The miRNAs were retro-transcribed and amplified (by PCR) using the multiplex RT TaqMan MicroRNA Low Density Array (LDA) Assays (Applied Biosystems, Foster City, CA, USA). For each sample, 160 ng of starting total RNA (20 ng for each of the eight RT-PCRs) was used. The global miRNAs profiling, for 365 human miRNAs, was then performed using the TaqMan LDA Human microRNA Panel v1.0 (Applied Biosystems). All the quality control tests were validated: blanks and reproducibility (standard deviation of cycle threshold (CT) <1) of the two small nucleolar house-keeping RNAs, RNU48 (SNORD48) and RNU44 (SNORD44). The amount of RNA from each sample was calibrated to the more stable (between the different arrays) small nucleolar house-keeping RNA, RNU48. This value gave a delta CT (ΔCT) value for each miRNA (miRNA CT value minus RNU48 CT value). The average ΔCTs were calculated for each group of six samples and the delta ΔCTs (ΔΔCTs), corresponding to the differences between two groups, were obtained by subtracting the average ΔCT of the second group from those of the first one. Fold differences were calculated as 2(−ΔΔCT), as a decrease in 1 CT value was equivalent to a two-fold increase in the starting amount of cDNA. For the down-regulated miRNAs, it was calculated by 1/2(−ΔΔCT).
Statistical assessment
First, unsupervised average linkage hierarchical clustering including all samples was performed to exclude any experimental bias. Secondly, unpaired t-tests were used to assess the statistical significance of the difference in miRNA expression between two groups. In addition, a correction for multiple testing (CMT) was performed according to Bonferroni's family-wise error rate. For both tests, the threshold for statistical significance was p<0.05. Finally, in order to assess the discriminating potential of miRNA expression profiles, clustering using centred correlation and average linkage was performed.
RESULTS
Differentially expressed miRNAs during squamous cell carcinogenesis
Comparisons of miRNAs between successive stages did not show any statistically significant differences, probably due to the small number of patients included at each stage. Therefore, the nine histological stages were combined in larger groups according to a molecular classification previously performed on the same biopsies. These groups are described in table 1⇓.
Composition of groups
First, miRNA expression profiles were compared between normal bronchial epithelium of nonsmokers (group 0) and other groups of biopsies at different stages of squamous carcinogenesis in smokers (groups 1–3). MiRNA expression profiles between successive groups and subgroups in smokers (groups 1, 2, 3, 3A and 3B; table 1⇑) were then successively compared, and comparisons were also made between CIS and SQCC. Eight comparisons were thus successively performed between groups 0 versus 1, 2 or 3; 1 versus 2; 2 versus 3 or 3A; 3A versus 3B and CIS versus 3B. In total, 69 miRNAs were found to be significantly differentially expressed in at least one of the eight comparisons. The names of these 69 miRNAs and the stage(s) at which their expression is significantly modified are described in table 2⇓. In the eight successive comparisons of miRNA expression profiles, six, 19, 19, 16, 30, 15, 20 and 10 miRNAs were found to be significantly differentially expressed. In particular, expression of miR-34c, miR-15a and miR-32 was progressively reduced from group 0 to groups 1, 2 and 3. Additionally, 16 of the 69 miRNAs remained significantly differentially expressed in one or more comparison(s) after CMT: zero, three, six, four, five, four, one and two miRNAs at each of the eight successive comparisons.
Statistically significantly evolving microRNAs during lung carcinogenesis
Two-step evolution of miRNA expression
As illustrated in table 3⇓, two successive waves were globally observed in miRNA expression profiles during lung squamous carcinogenesis: first, a downregulation and, secondly, an upregulation of the majority of miRNAs. In the first four comparisons of miRNA profiles, between groups 0 and 1, 0 and 2, 0 and 3, and 1 and 2, the majority of miRNAs were downregulated: 83, 84, 74 and 87.5%, respectively. In the later four comparisons, the proportion of downregulated miRNAs decreased and that of upregulated miRNAs increased. This increase in the proportion of upregulated miRNAs began in the comparison between groups 2 and 3, and became ≥50% in the comparisons between groups 2 and 3A, 3A and 3B, and CIS and 3B.
Two-step evolution of microRNAs in lung squamous carcinogenesis
While some miRNAs were modified at one specific step of lung squamous carcinogenesis only (e.g. miR-199a or miR139), others evolved across successive steps. Among them, some miRNAs showed a two-wave evolution, all being firstly downregulated and then upregulated (miR-142-3p, miR-142-5p, miR-9, miR-10b and miR-214).
Sample group classification by miRNA profiles
Using the miRNAs that were differentially expressed between groups 0, 1 and 2, it was not possible to segregate samples a posteriori according their respective groups. However, the miRNAs differentially expressed between group 2 and 3 (fig. 1⇓), as well as between groups 2 and 3A (fig. 2⇓), and CIS and invasive carcinomas (group 3B; fig. 3⇓), were very efficient in enabling the separation and classification of individual samples of these groups. The first and second clustering experiments, based on the differentially expressed miRNAs between groups 2 and 3 or 3A, show that these miRNAs segregate samples of groups 2 and 3 or 3A into two groups and that only three samples of group 3A (two CIS and one SQCC) were misclassified into group 2 (figs 1⇓ and 2⇓). Finally, the list of miRNAs found to be differentially expressed between CIS and invasive carcinoma perfectly discriminates the samples of these groups into two categories according to their histological classification (fig. 3⇓).
Classification of samples from groups 2 and 3, based on microRNAs differentially expressed between these groups. Dendrogram for clustering experiments, using centred correlation and average linkage. The dendrogram, based on the differentially expressed microRNAs between groups 2 (metaplasia, mild and moderate dysplasia) and 3 (severe dysplasia, carcinoma in situ (CIS) and invasive carcinoma), segregated samples of groups 2 and 3 into two groups. Only three samples of group 3 (two severe dysplasias and one CIS) were misclassified into group 2.
Classification of samples from groups 2 and 3A (A), based on microRNAs differentially expressed between these groups. Dendrogram for clustering experiments, using centred correlation and average linkage. The dendrogram, based on the differentially expressed microRNAs between groups 2 (metaplasia, mild and moderate dysplasia) and 3A (severe dysplasia and carcinoma in situ (CIS)), segregated samples of groups 2 and 3A into two groups. Only three samples of group 3A (two severe dysplasias and one CIS) were misclassified into group 2.
Classification of samples between carcinoma in situ (CIS; stage 8) and invasive squamous cell carcinoma (SQCC; stage 9), based on microRNAs differentially expressed between these groups. a) Dendrogram for clustering experiments, using centred correlation and average linkage. MicroRNAs differentially expressed between CIS and SQCC segregated the samples of these stages into two groups according to their histological classification: the six samples of stage 8 and of stage 9 were well separated into two different branches. b) The log intensities are directly proportional to the cycle threshold (CT), and the green (smaller CT) and the red (higher CT) represent up- and downregulated microRNAs, respectively. The first five microRNAs were upregulated in SQCC (the first six samples) compared with CIS (the last six samples). The last five microRNAs were downregulated from CIS to SQCC.
DISCUSSION
The results of the present study showed that several miRNAs were differentially expressed in the course of bronchial carcinogenesis and were very efficient tools to discriminate between the samples, both between low-grade lesions (metaplasia, mild and moderate dysplasia) and high-grade pre-neoplastic ones (severe dysplasia and CIS), and between in situ and invasive carcinomas. Moreover, two successive steps in miRNA expression evolution were found. Indeed, at the earliest stages of bronchial carcinogenesis, corresponding to the progression from normal epithelium of nonsmokers to phenotypically normal epithelium from smokers, to hyperplasia and until the group of relatively benign morphological abnormal bronchial tissue of smokers (metaplasia and mild and moderate dysplasia), a significant reduction of the expression of the vast majority of miRNAs was observed. However, during the later stages of carcinogenesis (severe dysplasia, CIS and SQCC), even if 74% of the altered miRNAs also remained downregulated compared with their level of expression in normal bronchial tissue of nonsmokers, the proportion of upregulated miRNAs (43%) increased compared with their level of expression in the previous group of lesions (metaplasia, mild and moderate dysplasia) and more so (80%) from that latter group to the group of severe dysplasia and CIS.
Bronchial carcinogenesis thus appears to be a process with a two-step evolution. This is reflected by the global decrease of miRNA expression observed at the earliest stages of smoke-induced morphological modifications of bronchial epithelium, followed by an increase in upregulated miRNAs at the latter stages of carcinogenesis. It has been postulated that miRNA expression is highly regulated according to the cell's developmental lineage and to its differentiation stage 3, 7, 13. In agreement with the hypothesis that an increased level of miRNA expression induces cell and tissue differentiation, miRNA downregulation often occurs in cancer where tissue is losing its normal differentiation 7, 15. Interestingly, the present data illustrate this concept. Globally, in evolution from groups 0 to 2, corresponding to morphological modifications following the loss of the normal differentiation of the normal ciliated and pseudostratified bronchial epithelium that is replaced by a keratinised and metaplastic one, the large majority of miRNAs were downregulated. Among them, the expression of several miRNAs presented an inverse evolution of their level of expression, compared with their evolution during lung embryogenesis. In particular, miR-34c and miR-15a, which are upregulated during lung development 16, were progressively downregulated during lung squamous carcinogenesis. MiR-99a, miR-142-3p and miR-142-5p, also upregulated in lung development 16, were first downregulated at the earliest steps of carcinogenesis and secondly upregulated at later ones. MiR-214 and miR-301, downregulated in lung embryogenesis 16, were upregulated at stages of severe dysplasia, CIS and SQCC.
Several miRNAs enlisted in the present study have been previously reported in lung cancer miRNA studies. However, none of the five miRNAs recently reported as prognostic factors by Yu et al. 10 was detected in the current study. Ethnic differences, the study by Yu et al. 10 involving exclusively Chinese patients, perhaps explain these discrepancies. Two miRNAs upregulated in the current study of lung carcinogenesis, miR-224 and miR-214, and one downregulated, miR-32, were reported in a study by Yanaihara et al. 9, but inversely for miR-224, which was downregulated at invasive stages 9, perhaps in accordance with the two-steps model of bronchial carcinogenesis. MiR-9, downregulated in a series of lung cancers with a higher proportion of adenocarcinomas 9, was also downregulated at early stages in the current study, from normal tissue of nonsmokers to group 2, and thereafter upregulated, from the latter stages to group 3A and again to group 3B. This pattern is another example of the two-step evolution observed for miRNA expression levels during bronchial carcinogenesis. MiR-203, the most significantly upregulated miRNA in group 2 compared with both normal bronchial tissues of nonsmokers and of smokers, was previously reported to be increased in lung cancer A549 cell lines 17 and in patients 9. Finally, miR-17-5p, upregulated at the level of group 3A, belongs to the miR-17-92 oncogenenic polycistron 18, which is overexpressed in cancer and, in particular, in lung cancer 19, where it enhances cell proliferation and inhibits cell differentiation.
Some miRNAs identified in the present study have been described as related to known oncogenes or tumour suppressor genes, and their expression level evolved through all successive stages from normal bronchi of nonsmokers to invasive SQCC in smokers. In particular, the expression level of miR-34c progressively decreased from normal epithelium of nonsmokers to invasive bronchial lesions of smokers. This miRNA was recently identified as a transcriptional target for the tumour suppressor p53 4. The progressive downregulation of miR-34c therefore reflects the progressive loss of the protective activity of p53 and is in accordance with previous studies reporting p53 alterations starting in histologically normal tissue and hyperplasia of smokers 20. MiR-15a was also progressively downregulated. First described in chronic lymphocytic leukaemia 21, miR-15a has tumour suppressor properties through the inhibition of the anti-apoptotic Bcl2 protein, which is involved in lung carcinogenesis 5.
Other miRNAs enlisted in the current study have not yet been reported to be related to lung cancer but have been described in other cancers: miR-196a in pancreatic carcinogenesis 22; miR-487b in Burkitt lymphoma 23; miR-204 as an inhibitor of apoptosis and reported in acute lymphocytic leukaemias 24; miR-106b in hepatocellular carcinomas 25; and miR-199a in ovarian cancer 26 and in acute myeloid leukaemia, where it is a signature for a bad prognosis 27. Conversely, several miRNAs identified in the present study, such as miR-452, miR-196b, miR-375, miR-193a and miR-139, have not been reported to be associated with any type of cancer until now, and their function needs to be determined.
Finally, the present study of the evolution of miRNA expression in bronchial pre-invasive lesions aimed to identify early biomarkers, especially ones detectable by noninvasive methods (serum, sputum, etc.) in high-risk patients, for early detection of lung cancer. To the current authors' knowledge, no tests or biomarkers are routinely available to solve this problem. The only attempt to identify a molecular pattern for pre-neoplastic bronchial lesions using a high-throughput technique was performed by Rahman et al. 28, who used a proteomic mass spectrometry approach called MALDI MS (matrix-assisted laser desorption ionisation mass spectrometry). Invasive lung cancer was correctly classified but they could not validate their results for pre-neoplastic lesions. Potential explanations given by the authors were the lack of true normal epithelium for reference, heterogeneity of the sources of analysed tissues, technical problems related to the MALDI MS method and a small number of available pre-neoplastic lesions. Most of these problems were avoided or solved in the present study. It was not possible to assess the prognostic value of miRNA expression in the biopsies because lesions were resected and cured by endoscopic procedure when necessary. However, previously published data indicate that severe dysplasia and CIS are at higher risk of pejorative evolution than lower-grade lesions (32–87% versus 2–9%, respectively) 29–31. Analysis of the current data revealed miRNA signatures differentiating high-grade lesions (severe dysplasia and higher grades) from lower-grade ones. If confirmed by internal and external validity studies and after assessment by more adequate techniques for screening, such as sputum examination, the current observations could have potential important practical implications for the early detection of bronchial lesions requiring treatment. This miRNA signature thus provides an original and promising tool for new lung cancer diagnosis research.
In conclusion, the present study is the first to assess microRNA expression during the successive steps of bronchial carcinogenesis on human biopsies. Most microRNA expression levels evolved in two steps during bronchial carcinogenesis. Initially decreasing during the earliest morphological modifications and the aberrant differentiation of normal bronchial epithelium into a keratinised metaplastic one, they were thereafter found to increase when malignant transformation was engaged at the stage of severe dysplasia. Several microRNAs involved in human lung development showed an inverse evolution of their level of expression during bronchial carcinogenesis, such as miR-34c, a transcriptional target of p53, which is upregulated in lung formation, but was inversely and progressively downregulated during the whole process of bronchial carcinogenesis. However, the functional roles of the majority of upregulated or downregulated microRNAs found in the current study remain unknown. These questions require further experiments. Finally, the microRNA patterns of expression were able, remarkably, to segregate high-grade pre-neoplastic lesions from the lower-grade ones, and also carcinoma in situ from invasive carcinoma. The microRNAs differentially expressed between these stages are therefore potential tools for early detection of lung cancer, at pre-invasive stages, before the spreading of micrometastases.
Support statement
C. Mascaux is a Research Fellow of the FNRS (National Fund of Scientific Research, Brussels, Belgium). This study was supported by grants from the FRSM (Fund of Medical Scientific Research), from the Télévie-FNRS and from “Les Amis de l'Institut Jules Bordet” (all Brussels).
Statement of interest
None declared.
Acknowledgments
The authors would like to thank P. Pierard and J. Faber (Institut Jules Bordet and CHU Saint-Pierre, Brussels, Belgium), M. Bruyneel, S. Bensliman and O. Bauwens (all CHU Saint-Pierre), who performed fluorescence bronchoscopy; and C. Decoster (Institut Jules Bordet and CHU Saint-Pierre), J-F. Godart (no affiliation), B. Martin and Y. Cleuter (both Institut Jules Bordet) for their help in this project. They also thank all patients who participated in this study.
- Received June 4, 2008.
- Accepted October 1, 2008.
- © ERS Journals Ltd