Skip to main content

Integrated fecal microbiota and metabolomics analysis of the orlistat intervention effect on polycystic ovary syndrome rats induced by letrozole combined with a high-fat diet

Abstract

Background

This study aimed to compare the characteristics of the gut microbiota and their metabolite profiles between polycystic ovary syndrome (PCOS) and orlistat-treated PCOS rats (ORL-PCOS), which could help to better understand the underlying mechanism of the effect of orlistat on PCOS.

Methods

PCOS rat models were established using letrozole combined with a high-fat diet. Ten rats were randomly selected as a PCOS control group (PCOS). The other three groups (n = 10/group) were additionally supplemented with different doses of orlistat (low, medium, high). Then, fecal samples of the PCOS and ORL-PCOS groups were analysed by 16S rRNA gene sequencing and untargeted metabolomics. Blood samples were collected to detect serum sex hormones and lipids.

Results

The results showed that orlistat attenuated the body weight gain, decreased the levels of T, LH, the LH/FSH ratio, TC, TG and LDL-C; increased the level of E2; and improved estrous cycle disorder in PCOS rats. The bacterial richness and diversity of the gut microbiota in the ORL-PCOS group were higher than those in the PCOS group. The ratio of Firmicutes to Bacteroidetes was decreased with orlistat treatment. Moreover, orlistat treatment led to a significant decrease in the relative abundance of Ruminococcaceae and Lactobacillaceae, and increases in the abundances of Muribaculaceae and Bacteroidaceae. Metabolic analysis identified 216 differential fecal metabolites in total and 6 enriched KEGG pathways between the two groups, including steroid hormone biosynthesis, neuroactive ligand-receptor interaction and vitamin digestion and absorption. Steroid hormone biosynthesis was the pathway with the most significant enrichment. The correlations between the gut microbiota and differential metabolites were calculated, which may provide a basis for understanding the composition and function of microbial communities.

Conclusions

Our data suggested that orlistat exerts a PCOS treatment effect, which may be mediated by modifying the structure and composition of the gut microbiota, as well as the metabolite profiles of PCOS rats.

Introduction

Polycystic ovary syndrome (PCOS), a complex gynaecological endocrine and metabolic disorder, is characterized by polycystic ovarian morphology (PCOM), infrequent ovulation or anovulation, clinical or biochemical hyperandrogenism, and at least 2 of 3 features are needed to establish a diagnosis [1, 2]. PCOS is also associated with several metabolic disturbances, such as insulin resistance (IR), increased risk of type 2 diabetes mellitus (T2DM) and obesity [1, 3]. The prevalence of PCOS is between 8 and 13% according to the population studied and the definitions used [4].

Studies have found that metabolic factors play important roles in PCOS pathological mechanisms. Some PCOS related genes were found to be related to carbohydrate metabolism and steroid synthesis pathway [5]. Metabolic imbalance could lead to PCOS development [6]. Several studies have shown that PCOS patients have abnormal metabolite composition, including bile acids, short-chain fatty acids and branched-chain amino acids [7, 8]. In addition, gut microbiota (GM) dysbiosis was also found in PCOS development. The GM plays a key role in human and animal health, and maintains a dynamic balance to prevent the development of various diseases. In recent years, evidence has been provided on the correlation between the GM and the development of metabolic diseases, such as obesity and T2DM [9], which accordingly leads to the hypothesis that the GM is closely associated with the aetiology and pathological mechanisms of PCOS [6, 10, 11]. Many studies have investigated this relationship. For example, compared with healthy women, PCOS patients had lower α diversity of the GM, which correlated with the increase in androgens [12]. Sun et al. found that plasma dimethylamine, a product of choline metabolism by the GM, was significantly increased in PCOS patients, demonstrating that GM growth was increased in the group of women with PCOS [13]. In addition, it seems that abnormal short-chain fatty acid (SCFA) metabolism caused by an abnormal GM is associated with IR and hyperandrogenaemia in PCOS patients [14]. These studies have confirmed the close relationship between the GM, metabolites and PCOS.

Over 50% of PCOS cases are overweight [15]. Evidence implies that obesity worsens some features of PCOS, such as infertility, hyperandrogenaemia and IR [16,17,18]. Obesity prevention and treatment will benefit patients with PCOS [18]. Weight loss is the first-line treatment for overweight women with PCOS [19]. Orlistat, an anti-obesity drug approved by the U.S. FDA, has been used to improve or reverse the pathological characteristics of PCOS [20]. Kumar et al. found that orlistat in PCOS is as effective as metformin in reducing weight with improvement in the lipid profile and pregnancy rates [21]. Orlistat could also combine with a low-calorie diet to improve IR and reduce circulating androgens [22]. There is relatively little research on the mechanism by which orlistat improves PCOS. Considering the relationship between the GM, metabolites and PCOS, in the current study, we therefore aimed to compare the characteristics of the gut microbiota and their metabolite profiles between the PCOS and ORL-PCOS groups, which could help to better understand the mechanism of the effect of orlistat on PCOS.

Materials and methods

This study was performed in accordance with the ethical standards laid down in the Declaration of Helsinki, and the Medical Ethics Committee of Lunan Pharmaceutical Group Co., Ltd approved all procedures (Permit No. AN-IACUC-2021–067).

Animals

