# Endemicity and genetic diversity of Hepatitis delta virus among Pygmies in Cameroon, Central Africa

Yacouba Foupouapouognigni, Jacques Delors Toumansie Mfonkou, Onana Boyomo, Antoine Gessain, Richard Njouom

<figure><img src="image_1.png" /></figure>

## To cite this version:

Yacouba Foupouapouognigni, Jacques Delors Toumansie Mfonkou, Onana Boyomo, Antoine Gessain, Richard Njouom. Endemicity and genetic diversity of Hepatitis delta virus among Pygmies in Cameroon, Central Africa. BMC Research Notes, 2022, 15 (1), pp.87. 10.1186/s13104-022-05976-7. pasteur-04079391

HAL Id: pasteur-04079391

https://pasteur.hal.science/pasteur-04079391

Submitted on 24 Apr 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH NOTE

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Endemicity and genetic diversity of Hepatitis delta virus among Pygmies in Cameroon, Central Africa

Yacouba Foupouapouognigni³, Jacques Delors Toumansie Mfonkou<fnref n="1" />,², Onana Boyomo², Antoine Gessain⁴ and Richard Njouom<fnref n="1" />* ID

## Abstract

**Objective:** A single study conducted about three decades ago on hepatitis D virus (HDV) infection among Baka pygmies in Cameroon reported a very high anti-HDV antibodies prevalence of 46%, but HDV genetic diversity has not been studied in this population. The genetic diversity of strains from endemic ancient populations may help to understand the origin and evolutionary history of viruses. This study aimed to investigate the HDV seroprevalence and the genetic diversity in three remote Cameroonian Pygmies with chronic HBV infection.

**Results:** An unusually high 69% (36/52) level of HDV infection was found among HBsAg-positive pygmies in Cameroon. HDV RNA was detected and sequenced in 38.8% (14/36). The phylogenetic analysis revealed that 9/14 strains (64.3%) were identified and classified as genotype 1 (HDV-1) and 5/14 (35.6%) as genotype 7 (HDV-7), respectively with a bootstrap value of 100%. The further analysis showed the co-circulation of highly diverse HDV genotypes HDV-1 and HDV-7 in this population. These results highlight the endemicity of HDV infection in Central Africa. The highly diverse HDV-1 and HDV-7 in pygmies suggest an African origin of HDV. However, further studies are needed with larger sample size.

**Keywords:** Hepatitis D, Genetic diversity, Pygmies, Cameroon

## Introduction

Hepatitis delta virus (HDV) is necessarily associated with hepatitis B virus (HBV) infection and tends to cause more severe disease than HBV infection alone [1]. HDV infection occurs worldwide, but its prevalence largely varies between geographic regions. It is highly endemic in the Middle East, the Mediterranean region, the Amazon zone bassin, and several African countries [2, 3]. Worldwide, over eight phylogenetically distinct genotypes have been revealed from varied regions. Except the ubiquitous genotype 1 (HDV-1), each virus clade has its distinct

geographic distribution: HDV-2 is found in Japan, Taiwan, and Russia; HDV-4 in Taiwan and Japan; HDV-3 in the Amazon zone bassin; and HDV-5, HDV-6, HDV-7 and HDV-8 in Africa [3, 4]. In a recent large representative sample of the adult population of Cameroon, we found an extreme degree of heterogeneity of HDV distribution across Cameroon, ranging from 1 to 54% according to regions, with the highest prevalence found in the Southern and Eastern regions close to the Equator [5]. A single study conducted about three decades ago on HDV infection among Baka Pygmies in Cameroon reported a very high anti-HDV antibodies prevalence of 46%, but HDV genetic diversity has not been studied in this population [6]. In a previous study conducted in 2011 in Bantou population of Cameroon, we found a high genetic diversity of HDV in Cameroon with the circulation of at

<footnote>*Correspondence: njouom@pasteur-yaounde.org; njouom@yahoo.com</footnote>

