It was a prospective interventional study conducted at a tertiary care centre over a period of 2 years (January 2019 to December 2020). Approval was obtained from the institutional ethical committee vide letter no IEC/2019/67/MA-06. The study was conducted in accordance with the ethical standards of the institutional and national research committee and with the 1964 Helsinki Declaration and its later amendments. Children between the ages of 2 and 18 years with CHB who had received antivirals for at least 24 months with (i) alanine aminotransferase (ALT) levels persistently < 1.5 times the upper limit of normal (ULN) (3 consecutive readings done 6 months apart in last 12 months); (ii) hepatitis B virus – deoxyribonucleic acid (HBV DNA) not detected on 3 consecutive occasions 6 months apart in last 12 months and (iii) whose parents/ guardians gave consent for the study, were included. Normal ALT levels were defined as per the cut-offs obtained from SAFETY study [11]. Children with family history of cirrhosis, HCC or co-infection with hepatitis C virus and human immuno-deficiency virus were excluded. Liver Biopsy was done to evaluate the histological activity index (HAI) and fibrosis as per Ishak staging [12]. Children with HAI < 5 and fibrosis < 3 in liver histology were offered the option of stopping antiviral therapy. Additionally, quantitative estimation of covalently closed circular DNA (ccc-DNA) was also done in the liver tissue. ccc-DNA levels were measured by real-time polymerase chain reaction (PCR) based in-house assay on lightcycler 480 (LC480) (Roche Diagnostics, Netherlands) instrument using following primers: Forward primer: 5’-CTCCCCGTCTGTGCCTTCT-3’; Reverse Primer: 5’-GCCCCAAAGCCACCCAAG-3’. This is a qualitative real-time PCR based assay using sybergreen probes for the target detection [13]. HBV-DNA was measured in the plasma samples by using automated real-time PCR kit, COBAS® AmpliPrep/COBAS® TaqMan®, v2.0 (Roche Diagnostics, GmbH, Mannheim, Germany). This is an automated hepatitis B viral load quantitative assay needing 650 µL of serum or plasma, with the linear range of quantification being 1.3-8.23 log10 IU/ml in all of them. HBsAg quantification was done by chemiluminiscent immunoassay (CLIA) method (Abbott Laboratories, Chicago, IL, USA) [14].
After stopping the antivirals, the patient was followed up at 1 month, 3 months and every 3 months thereafter. Workup at each follow up included: ALT, quantitative HBV-DNA (qHBV-DNA), quantitative HBsAg (qHBsAg), HBeAg, alfa-fetoprotein (AFP) and ultrasonography of the abdomen. Virological relapse was defined as any elevation of qHBV-DNA > 2000 IU/mL. Biochemical relapse was defined as elevation of ALT > 2 times ULN. Primary outcome of the study was the incidence of relapse after stopping antiviral therapy. Secondary outcome was identification of the predictors of relapse after stopping antiviral therapy. The predictors studied included clinical parameters, type of therapy, HBeAg seroconversion status, time taken for HBV-DNA to become undetected on therapy, HBV genotype, ccc-DNA levels, duration of normal ALT before stopping therapy and qHBsAg at the time of stopping antivirals. Those with only VR were continued on regular follow up whereas those with VR as well as BR were treated with combination therapy: 8 weeks of antivirals (entecavir in those younger than 12 years; tenofovir in those older than 12 years), followed by combination of peg-IFN alpha 2b (60 mcg/m2/week) for 48 weeks. The study design is depicted in figure 1.
Statistical analysis: All analyses were performed with SPSS 22.0 (SPSS, Chicago, IL). As this was a pilot study, all the children fulfilling the inclusion criteria during the study period were enrolled. Continuous variables were expressed as median and inter-quartile range (IQR), while categorical variables were expressed as proportions. Both categorical and continuous variables were analyzed using cox regression analysis. Cox regression model was used for univariate and multivariate analyses for predictors of relapse. Receiver operating characteristic curves were used to calculate the area under the receiver operating characteristic curve (AUROC).