Association between three genetic variants in kallikrein 3 and prostate cancer risk.

BACKGROUND: Epidemiological studies have assessed the association between kallikrein 3 (KLK3) polymorphisms and prostate cancer (PCa) susceptibility. However, published data on this association are somewhat inconclusive.
METHODS: Articles investigating the association between three KLK3 (rs1058205, rs2735839, and rs266882) variants and PCa susceptibility were searched from online databases, which included 35,838 patients and 36,369 control participants. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to demonstrate the strength of the association. We also utilized ELISA to detect serum expression of KLK3. In addition, in silico tools were adopted to evaluate the relationship of KLK3 expression and PCa survival time.
RESULTS: The overall results indicated that polymorphism T>C of rs1058205 was associated with decreased risk of PCa (allele contrast: OR = 0.75, 95% CI = 0.64-0.88, P heterogeneity < 0.001; homozygote comparison: OR = 0.58, 95% CI = 0.42-0.81, P heterogeneity < 0.001), particularly in Caucasian population (allele contrast: OR = 0.77, 95% CI = 0.65-0.91, P heterogeneity < 0.001; homozygote comparison: OR = 0.58, 95% CI = 0.41-0.82, P heterogeneity < 0.001). No association was observed between the polymorphism A>G of rs2735839 and risk of PCa. In addition, no association was observed between polymorphism A>G of rs266882 and risk of PCa. Serum KLK3 levels in PCa patients carrying CC/CT genotypes were statistically lower than those carrying TT genotypes.
Conclusion: This meta-analysis suggests that rs1058205 polymorphism of KLK3 is a risk factor for PCa development, polymorphism T>C of rs1058205 is associated with decreased susceptibility to PCa particularly in Caucasian population.

Introduction

Prostate cancer (PCa) is one of the most common cancers among men in U.S.A. [1]. However, the exact etiology underlying the development and progression of PCa is still poorly understood. Previous studies indicated that both environmental and genetic factors may be associated with the development and progression of PCa [2,3]. Single nucleotide polymorphism (SNP) is considered as the most common type of genetic variation and has been reported to be related to increased susceptibility to PCa [4,5]. In recent decades, genome-wide association studies (GWAS) in multiple ethnicities have identified some susceptibility genes in PCa. However, the molecular mechanisms of these SNPs remain poorly defined [6].

Serum levels of prostate-specific antigen (PSA) are measured to screen for PCa. Approximately 40–45% of variations in PSA serum levels are estimated to be influenced by genetic components [7,8]. A previous GWAS on European population identified six genetic loci that were significantly associated with PSA level [9]. Although genetic analyses of PSA levels from Asian countries are relatively limited compared with the European countries, GWAS from China has demonstrated that three loci, solute carrier family 45 member 3 (SLC45A3), microseminoprotein-β (MSMB) and kallikrein 3 (KLK3), were significantly associated with PSA levels [10]. The KLK locus is clustered in a tandem array of approximately 300 kb on chromosome 19q13.4 and contains the largest cluster of 15 homologous protease genes. In addition, the KLK3 gene is a serine protease kallikrein family member and may be involved in the occurrence and metastasis of PCa [11,12].

A number of epidemiological studies have explored the relationship between KLK3 polymorphisms and PCa risk. Previous studies have shown that the TT genotype of rs1058205 was associated with decreased PCa aggressiveness in the Caucasian population based on the Gleason score [13]. However, Zhang and colleagues found no correlation between any of the SNPs at the 19q13.33 locus and PCa risk in Han Chinese men [14]. Various case–control studies showed inconsistent results on the association between KLK3 polymorphism and PCa. Hence, in our meta-analysis, accumulated data from all eligible studies were utilized to summarize and enhance the statistical powers of these published studies for the three KLK3 variants. All the articles that met the selection criteria were included in the present analysis. The pooled results from all these studies allowed us to draw novel conclusions.

Materials and methods

Search strategy and identification of eligible studies

