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.