Lycopus lucidus Extract Improves Experimental Polycystic Ovary Syndrome with Gut Microbiota and Metabolome Remodeling

Summary

ZeLan extract improved endocrine and metabolic phenotypes in DHEA-induced PCOS rats. Integrated 16S rRNA sequencing and serum metabolomics showed partial restoration of gut microbial alpha diversity, altered microbial composition, and metabolite changes involving steroid hormones, bile acids, and short-chain fatty acids.

Abstract

This study investigated the therapeutic effects of ZeLan extract on polycystic ovary syndrome (PCOS) and evaluated its association with gut microbiota-metabolome remodeling. Twenty-four female Sprague-Dawley rats were randomly divided into Control, PCOS Model induced by dehydroepiandrosterone (DHEA), and ZeLan treatment groups (n = 8 each). After 28 days of intervention, serum hormone levels were measured by ELISA; gut microbiota composition was analyzed using 16S rRNA gene sequencing, and serum metabolic profiles were characterized by UHPLC-Q-TOF-MS. Spearman correlation analysis with multiple-testing control was performed to explore microbiota-metabolome associations. ZeLan significantly improved PCOS-associated phenotypes, reducing body weight, testosterone (0.92 vs. 1.85 ng/mL, P < 0.05), LH/FSH ratio, and insulin levels compared with the Model group. Sequencing showed that ZeLan partially restored microbial alpha diversity (Shannon index: 5.79 vs. 4.19) and decreased the Firmicutes-to-Bacteroidetes ratio. ZeLan was associated with enrichment of Lactobacillus and Ruminococcus and lower abundance of Bacteroides, Prevotella, and Escherichia. Metabolomics identified 52 altered metabolite features, with key changes involving steroid hormones, bile acids such as lithocholic acid, and short-chain fatty acids such as butyric acid. These findings suggest that ZeLan ameliorates PCOS-related endocrine and metabolic disturbances and is associated with remodeling of gut microbial and serum metabolic profiles.

Introduction

Polycystic ovary syndrome (PCOS) represents one of the most prevalent endocrine disorders affecting women of reproductive age, with a global prevalence estimated between 6% and 20% depending on the diagnostic criteria employed¹. This complex metabolic and reproductive disorder is characterized by a constellation of clinical manifestations, including hyperandrogenism, ovulatory dysfunction, and polycystic ovarian morphology, which collectively contribute to significant reproductive, metabolic, and psychological consequences². Beyond its immediate impact on fertility, PCOS is increasingly recognized as a systemic condition associated with long-term health complications, including type 2 diabetes mellitus, cardiovascular disease, and endometrial carcinoma, thereby imposing substantial burdens on both individual patients and healthcare systems worldwide³. The pathophysiology of PCOS remains incompletely understood, though it is generally acknowledged to involve intricate interactions between genetic predisposition, environmental factors, and lifestyle determinants that culminate in the characteristic hormonal and metabolic disturbances observed in affected individuals⁴.

Contemporary therapeutic approaches for PCOS primarily focus on symptomatic management through pharmacological interventions, including oral contraceptives for menstrual regulation, metformin for insulin sensitization, and anti-androgens for hirsutism control⁵. However, these conventional treatments often provide incomplete symptom relief and are frequently associated with adverse effects that limit their long-term applicability, particularly in women with fertility aspirations⁶. Consequently, there has been growing interest in exploring complementary and alternative medicine approaches, particularly traditional Chinese medicine (TCM), which has accumulated centuries of empirical experience in treating gynecological disorders through holistic regulatory mechanisms⁷. Among the numerous botanical medicines investigated for PCOS management, Lycopus lucidus Turcz. (commonly known as ZeLan in Chinese medicine) has attracted considerable attention due to its traditional applications in promoting blood circulation, resolving stasis, and regulating menstruation⁸. Published phytochemical studies indicate that L. lucidus contains polyphenols and flavonoid-rich fractions with antioxidant and anti-inflammatory activity^8,9. Therefore, standardized total flavonoid content was used as the chemical quality-control index for the extract in the present study, although identification of the specific active monomers responsible for microbiota and metabolic regulation requires targeted phytochemical analysis in the future.

