Distribution of Human Sapovirus Strain Genotypes over the Last Four Decades in Japan: a Global Perspective.

Sapovirus (SaV) infections are a public health problem because they cause acute gastroenteritis in humans of all ages, both sporadically and as outbreaks. However, only a limited amount of SaV sequence information, especially whole-genome sequences for all the SaV genotypes, is publicly available. Therefore, in this study, we determined the full/near-full-length genomic sequences of 138 SaVs from the 2001 to 2015 seasons in 13 prefectures across Japan. The genogroup GI was predominant (67%, n = 92), followed by genogroups GII (18%, n = 25), GIV (9%, n = 12), and GV (6%, n = 9). Within the GI genogroup, four different genotypes were identified: GI.1 (n = 44), GI.2 (n = 40), GI.3 (n = 7), and GI.5 (n = 1). We then compared these Japanese SaV sequences with 3,119 publicly available human SaV sequences collected from 49 countries over the last 46 years. The results indicated that GI.1, and GI.2 have been the predominant genotypes in Japan, as well as in other countries, over at least four decades. The 138 newly determined Japanese SaV sequences together with the currently available SaV sequences, could facilitate a better understanding of the evolutionary patterns of SaV genotypes.

Evolutionary Analysis of the VP1 and RNA-Dependent RNA Polymerase Regions of Human Norovirus GII.P17-GII.17 in 2013-2017.

Human norovirus (HuNoV) GII.P17-GII.17 (Kawasaki2014 variant) reportedly emerged in 2014 and caused gastroenteritis outbreaks worldwide. To clarify the evolution of both VP1 and RNA-dependent RNA polymerase (RdRp) regions of GII.P17-GII.17, we analyzed both global and novel Japanese strains detected during 2013-2017. Time-scaled phylogenetic trees revealed that the ancestral GII.17 VP1 region diverged around 1949, while the ancestral GII.P17 RdRp region diverged around 2010. The evolutionary rates of the VP1 and RdRp regions were estimated at ~2.7 × 10-3 and ~2.3 × 10-3 substitutions/site/year, respectively. The phylogenetic distances of the VP1 region exhibited no overlaps between intra-cluster and inter-cluster peaks in the GII.17 strains, whereas those of the RdRp region exhibited a unimodal distribution in the GII.P17 strains. Conformational epitope positions in the VP1 protein of the GII.P17-GII.17 strains were similar, although some substitutions, insertions and deletions had occurred. Strains belonging to the same cluster also harbored substitutions around the binding sites for the histo-blood group antigens of the VP1 protein. Moreover, some amino acid substitutions were estimated to be near the interface between monomers and the active site of the RdRp protein. These results suggest that the GII.P17-GII.17 virus has produced variants with the potential to alter viral antigenicity, host-binding capability, and replication property over the past 10 years.

Polio vaccination coverage and seroprevalence of poliovirus antibodies after the introduction of inactivated poliovirus vaccines for routine immunization in Japan.

In Japan, the oral poliovirus vaccine (OPV) was changed to 2 types of inactivated poliovirus vaccine (IPV), the standalone conventional IPV (cIPV) and the Sabin-derived IPV combined with diphtheria-tetanus-acellular pertussis vaccine (DTaP-sIPV), for routine immunization in 2012. We evaluated polio vaccination coverage and the seroprevalence of poliovirus antibodies using data from the National Epidemiological Surveillance of Vaccine-Preventable Diseases (NESVPD) from 2011 to 2015. Several years before the introduction of IPV in 2012, OPV administration for children was refused by some parents because of concerns about the risk of vaccine-associated paralytic poliomyelitis. Consequently, in children aged <1 years who were surveyed in 2011-2012, polio vaccination coverage (45.0-48.8%) and seropositivity rates for poliovirus (type 1: 51.7-65.9%, type 2: 48.3-53.7%, and type 3: 15.0-29.3%) were decreased compared to those surveyed in 2009. However, after IPV introduction, the vaccination coverage (95.5-100%) and seropositivity rates (type 1: 93.2-96.6%, type 2: 93.1-100%, and type 3: 88.6-93.9%) increased among children aged <1 years in 2013-2015. In particular, seropositivity rates and geometric mean titers (GMTs) for poliovirus type 3 in <5-year-old children who received 4 doses of IPV (98.5% and 247.4, respectively) were significantly higher than in those who received 2 doses of OPV (72.5% and 22.9, respectively). Furthermore, in <5-year-old children who received 4 doses of either DTaP-sIPV or cIPV, the seropositivity rates and the GMTs for all 3 types of poliovirus were similarly high (96.5-100% and 170.3-368.8, respectively). Our findings from the NESVPD demonstrate that both the vaccination coverage and seropositivity rates for polio remained high in children after IPV introduction.

Molecular Evolution of the RNA-Dependent RNA Polymerase and Capsid Genes of Human Norovirus Genotype GII.2 in Japan during 2004-2015.

