Characterisation of the environmental presence of hepatitis A virus in low-income and middle-income countries: a systematic review and meta-analysis

Abstract

Objective

To characterise the environmental presence of hepatitis A virus (HAV) in low- and middle-income countries (LMICs).

Design

Systematic review and meta-analysis.

Data sources

EBSCOhost, PubMed, Scopus, ScienceDirect, Clinical Key and Web of Science were searched. Grey literature was sourced by searching the following electronic databases: Open Grey, National Health Research Database and Mednar.

Eligibility criteria for including studies

Cross-sectional and ecological studies reporting HAV environmental presence and conducted in LMICs between January 2005 and May 2019, irrespective of language of publication.

Data extraction and data synthesis

Relevant data were extracted from articles meeting the inclusion criteria, and two reviewers independently assessed the studies for risk of bias. High heterogeneity of the extracted data led to the results being reported narratively.

Results

A total of 2092 records were retrieved, of which 33 met the inclusion criteria. 21 studies were conducted in Tunisia, India and South Africa, and the rest were from Philippines, Pakistan, Morocco, Chad, Mozambique, Kenya and Uganda. In Tunisian raw sewage samples, the prevalence of HAV ranged from 12% to 68%, with an estimated average detection rate of 50% (95% CI 25 to 75), whereas HAV detection in treated sewage in Tunisia ranged from 23% to 65%, with an estimated average detection rate of 38% (95% CI 20 to 57). The prevalence of HAV detection in South African treated sewage and surface water samples ranged from 4% to 37% and from 16% to 76%, with an estimated average detection rates of 15% (95% CI 1 to 29) and 51% (95% CI 21 to 80), respectively. Over the review period, the estimated average detection rate of environmental HAV presence appeared to have declined by 10%.

Conclusion

The quality of included studies was fair, but sampling issues and paucity of data limited the strength of the review findings.

PROSPERO registration number

CRD42019119592.

Correspondence to Paul Kuodi; [email protected]

Strengths and limitations of this study

• A comprehensive literature search from multiple databases was conducted.

• There were no language restrictions during the literature search.

• Eligible non-English language articles were translated and screened for inclusion.

• Limited number of studies met the inclusion criteria.

• Meta-analysis for important study outcomes was not possible due to high heterogeneity.

Background

Hepatitis A is caused by the hepatitis A virus (HAV), a small non-enveloped, single-stranded RNA virus belonging to the Picornaviridae family.¹ It is estimated that approximately 1.5 million cases of hepatitis A occur annually, which is likely an underestimate due to the under-reporting of asymptomatic HAV infections.² The transmission of HAV occurs mainly through the faecal-oral route via contact with an infectious person and/or through ingestion of contaminated food or water.³ Foods that are most commonly linked with HAV infection include unpeeled, uncooked produce and raw shellfish.⁴ The incidence of HAV infection is influenced by socioeconomic status (SES), as the risk of transmission is increased in areas with poor water quality, sanitation and hygiene.⁵ Hepatitis A disease in humans can be caused by three genotypes of HAV, genotypes I, II and III, which can be further classified into subgenotypes A and B.^{3 6} Globally, HAV genotype I is the most commonly reported.¹

The environmental presence of HAV is considered to be a driver in the force of HAV infection in a population. Thus, the characterisation of HAV environmental presence is crucial to understanding the epidemiology of the disease. The virus can be detected in environmental sources using qualitative or quantitative (real-time (RT)-PCR) molecular methods as well as nucleic acid sequence-based amplification techniques.⁷ In settings where water quality, sanitation and hygiene are poor, the environmental presence of HAV is assumed to be high, leading to high endemicity of hepatitis A in these areas.⁸ Given the risk factors for HAV infection, global reports indicate that increases in SES lead to a reduction in HAV transmission in the relevant populations and that these changes should be accompanied by the reassessment of public health strategies used to control hepatitis A.^{4 5}

Rationale

The WHO recommends routine immunisation against HAV in countries with intermediate HAV endemicity, that is, in settings presumed to have a mix of HAV environmental exposure levels. The WHO recommends that countries routinely collect and review essential information to estimate national burdens of hepatitis A (including disease transmission rates, clinical outcomes, cost of cases and so on) to inform the need for routine immunisation against the disease. In collecting information on the burden of hepatitis A, information on the environmental presence of HAV should be considered as it is a driver of hepatitis A epidemiology. Such synthesised information is largely missing in the field, and this review was conducted to add to the body of knowledge used to inform hepatitis A control strategies.

Methods

This systematic review was registered with PROSPERO (International Prospective Register of Systematic Reviews), and the results are reported using the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines⁹ (see online supplemental appendix 1).

Objective

The objective was to characterise the environmental presence of HAV in low- and middle-income countries (LMICs).

Primary objective

• To quantify the environmental presence of HAV in LMICs.

Secondary objectives

• To describe the seasonal variation of the environmental presence of HAV in LMICs.

• To describe the genotypes of HAV circulating in LMICs.

Study eligibility criteria

Cross-sectional and ecological studies conducted in LMICs between January 2005 and May 2019, irrespective of language, were eligible for inclusion in this review. The inclusion start date of January 2005 was guided by the latest global review period of HAV endemicity used by the WHO to classify countries’ HAV endemicity status.⁵ Eligible studies reported quantifiable concentrations of HAV in different environments, including raw sewage, treated sewage, sludge, surface water, groundwater, drinking water, ocean water, irrigation water, soil, crops, food and/or shellfish.

