A total of 90 patients were enrolled and randomized: FT group (n = 30; 17 female/13 male), FL group (n = 30; 16 female/14 male), and ZD group (n = 30; 15 female/15 male). Baseline characteristics were comparable across groups, with no significant differences in gender distribution (p = 0.875), mean age (FT: 25.3 ± 3.8 years; FL: 26.1 ± 3.6 years; ZD: 26.9 ± 3.7 years; p = 0.251), Body Mass Index (BMI, FT: 23.0 ± 3.4 kg/m2; FL: 22.4 ± 3.8 kg/m2; ZD: 21.9 ± 3.5 kg/m2; p = 0.518), or baseline AO severity grades (all patients Grade 2 or 3; p = 0.870). All 90 participants completed the 6-month follow-up (Table 3). Typical treatment procedure and recovery of MNRF for AO are provided (Figure 3A–C). This representative patient (baseline PGS grade 3) achieved a good response (grade 1 at 6 months), with complete resolution of post-treatment erythema and no scarring.
Efficacy rates
Clinical success rates differed significantly among the three groups at 6 months post-treatment (p < 0.001, Figure 4). The FL (93.3%) and ZD (90.0%) groups achieved markedly higher success rates than the FT group (53.3%). Pairwise comparisons confirmed significant differences between each active treatment group and the FT group (FL vs. FT: p < 0.001; ZD vs. FT: p = 0.002), while no difference was observed between FL and ZD (p = 1.000). Detailed efficacy grade distributions, RD, and RR are presented in Table 4. At 7 days post-treatment, efficacy rates were 90.0% (FT), 100% (FL), and 96.7% (ZD), with no significant intergroup difference (p = 0.108, Figure 4). Efficacy grade distributions were: FT group (excellent: 53.3%, good: 20.0%, fair: 16.7%, poor: 10.0%), FL group (excellent: 73.3%, good: 23.3%, fair: 3.3%, poor: 0%), and ZD group (excellent: 50.0%, good: 26.7%, fair: 20.0%, poor: 3.3%). While FL and ZD showed a higher proportion of favorable outcomes, the differences did not reach statistical significance (Table 4).
Sensitivity analysis
To assess the robustness of the primary endpoint definition, a sensitivity analysis was performed using a stricter definition of clinical success (excellent + good only, excluding fair). At 7 days post-treatment, the efficacy rates were 73.3% in the FT group, 96.7% in the FL group, and 76.7% in the ZD group (p = 0.038). Pairwise comparisons showed significant differences between FL and FT (p = 0.007) and between FL and ZD (p = 0.016). At 6 months, the efficacy rates were 33.3% (FT), 90.0% (FL), and 66.7% (ZD) (p < 0.001). Pairwise comparisons showed significant differences between FL and FT (p < 0.001), FL and ZD (p = 0.028), and FT and ZD (p = 0.010). These sensitivity analyses confirmed the robustness of the primary findings. Detailed results are presented in Supplementary Table 1.
Recurrence
Recurrence was defined as a return to a PGS grade equal to or higher than the baseline level, restricted to initial responders (patients with reduced grade at 7 days vs. baseline). Among initial responders, 6-month recurrence rates were 40.7% (11/27) in the FT group, 6.7% (2/30) in the FL group, and 6.9% (2/29) in the ZD group, with a significant overall difference (p = 0.009, Figure 4). Pairwise comparisons showed FL and ZD groups had significantly lower recurrence risks than FT (Table 4).
DLQI
Baseline DLQI scores were comparable across groups (FT: 8.0 [IQR: 6.0–16.0]; FL: 9.0 [IQR: 6.8–15.0]; ZD: 9.5 [IQR: 6.8–16.0]; p = 0.729). All groups demonstrated significant improvements in DLQI from baseline at both follow-up time points (all p < 0.001). At 7 days, median DLQI scores were 5.0 (IQR: 3.0–9.3) (FT), 3.0 (IQR: 1.8–6.3) (FL), and 3.4 (IQR: 2.0–6.3) (ZD). Median reductions from baseline were 4.0 (95% CI: 1.50–5.50) (FT), 5.0 (95% CI: 3.00–7.50) (FL), and 6.0 (95% CI: 4.00–8.00) (ZD), with no significant intergroup differences (p = 0.085, Figure 5 and Figure 6). At 6 months, sustained improvements were observed, with median DLQI scores of 5.5 (IQR: 2.0–8.0) (FT), 1.0 (IQR: 0.0–3.0) (FL), and 2.0 (IQR: 0.0–3.0) (ZD). Median reductions were 5.0 (95% CI: 0.50–6.00) (FT), 7.0 (95% CI: 5.50–11.00) (FL), and 7.5 (95% CI: 6.00–11.00) (ZD), with significant intergroup differences (p < 0.001). Post-hoc analysis confirmed that the FL and ZD groups had significantly lower DLQI scores than the FT group (both p < 0.001), while FL and ZD showed no difference (p = 0.948, Figure 5 and Figure 6).
Safety
Significant intergroup differences were observed in the incidence of pain, erythema, and swelling (all p < 0.05, Table 5).
Pain: The FL group had a higher pain incidence (60.0%, 18/30) than the FT group (20.0%, 6/30; RR = 3.00, 95% CI: 1.39–6.50; RD = 40.0%, 95% CI: 15.4%–64.6%; p = 0.002) and ZD group (23.3%, 7/30; RR = 2.57, 95% CI: 1.24–5.34; RD = 36.7%, 95% CI: 12.0%–61.4%; p = 0.004), with no difference between FT and ZD (p = 1.000, Figure 7). Median VAS pain score was higher in FL (1.5 [IQR: 0.0–2.3]) than in FT and ZD (both 0.0 [IQR: 0.0–0.0]; p < 0.001), with a median difference of 1.0 (95% CI: 0.00–2.00). Two FL patients required lidocaine patches, and all pain resolved within 7 days (Table 5).
Other complications: Incidence of erythema/swelling was 40.0% (FL), 16.7% (ZD), and 10.0% (FT) (p = 0.013), with FL showing a significant difference vs. FT (p = 0.007; RR = 4.00, 95% CI: 1.25–12.75; RD = 30.0%, 95% CI: 9.4–50.6%). Ecchymosis rates were 20.0% (FL), 6.7% (ZD), and 3.3% (FT) (p > 0.05). No hematoma was reported, and all mild complications resolved within 2 weeks with local heat application. Dysesthesia occurred only in FL (6.7%, 2/30) and resolved within 1 month with mecobalamin (Table 5, Figure 7).
Long-term safety (6 Months Post-treatment)
At the 6-month follow-up, all early complications (pain, erythema, swelling, and ecchymosis) had resolved. Scar formation was observed in 10.0% (3/30) of the FL group, consisting of small (< 1 cm), non-hypertrophic, asymptomatic lesions (Table 5, Figure 7). No severe adverse events (e.g., burns, infections, necrosis) were reported throughout the trial. All adverse events were mild and transient. In summary, the ZD needling strategy achieved 6 months efficacy (90.0%) comparable to FL (93.3%), with a safety profile closer to FT, while FT showed lower efficacy (53.3%) and higher recurrence (40.7%). The ZD strategy effectively balanced therapeutic efficacy and safety, offering a promising optimized protocol for MNRF treatment of axillary osmidrosis.
DATA AVAILABILITY:
The dataset supporting the findings of this study is publicly available via Science Data Bank at https://doi.org/10.57760/sciencedb.45827

