Phylogenetic and cluster analysis of human rhinovirus species A (HRV-A) isolated from children with acute respiratory infections in Yamagata, Japan

Longitudinal Epidemiology of Viral Infectious Diseases Combining Virus Isolation, Antigenic Analysis, and Phylogenetic Analysis as Well as Seroepidemiology in Yamagata, Japan, between 1999 and 2018

We introduced a microplate method for virus isolation in the Department of Microbiology, Yamagata Prefectural Institute of Public Health (YPIPH) in 1999 in Yamagata, Japan. We have since carried out longitudinal epidemiological studies on viral infectious diseases, particularly respiratory viruses, combining traditional technologies such as virus isolation and serological techniques and newly developed molecular methods. Here, we provide an overview of our activities at YPIPH between 1999 and 2018. During the study period, we observed emerging and re-merging diseases such as those caused by echovirus type 13, enterovirus D68, parechovirus-A3 (PeV-A3), and Saffold virus. With regard to PeV-A3, we proposed a new disease concept, "PeV-A3-associated myalgia/myositis." We also revealed the longitudinal epidemiologies of several viruses such as enterovirus A71 and coxsackievirus A16. To perform longitudinal epidemiological studies at any time in Yamagata, we established a system for stocking clinical specimens, viral isolates, complementary DNAs, and serum specimens. We have also pursued collaboration works with virology laboratories across Japan. We hope our experiences, findings, and research materials will further contribute to the development of countermeasures against viral infectious diseases and improvement in public health strategies in Yamagata, Japan, Asia, and around the world.

MIP-1α level in nasopharyngeal aspirates at the first wheezing episode predicts recurrent wheezing

BACKGROUND

Respiratory virus-induced wheezing, such as that induced by respiratory syncytial virus (RSV) and human rhinovirus, is an important risk factor for recurrent wheezing and childhood asthma. However, no biomarkers for predicting recurrent wheezing have been identified.

OBJECTIVE

We searched for predictors of recurrent wheezing using nasopharyngeal aspirates obtained from patients during the first wheezing episode who were hospitalized with an acute lower respiratory tract illness.

METHODS

We enrolled 82 infants during the first wheezing episode (median age, 5.0 months) who were hospitalized for acute lower respiratory tract illness between August 2009 and June 2012 and followed these patients for 2.5 years. Nasopharyngeal aspirates and blood samples were obtained on the first day of hospitalization. Viral genomes were identified by using RT-PCR and sequencing. Levels of 33 cytokines, tryptase, IgE, anti-RSV IgE, and anti-RSV IgG were measured by using ELISAs or the Bio-Plex multiplex assay. Predictors of recurrent wheezing were examined by using a stepwise logistic regression model with backward elimination.

RESULTS

Sixty percent of the patients experienced recurrent wheezing episodes. One or more viruses were detected in the nasopharynxes of 93% of the patients during the first wheezing episode. IFN-γ, IL-2, IL-9, MIP-1α, and MIP-1β levels were significantly higher among patients with recurrent wheezing than among those without recurrent wheezing (P < .05 or .01). The stepwise model demonstrated that the MIP-1α level (odds ratio, 7.72; 95% CI, 1.50-39.77; P = .015) was the strongest independent predictor of the occurrence of recurrent wheezing.

CONCLUSION

An increased MIP-1α level in nasopharyngeal aspirates from patients with acute respiratory symptoms during the first wheezing episode caused by viral infections might predict recurrent wheezing.

Molecular evolution of VP3, VP1, 3C(pro) and 3D(pol) coding regions in coxsackievirus group A type 24 variant isolates from acute hemorrhagic conjunctivitis in 2011 in Okinawa, Japan

