Praziquantel and Upper Gastrointestinal Bleeding in Hepatic Schistosomiasis: A Quasi-Experimental Study.

BACKGROUND: There is a general consensus that widespread use of praziquantel in populations where schistosomiasis is endemic prevents development of hepatic schistosomiasis and its complications. However, a few studies have reported discordant findings linking praziquantel to the occurrence of upper gastrointestinal bleeding (UGIB) in some patients with hepatic schistosomiasis and varices.
OBJECTIVE: We explored if there was any causal association between recent praziquantel use (rPZQ) and upper gastrointestinal bleeding in hepatic schistosomiasis in rural Africa. PATIENTS AND METHODS: A quasi-experimental, retrospective case-controlled study was performed. It involved adult patients with past or acute UGIB, varices, periportal fibrosis, and/or cirrhosis. Cases had acute variceal bleeding while controls did not. The outcome was the frequency of lifetime episodes of UGIB and exposure was rPZQ (received praziquantel in the last 11 months from the date of enrollment). The data analysis included 2 × 2 tables, logistic regression, and propensity-score matching. Odds ratios (ORs), average treatment effects (ATEs), and their 95% confidence intervals (CIs) were used for inference.
RESULTS: Over 6 weeks, we enrolled 19 cases with 92 lifetime episodes of UGIB, and 66 controls with 192 lifetime episodes of UGIB. Cases were more likely to experience UGIB than controls following rPZQ (92% vs. 62%; OR 7.6; 95% CI 3.4-17). Factors predictive of more lifetime episodes of UGIB at multivariable analysis included rPZQ (adjusted OR 13; 95% CI 2.9-53), relative leukocytosis (adjusted OR 26; 95% CI 7.6-89), large varices (adjusted OR 5.0; 95% CI 1.7-15), a family member with hepatosplenic schistosomiasis (adjusted OR 19; 95% CI 7.4-51), advanced periportal fibrosis (adjusted OR 8.0; 95% CI 2.6-22), ascites (adjusted OR 14; 95% CI 4.3-47), and jaundice (adjusted OR 32; 95% CI 7.8-128). While the ATE following rPZQ among the treated was 0.40 (95% CI 0.33-0.48).
CONCLUSIONS: Our findings suggest the presence of a plausible causal association between recent praziquantel use and increased frequency of UGIB in our study population.

Key Points

Introduction

Praziquantel (PZQ) is currently the most effective drug for the treatment of schistosomiasis. By 2018, it was estimated that over 235 million people worldwide received PZQ for schistosomiasis. This was mainly as preventive chemotherapy through mass drug administration (MDA) [1, 2]. In hepatic schistosomiasis (HS) due to Schistosoma mansoni, PZQ was shown to be effective in preventing, decreasing, and sometimes even reversing periportal fibrosis. This was reported to occur within 12 months of treatment. PZQ has also been shown to decrease portal vein diameter and to reduce the prevalence of esophageal varices [3-8]. Predictors of successful treatment of periportal fibrosis following PZQ use include female gender, young age, mild-to-moderate periportal fibrosis, hepatitis B surface antigen-negative status, low S. mansoni re-infection rates, and genetics or family pedigree [9-14].Whereas PZQ is considered very safe, over the last 15 years, a number of drug safety concerns have emerged [2, 7]. Most of these have been highlighted in a number of reports from clinical studies, anthropological studies, and drug information leaflets [5, 15–21]. However, a few reports linking PZQ use to the occurrence of hematemesis or upper gastrointestinal bleeding (UGIB) in hepatic schistosomiasis seem to have been overlooked. One community study reported two to 30 out of 1000 persons experienced hematemesis as a side effect of PZQ [5]. One cohort described a high variceal UGIB rate and death rate among those with advanced schistosoma periportal fibrosis, and previous UGIB at 1 year, despite recent PZQ and endoscopic sclerotherapy [22]. A second cohort study found prior use of PZQ in patients with hepatic schistosomiasis (HS) was associated with greater odds of experiencing acute variceal UGIB [23]. A third cohort looked at the long-term effects of PZQ among patients with HS. They showed that acute variceal bleeding was more frequent among those with a history of UGIB than those without any lifetime history of bleeding [24]. These reports suggest the presence of a drug safety signal (information on a new or known adverse event that is potentially caused by a medicine and that warrants further investigation). Yet, no further evaluation of this unanticipated and bizarre drug safety signal has been undertaken [25].

