$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
This clinical study was conducted in compliance with the ethical principles outlined in the Declaration of Helsinki and was approved by the Ethics Committee of the Haiyan County People's Hospital(Approval No.: 2025-17).
Study participant assessment and selection
Clinical data were extracted for all patients diagnosed with BPH from the hospital's electronic medical record system, covering the period from January 2020 to March 2021.
PSM was employed using the 1:1 nearest-neighbor method to mitigate potential confounding factors and selection bias inherent in retrospective designs.
This matching process yielded a final cohort of 120 patients, with 60 patients in each group. Patients with missing data at any scheduled time point, including baseline, 6, 12, 24, or 36 months, were excluded from the analysis. Consequently, all 120 matched patients had complete follow-up data. As shown in Figure 1, the study workflow included patient screening, PSM, and final cohort allocation.
The following inclusion criteria were applied to screen potential participants. The diagnosis met the criteria of the Chinese Guidelines for the Diagnosis and Treatment of BPH. The patient's age was verified to be between 50 and 80 years, with an IPSS ≥ 8, a PSA level > 1.4 ng/mL, and a PV ≥ 40 mL as measured by transrectal ultrasound. The Qmax was confirmed to be ≤ 15 mL/s. Voluntary participation and signed informed consent were also verified.
The following exclusion criteria were applied. Patients with comorbid prostate cancer, urinary tract infection, neurogenic bladder, or bladder stones were excluded. Patients with a history of prostate surgery or prior treatment with 5α-reductase inhibitors or α1-blockers were excluded. Patients with severe cardiac, hepatic, or renal insufficiency, such as serum creatinine > 2.0 mg/dL, were excluded. Patients with a known allergy to terazosin hydrochloride or finasteride, as well as those with comorbid severe psychiatric illness or anticipated poor compliance, were also excluded.
The following dropout criteria were defined and applied after enrollment. Patients who did not ultimately meet all inclusion criteria or had critically incomplete information were excluded from the final analysis. Patients who demonstrated poor medication compliance or voluntarily withdrew from the study were excluded. Patients who experienced a sudden onset of other serious illnesses during the study period were excluded. Patients who developed severe psychological or emotional problems that prevented continued participation were excluded. Patients who underwent a change in treatment, such as the addition of a new surgical procedure, were also excluded.
The required sample size was calculated prior to matching. An a priori sample size calculation was performed using dedicated software. An independent-samples t-test model was used. The effect size (d) was set to 0.5, the significance level (alpha) was set to 0.05 (two-tailed), and the desired statistical power (1-β) was set to 0.8. The calculated sample size was inflated to account for an estimated 10% data attrition rate and potential losses during the PSM process. The initial calculation indicated a need for 54 participants per group, for a total sample size of 108. The final cohort was expanded to 120 individuals to meet statistical requirements.
Treatment protocol assignment and execution
Patient-specific treatment details were identified and compiled from Electronic Medical Record (EMR) prescription data.
For patients assigned to the combination therapy group, the administration of terazosin hydrochloride (2 mg) orally once daily at bedtime and finasteride (5 mg) orally once daily in the morning was confirmed from the records.
For patients assigned to the monotherapy group, the administration of finasteride (5 mg) orally once daily in the morning was confirmed from the records.
The treatment period was defined as 12 months, and the total follow-up period was defined as 36 months. Patients were instructed not to take any other medications affecting urinary function or sex hormones during the study.
Data collection for outcome measures
All study data were collected retrospectively from institutional electronic health records, laboratory databases, and nursing documentation systems.
Data for the primary outcome measures were collected at baseline, 6, 12, 24, and 36 months. Sex hormone levels, including dihydrotestosterone, serum testosterone, and estradiol, were extracted and recorded. Urinary control function metrics, including the number of urinary incontinence episodes within 24 h and the ICIQ-SF score, were also extracted and recorded.
Data for the secondary outcome measures were collected at the same time points: baseline, 6, 12, 24, and 36 months. The IPSS, Qmax, PVR, and PV were extracted and recorded. All reported adverse events, such as dizziness, hypotension, and sexual dysfunction, were extracted and recorded from adverse event report forms, along with relevant liver and kidney function indicators from laboratory reports.
Statistical analysis
All statistical analyses were conducted using dedicated software. To mitigate selection bias, 1:1 nearest-neighbor PSM without replacement was performed. The propensity score was estimated by a logistic regression model, with treatment assignment, defined as combination therapy versus monotherapy, as the dependent variable. The following baseline covariates were forced into the model based on clinical plausibility and prior literature: age, BMI, IPSS, Qmax, PVR, PV, PSA, C-reactive protein, total testosterone, estradiol, DHT, 24 h incontinence episodes, and ICIQ-SF score. Matching was performed on the logit of the propensity score using a caliper of 0.0224. Balance after matching was assessed using standardized mean differences (SMDs), with SMD < 0.1 indicating negligible imbalance. As shown in Table 1, all SMDs after matching were < 0.1, confirming adequate confounding control.
Continuous variables were tested for normality using the Shapiro-Wilk test. For baseline comparisons between the matched groups, normally distributed continuous variables were presented as mean ± standard deviation (x̄ ± s) and were compared using the paired t-test to account for the matched pairs. Categorical variables were presented as numbers and percentages [n (%)] and were compared using McNemar's test for paired binary data.
For longitudinal outcomes measured at multiple time points, including baseline, 6, 12, 24, and 36 months, linear mixed-effects models were employed, with treatment group, time, and their interaction as fixed effects, and subjects as random intercepts. This approach accounted for the matched design and within-subject correlation. When appropriate, post-hoc comparisons with Bonferroni correction were performed.
All statistical tests were performed as two-sided. A P-value of less than 0.05 was defined as the threshold for statistical significance.