A large acute hemorrhagic conjunctivitis (AHC) outbreak occurred in 2011 in Okinawa Prefecture in Japan. Ten strains of coxsackievirus group A type 24 variant (CA24v) were isolated from patients with AHC and full sequence analysis of the VP3, VP1, 3C(pro) and 3D(pol) coding regions performed. To assess time-scale evolution, phylogenetic analysis was performed using the Bayesian Markov chain Monte Carlo method. In addition, similarity plots were constructed and pairwise distance (p-distance) and positive pressure analyses performed. A phylogenetic tree based on the VP1 coding region showed that the present strains belong to genotype 4 (G4). In addition, the present strains could have divided in about 2010 from the same lineages detected in other countries such as China, India and Australia. The mean rates of molecular evolution of four coding regions were estimated at about 6.15 to 7.86 × 10(-3) substitutions/site/year. Similarity plot analyses suggested that nucleotide similarities between the present strains and a prototype strain (EH24/70 strain) were 0.77-0.94. The p-distance of the present strains was relatively short (<0.01). Only one positive selected site (L25H) was identified in the VP1 protein. These findings suggest that the present CA24v strains causing AHC are genetically related to other AHC strains with rapid evolution and emerged in around 2010.

Sequencing human rhinoviruses: direct sequencing versus plasmid cloning

Human rhinoviruses (RV) are associated with the majority of viral respiratory illnesses in infants, children and adults. Over the last several years, researchers have begun to sequence the many different species and strains of RV in order to determine if certain species were associated with increased disease severity. There are a variety of techniques employed to prepare samples for sequencing. One method utilizes plasmid-cloning, which is expensive and takes several hours to complete. Recently, some investigators have instead used direct sequencing to sequence RV strains, allowing for omission of the time- and labor-intensive cloning step. This study formally compares and contrasts the sequencing results obtained from plasmid-cloning and direct Sanger sequencing of a 500 base pair PCR product covering the VP4/VP2 region of RV. A slightly longer sequence (by 65 base pairs on average) was obtained when specimens were plasmid-cloned, and the sequences were 86% similar. After trimming the extra base pairs from the cloned sequences, the sequences were 99.7% identical. Overall success of directly sequencing samples was similar to that of cloning, 5% on average failed for each technique. Therefore, in many instances, directly sequencing samples may be considered in lieu of the more expensive and time-consuming plasmid-cloning technique.

Severity of respiratory signs and symptoms and virus profiles in Japanese children with acute respiratory illness

Associations between the severity of respiratory signs and symptoms and the respiratory viruses identified in 214 Japanese children with acute respiratory illness (ARI) enrolled between January and December 2012 were studied. Respiratory rate, wheezing, cyanosis, and the use of accessory muscles were used as indices of respiratory severity and phylogenetic analysis of the viruses identified in these children was performed. Respiratory viruses such as respiratory syncytial virus (RSV), human rhinovirus (HRV), human parainfluenza virus (HPIV), and human metapneumovirus (HMPV) were prevalent, being detected in approximately 70% of the patients (151/214 patients). Co-detection of viruses occurred in about 9% of patients. RSV was identified more frequently in cases scored as moderate/severe than in those scored as mild (P < 0.05). Severity scores of patients with RSV were significantly higher than those of cases with HPIV. Moreover, severity scores in patients with mild disease and co-detections were higher than in those in whom only HPIV or adenovirus was detected. Phylogenetic analysis showed that many genotypes of HRV-A and -C with wide genetic divergence were associated with acute respiratory illness (ARI). On the other hand, only a limited number of genotypes of RSV were associated with ARI. HPIV and HMPV were associated with ARI at similar frequencies. These results suggest that different respiratory viruses with unique genetic characteristics can be found in patients with mild to severe ARI.

Molecular evolution of the haemagglutinin-neuraminidase gene in human parainfluenza virus type 3 isolates from children with acute respiratory illness in Yamagata prefecture, Japan

