key: cord-0011393-bzglwqep authors: Miyoshi, Masahiro; Goto, Akiko; Komagome, Rika; Yamaguchi, Hiroki; Maruo, Yuji; Nakanishi, Masanori; Ishida, Setsuko; Nagano, Hideki; Sugisawa, Takahisa; Okano, Motohiko title: Genetic characterization of a novel recombinant echovirus 30 strain causing a regional epidemic of aseptic meningitis in Hokkaido, Japan, 2017 date: 2019-12-11 journal: Arch Virol DOI: 10.1007/s00705-019-04484-7 sha: 336bac252b0759ff776fadae90d4c2cfe969bee0 doc_id: 11393 cord_uid: bzglwqep A regional epidemic of aseptic meningitis caused by echovirus 30 (E30) occurred in Hokkaido, Japan, during the period of August-December 2017. To investigate their phylogenetic relationship to other human enteroviruses, we determined the complete genomic nucleotide sequences of isolates from this outbreak. Phylogenetic analysis of the viral capsid protein 1 gene showed that the strains were most closely related to E30 strains detected in Germany, France, and Russia in 2013. In contrast, the region encoding the viral protease and the RNA-dependent RNA polymerase had a close phylogenetic relationship to non-E30 enteroviruses detected in the United Kingdom and Switzerland in 2015-2017, suggesting that a recombination event had occurred. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-019-04484-7) contains supplementary material, which is available to authorized users. Echovirus 30 (E30) is one of the serotypes belonging to the species Enterovirus B [1] and is a major pathogen causing aseptic meningitis (AM) in both children and adults. Large epidemics caused by E30 have been frequently reported globally [2] . We experienced a regional epidemic of AM caused by E30 in Hokkaido, Japan, during the period of August-December 2017. Details of the epidemic and the clinical features were presented previously [3] . In this study, we determined the complete nucleotide sequences of viruses from the outbreak to determine their molecular features and their phylogenetic relationship to other strains of enterovirus B (EV-B). Virus isolation, RNA extraction, and sequencing were performed as described earlier [4] . Genome amplification and the nucleotide sequencing were carried out using specific primers (Supplementary Table S1 ). We determined the complete nucleotide sequences of four isolates, named E30/Hokkaido.JPN/21208 (Hokkaido21208), 21214, 21246, and 21326, that were collected on August 31, September 2 and 18, and November 24, 2017 , and submitted them to the GenBank database under the accession numbers LC416533 to LC416536, respectively. To investigate the genetic relationships of these isolates to globally circulating strains, phylogenetic analysis based on the viral capsid protein 1 (VP1) gene, the viral protease gene, and the RNA-dependent RNA polymerase (3CD) gene were carried out. Phylogenetic analysis based on the VP1 gene (876 bases) was performed in comparison with 248 E30 strains isolated in 33 countries during 1957-2017, and analysis based on the 3CD gene (1938 bases) was performed with 243 EV-B strains isolated in 33 countries during 1948-2017 (Supplementary Tables S2 and S3) . With respect to choosing the strains, we took into account the genetic classifications of the strains that had been evaluated [2, [4] [5] [6] [7] . Phylogenetic trees were constructed using Molecular Evolutionary Genetic Analysis software (MEGA v.5) [8], using 1 3 the maximum-likelihood method based on the general timereversible model. The genetic lineages were confirmed by bootstrap values using 1000 pseudoreplicates [9] . Furthermore, the nucleotide sequences were aligned with those of 243 strains (Supplementary Table S3 ) of EV-B using MEGA v.5. The gap-opening and extension penalties were set to 15 and 6.66, respectively, and the alignments were corrected manually to match the open reading frame (ORF). In order to identify potential recombination sites in the genome, a similarity plot analysis was performed using Simplot software version 3.5.1 [10] with a sliding window of 400 nucleotides moving in steps of 50 nucleotides. The genomes of four E30 strains from Hokkaido consisted of 7428 nucleotides including a 5' untranslated region (UTR) of 742 nucleotides followed by an ORF encoding a polypeptide of 2194 amino acids and a 3' UTR of 101 nucleotides. The nucleotide and amino acid sequences of these strains were highly conserved (99.6-9% and 99.8-9% identity), suggesting that a single lineage of E30 had spread. The sequence of the 5' UTR and the P1 genomic region of the strain Hokkaido21208 showed the highest identity to those of the strain 13-759, which was isolated in Germany in 2013 [11] (Table 1) . While a high degree of nucleotide sequence similarity was also found in the 2A and the 2B genes, the sequence identity values from the 2C gene to the 3' UTR were considerably lower (ranging from 71.2% to 84.0%). However, despite this divergence, the amino acid identity values from the P1 region to the P3 region ranged from 90.9% to 99.6% (Table 1) . The result of similarity plot analysis correlated with the above results ( Fig. 1) . From the 5' UTR to the 2B gene, the strain Hokkaido21208 was most closely related to the strain 13-759. However, the similarity of the 2C gene dropped abruptly to the same level as with other serotypes of EV-B. By contrast, the genomic region between the 2C gene and the 3' UTR of the strain Hokkaido21208 did not show significant similarity to any of the strains compared, suggesting that a recombination event had occurred at the border between the 2B and 2C genes. This breakpoint has been previously identified as a recombination hot spot [12] . A previous study indicated that the VP1 genes of E30 detected in Japan before 1998 were divided into six genogroups [5] . These genogroups correspond to the lineages E30_d, e, and g based on an advanced report by Bailly et al. [6] . In this study, our analysis indicated that the Hokkaido strains belong to the lineage E30_f ( Fig. 2A) , forming a novel cluster with Japanese strains detected in Fukushima prefecture, and with German, French, and Russian strains detected during 2008-2014. In contrast, several Japanese strains detected during 2008-2014 in Toyama, Osaka, Oita, and Okinawa prefectures, located mainly in western Japan, were classified in lineage E30_b ( Fig. 2A) . These data suggest that the Hokkaido strains might have originated overseas, in particular, Europe. E30 strains have been detected almost every year in Japan [13] . However, due to insufficient genetic information of these isolates, the detailed phylogenetic and topological relationships among Japanese strains and the time point of importation of lineage E30_f could not be clarified. Amino acid substitutions in VP1 of the Hokkaido strains were found in at least 20 positions in comparison with the strain 25-0284/osaka.AB908259/JPN/2013, which belongs to lineage E30_b (Table 2) . Among these sites, 13 substitutions (56Y/F, 87E/D, 120I/V, 145V/I, 156R/K, 157S/G, 234N/S, 247K/R, 269V/I, 274E/D, 284S/T, 285T/N, and 289L/M) showed divergence from these lineages. Codons 87 and 156, which comprise a portion of the BC loop and the EF loop, respectively, have been predicted not only to be exposed on the surface but also to represent putative positive selection sites [6, 14] . Besides, several substitutions are in the C-terminal region of VP1 (amino acids 260-292), which has been recognized as one of the main antigenic regions of enteroviruses [15] [16] [17] . The emergence of an epidemic would implicitly indicate the lack of herd immunity in the community. Furthermore, the novel strains, which have different antigenic properties, even within the same serotype, could have escaped the host's immune system. Earlier studies have demonstrated that E30 is antigenically heterogeneous and divided into three major groups [18, 19] . The VP1 nucleotide sequences of these groups are classified into the lineage E30_a or c, while the differences in antigenic properties do not correlate with the comprehensive classification [6] . Under these circumstances, it is important to explore whether the strains are pathogenic and able to spread beyond a regional scale in the near future. However, the strain-specific features that might be involved in disease are still not well understood [11, 20] , and more-detailed study is needed to evaluate the characteristics of each lineage and to understand the potential endemic risk and the dynamics of virus transmission. The 3CD genes of EV-B are divided phylogenetically into three groups (1, 2, and 3) (Fig. 2B) [4, 7] . Our findings indicate that the Hokkaido strains belong to group 3 and form a cluster with coxsackievirus type B2 (CVB2), CVB3, CVB4, CVB5, E3, E6, E11, E13, enteroviruses 85, and B106 strains detected in European and Asian countries during the last 20 years. Among these strains, CVB4, CVB5, and E3 strains detected in the United Kingdom and Switzerland during 2015-2017 occupied the nearest branch (Fig. 2B) . No other E30 strain except for the Hokkaido strains was found in this cluster. In contrast, the E30 strain 13-759, whose VP1 gene showed a high degree of similarity to those of the Hokkaido strains, was classified into group 2. Group 2 included several EV-B strains detected intermittently in Europe, suggesting that this genome has been maintained there since the previous century (Fig. 2B) [6]. In addition to the data for the VP1 genes described above, chronological and phylogenetic information about the 3CD genes strongly suggest that the Hokkaido E30 strains originated in Europe. The emergence of E30 has been found to exhibit a cyclic pattern with dominance and disappearance of individual recombinants over a period of years [21] [22] [23] [24] [25] [26] [27] . The turnover period for such recombinants has been calculated to be from 1.8 to 5 years [28, 29] . Therefore, long-term surveillance based on detection of virus in clinical specimens and environmental materials would contribute to our understanding trends and epidemiological patterns [30] [31] [32] [33] [34] . Furthermore, investigation of the antigenic properties of circulating viruses is fundamental to understand their pathogenic features and, ultimately, to promote measures that are important for preventive public health. Lineage E30_h Lineage CF1590401.AM237326/FRA/2005 LP329 CVB2_NS13-006-Nagasaki13 CVB2_NS13-010-Nagasaki13.LC191416/JPN/2013 CVB2_08-824-Shimane08 Changes to taxonomy and the international code of virus classification and nomenclature ratified by the international committee on taxonomy of viruses 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 Strain-dependent effects of clinical echovirus 30 outbreak isolates at the blood-CSF barrier Recombination in circulating Human enterovirus B: independent evolution of structural and non-structural genome regions National Institute of Infectious Diseases and Tuberculosis and Infectious Diseases Control Division, Ministry of Health The structure of coxsackievirus B3 at 3.5 Å resolution Phylogenetic patterns of human coxsackievirus B5 arise from population dynamics between two genogroups and reveal evolutionary factors of molecular adaptation and transmission Molecular epidemiology of echovirus 30: temporal circulation and prevalence of single lineages Antigenic sites of coxackievirus A9 A comparative study of 63 strains of ECHO virus type 30 The antigenic heterogeneity of type 30 echoviruses Transmission of human enterovirus 85 recombinants containing new unknown serotype HEV-B donor sequences in Xinjiang Uighur Autonomous Region, China Molecular epidemiology and prevalence of echovirus 30 in Zhejiang Province Molecular epidemiological study of HEV-B enteroviruses involved in the increase in meningitis cases occurred in Spain during Molecular characterization of two novel echovirus 18 recombinants associated with handfoot-mouth disease Recombination analysis of coxsackievirus B5 genogroup C Echovirus plays major roles in the natural recombination of Coxsackievirus B3 Intraand inter-serotypic recombinations in the 5'UTR-VP4 region of Echovirus 30 strains Echovirus plays a major role in natural recombination in the coxsackievirus B group Transmission networks and population turnover of echovirus 30 Widespread circulation of a new echovirus 30 variant causing aseptic meningitis and non-specific viral illness A large Finnish echovirus 30 outbreak was preceded by silent circulation of the same genotype Molecular epidemiology of echoviruses 11 and 13, based on an environmental surveillance conducted in Toyama Prefecture Widespread circulation of echovirus type 13 demonstrated by increased seroprevalence in Toyama Endemic transmission of echovirus 30 in Toyama, Japan in 2010 is verified by environmental surveillance Assessment of an enterovirus sewage surveillance system by comparison of clinical isolates with sewage isolates from The authors are grateful to the staffs of local hospitals, the Kushiro Center of Public Health, and the Department of Health and Welfare, Hokkaido Government, for collecting samples and providing the epidemiological information.Funding This work was financially supported by a grant-in-aid for scientific research from the Hokkaido Government, Japan. The authors declare no conflict of interest.