Specific pathogen-free female rats (5–6 weeks, 145–165 g, Animal Qualification Certificate No.370726211100731686) were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd. (Certificate No. SCXK, 2016–0006; Beijing, China). Before the experiment, all the rats were fed adaptively and quarantined for 1 week. Rats were housed in a breeding room under standard conditions (20–26 °C, 40–70% relative humidity, centralized ventilation of central air conditioner ≥ 15 times per hour, 12-h dark/light cycle). During the experiment, ten rats were randomly selected as a control group (control) and were fed a high-fat diet with 45% of kcal from fat (OpenSource Diets™ Research Diets #D12451). The other rats had free access to water and were fed a high-fat diet (HFD) with 45% of kcal from fat combined with letrozole (1 mg/kg/d. Jiangsu Hengrui Pharmaceutical Co., Ltd; Jiangsu, China) for 28 days to generate a PCOS model [23]. Forty successfully letrozole-induced PCOS rats were randomly selected and divided into four groups (n = 10/group). Ten rats were randomly selected as a PCOS control group (PCOS). The remaining thirty rats were randomly allocated into three groups, which were then supplemented with different doses of orlistat (low-20 mg/kg/d; medium-40 mg/kg/d; high-80 mg/kg/d; Lunan Pharmaceutical Group Corporation; Shandong, Linyi, China.) over the next 12 weeks. The dose of orlistat was selected based on a previous study [24]. Letrozole was given throughout the experiment (12 weeks) to the forty PCOS rats. The control group and PCOS group were given the same volume of solvent without orlistat. Orlistat or solvent was administered by oral gavage. Vaginal epithelial cell smears were taken daily for 10 consecutive days to analyse the estrous-cycle stage. Schematic diagram of the animal experiment design and timeline is listed in Fig. 1.

Fig. 1
figure 1

Schematic diagram of the animal experiment design and timeline

Sample collection

After orlistat treatment, the rats were weighed. At the end of the experiment, Lee’s index was calculated as previously described [25]. Rats in each group fasted for 24 h after the last administration. Blood was taken from the abdominal aorta after sacrificing the rat. Blood samples were collected and then centrifuged at 4,000 rpm for 10 min at 4 °C. The levels of serum sex hormones including testosterone (T), oestradiol (E2), follicle stimulating hormone (FSH) and luteinizing hormone (LH) were tested using radioimmunoassay (Beijing North Institute of Biotechnology Co., Ltd. Beijing, China) according to the specifications for each kit. Biochemical assays were used to analyse lipids, including total cholesterol (TC), triglycerides (TGs), high- density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C). Data were analysed with SPSS 18.0 (IBM SPSS, Armonk, NY, USA) using the Mann–Whitney U test and expressed as the mean ± standard deviation (SD). A p value of < 0.05 was considered statistically significant. Fresh fecal samples (1–3 g) of the PCOS group and the high-doses orlistat treatment PCOS group (ORL-PCOS) were collected in sterile plastic tubes, and then stored at -80 °C for subsequent analyses.

DNA extraction and 16S rRNA gene sequencing

Total bacterial DNA was extracted using the CTAB/SDS method. DNA concentration and purity were monitored on 1% agarose gels. DNA was then diluted to 1 ng/μl using sterile water. 16S rRNA genes were amplified by PCR (98 °C for 1 min, followed by 30 cycles of denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, elongation at 72 °C for 60 s, and 72 °C for 5 min) using the specific primer (16S V3-V4: 341F-806R). All PCRs were carried out in 30 μL reactions with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs); 0.2 μM of forward and reverse primers, and approximately 10 ng of template DNA. PCR products were mixed with the same volume of 1X loading buffer containing SYBR Green), and electrophoresis was performed on a 2% agarose gel for detection. Samples with bright main bands between 400–450 bp were chosen for further experiments. PCR products were mixed in equidensity ratios. Then, the mixed PCR products were purified with an AxyPrepDNA Gel Extraction Kit (AXYGEN). Sequencing libraries were generated using the NEB Next®Ultra™DNA Library Prep Kit for Illumina (NEB, USA) following the manufacturer’s recommendations. The library was sequenced on an Illumina MiSeq/HiSeq2500 platform and 250 bp/300 bp paired-end reads were generated.

Analysis of the gut microbiota

Paired-end reads from the original DNA fragments were merged using FLASH and assigned to each sample according to the unique barcodes. Sequence analysis was performed by the UPARSE software package using the UPARSE-OTU and UPARSE-OTUref algorithms. Sequences with ≥ 97% similarity were assigned to the same OTUs. In-house Perl scripts were used to analyse alpha diversity. Two metrics were calculated: the Chao1 and Shannon indices. For beta diversity analysis, principal coordinate analysis (PCoA) between the two groups based on both weighted and unweighted UniFrac algorithms was calculated using QIIME. The linear discriminant analysis (LDA)-effect size (LEfSe) method was used for the quantitative analysis of biomarkers within different groups. An LDA threshold > 3 was used as a threshold to identify the most differentially abundant taxa. Correlations between fecal microbiota and metabolites were analysed with R language (R 3.4.2).

Metabolite extraction and LC–MS/MS analysis

Metabolites in the stool samples (60 mg) were extracted as previously described [26]. Briefly, sample metabolites were extracted using 400 μL of a methanol: water (4:1, v/v) solution. Then, the mixture was treated by a high-throughput tissue crusher at -20 °C, and 50 Hz for 6 min, vortex mixed for 30 s, ultrasonicated at 40 kHz for 30 min at 5 °C, and placed at -20 °C for 30 min to precipitate proteins. The mixture was centrifuged for 15 min (13,000 g, 4 °C). The supernatant was dried in a vacuum centrifuge. For LC–MS analysis, the samples were redissolved in 100 μL of acetonitrile/water (1:1, v/v) solvent. Analysis was performed using an UHPLC (1290 Infinity LC, Agilent Technologies) coupled to a quadrupole time-of-flight (AB Sciex TripleTOF 6600) at Shanghai Applied Protein Technology Co., Ltd. For HILIC separation, samples were analysed using a 2.1 mm × 100 mm ACQUITY UPLC BEH 1.7 µm column (Waters, Ireland). In both ESI positive and negative modes, the mobile phase contained A = 25 mM ammonium acetate and 25 mM ammonium hydroxide in water and B = acetonitrile. The gradient was 85% B for 1 min, linearly reduced to 65% in 11 min, reduced to 40% in 0.1 min, maintained for 4 min, and then increased to 85% in 0.1 min, with a 5 min re-equilibration period employed.

