Studies characteristics
Our initial search identified 349 potentially relevant publications, 36 of which were retrieved for detailed review. Twenty six studies met our inclusion criteria (Supplementary Fig. S1) [19-44]. The main characteristics of the included studies are presented in Tables 1 and 2. The sample size of these studies ranged 6–1075 participants, and there were a total of 3451 CHB patients. Among these patients, 2573 (74.6%) ceased NAs. Of the 2573 patients who ceased NAs, 1899 (73.8%) were male. At baseline, 337 (13.9%) of the 2433 patients with cirrhosis (information of cirrhosis was not reported in two studies [24,27], including 140 off-NAs patients), which was diagnosed by histology findings and/or imaging studies.
Sixteen studies included Asian patients [19,22-28,30,33,35,38-41,44], six studies involved Caucasian patients [20,21,29,32,36,42], and four studies included mixed patients population [31,34,37,43]. Lamivudine was used in one [32], entecavir in three [23,24,28], tenofovir in two [20,34], several NAs in fifteen [19,21,22,26,30,31, 33,35-41,44], and NAs was not reported in the remaining five studies [25,27,29,42,43] (Table 1).
Only HBeAg-positive patients were involved in three studies [25,26,43], only HBeAg-negative patients in ten [20,29,31,32,34-38,42], and both HBeAg-positive and HBeAg-negative patients in thirteen [19,21,23,24,26-28,30,33,39-41,44]. In total, 747 (29.0%) HBeAg-positive and 1826 (71.0%) HBeAg-negative patients were included (Table 1).
Great heterogeneity was found among the 26 studies in respect of age, serum alanine aminotransferase (ALT) levels, HBV DNA levels, HBsAg levels, duration of virological response before NAs cessation, duration of HBeAg seroconversion before NAs cessation, treatment duration and duration of follow up after NAs cessation (Table 2). The shortest duration of follow up after NAs cessation was 12 months, and the longest was 62 months.
The definition of VR after NAs cessation was based on HBV DNA levels >20 IU/mL in one [22], >120 IU/mL in one [21], >2000 IU/mL (or 10,000 copies/mL) in 19 [19,20,22-26,28-31,33-35,38-40,43,44], and >20,000 IU/mL (or 100,000 copies/mL) in one study [27]. One study provided two definitions of VR which were based on HBV DNA levels >2000 IU/mL and >20,000 IU/mL [32]. Two studies provided three definitions of VR which were based on HBV DNA levels >200 IU/mL, >2000 IU/mL and >20,000 IU/mL [36,37]. While it was not specified in one study [42]. When calculating the overall pooled rates of VR and VR in view of the duration of follow up after NAs cessation, we adopted the definition of VR based on HBV DNA level >2000 IU/mL, therefore, four studies [21,22,27,42] were not included.
BR after NAs cessation was defined as ALT above the upper limit of normal (ULN) in three [24,31,41], >1.2ⅹULN in one [20], and >2ⅹULN in fifteen studies [19,22,25-29,32-35,38,39,43,44]. One study specified two definitions of BR which were based on ALT >ULN and > 2ⅹULN [36]. One study specified five definitions of BR which were based on ALT >ULN, > 2ⅹULN, > 3ⅹULN, > 5ⅹULN and > 10ⅹULN [37]. While it was not specified in five studies [21,22,30,40,42]. When calculating the overall pooled rates of BR and BR in view of the duration of follow up after NAs cessation, we adopted the definition of BR based on ALT level >2×ULN, therefore, nine studies [20-22,24,30,31,40-42] were not included.
VR after NAs cessation
VR was observed in 1686 of 2459 patients (random effects pooled estimate = 0.64, 95% CI 0.57-0.72; P for heterogeneity < 0.001) despite of the duration of follow up after NAs cessation (Table 3). The overall pooled rate of VR in initially HBeAg-positive patients (237/456; random effects = 0.62, 95% CI 0.41-0.82; P for heterogeneity < 0.001) was lower than initially HBeAg-negative patients (975/1375; random effects = 0.63, 95% CI 0.54-0.71; P for heterogeneity < 0.001). However, there was no statistic difference between initially HBeAg-positive patients and HBeAg-negative patients (P = 0.731). The forest plots of the probabilities of VR in all, initially HBeAg-positive, and HBeAg-negative patients are presented in Fig. 1A-C.
In view of the duration of follow up after NAs cessation, the pooled rates (95% CI) of VR were 0.48 (0.39-0.58), 0.56 (0.47-0.65), 0.61 (0.50-0.71), 0.53 (0.33-0.70), 0.71 (0.62-0.82) and 0.65 (0.43-0.86) at 6 (13 studies included), 12 (17 studies included), 24 (13 studies included), 36 (12 studies included), 48 (9 studies included) and 60 (7 studies included) months after NAs cessation, respectively (Fig. 2A). The 6-month, 12-month, 24-month, 36-month, 48-month and 60-month pooled rates (95% CI) of VR were lower, but not significantly, in initially HBeAg-positive patients (0.36 [0.28-0.44], 0.50 [0.32-0.65], 0.58 [0.42-0.71], 0.49 [0.38-0.60], 0.50 [0.43-0.59], 0.60 [0.42-0.62]) than HBeAg-negative patients (0.51 [0.40-0.62], 0.56 [0.45-0.65], 0.67 [0.48-0.85], 0.59 [0.39-0.83], 0.53 [0.41-0.69], 0.65 [0.52-0.76]) (P = 0.400) (Fig. 2B-C).
In all patients, the pooled probabilities of VR after NAs cessation were getting lower and lower in studies defining VR by HBV DNA >200, >2000 and >20,000 IU/mL (VR [95% CI] at 12 months after NAs cessation: 0.80 [0.72-0.87] versus 0.57 [0.47-0.67] versus 0.39 [0.29-0.48]); Odds ratio [95% CI]: 0.89 [0.03-4.8] and 0.25 [0.12-3.25], P = 0.006) (Table 4). Likewise, the pooled rates of off-NAs VR were significantly different according to the VR definition in initially HBeAg-negative patients (P < 0.001). There were too few studies of initially HBeAg-positive patients to calculate the pooled probabilities of VR according to the different VR definitions.
BR after NAs cessation
BR was reported in 843 of 1778 patients (random effects pooled estimate = 0.45, 95% CI 0.36-0.53; P for heterogeneity < 0.001) despite of the duration of follow up after NAs cessation (Table 3). The overall pooled rate of BR in initially HBeAg-positive patients (153/407; random effects = 0.43, 95% CI 0.24-0.62; P for heterogeneity < 0.001) was lower than initially HBeAg-negative patients (622/1200; random effects = 0.51, 95% CI 0.41-0.61; P for heterogeneity < 0.001). However, there was no statistic difference between HBeAg-positive and HBeAg-negative patients (P = 0.572). The forest plots of the probabilities of BR in all, initially HBeAg-positive, and HBeAg-negative patients are presented in Supplementary Fig. S2A-C.
In view of the duration of follow up after NAs cessation, the pooled rates (95% CI) of BR were 0.30 (0.15-0.42), 0.35 (0.28-0.41), 0.44 (0.31-0.58), 0.45 (0.36-0.56), 0.49 (0.36-0.61) and 0.48 (0.37-0.61) at 6 (10 studies included), 12 (14 studies included), 24 (10 studies included), 36 (9 studies included), 48 (6 studies included) and 60 (5 studies included) months after NAs cessation, respectively (Supporting Fig. S3A). The 6-month, 12-month, 24-month, 36-month, 48-month and 60-month pooled rates (95% CI) of BR were lower, but not significantly, in initially HBeAg-positive patients (0.15 [0.07-0.22], 0.36 [0.22-0.50], 0.39 [0.32-0.45], 0.44 [0.32-0.56], 0.47 [0.24-0.70], 0.40 [0.19-0.62]) than HBeAg-negative patients (0.29 [0.10-0.45], 0.40 [0.22-0.58], 0.52 [0.26-0.73], 0.53 [0.32-0.74], 0.57 [0.37-0.78], 0.52 [0.35-0.73]) (P = 0.217) (Supporting Fig. S3B-C).
In all patients, the pooled probabilities of BR after NAs cessation were getting lower and lower in studies defining BR by ALT >1×ULN, >2×ULN and >5×ULN (BR [95% CI] at 12 months after NAs cessation: 0.48 [0.13-0.95] versus 0.35 [0.28-0.81]) versus 0.19 [0.02-0.58]); Odds ratio [95% CI]: 0.38 [0.19-1.4] and 0.34 [0.11-2.84], P = 0.013) (data not shown). There were too few studies of initially HBeAg-positive and HBeAg-negative patients to calculate the pooled probabilities of BR according to the different BR definitions.
Retreatment after NAs cessation
In ten studies supplying such data [19,24,26,31,33-35,37,39,42], retreatment was started in 447 (84%) of 534 patients with BR or 447 (52%) of 859 patients with VR.
Durability of HBeAg seroconversion after NAs cessation
As only three studies provided data about the durability of HBeAg seroconversion after NAs cessation [20,21,43], we did not calculate the pooled estimate.
HBsAg loss or seroconversion after NAs cessation
After NAs cessation, HBsAg loss was reported in 144 of 1721 patients (random effects pooled estimate = 0.09, 95% CI 0.06-0.12; P for heterogeneity < 0.001). There was no significant difference between initially HBeAg-positive patients (14/209; random effects = 0.06, 95% CI 0.02-0.10; P for heterogeneity = 0.222) and HBeAg-negative patients (99/1048; random effects = 0.14, 95% CI 0.08-0.20; P for heterogeneity < 0.001) (P = 0.131) (Table 3).
HBsAg seroconversion was reported in 65 of 1036 patients (random effects pooled estimate = 0.06, 95% CI 0.05-0.11; P for heterogeneity = 0.017) (Table 3).
Factors associated with VR after NAs cessation
Results of subgroup meta-analysis found that the duration of on-therapy virological response had a strong effect in the rates of VR after NAs cessation in all, initially HBeAg-positive patients or HBeAg-negative patients (Table 4). Especially, the pooled rates (95% CI) of VR at 12 months after NAs cessation were significant different between duration of on-NAs virological response ≤24 months and >24 months in all patients (0.55 [0.39-0.72] versus 0.41 [0.13-0.66], OR [95% CI] = 0.89 [0.63-3.56], P = 0.011), initially HBeAg-positive patients (0.53 [0.32-0.73] versus 0.37 [0.16-0.68], OR [95% CI] = 0.67 [0.01-3.17], P = 0.031), and initially HBeAg-negative patients (0.59 [0.21-0.84] versus 0.41 [0.21-0.53], OR [95% CI] = 0.56 [0.14-4.26], P = 0.025).
Results of subgroup meta-analysis found that the duration of treatment did not affect the probability of VR after NAs cessation in all or in initially HBeAg-negative patients. There was too few data from studies of initially HBeAg-positive patients to do the subgroup meta-analysis (Table 4).
In spite of 14 studies provided data about predictors of VR after NAs cessation [19,22-25,28,32,33,36,37,39-41,43], they often reported conflicting results. Off-NAs VR was found to be associated with higher baseline ALT and platelet in one [24], higher baseline gamma-glutamyltransferase (GGT) in one [32], high baseline HBV DNA (>106copies/ml) in one [24], baseline advanced fibrosis in one [36], and higher end of treatment (EOT) ALT in one study [23], but not in any other study included in this review. Male gender was associated with off-NAs VR only in one study [24]. In addition, older age was associated with off-NAs VR in seven studies (all with both HBeAg-positive and HBeAg-negative patients) [19,23,33,39-41,43], and two of these studies found that higher rates of off-NAs VR in patients >35 years old than patients ≤35 years old [33,44]. They did not found off-NAs VR was association with alcohol consumption in one [36], diabetes and hypertension in one [28], and body mass index (BMI) in two studies [29,36].
The rate of VR after NAs cessation was not related to the type of NAs in six studies using first-line (tenofovir or entecavir) and second-line NAs (adefovir, lamivudine or telbivudine) involving 683 patients (330 HBeAg-positive, 353 HBeAg-negative) [33,36,37,39,40,44]. Whereas, the rate of VR was higher in 22 tenofovir-treated than 113 entecavir-treated HBeAg-negative CHB patients in one study [41]. The rate of VR after NAs cessation was not related to HBV genotype (genotypes B, C and D) only in one study with HBeAg-negative patients reporting such data [37]. No relationship was found between the risk of off-NAs VR and EOT liver stiffness in three studies [19,32,33].
Higher EOT HBsAg levels were related to VR after NAs cessation in six studies with both HBeAg-positive and HBeAg-negative patients (269 HBeAg-positive, 327 HBeAg-negative) [19,23,28,39,41,44], but not in one study with 18 HBeAg-positive patients [22], two studies with 80 HBeAg-negative patiens [32,36], and in one study with both HBeAg-positive and HBeAg-negative patients (60 HBeAg-positive, 22 HBeAg-negative) [33]. Total serum level of EOT antibody against the HBV core protein (anti-HBc) [39] and EOT hepatitis B core-related antigen (HBcrAg) [41] were associated with risk of VR after NAs cessation in each one study.
When compared with initially HBeAg-positive and HBeAg-negative patients, only one study suggested initially HBeAg-positive patients with higher off-NAs VR [22], whereas no difference was found in other seven studies [23,25,28,33,39,41,44]. In the studies including HBeAg-positive patients at the treatment onset, the rate of VR after NAs cessation was associated with time to HBeAg seroconversion in one study with 138 patients [40], but not in another study with 58 patients [25]. Moreover, the rates of off-NAs VR were affected by consolidation duration in two studies including 196 patients [25,40], but not in another three studies including 152 patients [19,23,39].
Clinical outcomes after NAs cessation
Data of clinical outcomes after NAs cessation was provided in 11 studies including 1365 CHB patients [21,26,27,29,32-36,39,40].
During follow-up after NAs cessation, 7 of 337 patients with cirrhosis at baseline suffered from hepatic decompensation, all of whom were initially HBeAg-negative patients. They were retreated promptly but three patients died at 3, 15 and 32 months after NAs cessation, respectively.
Among 1166 patients, twenty-three patients developed HCC after NAs cessation, in whom, 19 patients were cirrhosis at baseline, and 2 patients were non-cirrhosis. Moreover, all of them were initially HBeAg-negative patients.
Sensitivity and publication bias
Sensitivity analysis was performed by sequential omission of individual studies. For analyses that included more than three pooled individual studies, the significance of the VR was not significantly altered by omitting any single study (data not shown).
Among 26 studies that estimated the rates of VR, no evidence of publication bias was found in the main analysis (p > 0.05) (Supplementary Fig. S4).
Study quality
There were only one randomized study [34] and the other 25 cohort studies including five retrospective [24-26,31,43], three retrospective-prospective [19,32,35] and 17 prospective [20-23,27-30,33,36-42,44]. According to the risk of bias assessments, four studies had high quality [19,22,23,32], nine had acceptable quality [25,28,29,32,34,36,39,40,44], and 13 had low quality [20,21,24,26,27,30,32,35,37,38, 41-43] (Supplementary Table S1).