Association Between Nonselective Beta-Blocker Use and Hepatocellular Carcinoma in Patients With Chronic Hepatitis B Without Cirrhosis and Decompensation

Background: Nonselective beta-blockers (NSBBs) can reduce the incidence or mortality of certain types of cancers, and NSBBs exert a protective effect on hepatocellular carcinoma (HCC) in patients with cirrhosis. However, the potential preventive effect of NSBBs has not yet been investigated in patients with chronic hepatitis B (CHB) who have a high HCC risk regardless of the presence of underlying cirrhosis.

Aim: This study evaluated the association between NSBB use and HCC incidence in patients with CHB without cirrhosis and decompensation.

Methods: From the 2000 Longitudinal Generation Tracking Database, we enrolled patients who were newly diagnosed as having CHB from January 2001 to December 2011 and then followed them up for at least 5 years. To estimate the causal effect of NSBBs on the time-to-event outcomes of HCC, a marginal Cox proportional hazards model was used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs).

Results: After adjustment, no significant benefit of HCC risk reduction was observed between the NSBB users and nonusers (adjusted HR, 0.82; 95% CI, 0.52–1.31). The cumulative defined daily dose (cDDD) analysis revealed no significant dose correlation among the three groups [adjusted HR (95% CI): 1.08, (0.56–2.05), 0.54 (0.17–1.77), and 0.76 (0.40–1.42) in the <90 cDDD, 90 to <180 cDDD, and ≥180 cDDD groups, respectively]. Duration-dependent associations were not observed. Multivariable stratified analysis results demonstrated that HCC risk markedly decreased in the patients aged >55 years (adjusted HR, 0.49; 95% CI, 0.25–0.96; p = 0.04).

Conclusion: NSBB did not significantly prevent HCC in the patients with CHB infection without cirrhosis and decompensation. This study provided one of valuable results that it is not clinically required to use NSBBs as recommended chemoprevention for HCC in high-risk patients who have CHB.

1 Introduction

Hepatocellular carcinoma (HCC) is the sixth most common cancer worldwide, accounting for >80% of primary liver malignancies (Yang et al., 2019). The incidence and mortality of HCC remain high, and HCC is estimated to cause >1 million deaths in 2030 (Villanueva, 2019). Hepatitis B virus (HBV) infection is a crucial risk factor for HCC (Thiele et al., 2014). Although HBV vaccination can reduce the risk of chronic hepatitis B (CHB), numerous unvaccinated individuals and HBV carriers still have a risk of HCC (Villanueva, 2019). In addition, patients who exhibit HBV surface antigen (HBsAg) clearance after receiving nucleotide analogue (NA) therapy still have a risk of HCC (Kim et al., 2015; Papatheodoridis et al., 2015). HCC is mostly diagnosed at a late stage, thus resulting in limited treatment options and poor prognosis. Hence, effective strategies must be developed to prevent HCC in patients with CHB.

Trends in drug repurposing have led to an increase in the number of studies exploring the efficacy of chemopreventive drugs in various cancers. Clinical and human studies have reported that nonselective beta-blockers (NSBBs) reduced the incidence or mortality of the cancers of the upper gastrointestinal tract and other types of solid cancers (Chang et al., 2015; Lin et al., 2015; Pantziarka et al., 2016) The anticancer effects of NSBBs were demonstrated by the finding that a long treatment with propranolol, an NSBB, reduced mortality in patients with unresectable and metastatic HCC in a duration-dependent manner (Chang et al., 2019). In human liver cancer cell lines, propranolol inhibited cancer cell proliferation by inducing apoptosis and S-phase arrest (Wang et al., 2018). This effect of propranolol may be attributable to the inhibition of beta-2 adrenergic receptors, reduction of bacterial translocation, and inhibition of fibrosis or angiogenesis through the downregulation of the vascular endothelial growth factor (Hamdy and El-Demerdash, 2012; Kassahun et al., 2012; Zhang et al., 2012; Moriya and Minamino, 2017). For patients with cirrhosis and oesophageal varices, NSBBs are recommended for the primary and secondary prevention of variceal bleeding (Tripathi, 2012; Garcia-Tsao et al., 2017). Some clinical studies have reported that NSBBs exert pleiotropic effects such as reducing the risk of HCC in patients with liver cirrhosis (Nkontchou et al., 2012; Thiele et al., 2015; Herrera et al., 2016; Wijarnpreecha et al., 2021).

Clinical studies have reported controversial findings regarding the role of NSBBs in HCC prevention, and most studies have recruited only patients with cirrhosis as the study population (Kim et al., 2012; Hagberg et al., 2016; Yeh et al., 2019; Wijarnpreecha et al., 2021). No study has examined the association between NSBBs and the incidence of HCC in patients with CHB who may develop HCC without cirrhosis and have a lower risk of HCC than do patients with existing cirrhosis. In addition, most studies have recruited a Western population, whose aetiology of HCC and responses to NSBBs differ from those of the East Asian population (Zhou et al., 1989; Hu et al., 2007). Therefore, this study investigated the association between NSBB use and the incidence of HCC in patients with CHB in the absence of liver cirrhosis and decompensation.

2 Methods

2.1 Data Source

Taiwan’s National Health Insurance (NHI) program covers 99.99% of Taiwan’s population. The claims data of NHI beneficiaries are collected and added to databases by the Health and Welfare Data Science Center (HWDC) of the Ministry of Health and Welfare. We used data from the 2000 Longitudinal Generation Tracking Database (LGTD 2000), which contains the information of 2 million beneficiaries randomly sampled from the NHI database. This subset contains comprehensive information regarding patients’ demographic variables such as age, sex, outpatient visits, and hospitalisations; disease diagnoses; procedure codes; prescription details; and healthcare item costs (Lin et al., 2018). This study also used the Registry for Catastrophic Illness Patients Database. To deduct certain NHI payments and copayments, patients with malignancies are required to apply for the catastrophic illness certification. All patient applications include complete histopathological or imaging confirmation from physicians and are formally reviewed by experts (Hsieh et al., 2019). Therefore, cancer diagnoses in the Registry for Catastrophic Illness Patients Database are highly accurate and can be used for validation.