Figure 1: CONSORT diagram showing the flow of participants during the study. Adapted from the CONSORT 2025 Flow Diagram template under a Creative Commons Attribution License10. Please click here to view a larger version of this figure.

Figure 2: Schematic of three needling strategies for MNRF treatment of axillary osmidrosis. (A, D) Full Transverse (FT) Group; (B, E) Full Longitudinal (FL) Group; (C, F) Zone-Directed (ZD) Group. (A–C) Intraoperative needling orientation; (D–F) Top-down view of needle trajectory distribution. This figure was created by the authors using Adobe Photoshop. No part of this figure has been reproduced or adapted from previously published sources. Please click here to view a larger version of this figure.

Figure 3: Typical treatment procedure and recovery of MNRF for axillary osmidrosis (AO). (A) Preoperative marking of the axillary treatment area (scale bar = 1 cm, grid size of 3 cm2 [approximately 2 cm × 1.5 cm]); (B) Immediately post-treatment, showing pinpoint bleeding and erythema; (C) At 6 months post-treatment, demonstrating complete resolution of erythema with no scar formation. Please click here to view a larger version of this figure.

Figure 4: Clinical efficacy and recurrence outcomes. (A) Clinical success rate at 7 days post-treatment; (B) Clinical success rate at 6 months post-treatment; (C) 6-month recurrence rate (restricted to initial responders, i.e., patients with reduced PGS at 7 days). Between-group comparisons were performed using the Chi-square test with Bonferroni correction. ***p < 0.001, **p < 0.01. Please click here to view a larger version of this figure.