<footnote n="1">¹Virology Service, Centre Pasteur of Cameroon, BP 1274, Yaoundé, Cameroon</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
least 4 genotypes including genotypes 1, 5, 6 and 7 [6]. In a recent study, the HDV genetic diversity was studied using 211 HDV-positive samples obtained from blood donors in retrospective studies conducted between 2010 and 2016 in Douala and Yaounde in Cameroon as well as patients with illness of unknown etiology, antenatal patients, and participants of door-to-door and voluntary testing campaigns recruited from 14 villages and towns in the South region of Cameroon [7]. HDV genotypes 1, 6, 7, and 8 were identified in the study population, indicating a high level of HDV diversity in Cameroon. Genotype 1 was predominant (65.4%), followed by genotype 7 (28.9%), genotype 6 (5.2%) and one sample of genotype 8 (0.5%). Potential recombinant breakpoints as depicted by sequence branching basal to references were assayed with the aid of SimPlot, however, none was found. The genetic diversity of strains from endemic ancient populations may help to understand the origin and evolutionary history of viruses [8]. In the present study, we investigated the HDV seroprevalence and the genotype diversity in three remote Cameroonian Pygmy populations with chronic HBV infection.

## Main text

### Methods

This study was part of a survey on emerging viruses in Pygmies from Cameroon conducted from 2005 to 2008 [9–11]. The geographic localization and the characteristics of the subjects included in this study are described elsewhere [9]. After obtained consent, whole blood was collected in EDTA tube and sample obtained was centrifuged at 2000g during 5 min. Plasma isolated was aliquoted and stored at –80 °C for subsequent handling. The Murex anti-delta assay (Abbott, Wiesbaden, Germany) was employed as described in the manufacturer’s protocol to determine the presence of HDV antibodies (HDV-Ab). In order to characterize HDV samples, RNA was extracted from 140 µl HDV-Ab-positive plasma using the QIAamp RNA minikit (Qiagen, Courtaboeuf, France) and HDV RNA was detected by reverse transcription (RT)-PCR amplification of a 360-bp fragment coding for the small hepatitis delta antigen (HDV-sHD) gene as described previously [12]. Amplicons were sequenced bidirectionally with the BigDye Terminator v3.1 Sequencing Kit (Applied Biosystems) using primers described previously [12]. Sequences editing was done using EDITSEQ (software Seqman™ II Lasergene). HDV reference sequences were retrieved from the GenBank database for multi-alignment using the CLUSTAL W method [13] in the MEGA software version 6 [14]. Originally, we used the neighbour-joining method with 1000-bootstrap replicates to generate an extended phylogenetic tree with all known HDV genotypes. Then,

the HDV sequence data sets were downsized by omitting most branches containing no studied sequences. For the analysis of the maximum likelihood of phylogenetic relationship, a reduced data set was used whilst considering the best-fit evolution model. In this case, the General Time Reversible with Gamma distribution GTR+G+I) was used with the assumption that a certain fraction of sites was invariable on an evolutionary level (GTR+G+I). The 1000 bootstrap resampling method was used to assess the reliability of the tree nodes. Kimura two-parameter approach was used to determine the genetic distances between the strains.

For this study, 52 of 62 HBsAg-positive pygmies previously detected were available [9].

### Results

HDV-Ab were found in 36 (69.2%) of the 52 tested plasmas. There was a significant association of HDV-Ab prevalence with age, increasing from 36.6% in children below 10 years of age to 46.2% in 10–20 years old adolescents, and >75% among the >21 years old (P < 0.001) (Table 1). RNA was detected and sequenced only in 14 out of the 36 HDV-Ab positive samples. The obtained sequences were compared to HDV reference sequences retrieved from Genbank corresponding to the eight HDV genotypes. The phylogenetic analysis revealed that 9 (64.3%) out of the 14 included in this study were classified as genotype 1 (HDV-1) with a bootstrap value of 100% (Fig. 1). The pairwise comparison of these nine strains revealed 81% to 87.1% similarity with a consensus HDV-1

<table><caption>**Table 1** Seroprevalence of HDV-Ab in three pygmy groups of Cameroon infected with HBV</caption><thead><tr><th rowspan="2">Age</th><th colspan="4">n/N<br>(%; IC95%)</th></tr><tr><th>Baka</th><th>Bakola</th><th>Bedzam</th><th>Total</th></tr></thead><tbody><tr><td>≤ 10</td><td>1/3<br>(33.3; -20-86.6)</td><td>3/8<br>(37.5; 4-71)</td><td>0</td><td>4/11<br>(36.6; 8.1-65.1)</td></tr><tr><td>11-20</td><td>4/7<br>(57.1; 20.5-93.7)</td><td>1/5<br>(20; -15.1-55.1)</td><td>1/1<br>(100)</td><td>6/13<br>(46.2; 19.1-73.3)</td></tr><tr><td>21-30</td><td>8/9<br>(88.8; 68.2-109.4)</td><td>0</td><td>2/2<br>(100)</td><td>10/11<br>(90.9; 73.9-107.9)</td></tr><tr><td>31-40</td><td>6/6<br>(100)</td><td>0</td><td>3/3<br>(100)</td><td>9/9<br>(100)</td></tr><tr><td>41-50</td><td>4/4<br>(100)</td><td>0</td><td>0</td><td>4/4<br>(100)</td></tr><tr><td>> 50</td><td>3/4<br>(75; 32.6-117.6)</td><td>0</td><td>0</td><td>3/4<br>(75; 32.6-117.6)</td></tr><tr><td>Total</td><td>26/33<br>(78.8; 64.9-92.7)</td><td>4/13<br>(30.8; 5.7-55.9)</td><td>6/6<br>(100)</td><td>36/52<br>(69.2; 56.2-81.7)</td></tr></tbody></table>

