Molecular Epidemiology of Hepatitis B Virus Variants in Nonhuman Primates

JOURNAL OF VIROLOGY, June 2000, p. 5377 5381 Vol. 74, No. 11 0022-538X/00/$04.00 0 Copyright 2000, American Society for Microbiology. All Rights Reserved. Molecular Epidemiology of Hepatitis B Virus Variants in Nonhuman Primates STEFANIE GRETHE, 1 JENS-OVE HECKEL, 2 WOLFRAM RIETSCHEL, 2 AND FRANK T. HUFERT 3 * Department of Medical Microbiology, University of Göttingen, D-37075 Göttingen, 1 Zoological and Botanical Garden Wilhelma, D-70342 Stuttgart, 2 and Department of Virology, Institute for Medical Microbiology and Hygiene, University of Freiburg, D-79104 Freiburg, 3 Germany Received 21 October 1999/Accepted 28 February 2000 We characterized hepatitis B virus (HBV) isolates from sera of 21 hepatitis B virus surface antigen-positive apes, members of the families Pongidae and Hylobatidae (19 gibbon spp., 1 chimpanzee, and 1 gorilla). Sera originate from German, French, Thai, and Vietnamese primate-keeping institutions. To estimate the phylogenetic relationships, we sequenced two genomic regions, one located within the pre-s1/pre-s2 region and one including parts of the polymerase and the X protein open reading frames. By comparison with published human and ape HBV isolates, the sequences could be classified into six genomic groups. Four of these represented new genomic groups of gibbon HBV variants. The gorilla HBV isolate was distantly related to the chimpanzee isolate described previously. To confirm these findings, the complete HBV genome from representatives of each genomic group was sequenced. The HBV isolates from gibbons living in different regions of Thailand and Vietnam could be classified into four different phylogenetically distinct genomic groups. The same genomic groups were found in animals from European zoos. Therefore, the HBV infections of these apes might have been introduced into European primate-keeping facilities by direct import of already infected animals from different regions in Thailand. Taken together, our data suggest that HBV infections are indigenous in the different apes. One event involving transmission between human and nonhuman primates in the Old World of a common ancestor of human HBV genotypes A to E and the ape HBV variants might have occurred. Human hepatitis B virus (HBV), the prototype member of the genus Hepadnaviridae, can infect different primates. In nature, HBV infections occur mainly in humans, although they have been also documented in chimpanzees (16), gibbons (8, 10), and orangutans (17, 18). Other members of the family Orthohepadnaviridae infect rodents, e.g., woodchucks (3) and ground squirrels (14, 15). A new hepadnavirus was recently isolated from a woolly monkey, a New World primate (woolly monkey HBV [WMHBV]) (7). This virus is distantly related to the human HBV group and is therefore considered to represent a progenitor of the human viruses. HBV shows a remarkable genetic variability, which is reflected in the occurrence of at least six human HBV genotypes and two nonhuman primate HBV genotypes, consisting of a chimpanzee HBV isolate and a gibbon HBV isolate (9, 10). All these primate HBV genotypes are defined by an intergenotype variation of 8% based on complete viral genomes (12). The different human HBV genotypes show characteristic geographical distributions (11). There are only a few reports about the occurrence of HBV infections in ape populations (6, 17). It is still a matter of debate whether these infections are of human origin or indigenous to the different apes (7, 17). At present, the complete HBV sequence from ape isolates has been analyzed only for one chimpanzee and one gibbon (10, 16). These viruses might as well represent human viruses introduced artificially into * Corresponding author. Mailing address: Department of Virology, Institute for Medical Microbiology and Hygiene, University of Freiburg, Hermann-Herder-Str. 11, D-79104 Freiburg, Germany. Phone: 49-761-203.6591. Fax: 49-761-203.6603. E-mail: hufert@ukl.uni-freiburg.de. Present address: Landeskriminalamt Rheinland-Pfalz, D-55118 Mainz, Germany. these animals (7). In order to answer this question, we analyzed in this study HBV isolates from HBV surface antigen (HBsAg)-positive apes kept in different primate-keeping facilities in Germany, France, Thailand, and Vietnam. A correlation between the geographic origins of the apes and their HBV sequences was investigated by analyzing parts of the viral genome. In addition, we tried to find chains of HBV transmission within the different primate-keeping institutions. The sequences of two short HBV genomic regions were analyzed to characterize phylogenetically closely related viral isolates. One region (A) encompasses the pre-s1 genomic region, and the second region (B) includes an area that codes for the viral polymerase and a part of the X protein. However, conclusions concerning possible new genotypes can only be drawn on the basis of complete genomes as defined by Okamoto et al. (12). Thus, the complete genomes of representatives of each genomic group were then fully sequenced in order to analyze the intergenotype variation. Sera from 19 gibbon spp., 1 chimpanzee, and 1 gorilla were part of the primate serum bank stored at the Zoological and Botanical Garden Wilhelma, Stuttgart, Germany. Sera were collected during routine veterinary procedures in different zoos and primate-keeping facilities in Europe, Thailand, and Vietnam (Table 2). All sera were positive for at least one of the HBV serologic parameters HBsAg, HBeAg, anti-hbs, and anti-hbc, as determined by microparticle enzyme immunoassays using commercial test kits (IMx-HBsAg, IMs-HBe2, IMx- AUSAB, and IMx-CORE; Abbott Laboratories). HBV DNA, detected by PCR, was present in all samples. HBV DNA was extracted from serum as described previously (5). Primers used for amplification of overlapping fragments spanning the HBV genome are given in Table 1. Nested PCR was carried out (5) for 35 cycles, with denaturation at 95 C for 15 s, annealing at 57 C for 15 s, and extension at 72 C 5377