The gut microbiota has emerged as a critical regulator of host metabolism and endocrine function, with mounting evidence supporting its involvement in the pathogenesis of various metabolic disorders, including PCOS¹⁰. The concept of the gut-ovary axis has been proposed to describe bidirectional communication between intestinal microorganisms and ovarian function, mediated by microbial metabolites, immune signaling, and neuroendocrine pathways¹¹. Women with PCOS have been consistently shown to exhibit gut dysbiosis characterized by reduced microbial diversity, altered Firmicutes-to-Bacteroidetes ratios, and diminished abundance of beneficial bacteria such as Lactobacillus and Bifidobacterium species¹². These microbial perturbations are thought to contribute to PCOS pathophysiology through multiple mechanisms, including increased intestinal permeability, systemic low-grade inflammation, and altered bile acid and short-chain fatty acid (SCFA) metabolism, which subsequently influence insulin sensitivity and androgen biosynthesis¹³. Importantly, interventions targeting the gut microbiota, including probiotics, prebiotics, and dietary modifications, have shown promising results in ameliorating PCOS symptoms, thereby supporting the gut-ovary axis as a therapeutic target¹⁴.

Metabolomics, the comprehensive analysis of small-molecule metabolites in biological systems, provides a powerful approach for characterizing the metabolic perturbations associated with disease states and for evaluating the biochemical effects of therapeutic interventions¹⁵. Integration of 16S rRNA gene sequencing for microbiome profiling with untargeted metabolomics enables systematic investigation of relationships between gut microbial communities and host metabolism, offering insights into microbiota-metabolome interactions that may accompany both disease pathogenesis and treatment responses. Despite the theoretical rationale for targeting the gut-ovary axis in PCOS management and ZeLan's traditional indications for gynecological disorders, no studies to date have systematically investigated whether ZeLan extract modulates the gut microbiota-metabolome interface in PCOS. This knowledge gap limits our understanding of ZeLan-associated biological responses and hinders rational optimization of this botanical intervention for PCOS treatment.

The present study was designed to investigate the therapeutic efficacy of ZeLan extract in a DHEA-induced rat model of PCOS and to evaluate associated gut microbiota-metabolome changes through integrated 16S rRNA sequencing and untargeted metabolomics analysis. Specifically, we evaluated the effects of ZeLan extract on PCOS-related phenotypes, including body weight, serum testosterone, luteinizing hormone-to-follicle-stimulating hormone (LH/FSH) ratio, and insulin levels; characterized alterations in gut microbiota composition and alpha diversity following ZeLan treatment; identified differential serum metabolite features and enriched metabolic pathways associated with ZeLan intervention; and explored correlations between specific gut microbial taxa and key metabolites. By addressing these objectives, we sought to provide associative multi-omics evidence for ZeLan-related modulation of the gut-ovary axis in PCOS.

Protocol

All experimental procedures in the present study were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee (Approval No. IACUC-2023-0156).

Animals

Female Sprague-Dawley rats (n = 24, 6 weeks old, body weight 180–200 g) were obtained from the Laboratory Animal Center of the affiliated institution and housed under controlled environmental conditions (temperature 22 ± 2°C, humidity 50 ± 10%, 12-hour light/dark cycle) with ad libitum access to standard rodent chow and purified water. Following a one-week acclimatization period, animals were randomly assigned to three experimental groups (n = 8 per group): Control, PCOS Model, and ZeLan treatment. Sample size was determined based on previous studies investigating gut microbiota alterations in PCOS models, with power analysis indicating that 8 animals per group would provide 80% power to detect a 1.5-fold difference in alpha diversity indices at α = 0.05.

PCOS model establishment and treatment protocol

