key: cord-0778210-leci3y54
authors: Kemenesi, Gábor; Gellért, Ákos; Dallos, Bianka; Görföl, Tamás; Boldogh, Sándor; Estók, Péter; Marton, Szilvia; Oldal, Miklós; Martella, Vito; Bányai, Krisztián; Jakab, Ferenc
title: Sequencing and molecular modeling identifies candidate members of Caliciviridae family in bats
date: 2016-07-31
journal: Infection, Genetics and Evolution
DOI: 10.1016/j.meegid.2016.04.004
sha: c9e0215e61b3b268bcba56c10f8363360c1f4d24
doc_id: 778210
cord_uid: leci3y54

Abstract Emerging viral diseases represent an ongoing challenge for globalized world and bats constitute an immense, partially explored, reservoir of potentially zoonotic viruses. Caliciviruses are important human and animal pathogens and, as observed for human noroviruses, they may impact on human health on a global scale. By screening fecal samples of bats in Hungary, calicivirus RNA was identified in the samples of Myotis daubentonii and Eptesicus serotinus bats. In order to characterize more in detail the bat caliciviruses, large portions of the genome sequence of the viruses were determined. Phylogenetic analyses and molecular modeling identified firmly the two viruses as candidate members within the Caliciviridae family, with one calicivirus strain resembling members of the Sapovirus genus and the other bat calicivirus being more related to porcine caliciviruses of the proposed genus Valovirus. This data serves the effort for detecting reservoir hosts for potential emerging viruses and recognize important evolutionary relationships.

