key: cord-287748-co9j3uig
authors: Kobayashi, Tomoya; Murakami, Shin; Yamamoto, Terumasa; Mineshita, Ko; Sakuyama, Muneki; Sasaki, Reiko; Maeda, Ken; Horimoto, Taisuke
title: Detection of bat hepatitis E virus RNA in microbats in Japan
date: 2018-05-29
journal: Virus Genes
DOI: 10.1007/s11262-018-1577-9
sha: 
doc_id: 287748
cord_uid: co9j3uig

Several recent studies have reported that various bat species harbor bat hepatitis E viruses (BatHEV) belonging to the family Hepeviridae, which also contains human hepatitis E virus (HEV). The distribution and ecology of BatHEV are not well known. Here, we collected and screened 81 bat fecal samples from nine bat species in Japan to detect BatHEV RNA by RT-PCR using HEV-specific primers, and detected three positive samples. Sequence and phylogenetic analyses indicated that these three viruses were BatHEVs belonging to genus Orthohepevirus D like other BatHEV strains reported earlier in various countries. These data support the first detection of BatHEVs in Japanese microbats, indicating their wide geographical distribution among multiple bat species.

Bats are known to be natural reservoirs of various zoonotic viruses such as rabies virus, Nipah virus, and severe acute respiratory syndrome (SARS) coronavirus [1] [2] [3] . In addition, hepatitis E virus (HEV)-like viruses have been detected in bats in several countries [4, 5] . Bat HEV (BatHEV) belongs to family Hepeviridae, genus Orthohepevirus and is a non-enveloped, positive-sense, single-stranded RNA virus. Orthohepevirus is divided into four species, Orthohepevirus A-D [6, 7] . Orthohepevirus A contains HEV, which is a causative agent for human acute hepatitis. Orthohepevirus B includes avian-HEV associated with hepatitis-splenomegaly syndrome. Orthohepevirus C is divided into C1 (rodent-HEV) and C2 (carnivore-associated HEV) genotypes. Orthohepevirus D comprises BatHEV, which has 57.4-64.8% identity with HEV [8] . Recently, BatHEVs were detected from macrobats and microbats in a variety of countries [4, 5] . However, limited information about the distribution of BatHEVs in other regions and their ecology is available. Here, we surveyed several Japanese microbat species to detect BatHEVs.

To examine whether BatHEVs exist in Japanese bats, we collected 81 bat fecal samples from nine bat species captured in five different prefectures of Japan in 2015, with permission from the Ministry of the Environment, Japan and the respective local government (Fig. 1a ). Bats were caught using a harp trap and kept in a pouch for an hour to check the signs of disease as well as to obtain fresh feces. All captured bats did not show any obvious symptoms. The feces were added in a medium containing antibiotics, and frozen in dry ice. We extracted RNA from the fecal samples and performed RT-PCR using a primer set ( Table 1) that was designed specifically against the conserved region of the RNA-dependent RNA polymerase (RdRp) of HEV, to screen the HEV genomes. Two samples (BtHEV-Ej1 and BtHEV-Ej2) from the Japanese short-tailed bat (Eptesicus japonensis) and 1 sample (BtHEV-Ps1) from the Brown long-eared bat (Plecotus sacrimontis) were found to be positive. Although we attempted to sequence the entire genome of the BatHEVs, we failed to amplify the whole genome using RT-PCR. Therefore, we amplified partial ORF1, which encodes non-structural protein including RdRp, and entire capsid coding region of BtHEV-Ej1, -Ej2, and -Ps1 (corresponding to nucleotides (nt) 2171-6690, nt 3585-6690, and nt 3847-6690 of BatHEV/BS7, respectively) using specific primers ( Table 1 ). The sequences of BtHEV-Ej1 and BtHEV-Ej2, which were collected in the same site on different days, were 99% identical (nt 2904 out of 2915). The identity between BtHEV-Ej1/-Ej2 and BtHEV-Ps1 using a part of the RdRp was about 75% at the nucleotide level, suggesting the presence of multiple BatHEVs in Japan. BLAST analysis indicated that BtHEV-Ej1/-Ej2 showed the highest sequence identities to BatHEV/BS7, a German strain detected from the Serotine bat (Eptesicus serotinus), among strains previously reported in other countries. On the other hand, BtHEV-Ps1 showed the highest sequence identity to BtHEV/NMS098B, a German strain detected from the Daubenton's bat (Myotis daubentonii). In particular, 82% identity observed between BtHEV-Ej1/-Ej2 and BatHEV/BS7 or 77% identity observed between BtHEV-Ps1 and BtHEV/NMS098B were greater than that observed between Japanese strains. These data suggest that similar viruses exist in geographically distant regions.