The PCOS model was established using a well-validated protocol involving subcutaneous injection of dehydroepiandrosterone (DHEA; Sigma-Aldrich, USA) dissolved in sesame oil into the dorsal neck region at a dose of 6 mg/100 g body weight daily for 21 consecutive days. Control animals received equivalent volumes of sesame oil vehicle. The full experimental timeline comprised 7 days of acclimatization, 21 days of model induction, and 28 days of intervention, with endpoint sampling on protocol day 56, including the acclimatization period (Figure 1). The ZeLan extract was prepared from authenticated Lycopus lucidus Turcz. herb through aqueous extraction and was standardized to contain total flavonoids at no less than 2.5% by high-performance liquid chromatography, which served as the chemical quality-control criterion for the botanical extract. Beginning after completion of DHEA-induced model establishment, rats in the ZeLan group received oral gavage of ZeLan extract at 200 mg/kg body weight daily for 28 consecutive days, while the Control and Model groups received equivalent volumes of distilled water. The 200 mg/kg dose was selected as an exploratory pharmacological dose based on prior animal-use experience with botanical extracts and preliminary tolerability considerations; formal dose-response evaluation was not part of the present design. Body weight was recorded weekly throughout the experimental period. In this experimental design, successful PCOS model establishment was evaluated using the DHEA induction protocol together with endpoint endocrine and metabolic phenotypes, including body weight, serum testosterone, LH/FSH ratio, and fasting insulin. Ovarian histopathological examination, follicle counting, corpus luteum quantification, and serial estrous-cycle monitoring were not included as model-confirmation endpoints, and the interpretation of ovarian dysfunction was therefore based on endocrine and metabolic readouts rather than direct tissue-level assessment.

Sample collection and processing

Samples (serum and feces) were collected from all three groups (n = 8 per group) for subsequent biochemical, microbiota, and metabolomics analyses. Endpoint collections were performed in the morning after overnight fasting to reduce circadian and feeding-related variability. Blood samples were collected from the abdominal aorta under anesthesia and centrifuged at 3,000 × g for 15 minutes at 4°C to obtain serum, which was aliquoted and stored at -80°C until hormone and metabolomic analyses. Fresh fecal samples were collected directly from the cecum under aseptic conditions into sterile cryotubes, kept on dry ice during transfer, snap-frozen in liquid nitrogen within 10 minutes of collection, and stored at -80°C until 16S rRNA gene sequencing.

Serum hormone measurements

Serum concentrations of testosterone, luteinizing hormone (LH), follicle-stimulating hormone (FSH), and insulin were quantified using commercially available enzyme-linked immunosorbent assay (ELISA) kits (Cusabio Biotech, Wuhan, China) according to the manufacturer's instructions. The LH/FSH ratio was calculated as a composite index of gonadotropin imbalance characteristic of PCOS. All samples were analyzed in duplicate, and measurements with duplicate coefficients of variation greater than 15% were repeated. The intra-assay and inter-assay coefficients of variation were required to remain below 10% and 15%, respectively. Potential outliers were screened against source records and assay performance information, and no value was excluded solely on statistical grounds without a documented technical reason.

16S rRNA gene sequencing and bioinformatic analysis

Total bacterial genomic DNA was extracted from cecal contents using the QIAamp DNA Stool Mini Kit (Qiagen, Germany) following the manufacturer's protocol. DNA quality and quantity were assessed using NanoDrop spectrophotometry and agarose gel electrophoresis, and DNA extracts were stored at -80°C before library construction. The V3-V4 hypervariable regions of the 16S rRNA gene were amplified using universal primers 338F (5'-ACTCCTACGGGAGGCAGCA-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3') with Illumina adapter sequences. Amplicon libraries were constructed, purified, quantified, pooled in equimolar concentrations, and sequenced on the Illumina NovaSeq 6000 platform (Novogene, Beijing, China) using paired-end 250 bp reads.

Raw sequencing reads were processed using the QIIME2 (version 2022.2) pipeline. Paired-end reads were demultiplexed, merged, quality-filtered using a Phred score threshold of at least 20, and processed with DADA2 to remove chimeric sequences and generate amplicon sequence variants (ASVs). Taxonomic classification was performed against the SILVA 138 reference database using a naive Bayes classifier. Alpha diversity metrics, including Shannon index, Simpson index, and Chao1 estimator, were calculated to assess within-sample diversity and richness. Sequencing-depth outputs, including read summaries, ASV abundance tables, Good's coverage, and rarefaction-curve files, were treated as part of the microbiome dataset described in the Data Availability section. Linear discriminant analysis effect size (LEfSe) was used to identify disease-associated taxa in the Control versus Model comparison using an LDA score threshold greater than 2.0 and P < 0.05; ZeLan-associated reversal was evaluated by comparing the direction of taxonomic abundance shifts in the ZeLan group relative to the Model group.

