key: cord-0830879-xo793bf5
authors: Karayel-Hacioglu, Ilke; Alkan, Feray
title: Molecular characterization of bovine noroviruses and neboviruses in Turkey: detection of recombinant strains
date: 2019-03-08
journal: Arch Virol
DOI: 10.1007/s00705-019-04186-0
sha: 9ca31b72b40475e937c13e6617140ecbd6060480
doc_id: 830879
cord_uid: xo793bf5

To investigate the molecular epidemiology and genetic diversity of bovine enteric caliciviruses, a total of 167 fecal samples from diarrheic calves were screened. Bovine noroviruses (BoNoVs) and neboviruses were detected in 56 (33.5%) and 37 (22.1%) fecal samples, respectively. Sequences of the RdRp and capsid gene of selected BoNoVs showed that the GIII.1 and GIII.2 genotypes were in circulation in Turkey. Two of the BoNoV strains were identified as recombinant strains (GIII.P1/GIII.2). All examined neboviruses possessed a Nebraska-like RdRp gene. The two nebovirus strains were classified into lineage 4 based on phylogenetic analysis of VP1 amino acid sequences. One of them showed evidence of a recombination event within the S domain. This study is thus the first to reveal the presence of the BoNoV GIII.1 genotype and recombinant strains of BoNoV and neboviruses in Turkey. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-019-04186-0) contains supplementary material, which is available to authorized users.

with the following prototypes: Bo/Newbury1/76/UK, Bo/ Nebraska/80/US and Bo/DijonA216/06/FR [1, 6, 9] . The capsid gene of neboviruses can be divided into two genotypes. One genotype consists of strains Bo/LZB-1/17/CH and Bo/DijonA216/06/FR, and the other can be divided into four capsid lineages [9] [10] [11] .

Recombination within the same genogroup or genotype, or between genogroups occurs frequently in caliciviruses, and this is recognized as an important mechanism for its continual emergence [12] . Recombination events in the norovirus genotype usually occur close to or within the ORF1/ ORF2 junction [13] . Molecular studies have revealed that the BoNoV recombinant types represent two RdRp/capsid combinations, GIII.P1/GIII.2 and GIII.P2/GIII.1 [13] . The first recombinant BoNoV strain, Bo/Thirsk10/00/UK, had a GIII.1-related RdRp but a GIII.2-related capsid [14] . Molecular investigations of isolates from various geographical regions have shown that most naturally occurring BoNoV recombinants share the same combination with Bo/Thirsk10/00/UK [12, 15] , whereas only one recombinant strain with a GIII.P2/GIII.1 combination has been described to date [13] . Recombination events can also occur in neboviruses. Indeed, two recombinant nebovirus strains (Bo/16TE/2010/ITA and Bo/LZB-1/17/CH) have been identified in recent studies [11, 16] .

Since bovine enteric caliciviruses (BECs) are routinely excluded from diagnostic procedures for identifying enteric diseases, their real impact on livestock production is unclear.

Unlike other enteric viruses, such as bovine rotavirus group A (BRVA), and bovine coronavirus (BCoV), there is limited knowledge of BECs in Turkey [17] [18] [19] [20] . Before the conception of this study, all reports on BEC infections only addressed the detection of BoNoVs in fecal samples from diarrheic calves [19] . Thus, to investigate the molecular epidemiology and genetic diversity of BECs, a total of 167 fecal samples from diarrheic calves (ages 1 day to 7 months) were collected from 48 family farms (up to 10 cattle) and 11 organized dairy cattle farms in 10 different provinces of Turkey between 2002 and 2016 (Table S1 and Table S2 ).

Coinfection is frequently observed in diarrheic calves, although a single primary pathogen can be the cause in some cases. Therefore, the fecal samples were also screened for other enteric pathogens frequently detected in Turkey, such as BRVA, BCoV, E. coli and Cryptosporidium spp., using a commercial ELISA kit (Bio-X Diagnostics S.A., Belgium) according to the manufacturer's instructions.

