In Tunisia, the prevalence of naturally occurring surface (S) gene variants of hepatitis B virus (HBV) has not been determined. In the present study, the prevalence of these variants was examined in terms of the clinical and viral state in a series of 99 Tunisian patients with HBV infection. The S genes were amplified and directly sequenced. Genotype D was predominant (98%), 40.4% isolates belonged to subgenotypes D7 and 1 to subgenotype D2. The most common subtype was ayw2 (95.9%). In total, 60.6% of the studied strains harbored S mutations. Several novel mutation patterns were detected. Interestingly, the presence of S mutations was significantly correlated with the D7 subgenotype, low HBV DNA and advancing age (≥35 years), and tended to be higher in liver cirrhosis than in chronic infection. The global prevalence of the major hydrophilic region variants was 12.1%, with substitution S143L/T as the most frequent (4%). Only 33.9% of S substitutions produced amino acid changes in the polymerase gene. In conclusion, a high prevalence of naturally occurring HBsAg variants was observed among Tunisian HBV carriers. Natural viral variability in a geographical region and duration of infection are among the major factors associated with the occurrence of S mutations.