Untargeted metabolomics analysis

Serum metabolomic profiling was performed using ultra-high-performance liquid chromatography coupled with quadrupole time-of-flight mass spectrometry (UHPLC-Q-TOF-MS; Agilent 1290 Infinity II/6550, USA). Serum samples (100 μL) were mixed with 400 μL of cold methanol: acetonitrile (1:1, v/v) containing internal standards, vortexed, and incubated at -20°C for 1 hour to precipitate proteins. Following centrifugation at 14,000 × g for 15 minutes, the supernatants were dried under a nitrogen stream and reconstituted in 100 μL of acetonitrile:water (1:1, v/v) for analysis. Chromatographic separation was achieved on an ACQUITY UPLC HSS T3 column (2.1 × 100 mm, 1.8 μm) maintained at 40°C, using water containing 0.1% formic acid as mobile phase A and acetonitrile containing 0.1% formic acid as mobile phase B at a flow rate of 0.30 mL/min. A gradient program increasing the organic phase from 5% to 95%, followed by column washing and re-equilibration, was used for metabolite separation. Mass spectrometric data were acquired in both positive and negative electrospray ionization modes with a scan range of m/z 50–1200; the source parameters included capillary voltage 3.5 kV, drying gas flow 10 L/min, drying gas temperature 325°C, nebulizer pressure 35 psi, and fragmentor voltage 120 V.

Raw data files were converted to mzXML format and processed using the XCMS package in R for peak detection, retention-time alignment, and peak-area integration. Quality-control samples prepared by pooling equal aliquots of all serum extracts were inserted throughout the analytical sequence to monitor instrument stability. Features with excessive missingness or unstable quality-control reproducibility were removed before statistical modeling. Metabolite identification was performed by matching accurate mass, retention time, and MS/MS fragmentation patterns against the Human Metabolome Database (HMDB), METLIN, and in-house spectral libraries. Multivariate statistical analyses were conducted using MetaboAnalyst 5.0. Metabolite features retained for downstream reporting were prioritized using a combined criterion that incorporated multivariate discriminatory contribution, fold-change direction, and statistical evidence; P values were adjusted using the Benjamini-Hochberg method to account for multiple testing. The complete annotated reporting table, containing 52 metabolite features, is provided as Supplementary Table 1. Pathway enrichment analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.

Statistical analysis

All statistical analyses were performed using SPSS version 26.0 (IBM Corporation, Armonk, NY, USA) and R software (version 4.2.1). Continuous data were expressed as mean ± standard deviation (SD) and were assessed for normality using the Shapiro-Wilk test before inferential testing. For normally distributed data, one-way analysis of variance (ANOVA) followed by Tukey's post hoc test was used for multiple-group comparisons; for non-normally distributed data, the Kruskal-Wallis test with Dunn's post hoc test was used. Correlations between gut microbial taxa and serum metabolites were evaluated using Spearman's rank correlation coefficient and visualized as hierarchically clustered heatmaps. Correlation P-values were adjusted for multiple testing using the Benjamini-Hochberg false-discovery rate method, and only correlations that were consistent in direction and statistically significant were interpreted. A two-tailed P-value < 0.05 was considered statistically significant unless an adjusted P-value threshold was specified.

Results

Effects of ZeLan extract on PCOS clinical phenotypes

Administration of DHEA induced PCOS-like endocrine and metabolic phenotypes in rats, as evidenced by significant elevations in body weight, serum testosterone concentration, LH/FSH ratio, and fasting insulin levels compared to the Control group (all P < 0.05). Distributional assessment using the Shapiro-Wilk test supported parametric analysis for the principal body weight and hormone endpoints summarized in Table ...

Discussion