To identify all potentially eligible studies on KLK3 polymorphisms and PCa risk, a systematic search was carried out on PubMed, Embase and SinoMed (from China) databases, covering all researches published literature worldwide due to July 30, 2017. The search terms “kallikrein 3 or KLK3,” “polymorphism or variant,” and “prostate cancer or tumor or adenocarcinoma” were used to identify suitable literature. Furthermore, review articles and bibliographies of other studies were identified through a manual search. All full-text case–control studies were included if they had PCa risk and KLK3 polymorphism, sufficient available data to estimate an odds ratio (OR) with 95% confidence interval (CI) and the frequencies of the allele or genotype available.

Data extraction and quality assessment

Two researchers independently extracted and reviewed the information complying with the selection criteria. In the case of a disagreement, the researchers discussed the studies until a consensus was reached. For each study, the following characteristics were summarized: first author’s last name, year of publication, sample size of cases and controls, country of origin, ethnicity, study-design (sources of samples), age range, the number of cases and controls with variant allele and wild type, and genotyping methods. Finally, nine articles were included [15-23] and the Newcastle–Ottawa scale (NOS) score was shown in Table 1.

Table 1

Basic information for included studies of the association between three KLK3 polymorphisms and prostate cancer susceptibility

Author [Ref.]	Year	Origin	Ethnicity	Source	Case	Control	Case			Control			HWE	Age range		Method	NOS
														Case	Control
rs2735839							GG	GA	AA	GG	GA	AA
Hu [15]	2014	China	Asian	PB	108	242	25	68	15	44	108	90	0.25	70.01 ± 7.26	70.52 ± 7.11	TaqMan	8
Wang [19]	2013	China	Asian	PB	285	280	102	126	57	98	144	38	0.19	NA	NA	HRM-PCR	7
Eeles [18]	2008	U.K.	Caucasian	PB	1851	1887	35	406	1410	103	591	1193	0.01	NA	NA	Illumina array	7
Parikh [22]	2011	U.S.A.	Caucasian	PB	3150	2980	260	299	2591	261	288	2431	<0.01	NA	NA	GWAS	8
rs266882							GG	GA	AA	GG	GA	AA
Cicek [16]	2005	U.S.A.	Caucasian	FB	439	479	119	196	124	119	236	124	0.75	NA	NA	PCR-RFLP	6
Lai [17]	2007	Australia	Caucasian	PB	209	223	35	111	63	71	107	45	0.68	68 (50–91)	50 (18–75)	PCR-RFLP	8
Penney [20]	2011	U.S.A.	Caucasian	PB	966	1279	255	483	228	346	622	311	0.34	59 ± 8.3	59.3 ± 8.0	Sequenom technology	8
rs1058205							CC	CT	TT	CC	CT	TT
Penney [20]	2011	U.S.A.	Caucasian	PB	896	1165	14	176	706	28	269	868	0.19	59 ± 8.3	59.3 ± 8.0	Sequenom technology	8
Chen [21]	2017	China	Asian	HB	268	298	4	55	209	7	89	202	0.44	58 (52–62)	62 (52–67)	HRM-PCR	7
Eeles [18]	2008	U.K.	Caucasian	PB	1850	1887	54	490	1306	136	674	1077	0.03	NA	NA	Illumina array	7
Stegeman [23]	2015	Australia	Caucasian	PB	22301	22320	446	5798	161057	670	6476	15177	0.53	64.8 ± 8.0	60.6 ± 10.7	Illumina array	8
Parikh [22]	2011	U.S.A.	Caucasian	PB	3515	3329	373	329	2813	410	298	2621	<0.01	NA	NA	GWAS	8

Abbreviations: FB, family-based; GWAS, Genome-Wide Association Study; HB, hospital-based; HRM-PCR, high-resolution melting curve polymerase chain reaction; HWE, Hardy–Weinberg equilibrium of controls; NA, not available; NOS, Newcastle–Ottawa scale; PB, population-based; PCR-FLIP, polymerase chain reaction and restrictive fragment length polymorphism.

Enzyme-linked immunosorbent assay (ELISA)

