Molecular characterization of hepatitis C virus in liver disease patients in Botswana: a retrospective cross-sectional study
BMC Infectious Diseases volume 19, Article number: 875 (2019)
Hepatitis C virus (HCV) infection is a major cause of chronic liver disease globally. Direct acting antivirals (DAAs) have proven effective in curing HCV. However, the current standard of care (SOC) in Botswana remains PEGylated interferon-α (IFN-α) with ribavirin. Several mutations have been reported to confer resistance to interferon-based treatments. Therefore, there is a need to determine HCV genotypes in Botswana, as these data will guide new treatment guidelines and understanding of HCV epidemiology in Botswana.
This was a retrospective cross-sectional pilot study utilizing plasma obtained from 55 participants from Princess Marina Hospital in Gaborone, Botswana. The partial core region of HCV was amplified, and genotypes were determined using phylogenetic analysis.
Four genotype 5a and two genotype 4v sequences were identified. Two significant mutations – K10Q and R70Q – were observed in genotype 5a sequences and have been associated with increased risk of hepatocellular carcinoma (HCC), while R70Q confers resistance to interferon-based treatments.
Genotypes 5a and 4v are circulating in Botswana. The presence of mutations in genotype 5 suggests that some patients may not respond to IFN-based regimens. The information obtained in this study, in addition to the World health organization (WHO) recommendations, can be utilized by policy makers to implement DAAs as the new SOC for HCV treatment in Botswana.
Viral hepatitis is the seventh leading cause of death globally , and hepatitis C virus (HCV) is one of the leading causes of liver failure. The prevalence of HCV differs by geographical location. In Africa, a range in prevalence from 1.7 to 14.7% has been reported [2,3,4,5]. In Botswana, preliminary studies reported a low HCV seroprevalence of below 1% [6,7,8] although Botswana is endemic for both human immunodeficiency virus (HIV) and hepatitis B virus (HBV) with the prevalence of HIV/HBV co-infection ranging from 4 to 10.6% [6, 7].
Treatment of HCV has evolved rapidly with the introduction of direct acting antivirals (DAA) in the last decade. DAAs can achieve sustained virologic responses (SVR) of 95 to 100% , have shortened treatment duration, and have overcome many challenges faced by administering interferon-based regimens. The World health organization (WHO) has recommended that DAA regimens be used for HCV treatment, and the choice of antiviral therapy and duration of treatment of the DAAs are dependent on HCV genotype. Therefore, studies of circulating HCV genotypes are important prior to widescale treatment rollout . In Botswana, there are no data on HCV genotypes. This study aimed to determine the HCV genotypes circulating in Botswana and to identify clinically relevant mutations within the HCV core region.
This is a retrospective cross-sectional pilot study utilizing 55 stored plasma samples collected between February 2015 and July 2016 from 55 liver disease patients at the liver clinic of Princess Marina Hospital, a referral hospital in Gaborone, Botswana. The study was conducted at Botswana Harvard AIDS Institute Partnership.
HCV ribonucleic acid (RNA) was extracted using the Qiagen Viral RNA kit using 140 μl plasma samples according to the manufacturer’s specifications (Qiagen, Hilden, Germany). The extracts were stored at − 80 °C prior to genotyping. Amplification targeted the partial core region with two primer sets – outer core primers (5′ – ACT GCC TGA TAG GGT GCT TGC – 3′, nt 288 → 308 and 5′ – ATG TAC CCC ATG AGG TCG GC – 3′, nt 732 ← 751) and inner core primers (5′ – AGG TCT CGT AGA CCG TGC A – 3′, nt 321 → 339 and 5′ – CAT GTG AGG GTA TCG ATG AC – 3′, nt 705 ← 724) previously described, [11, 12] relative to the H77 reference . Nested PCR was carried out using Superscript (III) one-step RT-PCR with Platinum Taq DNA polymerase high fidelity kit (Invitrogen, USA) according to the manufacturer’s instructions. PCR products were visualized on 1.5% agarose gel stained with ethidium bromide, and the positive amplicons were purified using the ZR DNA sequencing clean up kit (Zymo, Irvine, CA, USA) according to the manufacturer’s specifications. Population sequencing was conducted using BigDye Terminator version 3.0 kit (Applied Biosystems; Foster City, CA, USA) sequencing chemistry on a 3130XL ABI genetic analyser (ABI PRISM 3130XL; Applied Biosystems).
Sequence files were edited using Sequencher version 5.0 software (Gene Codes Corp., Ann Arbor, MI, USA). Phylogenetic analysis was used to evaluate HCV genotypes and included reference sequences from the Los Alamos HCV sequence database (https://hcv.lanl.gov/content/index). Sequences were aligned in Clustal X version 2.1. Additional phylogenetic inference was conducted using the Bayesian Markov chain Monte Carlo (MCMC) in the Bayesian Evolutionary Sampling Trees (BEAST) version 1.7.5 as described previously [8, 14]. Posterior probabilities > 80% were deemed statistically significant. Sequences were exported to BioEdit software where nucleotide sequences were translated to amino acids. H77 (accession number AF009606 ) served as the reference sequence, and mutations were visually analysed per amino acid position. Importantly, H77 is a genotype 1a reference; therefore, several genotype 1 references together with genotype 4 and 5 reference sequences were included in the analysis to distinguish between polymorphic regions that differ by genotype and true mutations. Sequences were submitted to National Center for Biotechnology Information (NCBI) GenBank under accession numbers MK392625 to MK392630. Statistical analysis was performed using R version 3.6.0 .
Study participants were Batswana adults and 63.3% were female. The age ranged from 16 years to 74 years with a median of 44 years and an interquartile range (IQR) of 32–55 years. From the 55 participant samples we had access to, 6 (10.9%) of those were RNA positive as shown in Table 1. The six participants had ages, ranging from 24 to 70 years with a median of 55.5 years as shown in Fig. 1a, and HCV was observed in older males and middle-aged women as shown in Fig. 1b.
Sequence analysis of 6 partial core regions revealed that four of the Botswana sequences belong to genotype 5a, as shown by the clustering with other genotype 5a strains from South Africa, Ethiopia, and Denmark. Two Botswana strains belong to genotype 4v as shown by the clustering pattern with other genotype 4v strains from Ethiopia, Cyprus, and the United Kingdom as shown in Fig. 2.
Further analysis was performed to identify core mutations within the viral sequences as shown in Table 2.
Most mutations were observed in genotype 5 samples (R9T, K10Q, T15I, L36 V and R70Q). Only one mutation (G95E) was observed in genotype 4.
This is the first report on the circulating HCV genotypes from Botswana. The finding of genotypes 5a and 4v agrees with previously reported data from southern Africa [5, 26, 27]. Interestingly, genotype 4v has not been observed previously in southern Africa, although other genotypes, such as 4c, 4 g, 4 k, 4q and 4r, have been reported [28, 29]. Genotype 4v has been reported in central Africa and the Middle East. In two studies, genotype 4v was observed in Rwanda and Ethiopia [30, 31]. To achieve elimination of HCV by 2030, DAAs should be introduced and administered to patients as per genotype . According to the observed genotypes, Botswana will need to procure Ledipasvir, Daclatasvir and Sofosbuvir as the DAA regimens recommended to treat genotypes 4 and 5 for patients with and without liver cirrhosis . In this study, the HCV diversity in Botswana is low as compared to South Africa, since only two subtypes were observed. This observation could be due to small sample size and/or varying transmission dynamics.
In the current study the core gene was amplified as the region of interest, despite the recommended NS5B region for HCV genotype classification . We selected a significant fragment of the core region to classify the HCV isolates, which also includes (but is not limited to) putative HCC core related mutations, since the core region has higher amplification rates compared to the NS5B region as previously reported .
The core mutations in genotype 5a – R9T, K10Q, T15I, L36 Vand R70Q – account for most of the observed mutations. The K10Q mutation is associated with increased risk of HCC . This mutation was the nucleotide substitution A28C in the core gene bringing about an amino acid substitution from lysine to glutamine. Further analysis from 100 genotype 5 sequences downloaded from the GenBank showed that only 5% of these references contained the K10Q mutation. Interestingly, all four genotype 5 samples from the Botswana sequences had that mutation.
Several studies in genotype 1 strains have indicated that the R70Q mutation increases the risk of HCC due to its oncogenic effect [18, 23,24,25] and also confers resistance to IFN-based treatments . This mutation occurs due to nucleotide substitution G209A resulting in an amino acid change from arginine to glutamine . From the analysis with the GenBank genotype 5 sequences, it was interesting to note that 81% of the references contained the R70Q mutation compared to 16% with the wildtype amino acid. In a study conducted in South Africa, the R70Q mutation was observed in 90% of genotype 5a blood donors . Whether the impact of this mutation in genotype 1 strains is the same in genotype 5 remains unclear. Presence of this mutation in only genotype 5 isolates could either be due to transmission of a mutated strain or a naturally occurring drug resistant mutation to the SOC . However, the latter assumption cannot be confirmed since there are no data on duration of treatment. Despite the scarce data on impact of this mutation in genotype 5, all four genotype 5 individuals may have a poor response to the current SOC. Thus, DAA regimens may improve treatment of individuals infected with HCV strains harbouring these mutations.
Only one mutation – G95E – was observed in one genotype 4v sample. However, this mutation has not been well characterized or associated with altering viral fitness or drug sensitivity. Most reported mutation analyses have been conducted on genotype 1 strains . Therefore, more longitudinal studies in mutation analysis based on other genotypes, such as 4 and 5, are warranted, as mutations differ by genotype.
A major limitation of the study was the modest population size evaluated and the low amplification success rate; therefore, there was a limited number of positive samples available for analysis. Furthermore, as no HCV antigen enzyme-linked immunosorbent assay (ELISA) was conducted due to low sample volumes, the low amplification rate may reflect a low proportion of current HCV infections versus seropositive but resolved infections. The low amplification rate may also be due to multiple freeze-thaw cycles for sample use in other studies prior to the current analysis. The cross-sectional nature of this study made following up on the patients with HCC-risk factors impossible. Therefore, we could not observe the oncogenic effect of the reported mutations in genotype 5 patients.
In summary, genotypes 5a and 4v are the circulating genotypes in Botswana. The mutations observed in this study confer resistance to the SOC, and, as per WHO recommendations, there is a need to introduce DAAs as the new SOC for Botswana in order to achieve elimination of HCV by 2030. The DAAs for Botswana HCV patients should include Ledipasvir, Daclatasvir and Sofosbuvir. A longitudinal study with a larger representative population is warranted to develop more understanding of the HCV epidemiology in Botswana.
Availability of data and materials
The datasets used to support the results of this study are available from the corresponding author upon request. Sequences are available in the National Center for Biotechnology Information (NCBI) GenBank under accession numbers MK392625 to MK392630.
Direct acting antivirals
Deoxyribose nucleic acid
Enzyme-linked immunosorbent assay
Hepatitis B virus
Hepatitis C virus
Human immune-deficiency virus
Markov chain Monte Carlo
National Center for Biotechnology Information
Reverse transcriptase polymerase chain reaction
Standard of care
Sustained virologic response
World health organisation
Stanaway JD, Flaxman AD, Naghavi M, Fitzmaurice C, Vos T, Abubakar I, et al. The global burden of viral hepatitis from 1990 to 2013: findings from the global burden of disease study 2013. Lancet. 2016;388(10049):1081–8. https://doi.org/10.1016/S0140-6736(16)30579-7.
El-Zanaty F, Way A. Egypt demographic and health survey, 2008. In: Demographic and health survey (EDHS). Cairo: Ministry of Health and Population; 2009.
Razavi H, Waked I, Sarrazin C, Myers RP, Idilman R, Calinas F, et al. The present and future disease burden of hepatitis C virus (HCV) infection with today’s treatment paradigm. J Viral Hepat. 2014;21:34–59.
Agyeman AA, Ofori-Asenso R, Mprah A, Ashiagbor G. Epidemiology of hepatitis C virus in Ghana: a systematic review and meta-analysis. BMC Infect Dis. 2016;16(1):391.
Gower E, Estes C, Blach S, Razavi-Shearer K, Razavi H. Global epidemiology and genotype distribution of the hepatitis C virus infection. J Hepatol. 2014;61(1):S45–57. https://doi.org/10.1016/j.jhep.2014.07.027.
Wester CW, Bussmann H, Moyo S, Avalos A, Gaolathe T, Ndwapi N, et al. Serological evidence of HIV-associated infection among HIV-1-infected adults in Botswana. Clin Infect Dis. 2006;43(12):1612–5.
Patel P, Davis S, Tolle M, Mabikwa V, Anabwani G. Prevalence of hepatitis B and hepatitis C coinfections in an adult HIV centre population in Gaborone, Botswana. Am J Trop Med Hyg. 2011;85(2):390–4. https://doi.org/10.4269/ajtmh.2011.10-0510.
Choga WT, Anderson M, Zumbika E, Moyo S, Mbangiwa T, Phinius BB, Melamu P, Kayembe MK, Kasvosve I, Sebunya TK, Blackard JT, Essex M, Musonda RM, Gaseitsiwe S. Molecular characterization of hepatitis B virus in blood donors in Botswana. Virus Genes. 2018. https://doi.org/10.1007/s11262-018-1610-z.
Van der Meer AJ. Value anti-hepatitis C virus therapy by its clinical efficacy. Hepatology. 2015;62:334–6.
WHO. Global hepatitis report. 2017. Available online: https://www.who.int/hepatitis/publications/global-hepatitis-report2017/en/. Accessed 8 Oct 2018.
Lole KS, Jha JA, Shrotri SP, Tandon BN, Prasad VGM, Arankalle VA. Comparison of hepatitis C virus genotyping by 5 Ј noncoding region- and core-based reverse transcriptase PCR assay with sequencing and use of the assay for determining subtype distribution in India. J Clin Microbiol. 2003;41(11):5240–4.
Cai Q, Zhao Z, Liu Y, Shao X, Gao Z. Comparison of three different HCV genotyping methods: Core, NS5B sequence analysis and line probe assay. Int J Mol Med. 2013;31:347–52. https://doi.org/10.3892/ijmm.2012.1209.
Kuiken C, Combet C, Bukh J, Shin-I T, Deleage G, Mizokami M, et al. A comprehensive system for consistent numbering of HCV sequences, proteins and epitopes. Hepatology. 2006;44(5):1355–61. https://doi.org/10.1002/hep.21377.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/.
Zhang J, Yamada O, Ito T, Akiyama M, Hashimoto Y, Yoshida H, et al. A single nucleotide insertion in the 5′-untranslated region of hepatitis C virus leads to enhanced cap-independent translation. Virology. 1999;261(2):263–70. https://doi.org/10.1006/viro.1999.9879.
Sugiyama K, Suzuki K, Nakazawa T, Funami K, Hishiki T, Ogawa K, et al. Genetic analysis of hepatitis C virus with defective genome and its infectivity in vitro. J Virol. 2009;83(13):6922–8. https://doi.org/10.1128/JVI.02674-08.
Hu Z, Muroyama R, Kowatari N, Chang J, Omata M, Kato N. Characteristic mutations in hepatitis C virus core gene related to the occurrence of hepatocellular carcinoma. Cancer Sci. 2009;100(12):2465–8.
Kuntzen T, Timm J, Berical A, Lennon N, Berlin AM, Young SK, et al. Naturally occurring dominant resistance mutations to hepatitis C virus protease and polymerase inhibitors in treatment-naïve patients. Hepatology. 2008;48(6):1769–78. https://doi.org/10.1002/hep.22549.
Bull RA, Eltahla AA, Rodrigo C, Koekkoek SM, Walker M, Pirozyan MR, Luciani F. A method for near full-length amplification and sequencing for six hepatitis C virus genotypes. BMC Genomics. 2016;17(1):247.
Sakamoto M, Akahane Y, Tsuda F, Tanaka T, Woodfield DG, Okamoto H. Entire nucleotide sequence and characterization of a hepatitis C virus of genotype V/3a. J Gen Virol. 1994;75(7):1761–8.
Akuta N, Suzuki F, Hirakawa M, Kawamura Y, Yatsuji H, Sezaki H, Suzuki Y, Hosaka T, Kobayashi M, Kobayashi M, Saitoh S, Arase Y, Ikeda K, Chayama K, Nakamura Y, Kumada H. Amino acid substitution in hepatitis C virus core region and genetic variation near the interleukin 28B gene predict viral response to telaprevir with peginterferon and ribavirin. Hepatology. 2010;52:421–9.
Nakamoto S, Imazeki F, Fukai K, Fujiwara K, Arai M, Kanda T, Yonemitsu Y, Yokosuka O, Nakamoto S, Imazeki F, Fukai K, Fujiwara K, Arai M, Kanda T, Yonemitsu Y, Yokosuka O. Association between mutations in the core region of hepatitis C virus genotype 1 and hepatocellular carcinoma development. J Hepatol. 2010;52(1):72–8 Erratum in: J Hepatol, 52(4):620.
Fishman SL, Factor SH, Balestrieri C, Fan X, Dibisceglie AM, Desai SM, Benson G, Branch AD. Mutations in the hepatitis C virus core gene are associated with advanced liver disease and hepatocellular carcinoma. Clin Cancer Res. 2009;15(9):3205–13.
El-Shamy A, Eng FJ, Doyle EH, Klepper AL, Sun X, Sangiovanni A, Iavarone M, Colombo M, Schwartz RE, Hoshinda Y, Branch AD. A cell culture system for distinguishing hepatitis C viruses with or without liver cancer-related mutations in the viral core gene. J Hepatol. 2015. https://doi.org/10.1016/j.jhep.2015.07.024.
Petruzziello A, Marigliano S, Loquercio G, Cozzolino A, Cacciapuoti C. Global epidemiology of hepatitis C virus infection: an up-date of the distribution and circulation of hepatitis C virus genotypes. World J Gastroenterol. 2016;22(34):7824–40. https://doi.org/10.3748/wjg.v22.i34.7824.
Gededzha MP, Selabe SG, Blackard JT, Kyaw T, Mphahlele MJ. Near full-length genome analysis of HCV genotype 5 strains from South Africa. Infect Genet Evol. 2014;21(November):118–23. https://doi.org/10.1016/j.meegid.2013.10.022.
Gededzha MP, Selabe SG, Kyaw T, Rakgole JN, Blackard JT, Mphahlele MJ. Introduction of new subtypes and variants of hepatitis C virus genotype 4 in South Africa. J Med Virol. 2012;84:601–7. https://doi.org/10.1002/jmv.23215.
Gededzha MP, Selabe SG, Kyaw T, Rakgole JN, Blackard JT, Mphahlele MJ. Introduction of new subtypes and variants of hepatitis C virus genotype 4 in South Africa. J Med Virol. 2012;84(4):601–7.
Iles JC, Raghwani J, Harrison GL, Pepin J, Djoko CF, Tamoufe D, LeBreton M, Schneider BS, Fair JN, Tshala FM, Kayembe PK, Muyembe JJ, Edidi-Basepeo S, Wolfe ND, Simmonds P, Klenerman P, Pybus OG. Phylogeography and epidemic history of hepatitis C virus genotype 4 in Africa. Virology. 2014;464-465:233–43.
Hundie GB, Raj VS, GebreMichael D, Pas SD, Haagmans BL. Genetic diversity of hepatitis C virus in Ethiopia. PLoS One. 2017;12(6):0179064.
Quer J, Gregori J, Rodríguez-Frias F, Buti M, Madejon A, Perez-del-Pulgar S, Garcia-Cehic D, Casillas R, Blasi M, Homs M, Tabernero D, Alvarez-Tejado M, Muñoz JM, Cubero M, Caballero A, delCampo JA, Domingo E, Belmonte I, Nieto L, Lens S, Muñoz-de-Rueda P, Sanz-Cameno P, Sauleda S, Bes M, Gomez J, Briones C, Perales C, Sheldon J, Castells L, Viladomiu L, Salmeron J, Ruiz-Extremera A, Quiles-Pérez R, Moreno-Otero R, López-Rodríguez R, Allende H, Romero-Gómez M, Guardia J, Esteban R, Garcia-Samaniego J, Forns X, Esteban JI. High- resolution hepatitis C virus subtyping using NS5B deep sequencing and phylogeny, an alternative to current methods. J Clin Microbiol. 2015;53:219–26. https://doi.org/10.1128/JCM.02093-14.
Prabdial-Sing N, Blackard JT, Puren AJ, Mahomed A, Abuelhassan W, Mahlangu J, Vermeulen M, Bowyer SM. Naturally occurring resistance mutations within the core and NS5B regions in hepatitis C genotypes, particularly genotype 5a, in South Africa. Antivir Res. 2016;127:90–8.
The authors would like to acknowledge the study participants. We also acknowledge the University of Botswana, Princess Marina Hospital, National University of Science and Technology in Zimbabwe and the Botswana Harvard AIDS Institute Partnership for their support and contributions to the study.
MA, WTG, SM and SG were funded by the Sub-Saharan African Network for TB/HIV Research Excellence (SANTHE), a DELTAS Africa Initiative [grant # DEL-15-006]. The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Wellcome Trust [grant # 107752/Z/15/Z] and the UK government. The views expressed in this publication are those of the author(s) and not necessarily those of AAS, NEPAD Agency, Wellcome Trust or the UK government. SM was partially supported was supported by the Fogarty International Center and National Institute of Mental Health, of the National Institutes of Health under Award Number D43 TW010543. The opinions, findings and conclusions expressed in this article are those of the authors and do not necessarily reflect the views of the sponsors.
Ethics approval and consent to participate
The study was approved by the Health Research and Development Committee (HRDC) at the Ministry of Health and Wellness (HPDME 13/18/1 X (897). Permission to access raw data was granted by the ethics committee of the university of Botswana and Princess Marina hospital. All participants consented prior to participation in the study.
Consent for publication
SM is an editorial board member for BMC Infectious disease. The remaining authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Bhebhe, L., Anderson, M., Souda, S. et al. Molecular characterization of hepatitis C virus in liver disease patients in Botswana: a retrospective cross-sectional study. BMC Infect Dis 19, 875 (2019). https://doi.org/10.1186/s12879-019-4514-1