Datasets released by the HWDC are anonymised and encrypted to protect patients’ privacy; thus, researchers cannot identify individuals. This study was approved by the Joint Institutional Review Board and Ethics Committee (JIRB) of Taipei Medical University, Taipei, Taiwan (TMU-JIRB No. 201909025). The requirement for patients’ informed consent was waived.

2.2 Study Design and Population

In this population-based retrospective cohort study, we enrolled patients who were newly diagnosed as having CHB [International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes: 070.2, 070.3, and V02.61] recorded at least three times in outpatient clinics or once during hospitalization from January 1, 2001, to December 31, 2011. The index date was the first prescription date of an NSBB in the NSBB use cohort and the matched date in the nonuse cohort (Supplementary Figure S1). The washout period for examining baseline characteristics was 1 year before the index date. To ensure the inclusion of the exposure time window and reduce immortal time bias, we included a follow-up period from the 180th day of the index date until HCC diagnosis, liver transplantation, death, or the end of the study period, whichever occurred first.

Figure 1 depicts the selection process of the study population. Patients aged <20 years were excluded. We also excluded individuals with a history of specific HCC-related diseases such as liver cirrhosis, liver decompensation, hepatitis C virus infection, other viral hepatitis, human immunodeficiency virus infection, alcoholic liver disease, haemochromatosis, biliary cirrhosis, Wilson’s disease, and alpha-1 antitrypsin deficiency. To apply a new user design, patients who were prescribed with NSBBs more than once before the CHB diagnosis were excluded. In addition, we excluded those who irregularly used NSBBs and nonusers who initiated NSBB treatment during the exposure time window. We excluded patients who had HCC before the index date, patients who met an endpoint (HCC, liver transplantation, or death) within the 180-days exposure period, and patients with missing medication records. A 1:2 propensity score matching was employed to ensure the comparability of the two study groups and reduce selection and confounding biases (Sturmer et al., 2014). Each NSBB user was randomly propensity score–matched at a ratio of 1:2 with NSBB nonusers on the basis of several sociodemographic variables related to the likelihood of NSBB treatment, namely age, sex, calendar year of the index date, hypertension, diabetes, mental disorders, cardiac diseases, asthma, and chronic obstructive pulmonary disease.

FIGURE 1

FIGURE 1. Flow chart of study population selection. HWDC, Health and Welfare Data Science Center; HCC, hepatocellular carcinoma; HCV, hepatitis C virus; HIV, human immunodeficiency virus; NSBB, nonselective beta blocker; CHB, chronic hepatitis B; COPD, chronic obstructive pulmonary disease.

2.3 Definition of Exposure

Patients who received NSBBs were identified from the database by using Anatomical Therapeutic Chemical (ATC) codes for nonselective beta-blocking agents (C07AA) and alpha- and beta-blocking agents (C07AG). Regular NSBB users were defined as those who received at least two prescriptions of NSBBs in the outpatient clinic and consumed the medication for >90 days during the 180-days exposure time window since the index date. Irregular users were defined as those who did not meet both the above-mentioned criteria. Patients who never received NSBBs during the enrolment period were included in the nonuser cohort. Alpha- and beta-blocking agents exert nitric oxide–generating and antioxidative effects that may be more beneficial to treating cancer than beta-blocking agents alone (Fonseca, 2010; Hickok and Thomas, 2010; Thyagarajan and Sahu, 2018). On the basis of pharmacological mechanisms, NSBBs were classified into only beta-blocking agents (alprenolol, bupranolol, carteolol, nadolol, oxprenolol, pindolol, propranolol, sotalol, and timolol) and alpha- and beta-blocking agents (carvedilol and labetalol). In terms of chemical properties, carvedilol has a longer half-life and lower lipophilicity than does propranolol. Because of variations in the properties of individual NSBBs, we performed another subgroup analysis based on the use of specific NSBBs (propranolol and carvedilol). Patients who received monotherapy with a certain drug throughout the enrolment period were included in the subgroups.

2.4 Covariates and Confounding Factors

Covariates related to HCC were comorbidities, lifestyle exposure, metabolic factors, and medication use. We identified the following comorbidities at the index date: hypertension, hyperlipidaemia, diabetes, and nonalcoholic liver diseases. In addition, some potential chemopreventive agents with robust evidence, namely antiviral therapy (NAs and interferons), aspirin, statin, and metformin, were considered covariates. Lifestyle factors, namely tobacco use, alcohol use, and obesity, were assessed and defined on the basis of physicians’ diagnoses in the databases. Diagnoses of these coexisting diseases for adjustment were considered if recorded more than two times during outpatient visits and once during hospitalization. Except for antiviral therapy, medication use was identified on the basis of a 28-days prescription within 1 year before the index date. The users of antiviral agents were defined as those who received antiviral agents for >1 week. In accordance with the reimbursement criteria of NAs and interferons in Taiwan, patients must reach a certain CHB severity level to be eligible for the reimbursement; for example, they must have an elevated alanine aminotransferase level of at least two times the two-fold upper limit of the normal level in combination with ≥2,000 copies/ml of HBV DNA. The prescription of these two drugs may reflect the disease severity of patients with CHB. Hence, the prescriptions of NAs and interferon were used as a surrogate variable for the adjustment of disease severity. To identify differences in covariates between the NSBB users and nonusers, the standardised mean difference (SMD) was used to balance diagnostics after propensity score matching (Austin, 2009). An absolute value of SMD greater than 0.10 (small effect size) indicated meaningful imbalance. Imbalanced variables were adjusted in the Cox regression model (Nguyen et al., 2017).

2.5 Outcome Measurement

2.5.1 Primary Outcome

