Role of microRNA-33a in regulating the expression of PD-1 in lung adenocarcinoma
- Laura Boldrini†1Email author,
- Mirella Giordano†1,
- Cristina Niccoli1,
- Franca Melfi1,
- Marco Lucchi1,
- Alfredo Mussi1 and
- Gabriella Fontanini1
© The Author(s) 2017
Received: 19 June 2017
Accepted: 9 November 2017
Published: 17 November 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–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–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–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 in 88 lung adenocarcinoma patients
miR33a expression (mean ± SD)
276.4 ± 24.9
175.2 ± 24.8
197.8 ± 22.6
274.8 ± 29.9
Adenocarcinoma prevalent pattern
242.7 ± 31.8
185.3 ± 33.6
272.9 ± 36.5
182.6 ± 51.6
634 ± 88.5
231.3 ± 20.1
168.6 ± 29.5
267.1 ± 26.7
191.3 ± 24.4
miR-33a and survival analysis
miR-33a target prediction
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).
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–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–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.
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.
Open AccessThis 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.
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Iorio MV, Croce CM. MicroRNAs in cancer: small molecules with a huge impact. J Clin Oncol. 2009;27:5848–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Inui M, Martello G, Piccolo S. MicroRNA control of signal transduction. Nat Rev Mol Cell Biol. 2010;11:252–63.View ArticlePubMedGoogle Scholar
- Iorio MV, Croce CM. MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 2012;4:143–59.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Smolle MA, Calin HN, Pichler M, Calin GA. Noncoding RNAs and immune checkpoints: clinical implications as cancer therapeutics. FEBS J. 2017;284:1952–66.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Boeri M, Pastorino U, Sozzi G. Role of microRNAs in lung cancer: microRNA signatures in cancer prognosis. Cancer J. 2012;18:268–74.View ArticlePubMedGoogle Scholar
- Krol J, Loedige I, Filipowicz W. The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010;11:597–610.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12:252–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Postow MA, Callahan MK, Wolchok JD. Immune checkpoint blockade in cancer therapy. J Clin Oncol. 2015;33:1974–82.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Duursma AM, Kedde M, Schrier M, Le Sage C, Agami R. miR-148 targets human DNMT3b protein coding region. RNA. 2008;14:872–7.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Rigoutsos I. New tricks for animal microRNAS: targeting of amino acid coding regions at conserved and nonconserved sites. Cancer Res. 2009;69:3245–8.View ArticlePubMedGoogle Scholar