Dysbiosis of Gut Microbiota Associated with Clinical Parameters in Polycystic Ovary Syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine and metabolic disorder in women. Gut microbiota has been implicated to play a critical role in metabolic diseases and may modulate the secretion of mediators of the brain–gut axis. Interaction between gut microbiota and the endocrine and biochemical disturbances in PCOS still remains elusive. Here, we showed an altered gut microbiota significantly correlated with PCOS phenotype. There were 33 patients with PCOS (non-obese PCOS individuals, PN, n = 12; obese PCOS individuals, PO, n = 21) as well as 15 control subjects (non-obese control individuals, CN, n = 9; obese control individuals, CO, n = 6) enrolled in our study. The plasma levels of serotonin, ghrelin, and peptide YY (PYY) were significantly decreased in patients with PCOS compared with controls, and have a significantly negative correlation with waist circumference and testosterone. Sequencing of the V3–V4 region of the 16S rRNA gene in fecal samples revealed the substantial differences of gut microbial species between the PCOS and non-obese controls. Bacterial species were clustered into 23 co-abundance groups (CAGs) based on the SparCC correlation coefficients of their relative abundance. The CAGs increased in PCOS, including the bacteria belonging to Bacteroides, Escherichia/Shigella and Streptococcus, were negatively correlated with ghrelin, and positively correlated with testosterone and BMI. Furthermore, the CAGs that were decreased in PCOS, including the bacteria from Akkermansia and Ruminococcaceae, showed opposite relationship with body-weight, sex-hormone, and brain–gut peptides. In conclusion, gut microbial dysbiosis in women with PCOS is associated with the disease phenotypes.

Introduction

Polycystic ovary syndrome is a common heterogeneous endocrine and chronic metabolic disease affecting 3–26% of reproductive-aged women by applying the different diagnostic criteria (Michelmore et al., 1999; March et al., 2010; Li et al., 2013). Women with PCOS are faced with increased risk of endometrial cancer, diabetes, metabolic syndrome and cardiovascular disease in the long term (Legro et al., 1999; Hardiman et al., 2003; Dumesic and Lobo, 2013). They showed hyperandrogenism, hyperinsulinemia, anovulation, and polycystic ovarian morphology (Landay et al., 2009), and had 2.8-fold higher prevalence of obesity compared to women without PCOS according to a recent meta-analysis (Lim et al., 2013). Although the genetic (Vink et al., 2006), gestational environment (Abbott et al., 1998) and lifestyle factors (Zhang J. et al., 2015) have been suggested to involve in the development of PCOS, the precise underlying triggers for these key biochemical and metabolic disturbances remain largely unclear.

Since high levels of circulating insulin stimulating the ovarian theca cells to produce androgen had been supported by various experiments, both in vitro and in vivo (Poretsky, 1991), a new microbiological hypothesis for the development of PCOS suggested that dysbiosis of gut microbiota can increase the ovaries production of androgens, and then interfere with normal follicle development via triggering a chronic inflammatory response and insulin resistance (Tremellen and Pearce, 2012). Further evidence suggests that gut microbiota and its metabolites have the ability to regulate inflammation pathway activation, brain–gut peptide secretion and islet β-cell proliferation, thus leading to abnormal or excessive fat accumulation, insulin resistance and compensatory hyperinsulinemia (Vrieze et al., 2010; Barber et al., 2015). The shifts of gut microbiota had been observed in the letrozole-induced PCOS mouse models (Guo et al., 2016; Kelley et al., 2016), and fecal microbiota transplantation from the letrozole-induced PCOS mice to germ-free mice indicated that gut microbiota also has a close relationship with its host’s sex hormone levels, estruscycles and ovarian morphological changes (Guo et al., 2016). In addition, according to another mouse model experiment, female recipients of male cecal microbiota displayed increased testosterone levels, compared with unmanipulated females and female recipients of female cecal micriobiota (Markle et al., 2013). On the other hand, some mediators of the brain–gut axis, such as serotonin, ghrelin and PYY, whose secretion were affected by gut microbiota (Samuel et al., 2008; Yano et al., 2015; Perry et al., 2016), were involved in appetite regulation, systemic energy homeostasis and LH secretion (Konturek et al., 2004; Ogata et al., 2009; Jun et al., 2015). But the experimental evidence of the relationship between gut microbiota, mediators of the brain–gut axis and metabolic phenotypes of PCOS in human is still absent.

In order to study the association of gut microbiota with obesity and non-obesity in relation to PCOS, we recruited non-obese or obese women with or without PCOS. We investigated the composition of gut microbiota in PCOS, and how the changes of gut microbiota correlated with the sex-hormones, brain–gut peptides, and various other metabolic parameters. Our study suggested an altered composition of gut microbiota of women with PCOS from controls, and identified important bacterial phylotypes which were associated with the clinical characteristics of PCOS.

Materials and Methods

Study Participants

