Abstract
The Von Hippel-Lindau (VHL) tumor suppressor gene is a crucial regulator of the hypoxia response pathway and plays an important role in tumorigenesis, particularly in tumor growth and vascularization. We hypothesize that polymorphisms in the functional region of VHL may influence susceptibility to prostate cancer (PCa). We genotyped a potentially functional polymorphism (rs779805) in 5′ UTR region of VHL in a case-control study of 665 PCa patients and 715 cancer-free controls in a Chinese population using the Taqman assay. The genetic associations between the incidence and progression of PCa were assessed by logistic regression. We observed that the rs779805 A>G polymorphism was significantly associated with risk for PCa. Compared with the AA genotype, the AG and AG/GG genotypes were associated with decreased risk of PCa (adjusted odds ratio [OR]=0.79, 95% confidence interval [CI]=0.62–0.99, and adjusted OR=0.76, 95% CI=0.61–0.95, respectively). Further, this decreased risk was more pronounced in the subgroups of nonsmokers (OR=0.73, 95% CI=0.54–0.98), nondrinkers (OR=0.70, 95% CI=0.54–0.91) and patients without family history of cancer (OR=0.72, 95% CI=0.57–0.92). In addition, the decreased risk associated with rs779805 variant genotypes (AG/GG) was more pronounced among the prostate specific antigen (PSA)>20 ng/mL subgroup (OR=0.68, 95% CI=0.49–0.95). Our findings suggest that the rs779805 A>G polymorphism in VHL may confer susceptibility to PCa in the Chinese population.
Introduction
Intratumoral hypoxia is a hallmark of the solid cancers (Hill et al., 2009). It has been shown to be correlated with increased tumor pathogenesis, invasiveness, and metastasis (Stewart et al., 2010). The VHL gene known to act is the causal gene for VHL syndrome disease (Latif et al., 1993), which is characterized by the risk factors of a number of different carcinomas in vivo and in vitro (Kim et al., 2010). The VHL gene product, VHL protein (pVHL), interacts with a number of cellular proteins and is implicated in the control of angiogenesis, extracellular matrix formation, cell metabolism, and mitogenesis (Kaelin, 2005). However, the most well-characterized function of pVHL is its role as the substrate recognition submit of an E3 ubiquitin ligase complex and as a master regulator of hypoxia-inducible factor (HIF) activity by targeting the HIF-alpha subunit for ubiquitination and rapid proteasomal degradation under normoxic conditions (Haase, 2009). HIF-1α is rapidly degraded by ubiquitination through the pVHL in the proteasome. Previous studies have demonstrated that the VHL/HIF-1 pathway plays a significant role in tumor progression and metastasis through activation of various genes linked to regulation of angiogenesis, cell survival, energy metabolism and apoptotic and proliferative responses in vivo and in vitro (Semenza, 2003; Vaupel, 2004). With regard to PCa, upregulation of HIF-1 activity has been shown to be an early event in prostate carcinogenesis in animal model studies (Huss et al., 2001). In PCa cell lines, expression of the HIF-1α protein is positively associated with cell growth rate and metastatic potential (Zhong et al., 1998, 2000). A recent immunohistochemical study illustrated that HIF-1α protein levels were upregulated in most prostatic intraepithelial neoplasia and further enhanced in primary and metastatic PCa (Du et al., 2003; Zhong et al., 2004).
Recently, several genetic association studies have suggested that polymorphisms in HIF1A (encoding HIF-1alpha) influence the susceptibility of PCa, although the results remained inconclusive (Chau et al., 2005; Li et al., 2007; Orr-Urtreger et al., 2007; Foley et al., 2009). We hypothesized that VHL, a crucial regulator of HIF-1, may have potentially functional polymorphisms that might play a role in the development of PCa. Many epidemiological studies have investigated the association between genetic variations in VHL and susceptibility and progression of various cancers, including pheochromocytoma (Waldmann et al., 2009), gastric cancer (Cao et al., 2008), and ovarian cancer (Goode et al., 2011). In our previous study, we observed no association between this SNP (rs779805) and renal cancer (Qin et al., 2011). To the best of our knowledge, there is no epidemiological report regarding the role of any of the VHL gene polymorphisms in PCa in the Chinese population. VHL SNPs were selected based on HapMap data (
Materials and Methods
Study subjects
The analysis included 665 patients with diagnosed and histologically confirmed PCa and 705 cancer-free control subjects. All subjects were genetically unrelated and age matched (mean±5 years) ethnic Han Chinese recruited in an ongoing study that began in September 2003 and lasted until January 2010 at The First Hospital of Nanjing Medical University, Nanjing, China. Before recruitment, all subjects were interviewed face-to-face to collect individual demographic data and exposure information, including age, race, history of tobacco and alcohol use, and family history of cancer. Those subjects who had been smoking daily for more than 1 year were defined as smokers and the rest as nonsmokers. Individuals who had imbibed alcohol at least three times per week for at least 6 months prior to the study were defined as drinkers, and the others were considered as nondrinkers. Family history of cancer was defined as any cancer in first-degree relative (parents, siblings, or children). Disease stage was determined by pathologic findings, pelvic computed tomography, magnetic resonance imaging, and radio-nucleotide bone scans. All cases were classified according to World Health Organization criteria. Clinical stage was divided into localized and advanced cancer based on the TNM classification system promulgated by the American Joint Committee on Cancer (Wittekind et al., 2002). Localized PCa was defined as cancer that was detectable upon clinical examination but had not spread outside the prostate (T1–2N0M0). Advanced cancer was defined as the cancer had spread through the prostatic capsule (T3–4NXMX or TXN1MX or TXNXM1). The Gleason score was evaluated by pathologists working at the hospital using the Gleason scoring system (Epstein et al., 2005). Prostate-specific antigen (PSA) is a protein produced by cells of the prostate gland. Based on the EAU Guidelines on PCa and D'Amico's Risk-Based management of PCa, serum PSA value was dimidiated into two groups PSA>20 ng/mL and PSA≤20 ng/mL (Heidenreich et al., 2008; D'Amico, 2011), predicting the likelihood of progression of PCa. After completion of the interview, each subject provided written informed consent and donated 5 mL of blood for genomic DNA extraction. The research protocol was approved by the institutional review board of Nanjing Medical University.
Polymorphism genotyping
The rs779805 A>G polymorphism was genotyped using predesigned TaqMan MGB SNP genotyping assays (Applied Biosystems, Foster City, CA). The primer sequences for the rs779805 A>G polymorphism was as follows: forward primer, 5′ GGAAATACAGTAACGAGTTGGCCTA 3′ and reverse primer, 5′ CGCTCGGTAGAGGATGGAAC 3′. Amplification was performed by using TaqMan Universal Master Mix, 80X SNP Genotyping AssayMix, Dnase-free water, and 10-ng genomic DNA in a final volume of 5 μL per reaction, under the following conditions: with an initial 2 min at 50°C, 10 min at 95°C, followed by 40 cycles of 95°C for 15 sec, and 60°C for 1 min. According to the manufacturer's instructions, amplifications and analysis were performed in the 384-well ABI 7900HT Real Time PCR System (Applied Biosystems) and the Sequence Detection Systems software (SDS 2.3; Applied Biosystems) were used to automatically collect and analyze the data and to generate the genotype calls. Four negative controls were included in each plate to ensure accuracy of the genotyping. For quality control, genotyping was performed by two people independently in a blinded fashion. About 5% of the samples were randomly selected for repeated genotyping, and the results were 100% concordant.
Statistical analysis
The χ2 test (for categorical variables) and the student's t-test (for continuous variables) were used to compare the differences in frequency distributions of selected demographic variables, smoking status, alcohol use, family history of cancer in addition to each allele and genotype of the VHL (rs779805) polymorphism between the cases and controls. The associations between the VHL polymorphism and risk of PCa were estimated by ORs at 95% CIs from unconditional univariate and multivariate logistic regression models with the adjustment for age, smoking status, drinking status, and family history of cancer. The statistical power was calculated by using the PS software (
Results
Characteristics of the study population
The frequency distributions of selected characteristics of the cases and controls are presented in Table 1. There was no significant age difference between the cases and controls (p=0.727). However, there were more smokers and drinkers among the cases than among the controls, and these differences were found to be statistically significant (p=0.023 and p=0.049 for smoking and drinking status, respectively). Further, the frequency of first degree relatives with cancer among the cases was higher than among the controls (p<0.001). Within the PCa group, approximately 58.8% of patients were at the localized stage and 41.2% of patients were in advanced stage; 394 (59.1%) of patients had a high serum PSA level (over 20 ng/mL).
t-test for age distributions between the cases and controls; two-sided χ2 test for others selected variables between the cases and controls.
Association between VHL rs779805 polymorphism and risk of PCa
Allele frequencies and genotype distributions of VHL rs779805 polymorphism are shown in Table 2. All observed genotype frequencies in both cases and controls conformed to Hardy-Weinberg equilibrium (p=0.36, and p=0.37 for cases and controls respectively). The frequencies of the AA, AG, and GG genotypes in cases were 68.4%, 28.0%, and 3.6%, respectively, and 62.2%, 32.6%, and 5.2%, respectively, in the controls (p=0.042). When we used the most common genotype (rs779805 AA) as a reference, we found that AG and AG/GG genotypes were associated with a statistically significantly decreased risk of PCa (adjusted odds ratio [OR]=0.79, 95% confidence interval [CI]=0.62–0.99 and adjusted OR=0.76, 95% CI=0.61–0.95, respectively). Further, the rs779805 G allele frequency was 17.6% in the cases and 21.6% in the controls. The difference was statistically significant (p=0.012).
Two-sided χ2 test for the either genotype distributions or allele frequencies between the cases and controls.
Adjusted for age, smoking status, drinking status, and family history of cancer in logistic regression model.
OR, odds ratio; CI, confidence interval.
Stratified analysis of rs779805 polymorphism and risk of PCa
We then evaluated the effect of the rs779805 polymorphism on PCa risk stratified by age, smoking status, drinking status, and family history of cancer. As shown in Table 3, the association between rs779805 polymorphism and the risk of PCa appeared stronger in the subgroups of nonsmokers (OR=0.73, 95% CI=0.54–0.98), nondrinkers (OR=0.70, 95% CI=0.50–0.91), and individuals without any family history of cancer (OR=0.72, 95% CI=0.57–0.92). However, the association between rs779805 polymorphism and PCa risk did not vary by age. Taken together, the rs779805 polymorphism is an independent risk factor for PCa in the Chinese population with the potential confounding factors having no significant modification.
Two-sided χ2 test for the distributions between the cases and controls.
Adjusted for age, smoking status, drinking status, and family history of cancer in logistic regression model.
Association between VHL polymorphism and clinicopathological characteristics of PCa patients
We also assessed the effects of the AG/GG genotypes on the severity of PCa, and investigated the association between the rs779805 polymorphism and clinicopathological characteristics of the PCa patients. As shown in Table 4, the genotypes of AG/GG occurred less frequently in patients with PSA≥20 ng/mL (adjusted OR=0.68, 95% CI=0.49–0.95). There were no other significant associations between the genotypes and clinical stage and Gleason score among the PCa patients.
The p value and ORs were calculated and adjusted for age, smoking, drinking, and family history of cancer in logistic regression model.
Localized: T1–2N0M0; Advanced: T3–4NxMx or TxN1Mx or TxNxM1. Clinical staging according to the international TNM system for PCa.
Discussion
In the current study, we assessed the association between the rs779805 A>G polymorphism and the risk of occurrence and progression of PCa in the Chinese population. Compared with individuals with the rs779805 A allele, we observed that individuals carrying G allele were significantly associated with the decreased risk for PCa, and individuals with variant genotypes were associated with low serum PSA (PSA≤20 ng/mL). To the best of our knowledge, this is the first study on the genetic association between the VHL polymorphism and the risk of PCa to confirm the initial hypothesis that genetic variations in VHL play a role in the pathogenesis of this malignancy.
Although hypoxia is emerging as a common feature of PCa and the development of hypoxia-targeted therapies is becoming more common, their clinical potential still appears to have been underestimated. The VHL/HIF-1α pathway plays a central role in adaptation to hypoxic environments. Previous investigators have observed that high HIF-1α levels were associated with tumor progression and poor prognosis (Lu and Kang, 2010). Zhong et al. (1999) suggested that HIF-1α was upregulated in primary PCa relative normal prostate epithelium. They confirmed that the overexpression of HIF-1α is an early event in prostate carcinogenesis (Zhong et al., 2004). High levels of HIF-1α protein have many possible causes. One possible explanation is aberration of VHL expression, which can lead to overexpression of HIF-1. Accumulation studies have suggested that SNPs in the functional region may influence the expression and functions of these genes (Putra et al., 2011; Theodoropoulos et al., 2011). Additionally, according to the web-based SNP analysis tool, TFSEARCH 1.3, we speculate that the rs779805 SNP might disrupt the predicted transcription factors binding site of c-Myb. Therefore, SNPs in functional region of VHL might also impact the expression of VHL. In the present study, we observed individuals carrying G allele were significantly associated with decreased risk for PCa. Given the location of this polymorphism and in silico predictions, the genetic association between rs779805 A>G and PCa is biologically plausible. One possible mechanism is that the A to G substitution of this polymorphism may increase the transcription and expression of VHL, which then plays a protective role in PCa carcinogenesis by downregulating expression of HIF-1. However, according to the HapMap database, it should be noted that the other two polymorphisms of VHL (rs779803 and rs1642742) are in complete linkage disequilibrium with rs779805. An in silico analysis of the VHL sequence predicted that rs779803 G>A polymorphism in the 5′ region of VHL generated a latent seed site for oct-1. These polymorphisms may also individually or jointly influence the expression of VHL by altering the affinity of transcription factor binding to the promoter or mRNA stability, which may further affect the expression level of HIF-1α.
In the current study, our results showed that smoking and drinking status were risk factors for PCa. In stratified analysis, we found that the decreased risk associated with VHL polymorphism was more pronounced in nonsmokers and nondrinkers, suggesting genetic effects might be overwhelmed by environmental effects. We also observed the effect of rs779805 polymorphism on the risk of PCa was more predominant among people without family history of cancer, suggesting that the polymorphism effect might be overwhelmed by the effect of other inherited genetic factors.
In addition, we did not find any association between polymorphism and pathological stage of PCa. However, we did observe that rs779805 AG./GG genotypes had a significant protective role in the patients with PSA over 20 ng/mL. The finding may arise from the significant relationship between serum PSA and hypoxia (Pipinikas et al., 2008). The clinical link between PSA levels and hypoxic fraction remains to be determined. One plausible mechanism is that hypoxia may enhance PSA expression, and increase the transcriptional activity of androgen receptors and PSA through HIF (Park et al., 2006; Horii et al., 2007; Mitani et al., 2010). Taken together, our results indicate that there might be different mechanisms underlying the early development and subsequent progression of PCa, and the SNP (rs779805) in the VHL gene may contribute to the progression of PCa.
Some limitations of the study should be addressed. First, our sample size is moderate and the statistical power of the study is limited, especially regarding stratification and interaction analysis. However, our post hoc power calculation suggests that we had an 81% power to detect an OR of 1.42 or greater and a minimal OR of 0.68 with an exposure frequency of 22% at the current sample size. Second, our controls were recruited from the hospital, so we could not exclude the inherent selection bias for any particular genotype. However, the genotype distributions of the controls in our study were similar to distributions reported to as typical for Asian population in other studies and conformed to Hardy-Weinberg Equilibrium (Droma et al., 2008; Qin et al., 2011). Third, ethnic variation is known to be a major risk for PCa. We found that genotype frequencies of rs779805 were different from the data population-wide genotype distribution reported in the NCBI dbSNP database. In our study with 715 control subjects, the frequencies of the AA, AG, and GG genotypes were 62.2%, 32.6%, and 5.2%, respectively, compared with 42.3%, 47.7%, and 9.9% in ss1596002 European samples (HapMap-CEU, n=222). Although the vital association between VHL and susceptibility to PCa came from a biologically based a priori hypothesis, our findings should be independently verified in other populations with a higher incidence, particularly in Caucasian. Last, the lack of detailed information on risk factors such as dietary habits, occupational exposures, and physical activity may limit the statistical power of our research. Therefore, a well-designed study involving the collection of more detailed information from more populations is warranted to confirm our finding.
Conclusions
Our present study revealed the association between the rs779805 polymorphism of VHL and susceptibility to PCa. We found that the variant G allele of rs779805 is associated with decreased risk of PCa. These findings suggest the rs779805 polymorphism may be a marker for genetic susceptibility to PCa. Further study is needed to investigate the relationship between VHL polymorphisms and the regulation of expression and function of VHL. A survey of the interactions between genetic and nongenetic risk factors including dietary habit and environmental factors must be performed to evaluate their contribution to the risk of PCa in the Chinese population.
Footnotes
Acknowledgments
We thank all the people who helped in this study. This study was partly supported by the Program for Development of Innovative Research Team in the First Affiliated Hospital of Nanjing Medical University, Provincial Initiative Program for Excellency Disciplines, Jiangsu Province and the National Natural Science Foundation of China (81171963).
Disclosure Statement
No competing financial interests exist.