5378 NOTES J. VIROL. TABLE 1. Primers used for amplification and sequencing a Fragments b First-round amplification Oligonucleotides used for: Second-round amplification Sequencing primers c 1 3207-21f/1276-22r 57-22f/1196-22r 57-22f, 433-22r, 412-22f, 720-21r, 703-20f, 1196-22r 2 2360-21f/477-23r 2382-21f/433-22r 2382-21f, 2850-22r, 57-22f, 433-22r 2724-25f, 3221-21r (region A) 3 1394-23f/2441-22r 1425-23f/2441-22r 1425-23f, 1789-20r, 1785-24f, 1959-23r, 2045-23f, 2441-22r 3.I 1099-23f/1789-20r 1253-21f/1789-20r 1253-21f, 1425-23f, 1789-20r 3.II 1785-24f/2850-22r 1785-24f/2441-22r 1785-24f, 1959-23r, 2045-23f, 2441-22r 4 685-20f/1701-20r 703-20f/1673-24r 703-20f, 1196-22r 1099-23f, 1673-24r (region B) a Sequences of oligonucleotides: 57-22f, CTGCTGGTGGCTCCAGTTCAGG; 703-20f, CGTAGGGCTTTCCCCCACTG; 1099-23f, TCGCCAACTTACAAGGC CTTTCT; 1196-22r, GTTGCGTCAGCAAACACTTGGC; 1253-21f, CTCCTCTGCCGATCCATACTG; 1785-24f, GGCATAAATTGGTCTGCGCACCAG; 1789-20r, TGCCTACAGCCTCCTAGTAC; 1959-23r, GGCAAAAAAGAGAGTAACTCCAC; 2441-22r, GAGATTGAGATCTTCTGCGACG. Sequences of the other oligonucleotides are given in reference 5. Regions 2 and 4 were amplified and regions A and B were analyzed for all HBV isolates. Amplification and sequence analysis of the whole genome were done with representatives of each genomic group (see Table 2). b Fragments 3.I and 3.II were analyzed when amplification of fragment 3 (spanning the nick of the HBV genome) failed because of low serum HBV DNA concentration. c Primers in boldface were used for amplification of regions A and B, respectively. for 70 s, as described previously. Fragments 2 and 4 were amplified from all samples to determine genomic groups of regions A and B. We were able to amplify these target sequences from almost all ape sera with the exception of that from gibbon 174. Here, we could amplify fragment 2 only. Fragments 1 and 3 were amplified additionally from one representative of each genomic group to sequence the complete genome. To determine the phylogenetic relationships, we sequenced two genomic regions, one located within the pre-s1/pre-s2 region (region A) and one including parts of the polymerase and the X protein open reading frame (region B). Nucleotide sequences were determined for both strands with the PRISM Ready Reaction DyeDeoxy Terminator Cycle sequencing kit on an ABI 373A DNA sequencer (Applied Biosystems, Weiterstadt, Germany) using PCR primers and internal oligonucleotides (Table 1). Sequence analysis was performed using the program package GCG (Wisconsin sequence analysis package, version 8.1, University of Wisconsin, Madison). Some phylogenetic analy- TABLE 2. Survey of the geographic origins of the HBV-infected apes and the phylogenetic classification of the corresponding virus isolates, based on sequence analysis of fragment A a Species Isolate b Genomic group Origin Western lowland gorilla (Gorilla gorilla gorilla) gor97 chimp (6% divergence) Leipzig, Germany (wild-caught animal from Cameroon) Chimpanzee (Pan troglodytes) chimp82 Welzheim, Germany White-handed gibbon (Hylobates lar) gib644 I c Wuppertal, Germany gib645 I c gib700 I c Silvery Javan gibbon (Hylobates moloch) gib156 V Berlin, Germany Capped gibbon (Hylobates pileatus) gib160 II Capped gibbon (Hylobates pileatus) gib731 II Mulhouse, France Gibbon sp. gib759 d Phuket, Thailand gib760 d gib761 d gib762 gib763 gib764 White-handed gibbon (Hylobates lar) gib174 II (region B not analyzed) Pata, Thailand Capped gibbon (Hylobates pileatus) gib157 II Chiang Mai, Thailand Black gibbon (Hylobates concolor) gib151 gibv Dusit, Thailand gib153 V gib154 V White-checked gibbon (Hylobates leucogenys) gib751 V White-checked gibbon (Hylobates leucogenys) gib824 gibv Cuc Phuong, Vietnam a Sequencing of the whole genome was done with representatives of each genomic group (boldface). b EMBL accession no.: AJ131538 to AJ131557 and AJ131567 to AJ131575. c Member of a family group and part of an infection chain. d Part of an infection chain.