The primary outcome was incident HCC. The diagnosis of HCC obtained from outpatient and inpatient records was identified on the basis of ICD-9-CM code 155 and International Classification of Diseases, Tenth Revision, Clinical Modification (ICD-10-CM) code C22. Patients with a diagnosis of HCC were also identified using the Registry for Catastrophic Illness Patient Database. Censored events in this study included liver transplantation and death. Liver transplantation was confirmed on the basis of the diagnosis of liver transplant status (ICD-9-CM code V42.7 and ICD-10-CM code Z94.4) or liver transplant surgery (procedure codes 505, 75020A, or 75020B). The Cause of Death Database, which contains information regarding the cause and date of death, was used to confirm deaths in the study population.

2.5.2 Secondary Outcome

Secondary outcomes included dose- and duration-dependent associations between NSBB use and HCC incidence, effect of pharmacological classes and individual NSBBs, and subgroups by baseline characteristics. We used the ATC and defined daily dose (DDD) to explore drug use in the population. The DDD of each NSBB in this study was based on the treatment of mild-to-moderate hypertension. In the subgroup analysis, the NSBB use group was categorised by their cumulative DDDs (cDDDs) into <90, 90–180, and ≥180 cDDD groups. The duration of medication use in the NSBB users was determined through the summation of all the intervals of NSBB use during the follow-up period. On the basis of pharmacological mechanism, NSBBs were classified as beta-blocking agents or alpha- and beta-blocking agents in the subgroup analysis. Because of variations in the properties of individual NSBBs, we examined propranolol and carvedilol use in one subgroup analysis.

2.6 Statistical Analysis

Baseline continuous variables are reported as the mean ± standard deviation. Categorical variables are reported as the percentage. The survival curves of cumulative HCC incidence were plotted using the Kaplan–Meier method. Differences between the two curves were examined using the log-rank test. To determine the causal effect of NSBBs on the time-to-event outcomes of HCC incidence, a marginal Cox proportional hazard model was used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs) (Chen et al., 2010; Austin, 2014). In the marginal Cox proportional hazard model, a robust sandwich covariance matrix estimator was used to account for clustering and produce unbiased HRs, which have a precise standard error. Univariable analysis was performed to examine potential risk factors for HCC in this population. Multivariable analysis was performed to determine the effect of NSBBs on HCC. Different variables were adjusted in the regression models: 1) age and sex; 2) age, sex, and comorbidities (namely hypertension, hyperlipidaemia, diabetes, and nonalcoholic liver diseases); 3) age, sex, comorbidities, and medication use (namely antiviral therapy, statin, metformin, and aspirin); 4) age, sex, comorbidities, medication use, and lifestyle (namely alcohol use, tobacco use, and obesity); and 5) imbalanced baseline variables that were tested on the basis of the SMD (Nguyen et al., 2017). To examine the subgroups of patients classified by baseline characteristics, we performed multivariable stratified analysis. To investigate the effect of pharmacological classes and individual NSBBs, the classified users were compared with their paired controls, and HRs were adjusted for unbalanced variables. To test the proportional hazards assumption, we determined the significance of the time-dependent explanatory variable and observed no violation of the proportional hazards assumption (Kleinbaum and Klein, 2012).

In the sensitivity analysis, we divided the study into three parts to evaluate the rigidity of the results. First, to reduce immortal time bias (from CHB diagnosis to NSBB prescription), we assigned the NSBB nonusers index dates by matching them on the basis of year of CHB diagnosis. Second, we excluded patients with prior selective beta-blocker use. Lastly, because the effect of NSBBs on HCC prevention may not persist for >10 years, we followed up with every participant for 10 years from their index dates. All data analyses were performed using SAS (version 9.4, SAS Institute Inc., Cary, NC, United States). Statistical significance was defined as a two-tailed p value of <0.05.

3 Results

3.1 Demographic Characteristics of the Study Population

Figure 1 presents the study population. We identified 1,062 patients with NSBB use and 32,491 patients without NSBB use between January 1, 2001, and December 31, 2011. After the propensity score matching of the 3,171 patients recruited in this study, we determined that the mean age of the patients was 51.35 (±13.36) years, and 1,878 (57.41%) of them were men. In total, 1,057 patients who used NSBBs were included in the NSBB use group, and 2,114 patients who did not use NSBBs were included in the non–NSBB use group. Table 1 lists the demographic characteristics of the participants. Overall, 83 (2.62%) patients received a diagnosis of HCC [25 (2.37%) and 58 (2.74%) in the NSBB user and nonuser cohorts, respectively]. The follow-up durations of the two groups were similar during a median follow-up period of 8.18 years (median, 8.18 years; range, 6.30–10.60 years). Most of the patients were followed until the end of the study period (2,826 cases, 89.12%), with few censored owing to death (262 cases, 8.26%). No patient was censored because of liver transplantation.

TABLE 1

TABLE 1. Demographic characteristics of participants.

3.2 Primary Outcome

Figure 2 presents the cumulative incidence of HCC. No significant difference was observed in the Kaplan–Meier curves between the two groups (p = 0.49). The 5-years cumulative incidence of HCC was 1.36% (95% CI, 0.81–2.29) and 1.37% (95% CI, 0.95–1.99) among the NSBB users and nonusers, respectively. The 10-years cumulative incidence of HCC was 3.55% (95% CI, 2.28–5.51) and 3.61% (95% CI, 2.71–4.79) among the NSBB users and nonusers, respectively.

FIGURE 2

FIGURE 2. Cumulative incidence of hepatocellular carcinoma. HCC, hepatocellular carcinoma; NSBB, nonselective beta blocker.

Potential risk factors for HCC were identified on the basis of the results of univariable analysis (Table 2). Old age (HR, 1.04; 95% CI, 1.02–1.05; p < 0.01), male sex (HR, 2.00; 95% CI, 1.20–3.34; p < 0.01), hypertension (HR, 1.88; 95% CI, 1.18–3.00; p < 0.01), diabetes (HR, 2.52; 95% CI, 1.56–4.06; p < 0.01), and metformin use (HR, 3.00; 95% CI, 1.80–5.00; p < 0.01) were associated with a high risk of HCC. After adjustment for covariates, old age and male sex remained independent risk factors.

TABLE 2