Viral RNA was extracted from 250 µl of fecal suspension (1:10, w/v) using the TRIzol LS procedure. All extracts were stored at -80 °C until tested. Reverse transcription was performed using a RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, USA) according to the manufacturer's instructions. BECs were detected by nested RT-PCR using two sets of external and internal primers targeting the RdRp gene. The following primer sets were used: for BoNoVs, CBECU-F/R (532 bp) and nCBECU-F/R (326 bp); for neboviruses, NBU-F/R (549 bp) and nNBU-F/R (194 bp) [21, 22] . Nested RT-PCR was performed as described previously [21, 22] . The capsid gene of the BoNoVs strains was amplified using two primer sets (NAcap-F/R and nCBECU-F/BoNoVcap-A) [12, 15, 22] . The primers NAcap-F/R (1615 bp) targeted the full-length capsid gene of BoNoVs, while the primer set nCBECU-F/BoNoVcap-A (781 bp) was used to amplify the 3′ end of RdRp and the 5′ end of the capsid gene. The primer set used for the amplification of the full-length nebovirus capsid gene (1692 bp) was NBcap-F3/R [12] . RT-PCR was performed using a Phusion High-Fidelity PCR Kit (Thermo Fisher Scientific, USA) as described previously [12, 15] with some modifications. The amplification products were analyzed by 1% agarose gel electrophoresis and visualized under UV light.

Some PCR products with the expected size were selected for sequencing based on the sampling year and the location. Purified PCR products were sequenced in both directions with the same primers used for amplification. Cognate sequences of reference BoNoVs and neboviruses representing different genotypes were retrieved from the GenBank database using the BLAST engine. Multiple sequence alignments were prepared by the MUSCLE algorithm as implemented in Aliview Software [23, 24] . MEGA 7.0 was used to build phylogenetic trees [25] , the statistical significance of which was estimated by bootstrap analysis (1000 replicates). Nucleotide (nt) and amino acid (aa) identity values were calculated in CLC Main Workbench (CLC Bio-QIAGEN, Aarhus, Denmark). Genetic recombination was analyzed using SimPlot software (version 3.5.1). All sequences (> 200 bp) detected in this study have been made available in the GenBank database under the following accession numbers: KX184721-KX184724; KX216596, KX216597; KX189067-KX189100.

Among the 167 samples from diarrheic calves, BECs were detected in 76 (45.5%). BoNoVs and neboviruses were detected in 56 (33.5%) and 37 (22.1%) fecal samples, respectively. Coinfection with BoNoVs and neboviruses was detected in 17 (10.1%) samples. BECs were detected on 20 family farms and 10 commercial dairy cattle farms, with 15 family farms and 8 commercial farms positive for BoNoVs, and 11 family farms and 7 commercial farms positive for neboviruses (Table S1 and Table S2 ). The results from the selected representatives for different farms and geographical regions showed that BECs are prevalent in Turkey (Table S3 and Table S4 ). Out of 76 BEC-positive samples, 18 and 14 were positive for the BoNoV and nebovirus alone, respectively, while the rest of the 44 fecal samples tested positive for other enteric pathogens. Of the concurrent infections of BECs with other pathogens, Cryptosporidium spp. was most commonly seen (46.5%), followed by BRVA (31.9%) and BCoV (20.8%). However, no E. coli was detected in any of the samples.

Molecular analysis of the partial RdRp gene sequences obtained from the selected 34 BoNoV strains revealed that nine BoNoV strains were clustered in GIII.1, whereas the remaining 25 BoNoV strains were located in GIII.2 (Fig. 1A) 1B) .

Considering that recombination frequently occurs at the ORF1-ORF2 junction of noroviruses, the capsid gene was amplified. Amplicons of the expected size for the full-length capsid gene for four strains (Aksaray/K59, Aksaray/K063, Aksaray/K056, and Aksaray/K060) and a partial capsid gene for one (Aksaray/K066) were successfully obtained. The sequences of all full-length capsid genes shared higher nt and aa sequence identity with Bo/Newbury2/76/UK (82.8-88.5% nt, 91.1-96.9% aa) than with Bo/Jena/80/DE (64.8-68.5% nt and 66.6-69.7% aa). Looking at the phylogenetic analysis, these four sequences were more closely related to GIII.2 ( Fig. 2A) , whereas Aksaray/K066 was clustered in GIII.1 (Fig. S1) .

Based on the phylogenetic analysis (ORF1/RdRp and ORF2/Capsid), the conflicting results obtained in grouping the RdRp and capsid genes of Aksaray/K056 and Aksaray/ K060 strains may indicate recombination within the ORF1 and ORF2 junction. Accordingly, both strains displayed the same RdRp type as GIII.1, whereas their capsid types were genetically closest to GIII.2. Based on a SimPlot analysis performed on the sequence of the partial ORF1 and the fulllength ORF2, a putative recombination breakpoint located within the ORF1/ORF2 junction (nt position 327) was identified (Fig. 2B) .