VOL. 74, 2000 NOTES 5379 ses were done with programs from the PHYLIP package, version 3.57c (J. Felsenstein, Department of Genetics, University of Washington, Seattle, 1993). The sequences obtained in our study (Table 2) were compared to published human, chimpanzee, and gibbon sequences. Sequence analysis of the pre-s1/pre-s2 region (region A). As described previously for the classification of human HBV isolates (4), the sequence analysis of region A allowed a rough estimate of the phylogenetic relationships among HBV isolates. A dendrogram of the sequences of region A including human and gibbon HBV isolates is given in Fig. 1. Here, human, chimpanzee, and gibbon HBV isolates clustered in different groups. All isolates obtained from gibbons were related to the gibbon HBV sequence (hbu46935; EMBL accession no. U46935) published by Norder et al. (10). Furthermore, the isolate from the gorilla (gor97a) was found to be distantly related to the chimpanzee HBV isolate (hpbvcg; EMBL accession no. D00220) described by Vaudin et al. (16). Surprisingly, the HBV isolate from the chimpanzee analyzed here (chimp82) was found within a gibbon HBV genomic group (). The robustness of this tree was confirmed by bootstrap resampling with 100 data sets. Here, the gibbon HBV sequences were found within five distinct genomic groups ( to gibv) with intergroup distances of 7 to 11 changes per 100 nucleotides, compared to intergenotype distances of 8 to 22 changes per 100 nucleotides for human HBV genotypes. Sequence analysis of the complete HBV genome from representatives of gibbon HBV genomic groups I to V. The complete genome of at least one representative of each group was sequenced. The results of phylogenetic analysis of complete viral genomes were similar to the ones obtained by analyzing region A (Fig. 2). There is one branch leading to the nonhuman primate HBV isolates, including the isolates from a naturally infected chimpanzee (hpbvcg) (16) and from a gorilla (whose HBV isolate was analyzed in this work). The other branch leads to the human HBV genotypes A to E. A third one leads to human HBV genotype F and to WMHBV (af046996). All gibbon HBV sequences are found within one monophyletic group, and within this gibbon HBV group, five genomic groups can be distinguished. Distances among these groups are be- FIG. 1. Phylogram based on nucleotide sequences of fragment A, encompassing a part of the pre-s1 gene, using the Kimura two-parameter matrix and the neighbor-joining method. The numbers represent the percentages of bootstrap replicates (of 100 total) for each node. Analysis of this region allows a rough estimate of phylogenetic relationships, as all human sequences cluster within their respective genomic groups. A to E, human HBV genotypes; chimp and to gibv, ape HBV genomic groups. FIG. 2. Phylogram based on complete nucleotide sequences, using the Kimura two-parameter matrix and the neighbor-joining method. The numbers represent the percentages of bootstrap replicates (of 100 total) for the nodes. A to F, human HBV genotypes; chimp and to gibv, ape HBV genomic groups.