The present study provides evidence that ZeLan extract ameliorates PCOS-related endocrine and metabolic phenotypes in a DHEA-induced rat model and that these effects are associated with favorable changes in gut microbiota composition and serum metabolic profiles. The integrated multi-omics approach showed that ZeLan treatment improved hyperandrogenemia, elevated LH/FSH ratio, and insulin resistance, while partially restoring gut microbial alpha diversity, rebalancing microbial composition, and normalizing several metabol...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Acetonitrile	Fisher Chemical / Thermo Fisher Scientific	A955-4	Optima LC-MS grade, 4 L; mobile phase and sample reconstitution
ACQUITY UPLC HSS T3 column	Waters Corporation	186003539	100 Å, 1.8 μm, 2.1 × 100 mm; UHPLC separation
Agarose gel electrophoresis system	Bio-Rad Laboratories	1704486	Mini-Sub Cell GT horizontal electrophoresis system; DNA quality assessment
Agilent 1290 Infinity II/6550 UHPLC-Q-TOF-MS	Agilent Technologies	G7120A/G7167B/G7116B/G6550B	1290 Infinity II UHPLC module set coupled to 6550 iFunnel Q-TOF MS; serum untargeted metabolomics platform
DADA2 algorithm	Bioconductor	Version 1.24.0	ASV generation and chimera removal; version inferred to match R 4.2.1/Bioconductor 3.15 workflow
Dehydroepiandrosterone (DHEA)	Sigma-Aldrich / MilliporeSigma	D4000	PCOS model induction; 6 mg/100 g body weight; CAS 53-43-0
Distilled water	Invitrogen / Thermo Fisher Scientific	10977015	UltraPure DNase/RNase-Free Distilled Water; vehicle for oral gavage in Control and Model groups
Formic acid	Fisher Chemical / Thermo Fisher Scientific	A117-50	Optima LC-MS grade, 50 mL; 0.1% additive for LC-MS mobile phases
Human Metabolome Database	University of Alberta / HMDB	HMDB 5.0	Metabolite annotation database
Illumina NovaSeq 6000 platform	Illumina	20012850	Sequencing system used by sequencing service provider for paired-end 250 bp 16S rRNA sequencing
Lycopus lucidus Turcz. herb	Beijing Tongrentang Chinese Medicine Co., Ltd.	Batch ZL-20230315	Authenticated botanical raw material for ZeLan aqueous extract
Methanol	Fisher Chemical / Thermo Fisher Scientific	A456-4	Optima LC-MS grade, 4 L; protein precipitation for serum metabolomics
METLIN database	Scripps Research	METLIN Classic 2023	Metabolite mass-spectral annotation database
NanoDrop spectrophotometer	Thermo Fisher Scientific	ND-ONE-W	NanoDrop One microvolume UV-Vis spectrophotometer; DNA quantity and purity assessment
QIAamp Fast DNA Stool Mini Kit	Qiagen	51604	DNA extraction from cecal contents; 50-prep stool DNA kit
QIIME2	QIIME2 development team	Version 2022.2	16S rRNA sequence processing
R software	R Foundation	Version 4.2.1	Statistical analysis and XCMS processing
Serum ELISA kits for testosterone, LH, FSH, and insulin	Cusabio Biotech	CSB-E05100r; CSB-E12654r; CSB-E06869r; CSB-E05070r	Rat testosterone, LH, FSH, and insulin ELISA kits for serum hormone and insulin quantification
Sesame oil	Sigma-Aldrich / MilliporeSigma	S3547	Vehicle for DHEA injection
Sodium pentobarbital	Sigma-Aldrich / MilliporeSigma	P3761	Anesthesia; 50 mg/kg intraperitoneal
SPSS software	IBM Corporation	Version 26.0	Statistical analysis
Universal 16S rRNA primers 338F/806R	Sangon Biotech (Shanghai) Co., Ltd.	Custom synthesis: 338F/806R	Amplification of V3-V4 regions; 338F ACTCCTACGGGAGGCAGCAG, 806R GGACTACHVGGGTWTCTAAT
XCMS package	Bioconductor / R	Version 3.18.0	Peak detection, retention-time alignment, and integration
ZeLan extract	In-house aqueous extract	Batch ZLE-20230402	Oral gavage treatment; standardized total flavonoids ≥2.5%