A total of 48 premenopausal women aged 17–45 years were enrolled from individuals who visited the department of endocrinology or gynecology at Shanghai General Hospital (Shanghai, China) between January 2014 and July 2014. Diagnosis criteria for PCOS established by Chinese Medical Association was based on Rotterdam criteria and the currently available evidences in Han Chinese. In our study, to avoid the confusion due to the different diagnostic criteria, each participant had met the three features based on Rotterdam criteria and the Chinese criteria as follows: oligo- or anovulation, clinical and/or biochemical signs of hyperandrogenism, and polycystic ovaries (Rotterdam Eshre/Asrm-Sponsored Pcos Consensus Workshop Group, 2004). Women in control group had no history of diagnosed PCOS and did not meet any of the Rotterdam Criteria. Individuals with one of the followings were excluded from the study: pregnancy, androgen-secreting tumors, adrenal disorders, thyroid dysfunction, Cushing’s syndrome, smoking, hypertension, and lipid-dysregulation. Included participants had also no administration of hormonal medication, insulin sensitizer or antibiotics within the preceding 3 months. After submitting written informed consent, all participants with PCOS (n = 33) or without PCOS (n = 15) were then stratified according to BMI into the following four subgroups: obese PCOS patients (PO group, n = 21), non-obese PCOS patients (PN group, n = 12), obese controls (CO group, n = 6) and non-obese controls (CN group, n = 9). Obesity was defined as BMI ≥ 25 kg/m² and the non-obesity otherwise according to the World Health Organization (World Health Organization, 2000). The study protocol was approved by the Human Research Ethics Committee of Shanghai General Hospital (No. 2014KY091) ahead of procedure of enrollment. The trial was registered in Chinese Clinical Trial Registry with the registration number ChiCTR-TRC-14005075.

Sampling

All subjects were examined in the morning after overnight fasting (≥8 h). 4-(2-aminoethyl)-benzenesulfonyl fluoride was immediately added to the blood samples for enzyme-linked immunosorbent assay (ELISA) of ghrelin at a final concentration of 1 mg/mL. After clotting at room temperature for 30 min, the blood samples were centrifuged at 3,000 × g for 15 min at 4°C, and were transferred into separate tubes and acidified by adding 5 N HCl to a final concentration of 0.05 N. The processed serum samples were stored at -20°C until analysis. For PYY, blood samples were immediately added DPP-IV inhibitor and aprotinin to a final concentration 500 KIU/mL. Fecal samples were collected on the same day. Feces was divided into aliquots and was frozen on dry ice immediately upon collection and stored at -80°C until analysis.

Anthropometric and Metabolic Parameters Measurements

Waist circumference was measured midway between the lower costal margin and the iliac crest, and hip circumference at the level of maximum extension of the buttocks. WHR was calculated as waist circumference divided by hip circumference. BMI was calculated as body weight in kilograms divided by body height squared in meters.

Biochemical indexes were measured on automatic biochemical analyzer (Beckman AU5800 Clinical Chemistry System, Beckman Coulter, Inc., South Kraemer Boulevard, Brea, CA, USA). HbA1c level was measured by high-performance liquid chromatography (Bio-Rad Variant II Turbo, Bio-Rad Laboratories, Co., Ltd, Germany). FINS, FSH, LH, and testosterone were tested using an automated immunoassay system (TOSOHAIA-1800ST, TOSOH Corporation, Tokyo, Japan). The HOMA-IR and beta cell function (HOMA-beta) were calculated using the following formula:

$$\text{HOMA-IR} \text{=} \text{fasting} \text{plasma} \text{glucose} (\text{FPG}) (\text{mM}) \text{*} \text{FINS }(\text{mIU/L})\text{/22}.\text{5};\left(1\right)$$

$$\text{HOMA-beta} \text{=} \text{20} \text{*} \text{FINS} (\text{mIU/l})\text{/}(\text{FPG }(\text{mM}) \text{-3}.\text{5})\text{*100}\%.\left(2\right)$$

Plasma serotonin, ghrelin, and PYY levels were determined with the commercially available ELISA kit according to manufacturer’s instructions. The kit for serotonin was purchased from EUROIMMUNUS, Inc. (EA602/96, Mountain Lakes, NJ, USA), and kits for ghrelin and PYY were bought from EMD Millipore Corporation (EZGRT-89K, EZHPYYT66K, Saint Charles, MO, USA, respectively). The intra-assay and inter-assay coefficients of variation were <5 and <10, respectively.

The investigation of mental psychological health was conducted with SDS and SAS.

Fecal DNA Extraction and Sequencing