For neboviruses, a full-length capsid gene for two strains was successfully recovered (AksarayK056 and (Fig. 3A) . The two strains with a lineage 4 capsid type, both of which were related to reference virus Bo/Nebraska/80/US (lineage 1) based on the RdRp gene, were further analyzed for possible recombination events. Recombination analysis did not show any recombination events in the genome of AksarayK056. Strain AnkaraA4 closely resembled Nebraska in RdRp but resembled Bo/Penrith143/00/UK, Bo/16TE/2010/ITA and Bo/DijonA386/08/FR (lineage 4) in the capsid region. The UK strain has a Newbury1-like RdRp, while Bo/16TE/2010/ ITA and Bo/DijonA386/08/FR have a Nebraska-like RdRp. Accordingly, AnkaraA4 was a possible recombinant strain, and a putative recombination breakpoint was identified within the S domain of the capsid protein (Fig. 3B) .

The involvement of BECs in calf diarrhea has become an important issue in recent years, and the prevalence of these viruses has therefore been investigated in many countries. In this study, we identified BoNoV (33.5%) and nebovirus (22.1%) infections in diarrheic calves from 48 family farms and 11 commercial dairy cattle farms in 10 different provinces of Turkey. Surprisingly, BECs were detected on commercial dairy farms at a much higher rate than on family farms. This could be due to the large number of animals on commercial farms and/or routine vaccination against the other enteropathogens. Further epidemiological studies with BoNoVs were detected in diarrheic calves at a much higher rate than previously seen in Turkish reports (1.7-8.5%) [18] [19] [20] but at a similar rate to that reported in European countries and North America [9, 21, 26] . The prevalence rate of nebovirus was higher than the reported rates in Tunisia [27] , Brazil [28] , Italy [16] , France [9] , South Korea [29] and the United Kingdom [30] , while it was similar to prevalence rates (21.6-39.5%) reported in North America and Iran [21, 26, 31] . There has been only one other study on the prevalence of neboviruses (25.19%) in Turkey [20] , and the similarity in its reported prevalence rates to those seen here is noteworthy. The higher prevalence rate of BoNoV compared to that of nebovirus in this study was an expected observation, as other investigators previously reported this in Tunisia and France [9, 27] . The results show that BECs are prevalent in Turkey and have been recorded in various cases of diarrheic calves, being commonly detected alone or in combination with other agents such as BRVA, BCoV and Cryptosporidium spp. The detected prevalence of BECs demonstrates their importance in calf diarrhea in Turkey, as in other countries [9, 15, 26, 27] . Therefore, BECs should be included in the routine diagnostic procedure for calf diarrhea to reveal the real causative agent(s), to determine their contributions to clinical manifestations and to take precautions to prevent calf diarrhea in Turkey.

This study not only identified GIII.2, which had been documented previously [18] [19] [20] , but was also able to identify the GIII.1 genotype for the first time ever in Turkey. Like this study, many studies in the USA, Europe, South Korea and Argentina have also detected both genotypes, although the GIII.2 genotype predominates [9, 16, 21, 32, 33] . Possible reasons for the low rate of detection of the GIII.1 genotype include the presence of GIII.P2/GIII.1 recombinant BoNoVs and the short duration of viral shedding of this genotype [2] .

Another interesting result was the detection of two potential recombinant BoNoV strains (Aksaray/K056 and Aksaray/K060). They are the first recombinant BoNoV strains described in Turkey and were identified as GIII. P1/GIII.2, similar to other recombinant strains reported in Europe [14, 15] . These recombinant strains were detected on the same farm (farm no. 2, Aksaray) sampled from diarrheic calves in 2008. Fecal samples were collected from this farm in three different years (2006, 2008 and 2016) . BoNoVs were detected in each sampling year while genotypes GIII.2 were identified in 2006 and 2016, and GIII.1 and GIII.2 were identified in 2008. These results support the argument that the presence of different genotypes concurrently is an important factor in increasing recombination [12] . As in other studies [13, 15] , these findings indicate that genotyping based only on the RdRp sequence may not be sufficient for characterizing BoNoVs due to the existence of a preferential recombination site in the ORF1-ORF2 junction region. Thus, sequence analysis of the capsid gene seems to be necessary for precise characterization [13, 15] .

Molecular and phylogenetic analyses of a portion of the RdRp region of the neboviruses in this study characterized them as Nebraska-like, similar to those reported in the USA, Italy, South Korea, Tunisia and Turkey [16, 20, 21, 27, 29] . Although there are only a few studies on the epidemiology of neboviruses, the results obtained so far indicate that Newbury1-like strains have been detected only in the United Kingdom, Brazil and Iran [28, 30, 31] .

Few studies have presented data about the capsid gene region of neboviruses. Our results show that nebovirus capsid lineage 4 is present in Turkey. AnkaraA4 is considered a potential recombinant strain due to the inconsistencies between RdRp and capsid clustering. To date, only two recombinant nebovirus strains have been reported, and their recombination sites are within the RdRp gene or at the end of the S domain [11, 16] . This study revealed that the AnkaraA4 strain has undergone a recombination event, with the breakpoint located within the S domain of the capsid protein, which is similar to another recombinant nebovirus strain, Bo/16TE/2010/ITA [16] . These results are also consistent with previous studies, suggesting that molecular analysis of the full-length capsid gene region is required to determine nebovirus lineages [10] .

In conclusion, the present study showed that both BoNoVs and neboviruses contribute to calf diarrhea, either alone or together with other enteropathogens. Genetic analysis revealed that both BoNoV GIII.1 and GIII.2 genotypes are circulating in Turkey, which allowed us to report the detection of GIII.1 in Turkey for the first time. This study also provides new genetic data for neboviruses in Turkey, which shows that they are genetically closely related to Bo/ Nebraska/80/US in the RdRp region and have a lineage 4 capsid type. Furthermore, this study is the first to report the presence of recombinant strains of BoNoVs and neboviruses in Turkey. Large-scale epidemiological studies would enhance our understanding of the genetic characteristics and evolution of BECs. Finally, considering their contribution to calf diarrhea, BECs might be incorporated into cattle vaccines in the future.

2214A) of the Scientific and Technological Research Council of Turkey (TUBITAK)

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

