Skip to main content

Main menu

  • Home
  • Current Issue
  • Archive
  • Info for
    • Authors
    • Editorial Policies
    • Advertisers
    • Editorial Board
  • Other Publications
    • Anticancer Research
    • In Vivo
    • Cancer Diagnosis & Prognosis
  • More
    • IIAR
    • Conferences
  • About Us
    • General Policy
    • Contact
  • Other Publications
    • Cancer Genomics & Proteomics
    • Anticancer Research
    • In Vivo

User menu

  • Register
  • Subscribe
  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Cancer Genomics & Proteomics
  • Other Publications
    • Cancer Genomics & Proteomics
    • Anticancer Research
    • In Vivo
  • Register
  • Subscribe
  • My alerts
  • Log in
  • My Cart
Cancer Genomics & Proteomics

Advanced Search

  • Home
  • Current Issue
  • Archive
  • Info for
    • Authors
    • Editorial Policies
    • Advertisers
    • Editorial Board
  • Other Publications
    • Anticancer Research
    • In Vivo
    • Cancer Diagnosis & Prognosis
  • More
    • IIAR
    • Conferences
  • About Us
    • General Policy
    • Contact
  • Visit iiar on Facebook
  • Follow us on Linkedin
Research ArticleArticles
Open Access

Combination of Urinary MiR-501 and MiR-335 With Current Clinical Diagnostic Parameters as Potential Predictive Factors of Prostate Biopsy Outcome

JAROSLAV JURACEK, MARIE MADRZYK, KAROLINA TRACHTOVA, MICHAELA RUCKOVA, JULIA BOHOSOVA, DOMINIK A. BARTH, MARTIN PICHLER, MICHAL STANIK and ONDREJ SLABY
Cancer Genomics & Proteomics May 2023, 20 (3) 308-316; DOI: https://doi.org/10.21873/cgp.20383
JAROSLAV JURACEK
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
MARIE MADRZYK
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
KAROLINA TRACHTOVA
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
MICHAELA RUCKOVA
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
JULIA BOHOSOVA
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOMINIK A. BARTH
2Research Unit of Non-Coding RNAs and Genome Editing in Cancer, Division of Clinical Oncology, Department of Medicine, Comprehensive Cancer Center Graz, Medical University of Graz, Graz, Austria;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
MARTIN PICHLER
2Research Unit of Non-Coding RNAs and Genome Editing in Cancer, Division of Clinical Oncology, Department of Medicine, Comprehensive Cancer Center Graz, Medical University of Graz, Graz, Austria;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
MICHAL STANIK
3Department of Urologic Oncology, Clinic of Surgical Oncology, Masaryk Memorial Cancer Institute, Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: stanik@mou.cz
ONDREJ SLABY
1Central European Institute of Technology, Masaryk University, Brno, Czech Republic;
4Department of Biology, Faculty of Medicine, Masaryk University, Brno, Czech Republic
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: ondrej.slaby@ceitec.muni.cz
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Background: The detection of prostate cancer (PCa) is currently based on prostate-specific antigen (PSA) quantification as an initial screening followed by ultrasound-guided transrectal biopsy. However, the high rate of false-negative biopsies often leads to inappropriate treatment. Therefore, new molecular biomarkers, such as urine microRNAs (miRNAs), are a possible way to redefine PCa diagnostics. Patients and Methods: Urine samples of 356 patients undergoing prostate biopsy (256 cases with confirmed prostate cancer, 100 cases with negative prostate biopsy) at the Masaryk Memorial Cancer Institute (Czech Republic) and additional 36 control subjects (healthy controls, benign prostatic hyperplasia – BPH) were divided into the discovery and validation cohorts and analyzed. In the discovery phase, small RNA sequencing was performed using the QIAseq miRNA Library Kit and the NextSeq 500 platform. Identified miRNA candidates were validated by the RT-qPCR method in the independent validation phase. Results: Using the small RNA sequencing method, we identified 12 urine miRNAs significantly dysregulated between PCa patients and controls. Furthermore, independent validation showed the ability of miR-501-3p and the quantitative miR-335:miR-501 ratio to distinguish between PCa patients and patients with negative prostate biopsy. The subsequent combination of the miR-335:miR-501 ratio with PSA and total prostate volume (TPV) using logistic regression exceeded the analytical accuracy of standalone parameters [area under curve (AUC)=0.75, positive predictive value (PPV)=0.85, negative predictive value (NPV)=0.51)] and discriminated patients according to biopsy outcome. Conclusion: Combination of miR-335:miR-501 ratio with PSA and total prostate volume was able to identify patients with negative prostate biopsy and could potentially streamline decision making for biopsy indication.

Key Words
  • Prostate cancer
  • diagnosis
  • prostate biopsy
  • microRNA
  • urine
  • non-invasive
  • NGS
  • miRNA ratio
  • regression model