TABLE 2. Univariable and multivariable Cox regression model analysis for HCC risk.

We assessed the association between NSBB use and incident HCC by performing multivariable analysis. No significant reduction in HCC risk was observed in the NSBB users compared with the NSBB nonusers, although the risk of HCC appeared to decrease (model 4, HR, 0.85; 95% CI, 0.52–1.38; p = 0.50; Table 2). We calculated HRs in different models as follows: Multivariable HRs were 0.80 (95% CI, 0.49–1.28; p = 0.35) in model 1, 0.80 (95% CI, 0.50–1.29; p = 0.36) in model 2, 0.84 (95% CI, 0.52–1.36; p = 0.48) in model 3, and 0.79 (95% CI, 0.50–1.27; p = 0.34) in model 5. Additional details regarding statistical analysis are provided in the Supplementary Appendix.

3.3 Secondary Outcome

The dose- and duration-dependent associations were not observed between NSBB use and HCC risk (Table 3). The 90–180 cDDD group had a lower adjusted HR (HR, 0.54; 95% CI, 0.17–1.77) than did the <90 cDDD group (HR, 1.08; 95% CI, 0.56–2.05) and ≥180 cDDD group (HR, 0.76; 95% CI, 0.40–1.42). No change in HCC risk was observed among the three groups. The median duration of NSBB use was 1.46 years (interquartile range, 0.61–4.55 years). In terms of duration, HCC risk in the NSBB users was similar to that in the NSBB nonusers (HR, 1.04; 95% CI, 0.54–1.99) among the short-term users (<1 year). Although the adjusted HR considerably decreased in the patients who used NSBBs for ≥5 years, no significant difference was observed between the users and nonusers (HR, 0.50; 95% CI, 0.19–1.37).

TABLE 3

TABLE 3. Dose and duration of NSBB use and risk of HCC.

3.4 Pharmacological Classes and Effect of Individual NSBBs

Properties such as pharmacological classes or effects of individual NSBBs did not affect HCC prevention (Table 4). The patients prescribed beta-blocking agents alone exhibited similar HCC risk to that of the nonusers (HR, 0.90; 95% CI, 0.49–1.67). Among those prescribed alpha- and beta-blocking agents, HCC risk did not significantly decrease, although their adjusted HR was lower than that of their matched counterparts (HR, 0.53; 95% CI, 0.21–1.36). In another subgroup analysis, we investigated the effect of major drugs, namely propranolol (n = 640, 60.6%) and carvedilol (n = 346, 32.7%). The risk of HCC incidence was similar between the propranolol users and their matched nonusers (HR, 0.88; 95% CI, 0.47–1.67). Similarly, the carvedilol users did not have a significantly lower incidence of HCC than did the nonusers (HR, 0.53; 95% CI, 0.19–1.53).

TABLE 4

TABLE 4. Pharmacological classes of NSBBs and risk of HCC.

3.5 Multivariable Stratified Analysis for Aspirin Therapy

In the multivariable stratified analysis, we stratified patients by several baseline characteristics that affect HCC development (Table 5). Most of the values were consistently <1.0 and did not reach statistical significance. The NSBB users had a 51% lower risk of HCC than did the nonusers among the patients aged >55 years (HR, 0.49; 95% CI, 0.25–0.96; p = 0.04). However, by contrast with the other subgroups, the NSBB users without hypertension had a slightly higher but not significant risk of HCC than that of the nonusers without hypertension (HR, 1.05; 95% CI, 0.50–2.20).

TABLE 5

TABLE 5. Multivariate stratified analyses of the association between NSBB use and risk of HCC.

3.6 Sensitivity Analysis

In our main analysis, although we began the follow-up from the 180th day from the index date to eliminate partial immortal time bias, the time between CHB diagnosis and the first NSBB use was suspected to be the immortal period. The results of sensitivity analysis revealed consistency in the findings even when we reduced immortal time bias by assigning matched dates to the NSBB nonusers on the basis of year of CHB diagnosis (Supplementary Table S1, HR, 0.85; 95% CI, 0.52–1.38). The patients who used selective beta-blockers were excluded because they might have been a potential confounder. The results are in accordance with our earlier observations indicating that HCC risk was not significantly reduced in the NSBB users (HR, 0.67; 95% CI, 0.37–1.20). Lastly, we repeated our analysis with 10 years of follow-up on the basis of the hypothesis that the effects of NSBBs would not last for >10 years. This analysis demonstrated that NSBB use did not reduce the risk of HCC even within 10 years (HR, 0.81; 95% CI, 0.50–1.32).

4 Discussion

This nationwide cohort study determined whether the protective effect of NSBBs can affect the risk of HCC in patients with CHB but without advanced liver diseases. Our cohort consisted of patients with CHB but without progress to cirrhosis and decompensation, with the crude HCC incidence rate being 0.31 per 100 person-years; this finding is consistent with that of a study of patients with CHB but not underlying cirrhosis (Do et al., 2014). The results of this study did not reveal an association between NSBB use and a significantly decreased HCC incidence in the patients with CHB.

Studies have indicated a significant association between NSBB use and a significantly lower risk of HCC in patients with cirrhosis (Nkontchou et al., 2012; Herrera et al., 2016; Wijarnpreecha et al., 2021). However, patients with cirrhosis who cannot use NSBBs may have more advanced liver diseases and thus a higher risk of HCC. Therefore, the higher HCC incidence observed in patients with cirrhosis in the non-NSBB group may be associated with severe advanced liver conditions (Huang et al., 2021). This confounding by different disease severity was lower in our study because we included only patients newly diagnosed as having CHB and excluded those with cirrhosis and liver compensation. Patients with cirrhosis had a higher risk of HCC than did those without cirrhosis (Thiele et al., 2014). The haemodynamic-related benefits of NSBBs, such as the reduction in portal hypertension, strongly affect patients with cirrhosis. Therefore, nonhaemodynamic effects, such as HCC prevention, are more prevalent in those with cirrhosis than in those with only CHB (Kim et al., 2012; Hagberg et al., 2016). The results of this study indicated that haemodynamic-related effects might mainly contribute to the protective effects of NSBBs against hepatocarcinogenesis.