Despite ongoing advances in biomedicine and molecular biology, infectious diseases remain a major threat to human health, economic sustainability and wildlife conservation. In the last decades bats became a possible source for several human infecting agents such as Ebola, SARS and MERS coronaviruses, Nipah virus, and Hendra virus (Wynne and Wang, 2013) . With over 1250 species distributed globally, bats represent the most diverse and species rich taxa among mammals (Teeling et al., 2005) . The zoonotic potential of numerous bat-harbored viruses has been described also in Europe (Kohl and Kurth, 2014) . Additionally, several novel viruses have been detected recently; however, there is no data on their role in association with bat or human diseases (Kemenesi et al., 2014 (Kemenesi et al., , 2015a .

The family Caliciviridae comprises small non-enveloped RNA viruses, which are classified in five genera: Vesivirus, Lagovirus, Norovirus, Sapovirus and Nebovirus. In addition, unclassified caliciviruses have been identified in monkeys, pigs, avian species and fishes (Farkas et al., 2008; L'Homme et al., 2009; Wolf et al., 2011; Liao et al., 2014; Mikalsen et al., 2014) . Calicivirus related sequences have been also reported from different bat species in Hong Kong, Hungary and New Zealand (Tse et al., 2012; Kemenesi et al., 2014; Wang et al., 2015) . The common characteristic of these bat-derived calicivirus sequences is the phylogenetic divergence from other known members within the family Caliciviridae. The generation of partial or full-length genome sequence data on these recently discovered viruses is crucial to improve viral taxonomy and classification, to develop reliable diagnostic assays and to better understand the molecular mechanisms driving the evolution and host species adaptation of caliciviruses. Human and animal caliciviruses may evolve by accumulation of punctate mutations and by exchange of genetic material via recombination (Giammanco et al., 2013) . Also, interspecies transmission has been documented repeatedly and the origin of some human calicivirus strains may have been triggered by a combination of various evolutionary mechanisms (Smith et al., 1998; Scheuer et al., 2013) . In this regard the possible reservoir role of bats has been investigated minimally. Here we report the Infection, Genetics and Evolution 41 (2016) 227-232 genomic characterization and molecular modeling based analysis for two novel calicivirus strains (BtCalV/BS58/HUN/2013 and BtCalV/ M63/HUN/2013) identified from Serotine bat (Eptesicus serotinus) and Daubenton's bat (Myotis daubentonii).

Bat guano samples were collected during the swarming period in 2013 as part of a countrywide survey program as described previously (Kemenesi et al., 2014) . Nucleic acid extractions were performed with DiaExtract Total RNA Isolation Kit (Diagon Ltd., Hungary). 3′ RACE PCRs were fulfilled on both samples with SuperScript III One-Step RT-PCR System with Platinum Taq DNA Polymerase (Invitrogen), according to the manufacturer's instructions with using the previously described forward primer by Jiang et al. (1999) . 5′ RACE PCR was also performed with 5′/3′ RACE Kit, 2nd Generation (Roche). In order to obtain sequences upstream the RNA dependent RNA polymerase (RdRp) region we performed a two-step RT-PCR protocol. cDNA synthesis was carried out with the Transcriptor First Strand cDNA Synthesis Kit (Roche) using random hexamer oligonucleotide mixture. PCR assays were implemented with Thermo Scientific Phusion High-Fidelity PCR Master Mix following the protocol provided by the manufacturer, using previously published generic calicivirus forward primer (L'Homme et al., 2009) with RdRp specific reverse primer (Jiang et al., 1999) .

High throughput sequencing was performed on PCR amplicons to obtain primary sequence data (Mihalov-Kovács et al., 2015) . Briefly, an Ion Torrent compatible library was prepared applying the NEBNext® Fast DNA Fragmentation & Library Prep Set for Ion Torrent (New England Biolabs) with the Ion Torrent Xpress barcode adapters (Life Technologies). The emulsion PCR to obtain clonally amplified fragments was carried out using the Ion OneTouch™ 200 Template Kit (Life Technologies) on a OneTouch version 2 equipment (Life Technologies) as recommended by the manufacturer. Templated beads were enriched using an Ion OneTouch™ ES pipetting robot (Life Technologies). The 200 bp sequencing protocol was performed on a 316 chip (Life Technologies) using the Ion Torrent PGM (Life Technologies) semiconductor sequencing equipment. Trimmed sequence reads were used for de novo assembly utilizing the MIRA (version 3.9.17) (Chevreux et al., 1999) . Additional bioinformatic analyses, validation with remapping were performed using the CLC Genomics Workbench (version 6.5.1; http://www.clcbio.com) and the DNAStar (version 12; http://www. dnastar.com). GenBank accession numbers for BtCalV/BS58/HUN/2013 and BtCalV/M63/HUN/2013 are KJ652318 and KJ652319, respectively.

Basic sequence alignments and sequence manipulations were fulfilled with ClustalX v2.1 and GeneDoc v2.7 softwares. Multiple evolutionary models were tested on our sequence alignments with PhyML v3.0 software in order to designate the best phylogenetic model for our data set. The phylogenetic trees were constructed with MEGA v5.0 software using the Maximum-Likelihood method, based on the General Time Reversible model with Gamma Distribution (GTR + G). Number of bootstraps for simulations was 1000.

The 3D structure models of mature BtCalV/M63/HUN/2013 (aa. 317 to 831) and BtCalV/BS58/HUN/2013 (aa. 268 to 746) bat calicivirus VP1 capsid protein (CP) structures were generated with I-TASSER (Zhang, 2008; Roy et al., 2010) . The models were built using the recently identified bat VP1 CP amino acid sequences. The following templates were used to thread the VP1 CP structure of the bat calicivirus BtCalV/M63/ HUN/2013: PDB ID codes 3M8L (Feline calicivirus (FCV) capsid protein) and 2GH8 (X-ray structure of a native calicivirus). The structure of Norwalk virus capsid (PDB ID code 1IHM) was used for modeling of the VP1 capsid protein of the bat calicivirus BtCalV/BS58/HUN/2013. The raw protein model structures were refined with the MacroModel energy minimization module of the Schrödinger Suite (Schrödinger, 2015) to eliminate the steric conflicts between the side-chain atoms. Pairwise protein sequence alignments were calculated with the NeedleP tool of the SRS bioinformatics software package. The T = 3 virion models were created with the Oligomer Generator application of VIPERdb (available at http://viperdb.scripps.edu/oligomer_multi.php). Prior to virion model generations the asymmetric units were constructed with Schrödinger Suite using the coordinates of subunit A, B and C of FCV and Norwalk virus in vdb convention format. Molecular graphics and sequence alignment visualization were prepared using VMD v1.9.1 (Humphrey et al., 1996) and the Multiple Sequence Viewer of the Schrödinger Suite, respectively.

Since 5′ RACE PCR protocols were not successful for the two samples, we obtained a 3278 and 3130 nt long partial genome sequences, respectively, for the strains BtCalV/M63/HUN/2013 and BtCalV/BS58/HUN/ 2013. The sequenced genomic fragment spanned the partial RNAdependent RNA polymerase (RdRp) gene (M63, 879 nt; BS58, 841 nt) and the whole coat protein coding region (1545 nt and 1436 nt) along with the complete VP2 protein coding region (656 nt and 668 nt). The calicivirus conserved aa motifs GLPSG and YGDD (Smiley et al., 2002) were identified in the RdRp region. A third calicivirus RdRp motif DYXXWDST, was also well conserved in both viruses, exhibiting the following aa sequence DFSKWDST (Wolf et al., 2011) . Strain BtCalV/M63/HUN/2013, detected in a M. daubentonii bat, was grouped with other bat-derived sapoviruses and showed the closest relationship with a Chinese sequence (GenBank: KJ641703), detected from the feces of a Myotis myotis bat (Fig. 2a) , with 67% nt and 67% aa identity in the RdRp and 66% nt and 61% aa identity in the capsid gene. Strain BtCalV/BS58/HUN/2013, derived from an E. serotinus bat species, was grouped with a bat calicivirus and with porcine caliciviruses of the proposed genus Valovirus. The closest phylogenetic relationship was described with a Chinese bat calicivirus strain, BtMspp.-CalV/GD2012 (GenBank: KJ641700), detected in an undefined Myotis species. Identity was 63 nt and 60% aa in the RdRp and 63% nt and 56% aa in the capsid gene between strain BS58/HUN/ 2013 and strain BtMspp.-CalV/GD2012 (Table 1) . This geographically distant evolutionary relationship among European and Asian origin bat viruses has already been described when picornaviruses and parvoviruses detected in European Miniopterus schreibersii bats were characterized (Kemenesi et al., 2015a,b) . Our new results with related caliciviruses in Myotis and Eptesicus bats, widens the list of bat species harboring genetically related viruses in large geographic distances, which can be a result of ancient evolutionary history of bats and their viruses and presumably could be a result of biogeographical history of these bat species.

In depth structural analysis of the capsid of caliciviruses has been reported in previous studies (Ossiboff et al., 2010; Prasad et al., 1999) . There was 27.3% aa identity, and 42.3% aa similarity between the template FCV capsid protein (PDB ID: 3M8L) and the mature CP of the bat calicivirus strain BtCalV/M63/HUN/2013. Likewise, there was 34.7% aa identity and 47.6% similarity, between the template Norwalk virus CP (PDB ID: 1IHM) and the CP of the bat calicivirus strain BtCalV/BS58/ HUN/2013. These identity and similarity values allowed reliable homology modeling and reconstruction of the capsid and virion structure of the two viruses to be made. The mature calicivirus capsid protein has three main domains: S, P1 and P2 ( Fig. 1c and d) . The S (i.e. shell) domain is composed of a canonical eight-stranded β-barrel structure. This domain forms the inner part of the virion and it has the most conserved structure within caliciviruses. The P1 and P2 subdomains build up the middle and the outer protruding part of the virion, respectively. The hypervariable regions, the primary targets of neutralizing antibodies, are located on the P2 subdomain (Ossiboff et al., 2010) . The basic structure together with the recognized main domains was predicted to be conserved in the two Hungarian bat caliciviruses. However in the P2 region we detected three deletions and an insertion (Fig. 1a and c) in the CP of the bat calicivirus BtCalV/M63/HUN/2013 when compared to the FCV template sequence. In the case of the strain BtCalV/ BS58/HUN/2013, we found two insertions and two deletions on the P2 subdomain ( Fig. 1b and d) with respect to the template structure. Insertions and/or deletions in the P2 domain are common even among caliciviruses of the same genus and/or genotype (Pinto et al., 2012; Medici et al., 2015) . These small structural and sequential differences may play important roles in the P2 subdomain dual functions. The protruding regions of the virion may be involved as antigens in the humoral immune response (Fig. 1e and f) . Furthermore, this variable part of the P2 subdomain may take part in protein-protein interactions between the CP and its functional cell receptor. Additionally, these regions are responsible for the fine tuning of cell-virus interaction during cell entry (Ossiboff et al., 2010) .

In conclusion, genome sequencing, phylogenetic analysis and homology modeling allowed a more precise characterization of two novel bat caliciviruses. One virus could be firmly assigned to the Sapovirus genus. Sapoviruses are enteric viruses associated with gastro-enteritis in humans and with enteritis in pigs (Oka et al., 2015) . These viruses have been also identified in diarrheic dogs and cats (Soma et al., 2015; Li et al., 2011) and display an impressive genetic heterogeneity (Oka et al., 2015) . Interestingly, genetic heterogeneity is also displayed by the few bat sapoviruses identified thus far, thus posing a challenge for the design of specific diagnostic tools. Based on our phylogenetic analyses species-specific patterns are clear among sapoviruses of different animal origin (Fig. 2b) . However, high phylogenetic diversity can be observed within these species-related clusters with several separate branches, which indicate a long-term evolutionary diversification. Neither amino acid and nucleic acid sequence homology patterns between the examined sequences (Table 1 .) nor the phylogenetic results indicate evidence for possible recombination event related to our sequences.

The other bat calicivirus was genetically more similar to porcine caliciviruses of the putative novel genus Valovirus. Valoviruses are orphan viruses, detected mostly in finisher pigs with a low prevalence (L'Homme et al., 2009; Wang et al., 2011; Di Martino et al., 2011) . Antibodies specific for valoviruses have been detected in 63 (10.3%) of 614 porcine serum samples in a previous study (Di Martino et al., 2012) , indicating that these viruses circulate in pigs, although their pathogenic role is unclear. Although these viruses were described in asymptomatic swine from several geographically distant countries, such as Japan and Canada (Sato et al., 2014; L'Homme et al., 2009) , the circumstances of their worldwide distribution should be extensively examined in future. The genomic organization of the bat calicivirus was not identical to the organization observed in valoviruses. The closest phylogenetic relatedness was observed with valoviruses, however the novel bat calicivirus sequence (BtCalV/BS58/HUN/2013) clearly formed a separate branch.

The limited number of studies describing bat caliciviruses is providing evidence that these viruses are highly heterogeneous genetically. It is possible that as data will accumulate, this picture will be confirmed and the interactions between caliciviruses from bats and from other animal species, including humans, will be more evident. Also, this will help clarify whether caliciviruses may have a pathogenic role in bats and will help reveal the evolutionary pathways of these viruses across the various bat species and geographic areas.

Genome sequence assembly using trace signals and additional sequence information

Detection of St-Valerien-like viruses in swine

Seroprevalence of St-Valerien-like caliciviruses in Italian swine

Characterization of a rhesus monkey calicivirus representing a new genus of Caliciviridae

Norovirus GII.4/Sydney/2012 in Italy

VMD -visual molecular dynamics

Design and evaluation of a primer pair that detects both Norwalk-and Sapporo-like caliciviruses by RT-PCR

Molecular survey of RNA viruses in Hungarian bats: discovering novel astroviruses, coronaviruses, and caliciviruses

Genetic diversity and recombination within bufaviruses: detection of a novel strain in Hungarian bats

Genetic characterization of a novel picornavirus detected in Miniopterus schreibersii bats

European bats as carriers of viruses with zoonotic potential

Genomic characterization of swine caliciviruses representing a new genus of Caliciviridae

Viruses in diarrhoeic dogs include novel kobuviruses and sapoviruses

Complete genome sequence of a novel calicivirus from a goose

Identification of the novel Kawasaki 2014 GII.17 human norovirus strain in Italy

Candidate new rotavirus species in sheltered dogs

Characterization of a novel calicivirus causing systemic infection in atlantic salmon (Salmo salar L.): proposal for a new genus of caliciviridae

Comprehensive review of human sapoviruses

Conformational changes in the capsid of a calicivirus upon interaction with its functional receptor

Discovery and genomic characterization of noroviruses from a gastroenteritis outbreak in domestic cats in the US

X-ray crystallographic structure of the Norwalk virus capsid

I-TASSER: a unified platform for automated protein structure and function prediction

Characterization of St-Valerien-like virus genome detected in Japan

Prevalence of porcine noroviruses, molecular characterization of emerging porcine sapoviruses from finisher swine in the United States, and unified classification scheme for sapoviruses

Schrödinger Suite, Schrödinger. LLC

Characterization of an enteropathogenic bovine calicivirus representing a potentially new calicivirus genus

Calicivirus emergence from ocean reservoirs: zoonotic and interspecies movements

Detection of Norovirus and Sapovirus from diarrheic dogs and cats in Japan

A molecular phylogeny for bats illuminates biogeography and the fossil record

Discovery and genomic characterization of a novel bat sapovirus with unusual genomic features and phylogenetic position

Characterization and prevalence of a new porcine calicivirus in swine, United States

Discovery of novel virus sequences in an isolated and threatened bat species, the New Zealand lesser short-tailed bat (Mystacina tuberculata)

Genetic characterization of a novel calicivirus from a chicken

Bats and viruses: friend or foe? PLoS Pathog. 10, e1003651. Zhang, Y., 2008. I-TASSER server for protein 3D structure prediction

Research activity of Ferenc Jakab was supported by the TÁMOP 4.2.4. A/2-11-1-2012 0001 -National Excellence Program elaborating and operating an inland student and researcher personal support system. The project was subsidized by the European Union and co-financed by the European Social Fund. Krisztián Bányai was supported by the "Momentum program". Ákos Gellért and Szilvia Marton were the recipients of a János Bolyai fellowship from the Hungarian Academy of Sciences. The licensing of the Schrödinger Suite software package was funded from the OTKA grant under agreement no.:108793. The present scientific contribution is dedicated to the 650th anniversary of the foundation of the University of Pécs, Hungary.