We conducted detailed genetic analyses of the haemagglutinin-neuraminidase (HN) gene in 272 human parainfluenza virus type 3 (HPIV3) isolates from children with acute respiratory illness during the period 2002-2009 in Yamagata prefecture, Japan. A phylogenetic tree reconstructed by the Bayesian Markov chain Monte Carlo method showed that the strains diversified at around 1946 and that the rate of molecular evolution was 1.10×10(-3) substitutions per site per year. Identity was high among the present strains (<90 %) and the pairwise-distances were short. Furthermore, we found four positive selection sites and some key amino acid substitutions in active/catalytic sites of the HN protein. The results suggest that the HN gene of HPIV3 in the present strains evolved rapidly, similarly to other virus genes such as the G gene of respiratory syncytial virus. However, the biological functions and detailed structures of the HN glycoprotein in some of these strains may have been altered.

Genetic analysis of human rhinovirus species A to C detected in patients with acute respiratory infection in Kumamoto prefecture, Japan 2011-2012

We performed detailed genetic analysis of the VP4/VP2 coding region in human rhinovirus species A to C (HRV-ABC) strains detected in patients with a variety of acute respiratory infections in Kumamoto, Japan in the period 2011-12. The phylogenetic tree and evolutionary timescale were obtained by the Bayesian Markov chain Monte Carlo method. Phylogenetic analyses showed that the present HRV-A, -B, and -C strains belonged to 25, 4, and 18 genotypes, respectively. Some new genotypes were confirmed as prevalent strains of HRV-C. An ancestor of the present HRV-ABCs could be dated back to about 20,000 years ago. The present HRV-A and -C strains have wide genetic divergence (pairwise distance >0.2) with rapid evolutionary rates (around 7 × 10(-4) to 4 × 10(-3)substitutions/site/year). Over 100 sites were found to be under negative selection, while no positively selected sites were found in the analyzed region. No evidence of recombination events was found in this region of the present strains. Our results indicate that the present HRV strains have rapidly evolved and subsequently diverged over a long period into multiple genotypes.

Molecular epidemiology of respiratory viruses in virus-induced asthma

Acute respiratory illness (ARI) due to various viruses is not only the most common cause of upper respiratory infection in humans but is also a major cause of morbidity and mortality, leading to diseases such as bronchiolitis and pneumonia. Previous studies have shown that respiratory syncytial virus (RSV), human rhinovirus (HRV), human metapneumovirus (HMPV), human parainfluenza virus (HPIV), and human enterovirus infections may be associated with virus-induced asthma. For example, it has been suggested that HRV infection is detected in the acute exacerbation of asthma and infection is prolonged. Thus it is believed that the main etiological cause of asthma is ARI viruses. Furthermore, the number of asthma patients in most industrial countries has greatly increased, resulting in a morbidity rate of around 10-15% of the population. However, the relationships between viral infections, host immune response, and host factors in the pathophysiology of asthma remain unclear. To gain a better understanding of the epidemiology of virus-induced asthma, it is important to assess both the characteristics of the viruses and the host defense mechanisms. Molecular epidemiology enables us to understand the pathogenesis of microorganisms by identifying specific pathways, molecules, and genes that influence the risk of developing a disease. However, the epidemiology of various respiratory viruses associated with virus-induced asthma is not fully understood. Therefore, in this article, we review molecular epidemiological studies of RSV, HRV, HPIV, and HMPV infection associated with virus-induced asthma.

Molecular evolution of attachment glycoprotein (G) gene in human respiratory syncytial virus detected in Japan 2008-2011

We investigated the evolution of the C-terminal 3rd hypervariable region of G gene in the prevalent human respiratory syncytial virus (RSV) subgroups A (RSV-A) and B (RSV-B) in Japan in 2008-2011. Phylogenetic analysis and the evolutionary timescale was obtained by the Bayesian Markov Chain Monte Carlo method. All 38 RSV-A strains detected were classified into genotype NA1 and the 17 RSV-B strains detected belonged to genotypes BA and GB2. NA1 subdivided around 1998 in the present phylogenetic tree. Genotype BA subdivided around 1994. The evolutionary rates for RSV-A and RSV-B were estimated at 3.63×10⁻³ and 4.56×10⁻³ substitutions/site/year, respectively. The mean evolutionary rate of RSV-B was significantly faster than that of RSV-A during all seasons. The pairwise distance was relatively short (less than 0.06). In addition, some unique sites under positive selection were found. The results suggested that this region of the RSV strains rapidly evolved with some unique amino acid substitutions due to positive pressure.

