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

miRNAs in lung cancer - Studying complex fingerprints in patient's blood cells by microarray experiments

Abstract

Background

Deregulated miRNAs are found in cancer cells and recently in blood cells of cancer patients. Due to their inherent stability miRNAs may offer themselves for blood based tumor diagnosis. Here we addressed the question whether there is a sufficient number of miRNAs deregulated in blood cells of cancer patients to be able to distinguish between cancer patients and controls.

Methods

We synthesized 866 human miRNAs and miRNA star sequences as annotated in the Sanger miRBase onto a microarray designed by febit biomed gmbh. Using the fully automated Geniom Real Time Analyzer platform, we analyzed the miRNA expression in 17 blood cell samples of patients with non-small cell lung carcinomas (NSCLC) and in 19 blood samples of healthy controls.

Results

Using t-test, we detected 27 miRNAs significantly deregulated in blood cells of lung cancer patients as compared to the controls. Some of these miRNAs were validated using qRT-PCR. To estimate the value of each deregulated miRNA, we grouped all miRNAs according to their diagnostic information that was measured by Mutual Information. Using a subset of 24 miRNAs, a radial basis function Support Vector Machine allowed for discriminating between blood cellsamples of tumor patients and controls with an accuracy of 95.4% [94.9%-95.9%], a specificity of 98.1% [97.3%-98.8%], and a sensitivity of 92.5% [91.8%-92.5%].

Conclusion

Our findings support the idea that neoplasia may lead to a deregulation of miRNA expression in blood cells of cancer patients compared to blood cells of healthy individuals. Furthermore, we provide evidence that miRNA patterns can be used to detect human cancers from blood cells.

Peer Review reports

Background

Lung cancer is the leading cause of cancer death worldwide [1]. Its five-year survival rate is among the lowest of all cancer types and is markedly correlated to the stage at the time of diagnosis [2]. Using currently existing techniques, more than two-thirds of lung cancers are diagnosed at late stages, when the relative survival rate is low [3]. This reality calls for the search of new biomarkers that are able to catch lung cancer while it is still small and locally defined.

MicroRNAs (miRNA) are a recently discovered class of small non-coding RNAs (17-24 nucleotides) [4]. Due to their function as regulators of gene expression they play a critical role both in physiological and in pathological processes, such as cancer [5–8]. This fact is also outlined by the "Human MiRNAs & Diseases" database, the most comprehensive resource on the web, containing hundreds of entries showing the deregulation of miRNAs in a manifold of human diseases [9].

There is increasing evidence that microRNAs are not only found in tissues but also in human blood cells both as free circulating nucleic acids and in mononuclear cells. A recent proof-of-principle study demonstrated miRNA expression pattern in pooled blood sera and pooled blood cells, both in healthy individuals and in cancer patients including patients with lung cancer [10]. In addition, a remarkable stability of miRNAs in human sera was recently demonstrated [10, 11]. These findings make miRNA a potential tool for a cancer diagnostics based on blood analysis. Since single biomarkers usually lack sufficient specificity and sensitivity, we set out to analyze complex miRNA expression pattern in blood cells of cancer patients. We synthesized 866 human miRNAs and miRNA star sequences as annotated in the Sanger miRBase ([12, 13], Version 12.0) on a microarray designed by febit biomed gmbh. This array combined with the fully automated Geniom Real Time Analyzer (GRTA) platform allows for measuring miRNA fingerprints and ensures a high degree of reproducibility. To identify miRNA expression pattern we analyzed the expression of 866 miRNAs in 17 blood samples of patients with non-small cell lung carcinomas and in 19 blood samples of healthy controls.

The aim of our study was to address the following questions: Is there a larger number of differentially regulated miRNAs in blood cells of lung cancer patients as compared to healthy controls? To what extend do miRNA expression profiles in blood cells allow for the discrimination of lung cancer patients from controls? What is the information content of single miRNAs for such discrimination? Does the platform used in these experiments offer the option of a highly reproducible and efficient large-scale diagnostic test? The answers to these questions will also lay the ground for the analysis of blood based miRNA expression profiles in other cancers.

Methods

Samples