Metabolomics analysis

The raw MS data were converted to MzXML files using ProteoWizard MSConvert before importing into freely available XCMS software. After normalization to the total peak intensity, the processed data were analysed by the R package, where they were subjected to multivariate data analysis, including orthogonal partial least-squares discriminant analysis (OPLS-DA). The sevenfold cross-validation and response permutation testing were used to evaluate the robustness of the model. The variable importance in the projection (VIP) value of each variable in the OPLS-DA model was calculated to indicate its contribution to the classification. Metabolites with VIP values > 1 were further subjected to Student’s t-test at the univariate level to measure the significance of each metabolite, and a p value of < 0.05 was considered statistically significant. For pathway enrichment analysis, the differential metabolites were annotated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.kegg.jp/kegg/pathway.html).

Statistical analysis

Values are expressed as the mean ± standard deviation (SD). All data were analyzed using SPSS 18.0 (International Business Machines Corporation, Armonk, New York, USA). The body weight, Lee’s index, serum sex hormones, and lipid profiles were analyzed by one-factor analysis of variance (ANOVA), followed by post-hoc the least significance difference (LSD)-t test. p < 0.05 was deemed statistically significant in all tests.

Results

Effects of orlistat on weight, sex hormones levels, lipids and the estrous cycle

Compared with the control (CON) group, the weight of the PCOS rats was significantly higher from week 4 to week 16. Meanwhile, the Lee’s index of the PCOS group was also significantly increased (Table 1). The weight and Lee’s indices of the rats in all the orlistat treatment groups were significantly decreased (p < 0.05) compared with those of the PCOS group.

Table 1 Effects of Orlistat on body weight and Lee’s index in PCOS rats

Compared with those in the CON group, the levels of serum T, and LH, and the ratio of LH/FSH were significantly higher (p < 0.05), and E2 was significantly lower (p < 0.01) in the PCOS group (Table 2). The rats in the high-orlistat group had significantly lower levels of T, LH, and the ratio of LH/FSH, and higher levels of E2 compared with those of the PCOS group. In addition, compared with the CON group, the PCOS rats had significantly higher concentrations of TC, TG and LDL-C (Table 2). The rats in the high- and medium-orlistat group had significantly lower levels of TC, TG and LDL-C than those in the PCOS group. No differences in HDL-C levels were found.

Table 2 Effect of supplementation with orlistat on serum sex hormone and lipids levels in PCOS rats

Estrous cycle assessment revealed that PCOS rats experienced a longer dioestrus and shorter proestrus and estrous than control rats (Fig. 2). Following orlistat treatment, the estrous cycle disorder was restored. The improvement effect on the estrous cycle disorder of the high-orlistat group was better than that of the low-orlistat group and medium-orlistat group.

Fig. 2
figure 2

Effects of orlistat treatment on the estrous cycle. Three rats were randomly selected in each group and representative estrous cycles of the different groups are shown. P: Proestrus, E: Estrous, M: Metoestrus, D: Dioestrus

Effects of orlistat on the gut microbiota in PCOS mice

A total of 3,208 OTUs were obtained in this study (Fig. 3A). Then the Shannon index and Chao1 index (Fig. 3B and C) were calculated. No significant differences were observed between the diversity of the gut microbiota in the ORL-PCOS and PCOS groups. PCoA between groups based on the weighted UniFrac (Fig. 3D) algorithms was performed and showed a separation between the samples in these two groups. LEfSe (Fig. 4A and B) was used to analyse biomarkers in the microbiota of each group. The PCOS group was significantly enriched in Lactobacillus, Ruminococcaceae_UCG_005 and Ruminococcaceae_UCG_014. The ORL-PCOS group was significantly enriched in Bacteroides and Muribaculaceae.

Fig. 3
figure 3

Characteristics of the gut microbiota after orlistat treatment. A Venn diagram of OTUs. The overlapping section represents the shared OTUs. There was no significant difference in the Shannon index. B and Chao1 (C) index between the ORL-PCOS and PCOS groups. D PCoA of the gut microbiota of ORL-PCOS and PCOS group based on weighted UniFrac algorithms

Fig. 4
figure 4

LEfSe of the gut microbiota of the ORL-PCOS and PCOS groups. A LEfSe taxonomic cladogram. B LDA chart (p = phylum, c = class, o = order, f = family, g = genus). LDA scores higher than 3.0 were used as a threshold for significance in LEfSe analyses

Effects of orlistat on fecal metabolomics in PCOS mice

