| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology, and Genetics |
Departments of 1 Molecular, Cellular and Developmental Biology, 2 Therapeutic Radiology, 3 Genetics, 4 Pathology, and 5 Epidemiology, Yale University, New Haven, Connecticut; and 6 Department of Environmental Health, Harvard School of Public Health and 7 Pulmonary and Critical Care Unit, Department of Medicine, Massachusetts General Hospital/Harvard Medical School, Boston, Massachusetts; and 8 Department of Internal Medicine, University of New Mexico and 9 Lung Cancer Program, Lovelace Respiratory Research Institute, Albuquerque, New Mexico
Requests for reprints: Joanne B. Weidhaas and Frank J. Slack, Yale University, 333 Cedar Street, New Haven, CT 06520. Phone: 203-737-4267; Fax: 203-785-6309; E-mail: joanne.weidhaas{at}yale.edu and frank.slack{at}yale.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The let-7 family of miRNAs seems to play a key role in lung cancer: they are at low levels in non–small cell lung cancer (NSCLC; refs. 12, 13); their lower levels are biomarkers of a poor outcome (14, 15); they regulate multiple important lung cancer oncogenes, including RAS (12, 16); and they inhibit growth of lung cancer cell lines in vitro (12, 14) and in vivo (17, 18). The purpose of this study was to identify SNPs that could potentially modify let-7 binding and to assess the effect of these SNPs on target gene regulation and the risk for NSCLC. Here, we show that a variant allele at a SNP in a LCS in the KRAS 3'UTR is associated with increased risk for NSCLC in moderate smokers. Furthermore, this variant allele leads to altered KRAS regulation in vitro, with higher KRAS expression in the presence of the variant allele. In addition, in tumors from patients harboring the variant allele, let-7s were lower than in nonvariant allele tumors, suggesting the SNP might be associated with NSCLC with a poor prognosis. The finding that a SNP disrupting let-7 miRNA regulation of a known oncogene can affect cancer predisposition to NSCLC is a new paradigm and supports further studies to identify similar SNPs in all cancer types.
| Materials and Methods |
|---|
|
|
|---|
To determine the frequency of the SNP alleles, 2,433 individuals were genotyped from a global sample of 46 populations. According to population ancestry and geographic locations, these 46 populations are categorized into 4 groups: European (including West Asia), African, Asian (including the Pacific), and Native American. Sample descriptions and samples sizes can be found in the Allele Frequency Database (ALFRED; 19) by searching for the population names.10 DNA samples were extracted from lymphoblastoid cell lines established and/or grown in the Yale University laboratory of K.K.K. The methods of transformation, cell culture, and DNA purification have been described. (20) All volunteers were apparently normal and otherwise healthy adult males or females and samples were collected after receipt of appropriate informed consent.
Lung cancer cases (n = 325) for the New Mexico case-control study were recruited beginning in 2004 from Albuquerque through two local hospitals, the Veterans hospital and the University of New Mexico hospital. All stages and histologic types of lung cancer were included. Controls (n = 325) with no history of any prior cancer were recruited from two ongoing local smoker cohorts, the Veterans Smokers Cohort (mainly veterans from Albuquerque), and the Lovelace Smokers Cohort (general residents in Albuquerque). Those two cohorts started to recruit participants in 2001 to conduct longitudinal studies on molecular markers of respiratory carcinogenesis in biological fluids such as sputum from people at risk for lung cancer. A standardized questionnaire was used to collect information on medical, family, and smoking exposure history, and quality of life for both lung cancer cases and control cohort members. Controls were randomly matched to lung cancer cases after categorization into different age groups (5-y differences) by sex and cohort (see Supplementary Table S2A). Cases with small cell lung cancer were excluded to more precisely assess the effect of the let-7 complementary site (LCS)6 SNP on risk for NSCLC. Cases over ages 82 y (the maximum age in the control group), cases with any prior cancer history, never smokers, or cases with missing data on smoking-related covariates were also excluded in the data analysis, resulting in 218 cases in the analysis.
A second lung cancer case-control study was conducted to validate findings from the New Mexico study. The study population was derived from a large ongoing molecular epidemiologic study in Boston, MA, that began in 1992 and now has >2,205 NSCLC patients. Details of this case-control population have been described previously (21–23). For this study >3,700 samples were analyzed, which included smokers and nonsmokers. This study was approved by the Human Subjects Committees of Massachusettts General Hospital and Harvard School of Public Health, Boston, MA. Briefly, all histologically confirmed, newly diagnosed patients with NSCLC at Massachusettts General Hospital were recruited between December 1992 and February 2006. Before 1997, only early stage (stage I and II) patients were recruited. After 1997, all stages of NSCLC cases were recruited in this study. Controls were recruited at Massachusettts General Hospital from healthy friends and nonblood-related family members (usually spouses) of several groups of hospital patients: (a) patients with cancer, whether related or not related to a case; or (b) patients with a cardiothoracic condition undergoing surgery. No matching was performed. Importantly, none of the controls were patients. Potential controls with a previous diagnosis of any cancer (other than nonmelanoma skin cancer) were excluded from participation. Over 85% of eligible cases and over 90% of controls participated in this study and provided blood samples. A research nurse administered questionnaires on demographic information and detailed smoking history (see Supplementary Table S2B). To reduce potential variation in allele frequency by ethnicity, only Caucasians were considered in the analysis.
Evaluation of 3'UTR sequences and the LCS6 SNP. DNA was isolated from frozen and FFPE lung tissue using the DNeasy Blood and Tissue kit (Qiagen). Segments of the KRAS 3'UTR were amplified using PfuTurbo DNA polymerase (Stratagene) and DNA primers specific to this sequence (see Supplementary Table S3). PCR products were purified using the QIAquick PCR Purification kit or 96 PCR Purification kit (Qiagen) and sequenced using the same primers. The NRAS 3'UTR was sequenced in the same manner.
For high-throughput genotyping, the DNA isolated from lymphocytes, blood, or tumor samples was amplified using TaqMan PCR assays designed specifically to identify the T or G allele of the LSC6 SNP (Applied Biosciences).
Determining the effect of the LCS6 variant allele on KRAS expression. We generated a pGL3 derivative containing almost the entire KRAS 3'UTR (KRAS wild-type) as follows. KRAS wild-type includes 3,910 bp of the KRAS 3'UTR, which was amplified from human genomic DNA using the forward primer SMJ104 and reverse primer LCJ5 (see Supplementary Table S3). NheI restriction sites were included on the 5' ends of the primers for convenient cloning. The product was first cloned into the TOPO cloning vector (Invitrogen) and then subcloned into pGL3 (Ambion) for use in subsequent luciferase assays. The luciferase reporter with the variant LCS6 KRAS 3'UTR (KRAS mLCS6) was constructed through site-directed mutagenesis of KRAS wild-type using GeneTailor (Invitrogen). A549 cells were cultured in DMEM with 10% fetal bovine serum and penicillin/streptomycin (Invitrogen). A549 cells were transfected with 500 ng KRAS wild-type or KRAS mLCS6 and 50 ng pRL-TK (Promega) using Lipofectamine 2000 (Invitrogen) for 24 h. Reporter expression was analyzed using the Dual-Luciferase Reporter Assay (Promega) and Wallac Victor2 1420 (PerkinElmer). Two-tailed t tests were done to verify statistical significance of differences in luciferase expression using GraphPad Prism.
Measuring let-7 in vivo. RNA was isolated from normal or tumor tissue specimens using Ambion isolation kits. RTPCR was done for let-7a, b, d, and g using the Ambion RTPCR specific primers for eight tumor samples with the variant allele and eight samples without the variant allele, using the ABI 7900.
Statistical analysis. To calculate significance, a
2 test was used for categorical variables, a t test was used for continuous variables, and in some cases, a two-sided Fisher's exact test was used. Two-sided two-sample t tests,
2 analyses, and two-sided Wilcoxon rank-sum tests were performed, as appropriate to compare the demographic variables between cases and controls. An unconditional logistic regression model was used to calculate odds ratio (OR) and confidence intervals (CI) for the KRAS LCS6 SNP in moderate and heavy smokers (defined by below and above the median pack-years for the study population, respectively) with adjustment for selected covariates. A likelihood ratio test was used to assess the association with the allele and the pack-year interaction for NSCLC. The dominant model was used for all genetic association analysis due to the low frequency of the rare allele.
| Results |
|---|
|
|
|---|
|
Frequency of the variant allele across world populations. To determine if the prevalence of the variant allele in our NSCLC population was higher than expected, we determined the allele frequencies of the LCS6 SNP in the general population with a unique human genetic resource found at Yale University, a collection of genomic DNA from 2,433 healthy individuals from a global set of 46 populations (26). An extensive data base of genetic variations in these samples exists and can be found, along with the population descriptions in ALFRED (19). Using a TaqMan assay, we found that <3% of the 4,866 chromosomes, or 5.8% of the people tested, had the G allele (variant) at the LCS6 SNP site (Fig. 2 ). The frequency of this allele varied across geographic populations, with "European" populations exhibiting the variant allele most frequently (7.6% of the chromosomes tested); African populations less frequently (<2.0% of chromosomes tested); and "Asian" and Native American populations infrequently (<0.4% of chromosomes tested). Of note, over 85% of the patients in our retrospective patient cohort were of European descent. Importantly, these findings indicate that the prevalence of the variant allele (20.3%) in our NSCLC patient cohort is significantly higher than expected in any existing geographic population, further supporting the hypothesis that that this variant allele is a marker of an increased risk to NSCLC.
|
|
40 pack-years. Our finding that the LCS6 SNP only affects cancer risk for those with less cigarette exposure agrees with other studies showing a dose-dependent gene-environment interaction for smoking-induced lung cancer risk (23, 27, 28); with higher smoking exposure, any genetic predisposition is hypothesized to be overwhelmed by the extent of smoking-related damage. The LCS6 variant allele affects KRAS regulation and is associated with lower let-7. To determine if LCS6 affects KRAS regulation, we used a previously described luciferase reporter for KRAS expression to determine the effect of the LCS6 variant on expression (12). We transfected a luciferase reporter with a full-length KRAS 3'UTR containing the variant allele at LCS6 (KRAS mLCS6) into A549 cells, a lung cancer cell line, with known low let-7 levels (14), and compared luciferase expression to A549 cells transfected with the luciferase reporter with a wild-type KRAS 3'UTR (KRAS wild-type). We found that there was a significant increase in luciferase activity in cells transfected with KRAS mLCS6 versus the KRAS wild-type (Fig. 3A and B ). These findings support the hypothesis that the variant allele has altered the dependence of the KRAS 3'UTR on LCS6, and allows increased KRAS expression in the presence of this variant allele.
|
We next measured the levels of let-7 in available patient tumors with or without the variant allele, using reverse transcription-PCR. On average, we found that the levels of let-7a, b, d, and g were lower in patients with the variant allele compared with patients without the variant allele (Fig. 3C). These findings suggest that the variant allele is associated with lower let-7, at least in the NSCLC tumors in this study.
| Discussion |
|---|
|
|
|---|
We find that this variant allele at the LCS6 site indeed affects regulation of KRAS expression in vitro, allowing increased levels of KRAS. Interestingly, tumors containing the variant allele had lower let-7 levels than tumors without the variant allele, and as low let-7 has been associated with a poor prognosis in NSCLC (14, 15, 29), these findings raise the possibility that the subset of NSCLC patients harboring the variant allele may be those with an especially poor prognosis, a hypothesis requiring future studies to validate.
Although it is not possible at this time to further define the mechanism of cancer predisposition for patients who harbor this variant allele, the finding that the variant allele seems to allow increased KRAS in vitro suggests that KRAS overexpression might be one plausible mechanism. An additional hypothesis could be that altered let-7 binding in the KRAS 3'UTR could somehow lower cellular levels of let-7, perhaps through a feedback loop, and lower let-7 levels could further result in increased cell growth, as let-7 is known to repress cell growth pathways (16). For example, recent evidence indicates that a negative feedback loop involving lin-28 specifically regulates cellular let-7 levels (30). In the background of the variant allele, a cycle creating lowered let-7 levels and increased KRAS expression could in concert act as the first steps in oncogenesis. These hypotheses need testing in future studies.
Because lung cancer is so deadly when caught at later stages, screening programs have been initiated in current and exsmokers: The Early Lung Cancer Action Project found that a chest computed tomography (CT) scan is thrice more sensitive in detecting early-stage lung cancer than a chest X-ray in "high-risk" populations (2.4% versus 0.7%; refs. 31, 32). Yet there remains considerable controversy over the use of lung CT scans as a global screening approach for lung cancer because of the expense (estimated cost 2 billion dollars yearly in the United States alone) and the very low yield of cancers detected yearly (1.2%; ref. 33). One of the primary problems is that although smoking is the number one risk factor for developing lung cancer, only 10% of smokers ever develop lung cancer. With 44.5 million current smokers in the United States (20.9% of the population) and over 1.3 billion smokers worldwide, there is a clear need to identify markers whose genetic variation is associated with lung cancer risk that would help better prioritize who should be screened clinically and be offered chemopreventive agents.
Although numerous studies of carcinogen metabolizing and detoxifying genes and DNA repair genes have identified sequence variations associated with risk for lung cancer, a meta-analysis of such polymorphisms in DNA repair pathways concluded that for any increased risk association between a single SNP and lung cancer, the risk fluctuation would likely be minimal, and only panels evaluating a collection of SNPs would ever successfully predict lung cancer risk (34). Supporting this hypothesis, 2 SNPs associated with lung cancer risk (OR, 1.19–1.8) were recently identified from three large studies examining over 300,000 SNPs. The mechanism of lung cancer predisposition caused by these SNPs is controversial but is hypothesized to be through affects on a nicotine receptor (35, 36).
In contrast, in this study, we have identified the first miRNA binding-site SNP that alone can predict a significant increase in NSCLC risk in people with a moderate smoking history. The mechanism may be due to altered regulation of the KRAS oncogene and possibly cellular let-7 levels as well. These findings give insight into a new paradigm and support the pursuit of 3'UTR sequencing for similar SNPs in all tumor-related genes to better understand their role in genetic cancer risk. This strategy can be a complementary, and likely productive, approach to enhance current efforts to define genetic cancer risk.
| Disclosure of Potential Conflicts of Interest |
|---|
|
|
|---|
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank William C. Speed and Dr. Judith R. Kidd for their assistance in setting up and performing the population genetic screen, Welela Tereffe and Lynn Wilson for their critical reading of the manuscript, and Drs. Kofi Asomaning and Yen-tsung Huang for the assistance with the statistical analysis for the Boston study.
| Footnotes |
|---|
F.J. Slack and J.B. Weidhaas contributed equally.
10 http://alfred.med.yale.edu/ ![]()
Received 6/ 6/08. Revised 7/28/08. Accepted 7/29/08.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. A. Edson, A. K. Nagaraja, and M. M. Matzuk The Mammalian Ovary from Genesis to Revelation Endocr. Rev., October 1, 2009; 30(6): 624 - 712. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Song, H. Zheng, B. Liu, S. Wei, H. Dai, L. Zhang, G. A. Calin, X. Hao, Q. Wei, W. Zhang, et al. An miR-502-Binding Site Single-Nucleotide Polymorphism in the 3'-Untranslated Region of the SET8 Gene Is Associated with Early Age of Breast Cancer Onset Clin. Cancer Res., October 1, 2009; 15(19): 6292 - 6300. [Abstract] [Full Text] [PDF] |
||||
![]() |
S M Metias, E Lianidou, and G M Yousef MicroRNAs in clinical oncology: at the crossroads between promises and problems J. Clin. Pathol., September 1, 2009; 62(9): 771 - 776. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Sotiropoulou, G. Pampalakis, E. Lianidou, and Z. Mourelatos Emerging roles of microRNAs as molecular switches in the integrated circuit of the cancer cell RNA, August 1, 2009; 15(8): 1443 - 1461. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Fabbri, N. Valeri, and G. A. Calin MicroRNAs and genomic variations: from Proteus tricks to Prometheus gift Carcinogenesis, June 1, 2009; 30(6): 912 - 917. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. C. Christensen, B. J. Moyer, M. Avissar, L. G. Ouellet, S. L. Plaza, M. D. McClean, C. J. Marsit, and K. T. Kelsey A let-7 microRNA-binding site polymorphism in the KRAS 3' UTR is associated with reduced survival in oral cancers Carcinogenesis, June 1, 2009; 30(6): 1003 - 1007. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |