Background: Glycogen storage diseases (GSDs) are known as a group of disorders characterized by genetic errors leading to accumulation of glycogen in various tissues. Since different types of GSD can sometimes be clinically indistinguishable, next generation sequencing is becoming a powerful tool for clinical diagnosis. Methods: 12 patients with suspected GSDs and their parents were enrolled in this study. The clinical and laboratory data of the patients were reviewed. Causative gene variants were identified in the patients using whole exome sequencing (WES) and verified by Sanger sequencing. Results: Genetic testing and analysis showed that 7 patients were diagnosed with GSD II (Pompe disease), 2 patients with GSD III, 1 patient with GSD VI, and 2 patients with GSD IXα. A total number of 18 variants were identified in 12 patients including 11 variants in GAA gene, 3 variants in AGL gene, 2 variants in PYGL gene and 2 variants in PHKA2 gene, of which 9 variants were reported and 9 variants were novel. SIFT, Polyphen-2, Mutation Taster, and REVEL predicted the novel variants (except GAA c.1052_1075 + 47del) to be disease-causing. The 3D structures of wild/mutant type GAA protein were predicted indicating that variants p. Trp621Gly, p. Pro541Leu, p. Ser800Ile and p. Gly293Trp might affect the proteins function via destroying hydrogen bonds or conformational constraints. Neither liver size nor laboratory findings allow for a differentiation among GSD III, GSD VI and GSD IXα. Conclusion: Our study expanded the variation spectrum of genes associated with GSDs. WES, in combination with clinical, biochemical, and pathological hallmarks, could provide accurate results for diagnosing and sub-typing GSD and related diseases in clinical setting.