Genetic Association of the KLK4 Locus with Risk of Prostate Cancer

Abstract

The Kallikrein-related peptidase, KLK4, has been shown to be significantly overexpressed in prostate tumours in numerous studies and is suggested to be a potential biomarker for prostate cancer. KLK4 may also play a role in prostate cancer progression through its involvement in epithelial-mesenchymal transition, a more aggressive phenotype, and metastases to bone. It is well known that genetic variation has the potential to affect gene expression and/or various protein characteristics and hence we sought to investigate the possible role of single nucleotide polymorphisms (SNPs) in the KLK4 gene in prostate cancer. Assessment of 61 SNPs in the KLK4 locus (±10 kb) in approximately 1300 prostate cancer cases and 1300 male controls for associations with prostate cancer risk and/or prostate tumour aggressiveness (Gleason score <7 versus ≥7) revealed 7 SNPs to be associated with a decreased risk of prostate cancer at the P_trend<0.05 significance level. Three of these SNPs, rs268923, rs56112930 and the HapMap tagSNP rs7248321, are located several kb upstream of KLK4; rs1654551 encodes a non-synonymous serine to alanine substitution at position 22 of the long isoform of the KLK4 protein, and the remaining 3 risk-associated SNPs, rs1701927, rs1090649 and rs806019, are located downstream of KLK4 and are in high linkage disequilibrium with each other (r²≥0.98). Our findings provide suggestive evidence of a role for genetic variation in the KLK4 locus in prostate cancer predisposition.

About the Authors:

Felicity Lose

Affiliations Molecular Cancer Epidemiology Group, Genetics and Population Health Division, Queensland Institute of Medical Research, 300 Herston Rd, Herston, Brisbane, Queensland, Australia, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Srilakshmi Srinivasan

Affiliation: Australian Prostate Cancer Research Centre – Queensland, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Tracy O’Mara

Affiliations Molecular Cancer Epidemiology Group, Genetics and Population Health Division, Queensland Institute of Medical Research, 300 Herston Rd, Herston, Brisbane, Queensland, Australia, Australian Prostate Cancer Research Centre – Queensland, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Louise Marquart

Affiliation: Statistics Unit, Queensland Institute of Medical Research, Brisbane, Queensland, Australia

Suzanne Chambers

Affiliations Griffith Health Institute, Griffith University, Brisbane, Queensland, Australia, Viertel Centre for Cancer Research, Cancer Council Queensland, Brisbane, Queensland, Australia, University of Queensland Centre for Clinical Research, Royal Brisbane Hospital, Brisbane, Queensland, Australia

Robert A. Gardiner

Affiliation: University of Queensland Centre for Clinical Research, Royal Brisbane Hospital, Brisbane, Queensland, Australia

Joanne F. Aitken

Affiliations Griffith Health Institute, Griffith University, Brisbane, Queensland, Australia, Viertel Centre for Cancer Research, Cancer Council Queensland, Brisbane, Queensland, Australia

the Australian Prostate Cancer BioResource7YeadonTrinaInstitute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, AustraliaSaundersPamelaInstitute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, AustraliaEckertAllisonInstitute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, AustraliaClementsJudithInstitute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, AustraliaHeathcotePeterBrisbane Urology Clinic, Brisbane Urology Clinic Central Queensland Urology Clinic, Wickham Terrace, Brisbane, Queensland, AustraliaWoodGlennBrisbane Urology Clinic, Brisbane Urology Clinic Central Queensland Urology Clinic, Wickham Terrace, Brisbane, Queensland, AustraliaMaloneGregBrisbane Urology Clinic, Brisbane Urology Clinic Central Queensland Urology Clinic, Wickham Terrace, Brisbane, Queensland, AustraliaSamaratungaHemaAquesta Pathology, Toowong, Queensland, AustraliaCollinsAngusSullivan and Nicolaides Pathology, Brisbane, Queensland, AustraliaTurnerMeganSullivan and Nicolaides Pathology, Brisbane, Queensland, AustraliaKerrKrisSullivan and Nicolaides Pathology, Brisbane, Queensland, Australia

¶Membership of the Australian Prostate Cancer BioResource is provided in the Acknowledgments.

Affiliation: Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Amanda B. Spurdle

Affiliations Molecular Cancer Epidemiology Group, Genetics and Population Health Division, Queensland Institute of Medical Research, 300 Herston Rd, Herston, Brisbane, Queensland, Australia, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Jyotsna Batra

Contributed equally to this work with: Jyotsna Batra, Judith A. Clements

* E-mail: [email protected]