Proposals for the classification of human rhinovirus species A, B and C into genotypically assigned types

Human rhinoviruses (HRVs) frequently cause mild upper respiratory tract infections and more severe disease manifestations such as bronchiolitis and asthma exacerbations. HRV is classified into three species within the genus Enterovirus of the family Picornaviridae. HRV species A and B contain 75 and 25 serotypes identified by cross-neutralization assays, although the use of such assays for routine HRV typing is hampered by the large number of serotypes, replacement of virus isolation by molecular methods in HRV diagnosis and the poor or absent replication of HRV species C in cell culture. To address these problems, we propose an alternative, genotypic classification of HRV-based genetic relatedness analogous to that used for enteroviruses. Nucleotide distances between 384 complete VP1 sequences of currently assigned HRV (sero)types identified divergence thresholds of 13, 12 and 13 % for species A, B and C, respectively, that divided inter- and intra-type comparisons. These were paralleled by 10, 9.5 and 10 % thresholds in the larger dataset of >3800 VP4 region sequences. Assignments based on VP1 sequences led to minor revisions of existing type designations (such as the reclassification of serotype pairs, e.g. A8/A95 and A29/A44, as single serotypes) and the designation of new HRV types A101–106, B101–103 and C34–C51. A protocol for assignment and numbering of new HRV types using VP1 sequences and the restriction of VP4 sequence comparisons to type identification and provisional type assignments is proposed. Genotypic assignment and identification of HRV types will be of considerable value in the future investigation of type-associated differences in disease outcomes, transmission and epidemiology.

Levofloxacin inhibits rhinovirus infection in primary cultures of human tracheal epithelial cells

Respiratory virus infections, including infections with rhinoviruses (RVs), are related to exacerbations of chronic obstructive pulmonary disease (COPD). A new quinolone antibiotic, levofloxacin (LVFX), has been used to treat bacterial infections that cause COPD exacerbations as well as bacterial infections that are secondary to viral infection in COPD patients. However, the inhibitory effects of LVFX on RV infection and RV infection-induced airway inflammation have not been studied. We examined the effects of LVFX on type 14 rhinovirus (RV14) (a major human RV) infection of human tracheal epithelial cells pretreated with LVFX. LVFX pretreatment reduced the RV14 titer, the level of cytokines in the supernatant, the amount of RV14 RNA in the cells after RV14 infection, and the cells' susceptibility to RV14 infection. LVFX pretreatment decreased the mRNA level of intercellular adhesion molecule 1 (ICAM-1), a receptor for RV14, in the cells and the concentration of the soluble form of ICAM-1 in the supernatant before RV14 infection. LVFX pretreatment also decreased the number and the fluorescence intensity of the acidic endosomes from which RV14 RNA enters the cytoplasm. LVFX pretreatment inhibited the activation of nuclear factor κB proteins, including p50 and p65, in nuclear extracts. LVFX pretreatment did not reduce the titers of RV2 (a minor human RV) but reduced the titers of RV15 (a major human RV). These results suggest that LVFX inhibits major-group rhinovirus infections in part by reducing ICAM-1 expression levels and the number of acidic endosomes. LVFX may also modulate airway inflammation in rhinoviral infections.

Molecular epidemiology of the attachment glycoprotein (G) gene in respiratory syncytial virus in children with acute respiratory infection in Japan in 2009/2010

This study performed a detailed genetic analysis of the glycoprotein (G) gene of respiratory syncytial virus (RSV) detected in 50 Japanese children with acute respiratory infection (ARI) in the 2009/2010 season. A phylogenetic tree constructed by the neighbour-joining method showed that 34 and 16 of the RSV strains could be classified into subgroups A and B, respectively. Strains belonging to subgroups A and B were further subdivided into GA2 and BA, respectively. The nucleotide and deduced amino acid sequence identities were relatively high among these strains (>90%). The deduced amino acid sequences implied that a relatively high frequency of amino acid substitutions occurred in the C-terminal 3rd hypervariable region of the G protein in these strains. In addition, some positively selected sites were estimated. The results suggest that RSV with genotypes GA2 and BA was associated with ARI in Japanese children in 2009/2010.

