To elucidate mutation spectrum and genotype-phenotype correlations in Japanese patients with OI, we conducted comprehensive genetic analyses using NGS, as this had not been analyzed comprehensively in this patient population. Most mutations were located on COL1A1 and COL1A2. Glycine substitutions in COL1A1 resulted in the severe phenotype. Introduction Most cases of osteogenesis imperfecta (OI) are caused by mutations in COL1A1 or COL1A2, which encode α chains of type I collagen. However, mutations in at least 16 other genes also cause OI. The mutation spectrum in Japanese patients with OI has not been comprehensively analyzed, as it is difficult to identify using classical Sanger sequencing. In this study, we aimed to reveal the mutation spectrum and genotype-phenotype correlations in Japanese patients with OI using next-generation sequencing (NGS). Methods We designed a capture panel for sequencing 15 candidate OI genes and 19 candidate genes that are associated with bone fragility or Wnt signaling. Using NGS, we examined 53 Japanese patients with OI from unrelated families. Results Pathogenic mutations were detected in 43 out of 53 individuals. All mutations were heterozygous. Among the 43 individuals, 40 variants were identified including 15 novel mutations. We found these mutations in COL1A1 (n = 30, 69.8%), COL1A2 (n = 12, 27.9%), and IFITM5 (n = 1, 2.3%). Patients with glycine substitution on COL1A1 had a higher frequency of fractures and were more severely short-statured. Although no significant genotype-phenotype correlation was observed for bone mineral density, the trabecular bone score was significantly lower in patients with glycine substitutions. Conclusion We identified pathogenic mutations in 81% of our Japanese patients with OI. Most mutations were located on COL1A1 and COL1A2. This study revealed that glycine substitutions on COL1A1 resulted in the severe phenotype among Japanese patients with OI.