The analysis of blood from lung cancer patients and healthy subjects has been approved by local ethics committee and participants have given written informed consent. Blood samples were obtained with patients' informed consent. The patient samples stem from 17 patients with non-small cell lung carcinoma and normal controls. Normal samples were obtained from 19 different volunteers. More detailed information of patients and controls is given in Tables 1 and 2.

Table 1 Information on lung cancer patients and healthy control subjects
Table 2 Staging of patients and tumor type

miRNA extraction and microarray screening

Blood of lung cancer patients and volunteers without known disease was extracted in PAXgene Blood RNA tubes (BD, Franklin Lakes, New Jersey USA). For each blood donor, 5 ml of peripheral blood were obtained. The content of blood sampling tubes of each patient and healthy donor was centrifuged at 5000 g for 10 min at room temperature. The supernatant was discarded while the pellet was resuspended in 10 ml RNase free water by vortex. The suspension is again centrifuged at 5000 g for 10 min at RT. Subsequently, the total RNA including the miRNA was isolated using the miRNeasy kit (Qiagen GmbH, Hilden). Therefore the pellet was resuspended in 700 μl QIAzol lysis reagent and incubated for 5 min at RT. Then, 140 μl chloroform were added, vortexed for 15 sec, and incubated for 2-3 min at RT. This was followed by a centrifugation at 14000 rpm and 4°C for 15 min. The upper, watery phase was removed and 1,5 times of its volume added in 100% ethanol (ca. 900 μl). Each 700 μl of this mixture were placed on a column and centrifuged at 13000 rpm at RT for 15 sec. After the mixture had completely passed the column, 700 μl of Buffer RWT were added to each column, and again centrifuged at 13000 rpm at RT for 15 sec. Following this, 500 μl Buffer RPL were added to the column and centrifuged at 13000 rpm at RT for 15 sec. After this, another 500 μl Buffer RPL were added to the column and centrifuged at 13000 rpm at RT for 2 min. The column was then centrifuged at 13000 rpm and RT for 1 min to dry it. The following was the elution step with 40 μl RNase free water, centrifuging at 13000 rpm at RT for 1 min. The eluted RNA was stored at -70°C.

Samples were analyzed with the Geniom Realtime Analyzer (GRTA, febit gmbh, Heidelberg, Germany) using the Geniom Biochip miRNA homo sapiens. Each array contains 7 replicates of 866 miRNAs and miRNA star sequences as annotated in the Sanger miRBase 12.0 [13, 14]. Sample labeling with Biotine has been carried out either by using the miRVANAâ„¢ miRNA Labeling Kit (Applied Biosystems Inc, Foster City, California USA) or by microfluidic-based enzymatic on-chip labeling of miRNAs (MPEA [15]).

Following hybridization for 16 hours at 42°C the biochip was washed automatically and a program for signal enhancement was processed with the GRTA. The resulting detection pictures were evaluated using the Geniom Wizard Software. For each array, the median signal intensity was extracted from the raw data file such that for each miRNA seven intensity values have been calculated corresponding to each replicate copy of miRBase on the array. Following background correction, the seven replicate intensity values of each miRNA were summarized by their median value. To normalize the data across different arrays, quantile normalization [16] was applied and all further analyses were carried out using the normalized and background subtracted intensity values.

We stored all microarray data in the miRDBXP (manuscript in preparation), a database which is designed to store any type of microRNA expression pattern. The database is freely accessible at http://64.119.137.93/fmi/iwp. In addition, the data are also available in GEO (GSE17681).

Statistical analysis

