Europe PMC
Nothing Special   »   [go: up one dir, main page]

Europe PMC requires Javascript to function effectively.

Either your web browser doesn't support Javascript or it is currently turned off. In the latter case, please turn on Javascript support in your web browser and reload this page.

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

Abstract 


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.

Free full text 


Logo of frontmicrobioLink to Publisher's site
Front Microbiol. 2017; 8: 324.
Published online 2017 Feb 28. https://doi.org/10.3389/fmicb.2017.00324
PMCID: PMC5328957
PMID: 28293234

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

Abstract

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.

Keywords: gut microbiota, polycystic ovary syndrome, obesity, testosterone, ghrelin, serotonin

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/m2 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:

HOMA-IR = fasting plasma glucose (FPG) (mM) * FINS (mIU/L)/22.5;
(1)
HOMA-beta = 20 * FINS (mIU/l)/(FPG (mM) -3.5)*100%.
(2)

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 pipeline1. 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’ package2 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 Table11. 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

The clinical characteristics and metabolic profiles of all the participants.

ParametersControls (n = 15)
PCOS (n = 33)
Non-obese, n = 9Obese, n = 6Non-obese, n = 12Obese, n = 21FP adjusted for age
Age (years)32.2 ± 5.933 ± 5.425.5 ± 4.329.3 ± 6.53.3/
Anthropometric parameters
BMI (kg/m2)21.9 ± 2.2b27.5 ± 3.3a21.6 ± 2.2b30 ± 3.6a46<0.0001
Waist circumference (cm)68.8 ± 5.7c84.7 ± 12.3b75 ± 7.8c95.1 ± 9.2a24<0.0001
WHR0.78 ± 0.05c0.8 ± 0.1b0.79 ± 0.06bc0.92 ± 0.05a20<0.0001
Hirsutism score8 ± 0b8 ± 0b17.5 ± 6.5a16.6 ± 5.7a31<0.0001
Blood RT
Leukocyte (10^9/L)6.8 ± 26.4 ± 17.1 ± 2.17.3 ± 1.70.50.49
Neutrocyte (10^9/L)4.1 ± 0.84.1 ± 1.24.5 ± 24.6 ± 1.20.50.6638
Lymphocyte (10^9/L)2.1 ± 0.92.3 ± 0.62 ± 0.42.3 ± 0.70.40.7929
Sex hormones
FSH (IU/L)5.3 ± 2.27 ± 0.97.2 ± 2.36.4 ± 2.51.60.2072
LH (IU/L)7.5 ± 4.1b11.3 ± 5.7ab15.7 ± 8.5a19.4 ± 31.2ab3.20.0335
Testosterone (nmol/L)0.79 ± 0.41b1 ± 0.4b4.5 ± 2a5.4 ± 2.7a14<0.0001
LH/FSH1.5 ± 0.71.6 ± 0.72.3 ± 1.12.5 ± 22.10.1117
Metabolic parameters
FPG (mmol/L)5 ± 0.45.2 ± 0.44.7 ± 0.55 ± 2.12.70.0591
PPG (mmol/L)6 ± 0.9b5.5 ± 1b6 ± 1.3b8.5 ± 2.6a8.80.0001
FINS (uU/mL)11.9 ± 4a15.6 ± 5a5.5 ± 2.4b14.3 ± 7.4a11<0.0001
P2hINS (uU/mL)54.1 ± 31.5b83.7 ± 52.4ab65.5 ± 53.3b163 ± 99.3a7.70.0003
HOMA-IR1.7 ± 1.1a3.5 ± 1a1.1 ± 0.5b3.3 ± 2.1a11<0.0001
HOMA-beta159.8 ± 32.7a203.7 ± 107a174.8 ± 284.2a339.8 ± 254.8a3.30.0279
HbA1c (%)5 ± 0.3b5.3 ± 0.3b5.2 ± 0.3b5.9 ± 1.1a9.4<0.0001
ALT (U/L)14 ± 9.6b16.2 ± 10.7b17.5 ± 8.9ab32.9 ± 28.3a4.10.0123
AST (U/L)15.9 ± 1.916.5 ± 5.419.4 ± 5.527.7 ± 23.62.10.1198
GGT (U/L)14.1 ± 3.6b19.2 ± 10.4ab13.8 ± 3.4b26.3 ± 23.6a4.80.0057
TCH (mmol/L)4.5 ± 0.94.8 ± 0.94.4 ± 1.14.6 ± 0.70.20.8762
TG (mmol/L)1 ± 0.6b1.3 ± 1ab0.83 ± 0.38b1.5 ± 0.7a7.10.0005
Mediators of the brain–gut axis
Serotonin (ng/mL)294.2 ± 91.8a199.6 ± 38.1b211.5 ± 69b193.6 ± 54.4b4.10.012
Ghrelin (ng/mL)0.51 ± 0.08a0.44 ± 0.04ab0.4 ± 0.1b0.24 ± 0.05c27<0.0001
PYY (pg/mL)79.8 ± 24.1a79.2 ± 17.5a56.8 ± 16.1ab62.1 ± 48.4b4.90.0052
Psychological scales
SDS value39.2 ± 8.1b42.4 ± 12.4ab36.9 ± 6.9b50.1 ± 9.4a6.60.0009
SAS value32.7 ± 5.5b34.2 ± 7.2ab35.9 ± 4.5ab42.8 ± 9.1a5.40.003
Data are presented as mean ± SD. Multiple comparison: means with the same letter (a, b, or c) are not significantly different, with different letters are significantly different (adjusted P < 0.05). P-value by ANOVA for normally distributed variables or non-parameter Wilcoxon test for non-normally distributed. The parameters of age, waist circumference, leukocyte, testosterone, AST, ALT, TCH, and SDS value were normally distributed. PCOS, polycystic ovary syndrome; 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.

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 Table22). 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

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