Prostate cancer (PCa) is the second most frequently diagnosed cancer and the fifth leading cause of cancer-related death among men (1). A crucial event for timely and effective therapy is the early detection of the disease which is maintained primarily using the measurement of prostate-specific antigen (PSA) levels in the blood serum. Based on increased PSA levels patients are further indicated for transrectal ultrasound-guided biopsies (TRUS) of the prostate (2). Unfortunately, the PSA test shows low specificity, as elevated levels may also reflect inflammation or trauma. This leads to a significant number of PSA-positive men who are unnecessarily subjected to biopsy. Furthermore, even if patients are correctly referred for biopsy due to abnormal screening results, there is a high rate of false-negative TRUS results (3). For this reason, there is an effort to refine or replace PSA testing with novel PCa-specific biomarkers (2, 3).

Modern tumor biomarkers are described as non-invasive, easily measurable, highly specific, sensitive, inexpensive, and producing rapid results. Within prostate cancer diagnostics, urine represents an ideal source of biomarkers as prostate cancer cells release proteins, DNA, RNAs, other small molecules, and exosomes in a detectable amount. In the context of urogenital malignancies, one of the most promising classes of potential biomarkers meeting the described criteria are microRNAs (miRNAs) (4), and more importantly urinary miRNAs. These cell-free molecules possess remarkable stability and good analytical performance and can reflect the clinical status of the malignancy or characteristics of the tumor (5, 6). Despite a considerable number of studies focused on the importance of urinary miRNAs within prostate cancer detection, there is, however, a low overlap in potential biomarker miRNAs.

In our current study, we performed high-capacity analysis of urine miRNAs to identify potential non-invasive biomarkers for the early detection of prostate cancer using a small RNA sequencing approach. Within the independent validation, we employed a more clinically relevant design and evaluated the diagnostic power of selected miRNAs in a group of prostate cancer patients and patients with negative prostate biopsy, where we compared the analytical performance of miRNAs with current clinical diagnostic parameters.

Patients and Methods

Patient characteristics. Between May 2016 and June 2020, 356 men undergoing prostate biopsy due to either abnormal PSA values or digital rectal examination were recruited to participate in this prospective study at the Department of Urologic Oncology, Masaryk Memorial Cancer Institute (MMCI). The study has been approved by the Ethics Committee of MMCI (the Ethics approval number MOU171816) and all participants signed an informed consent. Patients with urinary tract infection or a history of cancer were excluded from the control group. In total, 392 samples were used for analysis in this study.

Urine samples were obtained after digital rectal examination (DRE) before performing a biopsy. The prostate core biopsies were examined by an expert uropathologist; in patients where prostate cancer was detected, the grading was assigned according to the International Society of Urological Pathology 2014 consensus guidelines. Three subgroups were histologically defined based on the pathological evaluation of biopsy samples. The subgroups of indolent, intermediate, and aggressive PCa corresponding to ISUP grade 1, grade 2, and grade 3-5, respectively.

In the discovery phase, the control group included 18 healthy men who were seen in the urology outpatient clinic for benign urological conditions or as part of a cancer prevention program with long-term follow-up and low PSA values below 1.5 ng/ml. In addition, 18 patients with histologically confirmed benign prostatic hyperplasia (BPH) were included. In the validation phase, urine samples from 100 patients with negative prostate biopsy were used as a control group. The clinicopathological characteristics of the cohorts are summarized in Table I. A schematic overview of the workflow is depicted in Figure 1.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table I.

Clinicopathological characteristics of prostate cancer patients and control groups.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Study flowchart. PCa: Prostate cancer; BPH: benign prostatic hyperplasia; NPB: negative prostate biopsy.

Sample collection and RNA isolation and purification. For all patients, the morning voided urine after DRE was collected into 15 ml tubes with EDTA as a nucleic acid preservative agent. The whole urine was then centrifuged at 4°C at 2,000 × g for 15 min and the supernatant was collected and stored at −80°C until further analysis. Before RNA isolation, the urine supernatant was centrifuged again at 4°C at 12,000 g for 15 min to separate cell fragments and debris. Finally, total RNA including small RNA and circulating and exosomal RNA was isolated from 1 ml of cell-free urine supernatant using the Urine microRNA Purification Kit (Norgen Biotek Corp., Thorold, ON, Canada) according to the manufacturer’s instructions. The quality and quantity of RNA were determined using the NanoDrop™ 2000 spectrophotometer and/or the Qubit Fluorometer (both Thermo Fisher Scientific, Waltham, MA, USA).

Small RNA sequencing and data processing. Sequencing libraries were prepared using a QIAseq™ miRNA Library Kit and QIAseq™ miRNA NGS 96 Index IL (both Qiagen, Hilden, Germany) strictly according to the manufacturer’s instructions. Adapter dilution and cycle numbers were set according to the user manual. In more detail, for 3′ ligation adapter was diluted 10-fold and the starting volume of undiluted template RNA was 5 μl. For 5′ ligation adapter was diluted 5-fold. In the library amplification step, the reaction was set according to the manual with the number of cycles programmed to 22. For the listed reactions, a ProFlex PCR System (Thermo Fisher Scientific) was used.