n: number of HDV-Ab positive samples; N: number of samples tested; %: rate of seroprevalence; IC 95%: 95% confidence interval
<figure><img src="image_4.png" /></figure>

reference sequence. The percentage of similarity was less than 71.9% when compared with HDV-2 to HDV-8 reference sequences (Table 2). Thus, these 9 strains were

assigned to HDV-1 genotype. The phylogenetic analyses of the Cameroonian HDV-1 strains identified three clusters (Fig. 1). The first cluster containing three strains
<table><caption>**Table 2** Percent similarity between the R0 gene nucleotide sequences of the 14 HDV samples from Cameroonian Pygmy populations and the representatives of the 8 reference HDV genotypes</caption><thead><tr><th></th><th>HDV-1</th><th>HDV-2</th><th>HDV-3</th><th>HDV-4</th><th>HDV-5</th><th>HDV-6</th><th>HDV-7</th><th>HDV-8</th></tr></thead><tbody><tr><td>PYG70</td><td>82.6</td><td>71</td><td>60.2</td><td>71.1</td><td>65.9</td><td>66.4</td><td>66.4</td><td>67.1</td></tr><tr><td>R0Bobak-115</td><td>83</td><td>71.9</td><td>60.2</td><td>71.2</td><td>66.5</td><td>66.9</td><td>66.6</td><td>67.3</td></tr><tr><td>Mebak-58</td><td>81.6</td><td>70.7</td><td>60.4</td><td>70.3</td><td>67.5</td><td>67</td><td>69.5</td><td>68.4</td></tr><tr><td>Lobak-47</td><td>85</td><td>71.3</td><td>61.5</td><td>70.6</td><td>68</td><td>66.9</td><td>68</td><td>68.8</td></tr><tr><td>R0PYM16</td><td>87.1</td><td>72.6</td><td>61.8</td><td>71.4</td><td>69</td><td>68.3</td><td>68.1</td><td>68.6</td></tr><tr><td>R0PYM23</td><td>82.9</td><td>70.7</td><td>60.5</td><td>68.9</td><td>66.5</td><td>66.3</td><td>66</td><td>67.2</td></tr><tr><td>R0Bobak-118</td><td>85.5</td><td>71.4</td><td>61.1</td><td>71.7</td><td>67.4</td><td>67.7</td><td>67</td><td>68.1</td></tr><tr><td>R0Bobak-66</td><td>85.3</td><td>71.8</td><td>61.6</td><td>71.6</td><td>67.2</td><td>67.6</td><td>67</td><td>67.8</td></tr><tr><td>R0Mebak-61</td><td>84</td><td>70.4</td><td>60.3</td><td>70.1</td><td>66.3</td><td>66.6</td><td>66.2</td><td>67.2</td></tr><tr><td>R0PYG-56</td><td>68</td><td>67.9</td><td>59.6</td><td>69.6</td><td>69.9</td><td>67.7</td><td>86.8</td><td>73.6</td></tr><tr><td>Bobak95</td><td>69</td><td>67.9</td><td>59.5</td><td>69.9</td><td>68.5</td><td>69</td><td>81</td><td>74.6</td></tr><tr><td>Mebak57</td><td>65.5</td><td>64.8</td><td>57.4</td><td>66.4</td><td>67.2</td><td>66.5</td><td>77.4</td><td>70.7</td></tr><tr><td>R0PYG65</td><td>59.6</td><td>59.9</td><td>54.2</td><td>62.6</td><td>63</td><td>61.4</td><td>72.2</td><td>65.8</td></tr><tr><td>R0PYG48</td><td>59.8</td><td>60</td><td>55.3</td><td>62.4</td><td>64.4</td><td>62.5</td><td>73.9</td><td>67.4</td></tr></tbody></table>