OPLS-DA was performed to verify differential the metabolites between the two groups (Fig. 5A and C), and was validated by permutation analysis (Fig. 5B and D). The results showed clear separation between the ORL-PCOS and PCOS groups, suggesting a significant difference in fecal metabolites. In this study, VIP > 1.0 was used to screen metabolite differences, and only metabolites identified with a fold change (FC) > 1.5 or < 1/1.5, p < 0.05 and VIP > 1.0 were considered differential. Based on the FC and p values, volcanic diagrams were made to illustrate the up- and down-regulated metabolites between the two groups (Fig. 6A and B). All detailed information on the different metabolites is listed in Supplementary Table 1. In addition, the differential metabolites were enriched and analysed based on the KEGG database (Fig. 7A). Most of the differential metabolites were involved in steroid hormone biosynthesis (Fig. 7B), vitamin digestion and absorption (Fig. 7C) and neuroactive ligand-receptor interaction (Fig. 7D).

Fig. 5
figure 5

Multivariate statistical analyses of fecal metabolites in the ORL-PCOS and PCOS groups. A OPLS-DA negative ion mode scatter diagram. B Corresponding validation plots of the OPLS-DA negative ion mode obtained from the permutation test. C OPLS-DA positive ion mode scatter diagram. D Corresponding validation plots of the OPLS-DA positive ion mode obtained from the permutation test

Fig. 6
figure 6

Volcanic map analysis in negative (A) and positive (B) ion modes. According to the VIP value, p value and FC value, volcano maps were drawn to display the overall distribution of the differential metabolites. The blue dots represent significantly downregulated metabolites in the ORL-PCOS group, and the red dots represent significantly upregulated differential metabolites. The black dots represent nonsignificant differential metabolites

Fig. 7
figure 7

KEGG pathway enrichment analysis of differential metabolites. A The ordinate represents the name of the pathway, and the abscissa represents the number of differentially expressed metabolites included in each KEGG metabolic pathway. The colour indicates the p value of enrichment analysis. The number on the column represents the rich factor. B Heatmap for steroid hormone biosynthesis, C Heatmap for vitamin digestion and absorption, D Heatmap for neuroactive ligand-receptor interaction

Correlation analysis of the fecal microbiota with metabolites in orlistat treated PCOS mice

The correlation between the differential metabolites and the microbial diversity in orlistat treated PCOS mice was analysed to further explore the composition and function of microorganisms. As shown in Fig. 8, we found that steroid hormone biosynthesis pathway related differential metabolites, such as 3-dehydroepiandrosterone sulfate (DHEAS), and 5alpha-androstan-17-beta-ol-3-one (dihydrotestosterone; DHT), were negatively correlated with the abundance of Bacteroides, and positively correlated with the abundances of Ruminococcaceae_UCG_014, Ruminococcaceae_UCG_005 and Pseudomonas. Progesterone was positively correlated with the abundance of Bacteroides, and Desulfovibrio and negatively correlated with the abundance of Ruminococcaceae_UCG_014 and Ruminococcaceae_UCG_005 by correlation analysis.

Fig. 8
figure 8

Correlation heatmap analysis of the differential fecal metabolites and gut microbiota between the PCOS and ORL-PCOS groups. A significant positive correlation is shown in red; a significant negative correlation is shown in blue, with darker colours indicating stronger correlations. * p < 0.05, ** p < 0.01, *** p < 0.001

Discussion

Orlistat exerts an anti-obesity effect by inhibiting the absorption of triglycerides. Multiple clinical studies have proven its effect on weight loss [27, 28]. A recent study found that orlistat treatment has beneficial effects on body weight in HFD-induced obese mice by modifying the composition of the gut microbiota, presenting as reduced total microbial abundance and obviously upregulated bacteria [29]. In addition, 5% weight loss of initial body weight has been shown to increase ovulation frequency and fertility, and to improve testosterone and lipid levels in PCOS patients [30,31,32]. In this study, the gut microbiota and its metabolite profiles were first compared between PCOS mice and orlistat-treated PCOS mice, highlighting the critical role of microbiota and microbiota metabolites in controlling PCOS development, which could help to better understand the interaction between orlistat and PCOS.

Studies have shown that the gut microbiota is related to obesity or PCOS, and several phyla of microbiota may be involved in the occurrence and development of obesity or PCOS [14, 29, 33]. Firmicutes and Bacteroidetes are the two most abundant bacterial phyla in the human gut [34]. It has been found that the fecal community of obese mice/humans or PCOS patients/mouse models is characterized by an increased ratio of Firmicutes to Bacteroidetes (F/B) [35,36,37]. Changes in the F/B ratio were found during PCOS treatment. For example, Lin et al. observed a reduced abundance of Firmicutes and an elevated abundance of Bacteroidetes in response to sleeve gastrectomy in a dehydroepiandrosterone-induced PCOS rat model [38]. Additionally, flaxseed oil intervention modulated the gut microbiota and ameliorated PCOS in rats by decreasing the ratio of F/B [39]. A recent study found that orlistat could further decrease the F/B ratio in HFD-induced obese mice, exerting beneficial effects on body weight [29]. Consistently, in this study, we also found that the ORL-PCOS had significant changes in the composition of the gut microbiota, such as an increase in the abundance of Bacteroidetes and a decrease in Firmicutes. Bacteroidetes exert immunomodulatory effects on the host, and its elevation has been shown to be associated with weight loss in humans and animal models [40, 41]. Additionally, the Bacteroides genus belonging to the Bacteroidetes phylum was negatively associated with T2DM [42]. Firmicutes, which is closely related to obesity and metabolic syndrome, play a role in energy resorption and facilitate fat storage in the host body [11]. A previous study showed that orlistat is a reversible gastric and pancreatic lipase inhibitor and could promote weight loss by partially preventing intestinal fat absorption [43]. Weight loss is associated with changes in the gut microbiota [44]. Therefore, we speculated that orlistat may exert beneficial effects on body weight by modifying the gut microbiota, such as the decreased abundance of Firmicutes, which may help reduce fat storage and weight loss.