ParametersSerotonin
PYY
Ghrelin
rPrPrP
Anthropometric parameters
Waist circumference-0.3310.025-0.5070.000-0.2970.045
WHR-0.2790.061-0.3790.010-0.3360.023
Hirsutism score-0.3400.026-0.4210.005-0.3810.012
Blood RT
Leukocyte count0.2650.0820.1640.287-0.0080.959
Neutrocyte count0.3220.0330.3500.020-0.0620.689
Lymphocyte count0.2050.1810.0280.8560.0590.706
Sex hormones
FSH0.0530.7260.2370.113-0.0720.636
LH-0.1270.399-0.1360.369-0.2220.139
Testosterone-0.3120.035-0.3710.011-0.627<0.0001
LH/FSH-0.2210.140-0.2550.087-0.1870.213
Metabolic parameters
FPG-0.1570.298-0.0650.669-0.0120.939
PPG-0.0470.757-0.3960.007-0.4180.004
FINS0.1730.2510.0830.5830.0860.568
P2hINS0.0620.682-0.0240.874-0.3380.022
HOMA-beta0.3510.0170.3170.0320.1570.298
HbA1c0.0240.874-0.1640.276-0.1110.462
ALT0.0140.926-0.2910.050-0.4180.004
AST0.0420.784-0.3030.041-0.3010.042
TCH0.3120.035-0.1560.299-0.0610.688
TG0.1150.447-0.2530.090-0.1860.215
Mediators of the brain–gut axis
Serotonin1.000/0.2820.0580.2150.151
Ghrelin0.2150.1510.2380.1121.000/
PYY0.2820.0581.000/0.2380.112
Psychological scales
SDS value-0.1850.219-0.1350.370-0.0340.825
SAS value-0.3550.016-0.2310.122-0.0640.671
After taking account of age and BMI, bold values indicate statistical significance determined by Spearman partial correlation analysis. 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-beta, homeostasis model assessment for beta cell function; HbA1c, hemoglobin A1c; ALT, alanine aminotransferase; AST, aspartate transaminase; TCH, total cholesterol; TG, triglyceride; PYY, peptide YY; SDS, self-rating depression scale; SAS, self-rating anxiety scale; CN, control with non-obesity; CO, control with obesity; PN, PCOS with non-obesity; PO, PCOS with obesity.

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.

An external file that holds a picture, illustration, etc.
Object name is fmicb-08-00324-g001.jpg

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 Figure22). 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 Figure33).

An external file that holds a picture, illustration, etc.
Object name is fmicb-08-00324-g002.jpg

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.

An external file that holds a picture, illustration, etc.
Object name is fmicb-08-00324-g003.jpg

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 Figure4A4A; 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 Figure4B4B).

An external file that holds a picture, illustration, etc.
Object name is fmicb-08-00324-g004.jpg

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 Figure55). 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.

An external file that holds a picture, illustration, etc.
Object name is fmicb-08-00324-g005.jpg

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.

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.

Acknowledgments

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

Abbreviations