(PYG70, Mebak58, and R0Bobak115) was supported by a bootstrap value of 98% and a similarity score of more than 81% was observed between these strains (Table 2). This cluster grouped with two previously described HDV sub-genotype 1b strains from Cameroon and Central African Republic (Fig. 1). The second cluster containing three strains (Lobak47, R0PYM23, and R0PYM16) was supported by a bootstrap value of 85% (Fig. 1) and a similarity score of more than 82% was observed between these strains (Table 2). This second cluster grouped with five previously described HDV sub-genotype 1a strains from Cameroon (Fig. 1). Finally, the third cluster containing three strains (R0Bobak66, R0Bobak118, and R0Mebak61) was supported by a bootstrap value of 86% (Fig. 1) and a similarity score of more than 84% was observed between these strains (Table 2). This third cluster did not group with a significant bootstrap with any previously described HDV sub-genotype (Fig. 1). The average nucleotide similarity (83.24%) among the Cameroonian HDV-1 strains is in the same range with the 83% average reported among other worldwide HDV-1 strains [4]. These results also highlight the intra-genotypic diversity of the Cameroonian HDV-1 isolates with the presence of two clearly defined HDV sub-genotypes 1a and 1b and the circulation of three new possible HDV-1 sub-genotypes not yet assigned. The complete HDV nucleotide sequences of these three samples are thus required for a correct assignment. Five strains (35.71%) out of the 14 included in this study (R0PYG65, R0PYG48, R0PYG56, Mebak57, and Bobak95), clustered with previously described HDV genotype 7 (HDV-7) with a 100% bootstrap value (Fig. 1). The pairwise comparison of these five strains revealed a 72.2% to 86.8% similarity

