Some evidence suggest that bone health can be regulated by gut microbiota. To better understand this, we performed 16S ribosomal RNA sequencing to analyze the intestinal microbial diversity in primary osteoporosis (OP) patients, osteopenia (ON) patients and normal controls (NC). We observed an inverse correlation between the number of bacterial taxa and the value of bone mineral density. The diversity estimators in OP and ON groups were increased compared with that in NC group. Beta diversity analyses based on hierarchical clustering and principal coordinate analysis (PCoA) could discriminate the NC samples from OP and ON samples. Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria constituted the four dominant phyla in all samples. Proportion of Firmicutes was significantly higher and Bacteroidetes was significantly lower in OP samples than that in NC samples (p < 0.05), Gemmatimonadetes and Chloroflexi were significantly different between OP and NC group as well as between ON and NC group (p < 0.01). A total of 21 genera with proportions above 1% were detected and Bacteroides accounted for the largest proportion in all samples. The Blautia, Parabacteroides and Ruminococcaceae genera differed significantly between the OP and NC group (p < 0.05). Linear discriminant analysis (LDA) results showed 1 phylum community and 7 phylum communities were enriched in ON and OP, respectively. 35 genus communities, 5 genus communities and 2 genus communities were enriched in OP, ON and NC, respectively. The results of this study indicate that gut microbiota may be a critical factor in osteoporosis development, which can further help us search for novel biomarkers of gut microbiota in OP and understand the interaction between gut microbiota and bone health.