Hereditary spherocytosis (HS) is an inherited disorder characterized by anemia, splenomegaly, and spherical-shaped erythrocytes, caused by mutations in erythrocyte membrane Protein Genes (ANK1, SPTB, SLC4A1, SPTA1, and EPB42). We investigated molecular spectrum and genotype-phenotype correlation in HS patients in Hubei province, central China. Twenty-three patients with HS were included. A next-generation sequencing (NGS) panel targeting ANK1, SPTB, SLC4A1, SPTA1, and EPB42 genes was used to screen potential variants. Sanger sequencing was applied to validate variants. Of the twenty-three patients, thirteen patients carried ANK1 variants, and ten patients harbored SPTB variants, including ten non-sense, six indel, four splice site, one start-loss, and one missense variant. Four out of twenty-two variants in our study were known, and eighteen variants were novel. Most ANK1 and SPTB variants were indel (5/12) or non-sense (7/10), respectively. Family member analysis in thirteen families showed that six variants were de novo. Variable expressivities were observed in a pair of twins with ANK1 c.341C > T variant, and two unrelated patients both carried ANK1 c.2T > A variant. Genotype-phenotype analysis found no significant difference between ANK1 and SPTB regarding the levels of Hb, RBC, MCV, MCH, and MCHC. However, variants in the ANK1 death domain were associated with lower levels of MCV and MCH compared to other ANK1 domains. In conclusion, NGS is a fast way to provide a molecular HS diagnosis. We also identified unique genetic and clinical characteristics of patients with HS in Hubei Province, China. However, a large sample size is needed to further investigate the genotype-phenotype correlation.