Infection of calves with bovine norovirus GIII.1 strain jena virus: an experimental model to study the pathogenesis of norovirus ınfection

Advances in laboratory methods for detection and typing of norovirus

Genotype 1 and genotype 2 bovine noroviruses are antigenically distinct but share a cross-reactive epitope with human noroviruses

Molecular detection of norovirus in sheep and pigs in New Zealand farms

Isolation of small viruses resembling astroviruses and caliciviruses from acute enteritis of calves

Zackenvirus"(Jena agent 117/80)-a new diarrhea pathogen in calves

Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses

Possible novel nebovirus genotype in cattle

Heterogeneity in the capsid protein of bovine enteric caliciviruses belonging to a new genus

Genomic characterization of a RdRp-recombinat nebovirus strain with a novel VP1 genotype

Genetic recombination between two genotypes of genogroup III bovine noroviruses (BoNVs) and capsid sequence diversity among BoNVs and Nebraska-like bovine enteric caliciviruses

Norovirus recombination

A chimeric bovine enteric calicivirus: evidence for genomic recombination in genogroup III of the Norovirus genus of the Caliciviridae

Genetic heterogeneity of bovine noroviruses in Italy

Evidence for recombination in neboviruses

Identification of a bovine enteric calicivirus, Kırklareli virus, distantly related to neboviruses, in calves with enteritis in Turkey

Molecular characterization of the bovine noroviruses from diarrhoeic calves in Turkey 2

First report on the phylogeny of bovine norovirus in Turkey

Detection and molecular analysis of bovine enteric norovirus and nebovirus in Turkey

Reverse transcription-PCR assays for detection of bovine enteric caliciviruses (BEC) and analysis of the genetic relationships among BEC and human caliciviruses

Detection and characterization of bovine coronaviruses in fecal specimens of adult cattle with diarrhea during the warmer seasons

AliView: a fast and lightweight alignment viewer and editor for large datasets

MUSCLE: multiple sequence alignment with high accuracy and high throughput

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

Case-control study of microbiological etiology associated with calf diarrhea

Molecular prevalence of bovine noroviruses and neboviruses detected in central-eastern Tunisia

First detection and molecular characterization of Nebovirus in Brazil

Molecular detection and characterization of unclassified bovine enteric caliciviruses in South Korea

Genomic characterization of the unclassified bovine enteric virus Newbury agent-1 (Newbury1) endorses a new genus in the family Caliciviridae

Molecular prevalence of bovine noroviruses and neboviruses in newborn calves in Iran

Molecular detection of bovine Noroviruses in Argentinean dairy calves: circulation of a tentative new genotype

Molecular epidemiology of bovine noroviruses in South Korea

Conflict of interest The authors declare that they have no conflict of interest.Ethical approval This article does not contain any studies with animals performed by any of the authors.