Quantitative analysis of human kallikrein gene 14 expression in breast tumours indicates association with poor prognosis.

KLK14 (formerly known as KLK-L6) is a recently identified member of the human kallikrein gene family. This family harbours several genes aberrantly expressed in various cancers as well as established (PSA/hK3, hK2) and potential (hK6, hK10) cancer markers. Similar to other kallikrein genes, KLK14 was found to be regulated by steroid hormones, particularly androgens and progestins, in breast and ovarian cancer cell lines. Preliminary studies indicated that KLK14 is differentially expressed in breast, ovarian, prostatic and testicular tumours. Given the above, we determined the prognostic significance of KLK14 expression in breast cancer. We studied KLK14 expression in 178 histologically confirmed epithelial breast carcinomas by quantitative reverse transcription-polymerase chain reaction and correlated with clinicopathological variables (tumour stage, grade, histotype etc.) and with outcome (disease-free survival and overall survival), monitored over a median of 76 months. KLK14 mRNA levels ranged from 0 to 1,219 arbitrary units in breast cancer tissues, with a mean+/-s.e. of 136+/-22. An optimal cutoff value of 40.5 arbitrary units was selected, to categorise tumours as KLK14-positive or negative. Higher concentrations of KLK14 mRNA were more frequently found in patients with advanced stage (III) disease (P=0.032). No statistically significant association was found between KLK14 and the other clinicopathological variables. KLK14 overexpression was found to be a significant predictor of decreased disease-free survival (hazard ratio of 2.31, P=0.001) and overall survival (hazard ratio of 2.21, P=0.005). Cox multivariate analysis indicated that KLK14 was an independent prognostic indicator of disease-free survival and overall survival. KLK14 also has independent prognostic value in subgroups of patients with a tumour size </=2 cm and positive nodal, oestrogen receptor and progestin receptor status. We conclude that KLK14 expression, as assessed by quantitative reverse transcription-polymerase chain reaction, is an independent marker of unfavourable prognosis for breast cancer.

Breast cancer is the most prevalent malignancy among women, accounting for 21% of all female cancers and ranking third overall when both sexes are considered (Parkin ). Although the increased use of screening for early disease diagnosis and the widespread administration of systemic adjuvant therapies have lead to a decline in mortality rates, breast cancer is still the leading cause of death from cancer in women, causing over 39 500 deaths in the US annually (Peto ; Jemal ).

Given the heterogeneous nature of breast carcinomas, much attention has been focussed on the identification of tumour associated molecular markers that reveal the biological profile of each tumour and ultimately aid in determining cancer risk, diagnosis, screening, prognosis, monitoring, management and prediction of therapeutic response in breast cancer patients (Duffy, 2001). Among the multitude of markers discovered are serine proteases, which participate in many aspects of carcinogenesis, including stimulating cellular growth, angiogenesis and the degradation of the extracellular matrix (Gottesman, 1990; Duffy, 1991). These functions are in accord with clinical studies demonstrating that the aberrant expression of certain serine proteases correlates with the invasiveness and metastasis of cancer cells and predicts poor prognosis in various malignancies (Herszenyi ; Kuhn ).

Human kallikreins are a subset of secreted serine proteases found in a wide range of tissues and biological fluids and implicated in diverse physiological and pathological processes (Diamandis ; Yousef and Diamandis, 2001). The kallikrein genes, denoted KLK1–KLK15, are located on chromosome 19q13.4 and encode for corresponding kallikrein enzymes, hK1–hK15 (Diamandis ; Yousef ). Accumulating evidence indicates that many members of this family are differentially expressed in certain malignancies, including prostate (Rittenhouse ; Magklara ; Barry, 2001; Yousef ; Diamandis ), testicular (Luo ), breast (Yousef ,c) and ovarian (Anisowicz ; Diamandis ; Kim ; Luo ; Magklara ; Obiezu ; Yousef ) cancers. Also, many kallikrein genes examined thus far are under steroid hormone regulation, further suggesting a role for these enzymes in endocrine-related tissues (Yousef and Diamandis, 2002). Additionally, PSA/hK3, is the best tumour marker available in clinical medicine for diagnosing and managing prostate cancer (Diamandis, 1998; Barry, 2001), hK2 is useful for certain subgroups of prostate cancer patients (Magklara ), hK6, 10 and 11, have recently emerged as potential serological epithelial ovarian cancer markers (Luo ; Diamandis ), and several others possess clinical utility as prognostic/predictive markers (Diamandis and Yousef, 2001).