After having verified the normal distribution of the measured data, we carried out parametric t-test (unpaired, two-tailed) for each miRNA separately, to detect miRNAs that show a different behavior in different groups of blood donors. The resulting p-values were adjusted for multiple testing by Benjamini-Hochberg [17, 18] adjustment. Moreover, the Mutual Information (MI) [19] was computed as a measure to access the diagnostic value of single miRNA biomarkers. To this end, all biomarkers were transformed to z-scores and binned in three bins before the MI values of each biomarker, and the information whether the marker has been measured from a normal or lung cancer sample, was computed. In addition to the single biomarker analysis classification of samples using miRNA patterns was carried out using Support Vector Machines (SVM, [20]) as implemented in the R [21] e1071 package. In detail, different kernel (linear, polynomial, sigmoid, radial basis function) Support Vector Machines have been evaluated, where the cost parameter was sampled from 0.01 to 10 in decimal powers. The measured miRNA profiles were classified using 100 repetitions of standard 10-fold cross-validation. As a subset selection technique we applied a filter approach based on t-test. In detail, the s miRNAs with lowest p-values were computed on the training set in each fold of the cross validation, where s was sampled from 1 to 866. The respective subset was used to train the SVM and to carry out the prediction of the test samples. As result, the mean accuracy, specificity, and sensitivity were calculated together with the 95% Confidence Intervals (95% CI) for each subset size. To check for overtraining we applied permutation tests. Here we sampled the class labels randomly and carried out classifications using the permuted class labels. All statistical analyzes were performed using R [21–23].

Results

miRNA experiments

We analyzed the expression of 866 miRNAs and miRNA star sequences in blood cells of 17 patients with NSCLC. As a control we used blood cells of 19 volunteers without known disease (see also Materials and Methods).

Following RNA isolation and labeling by miRVANAâ„¢ miRNA Labeling Kit, the miRNA expression profiles were measured by the Geniom Bioship miRNA homo sapiens in the GRTA (febit gmbh, Heidelberg). Following intensity value computation and quantile normalization of the miRNA profiles [16], we determined a mean correlation value of 0.97 for technical replicates by using purchased total RNA from Ambion (four heart and four liver replicates). As comparison, we computed the mean correlation of all cancer samples with all other cancer samples and all normal samples with all other normal samples. Here, we still reached a high correlation of 0.87 and a variance of 0.009.

ruling out the influence of age and gender

To cross-check that age and gender do not have an influence on our analysis, we computed t-tests for the normal samples. In the case of males versus females we did not find any statistically significant deregulated miRNA. The most significant miRNA, hsa-miR-423, showed an adjusted significance level of 0.78.

To test for the influence of donor age we compared the profiles obtained from samples obtains from the oldest versus youngest patients by splitting the group in half based on age. Here the most significant miRNA, miR-890, obtained an adjusted p-value of 0.87. As for gender, we did not find any deregulated miRNAs, thus providing evidence that age and gender do not have a substantial influence on the miRNA profiles.

single deregulated miRNAs

We applied hypothesis testing to identify miRNAs deregulated in the blood cells of lung cancer patients as compared to the blood cells of the controls. Following verification of an approximately normal distribution, we performed two-tailed unpaired t-tests for each miRNA. The respective p-values were adjusted for multiple testing by the Benjamini-Hochberg approach [17, 18]. In total we detected 27 miRNAs significantly deregulated in blood cells of lung cancer patients as compared to the controls. A complete list of deregulated miRNAs is given in the Supplemental Material (see Additional file 1 - Table S1). The miRNAs that were most significantly deregulated included hsa-miR-126 with a p-value of 0.00003, hsa-let-7d with a p-value of 0.003, hsa-let-7i with a p-value of 0.003, and hsa-miR-423 with a p-value of 0.001 (Figure 1 and Figure 2). Other members of the let-7 family that were also found to be deregulated included hsa-let-7c, hsa-let-7e, hsa-let-7f, hsa-let-7g and hsa-let-7a. Besides miR-423, all above mentioned miRNAs were down-regulated in blood cells of lung cancer patients compared to blood cells of healthy subjects indicating an overall decreased miRNA repertoire.

Figure 1
figure 1

The bar-chart shows for 15 of the 27 deregulated miRNAs the median value of cancer samples and normal samples. Here, blue bars correspond to cancer samples while red bars to controls.

Figure 2
figure 2

Back to back histograms of two examples. The color-coding corresponds to the one used in Figure 1.

To validate our findings, we repeated the miRNA profiling using an enzymatic on-chip labeling technique termed MPEA (Microfluidic-based enzymatic on-chip labeling of miRNAs) [24]. For this control experiment, we used 4 out of the 17 lung cancer patients and 10 of the controls. Hereby, we detected 100 differentially regulated miRNAs. The miRNAs that were most significantly deregulated include hsa-miR-1253 with a p-value of 0.001, hsa-miR-126 with a p-value of 0.006, hsa-let-7d with a p-value of 0.006, and hsa-let-7f with a p-value of 0.006. Of the previously identified 27 miRNAs 12 were detected to be significant in the second experiment, while the remaining miRNAs showed increased p-values. The correlation of fold changes was 0.62. We also confirmed other members of the let-7 family as deregulated in blood cells of lung cancer patients. Furthermore, we confirmed that the majority of the deregulated miRNAs were down-regulated in patients' blood cells. In these Experiments, 62% of the significantly deregulated miRNAs showed decreased intensity values in lung cancer samples. In the Supplemental Material we provide a complete list of deregulated miRNAs identified by MPEA (see Additional file 1 - Table S1).

As a further control experiment we performed an expression analysis by qRT-PCR. As a test sample we analyzed the fold changes of hsa-miR-106b, miR-162, miR-98, miR-let7d, miR-22, and miR-140-3p in blood cells of eight tumor patients and five controls. Of these 6 miRNAs, 3 belong to the group of 12 miRNAs significant for miRVANA labeling and the MPEA method. The result of this validation experiment was a verification of the microarray screening. For example, in agreement with the up-regulation detected by the Geniom Biochip experiments, the qRT-PCR experiments showed on average an up-regulation of 1.4 fold in blood cells of lung cancer patients as compared to healthy controls for miR 106b. The excellent correlation of microarray and qRT-PCR data is demonstrated by an R2 value of 0.994 (see Figure 3).

Figure 3
figure 3

Scatterplot of fold quotients of qRT-PCR (x-axis) and microarray experiments (y-axis).

Diagnostic value of miRNA biomarkers

Mutual Information (MI) [19] is an adequate measure to estimate the overall diagnostic information content of single biomarkers [22]. In our study, Mutual Information is considered as the reduction in uncertainty about the class labels '0' for controls and '1' for tumor samples due to the knowledge of the miRNA expression. The higher the value of the MI of a miRNA, the higher is the diagnostic content of the respective miRNA.

We computed the MI of each miRNA with the class labels. First, we carried out a permutation test to determine the background noise of the miRNAs, e.g. the random information content of each miRNA. We randomly selected 1000 miRNAs (with replacements) and sampled the class labels for each miRNA. These permutation tests yielded a mean MI value of 0.029, a 95% quantile of 0.096 and a value of 0.217 for the highest random MI. Second, we calculated the MI values for the comparison between the miRNAs in blood cells of tumor patients and controls. The overall comparison of the 866 miRNAs yielded significantly increased MI values with a two-tailed p-value of ≤ 10-10 as shown by an unpaired Wilcoxon Mann-Whitney test [23, 24]. The miRNA hsa-miR-361-5p showed the highest MI with a value of 0.446. The miRNAs with the best significance values as computed by the t-test, namely hsa-miR-126 and hsa-miR-98, were also among the miRNAs showing the highest MI values. In total we detected 37 miRNAs with MI values higher than the highest of 1000 permuted miRNAs and 200 miRNAs with MI values higher than the 95% quantile (Figure 4). Of the initially detected 27 miRNAs by the t-test, 14 had values higher than the highest permutation test and all remaining 13 miRNAs belong to the group of miRNAs with MI values higher than the 95% quantile, demonstrating the overall well agreement of the two measures. Venn diagrams for the different analyses are provided as Venn diagrams (Figures 5A and 5B). A complete list of miRNAs, the respective MI and the enrichment compared to the background MI is provided in the supplemental material (see Additional file 2 - Table S2).

Figure 4
figure 4

The mutual information of all miRNAs that have higher information content than the best permutation test (upper red line). The middle red line denotes the 95% quantile of the 1000 permutation tests and the bottom red line the mean of the permutation experiments, corresponding to the background MI.

Figure 5
figure 5

Venn diagrams for different analyses. A) Venn diagram of the t-test mirVANA and MPEA labeling. B) Venn diagram of the t-test and the Mutual Information. The right circle shows significant MI miRNAs, separated in the miRNAs that are higher than the highest permutation test and the miRNAs higher than 95% of all 1000 permutation tests. The 27 miRNAs are all included in these miRNAs, indicated by the green circle. C) Venn diagram of t-test miRNAs and miRNAs used in the best classification subset. The 24 miRNAs used for classification are all contained in the set of 27 most significant mRNAs

