
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Oncology, Markers, Clinical Correlates |
Department of 1 Radiation Sciences, Oncology, and 2 Medical Biosciences, Pathology, University of Umeå, Umeå, Sweden; and 3 Center for Genomics and Bioinformatics, 4 Oncologic Center, Department of Surgical Sciences, and 5 Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Sweden
| ABSTRACT |
|---|
|
|
|---|
Experimental Design: Using 1624 prostate cancer cases and 801 unaffected controls, the truncating mutation E265X and five common sequence variants, including the two missense mutations R462Q and D541E, were evaluated for association between genotypes/haplotypes and prostate cancer risk.
Results: The prevalence of E265X carriers among unaffected controls and prostate cancer patients was almost identical (1.9 and 1.8% in controls and cases, respectively), and evidence for segregation of E265X with disease was not observed within any HPC family. Overall, the analyses of common sequence variants provided limited evidence for association with prostate cancer risk. We found a marginally significant inverse association between the missense mutation D541E and sporadic prostate cancer risk (odds ratio, 0.77; 95% confidence interval, 0.591.00) and reduced risk of prostate cancer in carriers of two different haplotypes being completely discordant.
Conclusions: Considering the high quality in genotyping and the size of this study, these results provide solid evidence against a major role of RNASEL in prostate cancer etiology in Sweden.
| INTRODUCTION |
|---|
|
|
|---|
RNASEL maps to HPC1 and is involved in the IFN-regulated 2'-5'linked oligoadenylates system that mediates cell proliferation and apoptosis (13 , 14) and has been suggested as a candidate tumor suppressor gene (15) . RNASEL was recently identified as a candidate gene for HPC1 through a positional cloning and candidate gene approach. A truncating mutation (E265X) and an initiation codon mutation (M1I) was reported to cosegregate within two hereditary prostate cancer (HPC; ref. 2 ) families linked to HPC1 (16) . Furthermore, loss of the wild-type allele and reduced RNASEL activity was observed in microdissected tumors carrying a germ-line mutation. In Finnish HPC families, RNASEL mutations did not explain disease segregation; however, the truncating mutation E265X was associated with prostate cancer risk (17) .
In addition to rare mutations, several polymorphisms in RNASEL have been reported. A common missense variant of RNASEL, R462Q, has been associated both with increased (17 , 18) and decreased (19 , 20) risk of HPC. Recently, Xiang et al. (21) reported that the R462Q variant reduces the ability of RNase L to cause apoptosis in response to activation by 2-5A. Hence, this polymorphism might have a functional role in prostate cancer development. Another missense variant, D541E, was associated with an increased risk for prostate cancer in a study from Japan (19) , although several other studies found no association (17 , 18 , 20) .
Additional evaluation of the impact that the RNASEL gene have on prostate cancer risk are required. Our aim was to assess the possible role of RNASEL sequence variants on prostate cancer risk in a Swedish population. We used a large population-based case-control study, composed of 1624 prostate cancer cases and 801 controls. The truncating mutation E265X and several common polymorphisms in RNASEL were evaluated, both for individual association and through haplotype analyses.
| MATERIALS AND METHODS |
|---|
|
|
|---|
79 years, diagnosed between July 1, 2001, and September 30, 2002, from the central and northern part of Sweden and all patients, ages
65 years, diagnosed in the same calendar period from the southeastern part of Sweden and the area of Stockholm. In total, 1961 prostate cancer patients were invited for participation and 1444 (74%) agreed to participate. All these cases donated a blood sample and answered a questionnaire, including questions regarding family history for prostate cancer. For all cases with at least one reported family member with prostate cancer, based on the initial questionnaire, a more detailed family history was obtained through a second questionnaire followed by a telephone interview.
Of the 1444 cases included in the study, 143 were classified as familial prostate cancer cases (two relatives with verified prostate cancer in a nuclear family) and 51 as HPC cases (three or more members with prostate cancer in a nuclear family). Detailed clinical data, including tumor-node-metastasis stage, Gleason score, prostate-specific antigen level at time of diagnosis, means of diagnosis and primary treatment were obtained through linkage to the National Prostate Cancer Registry. Cases were classified as either localized (tumor stage I-II and Gleason sum < 8/grade I-II and prostate-specific antigen < 100) or advanced (tumor stage III-IV or Gleason sum
8/grade III or prostate-specific antigen
100).
The continuously updated Swedish Population Registry encompasses the individually unique national registration number assigned to all residents. From this register, control subjects were randomly selected and frequency matched according to geographical origin (northern part of Sweden versus southeastern part of Sweden and the area of Stockholm) and to the expected age distribution of the cases (within 5-year age groups). The controls were recruited during the same calendar period as the cases. Of 1697 controls invited, 866 (51%) agreed to participate and to complete a blood draw and questionnaire (same as for cases). Subsequently, eight controls were excluded because linkage to the Swedish Cancer Registry revealed that they were diagnosed with a cancer of the prostate before inclusion. At the time of this study, DNA samples were available for 801 controls and 1427 cases (1247 sporadic prostate cancer cases, 131 familial prostate cancer cases, and 49 HPC cases).
Patients with Familial Prostate Cancer and HPC.
Recruitment of Swedish families with prostate cancer is an ongoing activity that started at the Department of Oncology, University of Umeå, in 1995. Identification of families has been based mainly on referrals by collaborating urologists and oncologists throughout Sweden but also on self-reported family history of prostate cancer from cases included into the Cancer Prostate in Sweden study, described above. Blood has been collected for as many family members as possible. For the present study, all affected family members with DNA information available (excluding cases from the Cancer Prostate in Sweden cohort already selected for this study) were included, resulting in 184 HPC cases and 25 familial prostate cancer cases. All diagnoses of prostate cancer in the families were confirmed both by reference to the Cancer Registry and by direct examination of medical records.
Total Study Population.
Our total study population is composed of 1636 prostate cancer patients and 801 unaffected controls. Of the 1636 prostate cancer patients, 156 (9.5%) and 233 (14.2%) were classified as familial prostate cancer and HPC cases, respectively. Pertinent characteristics of the study population are given in Table 1
. Written informed consent was obtained from each subject. The ethical committee at the Karolinska Institutet and University of Umeå approved the study.
|
13 kb. To comprehensively evaluate association between prostate cancer risk and common RNASEL sequence variants, we used the following approach. First, we defined the target region for selection of SNPs as 3 kb of the promoter, all exons, introns, and 3'-untranslated region. We then searched the public database National Center for Biotechnology Information (www.ncbi.nlm.nih.gov/SNP) within this target region, resulting in the identification of 34 SNPs. A subset of nine SNPs were then selected, with the criteria of a reported minor allele frequency of at least 5%, a density of one SNP per three kb, and as low as possible proportion of repetitive sequences. Preferentially validated SNPs were selected, and additional attention was paid to SNPs in the promoter region. Besides the two reported missense variants, R462Q and D541E, which were decided to be included in the analyses in advance, this subset consisted of four SNPs in the promoter region, one SNP in intron 4, one SNP in intron 5, and one SNP in 3'-untranslated region (Table 2)
|
Genotyping.
DNA samples were extracted from leukocytes with standard methods. Each DNA plate contained two Centre dEtude du Polymorphisme Humain controls, a water blank and blinded internal replicates. Two methods were used for genotyping. The stop mutation (E265X) and the two missense mutations (R462Q and D541E) were genotyped using a 5'-nuclease assay with TaqMan MGB probes. Primers and probes were designed with the Assay-by-Design service (Applied Biosystems, Inc., Foster City, CA). All reactions were performed in 10 µL of volume consisting of 10 ng of genomic DNA, 900 nmol/L of each primer, 200 nmol/L of each probe, and 5 µL of TaqMan universal master-mix. PCR cycling conditions were as follows: 50°C for 2 minutes, 95°C for 10 minutes, followed by 40 cycles of 92°C for 15 seconds and 60°C for 1 minute. The samples were analyzed on an ABI PRISM 7900 HT sequence detection system.
A dynamic allele-specific hybridization technique was adapted to genotype the remaining three SNPs. To thoroughly evaluate the genotype results, the missense variant R462Q was also genotyped with the dynamic allele-specific hybridization technique. Dynamic allele-specific hybridization was performed as described previously (24, 25, 26) . For this, two PCR primers and one dynamic allele-specific hybridization probe per target SNP were designed by custom DFold software (27) provided by DynaMetrix Ltd. (Hertfordshire, United Kingdom). These oligonucleotides were provided and HPLC purified, by Biomers GmbH (Ulm, Germany). The dynamic allele-specific hybridization PCRs entailed amplifying short genomic fragments spanning the variant of interest, with one of the primers carrying a 5'-biotin label. Amplifications were performed in 5 µL of volume, containing 1 to 2 ng of genomic DNA, 0.38 µmol/L biotinylated primer, 0.75 µmol/L nonbiotinylated primer, 0.03 units of AmpliTaq Gold (Applied Biosystems, Inc.), 10% dimethylsulphoxide, 1x AmpliTaq Gold Buffer, including 1.5 mmol/L MgCl2 (Applied Biosystems, Inc.), and 0.2 mmol/L each deoxynucleoside triphosphate. Thermal cycling was conducted on a MBS 384 device (Thermo-Hybaid, Thermo Electron Corporation, Milford, MA) as follows: 1x (10 minutes at 94°C) and 35x (15 seconds at 94°C, 30 seconds at annealing temperature). To verify successful amplification, 0.5 µL of several randomly chosen samples were examined on a 3.0% low-melt agarose gel.
Dynamic allele-specific hybridization analysis of the PCR product was conducted on membrane macroarrays, with the dynamic allele-specific hybridization-2 protocol (25) . Briefly, this entailed transferring samples to the membrane array by centrifugation (28) or robotic gridding. Resulting individual arrays with up to 9600 distinct features were rinsed in 0.1 mol/L NaOH to denature the PCR products. They were then exposed to 2 mL of HE buffer [0.1 mol/L HEPES and 10 mmol/L EDTA (pH 7.9)] containing 4 nmol of the appropriate ROX end-labeled probe. After heating to 85°C and cooling to room temperature, the membrane was briefly rinsed in HE buffer. The array was then soaked in 40 mL of HE buffer containing SYBR GreenI dye at 1:20,000 dilution for up to 3 hours. Using a dynamic allele-specific hybridization-2 device (DynaMetrix Ltd.), the membrane was taken through a dynamic allele-specific hybridization heating ramp (heating at 3°C/minute from room temperature to 85°C), whereas fluorescence from the ROX acceptor dye on the probe was monitored. Data were collected at intervals of 0.5°C. We used fluorescence changes with temperature (DNA melting profiles) to distinguish different alleles. This was done by the dynamic allele-specific hybridization-2 device software that uses negative derivatives of fluorescence against temperature to reveal peaks of denaturation rate (target-probe melting temperatures; Tm) and thereby automatically assign DNA samples into genotype groups.
Statistical Methods.
We tested Hardy-Weinberg equilibrium for each sequence variant and pairwise linkage disequilibrium between all sequence variants using a replication method (29)
. For each test, 10,000 permutations were performed, and the Fisher probability test statistic of each replicate was calculated from the new corresponding multilocus table. Empirical P values for each test were estimated as the proportion of replicate data sets found to be as probable as or less probable than the observed data, as implemented in the software package Genetic Data Analysis.
We used unconditional logistic regression to evaluate the association of each of the sequence variants with prostate cancer risk with the statistical software STATA. To adjust for the matching, we included indicator variables, representing each combination of age category (5-year age groups) and geographical region (northern part of Sweden versus southeastern part of Sweden and the area of Stockholm), in the regression model. A likelihood ratio test of a covariate equal to number of rare alleles (0, 1, or 2), corresponding to multiplicative allelic effects on the genotype relative risks, was used to test for association between genotypes and disease. To account for genotype correlations among individuals from the same family (in comparing familial prostate cancer or HPC cases with controls), a modified sandwich variance estimator in clustered logistic regression analysis was used (30) . Because the likelihood used for estimation is not a true likelihood, a Wald test was applied to test for association between genotypes and prostate cancer risk.
Association between haplotypes and prostate cancer risk was evaluated with a score test proposed by Schaid et al. (31) , as implemented in the software HAPLO.SCORE for the R programming language. On the basis of a generalized linear model, this method allows adjustment for possible confounding variables and provides both global and haplotype-specific tests. In these analyses, age and geographical region was adjusted for (as described above) and haplotypes with frequencies < 0.005 were pooled into a common group. Empirical P values, based on 10,000 simulations, were computed for the global score test and each of the haplotype-specific score tests.
No correction was made for multiple comparisons. All reported P values are two sided.
| RESULTS |
|---|
|
|
|---|
The truncating mutation, E265X, was observed in 15 (1.9%) controls, in 23 (1.9%) sporadic prostate cancer cases, and in five (1.4%) familial prostate cancer or HPC cases (Table 3)
. There was no significant difference in proportion of mutation carriers between controls and sporadic prostate cancer cases or between controls and familial prostate cancer or HPC cases. Among cases with HPC, four individuals (2.1%) were observed with the stop mutation. This proportion was not significantly different from the proportion observed among controls. The four HPC cases carrying the stop mutations are members of three different families. In one family, including five affected brothers with available DNA samples, two were mutation carriers, whereas the other three did not carry the mutation. In the other two HPC families, including one affected mutation carrier, no DNA samples were available from any other family member. The median age at sampling for the 15 control individuals carrying the mutation was 4 years younger than for the control individuals who were not carriers (P = 0.08, Mann-Whitney test). Stratified analyses based on age (<65 years or >65 years) and tumor aggressiveness (localized or locally advanced) did not reveal any significant differences in proportion of mutation carriers between controls and prostate cancer cases.
|
2 = 4.86, P = 0.03). However, between familial prostate cancer cases or HPC cases and controls, no significant differences were observed under a recessive model. In subgroup analyses based on age and tumor aggressiveness, no significant association between these variants and prostate cancer risk was observed.
Haplotype frequencies of the five common sequence variants were estimated in familial prostate cancer and HPC cases (selecting only the proband among related cases) and sporadic prostate cancer cases and controls (Table 4)
. A global score test provided no significant difference in haplotype frequency between sporadic prostate cancer cases and controls (P = 0.46), nor did a test of the maximum among all score-specific statistics (P = 0.10). This latter test should have greater statistical power compared with the global test if only a few haplotypes were associated with prostate cancer risk. Comparing familial prostate cancer and HPC cases with controls did not provide any support for association between the haplotypes and prostate cancer risk (global score test, P = 0.33; max-score test, P = 0.25). Individual haplotype analyses revealed two haplotypes significantly associated with prostate cancer risk. The frequency of the AAGAA haplotype of SNPs 2575 A>G, R462Q G>A, D541E G>T, +10721 G>A, and +13861 A>C was significantly decreased in sporadic prostate cancer cases compared with controls (2.9 versus 1.5%, P = 0.010). This haplotype also occurred at lower frequency among familial prostate cancer and HPC cases (1.6%) compared with controls, yet there was no significant difference (P = 0.181). The other haplotype significantly associated with prostate cancer risk was the GGTGC haplotype (containing the opposite allele of each SNP to the above associated haplotype) that also occurred at higher frequencies among controls (3.1%) compared with sporadic prostate cancer cases (2.5%, P = 0.415) and familial prostate cancer/HPC cases (0.8%, P = 0.020). The stop mutation E265X was in all cases observed on the most frequent haplotype AGTGC. Haplotype analyses, stratified by age and tumor aggressiveness, did not reveal any additional support for association between RNASEL sequence variants and prostate cancer risk (data not shown).
|
| DISCUSSION |
|---|
|
|
|---|
In the initial report of RNASEL as a candidate susceptibility gene for prostate cancer, an index case from each of 26 families at high risk for prostate cancer were screened for mutations in the RNASEL gene. Among these families, eight were linked to the HPC1 region with at least four affected individuals sharing an HPC1 haplotype. In this selected set of families, the stop mutation E265X, detected in one (3.8%) of the index cases, cosegregated with disease status within that family. Rökman et al. (17) reported the detection of E265X in five (4.3%) index patients from 116 Finnish families with HPC. In one of these five families, suggestive evidence of segregation between the mutation and prostate cancer was observed. In that study, they also compared the frequency of E265X mutations in unselected prostate cancer cases (n = 492) and control individuals (n = 566), without finding any significant association. In a study of 95 men with HPC, originating from 75 unrelated families, Chen et al. (32) found the E265X mutation in one family, with two of three affected brothers carrying the mutation. From these results, it is clear that only a small fraction of prostate cancer families can be explained by segregation of the E265X variant.
We observed virtually identical frequencies of E265X mutation carriers among sporadic prostate cancer cases and control individuals (1.9% in both groups). Among familial prostate cancer and HPC patients, expected to be enriched for genetic variants associated with prostate cancer risk, only 5 (1.4%) of 350 patients carried the E265X mutation; 4 of these had HPC and were members of three unrelated families. None of these families provided evidence for cosegregation between the mutation and disease status; in the only family with DNA samples available from other family members, two of five affected brothers were mutation carriers. Therefore, our results provide no support for the E265X mutation being a high-penetrance risk variant for prostate cancer. Furthermore, the complete lack of association between E265X and prostate cancer risk in our study provided no evidence for E265X being a modifier variant for prostate cancer risk.
The common missense mutation R462Q has been implicated in prostate cancer risk both from epidemiologic and functional studies. Using 423 prostate cancer cases and 454 sibling controls, Casey et al. (18) estimated that heterozygous carriers of this variant have 50% greater risk of prostate cancer than noncarriers, and homozygous carriers have more than double the risk. They also reported a reduction in enzymatic activity of the variant to one-third of normal activity. Additional functional evidence for this variants role in prostate cancer development comes from a recent study observing that the R462Q variant reduced the ability of RNase L to cause apoptosis in response to activation by 2-5A (21) . Rökman et al. (17) reported moderate evidence for an increased risk of prostate cancer for homozygous carriers of the R462Q variant in 66 index patients from HPC families. In addition, Wang et al. (20) reported a significant association between this variant and prostate cancer risk among 473 affected men from 181 families; however, they found an opposite trend with odds ratios of 0.8 for heterozygotes and 0.5 for homozygotes.
Despite the functional and epidemiologic findings supporting a role for R462Q in prostate cancer susceptibility, no significant associations between this variant and prostate cancer risk was observed in our study. Nor did we find any significant association with prostate cancer risk for any of the other four common sequence variants studied, except for the missense mutation D541E, which under a recessive genetic model, provided marginally significant evidence for association with sporadic prostate cancer (P = 0.03), with an estimated odds ratio of 0.77 (95% confidence interval, 0.591.00) for homozygous carriers compared with homozygous noncarriers. However, because of the large number of tests conducted in our study, this could represent a false positive finding. Moreover, no differences were observed between patients with familial prostate cancer or HPC and controls for this variant, irrespective of choice of genetic model (recessive, multiplicative, or dominant). This missense variant was significantly associated with an increased prostate cancer risk in a study from Japan based on 101 familial prostate cancer cases and 105 controls (19) , yet other studies have failed to find an association between this variant and prostate cancer risk (17 , 18 , 20) .
Haplotype analysis of the five SNPs provided weak support for RNASEL as a risk factor for prostate cancer. Global tests for association between haplotypes and prostate cancer risk were all nonsignificant; however, two specific haplotypes were significantly associated with prostate cancer risk. For both haplotypes, higher frequency was observed among controls compared with prostate cancer cases, indicating protective effects. It is therefore puzzling that these haplotypes were completely discordant for each of the five SNPs, and, as discussed above, due to the many statistical tests performed, the chance of false positive findings has to be considered.
Our study design has a certain limitation. There is the possibility of the existence of other common variants in RNASEL that are not well represented by the five SNPs genotyped in the present study. We used the public database National Center for Biotechnology Information for selection of SNPs, which is suboptimal compared with sequencing an adequate number of individuals from the study population to uncover all existing variants. However, considering the strong linkage disequilibrium between the studied SNPs, we believe that it is most likely that we were able to captured the major part of the genetic variation in the RNASEL gene in our study. Furthermore, because of the large number of individuals in this study, the possibility that our results represent false negative findings is of minimal concern. For example, comparing sporadic prostate cancer patients with controls, using the genotype frequencies observed for the R462Q variant among the controls and assuming a multiplicative genetic model, we have >87% power to detect an odds ratio of 1.5 for homozygotes at a significance level of 5% (two-sided test).
In summary, we failed to support a role of the rare truncating mutation E265X in Swedish prostate cancer. In addition, analyses of common sequence variants, both individually and through haplotypes, provided limited evidence for association with prostate cancer risk. Considering the large and well-characterized study population and the high quality in genotyping, these results provide strong evidence against a role of RNASEL in prostate cancer etiology in Sweden.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
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.
Requests for reprints: Fredrik Wiklund, Department of Radiation Sciences, Oncology, University of Umeå, 901 87 Umeå, Sweden. Phone: 46-90-7852270; Fax: 46-90-127464; E-mail: fredrik.wiklund{at}oc.umu.se
Received 5/19/04; accepted 7/26/04.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. J. Shook, J. Beuten, K. C. Torkko, T. L. Johnson-Pais, D. A. Troyer, I. M. Thompson, and R. J. Leach Association of RNASEL Variants with Prostate Cancer Risk in Hispanic Caucasians and African Americans Clin. Cancer Res., October 1, 2007; 13(19): 5959 - 5964. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Liu, F. R. Schumacher, S. J. Plummer, E. Jorgenson, G. Casey, and J. S. Witte trans-Fatty acid intake and increased risk of advanced prostate cancer: modification by RNASEL R462Q variant Carcinogenesis, June 1, 2007; 28(6): 1232 - 1236. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Fan A new human retrovirus associated with prostate cancer PNAS, January 30, 2007; 104(5): 1449 - 1450. [Full Text] [PDF] |
||||
![]() |
R. J. Molinaro, B. K. Jha, K. Malathi, S. Varambally, A. M. Chinnaiyan, and R. H. Silverman Selection and cloning of poly(rC)-binding protein 2 and Raf kinase inhibitor protein RNA activators of 2',5'-oligoadenylate synthetase from prostate cancer cells Nucleic Acids Res., December 2, 2006; 34(22): 6684 - 6695. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Li and B. C. Tai RNASEL Gene Polymorphisms and the Risk of Prostate Cancer: a Meta-analysis. Clin. Cancer Res., October 1, 2006; 12(19): 5713 - 5719. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Orr-Urtreger, A. Bar-Shira, D. Bercovich, N. Matarasso, U. Rozovsky, S. Rosner, S. Soloviov, G. Rennert, L. Kadouri, A. Hubert, et al. RNASEL Mutation Screening and Association Study in Ashkenazi and Non-Ashkenazi Prostate Cancer Patients. Cancer Epidemiol. Biomarkers Prev., March 1, 2006; 15(3): 474 - 479. [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 | Cell Growth & Differentiation |