Abstract
To elucidate the pathogenesis of prostate diseases, following in silico analysis, the LKB1 gene was selected for further investigation. The LKB1 gene has been associated with poor prognosis and is frequently mutated in different types of cancers. In this study, 50 benign prostatic hyperplasia (BPH) and 57 prostate cancer (PCa) tissues, including matched normal tissue for the patients, were analyzed by qRT-PCR and DNA sequencing for LKB1 expression and the mutation profile, respectively. Expression of LKB1 was increased in 60.7% of the PCa tissues compared with noncancerous tissue samples (p ≤ 0.001). However, LKB1 expression was lower when compared with normal tissues in BPH (p = 0.920). Four coding sequence alterations were detected in BPH. Three silent mutations were located in codons 9, 32, and 275 and a missense mutation was observed in codon 384. Six alterations were identified in the intronic regions of the LKB1 gene in both PCa and BPH. Five mutations were observed in both patient groups. A new alteration in intron 6 was observed in a patient with PCa. The LKB1 gene may be associated with benign transformations rather than the tumors in prostate pathogenesis when its expression and mutation status are considered. However, the mechanism of LKB1 in PCa needs further studies.
Introduction
Prostate cancer (PCa) is the second most common cancer in men and accounts for 15% of all cancers diagnosed (Ferlay et al., 2015). It is the second most common reason of death related to cancer after lung cancer with a rate of 10% (Siegel et al., 2014). Therefore, early diagnosis of PCa is of uttermost importance. An important marker used for screening of PCa is the prostate-specific antigen (PSA). Still, routine biopsy for levels of PSA ≥4 ng/mL in older men has a low diagnostic efficiency of 30% for PCa and provides an unreliable diagnostic marker (Catalona et al., 1994). PSA expression is increased during infections and inflammation in benign prostatic hyperplasia (BPH), and management of patients with high PSA levels has a higher rate of false negativity, resulting in repeated prostatic Tru-cut biopsies with transrectal ultrasonography and an increased rate of complications related to this invasive procedure (Loeb et al., 2013).
Some novel markers such as PCA3; prostate-specific membrane antigen (PSMA); and the TPMRSS2-ERG gene fusion have been suggested as possible biomarkers, but the utility of these as diagnostic and prognostic parameters is still not established and no single biomarker has ideal characteristics for detection and evaluation of PCa (Fiella et al., 2018). Therefore, for the diagnosis and management of prostatic diseases, new biomarkers that can selectively distinguish patients with BPH and PCa are essential. Identification of genes displaying differential expression in BPH may provide new insights into the mechanisms of prostatic disease and these genes may serve as potential candidates for future targets of therapy.
To figure out the principal genes and related pathways playing a role in the pathogenesis of prostate diseases, PubMed was text-mined using the Pathway Studio 9 program (Elsevier, Amsterdam, the Netherlands) and its internal ResNet 9 Mammalian database (Nikitin et al., 2003). Using gene set enrichment analyses, various interacting genes in curated pathways related to both PCa and benign prostate hyperplasia were determined. After the subnetwork enrichment analysis, genes were regrouped as those with known or unknown roles in PCa and benign prostate hyperplasia. Among these candidate genes with unclear regulatory function in prostate diseases, LKB1 was chosen for further analysis.
The LKB1 gene is frequently mutated in different types of cancers. Originally defined in patients with Peutz–Jeghers syndrome, LKB1 functions in multiple cellular processes, including regulation of the cell cycle, apoptosis, and metabolism (Hemminki et al., 1998a). Recent evidence indicates that LKB1 expression is associated with poor prognosis in cancer (He et al., 2014; Sun et al., 2016; Xiao et al., 2016). Although mutations of LKB1 have been investigated in different tumors, they have not been analyzed in human PCa. On the other hand, LKB1 expression in human prostate tumors has been investigated only in a single study (Lu et al., 2015).
In this study, we aimed to determine the differential expression and sequence alterations of the LKB1 gene in both PCa and BPH to clarify its role, if any, in the pathophysiology of PCa and to investigate its potential role as a diagnostic marker.
Methods
Samples of 50 BPH and 57 PCa tissues, including matched normal tissue controls, were obtained from BPH and clinically localized PCa patients. The study was approved by the Cerrahpasa Medical Faculty Ethics Committee of Istanbul University-Cerrahpasa (Approval No.: 83045809-604.01.02-161870), in accordance with ethical standards laid down in the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.
Tissue samples for BPH were taken from patients by either transurethral resection of the prostate (TUR-P) or open prostatectomy. Tissue samples for PCa were obtained by open radical retropubic prostatectomy. The mean ages of the patients with BPH and PCa at the time of surgery were 69.26 ± 8.93 and 59.96 ± 6.19, respectively. Bipolar TUR-P was performed in 40%, monopolar TUR-P in 30%, and open prostatectomy in 30% of the patients with BPH. In patients with PCa, the pathological T stages were pT2 in 95% and pT3 in 5%. Patients who underwent open prostatectomy had significantly higher prostate volumes and PSA levels compared with TUR-P and radical prostatectomy groups (p < 0.001). Diabetes mellitus (DM) was present in 48% and 51% of patients with BPH and PCa, respectively. The respective rates for intake of oral metformin were 61% and 59%. Seventeen percent and 23% of patients with BPH and PCa, respectively, were smokers. All patients with PCa had adenocarcinoma. The mean preoperative PSA level was 6.78 ± 2.97 ng/mL and the postoperative mean level was 0.065 ± 0.19 ng/mL (medians were 6 and 0.001, respectively). The mean volume ratio of cancerous tissue within the prostate gland was 24.96% ± 17.45% (median 23%) in patients with malignant tumors. Clinicopathological characteristics of patients with localized PCa are shown in Table 2.
RNA isolation and cDNA preparation
Total cellular RNA was obtained from both tumor and nontumor samples using the PureLink® RNA Mini Kit (Thermo Fisher Scientific, Inc., Waltham, MA). Nucleic acid yields were measured spectrophotometrically using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Inc.). cDNA was synthesized using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Inc.) and 300 ng of total RNA in a reaction volume of 20 μL. The G6PDH (glucose-6-phosphate dehydrogenase) gene was used as the reference gene for normalization. Five microliters of cDNA was used for each real-time PCR, which was performed in a 20-μL reaction volume containing 1 × PowerUp SYBR Green Master Mix (Thermo Fisher Scientific, Inc.) and gene-specific primers [LKB1 (forward 5′-AACCGGCCAAGAGGTTCT-3′ and reverse 5′-GATGGGCACTGGTGCTTC-3′) and G6PDH (forward 5′-ATGCCTTCCATCAGTCGGATACA-3′and reverse 5′-ATAGCCCACGATGAAGGTGTTTTC-3′)].
The target and reference genes were coamplified in the same reaction in different wells. Amplification was performed at 95°C for 10 s, 60°C for 30 s, and 72°C for 1 s for 40 cycles with an initial denaturation step at 95°C for 10 min, using the LightCycler 480 II platform (Roche Diagnostics, Basel, Switzerland). The expression levels were calculated by the 2−ΔΔCt method (Schmittgen and Livak, 2008).
LKB1 coding sequence analysis
Genomic DNA isolated from the tumor and normal tissues was used to amplify all exons and exon–intron boundaries of the LKB1 gene. Amplification reactions were performed using primers previously described (Onozato et al., 2007). The reaction contains 200–300 ng of genomic DNA, 10 pmol of each primer, 2.5 mM of MgCl2, 200 μM of dNTP mix, 1 × PCR buffer, and 1 U of Taq polymerase (Thermo Fisher Scientific, Inc.) in a final volume of 50 μL. Following an initial denaturation step at 95°C for 10 min, all exons of the LKB1 gene were amplified for 35 cycles at 95°C for 30 s, 61°C for 30 s, and 72°C for 30 s, with a final elongation step at 72°C for 10 min. After electrophoresis on 2% agarose gels and visualizing under ultraviolet light for confirmation of the products, sequencing was performed using an Applied Biosystems 3100 automated DNA sequencer (ABI, Pleasanton, CA) and the ABI Prism BigDye Terminator cycle sequencing ready reaction kit. Sequencing was performed twice in both forward and reverse directions for the mutation-positive samples.
Statistical analyses
Statistical analyses were performed using the IBM SPSS Statistics software (IBM Corp. Released 2012. IBM SPSS Statistics for Windows, version 21.0. Armonk, NY: IBM Corp.). We used the unpaired χ2 test to analyze the differences in expression levels in the tumor and normal tissues for both genes. The correlation of LKB1 expression with clinicopathologic parameters such as histologic grade, age, gender, and prostate volume was analyzed using Spearman's rho and chi-square tests. ROC analyses were performed for diagnostic discrimination of BPH and PCa. p < 0.05 was considered statistically significant.
Results
When we evaluated the expression level of the LKB1 gene in PCa tumor tissues compared with matched noncancerous tissues, we observed highly significant differences. The LKB1 gene expression was upregulated in tumor tissues compared with normal tissues (Table 1, p ≤ 0.001). LKB1 overexpression was observed in 60.7% of tumor tissues compared with noncancerous tissue samples. The average increase in the LKB1 mRNA level in tumor tissues was 40% compared with corresponding normal tissues. When the association of LKB1 expression levels with clinicopathological parameters was analyzed, the LKB1 gene expression levels did not show any significant difference according to age, DM, metformin use, smoking status, prostate volume, PSA levels, Gleason group, and lymphatic/vascular/perineural invasion (Table 2).
Expression Level of the LKB1 Gene in Tissue Samples with Prostate Cancer and Benign Prostatic Hyperplasia
Statistical analysis was performed by using the Wilcoxon signed rank test.
BPH, benign prostatic hyperplasia; PCa, prostate cancer; SD, standard deviation.
Clinicopathological Characteristics and LKB1 Expression in Patients with Prostate Cancer
Statistical analyses were performed by using the chi-square test.
DM, diabetes mellitus; PSA, prostate-specific antigen.
In contrast to PCa, expression of the LKB1 gene was lower in BPH tissues compared with their matched normal counterparts. The decrease in expression of LKB1 was 20%. Subgroup analyses of patients with BPH, considering the presence of age, DM, metformin use, smoking status, prostate volume, and PSA levels, did not reveal any significant difference with regard to the LKB1 gene expression (Table 3).
Clinicopathological Characteristics and LKB1 Expression in Patients with Benign Prostatic Hyperplasia
Statistical analyses were performed by using the chi-square test.
To determine the mutation status of the LKB1 gene in PCa or BPH patients, we screened the LKB1 gene by direct sequencing. Although we identified four different coding sequence alterations in samples from patients with BPH, no mutations were detected in any PCa tissue samples (Table 4). Three of these alterations were silent mutations occurring as a result of G to A, C to G, and G to A transitions at the third position of codons 9, 32, and 275, respectively. The remaining mutation was a C to G transition at the first nucleotide of codon 384 and resulted in a tryptophan to arginine substitution. In addition to coding sequence alterations, we detected six different alterations in intronic sequences of the LKB1 gene in PCa and BPH. Five of these alterations were common in both patient groups with different frequencies. Additionally, we observed a T to A transition in intron 6 of the gene in the tumor from a patient with PCa (Table 5).
Genetic Alterations in the Coding Regions of the LKB1 Gene in Benign Prostatic Hyperplasia
Genetic Alterations in the Noncoding Regions of the LKB1 Gene
Discussion
Early diagnosis of PCa is essential as the disease is the second most common cancer in men as well as the second most common cause of cancer deaths after lung cancer (Siegel et al., 2014; Ferlay et al., 2015). Currently, PSA is an important marker used for screening of PCa. However, especially for low PSA levels, clinically critical PCa can be missed when using the PSA strategy (Thompson et al., 2004). In addition, some tumors showing neuroendocrine differentiation produce low levels of PSA (Okotie et al., 2007). Thus, better and more specific markers other than PSA and PSA density are needed to identify PCa with low false negativity rates (Mettlin et al., 1994; Wever et al., 2010). In this study, we investigated the utility of the LKB1 gene in PCa and BPH.
One of the main characteristics of cancer cells is unlimited replication potential, which results in high energy requirement. The energy sensor of the cell is adenosine monophosphate-activated kinase (AMPK) and its activation protects cells against physiological and pathological stress. Upon activation, it switches on the metabolism to generate ATP while blocking energy expenditure (Oakhill et al., 2009). Therefore, AMPK plays an important role in regulating energy hemostasis, cell proliferation, senescence, and differentiation (Hardie et al., 2012). LKB1 activates AMPK through phosphorylation (Carling, 2006).
LKB1 has been accepted as a tumor suppressor gene because it is frequently mutated in the cancer-prone PJS (Hemminki et al., 1998b). LKB1 mutation frequency has been investigated in various types of cancers. However, the results are controversial and mutation types and rates are different in various cancers (Avizienyte et al., 1998; Su et al., 1999; Kim et al., 2004; Kenanli et al., 2010; Tigli et al., 2013; Yalniz et al., 2014). The mutation rate of the LKB1 gene in PCa has not yet been investigated. In a recent study using comprehensive genomic profiling, 1356 formalin-fixed, paraffin-embedded tissue sections of PCa samples were analyzed among other solid tumors. Alterations were identified in 18 PI3K pathway-associated genes, including LKB1. Although genetic alterations were detected in 46% of PCa samples, consistent with our results, none of these alterations were located in the LKB1 gene (Milis et al., 2019).
In another study, Ikediobi et al. (2006) investigated mutations of 24 known cancer genes in the NCI-60 cell line set and a frame-shift deletion (p.K178fs*86) has been identified in the LKB1 gene in the Du145 PCa cell line. Except those database analyses, there is no study in the literature analyzing mutations of the LKB1 gene in PCa. Therefore, we investigated the LKB1 mutation frequency in prostate tissue samples from patients with BPH or PCa. Although four coding sequence alterations were identified in BPH samples, we did not observe any alteration in PCa tumor samples. Previous reports indicate that LKB1 may function as a tumor promoter or tumor suppressor in a context and/or tissue-dependent manner (Baas et al., 2003; Boudeau et al., 2003; Dorfman and Macara, 2008; Zhan et al., 2012). In accordance with this, in previous studies, we have detected the M22I mutation in 20% of HNSCC tumor samples, while this mutation was never observed in bladder or renal cancers (Kenanli et al., 2010; Tigli et al., 2013; Yalniz et al., 2014).
The most frequent alteration was S19X in renal and bladder tumor samples with a frequency of 97.1% and 81%, respectively (Tigli et al., 2013; Yalniz et al., 2014). The absence of LKB1 coding sequence alterations in PCa tumor tissues suggests that it does not function as a tumor suppressor in prostate carcinogenesis. This is in accordance with a previous study that has associated LKB1 with Akt activation in head and neck cancers (Zhong et al., 2008). The PI3K/Akt pathway is known to promote tumorigenicity and progression of PCa, and more than 90% of prostate carcinoma cases display Akt activation (Graff et al., 2000; Malik et al., 2002; Kreisberg et al., 2004; Cao et al., 2013). On the other hand, alterations at the exon/intron boundaries are common in all tumor types with different frequencies.
The strength of the present study lies in simultaneous analysis of the mutations and expression of the LKB1 gene in PCa. Most of the mutations identified by our study are located within the catalytic domain or the binding region for SMARCA4/BRG1, which is an essential member of chromatin remodeling complexes and has been shown to play an important role in prostate carcinogenesis (Marignani et al., 2001; Muthuswami et al., 2019; Cyrta et al., 2020).
We also investigated LKB1 mRNA expression levels in association with LKB1 mutations in both PCa and BPH samples. Our results show that expression of LKB1 is repressed in similar fractions of the PCa and BPH tissues (40% and 42.5%, respectively). However, the mean mRNA expression level of LKB1 was higher in PCa tumor tissues than their matched normal counterparts, while lower in BPH samples. In accordance with the absence of mutation, high expression rates of LKB1 mRNA support the notion that LKB1 does not function as a tumor suppressor in PCa development and/or progression. In accordance with these results, recent evidence suggests that LKB1 may also exert an oncogenic effect (Martínez-López et al., 2012; Delgado et al., 2019). Activation of LKB1 has been found to be essential for cancer cell survival (Lee et al., 2015; Peart et al., 2015; Trapp et al., 2017). Knockdown of LKB1 leads to loss of cell viability and clonogenicity in breast cancer (Ng et al., 2012) and inhibits cell proliferation, while promoting apoptosis in hepatocellular cancer (Martínez-López et al., 2012). Upregulation of LKB1 has been reported in hepatocellular cancer and was associated with malignant characteristics, recurrence, and worse survival (Martínez-López et al., 2010). Similar data have also been reported for breast cancer (Chen et al., 2016), and LKB1 expression has been associated with poor prognostic and clinical factors (Syed et al., 2019). In line with these data in the literature, it has been suggested that LKB1 may exert different functions in different types of cancers (Tan et al., 2018).
Recently, a study investigating TCGA data from 498 PCa patients suggested that PCa can be divided into six subgroups displaying different types of genetic alterations involving the MSH1, FOXA1, ATM, BRCA1, RB1, LKB1, P53, and AR genes and leading to different metabolic expression profiles (Srihari et al., 2018). Among these, the C3 (LKB1, RB1, and HOXB13) subgroup displayed the highest deregulation for multiple pathways, was associated with worse prognosis, and had decreased LKB1 gene expression. However, this group did not represent the whole cohort with aberrant LKB1 gene expression, and no data were presented in the study on LKB1 gene expression levels in all patients with PCa. Interestingly, about one-third of the active surveillance patients were reassigned in subgroup C3 and this subgroup was found to have a 30–40% relapse rate within the first 72 months.
There is only one study in the literature investigating LKB1 mRNA expression in human prostate tumor tissues (Lu et al., 2015). In contrast to our results, this study reports a significant decrease in the LKB1 mRNA expression in PCa tissues. This may be due to using different calculation methods and ethnic origins of the patients. They also reported associations between low expression and creatinine level, advanced clinical stage, and PSA concentration. We were not able to detect any association between LKB1 expression and the clinicopathological characteristics.
On the other hand, there is another report in the literature investigating LKB1 expression in PCa at the protein level (Xu et al., 2014). According to this report, LKB1 protein levels were much lower in prostate tumor tissues than in normal tissue. The authors also reported that LKB1 mRNA was downregulated in two PCa cell lines (Du-145 and PC-3) compared with a control cell line (RWEP-1). However, it is known that the LKB1 gene is homozygously deleted in the Du-145 cell line. It is also well known that LKB1 shuttles between the nucleus and cytoplasm and its cellular location plays an important role in its catalytic activity. However, Xu et al.'s report evaluates the total cellular level of LKB1 protein. Therefore, it is not possible to evaluate its tumorigenic effect in that context.
At the same time, the LKB1 protein can be modified post-translationally by phosphorylation, acetylation, or binding to different isoforms of the pseudokinase STRADα, etc. (Marignani et al., 2007). Although investigating expression of the LKB1 protein to confirm mRNA expression levels would corroborate our results, unfortunately most of the samples made available for this study were not sufficient to analyze LKB1 protein expression by Western blotting in matched pairs of tissue specimens. However, for LKB1, a high degree of agreement between the mRNA and protein expression has been shown in PCa (Lu et al., 2015) as well as in cervical (Lao et al., 2014), pancreatic (Yang et al., 2015), and breast (Chen et al., 2016) cancers. A relevant study compared LKB1 protein expression in six benign prostatic tissue samples from cystoprostatectomy surgical specimens and 22 PCa tissue samples from radical retropubic prostatectomy surgical specimens (Grossi et al., 2015). Full-length LKB1 protein was shown in all benign samples and an immunoblotting analysis detected lower LKB1 expression in PCa. More than half of the PCa tissues showed no staining of LKB1, irrespective of grade (Grossi et al., 2015). On the other hand, we also investigated the LKB1 expression level in BPH patients. In contrast to PCa, the LKB1 mRNA expression was downregulated in these samples.
Conclusion
In conclusion, our data suggest that LKB1 is associated with a benign transformation rather than tumor formation. This is the first report on the LKB1 gene in PCa, which investigated its expression and alterations in association with clinicopathologic characteristics. We believe that our data may act as a reference, providing valuable information for future research. However, the activity of LKB1 is controlled by too many different factors, and there is a need for further functional studies to analyze the LKB1 gene in PCa in larger patient cohorts.
Footnotes
Authors' Contributions
H.K. carried out the expression analyses and participated in tissue collection and analysis and interpretation of data. A.C. and G.G. participated in the genotyping assays and expression analyses. N.D. contributed to the design and coordination of the study and to the final manuscript. H.O. provided tissue samples and clinical data. N.B. conceived the study, participated in its design and coordination, interpreted the data, and wrote the manuscript. All authors have read and approved the final manuscript.
Disclosure Statement
No competing financial interests exist.
Funding Information
This study was supported by the Istanbul University Research Fund (TDK-2016-22213).