Evaluating complex fingerprints

Even single miRNAs with highest MI values are not sufficient to differentiate between blood cells of tumor patients as compared to controls with high specificity. For example, the hsa-miR-126 separates blood cells of tumor patients from blood cells of healthy individuals with a specificity of 68%, only. In order to improve the classification accuracy we combined the predictive power of multiple miRNAs by using statistical learning techniques. In detail, we applied Support Vector Machines with different kernels (linear, polynomial, sigmoid, radial basis function) to the data and carried out a hypothesis test based subset selection as described in Material and Methods. To gain statistical significance we carried out 100 repetitions of 10-fold cross validation. Likewise, we computed 100 repetitions for the permutation tests.

The best results were obtained with radial basis function Support Vector Machines and a subset of 24 miRNAs (Table 3 and Venn Diagram in Figure 5C). These miRNAs allowed for the discrimination between blood cells of lung tumor patients and blood cells of controls with an accuracy of 95.4% [94.9%-95.9%], a specificity of 98.1% [97.3%-98.8%], and a sensitivity of 92.5% [91.8%-92.5%]. The permutation tests showed significantly decreased accuracy, specificity, and sensitivity with 49.2% [47.2%-51.3%], 56.9% [54.5%-59.3%] and 40.6% [37.9%-43.4%], respectively (Figure 6), providing evidence that the obtained results are not due to an overfit of the statistical model on the miRNA fingerprints.

Table 3 miRNAs that are in the best subset
Figure 6
figure 6

Boxplots of the classification accuracy, specificity, and sensitivity of the 100 repetitions of the 10-fold cross-validation (red boxes) and the same values for random permutation tests (blue boxes).

Discussion

While complex miRNA expression patterns have been reported for a huge variety of human tumors, far less information is available on miRNA expression patterns of blood cells derived from tumor patients. In the following we relate our miRNA expression profiling to both the miRNA expression in blood cells and in cancer cells of non-small cell lung cancer patients. We found significant down-regulation of hsa-miR-126 that was recently detected in blood cells of healthy individuals, but not in blood cells of lung cancer patients [10]. Down-regulation of hsa-miR-126 was also found in lung cancer tissue [25]. Functional studies on hsa-miR-126 revealed this miRNA as a regulator of the endothelial expression of vascular cell adhesion molecule 1 (VCAM-1), which is an intercellular adhesion molecule expressed by endothelial cells [26]. hsa-miR-126 is also reported to be an inhibitor of cell invasion in non-small cell lung cancer cell lines, and down-regulation of this miRNA might be a mechanism of lung cancer cells to evade these inhibitory effects [27]. Members of the hsa-let-7 family that were found down-regulated in our study were the first miRNAs reported as de-regulated in lung cancer [28]. This down-regulation of the let-7 family in lung cancer was confirmed by several independent studies [29–33]. Our data are also in agreement with a recent study showing the down-regulation of hsa-let-7a, hsa-let-7d, hsa-let-7f, hsa-let-7g, and hsa-let-7i in blood cells of lung cancer patients [10]. Notably, down-regulation of let-7 in lung cancer was strongly associated with poor clinical outcome [29]. The let-7 family members negatively regulate oncogene RAS [28]. The miRNA hsa-miR-22 that showed a high MI value and up-regulation in our study, was recently also reported to be up-regulated in blood cells of lung cancer patients [10]. The miRNA hsa-miR-19a that also showed a high MI value and up-regulation in our study was reported to be up-regulated in lung cancer tissue [34, 35]. In contrast, hsa-miR-20a, which is significantly down-regulated in our experiments, was reported as up-regulated in lung cancer tissue [34, 35]. The up-regulation of hsa-miR-20a was found in small-cell lung cancer cell lines, our study investigated only NSCLC. In summary, there is a high degree of consistency between miRNA expression found in the peripheral blood cells of lung cancer patients and miRNA expression in lung cancer tissue [25, 29–37].