At the family level, we observed that orlistat treatment led to a significant decrease in the relative abundance of Ruminococcaceae and Lactobacillaceae, and an elevation of Muribaculaceae and Bacteroidaceae. Previous studies have found that the relative abundance of Ruminococcaceae was higher in obese patients with PCOS and letrozole-induced PCOS rat models [37, 45, 46]. Furthermore, the Ruminococcaceae abundance was reported to be associated with diabetes and testosterone levels. Mokkala et al. revealed that the relative abundance of Ruminococcaceae increased twofold in women with gestational diabetes compared to healthy controls [47]. The Ruminococcus genus can increase inflammatory cytokine production, which is positively associated with T2DM [42]. Eyupoglu et al. reported that Ruminococcaceae appears to be associated with clinical androgen excess in patients with PCOS [45]. Lactobacillaceae can ferment carbohydrates, including pectin and glucose, to produce formic acid, lactic acid, acetic acid and ethanol. Although Lactobacillaceae is commonly used as a probiotic, studies have shown that Lactobacillaceae is positively correlated with weight gain, and a higher level of Lactobacillaceae was detected in obese children and adults than in lean control individuals [48, 49], suggesting that the effect of Lactobacillaceae on metabolism may be species and strain specific [50]. Moreover, a significantly increased relative abundance of the Lactobacillus genus was observed in PCOS with insulin resistance rat models [46]. A lower level of Lactobacillaceae could alleviate obesity in HFD-induced obese mice. Ye et al. suggested that ripened pu-erh tea serves as a great candidate to alleviate obesity by increasing the abundance of Bacteroidaceae, Muribaculaceae and decreasing the abundance of Lactobacillaceae in HFD-induced obese mice [51], which is consistent with the results of the present study. Bacteroidaceae could promote cytokine production and have a potential role in inhibiting autoimmune disease [51]. Obese adolescents with PCOS had a lower relative abundance of the family Bacteroidaceae, which conferred a 4.4-fold higher likehood ratio of taxa predictive of PCOS diagnosis [52]. Muribaculaceae could degrade dietary components and polysaccharides to produce short chain fatty acids (SCFAs) that play key roles in anti-inflammatory and glycolipid homeostasis balance [53]. PCOS has been proven to be a chronic inflammatory disease, and obesity, insulin resistance, T2DM and hyperandrogenemia are common symptoms [26]. In the present study, we speculated that orlistat may improve the above PCOS symptoms by changing the gut microbiota composition, indicating that these bacteria might be the most efficient taxa contributing to preventing the development of PCOS.

Lipid metabolism is dysregulated in women with PCOS. Lipids are involved in various metabolic pathways, such as steroid hormone biosynthesis, and fatty acid metabolism [54]. In this study, based on the KEGG enrichment analysis, steroid hormone biosynthesis was the pathway in which the differential metabolites were most significantly enriched, such as DHEAS, dihydrotestosterone (DHT), progesterone and tetrahydrocorticosterone. The levels of DHEAS and DHT were significantly downregulated 0.40-fold and 0.16-fold respectively in orlistat treated PCOS mice compared with the control group. Correlation analysis revealed that both were negatively correlated with the abundance of Bacteroides, and positively correlated with the abundance of Ruminococcaceae_UCG_014, Ruminococcaceae_UCG_005 and Pseudomonas. In a study on the fecal metabolites and gut microbiota in obese patients with PCOS, several fecal metabolites were used as characteristic metabolites, including DHEAS, which was significantly and positively correlated with the serum testosterone level and negatively correlated with body mass index (BMI) or fasting insulin among PCOS patients [26, 55]. DHEAS is a kind of androgen mainly secreted by the adrenal cortex, and found to be high in 22–25% of patients with PCOS [56]. Higher DHEAS could increase incidence of degenerated oocytes and early miscarriage rates in women with PCOS [57]. In addition, the prevalence of adrenal hyperandrogenaemia, which is defined as elevated circulating DHEAS levels, in women with PCOS is 15% to 45% [58]. Similar to DHEAS, the concentration of DHT, a potent androgen, was significantly higher in PCOS patients [59]. The DHT-induced PCOS model is one of the three typically used androgen-based PCOS models [60]. DHT exposure could cause the key reproductive characteristics (aperiodic, anovulatory and multifollicular ovaries) and metabolic characteristics (increased body weight and visceral adiposity) of PCOS [61]. Moreover, in this study, the level of progesterone was significantly upregulated by 6.74-fold in orlistat treated PCOS mice. Progesterone was positively correlated with the abundance of Bacteroides, and Desulfovibrio and negatively correlated with the abundance of Ruminococcaceae_UCG_014 and Ruminococcaceae_UCG_005 by correlation analysis. Progesterone is a key steroid hormone and is absolutely required for uterine implantation, decidualization, and maintenance of pregnancy [62]. PCOS women have lower progesterone levels during the luteal phase which overstimulates the immune system that produces more estrogen [63]. The chronic anovulation seen in PCOS implies long-term oestrogen excess or a lack of progesterone [64]. PCOS is a state of altered steroid hormone production and activity, and there is a close relationship between gut microbiota and sex hormones in PCOS [65]. Therefore, we speculated that orlistat may modify the gut microbiota of obese rats with PCOS, which has an impact on sex hormones, including decreased levels of DHEAS, and DHT and increased progesterone levels, to improve the steroid hormone state in PCOS.

Conclusions