Patient and public involvement

Patients or the general public were not involved in the design, planning and conduct of this study.

Search strategy

A comprehensive search strategy that included text words (all fields) and medical subject headings was used to search for eligible studies. The search strategy terms included hepatitis A virus, environmental presence and LMICs. The following electronic databases were searched for relevant published studies: EBSCOhost, PubMed, Scopus, ScienceDirect, Clinical Key and Web of Knowledge (Web of Science). Grey literature was sourced by searching the following electronic databases: Open Grey, National Health Research Database and Mednar.

An example of the search strategy applied to PubMed database is outlined in online supplemental appendix 2. This search strategy was adapted for use in the other databases with the relevant filters used to identify studies published between January 2005 and May 2019. The last literature search from all the databases was completed on 15 May 2019.

Study selection

All database search outputs were imported to Rayyan software for blinded screening.¹⁰ The first author (PK) and the second author (JP) independently screened 100% of titles and abstracts. Titles and abstracts of non-English outputs were translated into English using Google Translate.¹¹ Full-text articles of eligible studies were collected and independently screened by PK and JP. Where disagreements occurred regarding the inclusion, consensus was reached through consultation with the last author (BMK). The PRISMA diagram describing the study selection process is shown in figure 1.

Figure 1. PRISMA flow diagram for selection of studies. The figure details the flow of identification and selection of studies included in the review. Application of the search terms to different databases (PubMed, EBSCOhost, Web of Science, ScienceDirect, Scopus, Mednar, Open Grey and Clinical Key) yielded 2092 records. After the removal of duplicates, there were 1827 unique records eligible for screening. Title and abstract screening excluded 1504 records for various reasons: wrong outcome type (1311), wrong publication type (97), wrong study design (15) wrong population/setting (63), duplicates (8) and published outside the target study period (10). 323 records were assessed for eligibility by full-text screening. The full-text screening excluded 290 records because they were either reporting wrong outcome type or did not report quantifiable HAV environmental prevalence. 33 articles met the inclusion criteria. HAV, hepatitis A virus; n, number of records retrieved; PRISMA, Preferred Reporting Items for Systematic Reviews and Meta-Analyses.

Data extraction

Study characteristics and outcomes of interest were extracted from included studies by PK and JP on a predesigned, electronic data extraction form using KoBo Toolbox¹² (see online supplemental appendix 3). Prior to use of the data extraction form, the reliability of the extraction form was assessed by piloting 10 randomly selected articles that met the inclusion criteria. Any disagreements in data extraction were resolved through consultation and consensus with the last author.

Data management and analysis

Estimates of HAV environmental presence from different studies were pooled as weighted mean proportions by environmental source. Pooled estimates were calculated using the metaprop command, and a random effects model was fitted to the data after Freeman-Tukey double arcsine transformation to stabilise the variances.¹³ Inverse-variance weighting method was used in addition to the DerSimonian and Liard method to estimate the variance between studies during the pooling of proportions. Where heterogeneity between studies was found to be low (I² ≤60), pooled outcome measures were reported with 95% CIs for each respective outcome.¹⁴ Where heterogeneity was found to be high (I² ≤60), narrative reporting was used to describe the prevalence ranges for each respective outcome. All analyses were carried out using STATA V.15 software.¹⁵

Risk of bias assessment

Each included article was independently assessed for risk of bias by PK and JP using the Hoy et al¹⁶ tool (modified by Werfalli et al)¹⁷ (see online supplemental appendix 4). Using the checklist tool, a score of 0–5 points was categorised as ‘high risk’ of bias, 6–8 points as ‘moderate risk’ of bias and 9–10 points as ‘low risk’ of bias. Discrepancies in risk of bias assessments were resolved through consultation and consensus with BMK.

Results

Included studies

The database searches yielded 2092 records, from which 265 duplicates were removed (figure 1). A further 1504 records were excluded for failure to meet the inclusion criteria following the screening of titles and abstracts. Full texts of 323 records were screened, from which 33 studies met the final inclusion criteria.

Characteristics of included studies

The included 33 studies were conducted across 10 different LMICs. All the included studies were observational studies reporting environmental presence of HAV as either a primary or a secondary outcome. Twenty-two of the included studies were conducted in Africa, with Tunisia (nine studies) and South Africa (six studies) contributing to majority of studies from the region. The remaining 11 studies were conducted in Asia, with India (six studies) contributing to majority of the studies from this region.

Table 1 displays the characteristics of included studies and summary of study outcomes. All the included studies used RT-PCR methods for HAV detection, while the methods for environmental sample concentration among included studies differed (table 2). Detection of HAV in raw and treated sewage samples were the most frequently reported outcomes among the included studies.

Table 1

Characteristics of included studies and summary of outcomes from individual studies