Regarding the effect of the dose and duration on the enhancement of causality, only one study analysed the doses of NSBBs and indicated that the protective effect of NSBBs on HCC prevention was observed only in patients receiving a dose of >90 cDDD and not in low-cDDD groups (Yeh et al., 2019). Our findings revealed that the effects of NSBB use on HCC prevention cannot be extrapolated to patients with CHB even under high cumulative doses and increased durations. The results of the subgroup analysis revealed that the carvedilol monotherapy group did not have a significantly lower incidence of HCC than did the nonusers, although the adjusted HR for the carvedilol group was lower than that for the propranolol group. A study demonstrated a significantly lower risk of HCC in patients with cirrhosis by using carvedilol, propranolol, and nadolol (Wijarnpreecha et al., 2021). This finding differs from that of this study possibly because of differences in disease severity between the two studies.

The aetiology in our study population was HBV infection. The process of HCC manifestation in patients with HBV infection differs from that in patients with HCV infection, who mostly develop cirrhosis before HCC (Villanueva, 2019). Studies have reported a more significant protective effect on HCC incidence in patients with HCV- or alcohol-related cirrhosis than in those with HBV-related cirrhosis (Kim et al., 2012; Nkontchou et al., 2012; Herrera et al., 2016; Yeh et al., 2019). The results of our study are consistent with those of a stratified analysis of patients with HBV-associated cirrhosis (Wijarnpreecha et al., 2021). Tolerant doses to NSBBs in different ethnic groups may have resulted in no remarkable effect on HCC prevention in this study. Because the Asian population requires lower doses of NSBBs to reach the target blood pressure and heart rate than does the Caucasian population, the dose used on each day in this study was lower than that used in studies of the Caucasian population (Zhou et al., 1989; Hu et al., 2007). NSBB use appeared to benefit HCC prevention in the Caucasian population more than in the Asian population (Kim et al., 2012; Nkontchou et al., 2012; Herrera et al., 2016; Yeh et al., 2019; Wijarnpreecha et al., 2021).

The results of the multivariable stratified analysis indicated that HCC risk did not significantly decrease regardless of sex, underlying comorbidities, and baseline medication use. However, the patients aged ≥55 years and those using NSBBs had a substantially (51%) lower risk of HCC than did those who did not use NSBBs. By contrast, the patients aged <55 years exhibited no significant difference in HCC risk. This result is in accordance with that of a study conducted in the United States that reported that NSBB use exerted more significant protective effects for HCC in older patients (Wijarnpreecha et al., 2021) Aging is a crucial risk factor for HCC (Yip et al., 2017; Mittal et al., 2018). Thus, this result might be related to mechanisms underlying the effects of NSBB on HCC prevention, including the inhibition of angiogenesis and the beta-adrenergic signaling process, thus resulting in inflammation regulation. The inhibition of angiogenesis and inflammation exerted stronger effects on older patients whose immune system and angiogenic growth were less active than those in younger patients (Moriya and Minamino, 2017). Therefore, a significant reduction in HCC risk might be observed in older patients.

One study reported that nonselective beta blockers may inhibit hepatocarcinogenesis by inhibiting proliferation and inducing apoptosis and S-phase arrest in human liver cancer cell lines (Wang et al., 2018). In vivo and in vitro studies have demonstrated that beta-adrenergic signaling is associated with inflammation, angiogenesis, apoptosis, cellular replication, DNA damage repair, and cellular immune responses (Ming et al., 1995; Pantziarka et al., 2016; Wang et al., 2018). However, its protective effect was not observed in this study. The differences between this real-world study and experimental studies can be explained by the fact that the doses of NSBBs used in an animal or a cell line in these biological studies are not compatible with those used in humans. The optimal in vitro dose of propranolol for anticarcinogenesis in liver cells was 80 μmol/L (80 nmol/ml), which is considerably higher than the 5.3–300 ng/ml (0.02–1.16 nmol/ml) observed in human plasma (Wong et al., 1979; Wang et al., 2018). Therefore, the antitumor effect of NSBBs on liver cancer remains unclear in humans. Numerous pathogenic mechanisms are involved in HBV-associated HCC (Arzumanyan et al., 2013). Hepatocarcinogenesis is not only related to angiogenesis and inflammation, which are the main pathways of anti-HCC effects exerted by NSBBs. Although these processes are inhibited by NSBBs, other mechanisms may lead to the growth of cancer cells, such as the alteration of host gene expression (Nagaya et al., 1987; Tokino et al., 1991).

To the best of our knowledge, this is the first study to comprehensively investigate the effect of NSBBs on HCC prevention in patients with CHB but without cirrhosis and liver decompensation. The confounding by indication bias for advanced liver diseases was reduced in our study. This strong confounder was not allocated appropriately in studies on patients with cirrhosis who had a higher probability of receiving NSBBs for the prevention of oesophageal varices. In addition to avoiding cirrhosis as a confounding factor in our study, we used propensity score matching to ensure the comparability of the indications of NSBBs between the two cohorts. Another strength of this study is that the follow-up period was longer than that of other studies. Because hepatocarcinogenesis gradually advances, a short-term follow-up might not be sufficient for clinical applications. In addition, the LGTD 2000 used as the data source in this study consisted of updated records regarding prescriptions and diagnoses. Thus, we adjusted covariates for medication use, which had been ignored in some studies. Furthermore, our analyses included potentially chemopreventive drugs such as statins, aspirin, metformin, NA, and interferons. Lastly, we used an exposure window to minimise misclassification and ensure new user design. Several sensitivity analyses were performed to eliminate the epidemiological bias related to immortal time, potential confounding by selective beta-blockers, and unclear duration of the effect of NSBBs on cancer prevention.