The quality and size profiles of the prepared libraries were analyzed using High Sensitivity D1000 ScreenTape and 2200 TapeStation (both Agilent, Santa Clara, CA, USA). The library concentration was measured by Qubit 2.0 Fluorometer using reagents provided in Qubit™ dsDNA HS Assay Kit (both Thermo Fisher Scientific). Each miRNA sequencing library was diluted with nuclease-free water to 4 nM, mixed at an equimolar ratio to create a library pool, loaded onto the NextSeq 500/550 High Output v2 Kit (75 cycles) reagent cartridge, and sequenced on NextSeq 500 System (Illumina, San Diego, CA, USA).

The raw sequencing images from Illumina NextSeq 550 were converted to the fastq format using bcl2fastq (version 2.20.0). Raw reads were quality-checked with the FASTQC package (version 0.11.8). Adapter sequences were identified by Kraken package (15-065) and trimmed using Cutadapt (version 1.18). Subsequently, low-quality bases (Phred <10) were trimmed with Cutadapt and reads shorter than 15 nt were removed. Collapsing was performed utilizing unique molecular identifiers (UMIs) with FASTX-Toolkit (version 0.0.14). Reads originating from snoRNAs, snRNAs, rRNAs, tRNAs, piRNAs, and YRNAs (downloaded from Ensembl and RefSeq databases) were identified using Bowtie (version 1.2.2) and removed from the data. The remaining reads were mapped against the miRBase (version 21) and quantified using miraligner tool (version 1.2.4). Statistical analysis was carried out in R (version 3.4.3) with DESeq2 package (version 1.18.1).

RT-qPCR. The individual miRNAs were analyzed by standard TaqMan™ MicroRNA Assays (Thermo Fisher Scientific). First, the reverse transcription was performed using the TaqMan™ MicroRNA Reverse Transcription Kit, followed by real-time PCR with TaqMan probes and TaqMan® Universal Master Mix II, no UNG (all Thermo Fisher Scientific). All reactions and thermal cycling conditions were compliant with the standard protocol provided by the manufacturer. All reactions, including no-template controls and interplate controls, were run in duplicates. The RT-qPCR reaction was performed on the QuantStudio 12K Flex Real-Time PCR system (Thermo Fisher Scientific). The average Ct values of measured replicates were normalized using hsa-miR-101-3p as the most suitable reference gene identified with GenEx™ software 6.0 (MultiD Analyses AB, Göteborg, Sweden). Data normalization was performed using the 2−ΔCt method where the ΔCt value corresponds to the difference between the Ct of the detected miRNA and the Ct of the reference gene. As an alternative normalization approach, the ratio of potential biomarker miRNAs was calculated by dividing the relative expression values of hsa-miR-501-3p and hsa-miR-335-5p (7).

Statistical analysis. Within the validation, the differential level of selected miRNAs between the compared cohorts was evaluated using the Mann-Whitney U-test and GraphPad Prism version 8.0.0 for Windows (GraphPad Software, San Diego, CA, USA, www.graphpad.com). The discriminatory power of each biomarker and/or their combination, as well as the cut-off value for optimal analytical utilization, was established using ROC analysis (GraphPad Software).

For a combination of individual clinical biomarkers used for prostate cancer diagnosis, a stepwise logistic regression was performed using JMP® software (JMP®, Version 16, SAS Institute Inc., Cary, NC, 1989-2021). For all analyzes, a p-Value <0.05 was considered statistically significant.

Results

Small RNA sequencing. During the analysis of the small RNA sequencing data, samples with low miRNA coverage were discarded. Subsequently, the statistical analysis identified 12 miRNAs significantly dysregulated between a group of PCa patients (20 patients ISUP 1-5) vs a group of control samples (17 BPH patients and healthy individuals), thereof 10 miRNAs with higher levels and 2 miRNAs with lower levels in the urine of PCa patients (Table II and Figure 2A). For further validation in the independent cohort, we selected miR-335-5p and miR-501-3p based on their analytical parameters (fold change, p-value, and the average of normalized counts).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table II.

List of significantly dysregulated miRNAs in prostate cancer patients vs. control samples (BPH + healthy controls).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Hierarchical clustergram discriminating prostate cancer patients (20 patients ISUP 1-5) and controls (17 healthy controls and BPH patients) according to differentially expressed urine miRNAs. The heatmap shows 12 significantly dysregulated miRNAs (p<0.05) (A). Analytical performance of standalone parameters (miR-335:miR-501 ratio, PSA, and TPV) and the regression model combining the parameters within the discrimination of PCa patients and patients with negative prostate biopsy. The regression model shows superior discriminatory power (AUC=0.7478) (B).

RT-qPCR analysis and clinical diagnostic parameters. In the validation cohort, we determined the urinary levels of miR-335-5p and miR-501-3p in cell-free urine samples of overall 220 prostate cancer patients and 100 patients with negative prostate biopsy as a control group. Of these two potential biomarkers, only the level of miR-501-3p was significantly dysregulated (p=0.02) in the urine of PCa patients compared to patients with negative prostate biopsy.