The RNA-dependent RNA polymerase (RdRp) and capsid (VP1) genes of 51 GII.2 human norovirus (HuNoV) strains collected during the period of 2004-2015 in Japan were analyzed. Full-length analyses of the genes were performed using next-generation sequencing. Based on the gene sequences, we constructed the time-scale evolutionary trees by Bayesian Markov chain Monte Carlo methods. Time-scale phylogenies showed that the RdRp and VP1 genes evolved uniquely and independently. Four genotypes of GII.2 (major types: GII.P2-GII.2 and GII.P16-GII.2) were detected. A common ancestor of the GII.2 VP1 gene existed until about 1956. The evolutionary rates of the genes were high (over 10-3 substitutions/site/year). Moreover, the VP1 gene evolution may depend on the RdRp gene. Based on these results, we hypothesized that transfer of the RdRp gene accelerated the VP1 gene evolution of HuNoV genotype GII.2. Consequently, recombination between ORF1 (polymerase) and ORF2 (capsid) might promote changes of GII.2 antigenicity.

Molecular evolution of the capsid gene in human norovirus genogroup II.

Capsid protein of norovirus genogroup II (GII) plays crucial roles in host infection. Although studies on capsid gene evolution have been conducted for a few genotypes of norovirus, the molecular evolution of norovirus GII is not well understood. Here we report the molecular evolution of all GII genotypes, using various bioinformatics techniques. The time-scaled phylogenetic tree showed that the present GII strains diverged from GIV around 1630CE at a high evolutionary rate (around 10(-3) substitutions/site/year), resulting in three lineages. The GII capsid gene had large pairwise distances (maximum > 0.39). The effective population sizes of the present GII strains were large (>10(2)) for about 400 years. Positive (20) and negative (over 450) selection sites were estimated. Moreover, some linear and conformational B-cell epitopes were found in the deduced GII capsid protein. These results suggested that norovirus GII strains rapidly evolved with high divergence and adaptation to humans.

Genetic analysis of human rotavirus C: The appearance of Indian-Bangladeshi strain in Far East Asian countries.

Rotaviruses C (RVCs) circulate worldwide as an enteric pathogen in both humans and animals. Most studies of their genetic diversity focus on the VP7 and VP4 genes, but the complete genomes of 18 human RVCs have been described in independent studies. The genetic background of the Far East Asian RVCs is different than other human RVCs that were found in India and Bangladesh. Recently, a RVC detected in 2010 in South Korea had genetic background similar to the Indian-Bangladeshi RVCs. This study was undertaken to determine the whole genome of eight Japanese RVCs detected in 2005-2012, and to compare them with other human and animal global RVCs to better understand the genetic background of contemporary Far East Asian RVC. By phylogenetic analysis, the human RVCs appeared to be distinct from animal RVCs. Among human RVCs, three lineage constellations had prolonged circulation. The genetic background of the Far East Asian RVC was distinguished from Indian-Bangladeshi RVC as reported earlier. However, we found one Japanese RVC in 2012 that carried the genetic background of Indian-Bangladeshi RVC, whereas the remaining seven Japanese RVCs carried the typical genetic background of Far East Asian RVC. This is the first report of the Indian-Bangladeshi RVC in Japan. With that observation and the reassortment event of human RVCs in Hungary, our study indicates that the RVCs are spreading from one region to another.

Genetic analysis of attachment glycoprotein (G) gene in new genotype ON1 of human respiratory syncytial virus detected in Japan.

We studied the evolution of the G gene in the new genotype ON1 of RSV detected from patients with acute respiratory infection in Japan. Phylogenetic analyses and the evolutionary timescale were obtained by the Bayesian MCMC method. We also analyzed p-distance and positive selection sites. A new genotype ON1 emerged around 2001. The evolution rate was rapid (3.57 × 10(-3) substitutions/site per year). The p-distance was short and no positive selection site was found in the present strains. These results suggested that a new genotype ON1 of RSV-A emerged approximately10 years ago and spread to some countries with a high evolution rate.

Molecular epidemiological study of human rhinovirus species A, B and C from patients with acute respiratory illnesses in Japan.

Recent studies suggest that human rhinovirus species A, B and C (HRV-ABCs) may be associated with both the common cold and severe acute respiratory illnesses (ARIs) such as bronchiolitis, wheezy bronchiolitis and pneumonia. However, the state and molecular epidemiology of these viruses in Japan is not fully understood. This study detected the genomes of HRV-ABCs from Japanese patients (92 cases, 0-36 years old, mean±sd 3.5±5.0 years) with various ARIs including upper respiratory infection, bronchiolitis, wheezy bronchiolitis, croup and pneumonia between January and December 2010. HRV-ABCs were provisionally type assigned from the pairwise distances among the strains. On phylogenetic trees based on the nucleotide sequences of the VP4/VP2 coding region, HRV-A, -B and -C were provisionally assigned to 14, 2 and 12 types, respectively. The present HRV-A and -C strains had a wide genetic diversity (>30 % divergence). The interspecies distances were 0.230±0.063 (mean±sd, HRV-A), 0.218±0.048 (HRV-B) and 0.281±0.105 (HRV-C), based on nucleotide sequences, and 0.075±0.036 (HRV-A), 0.049±0.022 (HRV-B) and 0.141±0.064 (HRV-C) at the deduced amino acid level. Furthermore, HRV-A and -C were the predominant species and were detected throughout the seasons. The results suggested that HRV-A and -C strains have a wide genetic divergence and are associated with various ARIs in Japan.