We set out to investigate such phenomena in rural Africa were schistosomiasis is endemic, and mass drug administration of PZQ is frequent. Specifically, we explored if there was any causal association between PZQ use and occurrence of UGIB among patients with UGIB (past or present), varices, and HS that could not be explained by something else.

Patients and Methods

Ethics Approval

The Makerere University, College of Medicine, Institutional Research and Ethics Committee (#REC, REF2011-244) and Uganda National Council for Science and Technology (#HS-1620) approved the study. Written informed consent was obtained from each study participant. The study was conducted in accordance with the ethical standards as laid down in the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards.

Study Design and Setting

A quasi-experimental, retrospective case-controlled study design was performed. Participants were recruited and assessed at Pakwach Health Centre IV, a primary healthcare facility situated along the banks of the Albert Nile in rural northwestern Uganda. This area is a persistent hot spot for S. mansoni infection despite a reported annual mass drug administration coverage of 75% [26, 27].

Study Participants

Eligible participants included all patients attending Pakwach Health Centre IV, aged 12 years and above, who were able to provide written informed consent/assent, and who reported experiencing one or more episode(s) of upper gastrointestinal bleeding during their lifetime. Upper gastrointestinal bleeding was defined as any history of hematemesis (vomiting blood), melena (passing black or charcoal-colored stool), or hematochezia (fresh blood in stool). Any participant who was pregnant, HIV positive (as per medical records), or had any contraindication to upper gastrointestinal endoscopy was excluded.

Defining Clinical Outcome Assessments, Cases, and Controls

The main outcome was the sum of the frequency of UGIB (past episodes of UGIB and/or acute variceal bleeding) experienced by either cases or controls. Cases included inpatients admitted for acute variceal bleeding and liver disease (periportal fibrosis and/or cirrhosis). Controls included outpatients with prior history of upper gastrointestinal bleeding, ultrasound evidence of liver disease (periportal fibrosis and/or cirrhosis), esophageal varices at endoscopy, and did not have any evidence of acute variceal bleeding during the study period. We chose controls as patients who previously bled because of their high risk of re-bleeding [22, 24]. Cases and controls were recruited over the same period of 6 weeks. Figure 1 provides an overview of the recruitment process and outcome assessments.

Fig. 1

Flow diagram showing the recruitment of cases and controls with the total number of episodes of upper gastrointestinal bleeding experience by participants at each stage of the recruitment process

Study Procedures

All cases and controls were subjected to the same study procedures. These included written informed consent, relevant medical history (which included history of previous PZQ use, period between last PZQ use and recruitment into the study, history of blood transfusion, past admissions, and others), a physical examination, routine clinical blood draws (for complete blood count, hepatitis B and C serologies), and urine for circulating cathodic antigen (CCA), stool for S. mansoni ova, trans-abdominal ultrasound (for diagnosis of periportal fibrosis due to schistosomiasis and cirrhosis), and diagnostic esophagogastroduodenoscopy for detection of esophageal varices. These clinical findings were recorded in a clinical report form, entered into a Microsoft Access database, and later exported to Stata version 15 for statistical analysis. The details of these procedures are described elsewhere [26, 28]. Potential confounders included gender, age, clinical jaundice, grade of periportal fibrosis, splenomegaly at ultrasound, ascites at ultrasound, grade of varices, prior blood transfusion, presence of hepatitis B infection, current infection with S. mansoni (urine circulating cathodic antigen test or presence of S. mansoni ova in a Kato Katz stool smear), enlarged portal vein diameter, and recent treatment with PZQ. The treatment variable was any participant treated with PZQ in the last 11 months (recent PZQ use), and treatment-independent variables or covariates included clinical jaundice, large esophageal varices, a family member diagnosed with hepatosplenic schistosomiasis, advanced periportal fibrosis, and total white cell count ≥3 × 109/L.

Statistical Analysis

The initial analyses included 2 × 2 tables with frequencies, proportions, interval estimates, and 95% confidence intervals (CIs). Pearson’s chi-square test, Fisher’s exact test, and the rank-sum test were used for bivariate analysis. In addition, logistic regression with frequency weights of each covariate provided crude odds ratios (ORs), and their 95% CIs. Multivariable logistic regression with frequency weights was performed by purposeful selection of covariates to control for confounding. Post-estimation of average marginal effects was then implemented for the best-fit model and this is summarized in a table format.

We estimated the treatment effects of PZQ (recent PZQ use) on the frequency of UGIB from our observational data by propensity-score matching in Stata using “teffects” command [29-31].

This was estimated using the derived logistic regression model, in which treatment status was regressed on the measured pretreatment characteristics. Matching was 1 to 1. To assume reasonable estimation of experimental-type causal effects, we checked for adequate balance of matching between cases and controls in Stata. This is presented graphically. Both average treatment effects on the study population and average treatment effects on the treated patients with their CIs were generated. Average treatment effects were presented as coefficients. These represented the absolute difference in average outcomes between controls and those recently treated with PZQ. For ease of interpretation of average treatment effects, we calculated the number needed to harm from the average treatment effects (1 divide by-average treatment effects). All tests were two-sided and a p-level significance value of < 0.05 was considered. These results are presented in the text, tables, and as figures.

Results

Frequency of Upper Gastrointestinal Bleeding (UGIB) and Time Since Last Treatment with PZQ

During the months of July and August 2014, we studied 19 cases who experienced a total of 92 episodes of UGIB and 66 controls who experienced a total of 192 episodes of UGIB. Participants were grouped into five different period categories to determine a temporal relationship. These were defined by the period since a patient last used PZQ. Overall, the relative frequency of UGIB decreased with time since PZQ was last used.

Among the cases, 85 out of the 92 episodes of UGIB (93%; 95% CI 78–98) were experienced by those cases who received PZQ within the last 11 months. Among controls, 119 out of the 192 episodes of UGIB (62%; 95% CI 55–69) were experienced by those who received PZQ in the last 11 months. All cases received PZQ, while 13 (20%) controls with 20 lifetimes episodes of UGIB had never received PZQ during their lifetime. Figure 2 summarizes these findings.

Fig. 2

Column chart describing the relative frequency episodes of upper gastrointestinal bleeding by praziquantel (PZQ) use over time among cases and controls. Relative frequency = (Subgroup count /Total count) × 100. PZQ use was assessed as used PZQ < 2 months ago, used PZQ 2–11 months ago, used PZQ 12–59 months ago, used PZQ > 59 months ago, and never used PZQ

Demographics and Relevant Medical History of Cases and Controls

The median age (interquartile range) of cases and controls was 45 (11) years and 43 (14) years, respectively. The male-female ratio among controls was five to eight and among cases it was eight to 11. The lifetime frequencies of hematemesis, melena, and hematochezia among cases were 99%, 86%, and 83%, respectively. While the lifetime frequencies of hematemesis, melena, and hematochezia among controls was 99%, 32%, and 52%, respectively. Recurrent UGIB (two or more episodes of UGIB) was more frequent among cases 15 (79%) than controls 37 (56%). Among cases (19 participants with a total of 92 episodes of UGIB), the median number of episodes experience by each participant was three lifetime episodes. Among controls (66 participants with a total of 192 episodes of UGIB), the median number of episodes experience by each participant was two lifetime episodes. A comparison of these findings using the Mann-Whitney test or rank-sum test generated an exact p-value = 0.036 (see supplementary tables in Online Supplementary Material, OSM). All cases (100%) and almost all controls (99%) reported a past admission for UGIB. Only one case and three controls used alcohol. Only one case and four controls had a positive urine CCA test. None of the cases or controls had ever been subject to diagnostic endoscopy nor had they received any specific treatment for variceal UGIB (propranolol, endoscopic band ligation or sclerotherapy, shunt surgery). Cases were treated medically for acute variceal bleeding. No control developed acute variceal UGIB during the study period and neither did any case re-bleed during the study period. All participants were started on propranolol after endoscopy for secondary prevention of variceal UGIB according to the current standard of care and no participant experienced recurrent UGIB during the 6-week study period [32].

Contingency Table Analysis (Bi-Variable Analysis)

Cases were 3.9 (95% CI 1.1–14) times more likely to experience acute UGIB than controls if they received PZQ in the last 11 months (see OSM).

In addition, cases were 7.7 (95% CI 3.4–17) times more likely to experience more lifetime episodes of UGIB than controls if they received PZQ in the last 11 months (recent PZQ use). This was also true for several factor variables including female gender (crude OR 3.4; 95% CI 1.8–6.3), a family member diagnosed with hepatosplenic schistosomiasis (crude OR 6.2; 95% CI 3.5–11), a relative leukocytosis of ≥ 3 × 109/L (crude OR 3.7; 95% CI 2.0–6.5), large esophageal varices (crude OR 4.6; 95% CI 2.5–8.4), advanced periportal fibrosis grades E or F (crude OR 4; 95% CI 2.2–7.1), and jaundice (crude OR 25; 95% CI 8.6–75) (Table 1).

Table 1

The frequency of upper gastrointestinal bleeding (UGIB) among cases and controls by different factors and their crude odds ratios

Factor variables	Cases (n = 19)	Controls (n = 66)
Total UGIB episodes (%)	92 (100)	194 (100)
	UGIB episodes among cases n (%)	UGIB episodes among controls n (%)	Crude odd ratios (95% CI)
PZQ ≤ 11 months	85 (92)	119 (61)	7.7 (3.4–17)***
Age ≥ 40 years	79 (86)	152 (78)	1.7 (0.9–3.3)
Female	77 (84)	117 (60)	3.4 (1.8–6.3)***
Family member with HS	73 (79)	19 (10)	6.2 (3.5–11)***
Prior blood transfusion	92 (100)	164 (85)	–
Spleen size > 12 cm	92 (100)	178 (92)	–
Ascites at ultrasound	23 (25)	40 (21)	1.3 (0.71–2.3)
Liver flap	0	17 (9)	–
Clinical jaundice	32 (35)	4 (2)	25 (8.6–75)***
WBC ≥3 × 109/L	37 (40)	30 (16)	3.7 (2.0–6.5)***
HBsAg positive	0	10 (5)	–
Urine CCA positive	2 (2)	5 (3)	0.4 (0.09–1.9)
Stool mansoni ova	0	8 (4)	–
Liver fibrosis patterns E or F	36 (39)	27 (14)	4.0 (2.2–7.1)***
Large varices (F3)	76 (83)	99 (51)	4.6 (2.5–8.4)***

Percentages are rounded to the nearest whole number. The point estimates and their 95% CIs for crude odds ratios were obtained from logistic regression

CI confidence interval, PZQ praziquantel, WBC total white cell count, CCA circulating cathodic antigen, HBsAg hepatitis B surface antigen, HS hepatic schistosomiasis

***p value < 0.001

Multivariable Logistic Regression and Average Marginal Effects

Factors predictive of more lifetime episodes of UGIB after multivariable analysis included “recent PZQ use” (adjusted OR 12; 95% CI 2.9–53), relative leukocytosis (adjusted OR 26; 95% CI 7.7–89), large varices (adjusted OR 5.0; 95% CI 1.7–14.8), a family member with HS (adjusted OR 19; 95% CI 7.2–51), advanced periportal fibrosis (adjusted OR 7.6; 95% CI 2.6–22), ascites (adjusted OR 14.2; 95% CI 4.3–47), and jaundice (adjusted OR 31.5; 95% CI 7.8–128) (Table 2). The C-statistic for this best model was 0.94, n = 286. The effects of the above factors on the probability of experiencing more episodes of UGIB were 22 percentage points increase for recent PZQ use, 29 percentage points increase for relative leukocytosis ≥ 3 × 109/L, 14 percentage points increase for large varices, 26 percentage points increase for participants with a family history of hepatosplenic schistosomiasis, 18 percentage points increase for advanced periportal fibrosis (E, F), 23 percentage points increase for ascites, and 30 percentage points increase for jaundice. All these changes were statistically significant (Table 2).

Table 2

Multivariable regression and average marginal effects of various factors for prediction of upper gastrointestinal bleeding

N = 286	AOR (95% CI)	p value	Average marginal effectsDelta method dy/dx (95% CI)	p value
Family history of HS	19 (7.4–51)	0.001	0.26 (0.20–0.32)	0.001
Recent PZQ use (≤ 11months)	13 (2.9–53)	0.001	0.22 (0.10–0.34)	0.001
Large varices (F3 varices)	5 (1.7–15)	0.004	0.14 (0.05–0.23)	0.003
Advanced periportal fibrosis (patterns E or F)	8 (2.6–22)	0.001	0.18 (0.10–0.26)	0.001
WBC ≥ 3 × 109/L	26 (7.6–89)	0.001	0.29 (0.20–0.38)	0.001
Ascites at ultrasound	14 (4.3–47)	0.001	0.22 (0.15–0.32)	0.001
Clinical jaundice	32 (7.8–128)	0.001	0.30 (0.21–0.40)	0.001

AOR adjusted odds ratio, PZQ praziquantel, F3 varices large varices, HS hepatic schistosomiasis, WBC total white cell count, UGIB upper gastrointestinal bleeding

Estimated Average Treatment Effects and Number Needed to Harm

The adjusted estimated treatment effects of PZQ on the frequency of UGIB and on acute variceal bleeding in our study population are summarized in Table 3. The adjusted estimated treatment effect on the frequency of UGIB in the study population was 0.29 (95% CI 0.23–0.35), while the average treatment effect among the treated was 0.41 (95% CI 0.33–0.48). The number needed to harm was 3.3 (95% CI 2.7–4.3) for the average treatment effects for the study population and the number needed to harm was 2.3 (95% CI 2.1–2.9) for the average treatment effects on the treated cases (treated with PZQ in the last 11 months) (Table 3). The estimated average treatment effects of recent PZQ use on acute variceal bleeding and the number needed to harm following recent PZQ use are shown in Table 3.

Table 3

Treatment effects of recent use of PZQ on occurrence acute variceal bleeding and lifetime frequency of UGIB

Treatment effects	n	Coef	(95% CI) p value	NNH (95% CI)
Estimated average treatment effect on lifetime frequency of UGIB in the whole populationa	286	0.30	(0.23–0.36) p = 0.001	3.3 (2.8–4.4)
Estimated average treatment effect on lifetime frequency of UGIB in the treatedb	286	0.40	(0.33–0.48) p = 0.001	2.5 (2.1–3.0)
Estimated average treatment effect on occurrence of acute variceal bleeding in the study populationa	85	0.17	(0.04–0.29) p = 0.008	5 (3.5 to 25)
Estimated average treatment effect on occurrence of acute variceal bleeding in the treatedb	85	0.27	(0.13–0.40) p = 0.001	3.6 (2.4–6.7)

UGIB upper gastrointestinal bleeding, n = sample size, Coef coefficient (absolute difference in outcomes between controls and cases), NNH number needed to harm, 95% CI 95 % confidence intervals

aSupplementary Fig. 1

bSupplementary Fig. 2

Discussion

Our study looked at participants with a history of UGIB, ultrasound evidence of periportal fibrosis and/or cirrhosis, and esophageal varices. It showed participants who reported recent PZQ use had a greater probability of experiencing more episodes of UGIB than other participants (those who had never received PZQ and those who had not used PZQ recently). This probability was greatest among participants with advanced periportal fibrosis, large varices, ascites, jaundice, relative leukocytosis, and a family history of HS.

Chofle and others reported a history of previous PZQ use was associated with 4.4 times odds of developing acute variceal bleeding and a history of prior UGIB was associated with 2.7 times odds of developing acute variceal bleeding [23]. We found that recent PZQ use was associated with higher odds (OR 3.9; OSM) of experiencing acute variceal bleeding.

Richter and others found that recurrent variceal bleeding was more frequent among those with advanced fibrosis. From their study, most episodes of recurrent variceal bleeding occurred within the first 10 months of follow-up [22]. Advanced fibrosis was a predictor of UGIB in the same study. We found a significant association between the use of PZQ in 11 months and the occurrence of UGIB. We identified advanced periportal fibrosis, large varices, ascites, and jaundice as independent predictors of UGIB. Others have reported similar findings [14, 23, 33–35]. Our results also showed individuals with a close family member with HS were more likely to experience more episodes of UGIB. Though a family history of HS has not been directly associated with UGIB, progression of periportal fibrosis before and after PZQ use has been found to cluster in some families [9, 36–38]. From our study, participants who experienced more UGIB also had a relative leukocytosis. Leukocytosis is a frequent finding in UGIB and has been associated with poor outcomes [39, 40].

We explored the causal association between recent treatment with PZQ and UGIB by propensity-score matching [29]. Reasonable balance was achieved between cases and controls for the above factors. The effect of recent treatment with PZQ on the treated cases indicated a significant absolute risk increase (average treatment effect) on the occurrence of UGIB with a number needed to harm of 2.5. Similar observations were made when we looked at the effects of recent treatment with PZQ on the treated regarding the occurrence of acute variceal bleeding. There was a significant absolute risk increase in occurrence of acute variceal bleeding with a number needed to harm of 3.6. These findings are similar to what Domingues and Coutinho reported from their 5-year cohort [24]. Their cohort described the long-term effects of PZQ chemotherapy on HS and found one out of three persons (number needed to harm of 3) developed UGIB after treatment with PZQ. From the number need to harm of 2.4–6.7, we estimated 149–417 of 1000 persons from our study population developed acute variceal bleeding following PZQ. A study from the Albert Nile basin found two to 30 out of 1000 persons (0.9%) experienced hematemesis (a symptom of acute variceal bleeding) as a side effect of PZQ [5]. If we assume that about 10% of any population in an endemic region is likely to develop severe periportal fibrosis and hematemesis, a reasonable correction from our study would be 15–42 out of 1000 persons (149–417 divided by 10) [41]. This figure is similar to what Kabatereine and others reported [5].

We have no clear explanation for these unusual findings. Nonetheless, we postulate UGIB or hematemesis could be an indirect consequence of the primary action of PZQ on worms. We think PZQ causes the death of adult worms. These “dead adult worm lesions” may lead to intrahepatic thrombophlebitis, thrombosis, then a further increase in portal hypertension, and finally variceal bleeding. This phenomenon is probably limited to patients with a history of variceal UGIB who may not be able to compensate for changes in portal vein pressure resulting from “dead adult worm lesions” [42-45].

Our study has limitations We relied on patient reports to define PZQ use and the lifetime frequency of UGIB. This could have been subject to recall bias. For this reason, our results can be interpreted with caution. That said, patients’ reports are valid and reliable outcome measures in care or in estimating recent medication use [15, 46–49]. In addition, patients’ reports are a recognized and reliable source of drug safety signal in pharmacovigilance. Furthermore, our participants’ reports were consistent with their clinical findings and published reports on HS, varices, and UGIB [46, 50–54]. Unintentionally, we could have excluded patients who died in the past, indicating our conclusions are only valid for those who survived. Although we used propensity-score matching to mitigate selection bias, we concede that we were only able to adjust for known confounding factors. Finally, this was a single-center study meaning our treatment effects could have been overstated.

Conclusion

Our study provides evidence suggesting a plausible causal association between recent PZQ use and increased frequency of UGIB in our study population. Our findings also offers another reason supporting the current manufacturers recommendation on the careful use of PZQ in patients with decompensated liver disease and prior UGIB [19-21]. These findings are important for ongoing drug pharmacovigilance and PZQ mass drug-administration programs. The study findings also underscore the pressing need for proper management of recurrent UGIB in this population through use of drugs, bands, and shunts [32].

Below is the link to the electronic supplementary material.

Supplementary file1 (PDF 1192 KB)

Our study found that one out of four persons with varices and hepatosplenic schistosomiasis and/or cirrhosis who recently used PZQ (within the last 11 months) developed acute variceal bleeding.

We also demonstrated a strong causal association between recent praziquantel use and the occurrence of more episodes of UGIB in our study population even after accounting for other potential contributary factors.

These findings are significant in that they validate what others observed and re-emphasize the need for cautious use of praziquantel among patients with hepatic schistosomiasis and/or liver cirrhosis.