key: cord-0699624-efp5s20s authors: Yang, Ling'en; Wu, Jianmin; Hu, Tingsong; Qin, Shaomin; Deng, Bo; Liu, Jinfeng; Zhang, Fuqiang; He, Biao; Tu, Changchun title: Genetic diversity of bat orthohepadnaviruses in China and a proposed new nomenclature date: 2018-09-30 journal: Infection, Genetics and Evolution DOI: 10.1016/j.meegid.2018.05.024 sha: 2c569b07d9a88dd8e8a6b32178fafc89b9edc179 doc_id: 699624 cord_uid: efp5s20s Abstract The orthohepadnaviruses, which include the major human pathogen hepatitis B virus, exist in a wide range of hosts. Since 2013, a large group of orthohepadnaviruses has been identified in bats worldwide and classified as 4 species within the genus Orthohepadnavirus. To further investigate orthohepadnaviruses in the Chinese bat population, 554 archived bat samples from 20 colonies covering 3 southern provinces were screened with results showing that 9 (1.6%) were positive. A systematic phylogenetic analysis has indicated the need for a new nomenclature for bat hepatitis B virus-like viruses: BtHBV, with the addition of 3 new species, one being divided into 6 genotypes. Viruses identified here shared 9.0–19.2% full genome divergence and classified into 3 different genotypes. This study illustrates the genetic diversity of orthohepadnaviruses in the Chinese bat population, and emphasizes need for further investigation of their public health significance. Hepatitis B virus (HBV), the prototype member of the genus Orthohepadnavirus within the family Hepadnaviridae, is one of the most infectious human pathogens, responsible for > 250 million chronic infections worldwide and around 887,000 human deaths annually (ICTV, 2017; Mason et al., 2011; Ott et al., 2012; Seeger et al., 2013; WHO, 2017) . Orthohepadnaviruses have also been characterized in other primates and rodents (Seeger et al., 2013) . Additionally, 14 species of bats from the Old World (Myanmar, China and Gabon) and New World (Panama) have recently been found to harbor a variety of orthohepadnaviruses (Drexler et al., 2013; He et al., 2013 He et al., , 2015 Nie et al., 2018; Wang et al., 2017) , of which 4 have been approved in the latest report of the International Committee on Taxonomy of Viruses (ICTV) as new species within the genus Orthohepadnavirus: Long-fingered bat hepatitis B virus, Pomona bat hepatitis B virus, Roundleaf bat hepatitis B virus, and Tent-making bat hepatitis B virus (ICTV, 2017) . Molecular studies have shown that the Panamanian tent-making bat hepatitis B virus (TBHBV) is capable of infecting primary human hepatocytes through binding to the same receptor as human HBV, thereby suggesting that bat orthohepadnaviruses are ancestors of primate HBVs (Drexler et al., 2013; Rasche et al., 2016) . Orthohepadnaviruses are much more diverse in bats than in other host species. A single bat species can harbor diverse strains of these viruses, and a single orthohepadnavirus species can infect different bat species: e.g., Miniopterus schreibersi (M. schreibersi) bats have been found carrying both Neixiang-Ms69 and Anlong-Ms258 viruses, two viruses showing enough genetic diversity to be considered as two distinct species (Nie et al., 2018) , whereas Neixiang-Ms69 and Neixiang-Rpu92, sharing almost 100% full genome nucleotide (nt) similarity, have been found in two distinct bat species of two separate families, M. schreibersi and Rhinolophus pusillus (R. pusillus) (Nie et al., 2018) . Currently nomenclature of a bat orthohepadnavirus species uses HBV after the common name of the host in which they have been first detected; e.g., "Long-fingered bat hepatitis B virus" representing a virus species identified in Miniopterus bats, but this could refer to either strain 776 from Myanmar or Neixiang-Ms-69 from Chinatwo different viruses (He et al., 2013; Nie et al., 2018) . Such a mismatch can cause confusion. Accordingly, a new species nomenclature for bat HBV-like viruses (BtHBV) is necessary. With the exception of humans and possibly rodents, bats are the most abundant as well as the most widely distributed mammals (Wilson and Reeder, 2005) . They provide a huge natural virus bank from which > 200 viruses have been isolated or identified, including such lethal agents as lyssaviruses, henipaviruses, ebolaviruses, and SARSrelated coronaviruses (Moratelli and Calisher, 2015) . Apart from orthohepadnaviruses in bats in America, Asia and Africa, homologues of hepatitis A, C and E virus have also been sporadically identified in these animals (Drexler et al., 2015; Drexler et al., 2012; Quan et al., 2013) . To further understand the genetic diversity of HBV-like virus, we screened hepatitis viruses from bats collected in south China, and found another three variants. Systematic and phylogenetic analyses of the BtHBVs identified 7 species, for which we propose a new nomenclature that will be more helpful for understanding host species specificity and evolution of these viruses. The procedures for sampling of bats in this study were reviewed and approved by the Administrative Committee on Animal Welfare of the Institute of Military Veterinary, Academy of Military Medical Sciences, China (Laboratory Animal Care and Use Committee Authorization, permit numbers: JSY-DW-2010-02 and JSY-DW-2015-01). All live bats were maintained and handled according to the Principles and Guidelines for Laboratory Animal Medicine (2006), Ministry of Science and Technology, China. To investigate orthohepadnaviral ecology in Chinese bats, animals were collected between 2005 and 2015 with nets near or in human-inhabited communities in Yunnan, Guangxi and Guangdong provinces. All bats were apparently healthy at capture, and were initially identified by morphological examination by a trained field biologist and confirmed by sequence analysis of the mt-cyt b gene (Wang et al., 2003) . Following euthanization their livers and other organs were immediately collected and stored at −80°C ( Fig. 1 and Table 1 ). To screen for orthohepadnaviruses, viral DNA was extracted from liver tissue of each of the 554 bats using the QIAamp DNA Mini Kit (Qiagen) in a QIAcube (Qiagen).Virus detection was conducted using the TaKaRa PCR Kit (TaKaRa) as previously described (He et al., 2013) . Double-distilled water was used as a negative control. Two samples of each lineage were chosen for amplification of their complete genomes using LA Taq polymerase (TaKaRa) and previously described primers (He et al., 2013) . Positive PCR amplicons were ligated into pMD18T vector (TaKaRa) and used to transfect competent Escherichia coli DH5α cells (Tiangen). > 3 clones of each amplicon were randomly picked for sequencing by the Sanger method in an ABI 3730 sequencer (Coma-teBio). Overlapping amplicons were assembled with SeqMan v7.0 into full genomic sequences. Genomic structures were determined by Vector NTI v.8, followed by comparison with those of other orthohepadnaviruses. Full genomes and predicted open-reading frames (ORF) Pol, preS1/S2/S, preC/C, and X (only for orthohepadnaviruses) of representatives of hepadnaviruses, including unclassified ones, retrieved from Genbank were aligned with counterparts of viruses identified in this study using MAFFT v7.394 (Katoh and Standley, 2013) , with nt identities being calculated using MegAlign v7.1 (DNASTAR, Inc., Madison, WI). Evolutionary models were determined using ModelGenerator v0.85 with Akaike Information Criterion 1(AIC1) (Keane et al., 2006) . Phylogenetic and molecular evolutionary analyses were achieved using PhyML v3.3 by the maximum likelihood and best-fit substitution models under evaluation of 100 bootstraps (Guindon et al., 2010) . The generated tree files were visualized with FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/ figtree). Full genomic sequences obtained here were deposited in Genbank under the accession numbers KY905324-KY905329. During 2005-2015, 554 frugivorous and insectivorous bats representing 13 species were collected from 1 to 3 adjacent colonies in 11 locations in Yunnan, Guangxi and Guangdong provinces ( Fig. 1 ). Collection details are shown in Table 1 . Pan-orthohepadnavirus PCR screening gave positive results in 9 bat samples from 3 colonies. Of these, 10.0% (4/40) and 9.4% (3/32) were from Hipposideros pomona (H. pomona) bats in Zhenyuan, Yunnan province, 2013, and Baise, Guangxi province, 2014 respectively. The remainder (2/3) came from H. armiger bats in Hechi, Guangxi province, 2015 (Table 1 ). The amplicons were sequenced, and phylogenetic analysis revealed that they formed three distinct lineages. Of note is that the three lineages were from three different bat colonies. All BtHBVs in this study were named as described previously (He et al., 2014) , the first two letters representing the sampling location, with the remaining letters identifying the host species and numbers referring to the sampling order. Full genomes of six samples with two of each lineage were successfully sequenced. Analyses showed that the 6 full genomes ranged from 3254 to 3275 nt in length, consistent with previously published BtHBVs (3230-3287 nt) in Asia, including PEPRs in Pu′er and Rs3364 in Jinning, Yunnan, China, and 776 in Kachin state, Myanmar (He et al., 2013) . As in all orthohepadnaviruses, the viruses identified here contained four ORFs: polymerase (Pol), surface (S), core (C), and X. Positions of all ORFs in these bat orthohepadnaviruses were similar to ORFs of currently reported members of the Orthohepadnavirus genus but were clearly distinct from the ORFs of Avihepadnavirus. The S protein genes encoded in the large ORF contained preS1, preS2, and S domains. The preS1 domain contained a signal of N-myristoylation at glycine-2. Additional to their ORF organization, human HBV and the BtHBVs shared a similar location for the direct repeat (DR) sequences DR1 and DR2 involved in genome replication. Secondary structure prediction highlighted the structural similarities between BtHBVs and HBVs of other animal origin in their ε-loops, which serve as templates for the priming of reverse transcription of pre-genomic RNA in all hepadnaviruses (Seeger et al., 2013) . All representatives of hepadnaviruses available in Genbank were retrieved and aligned with the sequences obtained here. The best-fit substitution models, under AIC1, for full genomes, ORFs Pol, preS1/S2/ S, preC1/C and X were selected as GTR + G [Gamma distribution parameter alpha (a) =0.61], TVM + G (a = 0.61), GTR + G (a = 0.62), GTR + I + G (I = 0.05, a = 0.63), and TVM + G (a = 0.45), respectively. Phylogenetic trees based on the full genomic and ORF nucleotide sequences are shown in Fig. 2-4 , with the (Fig. 2, 5) . The larger group can be divided into 6 clades (BtHBVs 2-7) according to divergence limits used by ICTV as species demarcation criteria (Mason et al., 2011) (Fig. 2, 5) . BtHBV 2 consists of a lineage of viruses from Myanmar M. fuliginosus bats (He et al., 2013) (Fig. 2) . BtHBV 3 and 4 comprise viruses with a divergence of 21.1% isolated from Gabonese H. cf. ruber and R. alcyone bats, respectively (Drexler et al., 2013) (Fig. 2, 5) . BtHBV 5 consists of one isolate, Jiyuan-Rf-11 from Chinese R. ferrumequinum bats, showing 22.7-23.4% divergence with Gabonese viruses (Drexler et al., 2013; Wang et al., 2017) (Fig. 2, 5 ). BtHBV 6 is composed of viruses from a single location in China, but from two bat species, R. pearsonii and R. pusillus (Nie et al., 2018) (Fig. 2) . BtHBV 7 has diverse members from a variety of bat species in widely spread geographic locations in China, such as the PEPRs from H. pomona in Yunnan, Longquan-Ha-273 from H. armiger in Zhejiang, and Neixiang-Rm-103 from R. monoceros in Henan (He et al., 2015; Nie et al., 2018) (Fig. 2) .This clade can be further divided into 6 lineages (A to F) according to 6% divergence that has been used to differentiate genotypes of human HBVs (Norder et al., 1994) (Fig. 2, 5) . The six full genomes of bat viruses described here clustered within three different lineages (B to D) with < 06% intra-and 9.0-19.1% inter-lineages divergence within the clade BtHBV 7. The full genomes of the BSPRs (lineage B) from Baise, Guangxi, were much closer (~9.0% divergence) to the ZYPRs (lineage C) from Zhenyuan, Yunnan, than to the HCGRs (lineage D, 17.5% divergence) from Hechi, Guangxi, which showed~19% divergence from the ZYPRs (Fig. 5) (Fig. 5) , followed by 25.3-52.9% with members of BtHBV 1-6 ( Fig. 5) , then > 44.6% with rodent and primate orthohepadnaviruses (Fig. 5) . The HCGRs formed lineage D with Longquan-Ha-273, with the least divergence of 5% between each other (Fig. 2, 5) , and also shared the same host species, H. armiger, but in two locations > 1300 km apart (Nie et al., 2018) (Fig. 1) , then 16.8-19.8% with other lineages within the BtHBV7 (Fig. 2, 5) . The ORFs of Pol, preS1/S2/S, preC/C, and X were also used to validate the phylogenetic grouping. Generally, the Pol, preS1/S2/S, preC/C genes showed very similar topologies to the full genomes, and all BtHBVs can clearly be divided into 7 clades with BtHBV 7 containing 6 lineages ( Fig. 2, 3, 4A ). Due to the lack of X gene in other hepadnaviruses, the topology of orthohepadnaviruses differed from other trees, with BtHBV1s placing as an outgroup from other orthohepadnaviruses, and with the remaining BtHBVs forming 6 clades (Fig. 4B) . Unlike other trees within which BtHBV 4 and 5 shared an evolutionary root, BtHBV 5 in the X gene tree formed an outgroup to BtHBVs 3, 4, 6 and 7, indicating a different evolutionary rate of X gene (Fig. 4B) . Until very recently, the family Hepadnaviridae contains two genera, Avihepadnavirus and Orthohepadnavirus, the former comprising three species and the latter having 8 members (ICTV, 2017) . The species demarcation criteria in the genus Orthohepadnavirus are somewhat indefinite: for example, the full genome sequence divergence of WHV/ HBV is 40%, GHSV/WHV 15%, WMHBV/HBV 20%, and WMHBV/WHV 30%, but the genetic divergences between different species have tentatively been set at ≥20% (ICTV, 2017; Mason et al., 2011) . Recently a virus found in the white sucker fish (White sucker hepatitis B virus, WSHBV) has also been approved as a new species (White sucker hepatitis B virus), but unassigned to any genus (Hahn et al., 2015) . Additionally, two other hepadnaviruses have been identified in fish (Bluegill hepadnavirus, BGHBV) and amphibians (Tibetan frog hepadnavirus, TFHBV), with both showing a high genetic diversity compared with other prototypes, and might represent totally new genus candidates (Dill et al., 2016) . Currently viral nomenclature within the family follows the conventional tradition: hepatitis B virus (for avian, primate and chiropteran viruses) or hepatitis virus (for rodent viruses) prefixed by the host English name, e.g., duck hepatitis B virus (DHBV), woolly monkey hepatitis B virus (WMHBV), pomona bat hepatitis B virus (PBHBV), and woodchuck hepatitis virus (WHV). With increasing numbers of BtHBVs being reported, the current nomenclature cannot sustain the expansion of diversity of BtHBVs. Long-fingered bats and roundleaf bats are members of two genera, Miniopterus and Hipposideros, with each having > 10 species (Wilson and Reeder, 2005) . Within the Miniopterus spp. bats, not only do M. fuliginosus harbor orthohepadnaviruses (e.g. 776) (He et al., 2013) , but M. schreibersi from China also carry such a virus -Neixiang-Ms-69 (Nie et al., 2018) . However, Neixiang-Ms-69 only shows a distant relationship with isolate 776, with 35.3% divergence (Fig. 5) ; i.e., with a sufficient difference to separate them into two species. Similarly, GB09-303, the prototype of RBHBV, shares the same host genus (Hipposideros spp, roundleaf bats) as the HCGRs, but with~32% divergence -i.e., much lower than the species demarcation criterion (20%)indicating that they are members of two different species. The same situation has also been observed in Rhinolophus spp. bats (horseshoe bats) (Fig. 2) , which harbor diverse orthohepadnaviruses showing 22.7%-31.7% divergences that have been grouped into 4 clades, BtHBV 4-7. Accordingly, unlike non-bat HBVs that can use the nomenclature of HBV or HV specified by host common English name to clearly distinguish the virus species, use of names such as LBHBV, RBHBV or HBHBV cannot precisely provide a match to a single viral species, which results in confusion in our understanding of the evolution and taxonomy of BtHBVs. Here we propose and have useda new nomenclature: BtHBV for bat-associated hepatitis B virus-like viruses, with different species marked by number, and with the strain name containing the species abbreviation, identifier, host species abbreviation, country, and collection time, so that all essential information about a virus can be clearly read from its name. Accordingly, presently named species have been changed to BtHBVs 1-3 and 7 (Table 2, Fig. 2 ). With increasing numbers of orthohepadnaviruses being discovered in Chinese bats, 3 additional clades can be assigned showing > 20% divergence: namely BtHBV 4-6, containing new species candidates. The 3 lineages of viruses identified here all fell into BtHBV 7 as variants. Hepadnaviruses have long been limited to primates, rodents and birds, but diverse hepadnaviruses have recently been found in bats from both Old and New Worlds. With orthohepadnaviruses increasingly being reported in bats in China The increasing numbers of Chinese BtHBVs being reported (He et al., 2015; Nie et al., 2018; Wang et al., 2017) , are highly divergent genetically and distributed within all 3 of the proposed species, BtHBVs 5-7 (Fig. 2-5) . Such divergence is well reflected in BtHBV 7. To better understand the evolutionary dynamics of human HBVs, HBVs have been classified into genotypes A-J according to a divergence cutoff of 6% (Miyakawa and Mizokami, 2003; Olinger et al., 2008; Tatematsu et al., 2009) . Similarly, BtHBV 7 can also be divided into 6 lineages based on the same cutoff. Accordingly here we propose 6 genotypes within BtHBV 7, and viruses found here constitute genotypes B, C and D (Fig. 2) . All these findings indicate that orthohepadnaviruses in Chinese bats are much more diverse than has been expected and, undoubtedly, more new viruses or variants will be discovered in these creatures, which will be helpful in understanding Fig. 5 . Full genomic identity comparison of BtHBVs, and human-and rodent HBVs. Strains that showed little divergence (< 1%) with phenotypes of 4 bat HBV species (triangles) are excluded. Dark gray: species with divergences within the ≥20% cutoff. Light gray: species with divergences of between 6% and 20%. Circles: viruses identified in this study. ***: not available. their origin and evolution. HBVs are transmissible not only among humans but can also infect certain other primates, such as chimpanzees and macaques. Laboratory infection of the tree shrew (Tupaia belangeri) has also been reported (Seeger and Mason, 2015) , which suggests possible inter-species infection and transmission. It appears that this is also the case for BtHBVs, since some viruses within a species have been detected in distinctly different host species. For example, members of BtHBV 7 have been detected in bats of 10 species within 4 families (Fig. 2) , and even within BtHBV 6 there are 3 members with almost no difference but isolated from R. perarsonii and R. pusillus of the genus Rhinopophus, family Rhinolophidae, and M. schreibersi, genus Miniopterus, family Miniopteridae (Fig. 2) , which indicates that BtHBVs can jump the host barrier to circulate between different bat species. However, the HCGRs of H. armiger bats in Hechi do not appear to have infected cave-sharing H. pomona bats that are hosts to a number of BtHBVs, such as the PEPRs, BSPRs, and ZYPRs, even though both bat species are members of the genus Hipposideros. This may reflect potential host/genetic species barriers, although we cannot exclude the possibility that ecological reasons may prevent this cross-species transmission and/or that the limited sample size and survey duration has prevented the identification of positive individuals. All these observations suggested that the circulation of orthohepadnaviruses in bat population is very intricate. Understanding its complexity will require a study of viruses rescued from different infected hosts and epidemiological investigations of orthohepadnaviruses over a broader range of bat species, wider areas and of longer duration. Distinct viral lingeages from fish and amphilians reveal the complex evolutionary history of hepadnaviruses Bats worldwide carry hepatitis E virus-related viruses that form a putative novel genus within the family Hepeviridae Bats carry pathogenic hepadnaviruses antigenically related to hepatitis B virus and capable of infecting human hepatocytes The Hepatovirus Ecology Consortium, 2015. Evolutionary origins of hepatitis A virus in small mammals New algorithms and methods to estimates maximum-likelihood phylogenies: assessing the performance of PhyML 3.0 Characterization of a novel hepadnavirus in the white sucker (Catostomus commersonii) from the great lakes region of the United States Hepatitis virus in long-fingered bats Identification of diverse alphacoronaviruses and genomic characterization of a novel severe acute respiratory syndrome-like coronavirus from bats in China Identification of a novel Orthohepadnavirus in pomona roundleaf bats in China The 10th report of International Committee on Taxonomy of Viruses MAFFT multiple sequence alignment software version 7: improvements in performance and usability Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified Family Hepadnaviridae Classifying hepatitis B virus genotypes Bats and zoonotic viruses: can we confidentily link bats with emerging deadly viruses? Extensive diversity and evolution of hepadnaviruses in bats in China Complete genomes, phylogenetic relatedness, and structural proteins of six strains of the hepatitis B virus, four of which represent two new genotypes Possible new hepatitis B virus genotype, southeast Global epidemiology of hepatitis B virus infection: new estimates of age-specific HBsAg seroprevalence and endemicity Bats are a major natural reservoir for hepaciviruses and pegiviruses Bat hepadnaviruses and the origins of primate hepatitis B viruses Molecular biology of hepatitis B virus infection. Virology 479-480 Hepadnaviruses. In: Fields Virology A genetic variant of hepatitis B virus divergent from known human and ape genotypes isolated from a Japanese patient and provisionally assigned to new genotype Molecular phylogenetic of Hipposiderids (Chiroptera: Hipposideridae) and Rhinolophids (Chiroptera: Rhinolophidae) in China based on mitochondrial cytochrome b sequences Detection and genome characterization of four novel bat hepadnaviruses and a hepevirus in China The study was supported by the General Program of the National Natural Science Foundation of China (31572529), the National Key Basic Research and Development Program of China (2016YFC1200100, 2017YFD0501800).The funders played no role in study design, data collection and interpretation, or the decision to submit the work for publication.