This study has several limitations that should be addressed. First, this study was limited by the lack of information regarding laboratory data such as serological markers (e.g., HBV DNA levels, HBeAg, and HBsAg) and the genotype of HBV. However, we adjusted for antiviral therapy, which can serve as a surrogate variable for viral load. In addition, we adjusted for alcohol use, tobacco use, and obesity on the basis of related diagnoses; however, patients who were not hospitalised or did not consult a doctor for these problems were not identified. Nevertheless, patients with these diagnoses were considered to have more severe diseases than did those without related diagnoses. Regarding lifestyle factors, information regarding coffee consumption and aflatoxin exposure was lacking. Other information such as patients’ background and behaviour was not recorded, for example, family history and actual adherence to medication and screening. Lastly, the generalizability of these results is subject to certain limitations because most controls were not included in our studies after propensity score matching. However, the similarity of HCC incidence between the overall population in our studies and that in other studies partly supports generalizability.

In conclusion, this study revealed that NSBB use was not associated with decreased HCC occurrence in a nationwide population of patients with CHB. Our findings suggested that the effect of NSBBs on HCC chemoprevention cannot be extrapolated to patients with CHB, although the benefits were observed in patients with cirrhosis in some studies. The remarkable protective effect was only noted in the subgroup of patients aged >55 years. This study provided one of valuable results that it is not clinically required to use NSBBs as recommended chemoprevention for HCC in high-risk patients who have CHB. To evaluate the benefits of NSBBs for HCC prevention, additional studies should be conducted to investigate the anti-HCC effects of NSBBs on high-risk patients.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The data underlying this study are from National Health Insurance Research Database (NHIRD) and can be accessed by request to the Health and Welfare Data Science Center (HWDC) Health and Welfare, Taiwan (http://dep.mohw.gov.tw/DOS/np-2497-113.html) using the information outlined in the Methods section. All applications are reviewed for approval of data release and applicants must follow the Computer-Processed Personal Data Protection Law (http://www.winklerpartners.com/?p=987) and related regulations of National Health Insurance Administration. Requests to access these datasets should be directed to http://dep.mohw.gov.tw/DOS/np-2497-113.html.

Ethics Statement

The studies involving human participants were reviewed and approved by Datasets released by HWDC are anonymised and encrypted to protect patients’ privacy; thus, researchers cannot identify individuals. This study was approved by the Joint Institutional Review Board/Ethics Committee (JIRB) of Taipei Medical University, Taipei, Taiwan (TMU-JIRB No. 201909025). The requirement of patients’ informed consent was waived. The ethics committee waived the requirement of written informed consent for participation.

Author Contributions

H-YC contributed to conception and design, analysis and interpretation of data, drafting the article and revising it critically for important intellectual content, and final approval of the version to be published. HCL contributed to conception and design, analysis and interpretation of data, and final approval of the version to be published. HLL, YSU, and JJK partly contributed to acquisition of data, analysis and interpretation of data, revising it critically for important intellectual content and final approval of the version to be published. LHW contributed to conception and design, analysis and interpretation of data, drafting the article, revising it critically for important intellectual content and final approval of the version to be published. All authors have disclosed any potential competing financial interests regarding the submitted article. All authors have read the manuscript and approved the final version for submission to Frontiers in Pharmacology and those they accept responsibility for the manuscript’s contents.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

This research was supported by Taipei Medical University Hospital (109TMUH-P-13). This manuscript was edited by Wallace Academic Editing.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2021.805318/full#supplementary-material

References

Arzumanyan, A., Reis, H. M., and Feitelson, M. A. (2013). Pathogenic Mechanisms in HBV- and HCV-Associated Hepatocellular Carcinoma. Nat. Rev. Cancer 13, 123–135. doi:10.1038/nrc3449

PubMed Abstract | CrossRef Full Text | Google Scholar

Austin, P. C. (2009). Balance Diagnostics for Comparing the Distribution of Baseline Covariates between Treatment Groups in Propensity-Score Matched Samples. Stat. Med. 28, 3083–3107. doi:10.1002/sim.3697

PubMed Abstract | CrossRef Full Text | Google Scholar

Austin, P. C. (2014). The Use of Propensity Score Methods with Survival or Time-To-Event Outcomes: Reporting Measures of Effect Similar to Those Used in Randomized Experiments. Stat. Med. 33, 1242–1258. doi:10.1002/sim.5984

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, P. Y., Chung, C. H., Chang, W. C., Lin, C. S., Lin, H. H., Dai, M. S., et al. (2019). The Effect of Propranolol on the Prognosis of Hepatocellular Carcinoma: A Nationwide Population-Based Study. PloS one 14, e0216828. doi:10.1371/journal.pone.0216828

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, P. Y., Huang, W. Y., Lin, C. L., Huang, T. C., Wu, Y. Y., Chen, J. H., et al. (2015). Propranolol Reduces Cancer Risk: A Population-Based Cohort Study. Medicine (Baltimore) 94, e1097. doi:10.1097/md.0000000000001097

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Chen, K., and Ying, Z. (2010). Analysis of Multivariate Failure Time Data Using Marginal Proportional Hazards Model. Stat. Sin 20, 1025–1041.

PubMed Abstract | Google Scholar

Do, A. L., Wong, C. R., Nguyen, L. H., Nguyen, V. G., Trinh, H., and Nguyen, M. H. (2014). Hepatocellular Carcinoma Incidence in Noncirrhotic Patients with Chronic Hepatitis B and Patients with Cirrhosis of All Etiologies. J. Clin. Gastroenterol. 48, 644–649. doi:10.1097/mcg.0000000000000015

CrossRef Full Text | Google Scholar

Fonseca, V. A. (2010). Effects of Beta-Blockers on Glucose and Lipid Metabolism. Curr. Med. Res. Opin. 26, 615–629. doi:10.1185/03007990903533681

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia-Tsao, G., Abraldes, J. G., Berzigotti, A., and Bosch, J. (2017). Portal Hypertensive Bleeding in Cirrhosis: Risk Stratification, Diagnosis, and Management: 2016 Practice Guidance by the American Association for the Study of Liver Diseases. Hepatology 65, 310–335. doi:10.1002/hep.28906

PubMed Abstract | CrossRef Full Text | Google Scholar

Hagberg, K. W., Sahasrabuddhe, V. V., Mcglynn, K. A., and Jick, S. S. (2016). Does Angiotensin-Converting Enzyme Inhibitor and β-Blocker Use Reduce the Risk of Primary Liver Cancer? A Case-Control Study Using the U.K. Clinical Practice Research Datalink. Pharmacotherapy 36, 187–195. doi:10.1002/phar.1704

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamdy, N., and El-Demerdash, E. (2012). New Therapeutic Aspect for Carvedilol: Antifibrotic Effects of Carvedilol in Chronic Carbon Tetrachloride-Induced Liver Damage. Toxicol. Appl. Pharmacol. 261, 292–299. doi:10.1016/j.taap.2012.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Herrera, I., Pascual, S., Zapater, P., Carnicer, F., Bellot, P., and María Palazón, J. (2016). The Use of β-blockers Is Associated with a Lower Risk of Developing Hepatocellular Carcinoma in Patients with Cirrhosis. Eur. J. Gastroenterol. Hepatol. 28, 1194–1197. doi:10.1097/MEG.0000000000000677

CrossRef Full Text | Google Scholar

Hickok, J. R., and Thomas, D. D. (2010). Nitric Oxide and Cancer Therapy: the Emperor Has NO Clothes. Curr. Pharm. Des. 16, 381–391. doi:10.2174/138161210790232149

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsieh, C. Y., Su, C. C., Shao, S. C., Sung, S. F., Lin, S. J., Kao Yang, Y. H., et al. (2019). Taiwan's National Health Insurance Research Database: Past and Future. Clin. Epidemiol. 11, 349–358. doi:10.2147/CLEP.S196293

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, H., Jui, H. Y., Hu, F. C., Chen, Y. H., Lai, L. P., and Lee, C. M. (2007). Predictors of Therapeutic Response to Beta-Blockers in Patients with Heart Failure in Taiwan. J. Formos. Med. Assoc. 106, 641–648. doi:10.1016/s0929-6646(08)60021-2

CrossRef Full Text | Google Scholar

Huang, D. Q., Nguyen, M. H., and H, M. (2021). Letter: Hepatocellular Carcinoma Risk in Patients with Non‐selective Beta Blockers. Aliment. Pharmacol. Ther. 54, 1093–1094. doi:10.1111/apt.16576

PubMed Abstract | CrossRef Full Text | Google Scholar

Kassahun, W. T., Guenl, B., Ungemach, F. R., Jonas, S., and Abraham, G. (2012). Expression and Functional Coupling of Liver β2 - Adrenoceptors in the Human Hepatocellular Carcinoma. Pharmacology 89, 313–320. doi:10.1159/000337381

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, G. A., Lee, H. C., Kim, M. J., Ha, Y., Park, E. J., An, J., et al. (2015). Incidence of Hepatocellular Carcinoma after HBsAg Seroclearance in Chronic Hepatitis B Patients: a Need for Surveillance. J. Hepatol. 62, 1092–1099. doi:10.1016/j.jhep.2014.11.031

CrossRef Full Text | Google Scholar

Kim, T. W., Kim, H. J., Chon, C. U., Won, H. S., Park, J. H., Park, D. I., et al. (2012). Is There Any Vindication for Low Dose Nonselective β-blocker Medication in Patients with Liver Cirrhosis? Clin. Mol. Hepatol. 18, 203–212. doi:10.3350/cmh.2012.18.2.203

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleinbaum, D. G., and Klein, M. (2012). “Evaluating the Proportional Hazards assumption,” in Survival Analysis (New York: Springer), 161–200. doi:10.1007/978-1-4419-6646-9_4

CrossRef Full Text | Google Scholar

Lin, C. S., Lin, W. S., Lin, C. L., and Kao, C. H. (2015). Carvedilol Use Is Associated with Reduced Cancer Risk: A Nationwide Population-Based Cohort Study. Int. J. Cardiol. 184, 9–13. doi:10.1016/j.ijcard.2015.02.015

CrossRef Full Text | Google Scholar

Lin, L. Y., Warren-Gash, C., Smeeth, L., and Chen, P. C. (2018). Data Resource Profile: the National Health Insurance Research Database (NHIRD). Epidemiol. Health 40, e2018062. doi:10.4178/epih.e2018062

PubMed Abstract | CrossRef Full Text | Google Scholar

Ming, M., Wu, J., Lachance, S., Delalandre, A., Carrière, S., and Chan, J. S. (1995). Beta-adrenergic Receptors and Angiotensinogen Gene Expression in Mouse Hepatoma Cells In Vitro. Hypertension 25, 105–109. doi:10.1161/01.hyp.25.1.105

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittal, S., Kramer, J. R., Omino, R., Chayanupatkul, M., Richardson, P. A., El-Serag, H. B., et al. (2018). Role of Age and Race in the Risk of Hepatocellular Carcinoma in Veterans with Hepatitis B Virus Infection. Clin. Gastroenterol. Hepatol. 16, 252–259. doi:10.1016/j.cgh.2017.08.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Moriya, J., and Minamino, T. (2017). Angiogenesis, Cancer, and Vascular Aging. Front. Cardiovasc. Med. 4, 65. doi:10.3389/fcvm.2017.00065

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagaya, T., Nakamura, T., Tokino, T., Tsurimoto, T., Imai, M., Mayumi, T., et al. (1987). The Mode of Hepatitis B Virus DNA Integration in Chromosomes of Human Hepatocellular Carcinoma. Genes Dev. 1, 773–782. doi:10.1101/gad.1.8.773

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. L., Collins, G. S., Spence, J., Daurès, J. P., Devereaux, P. J., Landais, P., et al. (2017). Double-adjustment in Propensity Score Matching Analysis: Choosing a Threshold for Considering Residual Imbalance. BMC Med. Res. Methodol. 17, 78. doi:10.1186/s12874-017-0338-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Nkontchou, G., Aout, M., Mahmoudi, A., Roulot, D., Bourcier, V., Grando-Lemaire, V., et al. (2012). Effect of Long-Term Propranolol Treatment on Hepatocellular Carcinoma Incidence in Patients with HCV-Associated Cirrhosis. Cancer Prev. Res. (Phila) 5, 1007–1014. doi:10.1158/1940-6207.capr-11-0450

PubMed Abstract | CrossRef Full Text | Google Scholar

Pantziarka, P., Bouche, G., Sukhatme, V., Meheus, L., Rooman, I., and Sukhatme, V. P. (2016). Repurposing Drugs in Oncology (ReDO)-Propranolol as an Anti-cancer Agent. Ecancermedicalscience 10, 680. doi:10.3332/ecancer.2016.680

PubMed Abstract | CrossRef Full Text | Google Scholar

Papatheodoridis, G. V., Chan, H. L., Hansen, B. E., Janssen, H. L., and Lampertico, P. (2015). Risk of Hepatocellular Carcinoma in Chronic Hepatitis B: Assessment and Modification with Current Antiviral Therapy. J. Hepatol. 62, 956–967. doi:10.1016/j.jhep.2015.01.002

CrossRef Full Text | Google Scholar

Stürmer, T., Wyss, R., Glynn, R. J., and Brookhart, M. A. (2014). Propensity Scores for Confounder Adjustment when Assessing the Effects of Medical Interventions Using Nonexperimental Study Designs. J. Intern. Med. 275, 570–580. doi:10.1111/joim.12197

CrossRef Full Text | Google Scholar

Thiele, M., Albillos, A., Abazi, R., Wiest, R., Gluud, L. L., and Krag, A. (2015). Non-selective Beta-Blockers May Reduce Risk of Hepatocellular Carcinoma: a Meta-Analysis of Randomized Trials. Liver Int. 35, 2009–2016. doi:10.1111/liv.12782

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiele, M., Gluud, L. L., Fialla, A. D., Dahl, E. K., and Krag, A. (2014). Large Variations in Risk of Hepatocellular Carcinoma and Mortality in Treatment Naïve Hepatitis B Patients: Systematic Review with Meta-Analyses. PLoS One 9, e107177. doi:10.1371/journal.pone.0107177

PubMed Abstract | CrossRef Full Text | Google Scholar

Thyagarajan, A., and Sahu, R. P. (2018). Potential Contributions of Antioxidants to Cancer Therapy: Immunomodulation and Radiosensitization. Integr. Cancer Ther. 17, 210–216. doi:10.1177/1534735416681639

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokino, T., Tamura, H., Hori, N., and Matsubara, K. (1991). Chromosome Deletions Associated with Hepatitis B Virus Integration. Virology 185, 879–882. doi:10.1016/0042-6822(91)90564-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Tripathi, D. (2012). Drugs Used in Therapy of portal Hypertension. Clin. Liver Dis. (Hoboken) 1, 136–138. doi:10.1002/cld.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Villanueva, A. (2019). Hepatocellular Carcinoma. N. Engl. J. Med. 380, 1450–1462. doi:10.1056/NEJMra1713263

CrossRef Full Text | Google Scholar

Wang, F., Liu, H., Wang, F., Xu, R., Wang, P., Tang, F., et al. (2018). Propranolol Suppresses the Proliferation and Induces the Apoptosis of Liver Cancer Cells. Mol. Med. Rep. 17, 5213–5221. doi:10.3892/mmr.2018.8476

PubMed Abstract | CrossRef Full Text | Google Scholar

Wijarnpreecha, K., Li, F., Xiang, Y., Xu, X., Zhu, C., Maroufy, V., et al. (2021). Nonselective Beta-Blockers Are Associated with a Lower Risk of Hepatocellular Carcinoma Among Cirrhotic Patients in the United States. Aliment. Pharmacol. Ther. 54, 481–492. doi:10.1111/apt.16490

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, L., Nation, R. L., Chiou, W. L., and Mehta, P. K. (1979). Plasma Concentrations of Propranolol and 4-hydroxypropranolol during Chronic Oral Propranolol Therapy. Br. J. Clin. Pharmacol. 8, 163–167. doi:10.1111/j.1365-2125.1979.tb05815.x

CrossRef Full Text | Google Scholar

Yang, J. D., Hainaut, P., Gores, G. J., Amadou, A., Plymoth, A., and Roberts, L. R. (2019). A Global View of Hepatocellular Carcinoma: Trends, Risk, Prevention and Management. Nat. Rev. Gastroenterol. Hepatol. 16, 589–604. doi:10.1038/s41575-019-0186-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeh, W. S., Yang, S. C., Liang, C. M., Li, Y. C., Tai, W. C., Lee, C. H., et al. (2019). The Role of Non-selective β-Blockers in Compensated Cirrhotic Patients without Major Complications. Medicina (Kaunas) 56, 14. doi:10.3390/medicina56010014

PubMed Abstract | CrossRef Full Text | Google Scholar

Yip, T. C., Chan, H. L., Wong, V. W., Tse, Y. K., Lam, K. L., and Wong, G. L. (2017). Impact of Age and Gender on Risk of Hepatocellular Carcinoma after Hepatitis B Surface Antigen Seroclearance. J. Hepatol. 67, 902–908. doi:10.1016/j.jhep.2017.06.019

CrossRef Full Text | Google Scholar

Zhang, L., Wang, J. N., Tang, J. M., Kong, X., Yang, J. Y., Zheng, F., et al. (2012). VEGF Is Essential for the Growth and Migration of Human Hepatocellular Carcinoma Cells. Mol. Biol. Rep. 39, 5085–5093. doi:10.1007/s11033-011-1304-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H. H., Koshakji, R. P., Silberstein, D. J., Wilkinson, G. R., and Wood, A. J. (1989). Racial Differences in Drug Response. Altered Sensitivity to and Clearance of Propranolol in Men of Chinese Descent as Compared with American Whites. N. Engl. J. Med. 320, 565–570. doi:10.1056/nejm198903023200905

CrossRef Full Text | Google Scholar