| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Epidemiology and Prevention |
1 Cancer Research UK Department of Oncology and 2 Cancer Research UK Genetic Epidemiology Unit, Strangeways Research Laboratory, University of Cambridge, Cambridge, United Kingdom; 3 Translational Research Laboratory, Department of Gynaecological Oncology, Institute for Women's Health, University College London, London, United Kingdom; 4 Department of Viruses, Hormones, and Cancer, Institute of Cancer Epidemiology, Danish Cancer Society; 5 Gynaecologic Clinic, Juliane Marie Centre, Rigshospitalet, University of Copenhagen, Copenhagen, Denmark; 6 Department of Cancer Genetics, Roswell Park Cancer Institute, Buffalo, New York; and 7 Department of Health Research and Policy, Stanford University School of Medicine, Stanford, California.
Requests for reprints: Honglin Song, Cancer Research UK Department of Oncology, Strangeways Research Laboratory, University of Cambridge, Worts Causeway, Cambridge, United Kingdom CB1 8RN. Phone: 44-1223-740161; Fax: 44-1223-740147; E-mail: honglin{at}srl.cam.ac.uk.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The task of identifying the genes and genetic variants that might confer susceptibility to ovarian cancer is currently restricted, for the most part, to selecting functionally plausible candidates. For human cancers, pathways involved in the control of replication and cellular proliferation provide many possible candidate genes (3). During mitosis, faithful replication of DNA during S phase and equal distribution of the identical chromosomal copies to the daughter cells during M phase is essential (4). The G1-S phase of the cell cycle is crucial for the decision of cell whether to commit to growth arrest or proliferation. The function of the retinoblastoma protein, pRb (RB1 gene product), is to connect the cell cycle control with the transcriptional machinery of cell. pRb is a negative regulator of cellular proliferation, which is achieved by sequestering a variety of nuclear proteins involved in cellular growth. When in a hypophosphorylated state, pRb exerts its antiproliferative function by sequestering and altering the function of the E2F family of transcription factors that control the expression of a bank of genes essential for cells to progress from G1 into the S phase (5). Conditions that favor phosphorylation of pRb favor cell proliferation (6).
The RB1 gene spans
180 kb of genomic DNA on chromosome 13q14 and consists of 27 exons, encoding 928 amino acids. The largest intron (IVS 17) contains an open reading frame encoding the G proteincoupled receptor P2RY5 (purinergic receptor P2Y, G protein coupled, 5 NM_005767) in the reverse orientation relative to the transcription of RB1 (7). The importance of RB1 in human cancers and the role of the pRb pathway in cell proliferation and cancer development is well established. Most human cancers show somatic alteration, either in pRb or in one or other components of the pathway (8). Somatic inactivation of RB1 has been reported in retinoblastomas (9), breast cancer (10), small cell lung cancers (11), and many sarcomas and bladder carcinomas (12). Germ-line mutations in RB1 predispose to retinoblastoma (13). Very recently, polymorphisms of RB1, rs2854344 and rs198580, have been reported to be associated with altered breast cancer risk in the British population; carrying the minor allele of these single nucleotide polymorphisms (SNP) seems to reduce breast cancer risk (14).
Alterations in the pRb pathway are observed frequently in epithelial ovarian cancer (EOC);
30% of ovarian cancers exhibit loss of heterozygosity at the RB1 gene locus, but no mutations have been detected in the remaining RB1 allele (15). Systematic analysis of the genes involved in the entire pRb pathway (p16-cyclin-dependent kinase 4/cyclin D1-pRb) revealed that 80% of EOCs have abnormalities in this pathway (16). Recent research shows that the overall survival rate for EOC patients with a normal pRb pathway is significantly higher than for patients with an altered pRb pathway (17). Because RB1 is key in a regulatory pathway that suffers disruption during the pathogenesis of various human tumors and there is some evidence of its association with breast cancer, it is reasonable to consider that common variants in the gene might explain some of the interindividual variability in the risk of EOC.
Population-based genetic association studies using large sets of cases and controls are a powerful method for identifying common low/moderate penetrance disease susceptibility alleles. The aim of this study was to investigate whether common variants in RB1 gene were associated with ovarian cancer risk using a SNP tagging approach in a large, multicenter case-control study.
| Materials and Methods |
|---|
|
|
|---|
Tag SNP selection. Resequencing data from the National Institute of Environmental Health Sciences Environmental Genome Project (EGP) were available for RB1 to enable a systematic selection of tagging SNPs (stSNP).8 The EGP has been resequencing candidate genes for cancer across panels of individuals representative of U.S. ethnicities. The original panel (P1-PDR90) of 90 individuals consists of 24 European Americans, 24 African Americans, 12 Mexican Americans, 6 Native Americans, and 24 Asian Americans, but the ethnic group identifiers are not available. It is known that there is greater genetic and haplotype diversity in individuals of African origin and so we have identified and excluded 28 of the samples with the greatest African ancestry in this population by comparing the genotypes of the PDR90 sample with genotypes for the same SNPs from the National Heart, Lung, and Blood Institute Variation Discovery Resource Project African American panel.9 Data from the remaining 62 individuals were used to identify stSNPs. Exclusion of the samples from Native American, Hispanic American, and Asian American individuals would be ideal, but because there is less genetic diversity between these groups, this cannot be done with any certainty.
We used the program Tagger to select stSNPs. Tagger uses a strategy that combines the simplicity of pairwise methods with the potential efficiency of multimarker approaches (19).10 The aim of the SNP tagging was to identify a set of stSNPs that efficiently tags all the known common variants [minor allele frequency (MAF) >0.05] and is likely to tag most of the unknown common variants in the gene. The best measure of the extent to which one SNP tags another SNP is the pairwise correlation coefficient (rp2) because the loss in power incurred by using a marker SNP in place of a true causal SNP is directly related to this measure. We aimed to define a set of tagging SNPs such that all known common SNPs had an estimated rp2 of >0.8, with at least one tagging SNP. However, some SNPs are poorly correlated with other single SNPs but may be efficiently tagged by a haplotype defined by multiple SNPs, so called "aggressive tagging," thus reducing the number of tagging SNPs needed. As an alternative, therefore, we aimed for the correlation between each SNP and a haplotype of tagging SNPs (rS2) to be at >0.8.
Genotyping. All samples were genotyped using the Taqman 7900HT Sequence Detection System according to the manufacturer's instructions. Each assay was carried out using 10 ng DNA in a 5 µL reaction using Taqman Universal PCR Master Mix (Applied Biosystems, Warrington United Kingdom), forward and reverse primers, and FAM- and VIC-labeled probes designed by Applied Biosystems (ABI Assay-by-Design). Primer and probe sequences and assay conditions used for each polymorphism analyzed are available from the corresponding author on request. All assays were carried out in 384-well arrays with 12 duplicate samples in each plate for quality control. Where discordant genotypes were observed in duplicates, the genotyping was repeated. Genotypes were determined using Allelic Discrimination Sequence Detection software (Applied Biosystems). DNA samples that did not give a clear genotype result at the first attempt were not repeated. Hence, there are variations in the number of samples successfully genotyped for each polymorphism. Call rates ranged from 94.5% to 99.5% for all the studies and SNPs and overall concordance between duplicate samples was >99%.
Statistics. The three studies were treated as separate strata in the analyses. Deviations of the genotype frequencies in the controls from those expected under Hardy-Weinberg equilibrium (HWE) were assessed by
2 tests [1 degree of freedom (df)]. The primary tests of association were the univariate analyses between each tagging SNP and ovarian cancer. Genotype frequencies in cases and controls were compared for each study separately using
2 tests for heterogeneity (2 df). The data from the three studies were then pooled and genotype frequencies were compared in cases and controls using unconditional logistic regression with terms for genotype and study and an appropriate likelihood-ratio test. We tested for heterogeneity between study strata by comparing logistic regression models with and without a genotype-stratum interaction term using likelihood ratio tests. Genotypic specific risks with the common homozygote as the baseline comparator were estimated as odds ratios (OR) with associated 95% confidence intervals (95% CI) by unconditional logistic regression. We also tested for rare allele dose effect (assuming a multiplicative model) using
2 tests (1 df) for each study separately and unconditional logistic regression for the pooled data.
In addition to the univariate analyses, we carried out a global haplotype test for association of the common haplotypes with EOC. Haplotype frequencies and subject-specific expected haplotype indicators were calculated separately for each study using the program TagSNPs, which implements an expectation-substitution approach to account for haplotype uncertainty given unphased genotype data (20). Subjects missing >50% genotype data were excluded. We considered haplotypes with >2% frequency in at least one study to be "common." Rare haplotypes were pooled. We used unconditional logistic regression to test the global null hypothesis of no association between haplotype frequency and ovarian cancer stratified by study, by comparing a model with multiplicative effects for each common haplotype (treating the most common haplotype as the referent) to the intercept-only model. Haplotype-specific ORs were also estimated with their associated confidence intervals.
| Results |
|---|
|
|
|---|
Association analysis of SNPs with ovarian cancer risk. We have genotyped these 11 SNPs in three ovarian cancer case-control studies. Genotype distributions in controls were consistent with HWE, except for rs425834, in the MALOVA study (P = 0.02). This is likely to be a chance observation because the discrimination of genotypes for this assay was good and deviation from HWE was not seen in cases.
The genotype frequencies in controls were similar for all three study populations with minor differences for rs2854345, rs425834, and rs4151540 (Table 1
). The observed genotype frequencies by study and for the combined data set are presented in Table 1, which also shows the result of the test for the comparison of genotype frequencies (Pheterogeneity) between cases and controls. The genotypic specific risks for the individual studies and combined data and the result of the trend test for association are presented in Table 2
. There was no association in controls between age and genotype frequency for any of the SNPs, and age-adjusted genotype-specific ORs were similar to unadjusted ORs (data not shown). There was no significant heterogeneity between strata for any of the SNPs studied (data not shown). There was a significant difference in genotype frequency in ovarian cancer cases and controls for rs2854344 (Pheterogeneity = 0.0015) and rs4151620 (Pheterogeneity = 0.0001). There was no evidence for heterogeneity between studies for these two SNPs (P = 0.19 and 0.62, respectively). There was no difference in genotype frequencies for the other nine SNPs. The best fitting model for rs4151620 was recessive; rare allele homozygote being at
80% reduced risk of ovarian cancer (OR, 0.19; 95% CI, 0.07-0.53; P = 0.00005). For rs2854344, the best fitting model was dominant, this model fitting the data slightly better than a codominant one (P = 0.0009 versus P = 0.005), with 27% reduced risk for minor allele carriers (OR, 0.73; 95% CI, 0.61-0.89). The risks were similar after exclusion of the additional controls from the SEARCH breast cancer study but less significant as would be expected with a smaller sample size. The OR (95% CI) for rs4151620 rare allele homozygote compared with common homozygote was 0.22 (0.08-0.62), and the OR (95% CI) for rs2854344 rare allele carriers was 0.77 (0.62-0.95). The two associated SNPs are only weakly correlated with each other (r2 = 0.07) and their effects seem to be independent. In a stepwise logistic regression, including all 11 SNPs, only the terms for the minor allele homozygotes of rs4151620 (P = 0.003) and the rs2854344 heterozygote (P = 0.001) remained in the final model.
|
|
|
|
| Discussion |
|---|
|
|
|---|
The SNPs tested here were chosen to efficiently tag all the known common variants in the gene and not for their likely functional effect. Our stSNP selection was based on EGP resequencing data, but only 38% of the gene had been resequenced; therefore common SNPs in the nonsequenced regions will have been missed. However, the HapMap project has genotyped over 300 SNPs in RB1 (phase II release no. 20). Of these, only 22 SNPs in a single linkage disequilibrium (LD) block were common (MAF >0.05) in the samples of north-western European ancestry by the centre d'Etude du polymorphisme Humain (CEPH). Nine HapMap SNPs were also in the EGP data set (5 stSNPs). The 5 stSNPs (i.e., rs1981434, rs520342, rs399413, rs4151551, and rs3092904) tagged the 22 HapMap SNPs with mean rp2 of 0.83 and 17 SNPs were captured with r2 > 0.8. The full set of 11 stSNPs would be expected to do better. stSNPs selected from HapMap phase II data provide good power to detect all common variation (21). We are therefore confident that our set of stSNPs has adequately tagged the known and unknown common variants within the gene.
We found a significant association for the minor allele of two SNPs, rs2854344 and rs4151620; both were associated with reduced EOC risks. Similarly, rs2854344 has also been reported to be associated with reduced breast cancer risk (14). Carriers of the minor allele of rs2854344 seems to have reduced breast cancer risk compared with common allele homozygote (OR, 0.86; 95% CI, 0.76-0.96). Nevertheless, these results ought to be interpreted with caution. The P-value presented above have not been adjusted for multiple hypothesis testing. Because SNPs within the same gene may be in LD, the test statistics for the 11 SNPs are not independent, and standard methods for adjusting for multiple testing, such as the Bonferroni correction, would be too conservative. We therefore used a simulation to determine an empirical P for the most significant result (i.e., Precessive = 0.00005 for rs4151620). In this analysis, we randomly shuffled the case-control status among individuals multiple times and estimated how frequently a P < 0.00005 was obtained from the randomly permuted data. This method also accounts for the testing of multiple genetic models with each SNP. In 10,000 permutations, a P < 0.00005 was observed on 11 occasions, giving the most significant P corrected for multiple testing of 0.0011. Nevertheless, the possibility that this result is a false positive remains substantial. One way of assessing this probability is to use the false-positive report probability (FPRP). The FPRP is the estimated probability that a specific result is a false positive. It depends on the prior probability of a true association, the observed level of significance (
), and the statistical power to detect the OR of the alternative hypothesis at the given
(22). As there is a very large number of common SNPs in the genome, the overall prior probability of association is very low (say 1 in a million). However, the prior is likely to be more favorable because we selected RB1 as a candidate based on known biology and one of the two associated SNPs (rs2854344) has also been reported to be associated with breast cancer with similar protective effect. The FRPRs for the two associated SNPs under different prior probabilities are presented in Table 5
. Hidden population stratification is an alternative explanation for a spurious association. This occurs when allele frequencies differ between population subgroups and case and controls are drawn differentially from those subgroups. It seems unlikely that population stratification is relevant in this association study because the cases and controls in the three studies reported here were drawn from the same ethnic groups. Furthermore, if stratification were present, it is unlikely that the same degree of stratification would be seen in all three studies.
|
In conclusion, we have genotyped 11 stSNPs tagging the known common variants in RB1 in three ovarian cancer case-control studies. We found some evidence of association for rs2854344 and rs4151620 with invasive ovarian cancer risk. The observed associations with ovarian cancer risk warrant confirmation in independent studies before further functional studies.
| 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 Joan MacIntosh, Hannah Munday, Barbara Perkins, Clare Jordan, Kristy Driver, Mitul Shah, the local general practices and nurses, and the East Anglian Cancer Registry for recruitment of the United Kingdom cases; the EPIC-Norfolk investigators for recruitment of the United Kingdom controls; Craig Luccarini and Don Conroy for expert technical assistance; and all the study participants who contributed to this research.
| Footnotes |
|---|
9 http://pga.gs.washington.edu/finished_genes.html. ![]()
10 http://www.broad.mit.edu/mpg/tagger/. ![]()
11 http://genome.ucsc.edu/cgi-bin/hgGateway. ![]()
12 http://pupasnp.bioinfo.cnio.es. ![]()
Received 6/16/06. Revised 7/18/06. Accepted 8/ 7/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. L. Goode, B. L. Fridley, R. A. Vierkant, J. M. Cunningham, C. M. Phelan, S. Anderson, D. N. Rider, K. L. White, V. S. Pankratz, H. Song, et al. Candidate Gene Analysis Using Imputed Genotypes: Cell Cycle Single-Nucleotide Polymorphisms and Ovarian Cancer Risk Cancer Epidemiol. Biomarkers Prev., March 1, 2009; 18(3): 935 - 944. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Song, E. Hogdall, S. J. Ramus, R. A. DiCioccio, C. Hogdall, L. Quaye, V. McGuire, A. S. Whittemore, M. Shah, D. Greenberg, et al. Effects of Common Germ-Line Genetic Variation in Cell Cycle Genes on Ovarian Cancer Survival Clin. Cancer Res., February 15, 2008; 14(4): 1090 - 1095. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Beesley, S. J. Jordan, A. B. Spurdle, H. Song, S. J. Ramus, S. K. Kjaer, E. Hogdall, R. A. DiCioccio, V. McGuire, A. S. Whittemore, et al. Association Between Single-Nucleotide Polymorphisms in Hormone Metabolism and DNA Repair Genes and Epithelial Ovarian Cancer: Results from Two Australian Studies and an Additional Validation Set Cancer Epidemiol. Biomarkers Prev., December 1, 2007; 16(12): 2557 - 2565. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Gayther, H. Song, S. J. Ramus, S. K. Kjaer, A. S. Whittemore, L. Quaye, J. Tyrer, D. Shadforth, E. Hogdall, C. Hogdall, et al. Tagging Single Nucleotide Polymorphisms in Cell Cycle Control Genes and Susceptibility to Invasive Epithelial Ovarian Cancer Cancer Res., April 1, 2007; 67(7): 3027 - 3035. [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 |