5380 NOTES J. VIROL. Genotype TABLE 3. Unweighted distances between representatives of all known human HBV genotypes, the chimpanzee and gibbon HBV genomic groups, and WMHBV A B C D Unweighted distance a (exchanges/100 nucleotides) from genotype: E F Chimp b c (n 3) I II V gibv A 3.7 9.4 9.5 9.3 9.9 14.1 9.8 10.5 10.8 10.4 10.2 10.6 22.6 B 4.2 9.1 10.1 10.8 13.7 10.2 10.5 11 10.2 10.2 10 21.8 C 3.3 10.4 10.6 13.9 10.6 10.4 11 10.4 10.3 10.4 22 D 3.9 7.7 13.5 10.1 10.9 11.2 11 11 10.9 22.1 E 1.4 13.4 9 10.2 10.2 10.2 10.2 10.4 22 F 5 12.9 12.8 13.3 13.4 13.1 13.7 23.4 chimp 6 9 9.2 9.4 9.4 9.5 22.3 2.5 6.3 7.5 7.4 7.5 22.3 I NA d 7.8 7.5 7.9 22.5 II 3.8 7.3 7.5 22.8 V NA 6.8 22.2 gibv 4.6 22.4 a Intragroup differences are in boldface. b Includes the isolate from a chimpanzee (16) and the HBV isolate from a naturally infected gorilla. c Includes the gibbon HBV isolate described by Norder et al. (10) and the HBV sequence from a chimpanzee analyzed in this work. d NA, not applicable. The whole genome of only one isolate was analyzed so no intragroup distance can be given. WMHBV tween 6.3 and 7.9 changes per 100 nucleotides, thus slightly lower than distances among genotypes, which are defined by a distance of at least 8% (12). Intragroup distances were between 2.5 and 4.6%, comparable to intragroup distances within human genotypes. A summary of the intra- and intergroup distances of human HBV, chimpanzee HBV, WMHBV, and gibbon HBV variants is given in Table 3. Isolate gib731 (group II) had a 236-nucleotide deletion, from nucleotide 1965 to 2200, within the core gene, which should lead to the formation of an altered, truncated core protein. The HBV sequences obtained from gibbons 759, 760, and 761 were almost identical in both regions A and B. Short distances were due to ambiguities only. Assuming a hepadnaviral mutation rate for nucleotide substitutions of between 1.75 10 5 and 7.62 10 5 substitutions per site per year (13), this sequence identity indicates a recent transmission event among these apes. A second viral transmission was found in a gibbon family group living in a German zoo. The HBV isolates of all three individuals (gibbons 644, 645, and 700) were nearly identical in regions A and B and cluster within genomic group I. Low distances were again caused by ambiguous nucleotide positions. The transmission events are indicated in Table 2. There are some reports about the occurrence of HBV infections in ape populations (6, 17). The question of whether these infections were of human origin or indigenous to the different ape populations is still a matter of debate (7, 17, 18). We found strong evidence for the presumption that these HBV infections are indigenous to the different nonhuman primate populations by analyzing a large number of ape HBV sequences. All isolates analyzed in this work cluster within one monophyletic group, including the previously described chimpanzee and gibbon sequences. Phylogenetic analysis shows that there is one branch leading to human HBV genotypes A to E, a second leading to the ape HBV genomic groups, and a third one leading to human HBV genotype F and to the New World WMHBV. All ape isolates were from animals of different geographic origins and include different species. These facts are consistent with the presumption that there was a single transmission event between humans and apes in the Old World involving a common ancestor of human HBV genotypes A to E and the ape viruses. Prior to this, the branch leading to the New World HBV isolates, WMHBV, and human HBV genotype F might have split off. Norder et al. (11) found indications that group F represented the genomic group of HBV among populations with origins in the New World. With few exceptions, gibbon isolates from the different Thai and the Vietnamese locations could be classified into four genomic groups (, -III, -IV, and -V, with genetic distances between 6.3 and 7.9%), corresponding to the geographical origins of the apes. This finding is analogous to the finding by Norder et al. (12) of different geographical distributions of the different human HBV genotypes. Thus, genotype definition (12) should be discussed for the gibbon HBV variants. Genomic gibbon HBV group I represents nearly identical isolates from a family group living in a German zoo. The HBV isolates from the primates living in other European zoos could be integrated into the Thai genomic groups, II, and V. Therefore, the HBV infection of these apes might be introduced by captured wild apes which were already infected. Surprisingly, the sequence of the HBV isolate of the chimpanzee analyzed here showed a very close phylogenetic relationship to the gibbon isolates within group. Most likely, this reflects an infection of this chimpanzee by a gibbon HBV variant, as this animal was kept with some gibbons in another zoo. Unfortunately, we did not get serum samples from apes of this zoo. However, this is the first probable quasinatural interspecies transmission of HBV. The infectivity of a gibbon HBV variant to a chimpanzee was documented by Mimms et al. (8). In addition, different experimental studies revealed that human HBV isolates can infect gibbons and chimpanzees (1, 2). In general, these results imply that nonhuman primate HBV variants may be transmittable to humans. At least they infect representatives of the families Pongidae and Hylobatidae. However, at present it remains to be elucidated if nonhuman primate HBV variants are transmittable to humans and whether they cause disease in humans. We think that this is likely, because the transmission of human HBV variants to nonhuman primates is well documented. Thus, we strongly recommend routine screening for HBV markers and HBV vaccination of zoo-kept apes, primate-handling zoo keepers, and veterinarians.