We performed the same analysis for available clinical diagnostic parameters: total PSA, free PSA (fPSA), prostate-specific antigen density (PSAD), and total prostate volume measured by transrectal ultrasonography (TPV). For average values of clinical diagnostic parameters within control and testing groups, see Table III.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table III.

Average values of clinical diagnostic parameters within control and testing groups.

Next, to better discriminate the compared groups, we applied a quantitative ratiometric approach and calculated the ratio of miR-335-5p and miR-501-3p levels. This miRNA ratio outperformed individual miRNAs and distinguished PCa patients from patients with negative prostate biopsy with higher discriminatory power (p=0.0004, AUC=0.64). For comparison, only PSAD (p=0.019, AUC=0.62) and total prostate volume (p<0.0001, AUC=0.66) were able to distinguish prostate cancer patients from controls among the monitored clinical diagnostic parameters (Table IV). Since PCa patients and the control group were both subjected to the prostate biopsy based on increased PSA level, this parameter could not significantly discriminate the cohorts (p=0.595).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table IV.

Analytical performance of standalone urine miRNAs and clinical parameters within the comparison of patients with negative prostate biopsies with PCa patients (ISUP 1-5) (I), patients with indolent PCa (ISUP 1) (II), and patients with clinically insignificant PCa (ISUP 1-2) (III). Additionally, analytical parameters were obtained from comparison of group combining patients with negative prostate biopsy with indolent PCa (NPB + ISUP 1) from significant PCa (ISUP 2-5) (IV), comparison of patients with indolent (ISUP 1) and aggressive (ISUP 3-5) prostate cancer (V) and comparison among ISUP grade groups (VI).

Additionally, both miR-501-3p and the miR-335:miR-501 ratio were able to distinguish indolent (p=0.0009 and p=0.0002, respectively) and clinically insignificant (p=0.0036, and p<0.0001, respectively) prostate cancer from controls. Indolent prostate cancer (ISUP 1) was distinguished also by the use of total PSA level (p=0.0039) and total prostate volume (p<0.0001) which significantly discriminated also clinically insignificant PCa (p<0.0001). Neither miR-501-3p, miR-335, nor the miRNA ratio could not discriminate clinically insignificant group combining patients with negative prostate biopsy with indolent PCa from significant PCa (ISUP 2-5).

Furthermore, miR-501-3p and the miR-335:miR-501 ratio (p=0.0039 and p=0.0105, respectively) as well as total PSA level and PSAD could discriminate between indolent (ISUP 1) and aggressive PCa (ISUP 3-5), which was later confirmed with the Kruskal-Wallis test showing association of all these parameters/potential biomarkers with the ISUP grade. All statistical comparisons are shown in Table IV.

Combination of urine miRNAs and clinical diagnostic parameters. Since individual biomarkers did not show satisfactory discriminative power, we decided to use nominal logistic regression and apply a combination of all clinical diagnostic parameters with and without individual miRNAs and the miRNA ratio. ROC analysis and whole model test significance were performed for all conditions. The most significant combinations of miRNA biomarkers and clinical parameters for each condition are summarized in Table V. Standalone miRNAs or the miRNA ratio could not discriminate group combining patients with negative prostate biopsy with indolent PCa from significant PCa, therefore no regression models were calculated for this condition.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table V.

Most significant regression models combining urine miRNAs and clinical parameters within the comparison of patients with negative prostate biopsy with PCa patients, patients with indolent PCa, clinically insignificant PCa and the comparison of indolent and aggressive PCa. For each model we state the AUC, PPV, NPV, and the cut-off value.

Prostate cancer patients vs. patients with negative prostate biopsy. For comparison of prostate cancer patients and patients with negative prostate biopsy, the highest increase in discriminatory power for miR-335 was achieved in combination with PSAD, PSA, and fPSA (AUC=0.7325). Only a minor improvement for miR-501 was reached in combination with PSAD, PSA, and fPSA (AUC for the model 0.7215). The highest discriminatory power was reached for the miR-335:miR-501 ratio when combined with PSA + TPV [AUC=0.7478, PPV=0.848, NPV=0.507; DxScore=3.55618451 + (−0.6526161 * miR-335:miR-501 ratio) + (0.03623058 * PSA) + (−0.028844 * TPV); Figure 2B; Table V]. For the PSA grey zone (patients with total PSA level between 4-10 ng/ml) this regression model reached the AUC of 0.7178. Among other models with the highest discriminatory power, we detected the miR-335:miR-501 ratio + PSA + TPV + PSAD or the miR-335:miR-501 ratio + TPV + PSAD (both conditions AUC=0.7476) and PSAD + PSA + fPSA (AUC=0.7475) (Table VI).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table VI.

Discriminatory power of regression models combining clinical parameters and urinary miRNAs to compare patients with negative prostate biopsies with PCa patients (ISUP 1-5) (I), patients with indolent PCa (ISUP 1) (II), patients with clinically insignificant PCa (ISUP 1-2) (III), and to compare patients with indolent (ISUP 1) and aggressive (ISUP 3-5) prostate cancer.