Some of the deregulated miRNAs identified in our study are also reported as de-regulated in other cancer entities, e.g. hsa-miR-346 in gastritic cancer, hsa-miR-145 in bladder cancer, and hsa-miR-19a in hepatocellular carcinoma and B-cell leukemia [38–42]. We also found miRNAs with high diagnostic potential e.g. high MI value, that were not yet related to cancer as for example hsa-miR-527 or hsa-mir-361-5p that were both up-regulated in blood cells of lung cancer patients.

Besides the deregulation of single miRNAs, we analyzed the overall expression pattern of miRNAs in peripheral blood cells of lung cancer patients in comparison to the pattern in blood cells of healthy controls. Recently, Chen et al. [10] reported a high correlation of 0.9205 between miRNA profiles in serum and miRNA profiles in blood cells, both in healthy individuals. The correlation of the miRNA profiles between serum and blood cells in lung cancer patients were significantly lower (0.4492). These results are indicative of deregulated miRNAs in blood and/or serum of patients and are in agreement with our data that show the deregulation of miRNAs in the blood cells of lung carcinoma patients. These deregulated miRNAs can be used to differentiate patients with lung cancer from normal controls with high specificity and sensitivity.

Conclusion

This is the first evidence for the diagnostic potential of miRNA expression profiles in peripheral blood cells of cancer patients and healthy individuals. It remains to be seen whether other cancers and other non-cancer diseases also show a specific miRNA expression pattern that might be used to tell these diseases apart form controls and possibly apart from each other.

We are optimistic that the high specificity and sensitivity of the miRNAs deregulated in blood cells of lung carcinoma patients together with the high reproducibility of the applied technique will open avenues for applying this approach in prospective trials of a non-invasive diagnostic test.

Abbreviations

miRNA:

microRNA

NSCLC:

non-small cell lung cancer

qRT-PCR:

quantitative real time polymerase chain reaction

MI:

Mutual Information

GRTA:

Geniom Realtime Analyzer

MPEA:

microfluidic-based enzymatic on-chip labeling of miRNAs

SVM:

Support Vector Machines.

