The study protocol was reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Chongqing Medical University (No. 2024-020-01). The study was prospectively registered with the National Medical Research Registry, China (registration number: MR-50-24-036087; registration date: September 6, 2024). Written informed consent was obtained from all participants before enrollment, including consent for the publication of de-identified facial photographs. The study was conducted in accordance with the Declaration of Helsinki.
Study design
This was a 12‑week, split‑face, randomized, prospective clinical study. The compound acid peel used in this study consisted of AHA (4%), salicylic acid (2%), and azelaic acid (1.5%), with a total acid concentration of approximately 7.5%. The pH of the formulation was 3.6. Prior to initiating the study, the allocation of participants to the treatment and control arms was determined using a random sequence generated in a spreadsheet by a blinded statistician and implemented by a blinded doctor. The treatment side received compound acid peeling twice weekly for 3 weeks (six sessions total). The treatment area was gently cleansed with a mild, non-oily cleanser and dried completely, then the peel was applied in a thin, uniform layer (~2 mL per hemiface). The formulation was applied in a standardized sequence from the forehead and temples to the cheeks, nose, perioral region, and chin, while avoiding the eyelids, lips, mucosal surfaces, and damaged skin. The product was left in contact with the skin for 15 min, then neutralized with water. The clinical endpoint was mild, uniform erythema with tolerable tingling or warmth and without excessive burning, edema, blistering, or extensive frosting. After removal, a soothing and moisturizing mask was applied for 15 min.
Participants were instructed to avoid irritating topical products, exfoliation, picking, and direct sun exposure and to use broad-spectrum sunscreen during the daytime. This 3‑week treatment phase was followed by a 9‑week follow‑up period. The control side received no intervention during the first 3 weeks. After this initial waiting period, the control side was given the same 3‑week compound acid peeling course, followed by a 6‑week follow‑up. Five sessions per hemiface were administered by trained medical personnel at the hospital, whereas one session per hemiface was performed by the participant at home under standardized instructions provided by the investigators. The participant was provided with standardized written and verbal instructions regarding the amount and distribution of the product, the application sequence and treatment area, the prescribed exposure time, and the removal procedure. A schematic overview of the study design and experimental workflow is presented in Figure 1.
Patient recruitment
Thirty patients aged 18–60 years with a clinical diagnosis of moderate acne vulgaris (AV) were enrolled in the study. Eligible participants were of any sex and met the diagnostic criteria for moderate AV, corresponding to grade 2 (31–50 lesions, including comedones, papules, and pustules) or grade 3 (51–100 lesions, including numerous papules and pustules, occasional large inflammatory lesions, and fewer than three nodules). Patients were excluded if they had photosensitive skin or had used photosensitizing medications within the preceding 4 weeks; had a history of intense sun exposure within the preceding month or anticipated substantial sun exposure within 1 month after treatment; had scarring, active infection, or non-healing wounds at the treatment site; had a diagnosed immunodeficiency disorder; had severe cardiac, hepatic, or renal dysfunction; or had a documented psychiatric illness. All 30 eligible participants completed the study and were included in the final analysis.
Outcome assessments
Objective evaluations
Erythema and porphyrin analysis
Under standardized environmental conditions (fixed temperature, humidity, and lighting), a facial imaging system was used to capture images of the bilateral forehead, cheeks, jawline, and other affected areas at baseline and at each follow-up. This system allowed objective quantification of erythema intensity and porphyrin levels.
Microbiota sampling and immunohistochemistry
Skin microbiome profiling was performed on microbial swabs collected from the treated side before treatment and at week 3. Microbial DNA was extracted from the swabs with the specified extraction kit, following the manufacturer’s protocol, and the recovered DNA was then quantified. The bacterial 16S rRNA region was amplified with the universal primers 341F (5’-CCTACGGGNGGCWGCAG-3’) and 805R (5’-GACTACHVGGGTATCTAATCC-3’). Amplification began with denaturation at 98 °C for 30 s, followed by 32 cycles of 98 °C for 10 s, 54 °C for 30 s, and 72 °C for 45 s. A final extension was carried out at 72 °C for 10 min.
PCR amplicons were purified using the specified beads and quantified prior to library assessment. Product quality was evaluated using a bioanalyzer and quantitative PCR kit, after which eligible amplicons were pooled for sequencing. The demultiplexed raw reads were processed by first removing primer sequences with Cutadapt. Paired-end reads were assembled in FLASH, and fqtrim was used to exclude reads with quality scores below 20, lengths below 100 bp, or ambiguous bases accounting for more than 5% of the sequence. Clean tags that passed quality filtering were retained for downstream analysis. Chimeras were removed with Vsearch, and DADA2 was used to denoise the remaining reads and generate amplicon sequence variants (ASVs).
Species annotation was completed in QIIME2 with the feature-classifier plugin against the SILVA and NT-16S databases. QIIME2 was also used to calculate alpha and beta diversity. Bacterial composition was evaluated from relative-abundance profiles. Differences in genus-level abundance were tested using the Wilcoxon test, with P < 0.05 considered significant. Linear discriminant analysis effect size was performed with nsegata-lefse, using LDA ≥3.0 and P < 0.05 as the thresholds.
Among the enrolled participants, 10 patients were randomly selected for immunohistochemical analysis of inflammatory mediators, including IL-1α, IL-6, IL-17, TGF-β, and TLR2. Skin biopsy samples (2 mm) were collected from active acne lesions at baseline (week 0) and week 3 after disinfection with iodophor. Samples were fixed in 4% paraformaldehyde, embedded in paraffin, and cut into 4‑μm sections. Paraffin-embedded sections (4 μm) were deparaffinized and rehydrated. Antigen retrieval was performed using EDTA buffer (pH 9.0) with microwave heating. Endogenous peroxidase was quenched with 3% H2O2, after which sections were treated with 3% goat serum to reduce nonspecific antibody binding. The sections were then incubated with the primary antibodies overnight at 4 °C: IL-1α (1:400), IL-6 (1:400), IL-17A (1:400), TGF-β (1:400), and TLR2 (1:300). Subsequently, sections were incubated with a Polymer-HRP secondary antibody kit for 30 min at room temperature, followed by DAB chromogen detection and hematoxylin counterstaining. Negative controls were incubated with PBS instead of primary antibodies. Images were acquired using a digital scanner at 400× magnification. Five fields were randomly selected per section, and scoring was performed by two blinded dermatologists.
Subjective evaluations
Investigator global assessment
At baseline and week 12, a dermatologist blinded to treatment allocation assessed acne severity over the entire face using the Global Acne Grading System (GAGS).
Patient self‑assessment
At the end of the treatment and follow-up period, participants rated their perceived clinical improvement using a four-point self-assessment scale. The scores were defined as follows: 0, no improvement or worsening; 1, mild improvement; 2, moderate improvement; and 3, obvious improvement.
Blinded outcome assessment
To reduce assessment bias, a dermatologist blinded to treatment allocation independently evaluated clinical improvement using baseline photographs.
Safety monitoring
Adverse events were evaluated at every follow‑up visit to assess treatment safety.
Statistical analysis
Before study initiation, the sample size was estimated to detect a change in the whole-face GAGS score from baseline to week 12. The sample size calculation was based on a pilot study of compound acid peeling, assuming a standardized within-participant effect size of 0.60, a two-sided significance level of 0.05, and 80% power, 24 evaluable participants were required. Allowing for 20% attrition, the target enrollment was 30 participants.
Statistical analysis was conducted using statistical software and graphing software (see Table of Materials). Continuous data were described according to their distribution. Variables that were approximately normally distributed were reported as the mean ± standard deviation, whereas non-normally distributed variables were summarized as the median with interquartile range or range, as appropriate. Categorical data were expressed as counts and percentages.
Data distribution was evaluated using the Shapiro-Wilk test and supported by visual review of histograms and Q-Q plots. Because the sample size was small and several variables did not follow a normal distribution, the primary analyses were performed using nonparametric tests. Primary outcomes, including GAGS scores, red area, and porphyrin levels, were analyzed using the paired Wilcoxon signed-rank test.
Effect sizes and 95% confidence intervals were calculated for the primary clinical outcomes. For paired nonparametric comparisons, the matched-pairs rank-biserial correlation was used as the effect-size measure. Spearman's rank correlation was applied to evaluate relationships between continuous or ordinal variables. Statistical tests were two-sided, and results were considered significant at P < 0.05.
Microbiome analyses were performed on treatment-side samples collected from the same participants. Alpha diversity was evaluated using the Shannon and Simpson indices. Changes in alpha-diversity indices between baseline and post-treatment samples were assessed using the paired Wilcoxon signed-rank test.
The bacterial community composition was summarized by relative abundance at the phylum, family, and genus levels. Taxa with low abundance or low prevalence were retained for descriptive analysis. Differences in the relative abundance of specific taxa between baseline and post-treatment samples were evaluated using paired nonparametric tests. For microbiome analyses, taxon-level testing was exploratory, and unadjusted P values are reported. All tests were two-sided, with a significance threshold of P < 0.05.