Human kallikrein gene 14 (KLK14), formerly known as KLK-L6, is a recently identified member of the human kallikrein gene family (Yousef ). Structurally, this gene is formed of five coding exons and four intervening introns. The encoded protein, hK14, is a trypsin-like serine protease, translated as an inactive 251 amino acid preproenzyme precursor of about 27.5 kDa, of which 18 amino acids constitute the signal peptide and six amino acids the activation peptide. hK14 harbours the conserved catalytic triad characteristic of serine proteases and is highly homologous to other kallikreins, including PSA/hK3. KLK14 has a restricted tissue expression pattern and is found in the central nervous system as well as in endocrine-related tissues such as the uterus, ovary, thyroid and testis. Additionally, in situ hybridisation studies demonstrated that KLK14 is expressed by the secretory epithelial cells of benign prostate gland, prostatic intraepithelial neoplasia and malignant prostate cells (Hooper ). Preliminary studies have shown that KLK14 is down-regulated at the mRNA level in prostatic, testicular, ovarian and breast cancer tissues and in two breast cancer cell lines (Yousef ). Hormonal regulation studies in breast and ovarian cancer cell lines indicate that KLK14 expression is controlled by the androgen receptor in response to steroid hormones, particularly androgens and progestins (our unpublished data). Based on these collective findings, we hypothesised that KLK14 expression in malignant breast tissues may have prognostic/predictive value for patients with breast carcinomas.

MATERIALS AND METHODS

Study population

Tumour specimens from 178 consecutive patients undergoing surgical treatment for primary breast carcinoma at the Department of Gynecologic Oncology, University of Turin, Turin, Italy were analysed in this study. Patient age ranged from 25 to 87 years, with a median of 58 years (Table 1). Follow-up information (median follow-up period of 76 months) was available for 164 patients, among whom 60 (36%) had relapsed and 51 (31%) died. All tissue specimens were histologically confirmed and frozen in liquid nitrogen immediately after surgery.

Table 1

Distribution of numerical variables in the study

Clinical and pathological information documented at the time of surgery included clinical stage, grade, histology and size of the tumour, number of positive axillary nodes, steroid hormone receptor status and treatment strategy (Table 2). Tumour sizes ranged from 0.1 to 15 cm, with a median of 2.2 cm. Out of the 178 patients, 109 (61%) had ductal carcinoma, 28 (16%) lobular carcinoma while 37 (21%) possessed other histological types. The histotype was unknown for four patients. Patients with disease of all three Stages (I–III) and tumour grades (I–III) were represented in this study, with staging determined according to the TNM classification and grading in accordance to the Bloom-Richardson grading system (Bloom and Richardson, 1957). Oestrogen and progesterone receptor status was established as described by the European Organization for Research and Treatment of Cancer (Eortc, 1980). With respect to treatment, 33 (19%) received no adjuvant treatment, 83 (47%) received tamoxifen, while 59 (33%) received chemotherapy with or without tamoxifen. This study has been approved by the Institutional Review Board of the University of Turin.

Table 2

Associations between KLK14 statusa and other variables in 178 patients with epithelial breast carcinomas

Total RNA extraction and cDNA synthesis

Tumour tissues were minced with a scalpel, on dry ice, and transferred immediately to 2 ml polypropylene tubes. They were then homogenized and total RNA was extracted using Trizol™ reagent (Gibco–BRL) following the manufacturer's instructions. The concentration and purity of mRNA were determined spectrophotometrically. Two μg of total RNA was reverse-transcribed into first-strand cDNA using the Superscript™ pre-amplification system (Gibco–BRL). The final volume was 20 μl.

Quantitative real-time polymerase chain reaction (PCR) and continuous monitoring of PCR products

