Analysis of Prostate Cancer Susceptibility Variants in South African Men: Replicating Associations on Chromosomes 8q24 and 10q11
Genome-wide association studies (GWAS) have implicated single nucleotide polymorphisms (SNPs) on chromosomes 2p15, 6q25, 7p15.2, 7q21, 8q24, 10q11, 10q26, 11q13, 17q12, 17q24, 19q13, and Xp11, with prostate cancer (PCa) susceptibility and/or tumour aggressiveness, in populations of African, European, and Asian ancestry. The objective of this study was to confirm these associations in South African Mixed Ancestry and White men. We evaluated 17 prioritised GWAS SNPs in South African cases (331 Mixed Ancestry and 155 White) and controls (178 Mixed Ancestry and 145 White). The replicated SNP associations for the different South African ethnic groups were rs7008482 (8q24) (), rs6983267 (8q24) (), and rs10993994 (10q11) () in Mixed Ancestry men and rs10993994 () in White men. No significant associations were observed for the analyses stratified by disease aggressiveness in the individual and the combined population group analysis. The present study demonstrates that a number of known PCa susceptibility variants may contribute to disease susceptibility in South African men. Larger genetic investigations extended to other South African population groups are warranted to confirm the role of these and other SNPs in disease susceptibility.
Prostate cancer (PCa) is the most frequently diagnosed cancer in men in developed countries, with higher incidence rates in the United States of America (USA) in individuals of African ancestry . African American (AA) men have a PCa incidence rate that is 1.6 times that in European American (EA) men and the mortality rate is 2.5 times greater than EA men . In developing countries, recent estimates suggest that PCa is also the leading cancer in terms of incidence and mortality in men from Sub-Saharan Africa (SSA) and the Caribbean . Moreover, underdiagnoses and underreporting have been shown to be widespread in men from SSA and the Caribbean , further indicating that incidence and mortality rates may be substantially underestimated in developing countries.
Numerous genome-wide association studies (GWAS) have identified close to 100 PCa susceptibility variants, primarily in populations of European ancestry [5–12], to a lesser extent in AA men [8, 12, 13], and, recently, in African men [14, 15]. Even though many of these associations exceeded or were close to genome-wide significance levels, it remains necessary to confirm these findings in multiple independent populations to rule out false positive results, thereby improving the likelihood that they represent true associations. To this end, several studies have replicated GWAS associations in AA, EA, and Asian men [16–29]. Although GWAS have generally been undertaken to a lesser extent in African descent populations, a recent meta-analysis of 82 PCa susceptibility variants identified primarily in GWAS in European and Asian descent populations generally demonstrated consistent replication of the direction of effects in men with African ancestry .
In the present study, we investigated PCa susceptibility single nucleotide polymorphisms (SNPs) in men from two separate South African population groups, one group of mixed Khoisan (San), African, European, and Asian descent and the other of European descent. We report replicating associations with PCa susceptibility for three of the 17 investigated SNPs among Mixed Ancestry (also termed South African Coloured) and White men.
2. Materials and Methods
2.1. Study Population
The study population comprised 331 Mixed Ancestry and 155 White cases (ethnicity was self-reported for all participants, as previously described) . All cases had no prior cancer at any other site and were recruited from the Division of Urology, Tygerberg Hospital (Cape Town, South Africa). Clinical characteristics including prostate specific antigen (PSA), Gleason score, tumour node metastasis (TNM) stage, age at diagnosis, and other cancer diagnoses were obtained from medical records. All cases underwent radical prostatectomy, transurethral resection of the prostate, or prostatic biopsy and had histologically confirmed PCa. Cases were additionally stratified by risk of developing aggressive PCa using the D’Amico Cancer Risk Group Classification system. Probability of developing aggressive PCa in cases was categorized as “low” (tumour category < T1-2a or PSA ≤ 10 ng/mL or Gleason score < 7) or “high” (tumour category ≥ T2b or PSA > 10 or Gleason score ≥ 7) (“high” was defined here as the combination of intermediate + high risk).
Controls were selected among subjects admitted to Tygerberg Hospital for routine PSA examinations or benign prostatic hyperplasia (BPH) and comprised of 178 Mixed Ancestry men and 145 White men, matched for age and ethnicity (self-reported), and were from the same geographical region. Inclusion criteria for controls were PSA levels of ≤4.0 ng/mL and a normal digital rectal examination (DRE) or negative biopsy on histology.
All participants were informed and gave written consent to participate in the study and allow their biological samples to be genetically analysed, conducted in accordance with the Declaration of Helsinki. The study was approved by the Stellenbosch University Faculty of Medicine and Health Sciences Human Research Ethics Committee.
2.2. SNP Selection and Genotyping
Based on previously reported PCa GWAS and replication studies, we selected 17 SNPs on chromosomes 2p15, 6q25, 7p15.2, 7q21, 8q24, 10q11, 10q26, 11q13, 17q12, 17q24, 19q13, and Xp11. Given the diverse ancestral lineages of the South African population, the SNPs were prioritised for analyses primarily if they had been implicated in PCa susceptibility in GWAS and/or replication studies in European, African, and Asian descent populations. At the time of initial SNP selection, the novel PCa susceptibility variants identified in African men [14, 15] had not yet been reported and were thus not investigated in this study.
Genotyping was performed at the Stellenbosch University Central Analytical Facility (CAF) (Cape Town, South Africa) on the 17 selected SNPs using commercially available Taqman SNP Genotyping Assays (Life Technologies, USA). Quality control (QC) for exclusion from the study for individual samples was a genotype call rate <90%. For all SNPs, genotyping QC was assessed by including DNA duplicates (to an equivalent of 5% duplicates per 384-well plate), and the concordance genotype call rate between samples and duplicates was set at ≥95%. We excluded all SNPs whose genotype frequency departed from Hardy-Weinberg equilibrium (HWE) at in the controls and SNPs that were monomorphic (had a frequency <1%) in our study populations.
2.3. Statistical Analyses
Descriptive analyses for discrete traits were carried out using contingency table methods and Fisher’s exact tests (FET) or chi-square statistics. Medians were used to summarize continuously distributed traits. We tested the SNPs in the controls for departures of HWE and determined allele and genotype frequencies in cases and controls. For each SNP in each population group, logistic regression models were used to estimate PCa genotype and allelic odds ratios (OR), 95% confidence intervals (CI), and corresponding values adjusted for age (i.e., age at diagnosis for cases and age at ascertainment for controls). We adjusted for age and ethnicity in the combined population logistic regression analysis. Logistic regression was similarly used to estimate OR, 95% CI, and values for cases stratified by disease aggressiveness (“low” versus “high”). All analyses were performed in R, freely available from http://www.R-project.org/. In the present study, we adopted a conservative approach by only considering a Bonferroni corrected value threshold of statistically significant ( SNPs × two population groups). This correction may be considered to be overly conservative ; however, because the number of multiple comparisons we performed was relatively small, we opted to use a family-wise error rate approach instead of a less stringent approach .
The median age at diagnoses of PCa cases and the median age of controls were significantly different among the groups; PSA levels were significantly different between the groups, with Mixed Ancestry controls and cases showing higher median levels (Table 1). No significant differences were noted for Gleason score and tumour stage between the two groups, although Mixed Ancestry men showed significantly higher frequency of distant metastasis at the time of diagnosis (Table 1).
For the genotype QC analysis, two cases and one control Mixed Ancestry sample achieved a genotype call rate <90% and were thus excluded, leaving the analysed total of cases at 329 Mixed Ancestry and 155 White and a total number of controls at 177 Mixed Ancestry and 145 White. Additionally, four of the 17 selected SNPs were excluded from further study: rs721048 (2p15) and rs5945619 (Xp11) failed to meet the QC concordance genotype call rate of ≥95%, and Bd11934905 (8q24) and rs7210100 (17q21) were monomorphic in both population groups.
Of the 13 remaining prioritised SNPs, three were significantly associated with PCa after correction for multiple testing in Mixed Ancestry men: rs7008482 (8q24) (OR = 0.41, 95% CI: 0.26–0.62; ), rs6983267 (8q24) (OR = 0.35, 95% CI: 0.23–0.53; ), and rs10993994 (10q11) (OR = 1.56, 95% CI: 1.56–6.08; ) (Table 2). In White men, only rs10993994 (OR = 3.49, 95% CI: 2.28–5.54; ) was significantly associated with PCa after correction for multiple testing (Table 2). In the combined population analysis, the associations for rs7008482, rs6983267, and rs10993994 and PCa susceptibility were statistically significant (Table 2). No significant associations were observed for the case-only analyses stratified by disease aggressiveness (“low” versus “high” risk) in the individual and the combined population group analysis (data not shown).
The present study replicates previously reported associations for rs7008482, rs6983267, and rs10993994, in two different South African population groups. Men from the two population groups reported here have PCa age-standardized incidence rates that are between 1.4 and 2.4 times higher than men from other South African population groups (Mixed Ancestry: 39.1 per 100,000 and White: 52.0 per 100,000 compared to Asian: 28.6 per 100,000 and Black: 21.6 per 100,000) (National Cancer Registry, 2009) . When we compared the age at diagnosis between the two groups, Mixed Ancestry men were shown generally to be diagnosed with PCa at an earlier age than White men (67 versus 70 years). Additionally, Mixed Ancestry men, who have substantial proportions of African ancestry (combination of San, European, Black, and Asian ancestry ), showed significantly higher frequency of distant metastasis at the time of diagnosis. A few studies have suggested that men with African ancestry more often present with aggressive disease compared to men of European ancestry [29, 46]. Though highly speculative, our observations may support this view, although our genetic analysis stratified by disease aggressiveness failed to provide supporting evidence.
Most GWAS and replication studies have consistently identified genetic associations with SNPs on chromosome 8q24 [5–8, 10, 11, 16, 19, 20, 22, 23, 25, 26]; two 8q24 SNPs were associated in our study, whereas in another study South African Black men showed associations with the rs6983561 . It has been suggested that multiple, independent polymorphic variants on chromosome 8q24 may produce a common biological mechanism that contributes to disease, or, alternatively, the 8q24 regions may cumulatively influence the regulation of adjacent genes (cis-regulation) or genes on other chromosomes (trans-regulation) . Additionally, it has been suggested that the 8q24 regions harbor previously unannotated microRNAs (miRNAs) and that the miRNA elements in these regions are involved in cis-regulation of distal genes, thereby affecting RNA expression . However, it remains inconclusive why so many studies have shown PCa susceptibility associations to 8q24 given that it is a gene-poor region.
We showed that the SNP rs10993994 was associated with PCa susceptibility across both population groups, as well as in the combined populations. Pomerantz et al.  demonstrated that the T-allele of rs10993994 decreased expression in normal tissue of both isoforms of the beta-microseminoprotein gene (MSMB), as well as increasing the expression of the five isoforms of the nuclear receptor coactivator 4 gene (NCOA4). Furthermore, it was speculated that because rs10993994 is not a susceptibility allele for other cancers, such as colon or breast cancer, the associations observed with this SNP may thus be specific to prostate cancer .
In the present study, we opted to use a very conservative approach to reduce the chance of identifying false positive associations. However, we noted SNPs that would conventionally be considered associated with the disease if a less stringent statistical approach was employed. In the White group, rs1048657 (7p15), rs4430796 (17q12), and rs1859962 (17q24) were significant at but not after correction for multiple testing (Table 2). For the Mixed Ancestry, because of the effect of population stratification, we only consider a statistically significant (also applied when combining the Mixed Ancestry with other groups) . For the Mixed Ancestry and the combined group, only rs7008482, rs6983267, and rs10993994 were significantly associated after correction for population stratification as well as for multiple testing (Table 2). The individual samples sizes in our study were small; therefore, larger genetic investigations in the population groups investigated here are required to conclusively establish whether or not the selected SNPs play a role in PCa susceptibility.
Our study findings should be interpreted in light of some limitations. We acknowledge that our investigations were only performed on men from the two South African population groups with the highest incidence of PCa and recognize that significant underreporting and underdiagnosis of PCa have been documented in other groups, particularly in South African Black men ; thus, future studies would need to be extended to other population groups. Our study sample was relatively small and this may have disadvantaged our power to detect statistical significance. We excluded two SNPs (Bd11934905 and rs7210100) with minor allele frequencies (MAFs) < 1%, because there may not be sufficient power to detect the effect of rare variants due to our sample size; however, larger studies should be undertaken as we cannot disregard them as potential disease-causative SNPs. Population stratification is an additional confounder in case-control association analyses. Daya and colleagues  have shown using ancestry informative markers (AIMs) specific for the South African Mixed Ancestry population that genetic associations may change after adjustment for ancestry (previous associations may be abolished or associations may only become apparent after adjustment for ancestry). However, the authors also noted that association results are more likely to be unaffected by adjustment for ancestry if the allele frequencies in the stratified population do not differ markedly from the source populations . When we compared the allele frequencies of the selected SNPs in the Mixed Ancestry to likely source populations in the 1000 Genomes database (African, European, South, and East Asian), the allele frequencies we observed often differed by <5%.
We demonstrate that a number of reported susceptibility SNP variants may contribute to PCa susceptibility in South African Mixed Ancestry and White men. However, this was a relatively small study and only men from two population groups were investigated. Consequently, larger genetic investigations extended to different South African population groups are warranted to confirm the role of these and other SNPs in PCa susceptibility.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors would like to thank all the participants of the study and the urologists who facilitated the collection of samples. The study was supported by grants from the AORTIC/NCI BIG_CAT Programme and Cancer Association of South Africa (CANSA) (424366) (Pedro Fernandez).
R. A. Eeles, Z. Kote-Jarai, A. A. Al Olama et al., “Identification of seven new prostate cancer susceptibility loci through a genome-wide association study,” Nature Genetics, vol. 41, no. 10, pp. 1116–1121, 2009.View at: Google Scholar
F. R. Schumacher, S. I. Berndt, A. Siddiq et al., “Genome-wide association study identifies new prostate cancer susceptibility loci,” Human Molecular Genetics, vol. 20, no. 19, pp. 3867–3875, 2011.View at: Google Scholar
R. A. Eeles, A. A. Olama, S. Benlloch et al., “Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array,” Nature Genetics, vol. 45, no. 4, pp. 385–392, 2013.View at: Google Scholar
J. Shan, K. Al-Rumaihi, D. Rabah et al., “Genome scan study of prostate cancer in Arabs: identification of three genomic regions with multiple prostate cancer susceptibility loci in Tunisians,” Journal of Translational Medicine, vol. 11, no. 1, article 121, 2013.View at: Publisher Site | Google Scholar
N. J. Camp, J. M. Farnham, J. Wong, G. B. Christensen, A. Thomas, and L. A. Cannon-Albright, “Replication of the 10q11 and Xp11 prostate cancer susceptibility variants: results from a Utah pedigree-based study,” Cancer Epidemiology Biomarkers and Prevention, vol. 18, no. 4, pp. 1290–1294, 2009.View at: Publisher Site | Google Scholar
C. Robbins, J. B. Torres, S. Hooker et al., “Confirmation study of prostate cancer susceptibility variants at 8q24 in African Americans identifies a novel susceptibility locus,” Genome Research, vol. 17, no. 12, pp. 1717–1722, 2007.View at: Google Scholar
J. Xu, A. S. Kibel, J. J. Hu et al., “Prostate cancer susceptibility associated loci in African Americans,” Cancer Epidemiology Biomarkers and Prevention, vol. 18, no. 7, pp. 2145–2149, 2009.View at: Google Scholar
E. A. Tindall, M. R. Bornman, S. Van Zyl et al., “Addressing the contribution of previously described genetic and epidemiological risk factors associated with increased prostate cancer risk and aggressive disease within men from South Africa,” BMC Urology, vol. 13, article 74, 2013.View at: Publisher Site | Google Scholar
P. Fernandez, C. M. Zeigler-Johnson, E. Spangler et al., “Androgen metabolism gene polymorphisms, associations with prostate cancer risk and pathological characteristics: a comparative analysis between South African and Senegalese men,” Prostate Cancer, vol. 2012, Article ID 798634, 8 pages, 2012.View at: Publisher Site | Google Scholar
P. H. Westfall and S. S. Young, Resampling-Based Multiple Testing, Wiley, New York, NY, USA, 1993.
J. J. Chen, P. K. Roberson, and M. J. Schell, “The false discovery rate: a key concept in large-scale genetic studies,” Cancer Control, vol. 17, no. 1, pp. 58–62, 2010.View at: Google Scholar
H.-Y. Yeh, S.-W. Cheng, Y.-C. Lin, C.-Y. Yeh, S.-F. Lin, and V.-W. Soo, “Identifying significant genetic regulatory networks in the prostate cancer from microarray data based on transcription factor analysis and conditional independency,” BMC Medical Genomics, vol. 2, article 70, 2009.View at: Publisher Site | Google Scholar
C. Debiais-Delpech, J. Godet, N. Pedretti et al., “Expression patterns of candidate susceptibility genes HNF1β and CtBP2 in prostate cancer: association with tumor progression,” Urologic Oncology: Seminars and Original Investigations, vol. 32, no. 4, pp. 426–432, 2014.View at: Publisher Site | Google Scholar
M. Brecht, A. C. M. Steenvoorden, S. Luf, C. R. Bartram, and J. W. G. Janssen, “Rearrangement and expression of myeov and hst in NIH/3T3 transfectants: a caveat for the interpretation of DNA transfection analyses,” Oncology Reports, vol. 17, no. 5, pp. 1127–1131, 2007.View at: Google Scholar
South African National Cancer Registry (NCR) Reports, http://www.nioh.ac.za/assets/files/NCR_2009_FINAL.pdf.
C. R. Ritch, B. F. Morrison, G. Hruby et al., “Pathological outcome and biochemical recurrence-free survival after radical prostatectomy in African-American, Afro-Caribbean (Jamaican) and Caucasian-American men: an international comparison,” BJU International, vol. 111, no. 4, pp. E186–E190, 2013.View at: Publisher Site | Google Scholar
N. Ahmadiyeh, M. M. Pomerantz, C. Grisanzio et al., “8q24 prostate, breast, and colon cancer risk loci show tissue-specific long-range interaction with MYC,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 21, pp. 9742–9746, 2010.View at: Publisher Site | Google Scholar