We collected the blood of participants in standard cubes without anticoagulant according to the instructions. Patients were enrolled from the Affiliated Changzhou No.2 People’s Hospital of Nanjing Medical University and third Affiliated Hospital of Nantong University. Then, we centrifuged the blood sample immediately at 1000 rpm for 20 min by ELISA kit (Thermo Scientific Ltd.). A standard ELISA plate reader was used for measuring absorbance at 450 and 550 nm. We correlated the absorbance against a standard curve. Written informed consent by each participant was obtained before samples were collected. The present study was approved by Institution Review Board of these two hospitals.

Statistical analysis

The strength of the association between the three SNPs and PCa risk was evaluated by OR with its corresponding 95% CI. The statistical significance of the OR was evaluated with a Z-test. The heterogeneity assumption among studies was estimated utilizing a chi-square-based Q-test. The wild-type homozygote (WW), the risk of rare allele homozygote (RR), and heterozygous (WR) genotypes were evaluated. Then we assessed the risk of PCa under a dominant model (RR+WR vs. WW) and recessive model (RR vs. WR+WW). Subgroup analysis was stratified by ethnicity and source of control (hospital-based, HB, and population-based, PB). A P-value of greater than 0.10 for the Q-test indicated a lack of heterogeneity among the studies. If a significant heterogeneity was detected, the DerSimonian–Laird (random-effects model) was applied; otherwise the Mantel–Haenszel (fixed-effects model) was performed to calculate the pooled OR [24,25]. Hardy–Weinberg equilibrium (HWE) was calculated using the Pearson’s chi-square test. Publication bias was assessed using both Egger’s test and Begg’s test, and a P-value of less than 0.05 was considered significant [26]. All statistical tests for our meta-analysis were carried out using the Stata software version 11.0 (StataCorp LP, College Station, TX, U.S.A.).

Quality score assessment, sensitivity analysis, and in silico analysis