with a consensus HDV-7 reference sequence (Table 2). Phylogenetically, these five HDV Cameroonian strains formed three clusters (Fig. 1). The first contained one strain (R0PYG56) with a similarity score of 86.8% when compared to previously described HDV-7 genotype, suggesting that this strain could be assigned to HDV-7 genotype. The second contained two strains (R0PYG65 and R0PYG48), was supported by a bootstrap value of 100% and a similarity score of less than 74% when compared to all 8 previously described HDV genotypes, suggesting that they are distinct from any of the 8 previously described HDV genotypes. The third group contained two strains (Bobak95 and Mebak57) was supported by a bootstrap value of 99% and presented a similarity score of 77.4% and 81% with HDV-7 reference sequence while it was less than 74.6% when compared with the other 7 representatives reference sequences suggesting that they are different from all the 8 previously described HDV genotypes. Although these results do not fulfil the recommendations for a designation of a new genotype (that is, ≥ 3 distinct isolates with high bootstrap values and showing high scores of similarity [4], it is likely that the above strains belong to one or even 2 new HDV genotypes. Further studies in pygmy populations are necessary to get the complete HDV nucleotide sequences of these samples for a correct assignment and to confirm if they could be defined as new HDV genotypes.

## Discussion

Our results point to an endemicity of HDV infection in the Cameroonian pygmy populations which is in sharp accordance with data reported from the same population about 30 years ago and to the situation in neighbouring
Central African countries. Additionally, HDV infection in Cameroon is characterized by a wide genetic diversity with more than two different genotypes in circulation. This study allows us to highlight possible new HDV genotypes or sub-genotypes obtained from strains isolated from Cameroonian Pygmy populations.

## Conclusion

Our results show an endemic circulation of HDV within the pygmy populations of Cameroon. The description of highly diverse HDV-1 and HDV-7 (African genotype) in this ancient population suggest an African origin of HDV. However, further studies are needed in these groups with larger sample size.

## Limitations

The study of genetic diversity of HDV was carried out on a short sequence (RO). The complete genome would be necessary for assignment of sub-genotype or possibly recombinant forms.

## Abbreviations

HDV: Hepatitis D virus; HBsAg: Hepatitis B surface antigen; ELISA: Enzyme linked immunosorbent assay.

## Acknowledgements

The authors are grateful to the members of the three Pygmy ethnic groups who agreed to participate in this study.

## Authors' contributions

RN and AG designed and acquired funding for this study. YF performed the experiments, drafted, reviewed and edited the manuscript. JDTM performed genetic analysis, reviewed and edited the manuscript. RN, OB and AG drafted, reviewed and edited the manuscript. All authors read and approved the final manuscript.

## Funding

This study was funded by the Agence Nationale de Recherche sur le Sida et les Hépatites Virales (Grant ANRS 12366).

## Availability of data and materials

Nucleotide sequences obtained in the present study were submitted to Gen-Bank under the registration numbers MW773622 to MW773635.

## Declarations

### Ethics approval and consent to participate

This study was a part of survey which was approved by the Cameroonian National Ethics Committee and the Ministry of Public Health. Prior to blood specimen collection, written informed consent was obtained from each participating adult (or parents, in case of minors).

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Author details

I, BP 812, Yaoundé, Cameroon. ⁴Unit of Epidemiology and Physiopathology of Oncogenic Viruses, URA CNRS 3015, Department of Virology, Institut Pasteur, 25-28 rue du Dr Roux, 75724 Paris Cedex 15, France.

Received: 22 October 2021 Accepted: 16 February 2022

Published online: 03 March 2022

## References

1. Botelho-Souza LF, Vasconcelos MPA, de Dos Santos AO, et al. Hepatitis delta: virological and clinical aspects. Virol J. 2017;14:177. https://doi.org/10.1186/s12985-017-0845-y.

2. Miao Z, Zhang S, Ou X, et al. Estimating the global prevalence, disease progression, and clinical outcome of hepatitis delta virus infection. J Infect Dis. 2020;221:1677-87. https://doi.org/10.1093/infdis/jiz633.

3. Stockdale AJ, Kreuels B, Henrion MYR, et al. The global prevalence of hepatitis D virus infection: systematic review and meta-analysis. J Hepatol. 2020;73:523-32. https://doi.org/10.1016/j.jhep.2020.04.008.

4. Le Gal F, Brichler S, Drugan T, et al. Genetic diversity and worldwide distribution of the deltavirus genus: a study of 2,152 clinical strains. Hepatology. 2017;66:1826-41. https://doi.org/10.1002/hep.29574.

5. Besombes C, Njouom R, Paireau J, et al. The epidemiology of hepatitis delta virus infection in Cameroon. Gut. 2020;69:1294-300. https://doi.org/10.1136/gutjnl-2019-320027.

6. Foupouapouognigni Y, Noah DN, Sartre MT, Njouom R. High prevalence and predominance of hepatitis delta virus genotype 1 infection in Cameroon. J Clin Microbiol. 2011;49:1162-4. https://doi.org/10.1128/JCM.01822-10.

7. Butler EK, Rodgers MA, Coller KE, et al. High prevalence of hepatitis delta virus in Cameroon. Sci Rep. 2018;8:11617. https://doi.org/10.1038/s41598-018-30078-5.

8. Kramvis A, Kew MC. Epidemiology of hepatitis B virus in Africa, its genotypes and clinical associations of genotypes. Hepatol Res. 2007;37:59-19. https://doi.org/10.1111/j.1872-034X.2007.00098.x.

9. Foupouapouognigni Y, Serge Alain S, Betsem à Betsem E, et al. Hepatitis B and C virus infections in the three Pygmy groups in Cameroon. J Clin Microbiol. 2011;49:737-40. https://doi.org/10.1128/JCM.01475-10.

10. Calattini S, Betsem E, Bassot S, et al. New strain of human T lymphotropic virus (HTLV) type 3 in a Pygmy from Cameroon with peculiar HTLV serologic results. J Infect Dis. 2009;199:561-4. https://doi.org/10.1086/596206.

11. Calattini S, Betsem E, Bassot S, et al. Multiple retroviral infection by HTLV type 1, 2, 3 and simian foamy virus in a family of Pygmies from Cameroon. Virology. 2011;410:48-55. https://doi.org/10.1016/j.virology.2010.10.025.

12. Ivaniushina V, Radjef N, Alexeeva M, et al. Hepatitis delta virus genotypes I and II cocirculate in an endemic area of Yakutia, Russia. J Gen Virol. 2001;82:2709-18. https://doi.org/10.1099/0022-1317-82-11-2709.

13. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673-80. https://doi.org/10.1093/nar/22.22.4673.

14. Tamura K, Stecher G, Peterson D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725-9. https://doi.org/10.1093/molbev/mst197.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<footnote n="1">¹Virology Service, Centre Pasteur of Cameroon, BP 1274, Yaoundé, Cameroon.</footnote>

<footnote n="2">²Laboratory of Microbiology, University of Yaounde I, BP 812, Yaoundé, Cameroon.</footnote>

<footnote n="3">³Laboratory of Pharmacology and Toxicology, University of Yaoundé.</footnote>