Indolent prostate cancer vs. patients with negative prostate biopsy. Next, we compared the ability of computed regression models to distinguish indolent ISUP1 PCa from patients with negative prostate biopsy. We observed a significant increase in discriminatory power for miR-335 when combined with PSA + TPV + fPSA (AUC=0.7548). Interestingly, every combination of miR-335 involving total prostate volume allowed to detect indolent PCa with AUC greater than 0.74. On the contrary, regression models based on the urinary level of miR-501 did not show a significant improvement from the standalone miRNA. The best analytical parameters were observed for the miR-335:miR-501 ratio combined with PSA + TPV + PSAD [AUC=0.7717; DxScore=2.75372875 + (−0.8142955 * miR-335:miR-501 ratio) + (−0.0915265 * PSA) + (−0.0252965 * TPV) + (4.04533464 * PSAD)], although more regression models reached the AUC of 0.75. Furthermore, when we evaluated significant regression models within the PSA grey zone, the additional refinement of the models was observed with AUC increasing to around 0.8. All comparisons and results are listed in Table VI.

Clinically insignificant prostate cancer vs. patients with negative prostate biopsy. Interestingly, in the comparison of clinically insignificant prostate cancer and patients with negative prostate biopsy, we observed the most marked increase in discriminatory power in combinations with miR-335. Specifically, the highest analytical parameters were detected in the regression model combining the miR-335 level with PSA + TPV + fPSA [AUC=0.7516; DxScore=2.34331463 + (−0.0001012 * miR-335) + (0.0693638 * PSA) + (−0.0384259 * TPV) + (−0.4683634 * fPSA)]. Subsequent analysis of these conditions within the PSA grey zone did not improve discriminatory capabilities (Table VI).

Indolent vs. aggressive prostate cancer. Finally, we verified the ability of identified miRNAs and clinical diagnostic parameters to differentiate between indolent and aggressive prostate cancer. The highest predictive power was detected for miR-501 in combination with PSA (AUC=0.7565) and analogously for the miR-335:miR-501 ratio combined with PSA [AUC=0.7619; DxScore=−1.5082927 + (0.54160878 * miR-335:miR-501 ratio) + (0.03907273 * PSA)] (Table VI).

Discussion

Prostate cancer detection is dependent on PSA as an initial screening test and ultrasound-guided transrectal biopsy as the gold-standard diagnostic method. However, about 70% of PCa cases subjected to biopsy turn out negative (8) causing unnecessary physical and psychological harm to the patients. Moreover, a high rate of false-negative biopsies (15-30%) leads to commonly repeated biopsies (9, 10) and undertreatment. Thus, the legitimate indication for biopsy remains an important clinical issue. Despite several biomarkers available for PCa detection, they are not included in routine practice. Therefore, new molecular biomarkers are a possible way to improve and streamline PCa diagnostics.

In this study, we performed a high-capacity urinary miRNA analysis to identify a set of miRNAs significantly dysregulated between PCa patients and healthy controls and BPH patients. In subsequent validation, we employed a clinically more relevant design and evaluated the discrimination power of the selected miR-335-5p and miR-501-3p in a group of prostate cancer patients and patients with negative prostate biopsy. In this cohort, only miR-501-3p was able to significantly distinguish patients with negative prostate biopsy (p=0.02), however, with rather unsatisfactory analytical parameters. To further refine the analytical accuracy of identified miRNAs, we applied the ratiometric approach and calculated the miR-335:miR-501 ratio. This approach emerges to be widely accepted, especially in body fluids, where stable reference genes might be challenging to find. The quantitative ratio does not require additional normalization and should be independent of the detection system. In fact, in our validation cohort the miR-335:miR-501 ratio outperformed standalone miRNAs and most clinical parameters in terms of distinguishing between PCa patients and patients with negative prostate biopsy (p=0.0004).

Both miR-335-5p and miR-501-3p were already shown to be associated with prostate cancer. Previously, Haldrup et al. identified new potential serum miRNAs including miR-501-3p able to identify PCa patients (11). More importantly, miR-501-3p was among miRNAs significantly dysregulated in urine exosomes from prostate cancer patients detected by deep sequencing analysis (12). As we employed a similar study design, these results could serve as independent validation of our NGS analysis data. Whole-genome miRNA profiling studies also described upregulation of miR-335-5p in the tissue of high grade PCa tumors (13) and in exosomes secreted by human PCa cells under hypoxic conditions (14). The only exploratory study describing miR-335-5p as a potential biomarker of prostate cancer is the work by Koh et al., where this miRNA was differentially expressed in urine, plasma, and tissue samples of PCa patients (15).