Affiliation: Molecular Cancer Epidemiology Group, Genetics and Population Health Division, Queensland Institute of Medical Research, 300 Herston Rd, Herston, Brisbane, Queensland, Australia

Judith A. Clements

Contributed equally to this work with: Jyotsna Batra, Judith A. Clements

Affiliation: Australian Prostate Cancer Research Centre – Queensland, Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia

Introduction

The Kallikrein (KLK) gene family consists of 15 genes in a tightly clustered locus over 320 kilobases (kb) at 19 q13.4 [1]. Many of the KLKs display altered expression in disease, in particular hormone-dependent cancers [1], [2]. KLK4 is hormone-regulated and is expressed predominantly in the prostate [3], [4], and to a lesser extent in other tissues [4], [5]. KLK4 has gained support as a potential biomarker for several hormone-dependent cancers [2], and for prostate cancer specifically, in that numerous studies have found KLK4 to be significantly overexpressed in prostate carcinoma tissues compared to benign prostatic hyperplasia [6] and normal tissues [7]–[11]. Of note, KLK4 is known to be expressed as a variety of isoforms [12], with the full length protein (254 amino acids long) showing the potential to be a better biomarker of prostate tumour cells than the commonly expressed shorter isoform (205 amino acids) [11]. In addition, KLK4 has been proposed to play a role in prostate cancer progression through its involvement in epithelial-mesenchymal transition [13], a more aggressive phenotype, and metastases to bone [14]. KLK4 overexpression has been reported to be associated with prostate cancer stage, although the direction of effect differed for KLK4 mRNA (associated with advanced stage) [6] versus KLK4 protein (early stage tumours) [15].