We then phylogenetically analyzed the sequences by maximum-likelihood analysis using ClustalW and MEGA version 7.0 [9] . A phylogenetic tree constructed using the partial amino acid sequences of RdRp indicated that the Japanese viruses were included in Orthohepevirus D (Fig. 1b) , demonstrating that all BatHEVs are classified in this species. We also amplified the full-length capsid (ORF2) sequences by RT-PCR and analyzed them phylogenetically. The resulting tree confirmed that the novel Japanese viruses were included in Orthohepevirus D (Fig. 1c) . The sequences of BtHEV-Ej1/-Ej2 were placed in a position neighboring the German BatHEV/BS7 strain, confirming the phylogenetic similarity between these strains.

All bats captured in this study were insectivores and hibernate in winter. Although there is no information about migration of E. japonensis and P. sacrimontis, bat species closely related to them were reported to migrate only a few kilometers from their colonies per night [10, 11] . They had different ecology in the terms of the habitat. BtHEV-Ej1 and -Ej2 were detected on different days from two E. japonensis bats, both of which used eaves of the same house as night roost in Nagano. Although the E. japonensis bats formed mixed colony at the roost with Myotis ikkonikovi and Rhinolophus ferrumequinum bats, BatHEVs were only detected in E. japonensis, implying BtHEV-Ej1/-Ej2 might have a narrow host range. P. sacrimontis bats usually form a small colony without other species of bats. Indeed, P. sacrimontis bats, from which BtHEV-Ps1 was detected in this study, were captured near such small colony in a ruin in Aomori. Thus, BtHEV-Ps1 is likely circulating in the P. sacrimontis since the bats have low opportunity to come in contact with other species of bats.

The closely related BatHEVs (BtHEV-Ej1/-Ej2 and Bat HEV/BS7) have been detected in different species of Eptesicus bats (E. serotinus and E. japonensis). Since the distribution areas of these bats are not overlapping, viruses ancestral to BtHEV-Ej1/-Ej2 and BatHEV/BS7 might have infected ancestral Eptesicus and might have branched into different species in the process of evolution.

For virus isolation, we inoculated the RT-PCR-positive fecal samples not only into several bat cells (BKT, FBKT, and DemKT1 cells) but also into other mammalian cell lines (Madin-Darby canine kidney (MDCK), African green monkey VeroE6, human A549, Madin-Darby bovine kidney (MDBK), and swine PK15 cells) since we suspected that the bat fecal samples may contain several pathogens other than BatHEVs. All inoculated cells were incubated for 12-15 days with media changes at 2-3 days interval. After the incubation, cells were blindly passaged three times. However, we could neither recover any infectious viruses nor detect BatHEV RNA in the inoculated cells by RT-PCR.

In conclusion, the present study showed the presence of several BatHEV strains, which were independently classified To obtain fresh feces, bats were kept in a pouch for an hour, and the feces were then collected by a sterilized cotton bud and transferred to 1 mL of Dulbecco's modified medium Eagle's minimum essential medium (DMEM) supplemented with 100 U/mL of penicillin, 1 mg/mL of streptomycin, 100 µg/mL of gentamycin, and 2 µg/mL of amphotericin. The feces were suspended well and then centrifuged at 10,000×g for 15 min at 4 °C. The supernatants were used for RNA extraction with ISOGEN LS reagent (Nippon Gene). cDNA was synthesized using Prime Script RT reagent kit (Takara Bio) with a mixture of random hexamer and oligo dT primers. PCR amplifications were performed using the KOD FX Neo (Toyobo) with consensus HEV primer sets (PanHEV F and R), which were designed in this study to amplify a 191-bp fragment of the RNA-dependent into Orthohepevirus D, in Japanese bats, suggesting wide geographical distribution of BatHEV among multiple bat species. Although these data suggest limited transmissibility of BatHEV to other animals, further studies are needed to determine its zoonotic potential. 

Fields Virology

members of the International Committee on the Taxonomy of Viruses Hepeviridae Study Group

Acknowledgements We thank Mr. Mitsuru Mukohyama for helping us