Global impact of hepatitis B vaccination on the risk of hepatocellular carcinoma: a systematic review and meta-analysis
Highlight box
Key findings
• This systematic review and meta-analysis demonstrates that hepatitis B virus (HBV) vaccination is associated with a substantial and consistent reduction in hepatocellular carcinoma (HCC) incidence worldwide. Vaccinated individuals experienced a 65–67% lower risk of developing HCC, with protective effects observed across sexes, geographic regions, vaccine formulations, and dosing regimens.
What is known and what is new?
• Chronic HBV infection is a major risk factor for HCC, and HBV vaccination has long been recognized as effective in preventing HBV infection and its complications.
• This study provides comprehensive, global, and quantitatively robust evidence confirming that HBV vaccination confers durable protection against HCC across diverse populations. It demonstrates comparable effectiveness between plasma-derived and recombinant vaccines and across different dosing schedules using both unadjusted and adjusted risk estimates.
What is the implication, and what should change now?
• These findings reinforce HBV vaccination as a highly effective population-level cancer prevention strategy. Strengthening universal HBV immunization programs, improving vaccine coverage, and ensuring completion of recommended dosing regimens should remain global public health priorities to reduce the burden of HCC.
Introduction
Hepatocellular carcinoma (HCC) ranks among the most lethal malignancies globally, ranking as the third leading cause of cancer mortality and the sixth most commonly diagnosed cancer in 2020, with an estimated 906,000 new cases per year. Universal risk factors include chronic hepatitis B virus (HBV) infection, which contributes substantially to the global HCC burden, especially in endemic settings such as East Asia and sub-Saharan Africa (1-3).
Chronic HBV infection drives hepatocarcinogenesis via multifaceted mechanisms including ongoing inflammation, hepatocyte turnover, fibrosis, integration of viral DNA into the host genome, and eventual clonal expansion of transformed hepatocytes (2,4,5). Because HBV can integrate even in the absence of high-level replication, preventing initial infection is arguably the most effective strategy to block this cascade at its outset (1,2,5).
The advent of HBV vaccination in the 1980s marked a paradigmatic shift in liver disease prevention. The hepatitis B vaccine is widely recognized as the first vaccine with documented anticancer impact, given its success in preventing chronic HBV infection and, by extension, reducing HCC risk (1,2,6). Over four decades of implementation, global coverage of the three-dose HBV vaccination series has expanded into national immunization schedules of more than 190 countries (1,6). In many settings, inclusion of a timely birth dose has helped block mother-to-child transmission, a major route of HBV acquisition in high-prevalence areas (1,6).
Empirical observational and ecological studies have documented declines in HBsAg prevalence, chronic HBV carriage, cirrhosis, and HCC incidence following vaccine introduction. For example, universal childhood vaccination in Taiwan led to >90% reductions in HBsAg prevalence among vaccinated cohorts and progressively lower rates of childhood HCC over time (1,2,5). In Shanghai, joint-point regression analyses revealed accelerating decreases in HCC, particularly among younger male cohorts following implementation of universal HBV immunization (7). Nonetheless, interpretations of those trends may be complicated by secular changes in diagnostics, exposure to other risk factors (e.g., aflatoxin, alcohol, metabolic disease), and varying lag times between infection and cancer onset.
Despite the weight of these findings, prior evidence summaries on this topic have largely been narrative or descriptive in nature (1,2,6,8), with a paucity of formal quantitative syntheses pooling adjusted estimates of HCC across diverse populations. critical questions remain unresolved. First, no prior synthesis has robustly quantified incidence rate ratios (IRRs) or adjusted relative risks (aRRs) linking vaccination status to HCC across varied populations and vaccine formulations. Second, the durability and magnitude of protection—especially into adolescence and early adulthood—are still incompletely characterized. Third, differences by sex, vaccine type (plasma-derived vs. recombinant), and incomplete dosing regimens remain poorly understood at the level of cancer endpoints. Lastly, because many earlier studies report only ecological trends or unadjusted incidence declines, there is a need for rigorous meta-analysis combining both incidence-based and adjusted comparative data to strengthen causal inference.
To address these gaps, we conducted a systematic review and meta-analysis to estimate the pooled effect of HBV vaccination on incident HCC across global settings, synthesizing both occurrence-based IRRs and multivariable-adjusted risks. We performed subgroup analyses by sex, geographical region, vaccine formulation, and dosing regimen to explore heterogeneity and identify factors modifying vaccine effectiveness. Through this work, we aim to provide the most comprehensive quantitative estimate to date of how HBV vaccination contributes to reducing HCC incidence, thereby reinforcing the evidence base for vaccination as an oncologic prevention tool. We present this article in accordance with the PRISMA reporting checklist (available at https://jgo.amegroups.com/article/view/10.21037/jgo-2026-0202/rc) (9).
Methods
Study design and registration
This study was conducted as a systematic review and meta-analysis. The methodology followed a pre-defined protocol developed before data collection, and the final protocol was registered on PROSPERO (ID: 1355073). The research question was framed using the Population-Exposure-Outcome framework (10): among populations with exposure to the HBV vaccine, what is the risk of developing HCC compared with those who were unvaccinated or inadequately vaccinated?
Eligibility criteria
Studies were eligible for inclusion if they met the following criteria: (I) participants had received any formulation of the HBV vaccine (plasma-derived or recombinant), either as part of a national immunization program or within a clinical or observational setting; (II) the study reported HCC as an outcome, expressed in terms of incidence rate, prevalence rate, or RR metrics (e.g., risk ratio, odds ratio, IRR, or hazard ratio); (III) the number of HCC cases exceeded 20 (this threshold was adopted to ensure sufficient event counts for stable effect estimates and reliable variance calculation, given that HCC is a rare outcome and very small event counts produce unstable risk estimates and inflated standard errors when pooled); and (IV) the design was observational (cohort, case-control, or single-arm population study) or randomized controlled trial with extractable comparative or single-arm data. Studies reporting only mortality outcomes, reviews, editorials, conference abstracts, and non-English publications were excluded. When multiple reports described overlapping cohorts, the study with the most comprehensive or updated data was retained.
Search strategy
A comprehensive search of the electronic databases PubMed, Scopus, and Web of Science was conducted on July 24, 2025. The search was designed to identify all relevant studies evaluating the relationship between HBV vaccination and HCC incidence. The PubMed search strategy was developed using a combination of Medical Subject Headings (MeSH) and free-text terms related to “hepatitis B vaccine”, “hepatocellular carcinoma”, and risk metrics, and was adapted for use in the other databases. The full search syntax can be found in Table S1. All references were then imported into EndNote X8 (Clarivate Analytics) for duplicate removal and management.
Study selection
After removal of duplicate records, two reviewers independently screened all remaining titles and abstracts for relevance. Full-text articles of potentially eligible studies were then retrieved and assessed against the inclusion and exclusion criteria. Discrepancies were resolved by consensus or consultation with the senior reviewer.
Data extraction
Two investigators independently extracted data from each included study using a standardized data collection form. Extracted variables included first author, publication year, study design, country/region, study period, total sample size, number of vaccinated and unvaccinated participants, age range, sex distribution, diagnostic criteria for HCC, vaccine formulation, number of doses, and follow-up duration. Where reported, both unadjusted and age- and gender-adjusted effect estimates (IRRs, risk ratios, or hazard ratios) with corresponding 95% confidence intervals (CIs) were extracted. For studies reporting only incidence rates and person-years, IRRs and their standard errors were calculated using standard Poisson-based formula (11). When overlapping or repeated cohorts were identified, the most recent or complete dataset was prioritized.
Quality assessment
Methodological quality was evaluated using the National Institutes of Health (NIH) Quality Assessment Tool for Observational Cohort and Cross-Sectional Studies (12). Each study was assessed across 14 domains, including population definition, exposure measurement, confounding control, and outcome ascertainment. Each item was scored 0 (no), 1 (cannot determine/not reported/not applicable), or 2 (yes), giving a maximum possible total of 28. Scores ≥21 were considered good quality, 14–20 fair quality, and <14 poor quality. As for randomized trials, the revised Cochrane risk of bias (Rob-2) tool was used. Disagreements were resolved through consensus.
Statistical analysis
All quantitative analyses were performed using Stata version 18.0 (StataCorp, College Station, TX, USA), with confirmatory re-analysis of the adjusted pooled estimates performed in R (version 4.3.2; meta package) to verify the symmetry of 95% CIs on the natural log scale. For each study, log-transformed IRRs or RRs and their corresponding standard errors were derived from the reported estimates or calculated from event counts and person-years. Unadjusted IRRs and aRRs were analyzed in two separate, parallel meta-analyses rather than synthesized together. No hazard ratios or odds ratios were extracted or pooled. This separation preserves the conceptual distinction between incidence-rate comparisons (which inherently account for person-time) and adjusted risk-based comparisons, and avoids any implicit assumption of equivalence between effect-measure types. Pooled estimates were generated using the DerSimonian-Laird random-effects model, which remains the most widely applied estimator in meta-analyses and is the default in standard meta-analytic software; consistency with restricted maximum likelihood (REML) estimation was confirmed during R-based re-analysis (13), which account for between-study heterogeneity. Heterogeneity was quantified using the I2 statistic, with values of 25%, 50%, and 75% representing low, moderate, and high heterogeneity, respectively.
Subgroup analyses were conducted by country/region, sex, age group, vaccine type (plasma-derived versus recombinant), and number of doses (two versus three or more). When available, multivariable-adjusted estimates (aRRs) were synthesized separately from unadjusted rate ratios to evaluate robustness against confounding. Publication bias was assessed visually through funnel plots and, when ≥10 studies were available, quantitatively using Egger’s regression test. All statistical tests were two-sided, with P values <0.05 considered statistically significant.
Results
Study selection
A total of 1,365 records were retrieved through PubMed, Scopus, and Web of Science searches. After removal of 492 duplicates in EndNote, 873 unique records were screened by title and abstract. Of these, 791 were excluded for irrelevance. Full texts of 82 articles were sought for retrieval; 79 were successfully obtained and assessed. Ultimately, 16 studies met the eligibility criteria and were included in the quantitative synthesis (Figure 1) (7,14-28). Excluded reports primarily focused on HCC-related mortality rather than incidence (n=2), lacked raw or extractable data (n=1), examined non-relevant populations (n=32; comprising studies in patients with established chronic HBV infection or pre-existing cirrhosis evaluating treatment outcomes, hepatitis B immunoglobulin recipients without active vaccination, occupational sub-cohorts unrepresentative of the general population, and animal or in vitro studies), did not include HBV vaccination data (n=14), had no reportable outcomes (n=5), used ineligible study designs (n=8), or were non-English (n=1).
Study characteristics
Key characteristics of the included studies are summarized in Table 1. The 16 studies were published between 1997 and 2022 and collectively encompassed more than 5 million participants across following countries/regions: Taiwan (n=8), mainland of China (n=5), South Korea (n=1), The Gambia (n=1), and the United States (n=1); the U.S. study comprised an Alaska Native cohort with a high baseline prevalence of HBV genotype F, a unique subset not representative of the broader U.S. population. Most investigations were retrospective cohort studies (n=12), accompanied by two prospective cohorts and two randomized controlled trials. Observation periods spanned 1979 to 2018, covering populations born before and after the introduction of national HBV vaccination programs.
Table 1
| Study | Country/region | Study design | Year of investigation | Sample | Age (years) | Male, n (%) | Liver cancer diagnostic method | HBV vaccine | ||
|---|---|---|---|---|---|---|---|---|---|---|
| Total | Vaccinated | Non-vaccinated | ||||||||
| Chang (1997) (14) | Taiwan | Retrospective cohort | 1981–1994 | 18,077 | 2,308 | 15,769 | – | – | Capture-recapture | Hevac B, Institut Pasteur, Marnes-la-Coquette, 0.5 mL (145 IU) |
| Lee (1997) (21) | Taiwan | Retrospective cohort | 1984 | – | – | – | – | – | – | – |
| Lee (1998) (22) | South Korea | Retrospective cohort | 1984–1986 | 370,285 | 35,934 | 237,343 | – | 370,285 (100.0) | Hepa S-Ab TestR, Green Cross, Seoul | HepaVaxR, Green Cross, 20 mL |
| Chang (2000) (16) | Taiwan | Retrospective cohort | 1981–1996 | 271 | 76 | 68 | – | 201 (74.16) | Capture-recapture | – |
| Chang (2005) (15) | Taiwan | Retrospective cohort | 1981–2000 | 298 | 35 | 263 | – | – | Capture-recapture | – |
| Viviani (2008) (26) | The Gambia | Prospective cohort | 1990 | 124,577 | 61,065 | 63,512 | – | 401 (0.3) | – | Plasma derived HBV vaccine, France, 10 mL |
| Chang (2009) (18) | Taiwan | Retrospective cohort | 1983–2004 | – | – | – | 6 to 29 | – | – | Plasma-derived HBV vaccine/recombinant yeast-derived vaccine |
| McMahon (2011) (24) | USA | Retrospective cohort | 1985 | 73,798 | 52,000 | 21,798 | Below 20 | – | AFP result >15 ng/mL underwent a liver US or CT | – |
| Chien (2014) (20) | Taiwan | Prospective cohort | 1984–2000 | 4,132,960 | 3,386,988 | 745,972 | – | – | – | Plasma-derived vaccine, 0.5 mL |
| Qu (2014) (25) | China | RCT | 1985 and 1990 | 39,292 | 38,366 | 34,441 | – | 18,673 (47.5) | B-ultrasonic scan; CT; MRI | Plasma-derived HBV vaccine (5 mg) & 10 mg recombinant vaccine (GlaxoSmithKline) |
| Chang (2016) (17) | Taiwan | Retrospective cohort | 1983–2011 | 1,509 | 166 | 1,343 | 6 to 26 | 1,133 (75.1) | Serum AFP level of >400 ng/mL & imaging | Plasma-derived HBV vaccine/recombinant vaccine |
| Liao (2021) (23) | Taiwan | Retrospective cohort | 1979–2017 | 458,171 | 236,948 | 221,223 | 0 to 84 | 181,618 (39.6) | Based on the AJCC system from 2004 and on the Barcelona Clinic Liver Cancer (BCLC) system from 2010 | – |
| Cao (2022) (13) | China | RCT | 1980–1990 | 82,866 | 41,136 | 41,730 | 10 to 14 | 42,826 (51.7) | ELISA | Plasma-derived HBV vaccination, 5-μg, 3 doses |
| Chiang (2022) (19) | Taiwan | Retrospective cohort | 1979–2018 | 43,604 | – | – | 33 (6) | 37,755 (86.6) | International Classification of Diseases for Oncology databases | Recombinant HBV vaccines, 3 doses |
| Wong (2022) (27) | China | Retrospective cohort | 1988–2002 | 372,250 | 53,960 | 318,290 | 30.5 (8.1) | 118,410 (31.8) | ICD-9-CM diagnosis codes | – |
| Yu (2022) (3) | China | Retrospective cohort | 1973–2014 | 66,036 | 51,234 | 14,802 | – | 47,344 (71.69) | ELISA | Plasma-derived hepatitis B vaccine, 3 doses |
Age are presented as range or mean (SD). “–” denotes variables that were not reported in the source publication. AFP, α-fetoprotein; AJCC, American Joint Committee on Cancer; CT, computed tomography; ELISA, enzyme-linked immunosorbent assay; HBV, hepatitis B virus; ICD, International Classification of Diseases; MRI, magnetic resonance imaging; RCT, randomized controlled trial; SD, standard deviation; US, ultrasound.
The mean age at follow-up ranged from childhood to early adulthood, with male representation varying between 31.8% and 86.6%. HCC diagnosis methods included capture-recapture algorithms, α-fetoprotein screening combined with imaging, and registry-based ascertainment using AJCC, BCLC, or ICD coding systems. Both plasma-derived and recombinant yeast-derived vaccines were represented, typically administered as three-dose primary series.
Methodological quality
Quality appraisal using the NIH Quality Assessment Tool for Observational Cohort and Cross-Sectional Studies indicated generally robust study conduct (Table 2). Five studies (35.7%) were rated good quality, and nine (64.3%) were fair quality. Strengths included clearly defined objectives, standardized exposure/outcome measurement, and complete follow-up in registry-based cohorts. Common limitations involved incomplete control of residual confounding, absence of assessor blinding, and limited reporting of potential effect modifiers. No study was rated as poor quality. As for the two included randomized trials, one had some concerns (14) and the other had high risk of bias (26).
Table 2
| ID | Design | Q1 | Q2 | Q3 | Q4 | Q5 | Q6 | Q7 | Q8 | Q9 | Q10 | Q11 | Q12 | Q13 | Q14 | Total score | Quality |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Chang 1997 | RC | 1 | 2 | 0 | 1 | 0 | 1 | 2 | 2 | 2 | 2 | 2 | 1 | 0 | 0 | 16 | Fair |
| Lee 1997 | RC | 1 | 2 | 0 | 1 | 0 | 1 | 2 | 2 | 2 | 2 | 2 | 1 | 0 | 0 | 16 | Fair |
| Lee 1998 | RC | 2 | 2 | 1 | 2 | 0 | 1 | 2 | 2 | 2 | 1 | 2 | 1 | 2 | 0 | 20 | Fair |
| Chang 2000 | RC | 2 | 2 | 2 | 1 | 0 | 1 | 2 | 2 | 2 | 1 | 2 | 2 | 2 | 0 | 21 | Good |
| Chang 2005 | RC | 2 | 2 | 1 | 1 | 0 | 1 | 2 | 2 | 2 | 1 | 2 | 1 | 0 | 0 | 17 | Fair |
| Viviani 2008 | PC | 2 | 2 | 1 | 1 | 0 | 2 | 2 | 2 | 2 | 2 | 2 | 1 | 2 | 0 | 21 | Good |
| Chang 2009 | RC | 2 | 2 | 2 | 1 | 0 | 1 | 2 | 2 | 2 | 2 | 2 | 1 | 2 | 2 | 23 | Good |
| McMahon 2011 | RC | 2 | 2 | 2 | 1 | 0 | 1 | 2 | 1 | 2 | 1 | 2 | 1 | 0 | 0 | 17 | Fair |
| Chien 2014 | PC | 1 | 2 | 2 | 1 | 0 | 2 | 2 | 2 | 2 | 2 | 2 | 1 | 2 | 0 | 21 | Good |
| Chang 2016 | RC | 2 | 2 | 2 | 1 | 0 | 1 | 2 | 1 | 2 | 2 | 2 | 1 | 0 | 2 | 20 | Fair |
| Liao 2021 | RC | 2 | 0 | 0 | 1 | 0 | 1 | 2 | 2 | 2 | 2 | 2 | 1 | 0 | 0 | 15 | Fair |
| Chiang 2022 | RC | 1 | 2 | 0 | 1 | 0 | 1 | 1 | 2 | 2 | 2 | 2 | 1 | 0 | 0 | 15 | Fair |
| Wong 2022 | RC | 2 | 2 | 2 | 2 | 0 | 1 | 2 | 2 | 2 | 2 | 2 | 1 | 2 | 0 | 22 | Good |
| Yu 2022 | RC | 1 | 2 | 1 | 1 | 0 | 1 | 2 | 1 | 2 | 2 | 2 | 1 | 0 | 0 | 16 | Fair |
A modified 0/1/2-point scoring of the 14-item NIH Quality Assessment Tool was applied (maximum 28 points), as previously described in published meta-analyses, to allow finer discrimination between studies. Good quality was defined as an overall score of 21 or above. Fair quality was defined as an overall score of 14-20. Poor quality was defined as an overall score below 14. NIH Quality Assessment Tool items (Q1–Q14): Q1, Was the research question or objective in this paper clearly stated? Q2, Was the study population clearly specified and defined? Q3, Was the participation rate of eligible persons at least 50%? Q4, Were all the subjects selected or recruited from the same or similar populations (including the same time period)? Were inclusion and exclusion criteria for being in the study prespecified and applied uniformly to all participants? Q5, Was a sample size justification, power description, or variance and effect estimates provided? Q6, For the analyses in this paper, were the exposure(s) of interest measured prior to the outcome(s) being measured? Q7, Was the timeframe sufficient so that one could reasonably expect to see an association between exposure and outcome if it existed? Q8, For exposures that can vary in amount or level, did the study examine different levels of the exposure as related to the outcome (e.g., categories of exposure, or exposure measured as continuous variable)? Q9, Were the exposure measures (independent variables) clearly defined, valid, reliable, and implemented consistently across all study participants? Q10, Was the exposure(s) assessed more than once over time? Q11, Were the outcome measures (dependent variables) clearly defined, valid, reliable, and implemented consistently across all study participants? Q12, Were the outcome assessors blinded to the exposure status of participants? Q13, Was loss to follow-up after baseline 20% or less? Q14, Were key potential confounding variables measured and adjusted statistically for their impact on the relationship between exposure(s) and outcome(s)? PC, prospective cohort; RC, retrospective cohort.
Pooled effect of HBV vaccination on HCC incidence
Seven studies provided sufficient raw data to estimate IRRs between vaccinated and unvaccinated populations. The pooled analysis demonstrated a 65% reduction in HCC incidence among vaccinated individuals (IRR =0.35, 95% CI: 0.26–0.47; I2=76.5%), consistent with substantial between-study heterogeneity (Figure 2).
By country/region
When stratified by country/region, a consistent protective effect was observed across all settings. The reduction in HCC incidence was most pronounced in mainland of China (IRR =0.15, 95% CI: 0.03–0.66) and Taiwan (IRR =0.34, 95% CI: 0.24–0.46), while the single study from South Korea showed a nonsignificant association (IRR =0.59, 95% CI: 0.32–1.09).
By gender
Both sexes benefited from vaccination, though the magnitude of risk reduction differed. Among males, the pooled IRR was 0.25 (95% CI: 0.12–0.53; I2=94.9%), compared with 0.43 (95% CI: 0.20–0.94; I2=87.2%) among females, indicating a stronger relative effect in males.
By age group
Stratification by age at diagnosis revealed a clear, consistent protective trend across all age categories. The pooled IRR ranged from 0.28 (95% CI: 0.21–0.38) for ages 6–9 years to 0.42 (95% CI: 0.32–0.56) for ages 20–26 years, with negligible heterogeneity (I2=0%).
By combined age-gender subgroups
The combined age-gender analysis demonstrated sustained vaccine effectiveness across strata. Among males, IRRs ranged from 0.26 (95% CI: 0.19–0.37) in those aged 6–9 years to 0.46 (95% CI: 0.04–0.83) in those aged 20–26 years. Corresponding estimates for females were 0.39 (95% CI: 0.19–0.78), 0.38 (95% CI: 0.24–0.58), and 0.39 (95% CI: 0.23–0.68) across the same age strata, all indicating a substantial and consistent reduction in HCC risk. Heterogeneity within these subgroups was minimal (I2=0%), supporting the robustness of these associations.
Adjusted analyses
Twelve studies reported adjusted effect estimates accounting for major covariates such as sex, age, birth cohort, and socioeconomic or clinical factors. The pooled aRR of HCC among vaccinated individuals was 0.33 (95% CI: 0.24–0.46; I2=83.4%), indicating a 67% lower adjusted risk of HCC compared with non-vaccinated counterparts (Figure 3). The leave-one-out sensitivity analysis showed no change in the reported estimate (Figure S1). The funnel plot showed symmetry (Figure S2) and Egger’s regression test did not indicate statistically significant funnel-plot asymmetry (P=0.6879); however, this test is known to have limited statistical power when fewer than 10–15 studies are available, and its result should therefore be interpreted with caution.
By country
Subgroup analyses demonstrated consistent vaccine-associated reductions in HCC risk across regions. The strongest association was observed in mainland of China (aRR =0.20, 95% CI: 0.15–0.27; I2=0%), followed by Taiwan (aRR =0.36, 95% CI: 0.24–0.54; I2=87.8%) and South Korea (aRR =0.58, 95% CI: 0.31–1.09). The single study from the United States reported an extremely low aRR (0.01, 95% CI: 0.00–0.61), although the wide CI reflects the rarity of HCC events in that cohort and the small sample size.
By vaccine type
Plasma-derived vaccines were associated with a significantly reduced HCC risk (aRR =0.34, 95% CI: 0.25–0.47; I2=0%). The pooled estimate for recombinant vaccines was directionally consistent with a protective effect but did not reach statistical significance (aRR =0.16, 95% CI: 0.01–3.72; I2=76.0%), with the wide CI reflecting the limited number of contributing studies (n=2) and substantial between-study heterogeneity. Direct head-to-head comparison between vaccine types was therefore not possible.
By vaccine dose
Analyses stratified by vaccine completion status indicated that both two-dose and three-dose regimens substantially reduced HCC risk. The pooled aRRs were 0.23 (95% CI: 0.07–0.39; I2=21%) for individuals who received at least two doses and 0.28 (95% CI: 0.18–0.38; I2=46%) for those who completed three doses.
Discussion
In this systematic review and meta-analysis spanning >5 million individuals across multiple regions, HBV vaccination was associated with a large and durable reduction in incident HCC. The pooled unadjusted effect showed ~65% lower HCC incidence among vaccinated versus unvaccinated populations, and multivariable models confirmed an ~70% risk reduction after adjustment for confounders. These magnitudes are consistent with, and extend, seminal population studies that first demonstrated cancer prevention at scale through HBV immunization. In Taiwan, the introduction of universal infant vaccination in 1984 was followed by a pronounced decline in childhood HCC, providing early quasi-experimental evidence that interrupting HBV transmission translates into fewer cancers. Subsequent follow-up into adolescence and young adulthood showed that benefits persisted as vaccinated birth cohorts aged, reinforcing a causal link between vaccination and reduced carcinogenesis (15,19).
Our findings align with broader global syntheses documenting that universal HBV vaccination programs reduce HBsAg seroprevalence in children, curtail chronic HBV infection, and ultimately lower HCC burden. A recent narrative and data-driven review concluded that the scale-up of HepB3 has driven substantial declines in chronic infection and HCC across regions, with the strongest effects in highly endemic settings. Modelling from the Global Burden of Disease liver collaboration similarly suggests that the counterfactual prevalence of HBsAg would have been markedly higher in 2019 in the absence of vaccination, implying averted downstream cirrhosis and HCC over the life course (1,29).
Evidence from non-Asian settings corroborates generalizability. Among Alaska Natives, combined strategies—universal newborn vaccination, immunization of susceptibles, and surveillance—virtually eliminated symptomatic HBV and HCC in children, underscoring the impact of high coverage and programmatic completeness in a distinct genotype and epidemiologic context. Earlier program reports and subsequent updates from this population are concordant with our pooled effect sizes and support real-world effectiveness beyond East Asia (25,30).
Biological plausibility is strong. Chronic HBV drives HCC through persistent inflammation, integration of HBV DNA into the host genome, and resultant clonal expansion; preventing initial infection interrupts this cascade. Long-term cohort studies indicate that even when anti-HBs titers wane, immune memory persists and anamnestic responses are rapid upon antigenic challenge—helping to explain durable protection we and others observe years after primary series completion. Contemporary reviews of vaccine immunology concur that sustained memory underlies long-term effectiveness, mitigating concerns about routine boosting in the general population (6,18,31).
Our subgroup analyses help refine policy-relevant questions. We observed comparable risk reductions with plasma-derived and recombinant vaccines—mirroring programmatic transitions globally without loss of effectiveness. The protective association with ≥2 doses (with further benefit at ≥3 doses) aligns with the immunogenicity profiles of modern schedules, while reaffirming that timely completion—especially the birth dose—remains critical to block perinatal and early horizontal transmission. The World Health Organization and national technical groups continue to emphasize universal birth-dose within 24 hours followed by completion of the primary series to achieve elimination targets; our estimates provide outcome-level justification (HCC prevention) for these coverage priorities (8).
The geography- and sex-stratified patterns we found are also coherent with known epidemiology. Larger relative effects in mainland of China and Taiwan likely reflect higher baseline HBV prevalence and a greater proportion of HCC attributable to HBV, amplifying the observable benefit of vaccination. Sex differences—with somewhat stronger relative protection among males—are plausible given men’s higher baseline HCC risk from combined viral and non-viral cofactors; vaccination may therefore yield a larger RR reduction against the HBV-attributable component. These interpretations are consistent with regional time-trend analyses showing accelerated post-program declines among younger male cohorts and with population reviews documenting the largest absolute HCC reductions where HBV is the dominant driver (1,32).
Our study has important strengths. We synthesized both rate-based (person-time) and adjusted comparative estimates, demonstrating consistency across statistical frameworks and mitigating concerns that our results simply reflect differential follow-up or crude confounding. Inclusion of large registry-linked cohorts and randomized or quasi-experimental designs improved external validity and allowed informative subgroup analyses by country/region, sex, age at diagnosis, vaccine formulation, and dose.
Limitations warrant consideration. Most included studies were observational; although we pooled adjusted estimates, residual confounding (e.g., alcohol, metabolic syndrome, aflatoxin exposure, variable surveillance intensity) may remain. Specifically, the rising prevalence of metabolic dysfunction-associated steatotic liver disease (MASLD), variable alcohol-use patterns across populations, aflatoxin exposure in endemic regions, and host-viral factors such as HBV genotype—where genotypes C and F are associated with substantially higher HCC risk than genotypes B or D, as shown in the Alaska Native cohort—are likely contributors to residual confounding that could not be fully controlled in the source studies (33,34). Pooled adjusted estimates may therefore over- or under-estimate the true vaccine effect depending on the direction of unmeasured confounding. Between-study heterogeneity was high in some strata, likely reflecting differences in underlying HBV epidemiology, secular trends, diagnostic coding, and time since program rollout. The particularly high I2 values observed in the sex-stratified analyses (I2=87.2% for females and 94.9% for males) likely reflect variation in baseline HCC incidence rates between sexes, differing follow-up durations, and population-specific risk profiles. Evidence is still concentrated in East Asia (Taiwan, mainland of China, and South Korea contributing the majority of included studies); more long-term follow-up from Africa, South Asia, and the Americas—where genotype distributions, co-infections (HCV/HDV), and metabolic cofactors differ—is needed to confirm generalizability. Finally, lifetime protection into older adulthood remains incompletely characterized as vaccinated cohorts age; ongoing surveillance will clarify the durability of cancer prevention benefits. These caveats echo those raised by authoritative reviews and global policy assessments and point to clear research priorities (1,8).
The public-health implications are substantial. HBV vaccination is among the most effective cancer-prevention interventions ever implemented. Contemporary guidance from WHO and national advisory bodies prioritizes high coverage of the birth dose and completion of the infant series as essential steps toward the 2030 viral hepatitis elimination targets; our pooled estimates translate these programmatic inputs into downstream cancer outcomes. Strengthening perinatal prevention, ensuring cold-chain reliability, reducing missed opportunities, and integrating vaccination with screening and linkage to care for mothers living with HBV can compound benefits across the life course (1,2,35).
Conclusions
In conclusion, across diverse settings, vaccine platforms, and analytic approaches, HBV vaccination is associated with an approximately three-fold reduction in incident HCC. These findings, taken together with decades of mechanistic and programmatic evidence, support universal HBV immunization—including timely birth dose and series completion—as a cornerstone of global liver-cancer prevention strategies. Continued investment in coverage, equity, and long-term cohort surveillance will be pivotal to sustain and extend these gains.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the PRISMA reporting checklist. Available at https://jgo.amegroups.com/article/view/10.21037/jgo-2026-0202/rc
Peer Review File: Available at https://jgo.amegroups.com/article/view/10.21037/jgo-2026-0202/prf
Funding: None.
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jgo.amegroups.com/article/view/10.21037/jgo-2026-0202/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Flores JE, Thompson AJ, Ryan M, et al. The Global Impact of Hepatitis B Vaccination on Hepatocellular Carcinoma. Vaccines (Basel) 2022;10:793. [Crossref] [PubMed]
- Pattyn J, Hendrickx G, Vorsters A, et al. Hepatitis B Vaccines. J Infect Dis 2021;224:S343-51. [Crossref] [PubMed]
- Yu S, Zi X, Zhu Q, et al. Accelerating Decreases in the Incidences of Hepatocellular Carcinoma at a Younger Age in Shanghai Are Associated With Hepatitis B Virus Vaccination. Front Oncol 2022;12:855945. [Crossref] [PubMed]
- Sun HY, Sheng WH, Tsai MS, et al. Hepatitis B virus coinfection in human immunodeficiency virus-infected patients: a review. World J Gastroenterol 2014;20:14598-614. [Crossref] [PubMed]
- Wong VW, Chan HL. Prevention of hepatocellular carcinoma: a concise review of contemporary issues. Ann Hepatol 2012;11:284-93.
- Mironova M, Ghany MG, Hepatitis B. Vaccine: Four Decades on. Vaccines (Basel) 2024;12:439. [Crossref] [PubMed]
- Chang SH, Shie JJ, Yu NY. Enhancing Executive Functions and Handwriting with a Concentrative Coordination Exercise in Children with ADHD: A Randomized Clinical Trial. Percept Mot Skills 2022;129:1014-35. [Crossref] [PubMed]
- World Health Organization. Hepatitis B vaccines: WHO position paper, July 2017 - Recommendations. Vaccine 2019;37:223-5. [Crossref] [PubMed]
- Page MJ, McKenzie JE, Bossuyt PM, et al. The PRISMA 2020 statement: an updated guideline for reporting systematic reviews. BMJ 2021;372: [Crossref] [PubMed]
- Morgan RL, Whaley P, Thayer KA, et al. Identifying the PECO: A framework for formulating good questions to explore the association of environmental and other exposures with health outcomes. Environ Int 2018;121:1027-31. [Crossref] [PubMed]
- Khan MK, Hasan MT. A poisson regression approach to modeling traffic accident frequency in urban areas. American Journal of Interdisciplinary Studies 2022;3:117-56.
- Zeng X, Zhang Y, Kwong JS, et al. The methodological quality assessment tools for preclinical and clinical studies, systematic review and meta-analysis, and clinical practice guideline: a systematic review. Journal of Evidence-based Medicine 2015;8:2-10.
- Jackson D, Bowden J, Baker R. How does the DerSimonian and Laird procedure for random effects meta-analysis compare with its more efficient but harder to compute counterparts? Journal of Statistical Planning and Inference 2010;140:961-70.
- Cao M, Fan J, Lu L, et al. Long term outcome of prevention of liver cancer by hepatitis B vaccine: Results from an RCT with 37 years. Cancer Lett 2022;536:215652. [Crossref] [PubMed]
- Chang MH, Chen CJ, Lai MS, et al. Universal hepatitis B vaccination in Taiwan and the incidence of hepatocellular carcinoma in children. Taiwan Childhood Hepatoma Study Group. N Engl J Med 1997;336:1855-9.
- Chang MH, Chen TH, Hsu HM, et al. Prevention of hepatocellular carcinoma by universal vaccination against hepatitis B virus: the effect and problems. Clin Cancer Res 2005;11:7953-7. [Crossref] [PubMed]
- Chang MH, Shau WY, Chen CJ, et al. Hepatitis B vaccination and hepatocellular carcinoma rates in boys and girls. JAMA 2000;284:3040-2. [Crossref] [PubMed]
- Chang MH, You SL, Chen CJ, et al. Long-term Effects of Hepatitis B Immunization of Infants in Preventing Liver Cancer. Gastroenterology 2016;151:472-480.
- Chang MH, You SL, Chen CJ, et al. Decreased incidence of hepatocellular carcinoma in hepatitis B vaccinees: a 20-year follow-up study. J Natl Cancer Inst 2009;101:1348-55. [Crossref] [PubMed]
- Chiang CJ, Jhuang JR, Yang YW, et al. Association of Nationwide Hepatitis B Vaccination and Antiviral Therapy Programs With End-Stage Liver Disease Burden in Taiwan. JAMA Netw Open 2022;5:e2222367. [Crossref] [PubMed]
- Chien YC, Jan CF, Chiang CJ, et al. Incomplete hepatitis B immunization, maternal carrier status, and increased risk of liver diseases: a 20-year cohort study of 3.8 million vaccinees. Hepatology 2014;60:125-32. [Crossref] [PubMed]
- Lee CL, Ko YC. Hepatitis B vaccination and hepatocellular carcinoma in Taiwan. Pediatrics 1997;99:351-3. [Crossref] [PubMed]
- Lee MS, Kim DH, Kim H, et al. Hepatitis B vaccination and reduced risk of primary liver cancer among male adults: a cohort study in Korea. Int J Epidemiol 1998;27:316-9. [Crossref] [PubMed]
- Liao SH, Chen CL, Hsu CY, et al. Long-term effectiveness of population-wide multifaceted interventions for hepatocellular carcinoma in Taiwan. J Hepatol 2021;75:132-41. [Crossref] [PubMed]
- McMahon BJ, Bulkow LR, Singleton RJ, et al. Elimination of hepatocellular carcinoma and acute hepatitis B in children 25 years after a hepatitis B newborn and catch-up immunization program. Hepatology 2011;54:801-7. [Crossref] [PubMed]
- Qu C, Chen T, Fan C, et al. Efficacy of neonatal HBV vaccination on liver cancer and other liver diseases over 30-year follow-up of the Qidong hepatitis B intervention study: a cluster randomized controlled trial. PLoS Med 2014;11:e1001774. [Crossref] [PubMed]
- Viviani S, Carrieri P, Bah E, et al. 20 years into the Gambia Hepatitis Intervention Study: assessment of initial hypotheses and prospects for evaluation of protective effectiveness against liver cancer. Cancer Epidemiology Biomarkers & Prevention 2008;17:3216-23.
- Wong GLH, Hui VWK, Yip TCF, et al. Universal HBV vaccination dramatically reduces the prevalence of HBV infection and incidence of hepatocellular carcinoma. Alimentary Pharmacology & Therapeutics 2022;56:869-77.
- GBD 2019 Hepatitis B Collaborators. Global, regional, and national burden of hepatitis B, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet Gastroenterol Hepatol 2022;7:796-829. [Crossref] [PubMed]
- McMahon BJ, Rhoades ER, Heyward WL, et al. A comprehensive programme to reduce the incidence of hepatitis B virus infection and its sequelae in Alaskan natives. Lancet 1987;2:1134-6. [Crossref] [PubMed]
- Simons BC, Spradling PR, Bruden DJ, et al. A Longitudinal Hepatitis B Vaccine Cohort Demonstrates Long-lasting Hepatitis B Virus (HBV) Cellular Immunity Despite Loss of Antibody Against HBV Surface Antigen. J Infect Dis 2016;214:273-80. [Crossref] [PubMed]
- Al-Busafi SA, Alwassief A. Global Perspectives on the Hepatitis B Vaccination: Challenges, Achievements, and the Road to Elimination by 2030. Vaccines (Basel) 2024;12:288. [Crossref] [PubMed]
- Domínguez A, Avellón A, Hernando V, et al. Hepatitis B Virus-Related Cirrhosis and Hepatocellular Carcinoma Hospital Discharge Rates from 2005 to 2021 in Spain: Impact of Universal Vaccination. Vaccines (Basel) 2024;12:1254. [Crossref] [PubMed]
- McMahon BJ, Nolen LD, Snowball M, et al. HBV Genotype: A Significant Risk Factor in Determining Which Patients With Chronic HBV Infection Should Undergo Surveillance for HCC: The Hepatitis B Alaska Study. Hepatology 2021;74:2965-73. [Crossref] [PubMed]
- Noubiap JJ, Ndoula ST. Prevention of mother-to-child transmission of hepatitis B: birth-dose vaccination is not enough. Lancet Glob Health 2022;10:e455-6. [Crossref] [PubMed]

