Study design and participants
This randomized, patient- and assessor-blinded, parallel-group trial adhered to the CONSORT guidelines10. The checklist is provided in Supplementary File 1. The protocol was approved by the Ethics Committee of the First Affiliated Hospital of Chongqing Medical University (2023-435) and registered at ClinicalTrials.gov (NCT07094009). The study was conducted between October 2023 and April 2024 in the dermatology outpatient department. All procedures were performed in a dedicated treatment room maintained at a constant temperature of 22–24 °C. All participants provided written informed consent. Patients were randomly allocated in a 1:1:1 ratio to the FT, FL, or ZD group using a computer-generated randomization sequence. Allocation concealment was maintained using sequentially numbered opaque sealed envelopes prepared by an independent statistician. Participants and investigators, including outcome assessors and statisticians, remained blinded throughout the trial.
The primary outcome was the clinical success rate of different MNRF needling patterns in AO. Based on prior clinical observation, efficacy rates of 50% (FT), 90% (FL), and 90% (ZD) were assumed. Sample size calculation using PASS (two-tailed α = 0.05, power = 90%) determined 27 patients per group. To account for attrition, 30 patients were enrolled per group (total n = 90). Inclusion criteria: (1) age 18–50 years with AO; (2) AO severity Grade 2–3 per Park’s Grading System11; (3) willingness to provide informed consent. Exclusion criteria: (1) concurrent axillary dermatoses; (2) history of keloid or hypertrophic scarring; (3) pregnancy, lactation, or planned pregnancy; (4) coagulation disorders or immunodeficiency; (5) a previous history of invasive axillary odor treatment. Of 100 screened patients, 6 did not meet the inclusion criteria, and 4 declined, resulting in 90 randomized participants: 30 each in the FT, FL, and ZD groups (Figure 1).
Procedure of MNRF
On the day of treatment, the axillary region was thoroughly cleaned, and hair was shaved prior to the procedure. Each patient was positioned supine with the arm abducted at approximately 90° to expose the axillary vault. Following standard disinfection and sterile draping, tumescent anesthesia was administered at the dermal-subcutaneous junction. The anesthetic solution comprised 20 mL of 2% lidocaine, 0.5 mL of 0.1% epinephrine, and normal saline to a total volume of 200 mL. A waiting period of approximately 5 min was observed to ensure adequate anesthesia. The treatment area was defined as a square extending 1.0 cm beyond the hair-bearing region of the axilla and then subdivided into 3 cm2 grids. Each patient received a single treatment session. Treatment was delivered using an MNRF device equipped with a sterile, single-use tip. The device settings were fixed as follows: energy output at 15 W, pulse width at 300 ms, and needle depth at 4.0 mm. The needle penetration depth was adjusted to 4.0 mm to target the apocrine gland layer based on prior anatomical studies. The direction of needle insertion was accordingly defined as transverse (parallel to the basal skin folds) or longitudinal (perpendicular to them).
Patients were randomly assigned to one of three needling protocols: the full transverse (FT) group, in which transverse needling was applied across the entire treatment area (Figure 2A, D); the full longitudinal (FL) group, in which longitudinal needling was applied across the entire treatment area (Figure 2B, E); and the zone-directed (ZD) Group, in which a stratified approach was used (Figure 2C, F), with longitudinal needling restricted to the central region and transverse needling applied within a 1.5 cm peripheral band inside the marked border. All procedures were performed by a single experienced dermatologist to ensure consistency. Following treatment, the site was cleansed with normal saline, then a 30-minute ice compress was applied. The area was covered with sterile gauze and bandaged. Patients were instructed to keep the treatment area dry for 48 h and to avoid strenuous physical activity for 7 days. A detailed list of all materials, reagents, equipment, and software used in this study is provided in the Table of Materials. Generic institutional supplies (e.g., sterile gauze, disinfectant, and ice packs) were obtained from the hospital’s standard supply inventory.
Efficacy
Efficacy was assessed at 7 days and 6 months post-treatment. All evaluations were performed by two trained dermatologists who were blinded to group allocation. AO severity was graded on a 0–3 scale using the Park Grading System (PGS, Table 1). Treatment efficacy was determined by the change in the PGS score, which was classified as excellent, good, or fair based on a predefined AO elimination grading system (Table 2). Clinical success, defined as an efficacy grade of excellent, good, or fair, was based on published studies12,13. Recurrence was defined as a decrease in PGS score at 7 days followed by a return to baseline or higher at 6 months. The Dermatology Life Quality Index (DLQI, range 0–30) was recorded at baseline, 7 days, and 6 months.
Safety
Safety outcomes were documented at 7 days (for acute, transient events such as pain, erythema, swelling, and ecchymosis) and 6 months (for long-term sequelae such as scarring and pigmentation, as well as the status of any early complications that had resolved) according to predetermined checklists. Pain was assessed using a Visual Analog Scale (VAS, range 0–10, where 0 = no pain and 10 = worst imaginable pain).
Statistical analysis
Statistical analysis was conducted using SPSS. Normality of continuous variables was assessed with the Shapiro–Wilk test. Normally distributed data are expressed as mean ± standard deviation (SD) and were compared using one‑way ANOVA with Bonferroni post‑hoc tests. Non‑normally distributed data are presented as median (interquartile range, IQR) and analyzed using the Kruskal–Wallis H test followed by Dunn’s post‑hoc tests. Within‑group comparisons were performed using the paired t‑test (parametric) or the Wilcoxon signed‑rank test (non‑parametric). Categorical data are shown as frequencies (percentages) and compared with the Chi‑square or Fisher’s exact test, as appropriate. For multiple comparisons among the three groups, the Bonferroni correction was applied consistently to all primary and secondary endpoints, with a significance threshold of p < 0.017. All analyses were conducted on an intention‑to‑treat (ITT) basis. Treatment effects for categorical outcomes were expressed as Risk Difference (RD) and Relative Risk (RR), while the Hodges‑Lehmann median difference was used for non‑parametric continuous data. All effect sizes are reported with 95% confidence intervals (95% CIs). A two‑tailed p‑value < 0.05 was considered statistically significant for all other analyses. Exact p‑values are reported wherever feasible; p‑values below 0.001 are presented as p < 0.001. Post-hoc sensitivity analyses using a stricter efficacy definition (excellent + good only) were performed to test the robustness of the primary results.