VOL. 74, 2000 NOTES 5381 S.G. and J.-O.H. contributed equally to this work. We thank Rainer Thomssen for helpful discussion and personal support. J.-O.H. was supported by the Grimminger Foundation for Zoonosis Research, Stuttgart, Germany. REFERENCES 1. Alter, H. J., R. H. Purcell, J. L. Gerin, W. T. London, P. M. Kaplan, V. J. McAuliffe, J. Wagner, and P. V. Holland. 1977. Transmission of hepatitis B to chimpanzees by hepatitis B surface antigen-positive saliva and semen. Infect. Immun. 16:928 933. 2. Bancroft, W. H., R. Snitbhan, R. M. Scott, M. Tingpalapong, W. T. Watson, P. Tanticharoenyos, J. J. Karwacki, and S. Srimarut. 1977. Transmission of hepatitis B virus to gibbons by exposure to human saliva containing hepatitis B surface antigen. J. Infect. Dis. 135:79 85. 3. Galibert, F., T. N. Chen, and E. Mandart. 1982. Nucleotide sequence of a cloned woodchuck hepatitis virus genome: comparison with the hepatitis B virus sequence. J. Virol. 41:51 65. 4. Grethe, S. 1997. Neue genetische Varianten des Hepatitis B Virus Bedeutung für Pathogenese, Diagnostik, Infektkettenanalyse und Immunprophylaxe, vol. 1. Cuvillier, Göttingen, Germany. 5. Grethe, S., M. Monazahian, I. Bohme, and R. Thomssen. 1998. Characterization of unusual escape variants of hepatitis B virus isolated from a hepatitis B surface antigen-negative subject. J. Virol. 72:7692 7696. 6. Kalter, S. S., R. L. Heberling, A. W. Cooke, J. D. Barry, P. Y. Tian, and W. J. Northam. 1997. Viral infections of nonhuman primates. Lab. Anim. Sci. 47:461 467. 7. Lanford, R. E., D. Chavez, K. M. Brasky, R. B. Burns III, and R. Rico-Hesse. 1998. Isolation of a hepadnavirus from the woolly monkey, a New World primate. Proc. Natl. Acad. Sci. USA 95:5757 5761. 8. Mimms, L. T., L. R. Solomon, J. W. Ebert, and H. Fields. 1993. Unique pres sequence in a gibbon-derived hepatitis B virus variant. Biochem. Biophys. Res. Commun. 195:186 191. 9. Norder, H., A. M. Courouce, and L. O. Magnius. 1994. Complete genomes, phylogenetic relatedness, and structural proteins of six strains of the hepatitis B virus, four of which represent two new genotypes. Virology 198:489 503. 10. Norder, H., J. W. Ebert, H. A. Fields, I. K. Mushahwar, and L. O. Magnius. 1996. Complete sequencing of a gibbon hepatitis B virus genome reveals a unique genotype distantly related to the chimpanzee hepatitis B virus. Virology 218:214 223. 11. Norder, H., B. Hammas, S. D. Lee, K. Bile, A. M. Courouce, I. K. Mushahwar, and L. O. Magnius. 1993. Genetic relatedness of hepatitis B viral strains of diverse geographical origin and natural variations in the primary structure of the surface antigen. J. Gen. Virol. 74:1341 1348. 12. Okamoto, H., F. Tsuda, H. Sakugawa, R. I. Sastrosoewignjo, M. Imai, Y. Miyakawa, and M. Mayumi. 1988. Typing hepatitis B virus by homology in nucleotide sequence: comparison of surface antigen subtypes. J. Gen. Virol. 69:2575 2583. 13. Orito, E., M. Mizokami, Y. Ina, E. N. Moriyama, N. Kameshima, M. Yamamoto, and T. Gojobori. 1989. Host independent evolution and a genetic classification of the hepadnavirus family based on nucleotide sequences. Proc. Natl. Acad. Sci. USA 86:7059 7062. 14. Seeger, C., D. Ganem, and H. E. Varmus. 1984. Nucleotide sequence of an infectious molecularly cloned genome of ground squirrel hepatitis virus. J. Virol. 51:367 375. 15. Testut, P., C. A. Renard, O. Terradillos, L. Vitvitski-Trepo, F. Tekaia, C. Degott, J. Blake, B. Boyer, and M. A. Buendia. 1996. A new hepadnavirus endemic in arctic ground squirrels in Alaska. J. Virol. 70:4210 4219. 16. Vaudin, M., A. J. Wolstenholme, K. N. Tsiquaye, A. J. Zuckerman, and T. J. Harrison. 1988. The complete nucleotide sequence of the genome of a hepatitis B virus isolated from a naturally infected chimpanzee. J. Gen. Virol. 69:1383 1389. 17. Warren, K. S., H. Niphuis, Heriyanto, E. J. Verschoor, R. A. Swan, and J. L. Heeney. 1998. Seroprevalence of specific viral infections in confiscated orangutans (Pongo pygmaeus). J. Med. Primatol. 27:33 37. 18. Warren, K. S., J. L. Heeney, R. A. Swan, Heriyanto, and E. J. Verschoor. 1999. A new group of hepadnaviruses naturally infecting orangutans (Pongo pygmaeus). J. Virol. 73:7860 7865.