- Primary Research
- Open Access
Role of microRNA-33a in regulating the expression of PD-1 in lung adenocarcinoma
Cancer Cell International volume 17, Article number: 105 (2017)
MiRNAs are vital in functioning as either oncogenes or tumor suppressors in the cell cycle. Target transcripts for immune checkpoint molecules such as PD-1/PD-L1 and (programmed cell death-1/its ligand and cytotoxic T-lymphocyte antigen 4) have proven to be beneficial against several solid tumors, including lung adenocarcinoma.
Simultaneous quantification of the expression level of miR-33a and PD-1, PD-L1 and CTLA4 mRNAs with NanoString technology was performed in 88 lung adenocarcinoma specimens. A cohort of 323 lung adenocarcinoma patients from the cancer genome atlas (TCGA) database was further analyzed, in order to test our hypothesis. Potential interference of PD-1, PD-L1 and CTLA4 gene expression by miR-33a was predicted using the microRNA target prediction program RNA22.
High miR-33a expression was significantly associated with younger (p = 0.005), female (p = 0.04), patients with low grade (p < 0.0001), early stage (p = 0.03) tumors, and better survival. The hypothesis of the involvement of miR-33a in PD-1/PD-L1/CTLA4 mechanisms was corroborated by the finding of putative miR-33a binding sites in all three genes using the RNA22 method. We found an inverse correlation between miR-33a and PD-1 levels (p = 0.01), as well as for PD-L1 (p = 0.01) and CTLA4 (p = 0.03) expression, and a significant better prognosis for patients with high miR-33a/low PD-1. TCGA database analysis confirmed that miR-33a high levels were associated with low PD-1 expression and with longer survival on a larger population.
Our study emphasizes the notion of a potential value of miR-33a as a favorable prognostic marker through PD-1 regulation.
MicroRNAs (miRNAs) are small non-coding RNA molecules that function as indispensable regulators of an increasing number of cellular processes [1,2,3]. miRNAs are vital in regulating cell proliferation and apoptosis, in addition to functioning as either oncogenes or tumor suppressors in the cell cycle . In a recent study , miR–33a was identified as having potential tumor–suppressive activity. Overexpression of miR–33a was demonstrated to inhibit the growth of lung cancer cells, but the exact role of individual miRNAs strictly depends on their expression pattern and their targeted genes. Lung cancer is the most common cancer worldwide and despite recent progress with molecularly targeted agents, its prognosis is usually poor and new strategies need to be developed. Cancer immunotherapy is emerging as a very promising therapeutic strategy and recently various clinical trials exploring the use of anti-programmed cell death protein 1 (PD-1, also known as CD279), PD-L1 (programmed cell death-1 ligand, also known as CD274) and CTLA4 (cytotoxic T-lymphocyte antigen 4) inhibitors have successfully shown antitumor activities in lung cancer [6,7,8]. For most cancer biomarkers, determination of the levels of such immune function markers in formalin-fixed, paraffin-embedded (FFPE) samples has been generally performed by conventional immunohistochemistry using various antibodies with various levels of validation. Different thresholds have been used to define PD-1/PDL1 positivity, but cut-off points as the proportion of membrane-positive tumor cells were subjective in terms of estimated visual levels, such that reproducibility has not been formally assessed in the clinical setting.
Moreover, there are limited data on the prognostic significance of PD-1/PD-L1, with some showing poor , better prognosis  or no prognostic role . Another approach to regulating the immune response in the tumor microenvironment is by modulating the level of miRNAs. A recent review  emphasized the role of miRNAs as modulators of immune checkpoint molecules and consequently as potential cancer therapeutic targets.
Here, we investigated the potential role of miR-33a and demonstrated its potential value as prognostic marker by regulating PD-1/PD-L1 and CTLA4 expression in lung adenocarcinoma.
Materials and methods
Eighty-eight lung adenocarcinoma patients were retrospectively selected among those operated at the Unit of Thoracic Surgery of the A.O.U.P. Histological diagnoses were independently formulated according to the World Health Organization classification [13,14,15]. Clinical–pathological characteristics were collected whenever available for all the patients.
Total RNA were isolated from a representative area selected and marked on the surface of 10 micron sections of formalin-fixed, paraffin-embedded (FFPE) tissues using the miRNeasy FFPE Kit (Qiagen Inc., Hilden, Germany) according to the manufacturer’s instructions. The quality and concentration of RNA was assessed using a NanoDrop spectrophotometer (Thermo-Scientific, Wilmington, Del).
NanoString nCounter® assay, data normalization and analysis
Expression of the selected miR-33a, PD-1, PD-L1 and CTLA4 genes among a miRGE panel was measured using the NanoString nCounter technology (NanoString Technologies, Seattle, WA). The nCounter measures total counts of mRNAs/miRNA by a multiplex hybridization assay followed by scanning and digital readouts of fluorescent probes . Raw NanoString counts for each gene were subjected to a technical normalization, considering the counts obtained for positive control probe sets, followed by a biological normalization, using three reference genes included in the CodeSet, according to NanoString Technologies’ guidelines. The nCounter custom code set used in this study included three genes (PD-1, PD-L1 and CTLA4) with five housekeeping genes for reference [clathrin heavy chain 1 (CLTC), beta-glucuronidase (GUSB), tubulin beta (TUBB), hypoxanthine phosphoribosyltransferase (HPRT1), phosphoglycerate kinase 1 (PGK1)]. MiR-33a expression was tested simultaneously with other selected miRNAs and the normalization was performed using a scaling factor based on miRNAs with the lower variability coefficients.
The Cancer Genome Atlas (TCGA) database. From the TCGA data portal (http://tcga.cancer.gov/; accessed October 2017), we extracted PD-1 and miR-33a expression together with the corresponding clinical–pathological characteristics and survival data for 323 adenocarcinoma patients (LUAD).
Potential regulation of PD-1, PD-L1 and CTLA4 gene expression by miR-33a was predicted using the microRNA target prediction program RNA22, a pattern-based approach for the discovery of microRNA binding sites and their corresponding microRNA/mRNA complexes .
Once the RNA hybridization data had been correctly prepared, the data were subjected to 2-way hierarchical clustering analysis (HCA) using nSolver 2.5 Analysis Software. Differential expression was determined by applying the non-parametric Wilcoxon test in order to determine the association between miR-33a expression and the clinical–pathological parameters. Survival analyses were performed using the Kaplan–Meier method with log-rank test and the Cox proportional hazard model. Statistical analyses were performed using JMP10 software (SAS, Milan, Italy), and a two-tailed p value < 0.05 was considered significant.
This study included 88 patients with lung adenocarcinoma, 56 males and 32 females, with an age at diagnosis ranging from 30 to 81 years (mean 58.9 years, median 54.5). The predominant histologic patterns were characterized as follows: lepidic (29/88, 33%), solid (26/88, 29.5%), acinar (22/88, 25%), and papillar (11/88, 12.5%) variants. According to the degree of differentiation, three tumors were G1, whereas 58 and 27 were G2 and G3, respectively. The adenocarcinomas were all invasive. Their stages were classified as IA (17), IB (23), IIA (13), IIB (9), IIIA (23), IIIB (1), and IV (2). The survival data, with disease-free interval (DFI) and overall survival (OS), were available for all the patients and was last updated in March 2015. Concerning the smoking habits, there were 17 non-smokers, 16 former smokers, and 23 current smokers; for 33 patients, smoking history was unknown.
miR-33a expression and clinical–pathological characteristics
miR-33a expression was compared with the patients’ clinical–pathological characteristics. High miR-33a expression was significantly associated with younger (p = 0.005), female (p = 0.04) patients with low grade (p < 0.0001), early stage (p = 0.03) tumors (Table 1).
miR-33a and survival analysis
A Kaplan–Meier survival analysis using DFI (range 0–148, median 21 months) and OS (range 3–148, median 31.5 months) as endpoints and miR-33a expression as a dichotomous variable, distinguishing tumors with low from tumors with high expression, showed a significant better prognosis in patients with high miR-33a levels (Wilcoxon test, p = 0.02 and p = 0.008, for DFI and OS, respectively) compared to those with low levels (Fig. 1).
miR-33a target prediction
To investigate the putative involvement of miR-33a in the regulation of the most studied immune checkpoint molecules (PD-1, PD-L1, and CTLA4), we used RNA22 software for the discovery of microRNA binding sites and their corresponding microRNA/mRNA complexes. The RNA22 method suggests that the number of microRNA binding sites may be greater than hypothesized and that microRNA regulation may be effected through the 5′UTR and CDS of gene transcripts in animals, in addition to 3′UTRs. Specifically, the hypothesis of the involvement of miR-33a in PD-1/PD-L1/CTLA4 mechanisms was corroborated by the observation that the 3′ untranslated region (UTR) of the PD-1 mRNA carries a putative miR-33a binding site at position 1789. Two miR-33a binding sites for CTLA4, at positions 96 and 253 of the CDS region, and one in the 3′UTR of PD-L1 (position 1701) were also found (Fig. 2).
Correlation between miR-33a and PD-1 expression
We evaluated the abundance of PD-1, as well as of PDL-1 and CTLA4, mRNA by NanoString technology. The samples were divided into high and low expression groups based on the median fold-change value (1.76-fold change value ± 7.34 for PD-1; 2.79 ± 11.5 for PDL-1, and 30.23 ± 46.08 for CTLA4). Samples with low PD-1 mRNA levels also demonstrated low protein expression, assessed in FFPE tissue samples by immunohistochemistry using the anti-PD-1 mouse monoclonal antibody NAT105 (Ventana Medical Systems, Inc. USA) (data not shown).
We then analyzed their relationship with the expression level of miR-33a and found a statistically significant inverse correlation between miR-33a and PD-1/PD-L1 expression. MiR-33a levels were lower in patients with higher PD-1 (Chi square test, p = 0.01), PD-L1 and CTLA4 expression (p = 0.01 and p = 0.03, respectively) (Fig. 3). To evaluate possible relationships between patients’ survival, miR-33a, and PD-1, we grouped lung adenocarcinoma patients according to the co-expression of both factors (miR-33a and PD-1). We then compared the difference in survival between the two groups linked to a negative association, high miR-33a/low PD-1 and low miR-33a/high PD-1 expression. We found that the first group of patients, with high miR-33a and low PD-1, showed a better prognosis, either for DFI or OS, than the opposite group (p = 0.04 and p = 0.007, respectively) (Fig. 4), implicating miR-33a as a good prognostic marker as a consequence of PD-1 regulation.
TCGA data analysis
A cohort of 323 lung adenocarcinoma patients from TCGA database was further analyzed, in order to validate our findings on a larger population. The samples were divided into high and low miR-33a expression groups based on the median value, and statistical analysis revealed that high miR-33a levels were significantly associated with low PD-1 expression (t test; p = 0.03), as we found by NanoString methodology. Moreover, survival analysis confirmed the association between miR-33a high expression and a better survival (p = 0.004 for DFI, and p = 0.007 for OS) (data not shown).
Gaining insight into the molecular basis of lung cancers is of critical clinical relevance in order to identify subgroups of patients for a more accurate management. Several studies have identified miRNA signatures with diagnostic and prognostic relevance [18, 19]. MiRNAs potentially target hundreds of different mRNAs, thus regulating a wide variety of cellular processes . miR-33a has been shown to have potential tumor-suppressive activity by downregulating the expression of beta-catenin , but an improved deeper understanding of the role of this miRNA and its alternative antitumor mechanism is needed. In the current study, we showed that high miR-33a expression was associated with several favorable clinical–pathological characteristics, such as young age, female gender, low tumor grade, and early stage. Survival analysis also confirmed the good prognostic value of miR-33a. Our analysis next focused on the putative role of miR-33a in regulating the expression of immune markers in order to account for the antitumor function of miR-33a. The PD-1 receptor is a member of the immunoglobulin CD28 family, playing a crucial role in immune escape during tumor progression . PD-L1 is expressed in different cancer types, including lung cancer, and its interaction with PD-1 plays an important role in the blockade of the “cancer immunity cycle” . One of the most promising approaches in cancer, including lung adenocarcinoma, is antibody blockade of the PD-1/PDL-1 pathway [23,24,25]. This approach raises several questions, in particular whether PD-1/PD-L1 expression status is important in order to select patients eligible for these treatments, which technique is best suited in evaluating their levels, if these tumor immunity factors may be regulated by miRNAs, and how to define a threshold for positive PD-1/PD-L1 staining of tissue samples, considering that certain patients respond to treatment targeting PD-L1/PD-1, despite low or absent immunoreactivity for these biomarkers. Previously published studies suggest that tumor PD-L1 protein expression may be evaluated on human cancers using immunohistochemistry (IHC) in FFPE samples [6, 26]; however, several factors may impact the evaluation of PD-L1 positivity by conventional IHC, such as the specificity and reproducibility of the commercially available antibodies as well as the subjectivity of the staining interpretation. Velcheti et al.  reported a novel method of in situ measurement of PD-L1 mRNA NanoString nCounter technology, suggesting its utility for the accurate measurement of PD-1/PDL-1 levels in predicting the response of lung adenocarcinoma patients to targeted immunotherapy. NanoString employs unamplified nucleic acids and no cDNAs or enzymatic reaction step and is, therefore, less sensitive to tissue fixation effects; in our experience, this platform was extremely suitable to operate at very low levels of expression. Our findings provide a proof-of-principle that transcriptome-based analysis is a highly sensitive and effective means for PD-1/PD-L1 testing. Moreover, further analysis on a cohort of 323 lung adenocarcinoma patients from the TCGA database confirmed our findings on a larger population and also using a different transcriptome-based technology, such as Illumina HiSeq quantification.
The finding that targeting the PD-1/PD-L1 mechanism is of crucial importance in lung adenocarcinoma has led to increasing search for biomarkers predictive of response. Even if Nivolumab, an antibody targeting PD-1, has recently received Food and Drug Administration (FDA) approval for NSCLC therapy, regardless of the PDL1 expression status, a better prediction of which patients are more likely to respond to this cancer therapy may improve the treatment costs. Moreover, the effects of PD-1/PD-L1 signaling in the outcome of several neoplasms, including lung cancer, are not completely understood due to contradictory results [28, 29]. An alternative strategy in regulating immune response to tumors could be the expression of miRNAs that target immune checkpoint mRNAs. Using miRNA target prediction tools, we demonstrated a direct interaction of miR-33a and the 3′UTR of PD-1 and PD-L1 mRNAs, with consequent downregulation of the expression of both genes. In most circumstances, miRNAs bind to the target mRNA at the 3′ UTR region  and few of them were also reported to modulate genes through binding with the CDS region [31,32,33]. The target genes of miR-33a were surveyed globally by examining 3′UTR and CDS. Three putative binding sites of miR-33a-5p were found in the CDS region of the CTLA4 gene with similar miRanda scores.
Moreover, we identified a subset of lung cancer patients with high miR-33a levels and low PD-1 expression that had a favorable outcome, suggesting a better prognostic value of miR-33a via PD-1 regulation. Taken together, our findings identify a novel mechanism of tumor immune evasion regulated by miR-33a via PD-1/PD-L1, with potential application in clinical practice.
Our data present novel insights into PD-1/PD-L1 signaling, first, regarding a novel reliable method for their evaluation, alternative to immunohistochemical testing, second, concerning miRNAs as modulators of immune checkpoint mechanisms, and finally regarding the potential value of miR-33a as a favorable prognostic marker via PD-1 regulation. MiRNA analysis combined with tumor expression of immune-biomarkers will improve our ability to select the best candidates to receive immune-based therapies, with important clinical benefits.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Iorio MV, Croce CM. MicroRNAs in cancer: small molecules with a huge impact. J Clin Oncol. 2009;27:5848–56.
Inui M, Martello G, Piccolo S. MicroRNA control of signal transduction. Nat Rev Mol Cell Biol. 2010;11:252–63.
Iorio MV, Croce CM. MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 2012;4:143–59.
Zhu C, Zhao Y, Zhang Z, Ni Y, Li X, Yong H. MicroRNA-33a inhibits lung cancer cell proliferation and invasion by regulating the expression of beta-catenin. Mol Med Rep. 2015;11:3647–51.
Brahmer J, Reckamp KL, Baas P, Crinò L, Eberhardt WE, Poddubskaya E, Antonia S, Pluzanski A, Vokes EE, Holgado E, Waterhouse D, Ready N, Gainor J, et al. Nivolumab versus docetaxel in advanced squamous–cell non–small–cell lung cancer. N Engl J Med. 2015;373:123–35.
Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, Chow LQ, Vokes EE, Felip E, Holgado E, Barlesi F, Kohlhäufl M, Arrieta O, et al. Nivolumab versus docetaxel in advanced nonsquamous non–small–cell lung cancer. N Engl J Med. 2015;373:1627–39.
Gettinger SN, Horn L, Gandhi L, Spigel DR, Antonia SJ, Rizvi NA, Powderly JD, Heist RS, Carvajal RD, Jackman DM, Sequist LV, Smith DC, Leming P, et al. Overall survival and long–term safety of nivolumab (anti–programmed death 1 antibody, BMS–936558, ONO–4538) in patients with previously treated advanced non–small–cell lung cancer. J Clin Oncol. 2015;33:2004–12.
Mu CY, Huang JA, Chen Y, Chen C, Zhang XG. High expression of PD-L1 in lung cancer may contribute to poor prognosis and tumor cells immune escape through suppressing tumor infiltrating dendritic cells maturation. Med Oncol. 2011;28:682–8.
Lin C, Chen X, Li M, Liu J, Qi X, Yang W, Zhang H, Cai Z, Dai Y, Ouyang X. Programmed death-ligand 1 expression predicts tyrosine kinase inhibitor response and better prognosis in a cohort of patients with epidermal growth factor receptor mutation-positive lung adenocarcinoma. Clin Lung Cancer. 2015;16:e25–35.
Konishi J, Yamazaki K, Azuma M, Kinoshita I, Dosaka-Akita H, Nishimura M. B7-H1 expression on non-small cell lung cancer cells and its relationship with tumor-infiltrating lymphocytes and their PD-1 expression. Clin Cancer Res. 2004;10:5094–100.
Smolle MA, Calin HN, Pichler M, Calin GA. Noncoding RNAs and immune checkpoints: clinical implications as cancer therapeutics. FEBS J. 2017;284:1952–66.
Travis WD, Brambilla E, Muller-Hemerlink HK, Harris CC. World Health Organization classification of tumours. Pathology and genetics of tumours of the lung, pleura, thymus and heart. Lyon: IARC Press; 2004.
Travis WD, Brambilla E, Noguchi M, Nicholson AG, Geisinger KR, Yatabe Y, Beer DG, Powell CA, Riely GJ, Van Schil PE, Garg K, Austin JH, Asamura H, Rusch VW, Hirsch FR, Scagliotti G, Mitsudomi T, Huber RM, Ishikawa Y, Jett J, Sanchez-Cespedes M, Sculier JP, Takahashi T, Tsuboi M, Vansteenkiste J, Wistuba I, Yang PC, Aberle D, Brambilla C, Flieder D, Franklin W, Gazdar A, Gould M, Hasleton P, Henderson D, Johnson B, Johnson D, Kerr K, Kuriyama K, Lee JS, Miller VA, Petersen I, Roggli V, Rosell R, Saijo N, Thunnissen E, Tsao M, Yankelewitz D. International association for the study of lung cancer/American thoracic society/European respiratory society international multidisciplinary classification of lung adenocarcinoma. J Thorac Oncol. 2011;6:244–85.
Travis WD, Brambilla E, Noguchi M, Nicholson A, Geisinger K, Yatabe Y, Ishikawa Y, Wistuba I, Flieder DB, Franklin W, Gazdar A, Hasleton PS, Henderson DW, Kerr KM, Petersen I, Roggli V, Thunnissen E, Tsao M. Diagnosis of lung cancer in small biopsies and cytology: implications of the 2011 international association for the study of lung cancer/American thoracic society/European respiratory society classification. Arch Pathol Lab Med. 2013;137:668–84.
Waggott D, Chu K, Yin S, Wouters BG, Liu FF, Boutros PC. NanoStringNorm: an extensible R package for the pre-processing of NanoString mRNA and miRNA data. Bioinformatics. 2012;28:1546–8.
Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, Thomson AM, Lim B, Rigoutsos I. A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006;126:1203–17.
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:189–98.
Boeri M, Pastorino U, Sozzi G. Role of microRNAs in lung cancer: microRNA signatures in cancer prognosis. Cancer J. 2012;18:268–74.
Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010;11:597–610.
Dong H, Strome SE, Salomao DR, Tamura H, Hirano F, Flies DB, Roche PC, Lu J, Zhu G, Tamada K, lennon VA, Celis E, Chen L. Tumor-associated B7-H1 promotes T-cell apoptosis: a potential mechanism of immune evasion. Nat Med. 2002;8:793–800.
Herbst RS, Soria JC, Kowanetz M, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. 2014;515:563–7.
Taube JM, Klein A, Brahmer JR, Xu H, Pan X, Kim JH, Chen L, Pardoll DM, Topalian SL, Anders RA. Association of PD-1, PD-1 ligands, and other features of the tumor immune microenvironment with response to anti-PD-1 therapy. Clin Cancer Res. 2014;20:5064–74.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12:252–64.
Postow MA, Callahan MK, Wolchok JD. Immune checkpoint blockade in cancer therapy. J Clin Oncol. 2015;33:1974–82.
Sznol M, Chen L. Antagonist antibodies to PD-1 and B7-H1 (PD-L1) in the treatment of advanced human cancer. Clin Cancer Res. 2013;19:1021–34.
Velcheti V, Schalper KA, Carvajal DE, Anagnostou VK, Syrigos KN, Sznol M, Herbst RS, Gettinger SN, Chen L, Rimm DL. Programmed death ligand-1 expression in non-small cell lung cancer. Lab Invest. 2014;9:107–16.
Chen YB, Mu CY, Huang JA. Clinical significance of programmed death-1 ligand-1 expression in patients with non-small cell lung cancer: a 5-year-follow-up study. Tumori. 2012;98:751–5.
Boland JM, Kwon ED, Harrington SM, et al. Tumor B7-H1 and B7-H3 expression in squamous cell carcinoma of the lung. Clin Lung Cancer. 2013;14:157–63.
Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA target. Cell. 2005;120:15–20.
Duursma AM, Kedde M, Schrier M, Le Sage C, Agami R. miR-148 targets human DNMT3b protein coding region. RNA. 2008;14:872–7.
Tay Y, Zhang J, Thomson AM, Lim B, Rigoutsos I. MicroRNAs to Nanog, Oct4 and Sox2 coding regions modulate embryonic stem cell differentiation. Nature. 2008;455:1124–8.
Rigoutsos I. New tricks for animal microRNAS: targeting of amino acid coding regions at conserved and nonconserved sites. Cancer Res. 2009;69:3245–8.
LB, MG and GF conceived and designed the experiments. MG performed the experiments. LB wrote the paper. GF diagnosed lung cancer. CN performed FFPE sections. FM, ML, and AM performed lung surgery and follow-up. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The data supporting the conclusions of this paper are included within the article.
Consent of publication
All authors are responsible for the submission of this article and accept the conditions of submission.
Ethics approval and consent to participate
Our study was conducted in accordance with the ethical standards of our institutional research committee and with the 1964 Helsinki declaration; all the patients gave their informed consent to the molecular analyses.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Boldrini, L., Giordano, M., Niccoli, C. et al. Role of microRNA-33a in regulating the expression of PD-1 in lung adenocarcinoma. Cancer Cell Int 17, 105 (2017). https://doi.org/10.1186/s12935-017-0474-y
- Lung adenocarcinoma