Based on the published genomic sequence of KLK14 (GenBank accession #AF161221), two gene-specific primers were designed (6F5: 5′-AGT GGG TCA TCA CTG CTG CT-3′ and 6R5: 5′-TCG TTT CCT CAA TCC AGC TT-3′). These primers spanned more than two exons to avoid contamination by genomic DNA.

Real-time monitoring of PCR reaction was performed using the LightCycler™ system (Roche Molecular Systems, Indianapolis, USA) and the SYBR green I dye, which binds preferentially to double-stranded DNA. Fluorescence signals are proportional to the concentration of the product and are measured at the end of each cycle and immediately displayed on a computer screen, permitting real time monitoring of the PCR reaction. The reaction is characterised at the point during cycling when amplification of PCR products is first detected, rather than the amount of PCR product accumulated after a fixed number of cycles. The higher the starting quantity of the template, the earlier a significant increase in fluorescence is observed. The threshold cycle is defined as the fractional cycle number at which fluorescence passes a fixed threshold above baseline.

Endogenous control

For each sample, the amount of KLK14 cDNA and of an endogenous control (β actin, a housekeeping gene) were determined using a calibration curve (see below). The amount of KLK14 was then divided by the amount of the endogenous reference, to obtain a normalised KLK14 value.

Standard curve construction

The full-length mRNA sequence of the KLK14 gene was amplified by PCR using gene-specific primers, and the PCR product was cloned into a TOPO TA cloning vector (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. A plasmid containing β-actin cDNA, was similarly prepared. Plasmids were purified using a mini-prep kit (Qiagen Inc., Valencia, CA, USA). Different standard curves for actin and KLK14 were constructed using serial dilutions of the plasmid. These standards were included in each run. The LightCycler software automatically calculates the standard curve by plotting the starting dilution of each standard sample vs the threshold cycle, and the sample concentrations are then calculated accordingly. Standards for both KLK14 and actin RNAs were defined to contain an arbitrary starting concentration, since no primary preparations exist. Hence, all calculated concentrations are relative to the concentration of the selected standard. Each sample was repeated twice to ensure reproducibility.

PCR amplification

The PCR reaction was carried out on the LightCycle™ system. For each run, a master mixture was prepared on ice, containing 1 μl of cDNA, 2 μl of LC DNA Master SYBR Green 1 mix, 50 ng of primers and 1.2 μl of 25 mM MgCl2. After the reaction mixture was loaded into the glass capillary tube, the cycling conditions were carried out as shown in Table 3.

Table 3

Experimental protocol used for quantitative PCR amplification of the KLK14 gene

Melting curve

For distinguishing specific from non-specific products and primer dimers, a melting curve was obtained after amplification by holding the temperature at 70°C for 30 s followed by a gradual increase in temperature to 98°C at a rate of 0.2°C per s, with the signal acquisition mode set at step. To verify the melting curve results, representative samples of the PCR products were run on 1.5% agarose gels, purified, and cloned into the pCR 2.1-TOPO vector (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. The inserts were sequenced from both directions using vector-specific primers with an automated DNA sequencer.

Statistical analysis

Patients were subdivided into groups based on different clinical or pathologic parameters (Table 2) and statistical analyses were performed using SAS software (SAS Institute, Cary, NC, USA). An optimal cutoff value was defined by χ2 analysis based on the ability of KLK14 values to predict the disease-free survival (DFS) and overall survival (OS) of the study population. According to this value, tumours were categorised as KLK14-positive or KLK14-negative and associations between KLK14 status and other qualitative variables were analysed using the χ2 or the Fisher's Exact test, where appropriate. The cutoff value for tumour size was 2 cm. Lymph node status was either positive (any positive number of nodes) or negative. Age was categorised into three groups: less than 45 years, 45–55 years and greater than 55 years. Survival analyses were performed by constructing Kaplan–Meier disease free survival (DFS) and overall survival (OS) curves and differences between curves were evaluated by the log-rank test (Mantel, 1966), as well as by estimating the relative risks for relapse and death using the Cox proportional hazards regression model (Cox, 1972). Cox analysis was conducted at both univariate and multivariate levels. Only patients for whom the status of all variables was known were included in the multivariate regression models, which incorporated KLK14 and all other variables for which the patients were characterised. The multivariate models were adjusted for KLK14 expression in tumours, patient age, tumour size, stage, grade, histological type and oestrogen receptor (ER) and progestin receptor (PR) status.

RESULTS

Relationship between KLK14 expression and other parameters

KLK14 mRNA levels ranged from 0 to 1219 arbitrary units in breast cancer tissues, with a mean±s.e. of 136±22. The cutoff point (40.5 arbitrary units; 69th percentile) indicated that 55 (31%) of the 178 breast tumour tissues were positive for KLK14 expression (Figure 1). Table 2 depicts the distribution of KLK14 expression, positive or negative, in breast tumour tissues in relation to other established prognostic factors such as menopausal status, tumour size, stage, grade, histological type, nodal status, steroid receptor status and adjuvant therapy. Patients with KLK14-positive breast tumours more frequently had advanced stage (III) disease (P=0.032). Significant associations between KLK14 expression and tumour size, grade and histology, or menopausal, nodal and steroid receptor status were not observed.

Figure 1

Determination of the optimal cutoff value for KLK14 expression. For details, see text.

Univariate and multivariate survival analysis

The strength of association between each clinicopathological variable and DFS and OS is displayed in Table 4. In univariate analysis, patients with KLK14-positive tumours had a significantly increased risk of relapse (decreased DFS) and death (decreased OS) (hazards ratios of 2.31 and 2.21; P=0.001 and 0.005, respectively). Importantly, when treated as a continuous variable, KLK14 expression retained its statistically significant association with decreased DFS and OS (P<0.001). Further, when survival outcomes were adjusted for all other variables in the multivariate analysis (i.e. Cox proportional hazard regression model), the adverse effects of KLK14 positivity on DFS and OS were preserved (hazards ratios of 2.14 and 1.99; P=0.009 and 0.029, respectively), implying that KLK14 expression is an independent prognostic indicator. As expected, Kaplan–Meier survival curves (Figure 2) demonstrate that patients with KLK14-positive tumours have shorter DFS (P<0.001) and OS (P=0.004) compared to those who are KLK14-negative.

Table 4

Univariate and multivariate analysis of KLK14 expression in relation to DFS and OS

Figure 2

Kaplan–Meier survival curves for disease-free survival (A) and overall survival (B) in patients with KLK14-positive and negative breast tumours. n=number of samples.

Univariate and multivariate survival analysis insubgroups of patients

We further examined the associations between KLK14 expression levels and survival outcomes in subgroups of patients stratified by tumour size, nodal, OR and PR status (Table 5). KLK14 expression significantly impacted survival in subgroups of patients with a tumour size ⩽2 cm and positive nodal, OR and PR status. Univariate analysis revealed that KLK14-positive patients with a tumour size ⩽2 cm were about three times more likely to suffer disease progression and four times more likely to die than KLK14-negative patients (P=0.004 and 0.007, respectively). These survival differences remained significant after the data were subjected to multivariate analysis (P=0.021 and 0.014, respectively). Among patients with positive nodal status, high KLK14 expression was associated with approximately two-fold greater risk of relapse and death in both univariate (P=0.007 and 0.017, respectively) and multivariate analyses (P=0.010 and 0.025, respectively). Similarly, there was a tendency for an increased risk of disease progression and death in both OR and PR positive patients with KLK14-positive tumours. In univariate analysis, hazard ratios derived from the Cox regression model in relation to DFS and OS were 3.40 (P=0.001) and 3.47 (P=0.004) respectively, for ER positive and 3.61 (P=0.001) and 3.17 (P=0.014) for PR positive patients. Under multivariate analysis, KLK14 retained its independent prognostic value in these subgroups of patients.

Table 5

Associations between KLK14 and DFS and OS in subgroups of patients stratified by tumour size, nodal status, oestrogen and progesterone receptor status

DISCUSSION

The optimal management of women with breast cancer involves a multidisciplinary approach, including the use of biological markers. Ultimately, the goal is to select the best marker or panel of markers that are most informative in terms of their ability to predict relapse, disease progression, survival and response to therapy. Traditional prognostic/predictive factors in breast cancer include tumour size, grade, nodal status, hormone receptor status, vascular invasion and age (Denley ). However, only hormone receptor status has predictive value in selecting patients who are likely to respond to therapy, and is the only marker recommended for routine use by the American Society of Clinical Oncology (Smith ) and the College of American Pathologists Consensus Statement (Fitzgibbons ).

Many other potential biological markers have been identified as prognostic factors including p53, c-myc, BCL-2, HER-2, vascular endothelial growth factor (VEGF), urokinase plasminogen activator (uPA), CA 15-3, BR 27.29 and cathepsin B and D (Maguire ; Duffy ; Scorilas ; Hamilton and Piccart, 2000; Bundred, 2001). Some markers are also predictive. For instance, expression of HER-2 has value in determining response to treatment and in selecting patients for Herceptin therapy (Ross and Fletcher, 1998). The multitude of new candidate biomarkers will likely lead to insights into the cellular changes that correlate with, or determine the different biological properties and diverse behaviour of individual breast tumours. Furthermore, stratifying patients based on the presence or absence of such markers may help to tailor different therapeutic strategies to meet individual needs (Clark ).

In this paper, we identify KLK14 as a new independent marker of unfavourable prognosis in breast cancer. Patients with KLK14-positive breast tumours were more likely to have advanced stage (III) disease. When assessing KLK14 expression in terms of predicting survival outcomes, we found an increased risk of relapse and death for patients with KLK14-positive tumours. This was also observed in subgroups of patients with a tumour size ⩽2 cm, positive nodal, OR and PR status. Hence, KLK14 expression may aid in predicting relapse, disease progression and/or survival in breast cancer patients.

Interestingly, the results obtained for KLK14 in this study are comparable to those obtained for KLK3, KLK6 and KLK15 in breast cancer (Yu ; Anisowicz ) and KLK4, KLK5 and KLK10 in ovarian cancer (Kim ; Luo ; Obiezu ), in that high expression of these kallikrein genes also correlated with patient prognosis. Similar to KLK14, these genes are also under steroid hormone regulation (Yousef and Diamandis, 2001). These observations allow us to speculate that multiple kallikreins may participate in a common enzymatic pathway that plays a role in the normal physiology of the breast. This pathway may be deregulated in breast carcinogenesis. As is the case with many other serine proteases, certain kallikreins may degrade the extracellular matrix promoting tumour invasiveness and metastasis.

It was recently realised that KLK14 is under steroid hormone regulation, particularly androgens and progestins, and that these effects are mediated through the androgen receptor (our unpublished data). It is also known that breast cancer is a hormone-dependent malignancy (Russo and Russo, 1998) and that androgen receptors are present in 70–90% of primary breast tumours (Soreide ), and 75% of breast cancer metastases (Lea ). Thus, it is likely that the androgen receptor, acting as a ligand-activated transcription factor, upregulates KLK14 gene expression during breast carcinogenesis.

The role of androgens in the aetiology of breast cancer is ill-defined. In vitro, they both stimulate (MCF-7, MDA-453) or inhibit (T-47D, ZR-75-1, MFM-223) the growth of AR-positive breast cancer cell lines (Hackenberg ; Birrell ). Animal studies have demonstrated that androgens shorten the latency period, enhance tumour size and increase the incidence of breast tumours, by promoting rather than initiating carcinogenesis in rodents (Liao ; Xie ). These effects may be mediated by androgen-regulated genes that directly function in cell growth regulation and through interaction of the AR with other transcription regulators, allowing for cross-talk with other growth pathways (Brys, 2000). Thus, the identification of androgen regulated genes, such as KLK14, may help to define new targets for breast cancer treatment. Such new treatments may be particularly important in metastatic disease, where the AR is often the sole steroid receptor expressed.

In summary, we quantified KLK14 expression in breast tumours and found that high KLK14 expression is associated with decreased DFS and OS in both univariate and multivariate analysis. Additional basic and clinical studies are required to delineate the activity of KLK14 in both the normal and malignant breast and to further define the clinical value of this biomarker.