Citation	Period of study	Location(s) (setting)	Environmental sample	Positive samples, %	Mean concentration or range (GC/L)	Genotypes
South Africa
Rachida et al25	2013	Cato Ridge, KZN (Kwa Zulu-Natal (urban)	Surface water	72.7 (8/11)	NR	V
Osuolale and Okoh26	2012–2013	Eastern Cape (urban)	Treated sewage	8.6 (6/70)	NR	NR
Adefisoye et al27	2012–2013	Stutterheim (urban)	Raw and treated sewage	6.3 (3/48)	NR	NR
Saïd et al19	2012	Gauteng (urban)	Surface water and treated sewage	48.6 (35/70)	NR	IB
Chigor and Okoh20	2010–2011	Buffalo River (rural)	Surface water	41.3 (31/72)	1.5E+01 to 1.9E+05	NR
Venter et al21	1997–2000	Gauteng (urban)	Surface water	16.2 (50/308)	NR	NR
Mozambique
Nenonen et al28	2004	Maputo Bay (Urban)	Shellfish	33.3 (23/69)	NR	IB
Uganda
O’Brien et al29	2016	Kampala/Nakivubo (urban)	Raw and treated sewage, surface water (swamp)	46.7 (28/60)	5.88E+03 to 1.0E+04 (RS)	NR
					1.41E+03 to 1.11+04 (TS)
					2.73E+03 (SW)
Katukiza et al30	2011	Kampala (urban)	Treated sewage	17.1 (7/41)	NR	NR
Kenya
Kiulia et al31	2007–2008	Nairobi (urban)	Surface water (river), raw sewage	21.95 (9/41)	NR	NR
Tunisia
Jebri et al32	2012	North Tunis (urban)	Raw and treated sewage	0.0 (0/48)		NR
Ouardani et al33	2009–2010	Kairouan (urban)	Raw and treated sewage	53.9 (146/271)	6.7E+01 to 5.6E+07 (RS)	IA, IB
					7.2E+01 to 1.4E+08 (TS)
Ouardani et al34	2009–2010	Parkside (rural)	Raw and treated sewage	39.7 (129/325)	4.29E+03 (RS)	NR
					1.56E+03 (TS)
Gharbi-Khelifi et al35	2009	El-Mahres (urban)	Groundwater (well)	50 (5/10)	NR	IA
Béji-Hamza et al36	2007–2008	Mahdia (urban)	Raw and treated sewage	66.7 (100/150)	7.3E+01 to 4.33E+06 (RS)	IA, IB
					2.56E+03 to 2.32E+06 (TS)
Amri et al37	2007–2008	Lake Tunis (urban)	Shellfish	31.48 (17/54)	NR	NR
Hassine et al38	2005–2006	Monastir (urban)	Raw sewage	11.82 (11/93)	NR	NR
Gharbi-Khelifi et al39	2000–2001	Goz Amir (urban)	Raw and treated sewage	35.41 (17/48)	NR	NR
Elamri et al40	2000–2001	Northern Tunisia (urban)	Shellfish	26.1 (6/23)	NR	NR
Morocco
Amdiouni et al41	2009	Casablanca (urban)	Raw sewage	8.0 (4/50)	NR	NR
Benabbes et al42	2006–2010	Atlantic Coast (urban)	Shellfish	2.6 (2/77)	1.0E+02	NR
Chad
Guerrero-Latorre et al43	2009	Dogdore (urban)	Ground and surface water	0.0 (0/16)	NR	NR
Pakistan
Ahmad et al44	2016–2017	Lahore (urban)	Surface water	8.6 (6/70)	NR	IA
Ahmad et al45	2014	Peshawar (rural)	Raw sewage	12.63 (12/95)	NR	NR
Ahmad et al46	2013	Faisalabad (urban)	Raw sewage and surface water	4.0 (3/75)	NR	NR
Khan et al47	2011	Mardan (urban)	Raw and treated sewage and vegetables	10.3 (3/29)	NR	NR
India
Chobe and Arankalle48	2007	Shimla (urban)	Raw and treated sewage, surface water	80 (4/5)	NR	IIIA
Arankalle et al49	2005	Kottayam-Kerala (urban)	Raw sewage and surface water	17.64 (3/17)	1.0E+02 to 2.6E+03 (RS)	NR
Umesha et al50	2003–2005	Sasthan, Mulki (urban)	Shellfish	0.0 (0/194)	NR	IIIA
Phanuwan et al51	2005	Kampung Melayu (urban)	Surface water	70.58 (12/17)	NR	NR
Chadha et al52	2004	Daund, Pune (rural)	Raw sewage and groundwater	33.3 (1/3)	NR	IIIA
Chitambar et al53	2002	Pune (urban)	Raw sewage	90 (9/10)	NR	IB, IIIA
Philippines
Bai et al54	2012–2013	Metro Manila (urban)	Surface water	58.3 (7/12)	NR	IA, IIIA

GC/L, genome copies per litre; NR, not reported; RS, raw sewage; SW, surface water; TS, treated sewage.

Table 2

Summary of HAV detection and environmental samples concentration techniques used for each included study

Study	Concentration method	Detection method
Jebri et al,32 2012, Elamri et al,40 2006, Umesha et al,50 2007, Kiulia et al,31 2009, Mouna et al,55 2010, Be’ji-Hamza et al,36 2014, Amri et al,37 2009, Amdiouni et al,41 2017, Venter et al,21 2007, Saïd et al,19 2014, Rachida et al,25 2016, Ouardani et al,33 2016, Osuolale et al,26 2016, Katukiza et al,30 2013.	Protein precipitation with PEG	RT-PCR
Guerrero et al,43 2011, Khelifi et al,39 2006, Phanuwan et al,51 2006, O’Brien et al,29 2017, Khelifi et al,35 2011, Bai et al,54 2019, Ahmad et al,45 2016, Ahmad et al,46 2015, Ahmad et al,44 2018, Adefisoye et al,27 2015, Chigor et al,20 2012, Ouardani et al,34 2015.	Elution-based protocol	RT-PCR
Nenonen et al,28 2006.	Direct nucleic acid extraction	RT-PCR
Arankalle et al,49 2006, Chobe et al,48 2009.	Not concentrated	RT-PCR
Khan et al,47 2014, Chitambar et al,53 2007, Osuolale et al,26 2016.	Not reported	RT-PCR

HAV, hepatitis A virus; PEG, polyethylene glycol; RT-PCR, real-time PCR.

HAV genotypes

Twelve (36%) of the included studies reported data on phylogenetic analysis and genotypic classifications of HAV strains detected in environmental samples. HAV genotypes I, III and V were identified as circulating in LMICs. HAV genotypes I and V were reported in Southern Africa, while HAV subgenotypes IIIA, IB and IA were reported in Asia, and HAV subgenotypes IA and IB were reported in Northern Africa (figure 2).

Figure 2. Distribution of reported HAV genotypes. The map shows the geographical distribution of different HAV genotypes reported to be circulating in 10 low-income and middle-income countries by studies included in this review. Countries marked with blue colour scheme had studies reporting HAV environmental presence but did not report the specific genotypes circulating. HAV, hepatitis A virus.

Environmental HAV detection

Due to high heterogeneity among included studies, the environmental presence of HAV could not be pooled and thus is narratively reported by type of environmental source (table 1). As displayed in table 1, the detection of HAV varied greatly among included studies. When taken together, the overall detection of HAV in LMICs was found to be highest in raw and treated sewage samples and lowest in shellfish and groundwater (figure 3). These studies also reported HAV detection to be most frequent in summer (34%) and autumn (32%) seasons, with some detection remaining constant in all seasons (figure 4). Over the review period of 2005–2019, the environmental presence of HAV appears to have declined from an estimated average detection rate of 30% (95% CI 12 to 52) to an estimated detection rate of 20% (95% CI 8 to 35) (figure 5).

Figure 3. Proportion of HAV-positive samples. Boxplot diagram showing the proportions of HAV-positive samples in different environmental samples. A: raw sewage N=784, n=305, median 47, range 0-100. B: treated sewage N=622, n=184, median 28, range 0-64. C: shellfish N=407, n=48, median 26, range, 0-33 D: groundwater N=22, n=5, median, range 0-50. E: surface water N=706, n=147, median 22, range 0-76. N, number of samples, n, number of samples testing positive for HAV. HAV, hepatitis A virus.

Figure 4. Average seasonal HAV detection in low-income and middle-income countries. Mean proportions of HAV-positive samples per season. Summer: n=7, 34% of samples testing positive for HAV; autumn: n=6, 33% of samples testing positive for HAV; winter: n=8, 23% of samples testing positive for HAV; spring: n=4, 10% of samples testing positive for HAV. n, number of articles reporting. HAV, hepatitis A virus.

Figure 5. Forest plot showing reported environmental presence of HAV in LMICs by year of data collection (2000-2017). The forest plot shows reported environmental presence of HAV in LMICs by year of data collection. Estimates for proportions of HAV-positive samples by subgroup of 5-year intervals are shown. The study-specific proportions and 95% CIs are denoted by black boxes and black lines. The combined proportion estimates for the subgroups (5-year intervals) and the overall estimate are represented by black diamonds, where diamond width corresponds to 95% CI bounds. Box and diamond heights are inversely proportional to precision of the proportion estimate. The p value for heterogeneity of the overall proportion estimate is shown. HAV, hepatitis A virus; LMICs, low-income and middle-income countries.

HAV environmental detection in Africa

Given that 45% of included studies were conducted in Tunisia and South Africa, the review data presented a unique opportunity to analyse data for these two countries. Due to heterogeneity, however, these results are reported narratively. The prevalence of HAV detection in Tunisian raw sewage ranged from 12% to 68%, with an estimated average detection rate of 50% (95% CI 25 to 75) (figure 6). The prevalence of HAV detected in treated sewage in Tunisia ranged from 23% to 65%, with an estimated average detection rate of 38% (95% CI 20 to 57) (figure 6). The prevalence of HAV detection in South African treated sewage and surface water ranged from 4% to 37% and from 16% to 76%, with estimated average detection rates of 15% (95% CI 1 to 29) and 51% (95% CI 21 to 80), respectively (figure 7).

Figure 6. Proportion of HAV-positive samples from treated and raw sewage in Tunisia. The forest plot shows the presence of HAV in treated and raw sewage samples reported by studies from Tunisia. The estimate for proportions of HAV-positive samples for raw sewage was 50% (CI 25% to 75%) with heterogeneity of 97.45% (N=494), and the estimate for proportions of HAV-positive samples for treated sewage was 38% (20% to 57%) with heterogeneity of 93.3% (N=393). The study-specific proportions and 95% CIs are denoted by black boxes and black lines. The overall estimate is represented by the diamond, where diamond width corresponds to 95% CI bounds. Box and diamond heights are inversely proportional to precision of the proportion estimate. N, number of samples. HAV, hepatitis A virus.

Figure 7. Proportion of HAV-positive samples from surface water and treated sewage in South Africa. The forest plot shows the presence of HAV in surface water and treated sewage samples reported by studies from South Africa. The estimate for proportions of HAV-positive samples for surface water was 51% (CI 21% to 80%) with heterogeneity of 95.63% (N=412), and the estimate for proportions of HAV-positive samples for treated sewage was 15% (1% to 29%) with heterogeneity of 90.12% (N=169). The study-specific proportions and 95% CIs are denoted by black boxes and black lines. The overall estimate is represented by the diamond, where diamond width corresponds to 95% CI bounds. Box and diamond heights are inversely proportional to precision of the proportion estimate. HAV, hepatitis A virus; N, number of samples.

Risk of bias

Most of the included studies (28 of 33) were graded as having low risk of bias (see online supplementary appendix 5). Four studies were judged as having moderate risk of bias due to sampling from a single site. Table 2 displays the summary of methods used for concentration of environmental samples and HAV detection techniques employed for all included studies.

Discussion

This review synthesised available data on the environmental presence of HAV in LMICs. Data included in this review, however, only represent 10 of 109 (9.17%) LMICs. Given this paucity of data, our results cannot be generalised to all LMICs and indicate a need for more research to quantify the environmental presence of HAV in LMIC settings.

The included studies reported varying detection rates of HAV in a variety of environments, including raw sewage, treated sewage, groundwater, surface water, vegetables and shellfish. The high heterogeneity of HAV detection rates reported among these environmental types, however, was not surprising given the diversity of settings that the data were sourced from. The lack of obvious seasonality in HAV environmental presence reported in the review was also not surprising as HAV is a robust virus that is stable at low pH and at high temperatures of up to 80°C, which ensures the virus’s survival in environmental sources throughout the year.¹⁸

When assessing data from countries where heterogeneity allowed for the pooling of results, we were able to further explore the environmental presence of HAV in South Africa and Tunisia.

The prevalence of HAV detection rate in Tunisia was reportedly higher in raw than treated sewage, which is to be expected. In South Africa, however, HAV detection rates appeared to be higher in surface water than treated sewage. Furthermore, assessment of HAV presence in South African urban and rural surface waters was found to be comparable,^19–21 implying that rural and urban populations in the country may be experiencing similar environmental exposure risks for HAV infection. Interestingly, a study that evaluated HAV seroprevalence in South Africa in a similar period (2005–2019) reported HAV seroprevalence trends that would be expected in a population where urban and rural HAV exposure is equal.²² As of 2019, neither South Africa nor Tunisia had introduced routine HAV vaccination programmes.²³ Instead, in these two countries, HAV vaccination is available to few children in private healthcare and at an extra cost.

The HAV epidemiological transition has come into the public health spotlight globally, as some high-income settings are experiencing more frequent outbreaks of hepatitis A and some LMICs are reporting higher burdens of the disease (ie, an increase in hospitalisations and deaths due to HAV infection).²⁴ The phylogenetic analysis included in our review showed regional overlap of HAV genotypes. These results imply ongoing global transmission of HAV, further supporting Jacobsen and Koopman’s⁵ theory that globalisation (including migration and international food trade) may be significantly linked with hepatitis A outbreaks in high-income settings. Second, our results reported an overall decline in HAV environmental detection from 2005 to 2017. This overall decline in HAV environmental presence may be linked with the increase in the burden of hepatitis A observed in some LMICs as declining HAV transmission rates are often linked with an initial increase in the burden of hepatitis A disease in the corresponding population.⁵

Article summary

Strengths of the review

• We rigorously assessed the validity of included studies and quality of evidence was independently assessed by two reviewers.

• A comprehensive literature search was carried out including publications in different languages which were subsequently translated into English language.

Limitations of the review

• Limited number of studies met the inclusion criteria and hence paucity of the data.

• Due to high heterogeneity among included studies, meta-analysis for the main study outcomes was not possible.

Conclusion

This systematic review aimed to characterise the environmental presence of HAV in LMICs. We successfully addressed the aims of the study, although there was paucity of data to pool results. As mentioned previously, the findings of this review should be interpreted with caution as the results represent only 10 of 109 LMICs. Given the findings of this review, further research to provide quality evidence on the environmental presence of HAV in countries with indications of HAV epidemiological transitions should be prioritised. Such information will better the understanding of HAV transmission routes and inform hepatitis A control strategies, including the introduction of HAV immunisation.

The authors are grateful to the librarians at the School of Health Sciences (University of Cape Town) for their invaluable support during the development of the literature search strategy and the sourcing of grey literature.

Footnote

Contributors The original concept of the study was developed by JP, BK and GDH. PK and JP developed and implemented the systematic review protocol under the supervision of BK. Statistical analysis plan was guided by JP, SS and BK. PK and JP performed the literature search, screening of studies, data extraction and analysis under the guidance of BK, GDH and SS. PK wrote the first draft of the manuscript and all authors gave input to the final draft.

Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

Map disclaimer The depiction of boundaries on this map does not imply the expression of any opinion whatsoever on the part of BMJ (or any member of its group) concerning the legal status of any country, territory, jurisdiction or area or of its authorities. This map is provided without any warranty of any kind, either express or implied.

Competing interests None declared.

Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

Patient consent for publication Not required.

Provenance and peer review Not commissioned; externally peer reviewed.

Data availability statement Data are available upon reasonable request. Data used in our review were extracted from publicly available research articles. Steps for retrieval and analysis of the data are outlined in the systematic review methods.

References

¹ Costa-Mattioli M, Napoli AD, Ferré V, et al. Genetic variability of hepatitis A virus. J Gen Virol 2003; 84: 3191–201. doi:10.1099/vir.0.19532-0 http://www.ncbi.nlm.nih.gov/pubmed/14645901

² Franco E, Meleleo C, Serino L, et al. Hepatitis A: epidemiology and prevention in developing countries. World J Hepatol 2012; 4: 68. doi:10.4254/wjh.v4.i3.68 http://www.ncbi.nlm.nih.gov/pubmed/22489258

³ Vaughan G, Goncalves Rossi LM, Forbi JC, et al. Hepatitis A virus: host interactions, molecular epidemiology and evolution. Infect Genet Evol 2014; 21: 227–43. doi:10.1016/j.meegid.2013.10.023 http://www.ncbi.nlm.nih.gov/pubmed/24200587

⁴ Jacobsen KH. Globalization and the changing epidemiology of hepatitis A virus. Cold Spring Harb Perspect Med 2018; 8: a031716. doi:10.1101/cshperspect.a031716 http://www.ncbi.nlm.nih.gov/pubmed/29500305

⁵ Jacobsen KH, Koopman JS. The effects of socioeconomic development on worldwide hepatitis A virus seroprevalence patterns. Int J Epidemiol 2005; 34: 600–9. doi:10.1093/ije/dyi062 http://www.ncbi.nlm.nih.gov/pubmed/15831565

⁶ Endo K, Inoue J, Takahashi M, et al. Analysis of the full-length genome of a subgenotype IIIB hepatitis A virus isolate: primers for broadly reactive PCR and genotypic analysis. J Med Virol 2007; 79: 8–17. doi:10.1002/jmv.20757 http://www.ncbi.nlm.nih.gov/pubmed/17133545

⁷ Nainan OV, Xia G, Vaughan G, et al. Diagnosis of hepatitis A virus infection: a molecular approach. Clin Microbiol Rev 2006; 19: 63–79. doi:10.1128/CMR.19.1.63-79.2006 http://www.ncbi.nlm.nih.gov/pubmed/16418523

⁸ van der Poel W, Rzezutka A. Hepatitis A. In: Rose JB, Jiménez-Cisneros B, eds. Global Water Pathogen Project. http://www.waterpathogens.org/ (J.S Meschke, and R. Girones (eds) Part 3 Viruses) http://www.waterpathogens.org/book/hepatitisA. East Lansing: UNESCO, 2017.

⁹ Moher D, Shamseer L, Clarke M, et al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement. Syst Rev 2015; 4: 1. doi:10.1186/2046-4053-4-1 http://www.ncbi.nlm.nih.gov/pubmed/25554246

¹⁰ Ouzzani M, Hammady H, Fedorowicz Z, et al. Rayyan-a web and mobile APP for systematic reviews. Syst Rev 2016; 5: 210. doi:10.1186/s13643-016-0384-4 http://www.ncbi.nlm.nih.gov/pubmed/27919275

¹¹ Balk EM, Chung M, Chen ML, et al. Assessing the accuracy of Google translate to allow data extraction from trials published in non-English languages, 2013.

¹² Kobo toolbox: data collection for use in challenging environments, 2013. Available: https://www.kobotoolbox.org/

¹³ Freeman MF, Tukey JW. Transformations related to the angular and the square root. Ann Math Statist 1950; 21: 607–11. doi:10.1214/aoms/1177729756

¹⁴ Higgins JPT, Green S. Cochrane handbook for systematic reviews of interventions version 5.1.0 [updated March 2011]. The Cochrane collaboration, 2011, 2011. http://handbook.cochrane.org

¹⁵ StataCorp. Stata statistical software: release 15. College Station, TX: StataCorp LLC, 2017.

¹⁶ Hoy D, Brooks P, Woolf A, et al. Assessing risk of bias in prevalence studies: modification of an existing tool and evidence of interrater agreement. J Clin Epidemiol 2012; 65: 934–9. doi:10.1016/j.jclinepi.2011.11.014 http://www.ncbi.nlm.nih.gov/pubmed/22742910

¹⁷ Werfalli M, Musekiwa A, Engel ME, et al. The prevalence of type 2 diabetes mellitus among older people in Africa: a systematic review study protocol. BMJ Open 2014; 4: e004747.

¹⁸ Wang X, Ren J, Gao Q, et al. Hepatitis A virus and the origins of picornaviruses. Nature 2015; 517: 85–8. doi:10.1038/nature13806 http://www.ncbi.nlm.nih.gov/pubmed/25327248

¹⁹ Saïd R, Wolfaardt M, Taylor MB. Molecular characterisation of hepatitis A virus strains from water sources in South Africa. Water Sci Technol 2014; 69: 923–33. doi:10.2166/wst.2013.799 http://www.ncbi.nlm.nih.gov/pubmed/24622538

²⁰ Chigor VN, Okoh AI. Quantitative RT-PCR detection of hepatitis A virus, rotaviruses and enteroviruses in the buffalo river and source water dams in the Eastern Cape Province of South Africa. Int J Environ Res Public Health 2012; 9: 4017–32. doi:10.3390/ijerph9114017 http://www.ncbi.nlm.nih.gov/pubmed/23202829

²¹ Venter JME, van Heerden J, Vivier JC, et al. Hepatitis A virus in surface water in South Africa: what are the risks? J Water Health 2007; 5: 229–40. doi:10.2166/wh.2007.006b http://www.ncbi.nlm.nih.gov/pubmed/17674572

²² Patterson J, Abdullahi L, Hussey GD, et al. A systematic review of the epidemiology of hepatitis A in Africa. BMC Infect Dis 2019; 19: 651. doi:10.1186/s12879-019-4235-5 http://www.ncbi.nlm.nih.gov/pubmed/31331281

²³ World Health Organization. WHO vaccine-preventable diseases: monitoring system. 2019 global summary. 2019. Available: https://apps.who.int/immunization_monitoring/globalsummary [Accessed 20 May 2020 ].

²⁴ Jacobsen KH, Koopman JS. Declining hepatitis A seroprevalence: a global review and analysis. Epidemiol Infect 2004; 132: 1005–22. doi:10.1017/S0950268804002857 http://www.ncbi.nlm.nih.gov/pubmed/15635957

²⁵ Rachida S, Matsapola PN, Wolfaardt M, et al. Genetic characterization of a novel hepatitis A virus strain in irrigation water in South Africa. J Med Virol 2016; 88: 734–7. doi:10.1002/jmv.24370 http://www.ncbi.nlm.nih.gov/pubmed/26331799

²⁶ Osuolale O, Okoh A. Incidence of human adenoviruses and hepatitis A virus in the final effluent of selected wastewater treatment plants in Eastern Cape Province, South Africa. Virol J 2015; 12: 98. doi:10.1186/s12985-015-0327-z http://www.ncbi.nlm.nih.gov/pubmed/26104284

²⁷ Adefisoye MA, Nwodo UU, Green E, et al. Quantitative PCR detection and characterisation of human adenovirus, rotavirus and hepatitis A virus in discharged effluents of two wastewater treatment facilities in the eastern Cape, South Africa. Food Environ Virol 2016; 8: 262–74. doi:10.1007/s12560-016-9246-4 http://www.ncbi.nlm.nih.gov/pubmed/27236707

²⁸ Nenonen NP, Hernroth B, Chauque AA, et al. Detection of hepatitis A virus genotype IB variants in clams from Maputo Bay, Mozambique. J Med Virol 2006; 78: 896–905. doi:10.1002/jmv.20638 http://www.ncbi.nlm.nih.gov/pubmed/16721847

²⁹ O’Brien E, Nakyazze J, Wu H, et al. Viral diversity and abundance in polluted waters in Kampala, Uganda. Water Res 2017; 127: 41–9. doi:10.1016/j.watres.2017.09.063 29031798

³⁰ Katukiza AY, Temanu H, Chung JW, et al. Genomic copy concentrations of selected waterborne viruses in a slum environment in Kampala, Uganda. J Water Health 2013; 11: 358–70. doi:10.2166/wh.2013.184 http://www.ncbi.nlm.nih.gov/pubmed/23708582

³¹ Kiulia NM, Netshikweta R, Page NA, et al. The detection of enteric viruses in selected urban and rural river water and sewage in Kenya, with special reference to rotaviruses. J Appl Microbiol 2010; 109: 818–28. doi:10.1111/j.1365-2672.2010.04710.x http://www.ncbi.nlm.nih.gov/pubmed/20337763

³² Jebri S, Jofre J, Barkallah I, et al. Presence and fate of coliphages and enteric viruses in three wastewater treatment plants effluents and activated sludge from Tunisia. Environ Sci Pollut Res Int 2012; 19: 2195–201. doi:10.1007/s11356-011-0722-y http://www.ncbi.nlm.nih.gov/pubmed/22261843

³³ Ouardani I, Turki S, Aouni M, et al. Detection and molecular characterization of hepatitis A virus from Tunisian wastewater treatment plants with different secondary treatments. Appl Environ Microbiol 2016; 82: 3834–45. doi:10.1128/AEM.00619-16 http://www.ncbi.nlm.nih.gov/pubmed/27107113

³⁴ Ouardani I, Manso CF, Aouni M, et al. Efficiency of hepatitis A virus removal in six sewage treatment plants from central Tunisia. Appl Microbiol Biotechnol 2015; 99: 10759–69. doi:10.1007/s00253-015-6902-9 http://www.ncbi.nlm.nih.gov/pubmed/26286509

³⁵ Gharbi-Khelifi H, Ben Salem Abid N, Sdiri K, et al. Characterization of outbreak hepatitis A isolates in five Tunisian childcare centers. Braz J Microbiol 2011; 42: 1204–12. doi:10.1590/S1517-83822011000300046 http://www.ncbi.nlm.nih.gov/pubmed/24031743

³⁶ Béji-Hamza A, Khélifi-Gharbi H, Hassine-Zaafrane M, et al. Qualitative and quantitative assessment of hepatitis A virus in wastewaters in Tunisia. Food Environ Virol 2014; 6: 246–52. doi:10.1007/s12560-014-9163-3 http://www.ncbi.nlm.nih.gov/pubmed/25129101

³⁷ Amri I, Hmaïed F, Loisy F, et al. [Hepatitis A virus detection in shellfish from Tunisia by reverse transcription-nested PCR-investigation of a correlation between viral and bacterial contamination]. Pathol Biol 2011; 59: 217–21. doi:10.1016/j.patbio.2009.10.009 19944542

³⁸ Hassine M, Sdiri K, Riabi S, et al. [Detection of enteric viruses in wastewater of Monastir region by RT-PCR method]. Tunis Med 2010; 88: 70–5. http://www.ncbi.nlm.nih.gov/pubmed/20415161

³⁹ Gharbi-Khelifi H, Ferre V, Sdiri K, et al. Hepatitis A in Tunisia: phylogenetic analysis of hepatitis A virus from 2001 to 2004. J Virol Methods 2006; 138: 109–16. doi:10.1016/j.jviromet.2006.08.001 http://www.ncbi.nlm.nih.gov/pubmed/17005259

⁴⁰ Elamri DE, Aouni M, Parnaudeau S, et al. Detection of human enteric viruses in shellfish collected in Tunisia. Lett Appl Microbiol 2006; 43: 399–404. doi:10.1111/j.1472-765X.2006.01978.x http://www.ncbi.nlm.nih.gov/pubmed/16965370

⁴¹ Amdiouni HA, Maunula L, Al-Shuwaikh AM, et al. Comparison of two virus concentration methods for enteric viruses detection in Moroccan wastewater and treated effluent. Iraqi J Med Sci 2017; 15: 27–38.

⁴² Benabbes L, Ollivier J, Schaeffer J, et al. Norovirus and other human enteric viruses in Moroccan shellfish. Food Environ Virol 2013; 5: 35–40. doi:10.1007/s12560-012-9095-8 http://www.ncbi.nlm.nih.gov/pubmed/23412717

⁴³ Guerrero-Latorre L, Carratala A, Rodriguez-Manzano J, et al. Occurrence of water-borne enteric viruses in two settlements based in Eastern Chad: analysis of hepatitis E virus, hepatitis A virus and human adenovirus in water sources. J Water Health 2011; 9: 515–24. doi:10.2166/wh.2011.126 http://www.ncbi.nlm.nih.gov/pubmed/21976198

⁴⁴ Ahmad T, Adnan F, Nadeem M, et al. Assessment of the risk for human health of enterovirus and hepatitis A virus in clinical and water sources from three metropolitan cities of Pakistan. Ann Agric Environ Med 2018; 25: 708–13. doi:10.26444/aaem/99590 http://www.ncbi.nlm.nih.gov/pubmed/30586962

⁴⁵ Ahmad T, Arshad N, Adnan F, et al. Prevalence of rotavirus, adenovirus, hepatitis A virus and enterovirus in water samples collected from different region of Peshawar, Pakistan. Ann Agric Environ Med 2016; 23: 576–80. doi:10.5604/12321966.1226849 http://www.ncbi.nlm.nih.gov/pubmed/28030926

⁴⁶ Ahmad T, Anjum S, Sadaf Zaidi N-us-S, et al. Frequency of hepatitis E and Hepatitis A virus in water sample collected from Faisalabad, Pakistan. Ann Agric Environ Med 2015; 22: 661–4. doi:10.5604/12321966.1185903 26706972

⁴⁷ Khan W, Irshad M, Rehman G, et al. HAV in fresh vegetables: a hidden health risk in district Mardan, Pakistan. Springerplus 2014; 3: 675. doi:10.1186/2193-1801-3-675 http://www.ncbi.nlm.nih.gov/pubmed/26034681

⁴⁸ Chobe L, Arankalle V. Investigation of a hepatitis A outbreak from Shimla, Himachal Pradesh. Indian J Medical Research 2009; 130: 179.

⁴⁹ Arankalle VA, Sarada Devi KL, Lole KS, et al. Molecular characterization of hepatitis A virus from a large outbreak from Kerala, India. Indian J Med Res 2006; 123: 760. http://www.ncbi.nlm.nih.gov/pubmed/16885597

⁵⁰ Umesha KR, Bhavani NC, Venugopal MN, et al. Prevalence of human pathogenic enteric viruses in bivalve molluscan shellfish and cultured shrimp in South West coast of India. Int J Food Microbiol 2008; 122: 279–86. doi:10.1016/j.ijfoodmicro.2007.12.024 http://www.ncbi.nlm.nih.gov/pubmed/18279989

⁵¹ Phanuwan C, Takizawa S, Oguma K, et al. Monitoring of human enteric viruses and coliform bacteria in waters after urban flood in Jakarta, Indonesia. Water Sci Technol 2006; 54: 203–10. doi:10.2166/wst.2006.470 http://www.ncbi.nlm.nih.gov/pubmed/17037154

⁵² Chadha MS, Lole KS, Bora MH, et al. Outbreaks of hepatitis A among children in Western India. Trans R Soc Trop Med Hyg 2009; 103: 911–6. doi:10.1016/j.trstmh.2008.11.014 http://www.ncbi.nlm.nih.gov/pubmed/19155033

⁵³ Chitambar S, Joshi M, Lole K, et al. Cocirculation of and coinfections with hepatitis A virus subgenotypes IIIA and IB in patients from Pune, Western India. Hepatol Res 2007; 37: 85–93. doi:10.1111/j.1872-034X.2007.00025.x http://www.ncbi.nlm.nih.gov/pubmed/17300702

⁵⁴ Bai H, Shiota T, Yoshizaki S, et al. Detection of subgenotype IA and IIIA hepatitis A viruses in rivers flowing through Metro Manila, the Philippines. Jpn J Infect Dis 2019; 72: 53–5. doi:10.7883/yoken.JJID.2018.148 http://www.ncbi.nlm.nih.gov/pubmed/30175732

⁵⁵ Hassine M, Sdiri K, Riabi S, et al. Détection des virus entériques dans les eaux usées de la région de monastir par RT-PCR [Detection of enteric viruses in wastewater of Monastir region by RT-PCR method]. Tunis Med 2010; 88: 70–5.

⁵⁶ Shamliyan T, Kane RL, Dickinson S. A systematic review of tools used to assess the quality of observational studies that examine incidence or prevalence and risk factors for diseases. J Clin Epidemiol 2010; 63: 1061–70. doi:10.1016/j.jclinepi.2010.04.014 http://www.ncbi.nlm.nih.gov/pubmed/20728045