key: cord-0003005-2htnywef authors: Tsuchiaka, Shinobu; Naoi, Yuki; Imai, Ryo; Masuda, Tsuneyuki; Ito, Mika; Akagami, Masataka; Ouchi, Yoshinao; Ishii, Kazuo; Sakaguchi, Shoichi; Omatsu, Tsutomu; Katayama, Yukie; Oba, Mami; Shirai, Junsuke; Satani, Yuki; Takashima, Yasuhiro; Taniguchi, Yuji; Takasu, Masaki; Madarame, Hiroo; Sunaga, Fujiko; Aoki, Hiroshi; Makino, Shinji; Mizutani, Tetsuya; Nagai, Makoto title: Genetic diversity and recombination of enterovirus G strains in Japanese pigs: High prevalence of strains carrying a papain-like cysteine protease sequence in the enterovirus G population date: 2018-01-11 journal: PLoS One DOI: 10.1371/journal.pone.0190819 sha: eceb1667ac6722e9131eb7d0f560a7f76bfb25f1 doc_id: 3005 cord_uid: 2htnywef To study the genetic diversity of enterovirus G (EV-G) among Japanese pigs, metagenomics sequencing was performed on fecal samples from pigs with or without diarrhea, collected between 2014 and 2016. Fifty-nine EV-G sequences, which were >5,000 nucleotides long, were obtained. By complete VP1 sequence analysis, Japanese EV-G isolates were classified into G1 (17 strains), G2 (four strains), G3 (22 strains), G4 (two strains), G6 (two strains), G9 (six strains), G10 (five strains), and a new genotype (one strain). Remarkably, 16 G1 and one G2 strain identified in diarrheic (23.5%; four strains) or normal (76.5%; 13 strains) fecal samples possessed a papain-like cysteine protease (PL-CP) sequence, which was recently found in the USA and Belgium in the EV-G genome, at the 2C–3A junction site. This paper presents the first report of the high prevalence of viruses carrying PL-CP in the EV-G population. Furthermore, possible inter- and intragenotype recombination events were found among EV-G strains, including G1-PL-CP strains. Our findings may advance the understanding of the molecular epidemiology and genetic evolution of EV-Gs. a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 Next-generation sequencing was conducted on cDNA libraries constructed from total RNA of 222 fecal samples. Via a BLAST search, 59 EV-G-like contigs that were longer than 5,000 nt, including the entire VP1-coding sequence with more than threefold coverage of sequence reads, were identified in 35 (15.8%) normal fecal samples, five (2.3%) fecal samples for mild diarrhea, and 10 (4.5%) diarrheic fecal samples (Tables 1 and 2 ). Eight of 50 samples revealed more than two EV-G-like contigs (Table 1) . Apart from three samples (Ishi-Ya3, Ishi-Ka3, and Ishi-Ka7), 47 samples were found to contain other viruses: Rotaviruses A, C, or H; orthoreovirus; kobuvirus; picobirnavirus; astrovirus; porcine epidemic diarrhea virus (PEDV); posavirus; sapelovirus; St-Valerien virus; sapovirus; or teschovirus (Table 1 and S3 Table) . Because EV-G discrimination is based on sequence identities of the complete VP1 gene (http://www.picornaviridae.com/), a phylogenetic tree using the complete VP1 nt sequence was constructed. Japanese EV-G strains clustered with the reference strains of G1 (17 strains), G2 (four strains), G3 (22 strains), G4 (two strains), G6 (two strains), G9 (six strains), and G10 (five strains). One strain named Ishi-Ka2 branched independently and did not cluster with any reference strains (Fig 1) . Because >25% difference in complete VP1 nt sequences between isolates is a criterion for genotype classification [13] [14] , a pairwise comparison of complete VP1 nt sequences was conducted (Table 3 and S1 Table) . Bu6-5, Bu8-2, Bu8-4, and Ishi-Ya4-2 formed one cluster but were found to be slightly related to the G3 group. Although complete VP1 nt sequence identities of these strains to those of other G3 strains (except for Ishi-Ka3, Ishi-Ka3-1, Ishi-Ka4, Ishi-Ka5-1, Ishi-Ka6, and Ishi-Ka7) were <75.0% (69.8%-74.9%), these four strains shared !75.0% nt sequence identities with Ishi-Ka3, Ishi-Ka3-1, Ishi-Ka4, Ishi-Ka5-1, and Ishi-Ka6 (S1 Table) . Therefore, we tentatively classified these strains as G3-lineage 2 (Fig 1) . Ishi-Ka2 revealed low nt sequence identities (57.5% to 73.1%) with other genotypes and thus Ishi-Ka2 may represent a new serotype of EV-G (Table 3) . Complete or nearly complete aa sequences of the coding-sequence region (CDS) of 59 Japanese EV-G strains were aligned and compared. We found that 17 of the 59 strains contain extra 633 to 651 nt (211 to 217 aa) within the 2C-3A coding region. According to BLAST analysis, these sequences have sequence homology to the PL-CP sequence variants that were recently identified in EV-G strains from the USA and Belgium [24] [25] [26] . Each inserted sequence is located between the coding regions 2C and 3A as in the USA and Belgium strains. The insertion sequences were aligned and compared with those of the PL-CP of EV-G strains from the USA and Belgium and with the PL-VP sequences in the genome of nidoviruses including porcine and bovine toroviruses by phylogenetic analysis and pairwise sequence comparison (Fig 2 and S2 Table) . The PL-CP sequence of Japanese EV-G1 and that of Japanese EV-G2 revealed !74.0% nt and !74.6% aa sequence identities to each other and to the USA and Belgium EV-Gs and clustered in one group but are distantly related to those of porcine and bovine toroviruses, showing lower sequence identities (57.0% to 64.6% in the nt sequence and 49.6% to 58.7% in the aa sequence). EV-G HgYa2-1 and porcine torovirus HgYa2-2 were identified on the same farm at the same time; however, the nt and aa sequence identity between the PL-CP sequences of those strains was 62.3% and 54.3%, respectively (S2 Table) . Japanese EV-G strains carrying PL-CP were subdivided into G1-PL-CP lineage 1, G1-PL-CP lineage 2, and G2 in the VP1 phylogenetic tree (Fig 1) ; however, these groups were not clearly detectable in the PL-CP tree (Fig 2) To further investigate the genomic relations among EV-G strains, phylogenetic trees based on nt sequences of three regions (VP4-VP3, VP1, P2, and P3) were constructed. The tree for VP4-VP3 was similar to that of VP1, but the P2 and P3 trees showed topologies different from each other and no clear-cut EV-G types could be defined (Fig 3A) . EV-G1-PL-CP lineage 3 was found to be related to G1-PL-CP lineage 2 and 3 in the trees for VP4-VP3 and VP1, whereas G1-PL-CP lineage 3 was closely related to G3 and G9 strains in the P2 tree and to G3 and G1-PL-CP strains in the P3 tree. The G2-PL-CP strain HgYa2-1 showed high homology Genetic diversity and recombination of enterovirus G strains in Japanese pigs to G2 strains in the tree for VP4-VP3 and VP1, whereas HgYa2-1 showed high similarity with G1-PL-CP lineage 1 strains in regions P2 and P3 (Fig 3A) . By SimPlot analysis, the crossover site was mapped to the 2A region. The G2-PL-CP strain HgYa2-1 revealed that the downstream region of the crossover site has high similarity to the G1-PL-CP strain MoI2-2-1 ( Fig 3C) . To find a possible recombination breakpoint, a bootstrap scanning analysis was conducted. A possible recombination breakpoint was identified in the middle of the 2A region (Fig 3D) . Ishi-Ka2 branched independently in the trees for VP4-VP3 and VP1, whereas Ishi-Ka2 clustered with G3-lineage 1 and formed a cluster with G3 strains identified on the same farm (S1A Fig). SimPlot analysis suggested that Ishi-Ka2 has extremely high similarity to G3-lineage 1 strain Ishi-Ka7 in regions 2C and P3 (S1C Fig). G3-lineage 2 strains showed a topology similar to that of the VP1 tree in VP4-VP3; however, the G3-lineage 2 strains branched separately from the G3-lineage-1 strains in the P2 and P3 trees (S1A Fig) and were found to be closely related to each other throughout the genome (S1D Fig) . Although we did not initially aim to determine EV-G prevalence among pigs in Japan in this study, contigs that were longer than 5,000 nt were found in 22.5% (50/222) of pigs and on 23.4% (18/77) of farms, suggesting that EV-Gs are widespread among Japanese pigs. Fortyfour strains out of 59 (74.6%) were detected in healthy pigs, indicating that EV-Gs seem not to be associated with diarrhea in pigs, in accordance with other reports [13] [14] 38] . Because the detection limit of the method was not tested, a true prevalence study is needed in the future. EV-G genotyping is based on >25% divergence between VP1 nucleotide sequences [14, 39] . In the present study, according to the criteria, seven genotypes (G1, G2, G4, G6, G9, G10, and G?) were found in the feces samples from Japanese pigs, and the predominant genotypes were G3 (37.3%; 22/59) and G1 (28.8%; 17/59; Table 1 , Fig 1) . There are few studies on the genotyping of EV-Gs in pigs, and limited information is available from DDBJ/EMBL/GenBank databases. G1 and G6 types are predominant genotypes in Vietnam, whereas G3, G2, and G4 types appear to be common genotypes in Spain (however, that study did not analyze complete VP1 sequence) [40] . To date, G1-G16 genotypes and at least three EV-Gs with an unassigned genotype, including Ishi-Ka2 in this study, have been reported [13-14, 25, 41] . Owing to the limited number of reports on a specific geographic area, probably not all genotypes of EV-Gs are known at present. Further studies are needed for a comprehensive understanding of the genetic diversity of EV-Gs in other geographic areas. Picornaviruses show significant genetic variability driven by both mutations and recombination events [42] [43] . Ishi-Ka2 manifested >25% VP1 nucleotide sequence divergence from other strains; therefore, Ishi-Ka2 can be considered a new serotype of EV-Gs. Nonetheless, Ishi-Ka2 shares high sequence homology with the G3-lineage 1 strain Ishi-Ka7, which was identified in a pig kept on the same farm, except for the P1 region. It is likely that Ishi-Ka2 emerged by possible recombination events; however, the putative recombination points could not be identified, and the origins of these recombination events are unclear because the recombination counterparts of these strains could not be found in the DDBJ/EMBL/GenBank databases or our dataset. G3-lineage 2 strains have sequence homology to G3-lineage 1 in the P1 region, but they are distantly related to G3-lineage 1 on the basis of regions P2 and P3. Because VP1 induces a host immune response, serological properties can be hypothesized based on sequence homology of the VP1 gene. On the other hand, these results suggest that full genome analysis may be needed in addition to the genotyping approaches based solely on the VP1 gene for precise EV-G classification. RNA recombination events contribute to genetic diversity and may lead to changes in virulence, escape from host immunity, and adaptation to a new host [42] [43] [44] [45] [46] [47] [48] [49] [50] . EV-G strains carrying PL-CP in pigs with diarrhea have been reported in the USA and Belgium [24] [25] [26] . In these cases, EV-G-PL-CPs were detected solely or with low abundance of PEDV. Shang et al. constructed an infectious clone of the EV-G-PL-CP strain, 08/NC_USA/2015, and compared it with a PL-CP knockout recombinant virus. They found that the PL-CP knockout virus showed impaired growth and induced higher expression levels of innate-immunity genes, suggesting that EV-G-PL-CP strains acquire pathogenesis via a recombination event [25] . Four out of 17 Japanese EV-G-PL-CP strains were detected in diarrheic cases of pigs; however, 13 EV-G-PL-CP strains were isolated from healthy pigs. In all cases of detection of EV-G-PL-CP in Japan, EV-G-PL-CP strains were identified together with other enteric viruses, such as astrovirus, sapelovirus, posavirus, rotavirus, picobirnavirus, sapovirus, teschovirus, torovirus, PEDV, St-Valerien virus, or kobuvirus (Table 1) . Mixed infection with EV-G-PL-CP and other enteric viruses may influence the pathogenicity of EV-G-PL-CP strains. The sequences of PL-CP of Japanese EV-G PL-CP strains are distantly related to the sequences derived from ORF1 of toroviruses, even though they were simultaneously identified on the same farm, and they have homology to those of USA and Belgium strains (Fig 2) , suggesting that a recombination event between an EV-G and torovirus occurred in the past. By recombination analyses, possible recombination events between EV-G-PL-CP strains were uncovered and a recombination breakpoint was identified in the middle of region 2A (Fig 3) , in agreement with another report that describes a recombinant event between EV-G8 and EV-G9 [14] , suggesting that this point may be a hotspot of recombination events of EV-G. Furthermore, VP1-2A junction is a known recombination hot-spot in human enteroviruses and this was discussed in many papers [51] [52] [53] [54] [55] . The present recombination profile in EV-G described here apparently mirrors that in human enteroviruses. EV-Gs that received PL-CP have been evolving independently and gaining genetic diversity via recombination events. By a metagenomics approach, high genetic diversity of EV-Gs, including new genotypes and high prevalence of EV-Gs carrying PL-CP, was observed among EV-G isolates from the feces of Japanese pigs. EV-Gs comingle and pose a risk of coinfection in the current growing and high-density pig husbandry system of Japan. Coinfection of a single animal with multiple EV-Gs, including EV-G-PL-CP strains, may lead to recombination events, which may in turn promote genetic diversity of EV-Gs and EV-G-PL-CPs. These findings may improve our understanding of the molecular epidemiology and evolution of EV-Gs. Porcine enteric picornaviruses Classification of porcine enteroviruses by antigenic analysis and cytopathic effects in tissue culture: description of 3 new serotypes Grouping of porcine enteroviruses by indirect immunofluorescence and description of two new serotypes Genetic reclassification of porcine enteroviruses Sequencing of porcine enterovirus groups II and III reveals unique features of both virus groups Porcine teschoviruses comprise at least eleven distinct serotypes: molecular and evolutionary aspects The association of group III porcine enteroviruses with epithelial tissue Characterization of a novel porcine enterovirus in domestic pig in Hungary Characterization of a novel porcine enterovirus in wild boars in Hungary Natural interspecies recombinant bovine/porcine enterovirus in sheep First identification and characterization of porcine enterovirus G in the United States Complete genome sequence of a novel porcine enterovirus strain in China Prevalence, genetic diversity and recombination of species G enteroviruses infecting pigs in Vietnam Large-scale screening and characterization of enteroviruses and kobuviruses infecting pigs in Vietnam Characterization of a genetically heterogeneous porcine rotavirus C, and other viruses present in the fecal virome of a non-diarrheic Belgian piglet Development of Polioencephalomyelitis in Cesarean-Derived Colostrum-Deprived Pigs Following Experimental Inoculation with Either Teschovirus A Serotype 2 or Serotype 11 Detection of a novel sapelovirus in central nervous tissue of pigs with polioencephalomyelitis in the USA Experimental teschovirus encephalomyelitis in gnotobiotic pigs Multiple models of porcine teschovirus pathogenesis in endemically infected pigs The role of porcine teschovirus in causing diseases in endemically infected pigs Phylogeny and evolution of porcine teschovirus 8 isolated from pigs in China with reproductive failure Porcine teschovirus polioencephalomyelitis in western Canada Prevalence of porcine enterovirus 9 in pigs in middle and eastern China A porcine enterovirus G associated with enteric disease contains a novel papain-like cysteine protease A Naturally Occurring Recombinant Enterovirus Expresses a Torovirus Deubiquitinase Identification of an enterovirus recombinant with a torovirus-like gene insertion during a diarrhea outbreak in fattening pigs The papain like protease of severe acute respiratory syndrome coronavirus has deubiquitinating activity The papain-like protease from the severe acute respiratory syndrome coronavirus is a deubiquitinating enzyme Nidovirus papain-like proteases: multifunctional enzymes with protease, deubiquitinating and deISGylating activities Emerging role of ISG15 in antiviral immunity Interferon-induced ISG15 pathway: an ongoing virus-host battle Full genome analysis of bovine astrovirus from fecal samples of cattle in Japan: identification of possible interspecies transmission of bovine astrovirus The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods Confidence limits on phylogenies: an approach using the bootstrap Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination RDP4: Detection and analysis of recombination patterns in virus genomes The survey of porcine teschoviruses, sapeloviruses and enteroviruses B infecting domestic pigs and wild boars in the Czech Republic between Picornavirus and enterovirus diversity with associated human diseases Distribution and genetic characterization of Enterovirus G and Sapelovirus A in six Spanish swine herds A novel enterovirus species identified from severe diarrheal goats RNA structural elements determine frequency and sites of nonhomologous recombination in an animal plus-strand RNA virus Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses Frequency and dynamics of recombination within different species of human enteroviruses Role of recombination in evolution of enteroviruses Recombination among picornaviruses Mechanisms of genetic variation in polioviruses Recombination in viruses: mechanisms, methods of study, and evolutionary consequences Why do RNA viruses recombine? Nonhomologous recombination between defective poliovirus and coxsackievirus genomes suggests a new model of genetic plasticity for picornaviruses Recombination in circulating enteroviruses Recombination in uveitis-causing enterovirus strains Recombination in circulating Human enterovirus B: independent evolution of structural and non-structural genome regions Analysis of echovirus 30 isolates from Russia and new independent states revealing frequent recombination and reemergence of ancient lineages Recombination strategies and evolutionary dynamics of the Human enterovirus A global gene pool Mizutani. We received three funding including JSPS KAKENHI 15K07718 to Makoto Nagai, Global Innovation Research of Tokyo University of Agriculture and Technology to Tetsuya Mizutani, and the Research Project for Improving Food Safety and Animal Health of the Ministry of Agriculture, Forestry and Fisheries of Japan. We are very grateful to Ms. Kanako Satoh and Ms. Megumi Mori for technical assistance. Conceptualization: Tetsuya Mizutani, Makoto Nagai.