There are several limitations in this study. The number of samples was relatively small, a negative control group in the gut microbiome analysis was lacking, and the correlations were not verified. It is necessary to expand the sample size and verification in further research. This study was the first to use fecal microbiota combined with its metabolite profiles to explore the differences between PCOS mice and orlistat-treated PCOS mice, which may offer more specific mechanisms of orlistat treatment on PCOS and its related pathways. Based on 16S rRNA gene sequencing and untargeted metabolomics analysis, orlistat intervention modified the structure and composition of the gut microbiota, as well as the metabolite profiles of PCOS mice. Decreased abundance of Firmicutes and increased Bacteroidetes were observed in orlistat treated PCOS mice. Some families, such as Ruminococcaceae and Lactobacillaceae, which are associated with diabetes and testosterone levels or obesity, were decreased by orlistat intervention. Muribaculaceae and Bacteroidaceae which play anti-inflammatory roles, were notably enhanced by orlistat. The gut microbiota, fecal metabolites, and hormones were closely correlated. Steroid hormone biosynthesis was the pathway in which the differential metabolites were most significantly enriched.

Availability of data and materials

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI SRA (accession: SUB12040976).

References

  1. Teede HJ, Misso ML, Costello MF, et al. Recommendations from the international evidence-based guideline for the assessment and management of polycystic ovary syndrome. Hum Reprod. 2018;33:1602–18.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Rotterdam EA-SPcwg. Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome (PCOS). Hum Reprod. 2004;19:41–7.

    Article  Google Scholar 

  3. Li R, Zhang Q, Yang D, et al. Prevalence of polycystic ovary syndrome in women in China: a large community-based study. Hum Reprod. 2013;28:2562–9.

    Article  PubMed  Google Scholar 

  4. Azziz R, Carmina E, Chen Z, et al. Polycystic ovary syndrome. Nat Rev Dis Primers. 2016;2:16057.

    Article  PubMed  Google Scholar 

  5. Li T, Wu K, You L, et al. Common variant rs9939609 in gene FTO confers risk to polycystic ovary syndrome. PLoS one. 2013;8:e66250.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Giampaolino P, Foreste V, Di Filippo C, et al. Microbiome and PCOS: State-of-Art and Future Aspects. Int J Mol Sci. 2021;22:2048.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Qi X, Yun C, Sun L, et al. Gut microbiota-bile acid-interleukin-22 axis orchestrates polycystic ovary syndrome. Nat Med. 2019;25:1225–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Xu WL, Liu GY, Zhang N, et al. Untargeted metabolomics analysis of serum and follicular fluid samples from women with polycystic ovary syndrome. Minerva Endocrinol. 2023;48:160–71.

    Article  Google Scholar 

  9. Barko PC, McMichael MA, Swanson KS, et al. The gastrointestinal microbiome: a review. J Vet Intern Med. 2018;32:9–25.

    Article  CAS  PubMed  Google Scholar 

  10. Cani PD, Amar J, Iglesias MA, et al. Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes. 2007;56:1761–72.

    Article  CAS  PubMed  Google Scholar 

  11. Turnbaugh PJ, Ley RE, Mahowald MA, et al. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444:1027–31.

    Article  PubMed  Google Scholar 

  12. Torres PJ, Siakowska M, Banaszewska B, et al. Gut microbial diversity in women with polycystic ovary syndrome correlates with hyperandrogenism. J Clin Endocrinol Metab. 2018;103:1502–11.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Sun L, Hu W, Liu Q, et al. Metabonomics reveals plasma metabolic changes and inflammatory marker in polycystic ovary syndrome patients. J Proteome Res. 2012;11:2937–46.

    Article  CAS  PubMed  Google Scholar 

  14. Zhao X, Jiang Y, Xi H, et al. Exploration of the Relationship Between Gut Microbiota and Polycystic Ovary Syndrome (PCOS): a Review. Geburtshilfe Frauenheilkd. 2020;80:161–71.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Shi Y, Guo M, Yan J, et al. Analysis of clinical characteristics in large-scale Chinese women with polycystic ovary syndrome. Neuro Endocrinol Lett. 2007;28:807–10.

    CAS  PubMed  Google Scholar 

  16. Dahan MH, Reaven G. Relationship among obesity, insulin resistance, and hyperinsulinemia in the polycystic ovary syndrome. Endocrine. 2019;64:685–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Brassard M, AinMelk Y, Baillargeon JP. Basic infertility including polycystic ovary syndrome. Med Clin North Am. 2008;92(1163–1192):xi.

    PubMed  Google Scholar 

  18. Lim SS, Norman RJ, Davies MJ, et al. The effect of obesity on polycystic ovary syndrome: a systematic review and meta-analysis. Obes Rev. 2013;14:95–109.

    Article  CAS  PubMed  Google Scholar 

  19. Moran LJ, Pasquali R, Teede HJ, et al. Treatment of obesity in polycystic ovary syndrome: a position statement of the Androgen Excess and Polycystic Ovary Syndrome Society. Fertil Steril. 2009;92:1966–82.

    Article  PubMed  Google Scholar 

  20. Li Y, Chen C, Ma Y, et al. Multi-system reproductive metabolic disorder: significance for the pathogenesis and therapy of polycystic ovary syndrome (PCOS). Life Sci. 2019;228:167–75.

    Article  CAS  PubMed  Google Scholar 

  21. Kumar P, Arora S. Orlistat in polycystic ovarian syndrome reduces weight with improvement in lipid profile and pregnancy rates. J Hum Reprod Sci. 2014;7:255–61.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Panidis D, Farmakiotis D, Rousso D, et al. Obesity, weight loss, and the polycystic ovary syndrome: effect of treatment with diet and orlistat for 24 weeks on insulin resistance and androgen levels. Fertil Steril. 2008;89:899–906.

    Article  PubMed  Google Scholar 

  23. Wu YX, Yang XY, Han BS, et al. Naringenin regulates gut microbiota and SIRT1/ PGC-1a signaling pathway in rats with letrozole-induced polycystic ovary syndrome. Biomed Pharmacother. 2022;153: 113286.

    Article  CAS  PubMed  Google Scholar 

  24. Lucas N, Legrand R, Deroissart C, et al. Hafnia alvei HA4597 strain reduces food intake and body weight gain and improves body composition, glucose, and lipid metabolism in a mouse model of hyperphagic obesity. Microorganisms. 2019;8:35.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Peng MF, Tian S, Song YG, et al. Effects of total flavonoids from Eucommia ulmoides Oliv. leaves on polycystic ovary syndrome with insulin resistance model rats induced by letrozole combined with a high-fat diet. J Ethnopharmacol. 2021;273:113947.

    Article  CAS  PubMed  Google Scholar 

  26. Zhou L, Ni Z, Yu J, et al. Correlation between fecal metabolomics and gut microbiota in obesity and polycystic ovary syndrome. Front Endocrinol (Lausanne). 2020;11:628.

    Article  PubMed  Google Scholar 

  27. Sahebkar A, Simental-Mendia LE, Reiner Z, et al. Effect of orlistat on plasma lipids and body weight: a systematic review and meta-analysis of 33 randomized controlled trials. Pharmacol Res. 2017;122:53–65.

    Article  CAS  PubMed  Google Scholar 

  28. Aldekhail NM, Logue J, McLoone P, et al. Effect of orlistat on glycaemic control in overweight and obese patients with type 2 diabetes mellitus: a systematic review and meta-analysis of randomized controlled trials. Obes Rev. 2015;16:1071–80.

    Article  CAS  PubMed  Google Scholar 

  29. Ke J, An Y, Cao B, et al. Orlistat-induced gut microbiota modification in obese mice. Evid Based Complement Alternat Med. 2020;2020:9818349.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Legro RS. Obesity and PCOS: implications for diagnosis and treatment. Semin Reprod Med. 2012;30:496–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Pasquali R, Gambineri A, Pagotto U. The impact of obesity on reproduction in women with polycystic ovary syndrome. BJOG. 2006;113:1148–59.

    Article  CAS  PubMed  Google Scholar 

  32. Motta AB. The role of obesity in the development of polycystic ovary syndrome. Curr Pharm Des. 2012;18:2482–91.

    Article  CAS  PubMed  Google Scholar 

  33. Wang P, Gao J, Ke W, et al. Resveratrol reduces obesity in high-fat diet-fed mice via modulating the composition and metabolic function of the gut microbiota. Free Radic Biol Med. 2020;156:83–98.

    Article  CAS  PubMed  Google Scholar 

  34. John GK, Mullin GE. The gut microbiome and obesity. Curr Oncol Rep. 2016;18:45.

    Article  PubMed  Google Scholar 

  35. Gomes AC, Hoffmann C, Mota JF. The human gut microbiota: Metabolism and perspective in obesity. Gut Microbes. 2018;9:308–25.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Lindheim L, Bashir M, Munzker J, et al. Alterations in gut microbiome composition and barrier function are associated with reproductive and metabolic defects in women with Polycystic Ovary Syndrome (PCOS): a pilot study. PLoS one. 2017;12:e0168390.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Kelley ST, Skarra DV, Rivera AJ, et al. The gut microbiome is altered in a letrozole-induced mouse model of polycystic ovary syndrome. PLoS one. 2016;11:e0146509.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Lin W, Wen L, Wen J, et al. Effects of sleeve gastrectomy on fecal gut microbiota and short-chain fatty acid content in a rat model of polycystic ovary syndrome. Front Endocrinol (Lausanne). 2021;12:747888.

    Article  PubMed  Google Scholar 

  39. Wang T, Sha L, Li Y, et al. Dietary alpha-Linolenic acid-rich flaxseed oil exerts beneficial effects on polycystic ovary syndrome through sex steroid hormones-microbiota-inflammation axis in rats. Front Endocrinol (Lausanne). 2020;11:284.

    Article  PubMed  Google Scholar 

  40. Liu R, Hong J, Xu X, et al. Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat Med. 2017;23:859–68.

    Article  CAS  PubMed  Google Scholar 

  41. Turnbaugh PJ, Hamady M, Yatsunenko T, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.

    Article  CAS  PubMed  Google Scholar 

  42. Gurung M, Li Z, You H, et al. Role of gut microbiota in type 2 diabetes pathophysiology. EBioMedicine. 2020;51:102590.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Padwal RS, Majumdar SR. Drug treatments for obesity: orlistat, sibutramine, and rimonabant. Lancet. 2007;369:71–7.

    Article  CAS  PubMed  Google Scholar 

  44. Koutoukidis DA, Jebb SA, Zimmerman M, et al. The association of weight loss with changes in the gut microbiota diversity, composition, and intestinal permeability: a systematic review and meta-analysis. Gut Microbes. 2022;14:2020068.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Eyupoglu ND, Ergunay K, Acikgoz A, et al. Gut microbiota and oral contraceptive use in overweight and obese patients with polycystic ovary syndrome. J Clin Endocrinol Metab. 2020;105:dgaa600.

    Article  PubMed  Google Scholar 

  46. Zhu Y, Li Y, Liu M, et al. Guizhi Fuling Wan, Chinese Herbal Medicine, Ameliorates Insulin Sensitivity in PCOS Model Rats With Insulin Resistance via Remodeling Intestinal Homeostasis. Front Endocrinol (Lausanne). 2020;11:575.

    Article  PubMed  Google Scholar 

  47. Mokkala K, Houttu N, Vahlberg T, et al. Gut microbiota aberrations precede diagnosis of gestational diabetes mellitus. Acta Diabetol. 2017;54:1147–9.

    Article  PubMed  Google Scholar 

  48. Jin J, Cheng R, Ren Y, et al. Distinctive gut microbiota in patients with overweight and obesity with dyslipidemia and its responses to long-term orlistat and ezetimibe intervention: a randomized controlled open-label trial. Front Pharmacol. 2021;12:732541.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Bervoets L, Van Hoorenbeeck K, Kortleven I, et al. Differences in gut microbiota composition between obese and lean children: a cross-sectional study. Gut Pathog. 2013;5:10.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Million M, Angelakis E, Paul M, et al. Comparative meta-analysis of the effect of Lactobacillus species on weight gain in humans and animals. Microb Pathog. 2012;53:100–8.

    Article  PubMed  Google Scholar 

  51. Ye J, Zhao Y, Chen X, et al. Pu-erh tea ameliorates obesity and modulates gut microbiota in high fat diet fed mice. Food Res Int. 2021;144:110360.

    Article  CAS  PubMed  Google Scholar 

  52. Jobira B, Frank DN, Pyle L, et al. Obese adolescents with PCOS have altered biodiversity and relative abundance in gastrointestinal microbiota. J Clin Endocrinol Metab. 2020;105:e2134–44.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Li LL, Wang YT, Zhu LM, et al. Inulin with different degrees of polymerization protects against diet-induced endotoxemia and inflammation in association with gut microbiota regulation in mice. Sci Rep. 2020;10:978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Rajska A, Buszewska-Forajta M, Rachon D, et al. Metabolomic insight into polycystic ovary syndrome-an overview. Int J Mol Sci. 2020;21:4853.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Moran C, Arriaga M, Arechavaleta-Velasco F, et al. Adrenal androgen excess and body mass index in polycystic ovary syndrome. J Clin Endocrinol Metab. 2015;100:942–50.

    Article  PubMed  Google Scholar 

  56. Khan SH, Rizvi SA, Shahid R, et al. Dehydroepiandrosterone Sulfate (DHEAS) Levels in Polycystic Ovarian Syndrome (PCOS). J Coll Physicians Surg Pak. 2021;31:253–7.

    Article  PubMed  Google Scholar 

  57. Jimenez PT, Frolova AI, Chi MM, et al. DHEA-mediated inhibition of the pentose phosphate pathway alters oocyte lipid metabolism in mice. Endocrinology. 2013;154:4835–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Luque-Ramirez M, Escobar-Morreale HF. Adrenal hyperandrogenism and polycystic ovary syndrome. Curr Pharm Des. 2016;22:5588–602.

    Article  CAS  PubMed  Google Scholar 

  59. O’Reilly MW, Kempegowda P, Walsh M, et al. AKR1C3-mediated adipose androgen generation drives lipotoxicity in women with polycystic ovary syndrome. J Clin Endocrinol Metab. 2017;102:3327–39.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Walters KA, Allan CM, Handelsman DJ. Rodent models for human polycystic ovary syndrome. Biol Reprod. 2012;86(149):141–112.

    Google Scholar 

  61. Rodriguez Paris V, Edwards MC, Aflatounian A, et al. Pathogenesis of reproductive and metabolic PCOS traits in a mouse model. J Endocr Soc. 2021;5:bvab060.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Wetendorf M, DeMayo FJ. The progesterone receptor regulates implantation, decidualization, and glandular development via a complex paracrine signaling network. Mol Cell Endocrinol. 2012;357:108–18.

    Article  CAS  PubMed  Google Scholar 

  63. Luan YY, Zhang L, Peng YQ, et al. Immune regulation in polycystic ovary syndrome. Clin Chim Acta. 2022;531:265–72.

    Article  CAS  PubMed  Google Scholar 

  64. Li X, Feng Y, Lin JF, et al. Endometrial progesterone resistance and PCOS. J Biomed Sci. 2014;21:2.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Insenser M, Murri M, Del Campo R, et al. Gut microbiota and the polycystic ovary syndrome: influence of sex, sex hormones, and obesity. J Clin Endocrinol Metab. 2018;103:2552–62.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Prof Chen Chen (Endocrinology, SBMS, Faculty of Medicine, The University of Queensland, St Lucia, Qld 4072, Australia) for his assistance in discussing and editing the final version of the manuscript.

Funding

This work was financially supported by National Natural Science Foundation of China (82100920), Natural Science Foundation of Shandong Province (CN) (ZR2020MH107), China Postdoctoral Science Foundation (2021M691387) and Postdoctoral Creative Funding in Shangdong Province (202103062).

Author information

Authors and Affiliations

Authors

Contributions

Jingchun Yao conceived and designed the study. Enli Wang and Mingmin Jiang raised the animals. Jianmei Yang, Enli Wang, Yujun Tan, Fangfang Yao, Ling Gao, Chenghong Sun and Lihong Pan conducted the research. Jianmei Yang designed the statistical analyses and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jingchun Yao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, J., Wang, E., Jiang, M. et al. Integrated fecal microbiota and metabolomics analysis of the orlistat intervention effect on polycystic ovary syndrome rats induced by letrozole combined with a high-fat diet. J Ovarian Res 16, 109 (2023). https://doi.org/10.1186/s13048-023-01193-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13048-023-01193-3

Keywords