As for diagnostic tests or biomarkers, there are currently several alternatives to total PSA levels, such as free PSA, PSA density and velocity, PHI, 4Kscore, or PCA3. However, they are often of limited value and are not routinely used for economical or technical reasons. One possibility is to combine available parameters into one model which could significantly improve diagnostic performance. For example, Liu et al. used the multivariate model combining PSAD and multiparametric magnetic imaging (mpMRI), which performed significantly better than standalone parameters for the detection of PCa even within the PSA grey zone (3). Even more beneficial should be models that merge genomic biomarkers with clinical parameters. Adopting this approach, we used nominal logistic regression and tested a combination of available clinical diagnostic parameters (total PSA, free PSA, TPV, PSAD) with miR-335-5p and miR-501-3p and the miRNA ratio. Among all combinations, the highest discriminatory power was reached for the model based on the miR-335:miR-501 ratio with PSA and TPV which discriminated PCa patients according to the biopsy outcome with AUC of 0.7478, PPV=0.848, and NPV=0.507. Based on the NPV, therefore, our model can potentially reduce the number of unnecessary biopsies by 50%.

Additionally, we identified regression models discriminating patients with negative prostate biopsy from indolent (335:501 ratio + PSA + TPV + PSAD) and clinically insignificant PCa (miR-335 + PSA + TPV + fPSA). When comparing the patients with indolent PCa, the highest discriminating power reached combinations based on the miR-335:miR-501 ratio involving total prostate volume since this parameter alone allowed to distinguish compared groups (p<0.0001, AUC=0.69). However, the combination of markers significantly refined the analytical parameters (AUC=0.77). Moreover, an additional enhancement of the models was observed when evaluated within the PSA grey zone (AUC=0.8). For clinically insignificant prostate cancer, we observed the most marked increase in discriminatory power in combinations involving miR-335 and free PSA, despite the fact that these factors alone could not distinguish compared cohorts.

Lastly, several regression models based on the level of miR-501-3p or the miR-335:miR-501 ratio could differentiate between indolent and aggressive prostate cancer. A similar result was observed for miR-501-3p alone, whose level also decreased with the higher ISUP grade group. This is in agreement with the literature describing miR-501-3p as a potential tumor suppressor (16) that restricts prostate cancer cell growth by targeting CREPT to inhibit the expression of cyclin D1.

Our study has several potential limitations. For the identification of miRNA biomarkers, we used samples from healthy controls and BPH patients as a control group. Since a different control group was used in validation, this could lead to overlooking of miRNAs enabling better discrimination of patients with negative prostate biopsy. Given the high rate of false-negative biopsies, our data may be affected by this inaccuracy, as we do not possess detailed follow-up data showing whether the patient was indicated for re-biopsy or was diagnosed with PCa in time. The analytical power of identified regression models combining genomic biomarkers with clinical parameters could be improved by using other promising tests (PHI, PCA3, mpMRI), however, these are not routinely evaluated within the standard protocol at MMCI and could not be used. The weakness is also the one-level validation, which we preferred due to the increase in the homogeneity of the patient group. For these reasons, further independent studies are needed to confirm our results and demonstrate the potential clinical utility of the model. Another limitation of our study is a relatively high rate of positive biopsies (65%) in our patient cohort, which is higher when compared to the commonly observed positivity rates. This may reflect the specific pattern of patient referral to our tertiary hospital. Noteworthy, this fact may disturb the analytical outcomes for PSA in the prediction of positive prostate biopsy.

Conclusion

Current diagnostic procedures result in a large number of patients being misclassified for prostate biopsy. Using NGS analysis, we identified urine miRNAs for noninvasive detection of PCa. In independent validation, we successfully tested the ability of miR-335-5p and miR-501-3p and the miRNA ratio to discriminate PCa patients and patients with negative prostate biopsy. Next, we showed that combination of genomic biomarkers with clinical parameters could significantly improve diagnostic accuracy. Moreover, the combination of miR-335:miR-501 ratio + PSA + TPV was able to identify patients with negative prostate biopsy and could potentially improve decision making for biopsy indication.

Acknowledgements

This study was supported by Ministry of Health of the Czech Republic, grant no. NV18-03-00360 and by the project Excellence in research and development of non-coding RNA DIAGnostics in Oncology -RNADIAGON H2020-MSCA-RISE-2018- GA no. 824036. All rights reserved. The project National Institute for Cancer Research (Programme EXCELES, ID Project No. LX22NPO5102) – Funded by the European Union – Next Generation EU. The Core Facility Bioinformatics of CEITEC Masaryk University is gratefully acknowledged for obtaining of the scientific data presented in this paper.

