Polycystic ovary syndrome (PCOS) is a heterogeneous endocrine disorder, which affects 4-10 % women of reproductive age. Though accumulating scientific evidence, its pathogenesis remains unclear. In the current study, metabolic profiling as well as diagnostic biomarkers for different phenotypes of PCOS was investigated using non-invasive urinary GCMS based metabolomics. A total of 371 subjects were recruited for the study. They constituted the following groups: healthy women, those with hyperandrogenism (HA), women with insulin-resistance (IR) in PCOS. Two cross-comparisons with PCOS were performed to characterize metabolic disturbances. A total of 23 differential metabolites were found. The altered metabolic pathways included glyoxylate and dicarboxylate metabolism, pentose and glucuronate interconversions, and citrate cycle and butanoate metabolism. For differential diagnosis, a panel consisting of 9 biomarkers was found from the comparison of PCOS from healthy subjects. The area under the receiver operating characteristic (ROC) curve (AUC) was 0.8461 in the discovery phase. Predictive value of 89.17 % was found in the validation set. Besides, a panel of 8 biomarkers was discovered from PCOS with HA vs IR. The AUC for 8-biomarker panel was 0.8363, and a panel of clinical markers (homeostasis model assessment-insulin resistance and free androgen index) had 0.8327 in AUC. While these metabolites combined with clinical markers reached 0.9065 in AUC from the discovery phase, and 93.18 % in predictive value from the validation set. The result showed that differences of small-molecule metabolites in urine may reflect underlying pathogenesis of PCOS and serve as biomarkers for complementary diagnosis of the different phenotypes of PCOS.