The NOS was used to evaluate the quality of each study [27]. The NOS is a tool utilized for evaluating the quality of studies included in meta-analysis. Stars are awarded such that the highest quality studies are awarded up to nine stars. In the present analysis, studies with more than seven stars should be considered as high quality. We conducted leave-one-out sensitivity analysis to assess the stability of the overall results. The sensitivity was evaluated by removing each research once at a time in every genetic model for the three KLK3 variants. Furthermore, we used the online gene expression mini database to demonstrate the KLK3 expression in prostate tissues (http://gemini.cancer-pku.cn/).

Results

Study characteristics

A total of nine articles (including 12 case–control studies) met the above selection criteria and are included in our meta-analysis (Figure 1). The characteristic of these studies were summarized in Table 1. Furthermore, we checked the minor allele frequency (MAF) of three KLK3 variants reported for the main worldwide populations from 1000 genomes online. The MAFsof rs1058205 were as follows: American (AMR), 0.308; African (AFR), 0.477; East Asian (EAS), 0.442; South Asian (SAS), 0.395; and European (EUR), 0.171. Then, those of rs2735839 were as follows: AMR, 0.264; AFR, 0.391; EAS, 0.442; SAS, 0.347; and EUR, 0.140. At last, those of rs266882 were as follows: AMR, 0.290; AFR, 0.567; EAS, 0.144; SAS, 0.377; and EUR, 0.504 (Figure 2). Overall, 6872 PCa patients and 7740 control participants with three KLK3 polymorphisms were evaluated. In the subgroup of ethnicity, six studies were performed in individuals of European descent and three studies were in Asian descendants. Population-based controls were carried out in seven of the studies. Polymerase chain reaction-restriction fragment length polymorphism (RFLP), the classical genotyping method, was utilized in two of the comparisons. The distribution of genotypes in seven controls was consistent with HWE.

Figure 1

Flowchart illustrating the search strategy used to identify several studies investigating the association between three KLK3 variants and prostate cancer risk

Figure 2

Minor allele and major allele frequencies for KLK3 rs1058205 (A), rs2735839 (B), and rs266882 (C) polymorphisms in the controls stratified by ethnicity.

Vertical axis, allele frequency; horizontal axis, allele type. AMR: American; AFR: African; EAS: East Asian; SAS: South Asian; EUR: European.

Quantitative synthesis

The overall results indicated a significant association between the rs1058205 T>C polymorphism and PCa (allele contrast: OR = 0.75, 95% CI = 0.64–0.88, Pheterogeneity < 0.001; heterozygote comparison: OR = 0.787, 95% CI = 0.65–0.93, Pheterogeneity < 0.001; homozygote comparison: OR = 0.58, 95% CI = 0.42–0.81, Pheterogeneity < 0.001, Figure 3; dominant genetic model: OR = 0.74, 95% CI = 0.62–0.88, Pheterogeneity < 0.001; and recessive genetic model: OR = 0.63, 95% CI = 0.48–0.82, Pheterogeneity < 0.001) (Table 2), especially in the Caucasian population (allele contrast: OR = 0.77, 95% CI = 0.65–0.91, Pheterogeneity < 0.001; heterozygote comparison: OR = 0.80, 95% CI = 0.66–0.98, Pheterogeneity < 0.001; homozygote comparison: OR = 0.58, 95% CI = 0.41–0.82, Pheterogeneity < 0.001; dominant genetic model: OR = 0.76, 95% CI = 0.63–0.92, Pheterogeneity < 0.001; and recessive genetic model: OR = 0.63, 95% CI = 0.47–0.83, Pheterogeneity < 0.001. However, no association was observed between the rs2735839 A>G polymorphism and PCa risk in allele contrast (OR = 0.93, 95% CI = 0.63–1.39, Pheterogeneity < 0.001), heterozygote comparison (OR = 0.98, 95% CI = 0.59–1.64, Pheterogeneity < 0.001); homozygote comparison (OR = 0.85, 95% CI = 0.40–1.82, Pheterogeneity < 0.001, Figure 4); dominant genetic model (OR = 0.97, 95% CI = 0.58–1.61, Pheterogeneity < 0.001); and recessive genetic model (OR = 0.81, 95% CI = 0.48–1.34, Pheterogeneity < 0.001). In addition, no association was observed between rs266882 A>G variant and PCa risk in allele contrast (OR = 0.87, 95% CI = 0.67–1.12, Pheterogeneity = 0.003); heterozygote comparison (OR = 0.95, 95% CI = 0.81–1.12, Pheterogeneity = 0.240); homozygote comparison (OR = 0.75, 95% CI = 0.45–1.26, Pheterogeneity = 0.002, Figure 5); dominant genetic model (OR = 0.86, 95% CI = 0.65–1.15, Pheterogeneity = 0.060); and recessive genetic model (OR = 0.66, 95% CI = 0.30–1.46, Pheterogeneity = 0.001).

Figure 3

Forest plot of CC vs. TT genetic model of KLK3 rs1058205 polymorphism

Table 2

Stratified analyses of three KLK3 polymorphisms on prostate cancer risk

			R-allele vs. W-allele			WR vs. WW			RR vs. WW			RR+WR vs. WW			RR vs. WR+WW
Variables	N	Case/Control	OR (95%CI)	Pheter*	P	OR (95%CI)	Pheter*	P	OR (95%CI)	Pheter*	P	OR (95%CI)	Pheter*	P	OR (95%CI)	Pheter*	P
rs2735839	4	5394/5389	0.93 (0.63–1.39)	<0.001	0.736	0.98 (0.59–1.64)	<0.001	0.947	0.85 (0.40–1.82)	<0.001	0.678	0.97 (0.58–1.61)	<0.001	0.906	0.81 (0.48–1.34)	<0.001	0.407
Ethnicity
Asian	2	393/522	1.24 (0.63–2.44)	0.001	0.526	1.47 (0.23–9.20)	<0.001	0.682	1.50 (0.32–7.14)	<0.001	0.609	1.50 (0.26–8.56)	<0.001	0.647	1.11 (0.83–1.49)	0.418	0.469
Caucasian	2	5001/4867	0.72 (0.42–1.24)	<0.001	0.236	0.75 (0.45–1.25)	<0.001	0.267	0.53 (0.17–1.67)	<0.001	0.276	0.72 (0.41–1.26)	<0.001	0.248	0.57 (0.21–1.56)	<0.001	0.274
rs266882	3	1614/1981	0.87 (0.67–1.12)	0.003	0.276	0.95 (0.81–1.12)	0.240	0.533	0.75 (0.45–1.26)	0.002	0.274	0.86 (0.65–1.15)	0.062	0.305	0.81 (0.53–1.25)	0.002	0.347
Source of control
PB	2	1175/1502	0.79 (0.48–1.29)	0.001	0.349	1.00 (0.83–1.21)	0.170	0.981	0.61 (0.22–1.71)	0.001	0.352	0.81 (0.46–1.41)	0.060	0.456	0.66 (0.30–1.46)	0.001	0.310
rs1058205	5	28830/28999	0.75 (0.64–0.88)	<0.001	<0.001	0.78 (0.65–0.93)	<0.001	0.007	0.58 (0.42–0.81)	<0.001	0.001	0.74 (0.62–0.88)	<0.001	0.001	0.63 (0.48–0.82)	<0.001	0.001
Ethnicity/source of control
Caucasian/PB	4	28562/28701	0.77 (0.65–0.91)	<0.001	0.002	0.80 (0.66–0.98)	<0.001	0.027	0.58 (0.41-0.82)	<0.001	0.002	0.76 (0.63–0.92)	<0.001	0.005	0.63 (0.47–0.83)	<0.001	0.001

*P value of Q-test for heterogeneity test (Pheter)

Figure 4

Forest plot of GG vs. AA genetic model of KLK3 rs2735839 polymorphism

Figure 5

Forest plot of GG vs. AA genetic model of KLK3 rs266882 polymorphism

Serum expression of KLK3, sensitivity analysis, and in silico analysis

Serum expression of KLK3 was investigated by ELISA. Our research show evidence that serum KLK3 levels in PCa participants with CC/CT genotypes were relatively lower than those with TT genotypes (P = 0.021, Figure 6). We measured sensitivity by removing each study once a time in every genetic model for KLK3 rs1058205, rs2735839, and rs266882 variants. Results of sensitivity analysis indicated that no individual study could influence the pooled ORs (Figure 7). In silico analysis showed that KLK3 expression in PCa tissue was relatively higher than that in control tissue (Figure 8A). Furthermore, we assessed the correlation between KLK3 expression and overall survival (OS) or disease-free survival (DFS) time of PCa patients. Kaplan–Meier estimate identified no significant difference for either OS time (P = 0.67, Figure 8B) or DFS time (P = 0.094, Figure 8C).

Figure 6

Analysis of serum KLK3 levels in rs1058205 T>C genotype of PCa cases with mean values by ELISA (horizontal lines, mean values)

Figure 7

Sensitivity analysis about KLK3 rs1058205 (A), rs2735839 (B), and rs266882 (C) polymorphisms and PCa risk

Figure 8

In silico analysis of KLK3 expression

(A) Relative expression of KLK3 in PCa tissue and paracancerous tissue. The association between KLK3 expression and overall survival time (B) or disease free survival time (C) in prostate cancer patients.

Publication bias

The Egger’s test and Begg’s funnel plot were used to assess the publication bias within the included studies. No obvious evidence of publication bias was indicated in rs1058205 (C-allele vs. T-allele, t = −0.13, P = 0.909; CC vs. TT, t = −0.13, P = 0.912; CT vs. TT, t = −0.12, P = 0.916; CC + CT vs. TT, t = −0.12, P = 0.916; CC vs. CT + TT, t = −0.13, P = 0.912), rs2735839 (G-allele vs. A-allele, t = 0.57, P = 0.626; GG vs. AA, t = −0.05, P = 0.964; GA vs. AA, t = 0.92, P = 0.456; GG + GA vs. AA, t = 0.47, P = 0.682; GG vs. GA + AA, t = 0.01, P = 0.992), and rs266882 (G-allele vs. A-allele, t = −1.57, P = 0.361; GG vs. AA, t = −1.74, P = 0.332; GA vs. AA, t = −1.66, P = 0.345; GG + GA vs. AA, t = −1.71, P = 0.336; GG vs. GA + AA, t = −1.65, P = 0.346). The shape of the funnel plots seemed asymmetrical in the homozygote comparison for the three KLK3 polymorphisms, which means no publication bias was found (Figure 9).

Figure 9

Begg’s funnel plot of standard error for assessing publication bias of rs1058205 (A), rs2735839 (B), and rs266882 (C) under RR vs. WW genetic model

Discussion

To our knowledge, the pathogenesis and etiology of PCa is still inconclusive. However, the continuous development in molecular biology and immunology methods can provide a basis for the diagnosis of PCa and evaluation of its prognosis [28]. The KLK gene family, consisting of 15 genes spanning in a region of approximately 300 kb on 19q13.4, encodes a cluster of serine-proteases in the human genome [29,30]. KLK3 encodes a PSA, which is mainly expressed in the prostate tissue and is useful in the diagnosis and monitoring of PCa [31,32]. Previous studies have demonstrated the association between KLK3 rs2735839 G allele and higher PSA levels in various ethnicities. Nevertheless, the association between KLK3 polymorphism and PCa susceptibility remains controversial [33,34]. Our meta-analysis primarily aimed to combine previous studies to provide more clear conclusions.

The results of the overall analysis indicated that KLK3 rs1058205 T>C polymorphism was associated with decreased PCa risk, especially in the Caucasians. No association was detected between the KLK3 rs2735839 A>G polymorphism and PCa risk in all comparisons. In addition, no association was observed between rs266882 A>G variant and PCa risk. In population-based studies stratified by source of control, a significant association between the rs1058205 T>C polymorphism and PCa was observed. Furthermore, we assessed the serum expression of KLK3 by ELISA and indicated that individuals who were CC/CT carriers of KLK3 rs1058205 variant had a lower serum expression level than TT carriers. This result shows evidence that mutation from T to C at KLK3 rs1058205 locus may cause lower expression of KLK3. To further demonstrate the influence of KLK3 variants on PCa susceptibility, we used online gene expression mini database to assess the expression of KLK3 in PCa tissues. It indicated that KLK3 expression in PCa tissue was higher than that in normal tissue, which was in line with the results of our study.

Some factors can limit the generalization of the above results. The insufficient number of reported cases on KLK3 variations when specifying different ethnic backgrounds of PCa is one of the main limitations of the present study. Although all eligible studies based on the inclusion criteria were collected, the sample sizes of these studies were not large enough to provide adequate statistical power to assess the association between these polymorphisms and PCa susceptibility in further stratified analysis. The quantitative synthesis of several subgroups may have no sufficient testing power to precisely evaluate the true association. Histological types and Gleason score may also have varied between the evaluated studies. Moreover, the impact of environmental factors such as radiation, toxin, diet, and infectious agent were not evaluated in the analyzed studies. The impact of these factors on susceptibility to PCa should be further evaluated in further experiments. Due to the complex effects of the polygenic nature of complex diseases, including PCa, it is possible that gene-to-gene interactions play a more vital role in an individual’s susceptibility to this disease than single polymorphisms [35]. Furthermore, previous studies showed evidence that manuscript with ‘positive’ results was prone to be accepted faster than that with ‘negative’ findings, which could take longer time to be published in time-lag bias [36].

In spite of these limitations, some key advantages in the present analysis should also be acknowledged. Five case–control studies with 28,830 PCa cases and 28,999 control subjects related to the KLK3 rs1058205 variant, four case–control studies with 5394 cases and 5389 controls concerning the rs2735839 polymorphism, and three case–control studies with 1614 cases and 1981 controls related to rs266882 polymorphism were pooled in our present analysis. Statistical powers of these studies were strongly increased. Furthermore, publication bias was not identified in any of the genetic models, indicating that the results are relatively dependable and stable. We performed stratification analysis and sensitivity analysis to explore the sources of heterogeneity in the present study. Our data indicated that neither ethnicity nor source of controls was the sources of heterogeneity. Sensitivity analysis showed evidence that our data regarding these KLK3 polymorphisms were generalizable and credible.

In summary, our present study showed that rs1058205 polymorphism of KLK3 was a risk factor for PCa development, polymorphism T>C of rs1058205 was associated with decreased susceptibility to PCa particularly in the Caucasian population. Further researches on additional gene–environment interactions are warranted to provide a better understanding on the association between KLK3 polymorphisms and PCa risk.