key: cord-299827-lxyo3z3u authors: Di Martino, Barbara; Di Profio, Federica; Melegari, Irene; Sarchese, Vittorio; Cafiero, Maria Assunta; Robetto, Serena; Aste, Giovanni; Lanave, Gianvito; Marsilio, Fulvio; Martella, Vito title: A novel feline norovirus in diarrheic cats date: 2015-12-29 journal: Infect Genet Evol DOI: 10.1016/j.meegid.2015.12.019 sha: doc_id: 299827 cord_uid: lxyo3z3u By screening a collection of fecal samples from young cats housed in three different shelters in South Italy, noroviruses (NoVs) were found in 3/48 (6.2%) specimens of animals with enteritis signs while they were not detected in samples collected from healthy cats (0/57). Upon sequence analysis of the short RNA-dependent RNA polymerase (RdRp) region, the three strains displayed the highest nucleotide (nt) and amino acid (aa) identities to the prototype GIV.2 strain lion/Pistoia/387/06/ITA (91.0–93.0% nt and 97.0–98.0% aa). The sequence of ~ 3.4-kb portion at the 3′ end of the genome of a NoV strain, TE/77-13/ITA, was determined. In the full-length ORF2, encoding the VP1 capsid protein, the virus was genetically closest to the canine GVI.2 NoV strains C33/Viseu/2007/PRT and FD53/2007/ITA (81.0–84.0% nt and 93.0–94.0% aa identities), suggesting a recombination nature, with the cross-over site being mapped to the ORF1-ORF2 junction. Based on the full-length VP1 amino acid sequence, we classified the novel feline NoV, together with the canine strains Viseu and FD53, as a genotype 2, within the genogroup GVI. These findings indicate that, as observed for GIV NoV, GVI strains may infect both the canine and feline host. Unrestricted circulation of NoV strains in small carnivores may provide the basis for quick genetic diversification of these viruses by recombination. Interspecies circulation of NoVs in pets must also be considered when facing outbreaks of enteric diseases in these animals. Noroviruses (NoVs), Caliciviridae family, have been identified as the most common cause of viral gastroenteritis in humans. NoV infections affect persons of all age groups and are predominantly transmitted through the fecal-oral route, either indirectly through contaminated food, water or surfaces or directly from person to person (Patel et al., 2008) . Drop virions are nonenveloped and approximately 30 to 35 nm in diameter. The RNA genome is organized into three open reading frames (ORFs) (Green, 2007) . ORF1 encodes a polyprotein that is cleaved by the virus-encoded protease to produce several nonstructural proteins, including the RNA dependent RNA polymerase (RdRp), ORF2 encodes a major capsid protein (VP1) and ORF3 encodes a small basic protein (VP2) that has been associated with the capsid stability (Bertolotti-Ciarlet et al., 2003) . Based on the full-length VP1 amino acid sequence, NoVs have been divided into six genogroups (GI to GVI) and several genotypes (Zheng et al., 2006; Martella et al., 2009; Green, 2013) . Only GI, GII, and GIV NoVs infect humans, with GII strains being the most prevalent worldwide (Green, 2007) . NoVs genetically similar to human NoVs have been recently found in dogs and cats (Martella et al., 2007 (Martella et al., , 2008 Summa et al., 2012; Pinto et al., 2012; Soma et al., 2014) , raising public health concerns of potential cross-species transmission due to the strict social interaction between humans and pets. Feline NoVs were first detected in the stools of 8-12-week-old kittens from a feline shelter with an outbreak of diarrhea in New York State (Pinto et al., 2012) . In the VP1 encoding gene, the feline NoVs displayed the highest amino acid (aa) identity (97.9%) to the prototype NoV strain GIV.2/Pistoia/387/06/ITA, detected in a captive lion cub with severe hemorrhagic enteritis (Martella et al., 2007) and to the canine strain GIV.2/Bari-170/07/IT (90.4% aa), detected in a young dog with diarrhea (Martella et al., 2008) . Using baculovirus-expressed VP1 of the lion NoV strain GIV.2/Pistoia/387/06/ITA, antibodies specific for GIV NoVs have been identified in 16.1% of cats in Italy (Di Martino et al., 2010) , providing indirect evidence for the circulation of these NoVs in felines. In addition, the RNA of GIV.2 NoVs has been detected in 1.2% of fecal samples of cats with enteritis in Japan (Soma et al., 2014) . Upon genome sequencing, the feline NoV strain cat/GVI.1/JPN/2012/M49 (Takano et al., 2015) was found to be more similar (87.0% aa identity) in the full-length ORF2 to the canine NoV GVI.1/Bari/91/07/IT (Martella et al., 2009) . Infection, Genetics and Evolution 38 (2016) Infection, Genetics and Evolution j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / m e e g i d Altogether these findings indicate that diverse NoV strains may infect cats, as observed in dogs, and that the feline and canine host may be infected by the same NoV strains, thus constituting an enlarged host reservoir for these animal NoVs. In order to draw a more complete picture of NoVs molecular epidemiology in cats, in this study a collection of fecal specimens from diarrheic and healthy animals was screened using either broadly-reactive primers for caliciviruses and primers specific for NoVs. A total of 105 stool samples from domestic cats aged 2-12 months were collected from April to July 2013 in three different shelters located in South Italy. The fecal panel consisted of 48 samples from cats with signs of mild to severe gastroenteritis and 57 samples from asymptomatic animals. All the samples were stored at −80°C until use. Fecal specimens (10%) were re-suspended in phosphate-buffered saline pH 7.2, and the debris were removed by centrifugation at 8000 × g for 5 min. DNA and RNA extracts were prepared using the DNAeasy® and QIAamp® viral RNA kit (Qiagen GmbH, Hilden, Germany), according to the manufacturer's instructions and stored at −80°C until use. To assess the presence of NoV RNA, the samples were screened using a broadly reactive primer pair, p289-p290, targeted to highly conserved motifs DYSKWDST and YGDD of the RNA-dependent RNA polymerase (RdRp) region of the polymerase complex (Jiang et al., 1999) . In the samples yielding amplicons of the expected sizes, the presence of NoV was confirmed using norovirus-specific primer pair JV12Y-JV13I (Vennema et al., 2002) . All the fecal samples were also tested by PCR or RT-PCR for feline parvovirus (FPV) (Buonavoglia et al., 2001) , feline enteric coronavirus (FECV) (Gunn-Moore et al., 1998) and feline kobuviruses (FeKoV) (Di Martino et al., 2015) . The amplicons were excised from the gel and purified using a QIAquick gel extraction kit (Qiagen GmbH, Hilden, Germany). The fragment was then subjected to direct sequencing using BigDye Terminator Cycle chemistry and 3730 DNA Analyzer (Applied Biosystems, Foster, CA). Basic Local Alignment Search Tool (BLAST; http://www.ncbi.nlm. nih.gov) and FASTA (http://www.ebi.ac.uk/fasta33) with default values were used to find homologous hits. The sequence of~3.4-kb fragment at the 3′ end of the genome of one such strain, TE/77-13/ITA, including the partial RdRp gene and the complete ORF2 and ORF3 genes, was determined by 3′ RACE protocol, as previously described (Scotto-Lavino et al., 2006) . cDNA was synthesized by SuperScript III First-Strand cDNA synthesis kit (Invitrogen Ltd., Milan, Italy) with primer QT. PCR was then performed with TaKaRa La Taq polymerase (Takara Bio Europe S.A.S. Saint-Germain-en-Laye, France) with forward primer p290 and reverse primers QO and QI. Finally, the amplicons were purified and cloned by using TOPO® XL Cloning Kit (Invitrogen Ltd., Milan, Italy). Additional primers were designed to determine the complete 3.4 kb sequences by an overlapping strategy (Table 1) . Sequence editing and multiple alignments were performed with the BioEdit software package, version 2.1 (Hall, 1999) . Phylogenetic trees were generated using Bayesian analysis with MrBayes (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003) . The appropriate substitution model settings were derived using jModelTest (Posada, 2009 ). The sequence obtained was analyzed with Simplot (Lole et al., 1999) using a window size of 200 and step size of 20, with gap strip off and Hamming correction on. Additionally, recombination analysis was carried out with different algorithms implemented in the Recombination Detection Program v.4.43 (RDP4) (Martin et al., 2010) , with default settings. Out of 105 samples, 3 (2.8%) contained NoVs RNA, either alone (0.9%, 1/105) or in mixed infections with FeKoV or FECV (1.9%, 2/105). Sixteen samples (15.2%) were found to contain FPV DNA alone. All the NoV positive samples were identified from diarrheic cats with a prevalence rate of 6.2% (3/48), while they were not detected from asymptomatic animals (0/57). By sequence comparison in the short RdRp fragment, the viruses TE/68-13, ME/78-13 and TE/77-13/ITA shared 89.9-93.2% nt and 90.8-95.3% aa identities to each other and displayed the highest identity (91.0-93.0% nt and 97.0-98.0% aa) to the prototype GIV.2 strain lion/Pistoia/387/06/ITA. For the strain TE/77-13/ITA the sequence of~3.4-kb fragment at the 3′ end of the genome, including the partial RdRp (0.8 kb) and the complete ORF2 and ORF3 (GenBank accession number KT245136), was sequenced and the genome organization was determined. Phylogenetic analysis was based on the 750-nt sequence of the COOH terminus of the polymerase complex of the carnivore NoV strains available in the databases. Also, RdRp sequences of human GIV.1 NoVs were included in the analysis and used to calculate a nt identity matrix. By visual inspection of the tree, the carnivore NoVs segregated in at least three different genetic clusters (Fig. 1) . The strain cat/TE/77-13/ITA was grouped with the feline NoV strains lion/GIV.2/Pistoia/387/06/ITA, cat/GIV.2/ CU081210E/USA/2010 and cat/GVI.1/JPN/2012/M49 (Martella et al., 2007; Pinto et al., 2012; Takano et al., 2015) , with a nt identity of 90.8-94.0%. This group shared identity of 82.3-85.4% to the recombinant NoV dog/GVI.1/91/2007/ITA and to the GIV.2 strains dog/170/07/ ITA and dog/Thessaloniki/30/08/GR (Martella et al., 2008 (Martella et al., , 2009 Ntafis et al., 2010) , which in turn segregated in a second cluster (91.6-97.8% nt identity). A minor group, distantly related to the feline NoVs (78.2-79.6% nt identity) and to the canine NoVs GVI. Two additional clusters were resolved in the tree that included, respectively, the human GIV.1 strains detected from stool and sewage samples in different geographic settings (Fankhauser et al., 2002; La Rosa et al., 2008; Eden et al., 2012; Ao et al., 2014; Han et al., 2014) and GIV.1 NoVs found in sewage samples in Italy in -2009 (La Rosa et al., 2010 . The nt identity between this two groups was 87.0-91.5%. The ORF2 of strain TE/77-13/ITA was 1719 nt in length and encoded a VP1 capsid protein with a predicted size of 572 aa. ORF3 was 795 nt long and encoded a VP2 protein of 264 aa. A 16-nucleotide (nt) overlap was present in the ORF1-ORF2 junction region. In the complete VP1, the strain TE/77-13/ITA was most closely related (81.0-84.0% nt and 93.0-94.0% aa) to the strains dog/GVI.2/C33/ Viseu/2007/PRT and dog/GVI.2/FD53/2007/ITA, while identity to the feline strain cat/NoV/GVI.1/JPN/2012/M49 and to the canine strain GVI.1/ Bari/91/2007/ITA (Martella et al., 2009 ) was 70.0-71.0% nt and 80.0-81.0% aa. Strain TE/77-13/ITA displayed b54.0% aa identity to animal and human GIV NoVs. Phylogenetic analysis was performed with a selection of complete capsid sequences representative of the Norovirus genus. In the VP1-based tree (Fig. 2) , strain TE/77-13/ITA segregated with the canine NoVs GVI.2/C33/Viseu/2007/PRT and GVI.2/FD53/ 2007/ITA into genogroup GVI, genotype 2. A nucleotide identity plot of the genome of strain TE/77-13/ITA was elaborated, in comparison with the canine strain dog/NoV/GVI.2/C33/Viseu/2007/PRT and the feline strain cat/NoV/GIV.2/CU081210E/USA/2010. By SimPlot (Fig. 3) and RDP (Fig. S1 ) analyses, a putative recombination break-point event was mapped to ORF1-ORF2 junction region at nt 703 with a significant statistical support (P b 0.05). In this study direct evidence was collected for circulation of NoVs in cats. NoVs were detected in cats with enteric signs while they were not identified in samples collected from healthy animals used as control study group. Experimental inoculation of specific pathogen free cats with the feline GVI.1 strain JPN/2012/M49 can induce enteritis signs, diarrhea and vomiting (Takano et al., 2015) . Although the pathogenic role of NoV in cats should be confirmed in larger epidemiological studies 2007/ITA. Following strictly the outlines of Zheng's classification (Zheng et al., 2006) , we classified the novel feline NoV, together with the canine strains Viseu and FD53, as a genotype 2 (N 85% pairwise aa identity intergenotypes), within the genogroup GVI (N 55% pairwise aa identity intergenogroups). Accordingly, cats and dogs may harbor NoVs of the same genotypes, GIV.2, GVI.1 and GVI.2. Circulation of NoVs genetically related in different host species has been already demonstrated. Porcine NoVs cluster in GII (Wang et al., 2005) , but within different genotypes (GII.11, GII.18, and GII.19 ) from those infecting humans. GIII NoVs have been detected in large and small ruminants, with GIII.1 and GIII.2 (Liu et al., 1999; Oliver et al., 2007) found in cattle and GIII.3 in sheep (Wolf et al., 2009) . However, unlike small carnivores, circulation of NoV strains belonging to the same genotypes in heterologous species has not been reported thus far. This intriguing finding poses several questions. Binding of GVI.2 and GVI NoVs in dog tissues seems to be mediated by the presence of the H and A antigens of the histo-blood group antigen (HBGA) family (Caddy et al., 2014) and therefore to be genetically determined, as observed in humans (Marionneau et al., 2002) . This may suggests that dogs and cats share a similar pattern of HBGAs as attachment factor for NoV infections. Also, virus-like particles of seven different human NoV genotypes (GI.1, GI.2, GI.3, GII.3, GII.4, GII.6, and GII.12) have been shown to be able to bind to canine gastrointestinal tissues (Caddy et al., 2015) . It will be interesting to assess whether cats may also be infected by human NoVs, as observed in dogs (Summa et al., 2012; Caddy et al., 2015) , as this may have implications for the transmission of human NoVs. Recombination among NoVs of domestic carnivores has been already described. The canine NoV strain Bari/91/07/ITA, resembles GIV.2 NoVs in its polymerase gene while it is genetically unrelated in the VP1 gene to GIV NoV (Martella et al., 2009) . The feline NoV strain, JPN/2012/M49 (Takano et al., 2015) , possesses a GIV.2 RdRp region, and a GVI.1 ORF2 related to the canine virus Bari/91/07/ITA (87.0% aa identity). In all the cases, the site of recombination was mapped to the ORF1/ORF2 junction region. This part of NoV genome is highly conserved and has been individuated as a preferential recombination site (Bull et al., 2005) . Recombination has been shown to strongly influence the evolution and epidemiology of human NoVs (Ambert-Balay et al., 2005; Reuter et al., 2006) and surely poses a challenge for the development of specific diagnostic tools for NoV of carnivores. Analysis of the RdRp fragment cannot be used to characterize unequivocally these animal NoVs and a definitive characterization should rely on the ORF2. The development of molecular assays for caliciviruses and for NoVs has allowed gathering epidemiological information about these viruses in several animal species, including domestic carnivores. It is now clear that cats and dogs may harbor NoVs of several genotypes and genogroups, although the clinical relevance of these viruses remains to be investigated. Gathering information on the genetic diversity of animal NoVs will be useful to optimize/develop direct and indirect diagnostic tools, and to investigate more effectively the epidemiology of NoVs in carnivores. In addition, as NoVs of carnivores are suspected to have a zoonotic relevance (Peasey et al., 2004; Mesquita et al., 2013; Di Martino et al., 2014; Caddy et al., 2015) , this will be useful to understand the extent of inter-species transmission from cat to dogs, and vice versa, from carnivores to humans. Supplementary data to this article can be found online at http://dx. doi.org/10.1016/j.meegid.2015.12.019. All Authors declare that there are no financial or other relationships that might lead to a conflict of interest. All authors have seen and approved the manuscript and have contributed significantly to the work. Characterization of new recombinant noroviruses Detection of human norovirus GIV.1 in China: a case report The 3′ end of Norwalk virus mRNA contains determinants that regulate the expression and stability of the viral capsid protein VP1: a novel function for the VP2 protein Norovirus recombination in ORF1/ORF2 overlap Evidence for evolution of canine parvovirus type 2 in Italy Genogroup IV and VI canine noroviruses interact with histo-blood group antigens Evidence for human norovirus infection of dogs in the United Kingdom Seroprevalence of norovirus genogroup IV antibodies among humans Detection of feline kobuviruses in diarrhoeic cats Detection of antibodies against norovirus genogroup GIV in carnivores Complete genome of the human norovirus GIV.1 strain Lake Macquarie Epidemiologic and molecular trends of "Norwalk-like viruses" associated with outbreaks of gastroenteritis in the United States Caliciviridae: the noroviruses Caliciviridae: the noroviruses Detection of feline coronaviruses by culture and reverse transcriptase-polymerase chain reaction of blood samples from healthy cats and cats with clinical feline infectious peritonitis BioEdit: a user-friendly biological sequence alignment and analysis program for Windows 95/98/NT Detection of norovirus genogroup IV, klassevirus, and pepper mild mottle virus in sewage samples in South Korea MRBAYES: Bayesian inference of phylogeny Design and evaluation of a primer pair that detects both Norwalk-and Sapporo-like caliciviruses by RT-PCR Molecular detection and genetic diversity of norovirus genogroup IV: a yearlong monitoring of sewage throughout Italy Detection of genogroup IV noroviruses in environmental and clinical samples and partial sequencing through rapid amplification of cDNA ends Molecular characterization of a bovine enteric calicivirus: relationship to the Norwalk-like viruses Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination Norwalk virus binds to histo-blood group antigens present on gastroduodenal epithelial cells of secretor individuals Norovirus in captive lion cub Genetic heterogeneity and recombination in canine noroviruses Detection and molecular characterization of a canine norovirus RDP3: a flexible and fast computer program for analyzing recombination Novel norovirus in dogs with diarrhea Presence of antibodies against genogroup VI norovirus in humans Outbreak of canine norovirus infection in young dogs Complete genomic characterization and antigenic relatedness of genogroup III, genotype 2 bovine noroviruses Systematic literature review of role of noroviruses in sporadic gastroenteritis Seroepidemiology and risk factors for sporadic norovirus/Mexico strain Discovery and genomic characterization of noroviruses from a gastroenteritis outbreak in domestic cats in the US Selection of models of DNA evolution with jModelTest Epidemic spread of recombinant noroviruses with four capsid types in Hungary MrBayes 3: Bayesian phylogenetic inference under mixed models 3′ end cDNA amplification using classic RACE Detection of norovirus and sapovirus from diarrheic dogs and cats in Japan Pet dogs-a transmission route for human noroviruses? Molecular characterization and pathogenicity of a genogroup GVI feline norovirus Rational optimization of generic primers used for Norwalk-like virus detection by reverse transcriptase polymerase chain reaction Porcine noroviruses related to human noroviruses Molecular detection of norovirus in sheep and pigs in New Zealand farms Norovirus classification and proposed strain nomenclature This work was supported with funds from the grant "Calicivirus nei carnivori e nell'uomo: caratterizzazione molecolare, epidemiologia, implicazioni zoonosiche -PRIN 2008" (20084F4P7X_003).