References

  1. Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics, 2008. CA Cancer J Clin. 2008, 58 (2): 71-96. 10.3322/CA.2007.0010.

    Article  PubMed  Google Scholar 

  2. Scott WJ, Howington J, Feigenberg S, Movsas B, Pisters K: Treatment of non-small cell lung cancer stage I and stage II: ACCP evidence-based clinical practice guidelines (2nd edition). Chest. 2007, 132 (3 Suppl): 234S-242S. 10.1378/chest.07-1378.

    Article  PubMed  Google Scholar 

  3. Henschke CI, Yankelevitz DF: CT screening for lung cancer: update 2007. Oncologist. 2008, 13 (1): 65-78. 10.1634/theoncologist.2007-0153.

    Article  PubMed  Google Scholar 

  4. Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.

    Article  CAS  PubMed  Google Scholar 

  5. Calin GA, Croce CM: MicroRNA signatures in human cancers. Nat Rev Cancer. 2006, 6 (11): 857-866. 10.1038/nrc1997.

    Article  CAS  PubMed  Google Scholar 

  6. Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006, 6 (4): 259-269. 10.1038/nrc1840.

    Article  CAS  PubMed  Google Scholar 

  7. Sassen S, Miska EA, Caldas C: MicroRNA: implications for cancer. Virchows Arch. 2008, 452 (1): 1-10. 10.1007/s00428-007-0532-2.

    Article  CAS  PubMed  Google Scholar 

  8. Zhang B, Pan X, Cobb GP, Anderson TA: microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007, 302 (1): 1-12. 10.1016/j.ydbio.2006.08.028.

    Article  CAS  PubMed  Google Scholar 

  9. Lu M, Zhang Q, Deng M, Miao J, Guo Y, Gao W, Cui Q: An analysis of human microRNA and disease associations. PLoS ONE. 2008, 3 (10): e3420-10.1371/journal.pone.0003420.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, Guo J, Zhang Y, Chen J, Guo X, et al: Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008, 18 (10): 997-1006. 10.1038/cr.2008.282.

    Article  CAS  PubMed  Google Scholar 

  11. Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, Benjamin H, Kushnir M, Cholakh H, Melamed N, et al: Serum microRNAs are promising novel biomarkers. PLoS ONE. 2008, 3 (9): e3148-10.1371/journal.pone.0003148.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, D140-144. 10.1093/nar/gkj112. 34 Database

  13. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-158. 36 Database

  14. Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Eddy SR, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, D121-124. 33 Database

  15. Vorwerk S, Ganter K, Cheng Y, Hoheisel J, Stahler PF, Beier M: Microfluidic-based enzymatic on-chip labeling of miRNAs. N Biotechnol. 2008, 25 (2-3): 142-149. 10.1016/j.nbt.2008.08.005.

    Article  CAS  PubMed  Google Scholar 

  16. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.

    Article  CAS  PubMed  Google Scholar 

  17. Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.

    Google Scholar 

  18. Hochberg Y: A sharper bonferroni procedure for multiple tests of significance. Biometrica. 1988, 75: 185-193. 10.1093/biomet/75.4.800.

    Article  Google Scholar 

  19. Shannon C: A mathematical theory of communication. The Bell System Technical Journal. 1984, 27: 623-656.

    Article  Google Scholar 

  20. Vapnik V: The Nature of Statistical Learning Theory. 2000, Springer, ISBN 978-0-387-98780-4, Second

    Chapter  Google Scholar 

  21. Team RDC: R: A Language and Environment for Statistical Computing. 2008, Vienna, Austria: R Foundation for Statistical Computing

    Google Scholar 

  22. Keller A, Ludwig N, Comtesse N, Hildebrandt A, Meese E, Lenhof HP: A minimally invasive multiple marker approach allows highly efficient detection of meningioma tumors. BMC Bioinformatics. 2006, 7: 539-10.1186/1471-2105-7-539.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Wilcoxon F: Individual comparisons by ranking methods. Biometric Bull. 1945, 1: 80-83. 10.2307/3001968.

    Article  Google Scholar 

  24. Mann H, Wilcoxon F: On a test whether one of two random variables is stochastically larger than the other. Ann Mat Stat. 1947, 18: 50-60. 10.1214/aoms/1177730491.

    Article  Google Scholar 

  25. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, et al: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006, 9 (3): 189-198. 10.1016/j.ccr.2006.01.025.

    Article  CAS  PubMed  Google Scholar 

  26. Harris TA, Yamakuchi M, Ferlito M, Mendell JT, Lowenstein CJ: MicroRNA-126 regulates endothelial expression of vascular cell adhesion molecule 1. Proc Natl Acad Sci USA. 2008, 105 (5): 1516-1521. 10.1073/pnas.0707493105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Crawford M, Brawner E, Batte K, Yu L, Hunter MG, Otterson GA, Nuovo G, Marsh CB, Nana-Sinkam SP: MicroRNA-126 inhibits invasion in non-small cell lung carcinoma cell lines. Biochem Biophys Res Commun. 2008, 373 (4): 607-612. 10.1016/j.bbrc.2008.06.090.

    Article  CAS  PubMed  Google Scholar 

  28. Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120 (5): 635-647. 10.1016/j.cell.2005.01.014.

    Article  CAS  PubMed  Google Scholar 

  29. Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, et al: Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004, 64 (11): 3753-3756. 10.1158/0008-5472.CAN-04-0637.

    Article  CAS  PubMed  Google Scholar 

  30. Tong AW: Small RNAs and non-small cell lung cancer. Curr Mol Med. 2006, 6 (3): 339-349. 10.2174/156652406776894554.

    Article  CAS  PubMed  Google Scholar 

  31. Stahlhut Espinosa CE, Slack FJ: The role of microRNAs in cancer. Yale J Biol Med. 2006, 79 (3-4): 131-140.

    PubMed  Google Scholar 

  32. Williams AE: Functional aspects of animal microRNAs. Cell Mol Life Sci. 2008, 65 (4): 545-562. 10.1007/s00018-007-7355-9.

    Article  CAS  PubMed  Google Scholar 

  33. Zhang B, Wang Q, Pan X: MicroRNAs and their regulatory roles in animals and plants. J Cell Physiol. 2007, 210 (2): 279-289. 10.1002/jcp.20869.

    Article  CAS  PubMed  Google Scholar 

  34. Calin GA, Croce CM: MicroRNA-cancer connection: the beginning of a new tale. Cancer Res. 2006, 66 (15): 7390-7394. 10.1158/0008-5472.CAN-06-0800.

    Article  CAS  PubMed  Google Scholar 

  35. Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65 (21): 9628-9632. 10.1158/0008-5472.CAN-05-2352.

    Article  CAS  PubMed  Google Scholar 

  36. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al: MicroRNA expression profiles classify human cancers. Nature. 2005, 435 (7043): 834-838. 10.1038/nature03702.

    Article  CAS  PubMed  Google Scholar 

  37. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, et al: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006, 103 (7): 2257-2261. 10.1073/pnas.0510565103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Guo L, Huang ZX, Chen XW, Deng QK, Yan W, Zhou MJ, Ou CS, Ding ZH: Differential Expression Profiles of microRNAs in NIH3T3 Cells in Response to UVB Irradiation. Photochem Photobiol. 2009, 85 (3): 765-73. 10.1111/j.1751-1097.2008.00482.x.

    Article  CAS  PubMed  Google Scholar 

  39. Ichimi T, Enokida H, Okuno Y, Kunimoto R, Chiyomaru T, Kawamoto K, Kawahara K, Toki K, Kawakami K, Nishiyama K, et al: Identification of novel microRNA targets based on microRNA signatures in bladder cancer. Int J Cancer. 2009, 125 (2): 345-52. 10.1002/ijc.24390.

    Article  CAS  PubMed  Google Scholar 

  40. He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, et al: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435 (7043): 828-833. 10.1038/nature03552.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Alvarez-Garcia I, Miska EA: MicroRNA functions in animal development and human disease. Development. 2005, 132 (21): 4653-4662. 10.1242/dev.02073.

    Article  CAS  PubMed  Google Scholar 

  42. Feitelson MA, Lee J: Hepatitis B virus integration, fragile sites, and hepatocarcinogenesis. Cancer Lett. 2007, 252 (2): 157-170. 10.1016/j.canlet.2006.11.010.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We want to thank the healthy individuals and lung cancer patients for participating in this study. This Work was supported by funding of the German Ministry of Research Education (BMBF) under contract 01EX0806 and by Hedwig-Stalter foundation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Keller.