ALTalanine aminotransferase
ASTaspartate transaminase
Blood RTblood routine test
BMIbody mass index
CAGco-abundance group
CAPcanonical analysis of principal coordinates
CNcontrol with non-obesity
COcontrol with obesity
FINSfasting plasma insulin
FSHfollicular stimulating hormone
FPGfasting plasma glucose
GGTγ-glutamyltransferase
HbA1chemoglobin A1c
HOMA-betahomeostasis model assessment for beta cell function
HOMA-IRhomeostasis model assessment for insulin resistance index
LHluteinizing hormone
LPSlipopolysaccharide
MANOVAmultivariate analysis of variance
OTUoperational taxonomic units
P2hINS2h postprandial plasma insulin
PCOSpolycystic ovary syndrome
PNPCOS with non-obesity
POPCOS with obesity
PPG2h postprandial plasma glucose
PYYpeptide YY
RDAredundancy analysis
SASself-rating anxiety scale
SDSself-rating depression scale
TCHtotal cholesterol
TGtriglyceride
WHRwaist hip ratio

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).

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

References

  • Abbott D. H., Dumesic D. A., Eisner J. R., Colman R. J., 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. 10.1016/S1043-2760(98)00019-8 [Abstract] [CrossRef] [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., Willis T. J. (2003). Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology 84 511–525. 10.1890/0012-9658(2003)084[0511:CAOPCA]2.0.CO;2 [CrossRef] [Google Scholar]
  • Barber T. M., Dimitriadis G. K., Andreou A., Franks S. (2015). Polycystic ovary syndrome: insight into pathogenesis and a common association with insulin resistance. Clin. Med. 15 s72–s76. 10.7861/clinmedicine.15-6-s72 [Abstract] [CrossRef] [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. 10.3109/09513590.2015.1117069 [Abstract] [CrossRef] [Google Scholar]
  • Benjamini Y., 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., Curtis J. T. (1957). An ordination of the upland forest communities of southern Wisconsin. Ecol. Monogr. 27 325–349. 10.2307/1942268 [CrossRef] [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. 10.2337/db06-1491 [Abstract] [CrossRef] [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. 10.1038/nmeth.f.303 [Europe PMC free article] [Abstract] [CrossRef] [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. 10.1093/nar/gkt1244 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Dumesic D. A., Lobo R. A. (2013). Cancer risk and PCOS. Steroids 78 782–785. 10.1016/j.steroids.2013.04.004 [Abstract] [CrossRef] [Google Scholar]
  • Edgar R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10 996–998. 10.1038/nmeth.2604 [Abstract] [CrossRef] [Google Scholar]
  • Edgar R. C., Haas B. J., Clemente J. C., Quince C., Knight R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27 2194–2200. 10.1093/bioinformatics/btr381 [Europe PMC free article] [Abstract] [CrossRef] [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. 10.1073/pnas.1219451110 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Friedman J., Alm E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687 10.1371/journal.pcbi.1002687 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Godon J.-J., Zumstein E., Dabert P., Habouzit F., Moletta R. (1997). Molecular microbial diversity of an anaerobic digestor as determined by small-subunit rDNA sequence analysis. Appl. Environ. Microbiol. 63 2802–2813. [Europe PMC free article] [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 10.1371/journal.pone.0153196 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Hardiman P., Pillay O. S., Atiomo W. (2003). Polycystic ovary syndrome and endometrial carcinoma. Lancet 361 1810–1812. 10.1016/S0140-6736(03)13409-5 [Abstract] [CrossRef] [Google Scholar]
  • Jun N., Hail K., Sangkyu P. (2015). Peripheral serotonin: a new player in systemic energy homeostasis. Mol. Cells 38 1023–1028. 10.14348/molcells.2015.0258 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Kelley S. T., Skarra D. V., Rivera A. J., Thackray V. G. (2016). The gut microbiome is altered in a letrozole-induced mouse model of polycystic ovary syndrome. PLoS ONE 11:e0146509 10.1371/journal.pone.0146509 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Konturek S. J., Konturek J. W., Pawlik T., Brzozowski T. (2004). Brain-gut axis and its role in the control of food intake. J. Physiol. Pharmacol. 55 137–154. [Abstract] [Google Scholar]
  • Landay M., Huang A., 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. 10.1016/j.fertnstert.2008.06.021 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Lang U. E., Beglinger C., Schweinfurth N., Walter M., Borgwardt S. (2015). Nutritional aspects of depression. Cell. Physiol. Biochem. 37 1029–1043. 10.1159/000430229 [Abstract] [CrossRef] [Google Scholar]
  • Legro R. S., Kunselman A. R., Dodson W. C., 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. 10.1210/jc.84.1.165 [Abstract] [CrossRef] [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. 10.1093/humrep/det262 [Abstract] [CrossRef] [Google Scholar]
  • Lim S. S., Norman R., Davies M., Moran L. (2013). The effect of obesity on polycystic ovary syndrome: a systematic review and meta-analysis. Obes. Rev. 14 95–109. 10.1111/j.1467-789X.2012.01053.x [Abstract] [CrossRef] [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. 10.1093/humrep/dev231 [Abstract] [CrossRef] [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 10.3389/fmicb.2016.01270 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Lozupone C., Knight R. (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71 8228–8235. 10.1128/AEM.71.12.8228-8235.2005 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • March W. A., Moore V. M., Willson K. J., Phillips D. I., Norman R. J., Davies M. J. (2010). The prevalence of polycystic ovary syndrome in a community sample assessed under contrasting diagnostic criteria. Hum. Reprod. 25 544–551. 10.1093/humrep/dep399 [Abstract] [CrossRef] [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. 10.1126/science.1233521 [Abstract] [CrossRef] [Google Scholar]
  • Michelmore K., Balen A., Dunger D., Vessey M. (1999). Polycystic ovaries and associated clinical and biochemical features in young women. Clin. Endocrinol. 51 779–786. 10.1046/j.1365-2265.1999.00886.x [Abstract] [CrossRef] [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. 10.1007/BF03349222 [Abstract] [CrossRef] [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. 10.1159/000257421 [Abstract] [CrossRef] [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. 10.1038/nature18309 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Poretsky L. (1991). On the paradox of insulin-induced hyperandrogenism in insulin-resistant states. Endocr. Rev. 12 3–13. 10.1210/edrv-12-1-3 [Abstract] [CrossRef] [Google Scholar]
  • Price M. N., Dehal P. S., Arkin A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26 1641–1650. 10.1093/molbev/msp077 [Europe PMC free article] [Abstract] [CrossRef] [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. 10.1016/j.fertnstert.2003.10.004 [Abstract] [CrossRef] [Google Scholar]
  • Samantha C., Teede H. J., Harrison C. L., Joham A. E., Moran L. J., Stepto N. K. (2014). Biomarkers and insulin sensitivity in women with polycystic ovary syndrome: characteristics and predictive capacity. Clin. Endocrinol. 83 88–92. 10.1111/cen.12619 [Abstract] [CrossRef] [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. 10.1073/pnas.0808567105 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • Tremaroli V., Bäckhed F. (2012). Functional interactions between the gut microbiota and host metabolism. Nature 489 242–249. 10.1038/nature11552 [Abstract] [CrossRef] [Google Scholar]
  • Tremellen K., Pearce K. (2012). Dysbiosis of gut microbiota (DOGMA)–a novel theory for the development of polycystic ovarian syndrome. Med. Hypotheses 79 104–112. 10.1016/j.mehy.2012.04.016 [Abstract] [CrossRef] [Google Scholar]
  • Vink J. M., Sadrzadeh S., Lambalk C. B., Boomsma D. I. (2006). Heritability of polycystic ovary syndrome in a Dutch twin-family study. J. Clin. Endocrinol. Metab. 91 2100–2104. 10.1210/jc.2005-1494 [Abstract] [CrossRef] [Google Scholar]
  • Vrieze A., Holleman F., Zoetendal E. G., de Vos W. M., Hoekstra J. B., Nieuwdorp M. (2010). The environment within: how gut microbiota may influence metabolism and body composition. Diabetologia 53 606–613. 10.1007/s00125-010-1662-7 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
  • World Health Organization (2000). Obesity: Preventing and Managing the Global Epidemic. Geneva: World Health Organization; 10.1016/j.cell.2015.02.047 [Abstract] [CrossRef] [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. 10.1016/j.cell.2015.02.047 [Europe PMC free article] [Abstract] [CrossRef] [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. 10.1016/j.ebiom.2015.07.007 [Europe PMC free article] [Abstract] [CrossRef] [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 10.1371/journal.pone.0127094 [Europe PMC free article] [Abstract] [CrossRef] [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 10.1038/srep27572 [Europe PMC free article] [Abstract] [CrossRef] [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. [Abstract] [Google Scholar]

Articles from Frontiers in Microbiology are provided here courtesy of Frontiers Media SA

Citations & impact 


Impact metrics

Jump to Citations

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/20145696
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/20145696

Smart citations by scite.ai
Smart citations by scite.ai include citation statements extracted from the full text of the citing article. The number of the statements may be higher than the number of citations provided by EuropePMC if one paper cites another multiple times or lower if scite has not yet processed some of the citing articles.
Explore citation contexts and check if this article has been supported or disputed.
https://scite.ai/reports/10.3389/fmicb.2017.00324

Supporting
Mentioning
Contrasting
14
301
3

Article citations


Go to all (179) article citations

Data 


Data behind the article

This data has been text mined from the article, or deposited into data resources.