DNA extraction from each frozen fecal sample was conducted as previously described (Godon et al., 1997). All samples were sequenced on the Illumina Miseq platform (Illumina, Inc., USA) with Miseq reagent kit v3 (600-cycle) (MS-102-3033, Illumina, USA). Sequencing library of the V3–V4 regions of the 16S rRNA gene was constructed based on the manufacturer’s instruction (Part # 15044223 Rev. B, Illumina, USA), with some modifications as previously published (Zhang et al., 2016).

Statistical Analysis

Clinical Data Statistical Analysis

Statistical analyses were performed using SAS statistical software 9.2 (SAS Institute, Inc., Cary, NC, USA). Normal distribution of the data was calculated with the Kolmogorov–Smirnov test. Continuous variables were expressed as mean ± standard deviation (SD). Variables that were not normally distributed were rank transformed before analysis. Differences on clinical characteristics and metabolic biomarkers were evaluated using t-test or a one-way analysis of variance (ANOVA) for normally distributed variables or non-parameter Wilcoxon test for non-normally distributed, with multiple comparisons performed using SNK (Student-Newman-Keuls) correction. Partial correlation analysis was used to study the relationships between brain–gut peptides and other metabolic parameters. A multiple stepwise regression analysis was employed to identify influence factors of brain–gut peptides, using a stepwise selection method with a P-value of 0.15 at entry and a P-value of 0.15 at stay. All independent variables in the multiple regression were tested for multicollinearity. Statistical tests were based on two tailed probability. A P-value less than 0.05 was considered statistically significant. Power calculation was performed using PASS software 11.

Bioinformatics and Statistical Analysis

Quality control of the raw sequencing data was performed, as previously described (Zhang et al., 2016), to reduce spurious OTUs. OTUs were clustered and mapped using Uparse (Edgar, 2013) according to Uparse OTU analysis pipeline¹. Quality-filtered sequences were clustered into unique sequences and sorted in order of decreasing abundance to identify representative sequences, and singletons were omitted in this step. Chimera detection was performed using UCHIME (Edgar et al., 2011) against the Ribosomal Database Project (RDP) classifier (Cole et al., 2013). OTUs were clustered using a ≥97% similarity threshold. The representative sequences of each OTU were used to generate a phylogenetic tree using FastTree (Price et al., 2009), and were subjected to RDP for taxonomical assignments with a bootstrap cutoff of 80%.

The sequences of all the samples were downsized to 7500 (1000 permutations) to equal the difference in sequencing depth. The rarefaction estimates, the Shannon index, the Simpson’s diversity index and the Chao1 index were calculated using QIIME v1.2.1 (Caporaso et al., 2010), and the difference of alpha-diversity was evaluated by Kruskal–Wallis test using IBM SPSS Statistics 19. To evaluate the structural segregation of gut microbiota between individuals with/without PCOS, we performed CAP (Anderson and Willis, 2003) using FORTRAN program (Anderson, 2004). And MANOVA was estimated using MATLAB R2014a (The MathWorks, Natick, MA, USA). Both of CAP and MANOVA analysis were calculated based on the Bray–Curtis (Bray and Curtis, 1957) and the UniFrac distance (Lozupone and Knight, 2005) performed by QIIME.

To pick key OTUs that might be responsible for the segregation, RDA was performed using the R ‘vegan’ package. Statistical significance was assessed by Monte Carlo test with 9999 permutations.

The OTUs shared by at least 20% among all the samples were considered as prevalent OTUs. Correlations between prevalent OTUs were calculated using the SparCC algorithm (Friedman and Alm, 2012) based on their abundance. The correlation values were converted to a correlation distance (1-correlation value), and the OTUs were clustered using the Ward clustering algorithm via MATLAB. Permutational MANOVA was used to determine CAGs via dividing the cluster tree (9999 permutations, P < 0.001). Permutational MANOVA was performed using the R ‘vegan’ package² and the CAG network was visualized in Cytoscape. Spearman correlation between CAGs and clinical parameters was calculated using MATLAB, and Benjamini-Hochberg method was used to control false discovery rate (Benjamini and Hochberg, 1995). Kruskal–Wallis test was performed for comparing the difference of relative abundance at CAG level between the four groups using IBM SPSS Statistics 19.

Raw sequencing data of the 16S rRNA gene V3–V4 regions and accompanying information are available in Sequence Read Archive database under accession number SRP085887.

Results

General Anthropometric and Metabolic Parameters of All the Participants

The clinical data of all the 48 participants was summarized in Table 1. After adjusted P value for age, the differences between groups reflected in the following aspects: anthropometric parameters, sex hormones, metabolic parameters, mediators of the brain–gut axis and psychological scales. Compared to women without PCOS, women with PCOS showed significantly higher levels of testosterone and hirsutism score. PO group also had increased levels of metabolic parameters and SDS value than PN group, indicating that excess weight aggravated the metabolic abnormalities in women with PCOS. The level of serotonin was significantly higher in CN group than the other three groups. The level of ghrelin was lower in women with PCOS compared with women without PCOS, and it was also lower in obese women with PCOS compared with non-obese women with PCOS. PYY showed a significant decrease in PO group compared with CN and CO groups, while PN group was also decreased compared to CN and CO and it was higher than the PO group. Both PCOS patients and obese individuals possessed abnormal levels of the mediators of the brain–gut axis.

TABLE 1

TABLE 1. The clinical characteristics and metabolic profiles of all the participants.

In order to evaluate the association between the mediators of the brain–gut axis and the various metabolic characteristics of PCOS, we calculated the partial correlation between serotonin, PYY, ghrelin and the clinical parameters of all subjects after adjusting for age and BMI (Table 2). Plasma serotonin level was significantly negatively correlated with waist circumference, hirsutism score, testosterone and SAS value, whereas, significantly positively correlated with neutrocyte count, HOMA-beta and TCH. Plasma PYY level was significantly negatively correlated with waist circumference, WHR, hirsutism score, testosterone, PPG, ALT and AST, whereas, significantly positively correlated with neutrocyte count and HOMA-beta. Plasma ghrelin level was significantly negatively correlated with waist circumference, WHR, hirsutism score, testosterone, PPG, P2hINS, ALT, and AST. To explore whether plasma brain–gut peptide/indole levels were independently associated with the clinical parameters, multiple stepwise regression analysis involving all the parameters with significant correlations with the mediators of the brain–gut axis was performed (with sle = 0.15, and sls = 0.15). Adjusting for age and BMI, it revealed that the testosterone level (β, -2.5296; P < 0.001), neutrocyte count (β, 0.38753; P < 0.01) and waist circumference (β, -0.29091; P < 0.05) were independently associated with plasma serotonin levels; the waist circumference (β, -0.27497; P < 0.05), the PCOS group (β, -9.447; P < 0.05) and neutrocyte count (β, 0.41872; P < 0.01) were independently associated with plasma PYY levels; while the testosterone level (β, -1.3728; P < 0.05), the PCOS group (β, -12.6121; P < 0.001), BMI (β, -0.889; P < 0.01), and P2hINS (β, -0.3382; P < 0.01) were independently associated with plasma ghrelin levels in all participants. Thus, the levels of the brain–gut peptides and the brain–gut indole were correlated with the PCOS related phenotypes.

TABLE 2

TABLE 2. Partial correlations between the brain–gut peptide/indole levels and various clinical parameters after adjusted for age and BMI.

Dysbiosis of Gut Microbiota in Women with PCOS

We performed gene sequencing on the V3–V4 regions of the 16S rRNA to evaluate the dysbiosis of gut microbiota in PCOS. In total, 852,886 usable raw reads (264,105 unique sequences) were obtained from all 48 samples. 770,892 high-quality reads (14,480 ± 4,325 reads per sample) were clustered into 567 OTUs at 97% similarity level. According to observed OTUs and the Chao1 index, CN group owned the highest richness, followed by CO and PN group, and PO group showed the lowest richness of gut microbiota (Figures 1A,B), while no difference was found for the Shannon index and the Simpson’s diversity index. CAP and MANOVA based on the Bray–Curtis distance revealed a separation between CN group and the other three groups (Figures 1C,D). Notably, the gut microbiota composition of CO group was more similar with that of PN and PO groups, rather than CN group. CAP and MANOVA based on the UniFrac distance also showed similar trends (Supplementary Figure S1). Thus, the overall composition of the gut microbiota was disrupted in PCOS.

FIGURE 1

FIGURE 1. Overall structural differentiation of gut microbiota based on UniFrac distance between four groups. (A) OTU-level rarefaction (observed OTUs). (B) Chao 1 index. Values are shown by box-plot. Box represents the interquartile range. The line inside the box represents the median. And whiskers denote the minimum and maximum value. ^∗∗Adjusted P < 0.01 (Kruskal–Wallis test). (C) Bray–Curtis CAP using the first 13 PCs (accounting for 80.67% of the total variation). (D) Clustering of gut microbiota based on Bray–Curtis distance calculated with MANOVA using the first 13 PCs. ^∗∗P < 0.01. CN: non-obese control group, n = 9. CO: obese control group, n = 6. PN: non-obese PCOS group, n = 12. PO: obese PCOS group, n = 21.

We then used RDA to identify the gut bacteria responding to PCOS and obesity. RDA confirmed that both PCOS and obesity were significantly correlated to the dysbiosis of gut microbiota (Monte Carlo Permutation Procedure, P = 0.0002 and 0.0490, respectively). We identified 28 key OTUs that responded to the differentiation between the four groups according to RDA (Figure 2). OTU4 belonging to genus Bacteroides was enriched in PO group than CN and CO groups (P < 0.01 and P < 0.05, respectively), with PN group falling in between them. The average relative abundances of 21 OTUs were higher in CN group than PO group (P < 0.05), which mostly belonged to the genera Akkermansia, Bacteroides, Clostridium IV, Lactobacillus, Oscillibacter, or unclassified genera from family Ruminococcaceae, while in CO and PN groups, the abundances of these OTUs were lower than CN group but were higher than PO group. Five OTUs were enriched in CO group, which belonged to the genera Collinsella, Paraprevotella, Slackia, and unclassified genera from class Clostridia and kingdom Bacteria. OTU7 belonging to Bacteroides showed decreased abundance in CO group, compared with the other three groups (Figure 3).

FIGURE 2

FIGURE 2. The 28 key OTUs selected by RDA responding to PCOS or obesity. PCOS and obesity were conducted as environmental variables. Relative abundances (after log transformation) of all OTUs were used as response variables, and responding OTUs that had at least 20% of the variability explained by all canonical axes were selected and exhibited by black arrows. Statistical significance was assessed by Monte Carlo test with 9999 permutations.

FIGURE 3

FIGURE 3. Heat map of the relative abundances of 28 key OTUs related to the alteration of gut microbiota between the four groups. The color of the spots represents the relative abundance (normalized and log-transformed) of the OTUs. OTUs are ordered by spearman correlation analysis based on relative abundance. The genus-level taxonomic assignment is shown on the right. For CN, n = 9; for PN, n = 12; for CO, n = 6; and for PO, n = 21.

Association between Gut Microbiota and Clinical Parameters in PCOS and Control Subjects

As the bacteria act as interdependent functional groups (guilds) in the gut ecosystem (Zhang C. et al., 2015), we constructed a co-abundance network at OTU level. There were 225 OTUs shared by at least 20% of the samples were clustered into 23 CAGs based on SparCC correlation coefficients (Figure 4A; Supplementary Table S1). The results showed bacteria species in the same genus could have different response to the environment perturbations. For example, OTU5, 86, 506, 393, 491, and 336 belonging to genus Blautia were divided into CAG1 (both OTU5 and 506), 3, 15, 22, and 23, respectively. According to the Kruskal–Wallis test, 14 CAGs showed no significant difference between the four groups. CAG1, 4, and 7 were enriched in PO group rather than CN group. CAG 1 contained 14 OTUs, which belonged to genera Bacteroides, Escherichia/Shigella, Streptococcus, Blautia, and Parabacteroides. CAG4 contained 5 OTUs, which mainly belonged to Clostridium XlVa, Alistipes and Weissella. CAG7 contained OTUs belonging to Granulicatella, Streptococcus, Peptostreptococcus, and Rothia. While, CAG10, 11, 12, and 13 decreased in PO group rather than CN group. The dominant OTUs in CAG10-13 mostly belonged to Akkermansia, Alistipes, Coprococcus, and Ruminococcus. In CO and PN groups, the relative abundances of all these above CAGs fell in between those in PO and CN groups. In addition, CAG17 and 18 showed increased trends in CO group compared with the other groups (Figure 4B).

FIGURE 4

FIGURE 4. Bacterial correlation based on relative abundance in women with or without PCOS. (A) Co-abundance groups interaction network. Network shows correlation relationships between 23 CAGs of 225 OTUs from all samples. Node size represents the average abundance of each OTU. Lines between nodes represent correlations of each other, with the line width representing the correlation magnitude. The red ones represent positive correlations, and the blue ones represent negative correlations. Only lines whose absolute value of correlation coefficient greater than 0.50 are drawn, and unconnected nodes are omitted. (B) Group-level abundance differentiation of CAGs. Data are visualized by box-plot. Box represents the interquartile range. The line inside the box represents the median. And whiskers denote the minimum and maximum value. ^∗adjusted P < 0.05, ^∗∗adjusted P < 0.01, and ^∗∗∗adjusted P < 0.001 (multiple comparisons of Kruskal–Wallis test). For CN, n = 9; For PN, n = 12; For CO, n = 6; For PO, n = 21.

To explore the relationship between gut microbiota and host phenotype, Spearman correlation between 27 clinical parameters and 9 CAGs that showed significant difference among the four groups (Figure 5). The results showed significant correlation between gut microbiota and PCOS related symptoms, including obesity, inflammation, hyperinsulinmia, hyperandrogenism, abnormal ghrelin level, and psychological state. CAG1, 4 and 7, which were enriched in PO group, had a positive correlation with the disease phenotypes; while CAG10–13 and 18, which were enriched in CN group, had an opposite trend.

FIGURE 5

FIGURE 5. Associations between clinical parameters and gut microbiota. The color of spots represents R value of Spearman correlation between each CAG and clinical parameter. +FDR < 0.05, ++FDR < 0.01. For leucocyte, neutrocyte, lymphocyte and hirsutism, n = 46; for the other parameters, n = 48. BMI, body mass index; WHR, Waist hip ratio; Blood RT, blood routine test; FSH, follicular stimulating hormone; LH, luteinizing hormone; FPG, fasting plasma glucose; PPG, 2h postprandial plasma glucose; FINS, fasting plasma insulin; P2hINS, 2h postprandial plasma insulin; HOMA-IR, homeostasis model assessment for insulin resistance index; HOMA-beta, homeostasis model assessment for beta cell function; HbA1c, hemoglobin A1c; ALT, alanine aminotransferase; AST, aspartate transaminase; GGT, γ-glutamyltransferase; TCH, total cholesterol; TG, triglyceride; PYY, peptide YY; SDS, self-rating depression scale; SAS, self-rating anxiety scale.

Discussion

Our current study showed a broad view of the association between PCOS and gut microbiota in human. Our study, to the best of our knowledge, is the first clinical study that revealed the overall composition of gut microbiota in women with PCOS. In this study, we observed an altered gut microbial pattern in both obese and non-obese women with PCOS against non-obese control individuals, and found that the gut microbiota of PCOS shares similar compositional dysbiosis with that of obesity. Furthermore, the altered gut microbiota of PCOS was associated with the metabolic parameters, the sex-hormones and the mediators of the brain–gut axis.

Polycystic ovary syndrome leads to several complications, which include hyperandrogenism, obesity, and metabolic syndrome. In the current work, compared to women with obesity or PCOS alone, obese women with PCOS had more deleterious metabolic disorder, including insulin resistance and dyslipidemia, in accord with previous studies (Lim et al., 2013; Behboudigandevani et al., 2016). Based on our data, changes of gut microbiota also had trends similar to the variations of metabolic symptoms. Women with obesity and PCOS shared similar dysbiosis of gut microbiota. When compared with women with obesity or PCOS alone, obese women with PCOS possessed a more serious dysbiosis of gut microbiota. We observed that some gram-negative bacteria belonging to the genera Bacteroides and Escherichia/Shigella significantly increased in the gut of PCOS women with obesity. Bacterial LPS produced by gram-negative bacteria has been demonstrate to induce chronic inflammation, obesity, and insulin resistance in LPS-infused mice (Cani et al., 2007). While various experiments, both in vitro and in vivo, have verified that high levels of circulating insulin stimulates the ovarian theca cells to produce androgen (Poretsky, 1991). It suggests that gut microbiota antigen might also associate to the metabolic disorder in PCOS. On the contrary, we found that Akkermensia decreased in obesity and PCOS, which can restore host mucus layer and has been considered to contribute to the reduction in metabolic endotoxemia (Everard et al., 2013). Hence, the dysbiosis of gut microbiota, as increasing LPS-producing bacteria and decreasing protective bacteria, may correlate with the development of metabolic disorder in PCOS.

Some mediators of the brain–gut axis, such as serotonin and ghrelin were involved in appetite regulation (Konturek et al., 2004; Jun et al., 2015) as well as psychological well-being in PCOS (Lang et al., 2015). Moreover, intracerebroventricular injection of ghrelin significantly reduced serum LH level and pulsatile frequency in rats (Ogata et al., 2009). However, the relationship between the levels of serum mediators of the brain–gut axis and PCOS remains controversial in human studies. For example, the ghrelin level had significant difference between PCOS and BMI-matched controls in a Spanish study (Micic et al., 2007), but this orexigenic hormone showed association only with body weight in an Australian PCOS cohort (Samantha et al., 2014). The fasting PYY level also had conflicting results from different cohort (Zwirska-Korczala et al., 2008; Lin et al., 2015). Our study showed that women with PCOS had a significant decrease of serotonin, ghrelin, and PYY level compared with controls, and obesity enhanced the difference of serotonin and ghrelin in PCOS. Moreover, we find that the levels of serotonin, PYY, and ghrelin had a negative correlation with PCOS related parameters, such as waist circumference and testosterone. Our results indicated that the mediators of the brain–gut axis are associated with PCOS.

The secretion of some mediators of the brain–gut axis is regulated by gut microbiota. It has been shown that serotonin biosynthesis is promoted by spore-forming bacteria isolated from the mouse and human gut microbiota, known to be dominated by Clostridial species, while species belonging to Bacteroides didn’t show this function (Yano et al., 2015). In our study, most of the decreased key OTUs were Clostridial species, and the increased key OTU was Bacteroides spp. in both PCOS and obesity cohorts. PYY secretion is influenced by short chain fatty acid (SCFA) produced by gut microbiota (Tremaroli and Bäckhed, 2012). In the current, we also found the association between ghrelin and gut microbiota. For example, CAG1, which includes Bacteroides, Escherichia/Shigella and Blautia, was negatively correlated with ghrelin, while CAG10, including Akkermansia, was positively correlated with ghrelin. But we did not find a significant association between serotonin, PYY and the gut microbiota, which might be explained by the limitation of small sample size and require further studies. Regardless, according to existing evidences, regulating the mediators of the brain–gut axis might be one of the pathways for gut microbiota to effect on PCOS.

This cross-sectional study revealed the association between gut microbiota and PCOS. The dysbiosis of gut microbiota might contribute to the development of PCOS through the possible mechanisms discussed above. Due to the limited sample size (power calculation: the power values for BMI, waist circumference, testosterone level and fasting insulin level were 99.0, 96.9, 89.5, and 70.4% respectively, α = 0.0083) and the observation of association but not causality, further studies are needed to enlarge sample size and explore the potential causal mechanisms between gut microbiota and PCOS, such as the metabolites produced by gut microbiota and fecal microbiota transplantation.

Taken together, our work suggested an altered composition of gut microbiota, represented by the reduction of alpha diversity, the increase of LPS-producing bacteria, and the decrease of spore-forming species. This study also demonstrated the association between gut microbiota and PCOS-related clinical parameters. The previous study did not find any association of the saliva microbiota and PCOS phenotypes (Lindheim et al., 2016), indicating that the gut microbiota has a closer relationship with PCOS. It laid a foundation of research on the interaction between gut microbiota and its host with PCOS.

Ethics Statement

This study was carried out in accordance with the recommendations of Guideline on Informed Consent, Human Research Ethics Committee of Shanghai General Hospital with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the Human Research Ethics Committee of Shanghai General Hospital.

Author Contributions

The authors XD, YP, CZ, and RL conceived the study and designed the experiments. YS, JS, FZ, LL, XW, YL, HF, and WD performed the technical procedures. RL and XD provided the interpretation of data and wrote the paper. AG, AR, LZ, and XD revised the final draft of the paper. All authors read and approved the final version of the paper.

Funding

The study was supported by the Key Projects of Shanghai Municipal Health Bureau Research Fund (201440033), Shanghai Jiao Tong University Research Funding on Medical, Engineering Interdisciplinary Projects (YG2015MS30), the National Natural Science Foundation of China (81401141), the Science and Technology Commission of Shanghai Municipality (14YF1402200).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgment

We are grateful to all the participants and the medical staffs in the study.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00324/full#supplementary-material

Abbreviations

ALT, alanine aminotransferase; AST, aspartate transaminase; Blood RT, blood routine test; BMI, body mass index; CAG, co-abundance group; CAP, canonical analysis of principal coordinates; CN, control with non-obesity; CO, control with obesity; FINS, fasting plasma insulin; FSH, follicular stimulating hormone; FPG, fasting plasma glucose; GGT, γ-glutamyltransferase; HbA1c, hemoglobin A1c; HOMA-beta, homeostasis model assessment for beta cell function; HOMA-IR, homeostasis model assessment for insulin resistance index; LH, luteinizing hormone; LPS, lipopolysaccharide; MANOVA, multivariate analysis of variance; OTU, operational taxonomic units; P2hINS, 2h postprandial plasma insulin; PCOS, polycystic ovary syndrome; PN, PCOS with non-obesity; PO, PCOS with obesity; PPG, 2h postprandial plasma glucose; PYY, peptide YY; RDA, redundancy analysis; SAS, self-rating anxiety scale; SDS, self-rating depression scale; TCH, total cholesterol; TG, triglyceride; WHR, waist hip ratio.

Footnotes

1. ^ http://drive5.com/usearch/manual/uparse_pipeline.html
2. ^ http://www.r-project.org/

References

Abbott, D. H., Dumesic, D. A., Eisner, J. R., Colman, R. J., and Kemnitz, J. W. (1998). Insights into the development of polycystic ovary syndrome (PCOS) from studies of prenatally androgenized female rhesus monkeys. Trends Endocrinol. Metab. 9, 62–67. doi: 10.1016/S1043-2760(98)00019-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, M. (2004). CAP: A FORTRAN Computer Program for Canonical Analysis of Principal Coordinates. Auckland: University of Auckland, 14.

Google Scholar

Anderson, M. J., and Willis, T. J. (2003). Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology 84, 511–525. doi: 10.1890/0012-9658(2003)084[0511:CAOPCA]2.0.CO;2

CrossRef Full Text | Google Scholar

Barber, T. M., Dimitriadis, G. K., Andreou, A., and Franks, S. (2015). Polycystic ovary syndrome: insight into pathogenesis and a common association with insulin resistance. Clin. Med. 15, s72–s76. doi: 10.7861/clinmedicine.15-6-s72

PubMed Abstract | CrossRef Full Text | Google Scholar

Behboudigandevani, S., Tehrani, F. R., Dovom, M. R., Farahmand, M., Khomami, M. B., Noroozzadeh, M., et al. (2016). Insulin resistance in obesity and polycystic ovary syndrome: systematic review and meta-analysis of observational studies. Gynecol. Endocrinol. 32, 1–11. doi: 10.3109/09513590.2015.1117069

PubMed Abstract | CrossRef Full Text

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate - a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300.

Google Scholar

Bray, J. R., and Curtis, J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27, 325–349. doi: 10.2307/1942268

CrossRef Full Text | Google Scholar

Cani, P. D., Amar, J., Iglesias, M. A., Poggi, M., Knauf, C., Bastelica, D., et al. (2007). Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes Metab. Res. Rev. 56, 1761–1772. doi: 10.2337/db06-1491

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Cole, J. R., Wang, Q., Fish, J. A., Chai, B., McGarrell, D. M., Sun, Y., et al. (2013). Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res. 42, D633–D642. doi: 10.1093/nar/gkt1244

PubMed Abstract | CrossRef Full Text | Google Scholar

Dumesic, D. A., and Lobo, R. A. (2013). Cancer risk and PCOS. Steroids 78, 782–785. doi: 10.1016/j.steroids.2013.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | CrossRef Full Text | Google Scholar

Everard, A., Belzer, C., Geurts, L., Ouwerkerk, J. P., Druart, C., Bindels, L. B., et al. (2013). Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc. Natl. Acad. Sci. U.S.A. 110, 9066–9071. doi: 10.1073/pnas.1219451110

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedman, J., and Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687. doi: 10.1371/journal.pcbi.1002687

PubMed Abstract | CrossRef Full Text | Google Scholar

Godon, J.-J., Zumstein, E., Dabert, P., Habouzit, F., and Moletta, R. (1997). Molecular microbial diversity of an anaerobic digestor as determined by small-subunit rDNA sequence analysis. Appl. Environ. Microbiol. 63, 2802–2813.

PubMed Abstract | Google Scholar

Guo, Y., Qi, Y., Yang, X., Zhao, L., Wen, S., Liu, Y., et al. (2016). Association between polycystic ovary syndrome and gut microbiota. PLoS ONE 11:e0153196. doi: 10.1371/journal.pone.0153196

PubMed Abstract | CrossRef Full Text | Google Scholar

Hardiman, P., Pillay, O. S., and Atiomo, W. (2003). Polycystic ovary syndrome and endometrial carcinoma. Lancet 361, 1810–1812. doi: 10.1016/S0140-6736(03)13409-5

CrossRef Full Text | Google Scholar

Jun, N., Hail, K., and Sangkyu, P. (2015). Peripheral serotonin: a new player in systemic energy homeostasis. Mol. Cells 38, 1023–1028. doi: 10.14348/molcells.2015.0258

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, S. T., Skarra, D. V., Rivera, A. J., and Thackray, V. G. (2016). The gut microbiome is altered in a letrozole-induced mouse model of polycystic ovary syndrome. PLoS ONE 11:e0146509. doi: 10.1371/journal.pone.0146509

PubMed Abstract | CrossRef Full Text | Google Scholar

Konturek, S. J., Konturek, J. W., Pawlik, T., and Brzozowski, T. (2004). Brain-gut axis and its role in the control of food intake. J. Physiol. Pharmacol. 55, 137–154.

Google Scholar

Landay, M., Huang, A., and Azziz, R. (2009). Degree of hyperinsulinemia, independent of androgen levels, is an important determinant of the severity of hirsutism in PCOS. Fertil. Steril. 92, 643–647. doi: 10.1016/j.fertnstert.2008.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, U. E., Beglinger, C., Schweinfurth, N., Walter, M., and Borgwardt, S. (2015). Nutritional aspects of depression. Cell. Physiol. Biochem. 37, 1029–1043. doi: 10.1159/000430229

PubMed Abstract | CrossRef Full Text | Google Scholar

Legro, R. S., Kunselman, A. R., Dodson, W. C., and Dunaif, A. (1999). Prevalence and predictors of risk for type 2 diabetes mellitus and impaired glucose tolerance in polycystic ovary syndrome: a prospective, controlled study in 254 affected women 1. J. Clin. Endocrinol. Metab. 84, 165–169. doi: 10.1210/jc.84.1.165

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R., Zhang, Q., Yang, D., Li, S., Lu, S., Wu, X., et al. (2013). Prevalence of polycystic ovary syndrome in women in China: a large community-based study. Hum. Reprod. 28, 2562–2569. doi: 10.1093/humrep/det262

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, S. S., Norman, R., Davies, M., and Moran, L. (2013). The effect of obesity on polycystic ovary syndrome: a systematic review and meta-analysis. Obes. Rev. 14, 95–109. doi: 10.1111/j.1467-789X.2012.01053.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, T., Li, S., Xu, H., Zhou, H., Feng, R., Liu, W., et al. (2015). Gastrointestinal hormone secretion in women with polycystic ovary syndrome: an observational study. Hum. Reprod. 30, 2639–2644. doi: 10.1093/humrep/dev231

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindheim, L., Bashir, M., Munzker, J., Trummer, C., Zachhuber, V., Pieber, T. R., et al. (2016). The salivary microbiome in polycystic ovary syndrome (pcos) and its association with disease-related parameters: a pilot study. Front. Microbiol. 7:1270. doi: 10.3389/fmicb.2016.01270

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozupone, C., and Knight, R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235. doi: 10.1128/AEM.71.12.8228-8235.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

March, W. A., Moore, V. M., Willson, K. J., Phillips, D. I., Norman, R. J., and Davies, M. J. (2010). The prevalence of polycystic ovary syndrome in a community sample assessed under contrasting diagnostic criteria. Hum. Reprod. 25, 544–551. doi: 10.1093/humrep/dep399

PubMed Abstract | CrossRef Full Text | Google Scholar

Markle, J. G., Frank, D. N., Mortin-Toth, S., Robertson, C. E., Feazel, L. M., Rolle-Kampczyk, U., et al. (2013). Sex differences in the gut microbiome drive hormone-dependent regulation of autoimmunity. Science 339, 1084–1088. doi: 10.1126/science.1233521

PubMed Abstract | CrossRef Full Text | Google Scholar

Michelmore, K., Balen, A., Dunger, D., and Vessey, M. (1999). Polycystic ovaries and associated clinical and biochemical features in young women. Clin. Endocrinol. 51, 779–786. doi: 10.1046/j.1365-2265.1999.00886.x

CrossRef Full Text | Google Scholar

Micic, D., Sumaracdumanovic, M., Kendereski, A., Cvijovic, G., Zoric, S., Pejkovic, D., et al. (2007). Total ghrelin levels during acute insulin infusion in patients with polycystic ovary syndrome. J. Endocrinol. Invest. 30, 820–827. doi: 10.1007/BF03349222

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogata, R., Matsuzaki, T., Iwasa, T., Kiyokawa, M., Tanaka, N., Kuwahara, A., et al. (2009). Hypothalamic Ghrelin suppresses pulsatile secretion of luteinizing hormone via beta-endorphin in ovariectomized rats. Neuroendocrinology 90, 364–370. doi: 10.1159/000257421

PubMed Abstract | CrossRef Full Text | Google Scholar

Perry, R. J., Peng, L., Barry, N. A., Cline, G. W., Zhang, D., Cardone, R. L., et al. (2016). Acetate mediates a microbiome-brain-β-cell axis to promote metabolic syndrome. Nature 534, 213–217. doi: 10.1038/nature18309

PubMed Abstract | CrossRef Full Text | Google Scholar

Poretsky, L. (1991). On the paradox of insulin-induced hyperandrogenism in insulin-resistant states. Endocr. Rev. 12, 3–13. doi: 10.1210/edrv-12-1-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | CrossRef Full Text | Google Scholar

Rotterdam Eshre/Asrm-Sponsored Pcos Consensus Workshop Group (2004). Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome. Fertil. Steril. 81, 19–25. doi: 10.1016/j.fertnstert.2003.10.004

CrossRef Full Text | Google Scholar

Samantha, C., Teede, H. J., Harrison, C. L., Joham, A. E., Moran, L. J., and Stepto, N. K. (2014). Biomarkers and insulin sensitivity in women with polycystic ovary syndrome: characteristics and predictive capacity. Clin. Endocrinol. 83, 88–92. doi: 10.1111/cen.12619

PubMed Abstract | CrossRef Full Text | Google Scholar

Samuel, B. S., Shaito, A., Motoike, T., Rey, F. E., Backhed, F., Manchester, J. K., et al. (2008). Effects of the gut microbiota on host adiposity are modulated by the short-chain fatty-acid binding G protein-coupled receptor, Gpr41. Proc. Natl. Acad. Sci. U.S.A. 105, 16767–16772. doi: 10.1073/pnas.0808567105

PubMed Abstract | CrossRef Full Text | Google Scholar

Tremaroli, V., and Bäckhed, F. (2012). Functional interactions between the gut microbiota and host metabolism. Nature 489, 242–249. doi: 10.1038/nature11552

PubMed Abstract | CrossRef Full Text | Google Scholar

Tremellen, K., and Pearce, K. (2012). Dysbiosis of gut microbiota (DOGMA)–a novel theory for the development of polycystic ovarian syndrome. Med. Hypotheses 79, 104–112. doi: 10.1016/j.mehy.2012.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Vink, J. M., Sadrzadeh, S., Lambalk, C. B., and Boomsma, D. I. (2006). Heritability of polycystic ovary syndrome in a Dutch twin-family study. J. Clin. Endocrinol. Metab. 91, 2100–2104. doi: 10.1210/jc.2005-1494

PubMed Abstract | CrossRef Full Text | Google Scholar

Vrieze, A., Holleman, F., Zoetendal, E. G., de Vos, W. M., Hoekstra, J. B., and Nieuwdorp, M. (2010). The environment within: how gut microbiota may influence metabolism and body composition. Diabetologia 53, 606–613. doi: 10.1007/s00125-010-1662-7

PubMed Abstract | CrossRef Full Text | Google Scholar

World Health Organization (2000). Obesity: Preventing and Managing the Global Epidemic. Geneva: World Health Organization. doi: 10.1016/j.cell.2015.02.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Yano, J. M., Yu, K., Donaldson, G. P., Shastri, G. G., Ann, P., Ma, L., et al. (2015). Indigenous bacteria from the gut microbiota regulate host serotonin biosynthesis. Cell 161, 264–276. doi: 10.1016/j.cell.2015.02.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Yin, A., Li, H., Wang, R., Wu, G., Shen, J., et al. (2015). Dietary modulation of gut microbiota contributes to alleviation of both genetic and simple obesity in children. EBioMedicine 2, 968–984. doi: 10.1016/j.ebiom.2015.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Liu, Y., Liu, X., Xu, L., Zhou, L., Tang, L., et al. (2015). High intake of energy and fat in Southwest Chinese women with PCOS: a population-based case-control study. PLoS ONE 10:e0127094. doi: 10.1371/journal.pone.0127094

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Wu, Y., Wang, J., Wu, G., Long, W., Xue, Z., et al. (2016). Accelerated dysbiosis of gut microbiota during aggravation of DSS-induced colitis by a butyrate-producing bacterium. Sci. Rep. 6:27572. doi: 10.1038/srep27572

PubMed Abstract | CrossRef Full Text | Google Scholar

Zwirska-Korczala, K., Sodowski, K., Konturek, S., Kuka, D., Kukla, M., Brzozowski, T., et al. (2008). Postprandial response of ghrelin and PYY and indices of low-grade chronic inflammation in lean young women with polycystic ovary syndrome. J. Physiol. Pharmacol. 59(Suppl. 2), 161–178.

PubMed Abstract | Google Scholar