Additional information

Competing interests

AK: salary (febit), patent (febit); AB salary (febit), patent (febit); MS salary (febit); FW salary (febit); AW salary (febit);

Authors' contributions

AK contributed to study design, carried out the statistical analyses, initiated the study and contributed in paper writing. PL developed the protocol for miRNA extraction, extracted the miRNAs from blood cells, and contributed in paper writing. AB performed the miRNA screening. AW added the data into the miRNA database and contributed in the statistical evaluation. FW assisted the evaluation of the miRNA microarrays. MS contributed to study design and writing of the paper. HH contributed in study design, collected the blood samples and contributed in paper writing. HPL contributed in the bio-statistical analysis. EM contributed in study design, development of miRNA extraction, supervised the study and contributed in paper writing. All authors read and approved the final manuscript.

Andreas Keller, Petra Leidinger contributed equally to this work.

Electronic supplementary material

12885_2009_1687_MOESM1_ESM.XLS

Additional file 1: Table S1. Details of all miRNAs for the miRVANA screening (columns A-I), MPEA Assay (columns K-S) and the MI (columns U-V). Significant miRNAs are colored. (XLS 352 KB)

12885_2009_1687_MOESM2_ESM.XLS

Additional file 2: Table S2. Mutual Information of biomarkers. This table contains for each miRNA biomarker the Mutual Information (MI) together with the enrichment as compared to the background MI. The miRNAs with MI higher than the highest background MI are colored red, the miRNAs with MI values higher than highest 5% background MI's in blue. (XLS 140 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Keller, A., Leidinger, P., Borries, A. et al. miRNAs in lung cancer - Studying complex fingerprints in patient's blood cells by microarray experiments. BMC Cancer 9, 353 (2009). https://doi.org/10.1186/1471-2407-9-353

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-9-353

Keywords