Footnotes

  • Conflicts of Interest

    The Authors declare no conflicts of interest.

  • Authors’ Contributions

    Study conceptualization, J.J., M.S. and O.S.; Patients enrollment, samples characterization, M.S.; Laboratory analysis, J.J., M.M., M.R., and J.B.; Data analysis and statistical analysis, K.T., M.M. and J.J.; Original draft preparation and visualization, J.J.; Review and editing, all authors; Supervision, O.S. All Authors have read and agreed to the published version of the manuscript.

  • Received January 25, 2023.
  • Revision received February 19, 2023.
  • Accepted March 6, 2023.
  • Copyright © 2023, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY-NC-ND) 4.0 international license (https://creativecommons.org/licenses/by-nc-nd/4.0).

References

  1. ↵
    1. Sung H,
    2. Ferlay J,
    3. Siegel RL,
    4. Laversanne M,
    5. Soerjomataram I,
    6. Jemal A and
    7. Bray F
    : Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 71(3): 209-249, 2021. PMID: 33538338. DOI: 10.3322/caac.21660
    OpenUrlCrossRefPubMed
  2. ↵
    1. Descotes JL
    : Diagnosis of prostate cancer. Asian J Urol 6(2): 129-136, 2019. PMID: 31061798. DOI: 10.1016/j.ajur.2018.11.007
    OpenUrlCrossRefPubMed
  3. ↵
    1. Liu J,
    2. Dong B,
    3. Qu W,
    4. Wang J,
    5. Xu Y,
    6. Yu S and
    7. Zhang X
    : Using clinical parameters to predict prostate cancer and reduce the unnecessary biopsy among patients with PSA in the gray zone. Sci Rep 10(1): 5157, 2020. PMID: 32198373. DOI: 10.1038/s41598-020-62015-w
    OpenUrlCrossRefPubMed
  4. ↵
    1. Nam RK,
    2. Benatar T,
    3. Amemiya Y and
    4. Seth A
    : MiR-139 induces an Interferon-β response in prostate cancer cells by binding to RIG-1. Cancer Genomics Proteomics 18(3): 197-206, 2021. PMID: 33893074. DOI: 10.21873/cgp.20252
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Juracek J,
    2. Peltanova B,
    3. Dolezel J,
    4. Fedorko M,
    5. Pacik D,
    6. Radova L,
    7. Vesela P,
    8. Svoboda M,
    9. Slaby O and
    10. Stanik M
    : Genome-wide identification of urinary cell-free microRNAs for non-invasive detection of bladder cancer. J Cell Mol Med 22(3): 2033-2038, 2018. PMID: 29363887. DOI: 10.1111/jcmm.13487
    OpenUrlCrossRefPubMed
  6. ↵
    1. Zhou J,
    2. Gong G,
    3. Tan H,
    4. Dai F,
    5. Zhu X,
    6. Chen Y,
    7. Wang J,
    8. Liu Y,
    9. Chen P,
    10. Wu X and
    11. Wen J
    : Urinary microRNA-30a-5p is a potential biomarker for ovarian serous adenocarcinoma. Oncol Rep 33(6): 2915-2923, 2015. PMID: 25962395. DOI: 10.3892/or.2015.3937
    OpenUrlCrossRefPubMed
  7. ↵
    1. Avissar M,
    2. Christensen BC,
    3. Kelsey KT and
    4. Marsit CJ
    : MicroRNA expression ratio is predictive of head and neck squamous cell carcinoma. Clin Cancer Res 15(8): 2850-2855, 2009. PMID: 19351747. DOI: 10.1158/1078-0432.CCR-08-3131
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Lewicki P,
    2. Shoag J,
    3. Golombos DM,
    4. Oromendia C,
    5. Ballman KV,
    6. Halpern JA,
    7. Stone BV,
    8. O’Malley P,
    9. Barbieri CE and
    10. Scherr DS
    : Prognostic significance of a negative prostate biopsy: an analysis of subjects enrolled in a prostate cancer screening trial. J Urol 197(4): 1014-1019, 2017. PMID: 27836710. DOI: 10.1016/j.juro.2016.11.002
    OpenUrlCrossRefPubMed
  9. ↵
    1. Serefoglu EC,
    2. Altinova S,
    3. Ugras NS,
    4. Akincioglu E,
    5. Asil E and
    6. Balbay MD
    : How reliable is 12-core prostate biopsy procedure in the detection of prostate cancer? Can Urol Assoc J 7(5-6): E293-E298, 2013. PMID: 22398204. DOI: 10.5489/cuaj.11224
    OpenUrlCrossRefPubMed
  10. ↵
    1. Wolters T,
    2. van der Kwast TH,
    3. Vissers CJ,
    4. Bangma CH,
    5. Roobol M,
    6. Schröder FH and
    7. van Leenders GJ
    : False-negative prostate needle biopsies: frequency, histopathologic features, and follow-up. Am J Surg Pathol 34(1): 35-43, 2010. PMID: 19935058. DOI: 10.1097/PAS.0b013e3181c3ece9
    OpenUrlCrossRefPubMed
  11. ↵
    1. Haldrup C,
    2. Kosaka N,
    3. Ochiya T,
    4. Borre M,
    5. Høyer S,
    6. Orntoft TF and
    7. Sorensen KD
    : Profiling of circulating microRNAs for prostate cancer biomarker discovery. Drug Deliv Transl Res 4(1): 19-30, 2014. PMID: 25786615. DOI: 10.1007/s13346-013-0169-4
    OpenUrlCrossRefPubMed
  12. ↵
    1. Rodríguez M,
    2. Bajo-Santos C,
    3. Hessvik NP,
    4. Lorenz S,
    5. Fromm B,
    6. Berge V,
    7. Sandvig K,
    8. Linē A and
    9. Llorente A
    : Identification of non-invasive miRNAs biomarkers for prostate cancer by deep sequencing analysis of urinary exosomes. Mol Cancer 16(1): 156, 2017. PMID: 28982366. DOI: 10.1186/s12943-017-0726-4
    OpenUrlCrossRefPubMed
  13. ↵
    1. Walter BA,
    2. Valera VA,
    3. Pinto PA and
    4. Merino MJ
    : Comprehensive microRNA Profiling of Prostate Cancer. J Cancer 4(5): 350-357, 2013. PMID: 23781281. DOI: 10.7150/jca.6394
    OpenUrlCrossRefPubMed
  14. ↵
    1. Panigrahi GK,
    2. Ramteke A,
    3. Birks D,
    4. Abouzeid Ali HE,
    5. Venkataraman S,
    6. Agarwal C,
    7. Vibhakar R,
    8. Miller LD,
    9. Agarwal R,
    10. Abd Elmageed ZY and
    11. Deep G
    : Exosomal microRNA profiling to identify hypoxia-related biomarkers in prostate cancer. Oncotarget 9(17): 13894-13910, 2018. PMID: 29568403. DOI: 10.18632/oncotarget.24532
    OpenUrlCrossRefPubMed
  15. ↵
    1. Koh Y,
    2. Bustos MA,
    3. Moon J,
    4. Gross R,
    5. Ramos RI,
    6. Ryu S,
    7. Choe J,
    8. Lin SY,
    9. Allen WM,
    10. Krasne DL,
    11. Wilson TG and
    12. Hoon DSB
    : Urine cell-free microRNAs in localized prostate cancer patients. Cancers (Basel) 14(10): 2388, 2022. PMID: 35625992. DOI: 10.3390/cancers14102388
    OpenUrlCrossRefPubMed
  16. ↵
    1. Zhang Z,
    2. Shao L,
    3. Wang Y and
    4. Luo X
    : MicroRNA-501-3p restricts prostate cancer growth through regulating cell cycle-related and expression-elevated protein in tumor/cyclin D1 signaling. Biochem Biophys Res Commun 509(3): 746-752, 2019. PMID: 30621914. DOI: 10.1016/j.bbrc.2018.12.176
    OpenUrlCrossRefPubMed
Previous
Back to top

In this issue

Cancer Genomics - Proteomics: 20 (3)
Cancer Genomics & Proteomics
Vol. 20, Issue 3
May-June 2023
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Back Matter (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Cancer Genomics & Proteomics.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Combination of Urinary MiR-501 and MiR-335 With Current Clinical Diagnostic Parameters as Potential Predictive Factors of Prostate Biopsy Outcome
(Your Name) has sent you a message from Cancer Genomics & Proteomics
(Your Name) thought you would like to see the Cancer Genomics & Proteomics web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
2 + 8 =
Solve this simple math problem and enter the result. E.g. for 1+3, enter 4.
Citation Tools
Combination of Urinary MiR-501 and MiR-335 With Current Clinical Diagnostic Parameters as Potential Predictive Factors of Prostate Biopsy Outcome
JAROSLAV JURACEK, MARIE MADRZYK, KAROLINA TRACHTOVA, MICHAELA RUCKOVA, JULIA BOHOSOVA, DOMINIK A. BARTH, MARTIN PICHLER, MICHAL STANIK, ONDREJ SLABY
Cancer Genomics & Proteomics May 2023, 20 (3) 308-316; DOI: 10.21873/cgp.20383

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Reprints and Permissions
Share
Combination of Urinary MiR-501 and MiR-335 With Current Clinical Diagnostic Parameters as Potential Predictive Factors of Prostate Biopsy Outcome
JAROSLAV JURACEK, MARIE MADRZYK, KAROLINA TRACHTOVA, MICHAELA RUCKOVA, JULIA BOHOSOVA, DOMINIK A. BARTH, MARTIN PICHLER, MICHAL STANIK, ONDREJ SLABY
Cancer Genomics & Proteomics May 2023, 20 (3) 308-316; DOI: 10.21873/cgp.20383
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Patients and Methods
    • Results
    • Discussion
    • Conclusion
    • Acknowledgements
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

  • No related articles found.
  • PubMed
  • Google Scholar

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • Knockdown of G Protein-coupled Estrogen Receptor 1 (GPER1) Enhances Tumor-supportive Properties in Cervical Carcinoma Cells
  • The Role of Apoptotic Genes and Protein-Protein Interactions in Triple-negative Breast Cancer
  • A Multiplex Biomarker Assay Improves the Prediction of Survival in Epithelial Ovarian Cancer
Show more Articles

Similar Articles

Keywords

  • Prostate cancer
  • diagnosis
  • prostate biopsy
  • microRNA
  • urine
  • non-invasive
  • NGS
  • miRNA ratio
  • regression model
Cancer & Genome Proteomics

© 2023 Cancer Genomics & Proteomics

Powered by HighWire