Approximately 40% of prostate cancer is estimated to have a genetic component (http://www.genome.gov/gwastudies/) [16], and to date single nucleotide polymorphisms (SNPs) in over 40 loci have been identified by genome-wide association studies (GWAS) to be associated with prostate cancer risk [17]. One of these SNPs is located in the KLK locus, downstream of the KLK3 gene [18], [19], [20], and is thought to be a marker for a potentially functional non-synonymous SNP within the KLK3 gene [21]. Although no SNPs in KLK4 have been reported by GWAS to be associated with prostate cancer at genome-wide significance levels to date, commonly used GWAS chips only capture 22% [22] - 44% [23] of validated genetic variation in the locus with r²≥0.80. Hence we sought to comprehensively investigate the role of KLK4 in prostate cancer risk and tumour aggressiveness by genotyping the majority of validated genetic variation (±10 kb) around the KLK4 locus in a large prostate cancer study group and male controls not screened for PSA levels.

Materials and Methods

Study Subjects

Study subjects have been described elsewhere [24], [25]. Briefly, from 2004 onwards, 1349 histopathologically-confirmed prostate cancer cases were recruited through private and public urologists in Queensland, Australia via three prostate cancer studies or resources: the Retrospective Queensland Study (N = 154; [26]), the Prostate Cancer Supportive Care and Patient Outcomes Project (ProsCan, N = 857; [25]) and from the Australian Prostate Cancer BioResource (APCB, N = 338; http://www.apccbioresource.org.au/index.html). Men presented to urologists with lower urinary tract symptoms and/or abnormal serum Prostate Specific Antigen (PSA), and 72% of cases possessed prostate tumours of Gleason score 7 or above. Cases ranged in age at diagnosis from 40–88 years (median 63 years). Male controls (N = 1405) with no self-reported personal history of prostate cancer were randomly selected from the Australian Electoral Roll and age-matched (in 5 year groups) and post-code matched to cases (N = 569), or recruited through the Australian Red Cross Blood Services in Brisbane (N = 836). Controls were not screened for PSA levels and analyses excluded 50 controls with age at interview <40 years (the age of the youngest case); included controls ranged in age at interview from 40–89 years of age (median 62 years). All participants had self-reported European ethnicity and gave written informed consent. The study protocol was approved by the Human Research Ethics Committees of the Queensland University of Technology, Queensland Institute of Medical Research, the Mater Hospital (for Brisbane Private Hospital), the Royal Brisbane Hospital, Princess Alexandra Hospital and the Cancer Council Queensland.

SNP Selection and Genotyping

The KLK4 gene region used for SNP selection was chr19∶56091420…56115806 (hg18), which encompasses the longest KLK4 isoform ±10 kb. All SNPs in this region were extracted from the National Center for Biotechnology Information (NCBI) dbSNP build 130 [27], CHIP SNPper [28] and the “ParSNPs” database [29] and duplicates removed. SNPs not classified as validated were removed and validated SNPs were further investigated for occurrence in Europeans using SPSmart [30] and 1000 Genomes [23]. Additional SNPs excluded from investigation included all SNPs on the Illumina 550 K, 610 K and Omni1 genome-wide genotyping chips and SNPs assessed in the Cancer Genetics Markers of Susceptibility (CGEMS) project [31], unless there was evidence of association with prostate cancer by CGEMS (P<0.05). SNPs in high linkage disequilibrium (LD; r²≥0.80) with these excluded Illumina and CGEMS SNPs were also removed, determined by the SNP Annotation and Proxy Search program (SNAP) version 2.1 [32] using HapMap release 22 (1000 Genomes data was not available at the time of initiation of this study). We then prioritised for genotyping all independent SNPs (r²<0.80) according to SNAP using HapMap release 22 data (N = 74). An additional 8 KLK4 tagSNPs (selected using HapMap data release 24/phase II, Nov 2008, NCBI build 36, dbSNP b126, using the Tagger program within Haploview v4.1 [33]), genotyped as part of a previous study, were also included (N = 82 overall).

SNPs were genotyped using iPLEX Gold assays on the Sequenom MassARRAY platform (Sequenom, San Diego, CA), as described previously [34]. There were 4 negative (H₂O) controls per 384-well plate, and quality control parameters included genotype call rates >95%, a combination of cases and controls on each plate, inclusion of 20 duplicate samples per 384-well plate (>5% of samples) with ≥98% concordance between duplicates and Hardy-Weinberg Equilibrium P values >0.05. Of a total of 82 KLK4 SNPs selected for investigation, 11 could not be designed for Sequenom assays, and after application of quality control parameters, 61 SNPs were successfully genotyped. After the study was completed, 1000 Genomes data became available and revealed that 6 KLK4 SNPs not genotyped directly in our study (rs2659108, rs1654556, rs1090648, rs11881373, rs2569531 and rs73598979) were actually tagged by our genotyped SNPs (r²>0.80).

Statistical Methods

Predictive Analytics Software (PASW) Statistics version 17.0.2 (SPSS Inc., Chicago, IL) was used for all analyses. Genotype and allele frequencies were calculated for the patient and control groups. SNP allele and genotype distributions were compared using χ² and their association with prostate cancer susceptibility and clinical data were performed under codominant and linear models using logistic regression analysis. Prostate cancer cases with tumour Gleason scores ≥7 were classified as aggressive. All analyses were adjusted for age (as a continuous variable).

SNP Function Prediction

Alibaba (http://labmom.com/link/alibaba_2_1_tf_binding_prediction), TFsearch (http://www.cbrc.jp/research/db/TFSEARCH.html) and MatInspector (http://www.genomatix.de/online_help/help_matinspector/matinspector_help.html) were used to predict the transcription factor binding sites. The program SignalP was used to predict the signal peptide (http://www.cbs.dtu.dk/services/SignalP/). miRNADA (http://www.microrna.org/microrna/home.do), Patrocles and (http://www.patrocles.org/) miRBase (http://www.mirbase.org/) were used to determine the effect of the SNP alleles on miRNA binding. JASPAR (http://jaspar.binf.ku.dk/), CISTER (http://zlab.bu.edu/~mfrith/cister.shtml) and NHRScan were used for the prediction of nuclear hormone receptor response elements. Splicing effects (using splice-finder), protein structure and stability (Polyphen, SIFT, SNP3D) were determined through the SNPinfo web server (http://snpinfo.niehs.nih.gov). Histone marks, DNAse hypersensitive sites and conservation scores were obtained from HaploReg (http://www.broadinstitute.org/mammals/haploreg/haploreg.php), which extracts data from the UCSC Browser (http://genome.ucsc.edu/). The F-SNP web server (http://compbio.cs.queensu.ca/F-SNP) was used to determine the functional score and putative effect of each SNP.

Results

Seven SNPs were found to be monomorphic in our sample group (Table S1). Results of analyses of the remaining 54 KLK4 SNPs and risk of prostate cancer are displayed in Table 1. Although no KLK4 SNPs were statistically significantly associated with prostate cancer risk after Bonferroni correction (P<9×10⁻⁴), 7 SNPs were associated at the P_trend<0.05 significance level and the majority of these displayed a modest decrease in prostate cancer risk of around 20%. Two of these SNPs, rs268923 (Odds Ratio (OR) 0.89, 95% Confidence Interval (CI) 0.79–1.00, P_trend = 0.045) and rs56112930 (OR 0.37, 95% CI 0.14–0.96, P_trend = 0.040; Minor Allele Frequency (MAF) 0.006), are located several kilobases upstream of the long isoform of KLK4, 8.2 kb and 6.5 kb, respectively. The KLK4 tagSNP rs7248321 (OR 0.77, 95% CI 0.60–0.98, P_trend = 0.033; MAF 0.060), also represented on several of the genome-wide chips including the Illumina 550 K, 610 K and Omni1 chips, is located ∼4.5 kb upstream of KLK4, and rs7248321 tags two nearby SNPs rs13345980 and rs7246794 with r² 1.00. rs1654551 (OR 0.79, 95% CI 0.65–0.97, P_trend = 0.023; MAF 0.093) is a non-synonymous SNP in the full length KLK4 protein coding for a serine to alanine amino acid (aa) substitution at position 22. In the more commonly expressed 205 aa KLK4 isoform, this SNP is located in the 5′ untranslated region [11]. The remaining three risk-associated SNPs, rs1701927, rs1090649 and rs806019, we determined to be in high LD with each other (r²≥0.98) and accordingly all display ORs of around 0.85 (95% CI range 0.73–1.00, P_trend range 0.030–0.044; MAF 0.175). Results were similar for rs1701926 that is also part of this high LD block (OR 0.86, 95% CI 0.74–1.00, P_trend = 0.058). All four SNPs are located downstream of the KLK4 gene from 750 base pairs (bp) to 3.6 kb past the 3′ untranslated region (UTR).

[Figure omitted. See PDF.]

Table 1. Association of KLK4 SNPs and prostate cancer risk.

https://doi.org/10.1371/journal.pone.0044520.t001

Only one KLK4 SNP, rs198968, was associated with prostate tumour aggressiveness (Gleason score <7 vs. ≥7: OR 0.76 (95% CI 0.60–0.95, P_trend = 0.016; Table 2). However, this result was not reflected in a more robust Gleason score analysis comparing “extreme” Gleason categories, ≤6 (N = 329) vs. ≥8 (N = 173), with an odds ratio of 0.95 (95% CI 0.68–1.32; P_trend = 0.752).

[Figure omitted. See PDF.]

Table 2. Association of KLK4 SNPs and prostate tumour aggressiveness.

https://doi.org/10.1371/journal.pone.0044520.t002

Results of bioinformatic prediction of functions of the associated SNPs are provided in Table S2. SNPs rs268923, rs198968, rs1654551, rs1701926, rs1090649 and rs806019 were found to alter transcription factor binding sites as predicted by at least one prediction tool. rs198968 and rs1654551 also lie within promoter histone marks as well as DNAse hypersensitive sites (Table S2), and hence are better candidate for functional follow up studies.

SNP rs1654551 leads to a serine to alanine amino acid change, but is predicted to be benign using the FASTSNP web server, although the SNP is predicted to effect o-glycosylation. In addition, PsortII prediction (http://urgi.versailles.inra.fr/Tools/PsortII) predicted the serine variant to be only 44.4% extracellularly localised as compared to the alanine variant which is predicted to be 55.6% extracellular. This is backed by SignalP predicting alteration of the KLK4 signal peptide sequence for the serine variant. Further, this SNP is also predicted to be involved in differential splicing.

Discussion

We performed a comprehensive investigation of the role of variation in the KLK4 gene in prostate cancer risk and/or tumour aggressiveness by assessing the majority of SNPs that have not been covered by previously performed GWA studies. Our study of approximately 1,300 cases and 1,300 male controls provided suggestive evidence that several KLK4 SNPs may be associated with decreased risk of prostate cancer, and bioinformatic analysis provides evidence that some of these have potential biological relevance in prostate cancer.

None of the nominally risk-associated SNPs were located in known KLK4 hormone response elements [35]. Three SNPs lay several kb upstream of the KLK4 gene. rs7248321 is a tagSNP that has not been previously reported to be associated with prostate cancer risk in any GWAS, including CGEMS, and given the large numbers of samples assessed in previous studies [17], it is likely to be a false-positive result. Bioinformatic analyses of the rare rs56112930 SNP did not reveal any predicted effects on transcription factor binding sites [36], [37], [38]. SNP rs268923 was calculated by three different transcription factor binding site prediction programs to possibly have an effect [36]–[39], and although each program predicts different transcription factor binding sites to be altered by the SNP; one example of a prostate cancer relevant result is the predicted gain of an Oct-1 site [36]. Oct-1 is a known co-regulator of the androgen receptor [40], regulates growth of prostate cancer cells and is associated with poor prognosis [41].

The only SNP located in the KLK4 coding region found to be marginally associated with prostate cancer risk was rs1654551. Since splicing of the KLK4 locus is complex and results in several KLK4 mRNA forms being produced [12], there are several possible functional consequences of this substitution. The two protein isoforms identified to date, in order of expression in normal prostate, are an intracellular 205 amino acid (aa) protein which lacks the classical KLK signal peptide and is localised to the nucleus (“short” isoform) [9], [11], and a secreted 254 aa protein that is cytoplasmically localised [11], [13]. rs1654551 codes for a serine to alanine substitution at amino acid 22 of the long isoform, or is located in the 5′ UTR of the 205 aa KLK4 protein. Although both the short and long isoforms have been found to be overexpressed in prostate cancer cells, the “long” 254 aa KLK4 protein is better able to discriminate between tumour and normal cells [11] and hence may be the more biologically relevant isoform in prostate cancer. Amino acid 22 is located within the signal peptide region of KLK4, which is cleaved off between aa 26 and 27 to result in secretion. It is unknown what the potential functional effects of an amino acid substitution are within the signal peptide. However, a recent study has shown that this cleaved peptide may be a useful target in prostate cancer immunotherapy, with the KLK4 signal peptide successfully inducing and expanding the cytotoxic T lymphocyte response more readily than PSA or Prostatic Acid Phosphotase (PAP) [42]. In addition, in silico analysis using the signal peptide prediction program SignalP [43] predicted a Serine22Alanine substitution to alter the cleavage site from aa 26/27 to aa 21/22. This would result in a KLK4 pro-protein with an additional 5 aa, which could potentially affect localisation or possibly even activation of the KLK4 proenzyme. Of relevance, a form of PSA has been reported that has an altered signal/pro-peptide and, although the pro-PSA sequence is truncated (not lengthened as is predicted for KLK4), the signal peptide alteration does result in an isoform of PSA that is unable to be activated [44]. This [-2]pro-PSA isoform is now also the basis of a commercially available prostate cancer serum test [45].

Attempting to predict the possible functional effects of the four associated SNPs located downstream of KLK4, rs1701927, rs1701926, rs1090649 and rs806019, is not as clearly directed. It is possible that some or all of these SNPs might alter enhancer/silencer binding sites, affecting expression of KLK4. In silico transcription factor binding site analysis predicts that rs1701926, rs1090649 and rs806019 may alter transcription factor binding sites [36]–[39] relevant to prostate cancer. The closest validated gene to the KLK4 3′ end is the Kallikrein family pseudogene KLKP1, thus these SNPs may regulate the activity of this pseudogene or expression of KLKP1 transcripts, which have been shown to be down-regulated in prostate cancer tissues [46]. In addition, one other SNP not genotyped in this study, rs1654556, is in high LD (r²≥0.80) with these SNPs [23] and is predicted to alter mRNA folding [47] and miRNA binding (Table S2).

To the best of our knowledge, only two other studies have examined the role of KLK4 SNPs in cancer (aside from genome-wide investigations). Recently Klein et al. investigated the effect of common variation in the exons and putative promoter regions of all 15 KLK genes on prostate cancer risk and levels of PSA forms and KLK2 [48]. Five KLK4 SNPs – rs198969, rs198968, rs1654552, rs1654551 and rs1654553 - were assessed for association with prostate cancer risk in the Cancer Prostate in Sweden (CAPS) 1 sample set of over 1,400 cases and 700 controls and none were found to be associated. A second study in a small Korean sample set of 117 breast cancer cases and 194 controls found KLK4 SNP rs806019 to be associated with a decreased risk of breast cancer (Odds Ratio 0.53; 95% Confidence Interval 0.33–0.85; P = 0.007) [49], a finding of similar magnitude and direction to that observed in our study of prostate cancer. Part of our study design was to exclude KLK4 SNPs already assessed in GWAS, except for those reported to be associated with prostate cancer at the P<0.05 level in CGEMS. Four SNPs in CGEMS gave evidence of association with prostate cancer - rs17714461, rs8101572, rs8100631 and rs10427094 [31]. All four were genotyped in this study, but none were associated in our sample set. Of note, rs17714461 was recently reported to interact with the GWAS-detected KLK2/3 SNP rs2735839 in CGEMS [50]. The authors mention that this result is notable considering KLK4 and KLK2 collaborate to stimulate cellular proliferation in prostate cancer [51]. rs2735839 genotype data was available for only a small proportion of our samples and hence we did not investigate this interaction further.

Our well-sized study indicates a possible contribution of SNPs in the KLK4 gene to decreased risk of prostate cancer. However, these results should be interpreted cautiously considering the number of tests performed, and validation in much larger sample sets such as those of the PRACTICAL consortium is necessary.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

rs IDs found to be monomorphic in this study.

https://doi.org/10.1371/journal.pone.0044520.s001

(DOC)

Table S2.

In silico function prediction of prostate cancer associated KLK4 gene variants.

https://doi.org/10.1371/journal.pone.0044520.s002

(XLSX)

Acknowledgments

The authors thank the many patients and control subjects who participated so willingly in this study, and the numerous institutions and their staff who have supported recruitment. The authors are very grateful to staff at the Australian Red Cross Blood Services for their assistance with the collection of risk factor information and blood samples of healthy donor controls; members of the Cancer Council Queensland for ProsCan participant information, including Megan Ferguson and Andrea Kittila; the hospitals that participated in recruitment for the ProsCan study: Greenslopes Private, Royal Brisbane, Mater Adults, Princess Alexandra, Ipswich, QEII, Redlands and Redcliffe Hospitals, Townsville General Hospital, and Mackay Base Hospital; Drs. John Yaxley and David Nicol for recruitment of patients into the Retrospective Queensland Study; and the Urological Society of Australia and New Zealand. Thank you to members of the Australian Prostate Cancer Research Centre-Queensland at QUT and the QIMR Molecular Cancer Epidemiology Laboratory for their assistance with recruitment and biospecimen processing, and Dr John Lai, XiaoQing Chen and Dr Jonathan Beesley for technical advice.

Contributors

Membership of the Australian Prostate Cancer BioResource denotes the Queensland node participants in this study and includes: Trina Yeadon, Pamela Saunders, Allison Eckert and Judith Clements (Institute of Health and Biomedical Innovation, Queensland University of Technology, Brisbane, Queensland, Australia); Peter Heathcote, Glenn Wood, Greg Malone (Brisbane Urology Clinic, Brisbane Urology Clinic Central Queensland Urology Clinic, Wickham Terrace, Brisbane, Queensland, Australia); Hema Samaratunga (Aquesta Pathology, Toowong, Queensland, Australia); Angus Collins, Megan Turner and Kris Kerr (Sullivan and Nicolaides Pathology, Brisbane, Queensland, Australia).

Author Contributions

Conceived and designed the experiments: JB JAC ABS. Performed the experiments: FL JB SS TO. Analyzed the data: LM. Contributed reagents/materials/analysis tools: FL SC RAG SC JFA Australian Prostate Cancer BioResource ABS JB JAC. Wrote the paper: FL JB.

Citation: Lose F, Srinivasan S, O’Mara T, Marquart L, Chambers S, Gardiner RA, et al. (2012) Genetic Association of the KLK4 Locus with Risk of Prostate Cancer. PLoS ONE 7(9): e44520. https://doi.org/10.1371/journal.pone.0044520

References

1. Lawrence MG, Lai J, Clements JA (2010) Kallikreins on steroids: structure, function, and hormonal regulation of prostate-specific antigen and the extended kallikrein locus. Endocr Rev 31: 407–446.

2. Mavridis K, Scorilas A (2010) Prognostic value and biological role of the kallikrein-related peptidases in human malignancies. Future Oncol 6: 269–285.

3. Nelson PS, Gan L, Ferguson C, Moss P, Gelinas R, et al. (1999) Molecular cloning and characterization of prostase, an androgen-regulated serine protease with prostate-restricted expression. Proc Natl Acad Sci U S A 96: 3114–3119.

4. Yousef GM, Obiezu CV, Luo LY, Black MH, Diamandis EP (1999) Prostase/KLK-L1 is a new member of the human kallikrein gene family, is expressed in prostate and breast tissues, and is hormonally regulated. Cancer Res 59: 4252–4256.

5. Obiezu CV, Shan SJ, Soosaipillai A, Luo LY, Grass L, et al. (2005) Human kallikrein 4: quantitative study in tissues and evidence for its secretion into biological fluids. Clin Chem 51: 1432–1442.

6. Avgeris M, Stravodimos K, Scorilas A (2011) Kallikrein-related peptidase 4 gene (KLK4) in prostate tumors: quantitative expression analysis and evaluation of its clinical significance. Prostate 71: 1780–1789.

7. Day CH, Fanger GR, Retter MW, Hylander BL, Penetrante RB, et al. (2002) Characterization of KLK4 expression and detection of KLK4-specific antibody in prostate cancer patient sera. Oncogene 21: 7114–7120.

8. Obiezu CV, Soosaipillai A, Jung K, Stephan C, Scorilas A, et al. (2002) Detection of human kallikrein 4 in healthy and cancerous prostatic tissues by immunofluorometry and immunohistochemistry. Clin Chem 48: 1232–1240.

9. Xi Z, Klokk TI, Korkmaz K, Kurys P, Elbi C, et al. (2004) Kallikrein 4 is a predominantly nuclear protein and is overexpressed in prostate cancer. Cancer Res 64: 2365–2370.

10. Klokk TI, Kilander A, Xi Z, Waehre H, Risberg B, et al. (2007) Kallikrein 4 is a proliferative factor that is overexpressed in prostate cancer. Cancer Res 67: 5221–5230.

11. Dong Y, Bui LT, Odorico DM, Tan OL, Myers SA, et al. (2005) Compartmentalized expression of kallikrein 4 (KLK4/hK4) isoforms in prostate cancer: nuclear, cytoplasmic and secreted forms. Endocr Relat Cancer 12: 875–889.

12. Kurlender L, Borgono C, Michael IP, Obiezu C, Elliott MB, et al. (2005) A survey of alternative transcripts of human tissue kallikrein genes. Biochim Biophys Acta 1755: 1–14.

13. Veveris-Lowe TL, Lawrence MG, Collard RL, Bui L, Herington AC, et al. (2005) Kallikrein 4 (hK4) and prostate-specific antigen (PSA) are associated with the loss of E-cadherin and an epithelial-mesenchymal transition (EMT)-like effect in prostate cancer cells. Endocr Relat Cancer 12: 631–643.

14. Gao J, Collard RL, Bui L, Herington AC, Nicol DL, et al. (2007) Kallikrein 4 is a potential mediator of cellular interactions between cancer cells and osteoblasts in metastatic prostate cancer. Prostate 67: 348–360.

15. Seiz L, Kotzsch M, Grebenchtchikov NI, Geurts-Moespot AJ, Fuessel S, et al. (2010) Polyclonal antibodies against kallikrein-related peptidase 4 (KLK4): immunohistochemical assessment of KLK4 expression in healthy tissues and prostate cancer. Biol Chem 391: 391–401.

16. Gronberg H (2003) Prostate cancer epidemiology. Lancet 361: 859–864.

17. Hindorff LA, Junkins HA, Hall PN, Mehta JP, Manolio TA (2010) A Catalog of Published Genome-Wide Association Studies.

18. Eeles RA, Kote-Jarai Z, Giles GG, Olama AA, Guy M, et al. (2008) Multiple newly identified loci associated with prostate cancer susceptibility. Nat Genet 40: 316–321.

19. Kote-Jarai Z, Easton DF, Stanford JL, Ostrander EA, Schleutker J, et al. (2008) Multiple novel prostate cancer predisposition loci confirmed by an international study: the PRACTICAL Consortium. Cancer Epidemiol Biomarkers Prev 17: 2052–2061.

20. Eeles RA, Giles GG, Neal DE, Muir K, Easton DF (2008) Reply to "Variation in KLK genes, prostate-specific antigen and risk of prostate cancer". Nat Genet 40: 1035–1036.

21. Kote-Jarai Z, Amin Al Olama A, Leongamornlert D, Tymrakiewicz M, Saunders E, et al. (2011) Identification of a novel prostate cancer susceptibility variant in the KLK3 gene transcript. Hum Genet 129: 687–694.

22. The International HapMap Consortium (2003) The International HapMap Project. Nature 426: 789–796.

23. 1000 Genome Project Consortium (2010) A map of human genome variation from population-scale sequencing. Nature 467: 1061–1073.

24. Lose F, Batra J, O'Mara T, Fahey P, Marquart L, et al.. (2011) Common variation in Kallikrein genes KLK5, KLK6, KLK12, and KLK13 and risk of prostate cancer and tumor aggressiveness. Urol Oncol.

25. Baade PD, Aitken JF, Ferguson M, Gardiner RA, Chambers SK (2010) Diagnostic and treatment pathways for men with prostate cancer in Queensland: investigating spatial and demographic inequalities. BMC Cancer 10: 452.

26. Lai J, Kedda MA, Hinze K, Smith RL, Yaxley J, et al. (2007) PSA/KLK3 AREI promoter polymorphism alters androgen receptor binding and is associated with prostate cancer susceptibility. Carcinogenesis 28: 1032–1039.

27. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, et al. (2001) dbSNP: the NCBI database of genetic variation. Nucleic Acids Res 29: 308–311.

28. Riva A, Kohane IS (2004) A SNP-centric database for the investigation of the human genome. BMC Bioinformatics 5: 33.

29. Goard CA, Bromberg IL, Elliott MB, Diamandis EP (2007) A consolidated catalogue and graphical annotation of dbSNP polymorphisms in the human tissue kallikrein (KLK) locus. Mol Oncol 1: 303–312.

30. Amigo J, Salas A, Phillips C, Carracedo A (2008) SPSmart: adapting population based SNP genotype databases for fast and comprehensive web access. BMC Bioinformatics 9: 428.

31. Yeager M, Orr N, Hayes RB, Jacobs KB, Kraft P, et al. (2007) Genome-wide association study of prostate cancer identifies a second risk locus at 8q24. Nat Genet 39: 645–649.

32. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, et al. (2008) SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics 24: 2938–2939.

33. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.

34. Lose F, Nagle CM, O'Mara T, Batra J, Bolton KL, et al. (2010) Vascular endothelial growth factor gene polymorphisms and ovarian cancer survival. Gynecol Oncol 119: 479–483.

35. Lai J, Myers SA, Lawrence MG, Odorico DM, Clements JA (2009) Direct progesterone receptor and indirect androgen receptor interactions with the kallikrein-related peptidase 4 gene promoter in breast and prostate cancer. Mol Cancer Res 7: 129–141.

36. Grabe N (2000) AliBaba2.1.

37. Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, et al. (1998) Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res 26: 362–367.

38. Akiyama Y (2002) TFSEARCH: Searching Transcription Factor Binding Sites.

39. Cartharius K, Frech K, Grote K, Klocke B, Haltmeier M, et al. (2005) MatInspector and beyond: promoter analysis based on transcription factor binding sites. Bioinformatics 21: 2933–2942.

40. Gonzalez MI, Robins DM (2001) Oct-1 preferentially interacts with androgen receptor in a DNA-dependent manner that facilitates recruitment of SRC-1. J Biol Chem 276: 6420–6428.

41. Obinata D, Takayama K, Urano T, Murata T, Kumagai J, et al. (2012) Oct1 regulates cell growth of LNCaP cells and is a prognostic factor for prostate cancer. Int J Cancer 130: 1021–1028.

42. Wilkinson R, Woods K, D'Rozario R, Prue R, Vari F, et al. (2012) Human kallikrein 4 signal peptide induces cytotoxic T cell responses in healthy donors and prostate cancer patients. Cancer Immunol Immunother 61: 169–179.

43. Bendtsen JD, Nielsen H, von Heijne G, Brunak S (2004) Improved prediction of signal peptides: SignalP 3.0. J Mol Biol 340: 783–795.

44. Mikolajczyk SD, Marker KM, Millar LS, Kumar A, Saedi MS, et al. (2001) A truncated precursor form of prostate-specific antigen is a more specific serum marker of prostate cancer. Cancer Res 61: 6958–6963.

45. Catalona WJ, Partin AW, Sanda MG, Wei JT, Klee GG, et al. (2011) A multicenter study of [-2]pro-prostate specific antigen combined with prostate specific antigen and free prostate specific antigen for prostate cancer detection in the 2.0 to 10.0 ng/ml prostate specific antigen range. J Urol 185: 1650–1655.

46. Kaushal A, Myers SA, Dong Y, Lai J, Tan OL, et al. (2008) A novel transcript from the KLKP1 gene is androgen regulated, down-regulated during prostate cancer progression and encodes the first non-serine protease identified from the human kallikrein gene locus. Prostate 68: 381–399.

47. Zuker M (2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 31: 3406–3415.

48. Klein RJ, Hallden C, Cronin AM, Ploner A, Wiklund F, et al. (2010) Blood biomarker levels to aid discovery of cancer-related single-nucleotide polymorphisms: kallikreins and prostate cancer. Cancer Prev Res (Phila) 3: 611–619.

49. Lee JY, Park AK, Lee KM, Park SK, Han S, et al. (2009) Candidate gene approach evaluates association between innate immunity genes and breast cancer risk in Korean women. Carcinogenesis 30: 1528–1531.

50. Ciampa J, Yeager M, Amundadottir L, Jacobs K, Kraft P, et al. (2011) Large-scale exploration of gene-gene interactions in prostate cancer using a multistage genome-wide association study. Cancer Res 71: 3287–3295.

51. Mize GJ, Wang W, Takayama TK (2008) Prostate-specific kallikreins-2 and -4 enhance the proliferation of DU-145 prostate cancer cells through protease-activated receptors-1 and -2. Mol Cancer Res 6: 1043–1051.