Clinical characteristics and genetic variability of human rhinovirus in Mexico

Human rhinovirus (HRV) is a leading cause of acute respiratory infection (ARI) in young children and infants worldwide and has a high impact on morbidity and mortality in this population. Initially, HRV was classified into two species: HRV-A and HRV-B. Recently, a species called HRV-C and possibly another species, HRV-D, were identified. In Mexico, there is little information about the role of HRV as a cause of ARI, and the presence and importance of species such as HRV-C are not known. The aim of this study was to determine the clinical characteristics and genetic variability of HRV in Mexican children. Genetic characterization was carried out by phylogenetic analysis of the 5'-nontranslated region (5'-NTR) of the HRV genome. The results show that the newly identified HRV-C is circulating in Mexican children more frequently than HRV-B but not as frequently as HRV-A, which was the most frequent species. Most of the cases of the three species of HRV were in children under 2 years of age, and all species were associated with very mild and moderate ARI.

Detailed genetic analysis of hemagglutinin-neuraminidase glycoprotein gene in human parainfluenza virus type 1 isolates from patients with acute respiratory infection between 2002 and 2009 in Yamagata prefecture, Japan

Background

Human parainfluenza virus type 1 (HPIV1) causes various acute respiratory infections (ARI). Hemagglutinin-neuraminidase (HN) glycoprotein of HPIV1 is a major antigen. However, the molecular epidemiology and genetic characteristics of such ARI are not exactly known. Recent studies suggested that a phylogenetic analysis tool, namely the maximum likelihood (ML) method, may be applied to estimate the evolutionary time scale of various viruses. Thus, we conducted detailed genetic analyses including homology analysis, phylogenetic analysis (using both the neighbor joining (NJ) and ML methods), and analysis of the pairwise distances of HN gene in HPIV1 isolated from patients with ARI in Yamagata prefecture, Japan.

Results

A few substitutions of nucleotides in the second binding site of HN gene were observed among the present isolates. The strains were classified into two major clusters in the phylogenetic tree by the NJ method. Another phylogenetic tree constructed by the ML method showed that the strains diversified in the late 1980s. No positively selected sites were found in the present strains. Moreover, the pairwise distance among the present isolates was relatively short.

Conclusions

The evolution of HN gene in the present HPIV1 isolates was relatively slow. The ML method may be a useful phylogenetic method to estimate the evolutionary time scale of HPIV and other viruses.

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.

A molecular epidemiological study of respiratory viruses detected in Japanese children with acute wheezing illness

Background

Recent studies strongly suggest that some respiratory viruses are associated with the induction of acute wheezing and/or exacerbation of bronchial asthma. However, molecular epidemiology of these viruses is not exactly known.

Methods

Using PCR technology, we attempted to detect various respiratory viruses from 115 Japanese children. Furthermore, the detected viruses were subjected to homology, pairwise distance, and phylogenetic analysis.

Results

Viruses were detected from 99 (86.1%) patients. Respiratory syncytial virus (RSV) alone and human rhinovirus (HRV) alone were detected in 47 (40.9%) and 36 (31.3%) patients, respectively. Both RSV and HRV were detected in 14 (12.2%) patients. Human metapneumovirus (HMPV) alone and human parainfluenza virus (HPIV) alone were detected in 1 (0.9%) patient each, respectively. Homology and phylogenetic analyses showed that the RSV and HRV strains were classified into genetically diverse species or subgroups. In addition, RSV was the dominant virus detected in patients with no history of wheezing, whereas HRV was dominant in patients with a history of wheezing.

Conclusions

The results suggested that these genetically diverse respiratory viruses, especially RSV and HRV, might be associated with wheezing in Japanese children.