This study was aimed to construct classification and regression tree (CART) model of glycosaminoglycans (GAGs) for the differential diagnosis of Mucopolysaccharidoses (MPS). Two-dimensional electrophoresis and liquid chromatography-tandem mass spectrometry (LC-MS/MS) were used for the qualitative and quantitative analysis of GAGs. Specific enzyme assays and targeted gene sequencing were performed to confirm the diagnosis. Machine learning tools were used to develop CART model based on GAG profile. Qualitative and quantitative CART models showed 96.3% and 98.3% accuracy, respectively, in the differential diagnosis of MPS. The thresholds of different GAGs diagnostic of specific MPS types were established. In 60 MPS positive cases, 46 different mutations were identified in six specific genes. Among 31 different mutations identified in IDUA, nine were nonsense mutations and two were gross deletions while the remaining were missense mutations. In IDS gene, four missense, two frameshift, and one deletion were identified. In NAGLU gene, c.1693C > T and c.1914_1914insT were the most common mutations. Two ARSB, one case each of SGSH and GALNS mutations were observed. LC-MS/MS-based GAG pattern showed higher accuracy in the differential diagnosis of MPS. The mutation spectrum of MPS, specifically in IDUA and IDS genes, is highly heterogeneous among the cases studied.