Figure 5: Individual changes in DLQI scores from baseline. (A) At 7 days and (B) 6 months. Each bar represents the DLQI reduction for an individual patient, ordered by decreasing value. The dashed line indicates the minimally clinically important difference (MCID) for DLQI, above which the improvement is considered clinically meaningful. Please click here to view a larger version of this figure.

Figure 6: DLQI scores and improvements in the FT, FL, and ZD groups at baseline, 7 days, and 6 months. (A) DLQI scores at each time point; (B) DLQI improvement (ΔDLQI) from baseline to 7 days; (C) DLQI improvement (ΔDLQI) from baseline to 6 months. Between-group comparisons were performed using the Kruskal-Wallis test with Dunn’s post hoc test: *p < 0.05, **p < 0.01, ***p < 0.001. Please click here to view a larger version of this figure.

Figure 7: Incidence of complications at 7 days post-treatment in the FT, FL, and ZD groups. Data are presented as percentages. **p < 0.01. Please click here to view a larger version of this figure.
Table 1: The Park grading system of axillary osmidrosis. This grading system classifies the severity of axillary osmidrosis on a 0–3 scale. Higher grades indicate greater severity. Please click here to download this file.
Table 2: The efficacy grade of the axillary osmidrosis elimination grading system. Abbreviations: PGS = the Park grading system of axillary osmidrosis. Please click here to download this file.
Table 3: Baseline demographic and clinical characteristics of the participants. Intergroup comparisons for continuous variables (age, BMI) were performed using one‑way ANOVA, and for categorical variables (gender, PGS grade) using the χ2 test. A p value < 0.05 was considered statistically significant. Abbreviations: FT = full transverse; FL = full longitudinal; ZD = zone-directed; SD = standard deviation; PGS = the Park grading system for axillary osmidrosis; BMI = body mass index. Please click here to download this file.
Table 4: Comparison of short- and long-term efficacy, and recurrence among three groups.*p < 0.05 was considered statistically significant for overall comparisons among the three groups (Kruskal-Wallis H test/χ2 test). Post hoc pairwise comparisons (Dunn’s test/χ2 test) were performed when p < 0.05, with significance set at p’ < 0.017 (Bonferroni correction, α = 0.05/3). Abbreviations: FT = full transverse; FL = full longitudinal; ZD = zone-directed; RD = risk difference; RR = relative risk; CI = confidence interval; DLQI = dermatology life quality index; IQR = interquartile range. Please click here to download this file.
Table 5: Comparison of complications among different groups. *p < 0.05 was considered statistically significant for overall comparisons among the three groups (Kruskal-Wallis H test/χ2 test). Post hoc pairwise comparisons (Dunn’s test/χ2 test) were performed when p < 0.05, with significance set at p’ < 0.017 (Bonferroni correction, α = 0.05/3). Abbreviations: FT = full transverse; FL = full longitudinal; ZD = zone-directed; RD = risk difference; RR = relative risk; CI = confidence interval; IQR = interquartile range. Please click here to download this file.
Supplementary Table 1: Sensitivity analysis of short- and long-term efficacy among the three groups using a stricter efficacy definition (excellent + good only, excluding fair). *p < 0.05 was considered statistically significant for overall comparisons among the three groups (Kruskal‑Wallis H test/χ2 test). Post hoc pairwise comparisons (Dunn’s test/χ2 test) were performed when p < 0.05, with significance set at p' < 0.017 (Bonferroni correction, α = 0.05/3). Bold values indicate statistical significance. Abbreviations: FT = full transverse; FL = full longitudinal; ZD = zone-directed. Please click here to download this file.
Supplementary File 1: CONSORT checklist. This file contains the CONSORT checklist for the study. Please click here to download this file.