








key: cord-326816-a7yovhvk
authors: Piralla, Antonio; Principi, Nicola; Ruggiero, Luca; Girello, Alessia; Giardina, Federica; De Sando, Elisabetta; Caimmi, Silvia; Bianchini, Sonia; Marseglia, Gian Luigi; Lunghi, Giovanna; Baldanti, Fausto; Esposito, Susanna
title: Enterovirus-D68 (EV-D68) in pediatric patients with respiratory infection: The circulation of a new B3 clade in Italy
date: 2018-01-12
journal: J Clin Virol
DOI: 10.1016/j.jcv.2018.01.005
sha: 
doc_id: 326816
cord_uid: a7yovhvk

BACKGROUND: In recent years, several outbreaks due to Enterovirus D-68 (EV-D68) have been reported, and it was confirmed that the virus can cause upper and lower respiratory tract diseases and be associated with the development of neurological problems. OBJECTIVES: The main aim of this research was to study the genetic characteristics of EV-D68 strains that were circulating in Italy identified during an outbreak of an EV-D68 infection that occurred in Italy during the period March-October 2016. STUDY DESIGN: A retrospective study of the circulation of different types and subtypes of EV-D68 was performed. Nasopharyngeal swabs were collected from March 2016 through October 2016 in children admitted to the Emergency Room with respiratory diseases. RESULTS: Among 390 children, 22 (59.1% males; mean age 47 months) were found to be infected by EV-D68 and most of them were immunocompetent (72.7%). Pneumonia was diagnosed in 12 (54.5%) children. Phylogenetic analysis of the VP1 region showed that all the strains identified in this study belonged to clade B3. Within B3 subclade, the Italian EV-D68 strains were most closely related to strains detected in Southern China in 2015 as well as to strains detected in US and the Netherlands in 2016. CONCLUSIONS: These results showed that EV-D68 infections are a common cause of lower respiratory illness in pediatric age. The circulation of one EV-D68 lineage has been proven in Italy and in the European region during 2016. However, further studies are required to investigate whether some strains or lineages may possess a higher affinity for the lower airway or central nervous system.

Human enterovirus D68 (EV-D68) was originally isolated in 1962, and for many years it was considered a rare cause of respiratory disease because it was identified in isolated cases and in small outbreaks [1] . Although outbreaks of EV-D68 infection occurred in the Philippines [2] , Italy [3, 4] , the Netherlands [5] and Japan [6] between 2005 and 2012, it was only in 2014 that the virus gained epidemiological and clinical relevance. In that year, a large-scale outbreak of severe respiratory infections due to this agent, in some cases associated with central nervous system diseases, occurred in the USA and Canada [7] , with subsequent cases in other countries [8] [9] [10] [11] [12] [13] [14] [15] [16] . Since that period, several outbreaks in different continents have been reported, and it was confirmed that the virus can cause upper and lower respiratory tract diseases and be associated with the development of neurological problems, including acute flaccid paralysis (AFP), encephalitis and aseptic meningitis [17] [18] [19] .

Phylogenetic analysis of EV-D68 strains detected worldwide has led to the identification of four distinct viral clades that are designated groups A-D according to the characteristics of the structural protein VP1, the most variable region of the EV genome [20] . Moreover, within clades, subgroups were identified according to amino acid substitutions https://doi.org/10.1016/j.jcv.2018.01.005 in structural and nonstructural proteins. Several studies have shown that these different clades and subgroups can circulate or co-circulate during different periods and be associated with both mild and severe diseases [21] . However, little is known about the molecular evolution of this virus during outbreaks, and it is not precisely defined whether specific genetic variances are responsible for the most severe cases. Studies regarding the molecular phylogeny, diversity, and epidemiology of EV-D68 are essential to understanding the evolution of this virus and its adaptation to the human host. Moreover, they can contribute to the development of effective preventive and therapeutic measures.

The main aim of this research was to study the genetic characteristics of EV-D68 strains that were circulating in Italy during the period March-October 2016.

A retrospective study of the circulation of different types and subtypes of EV-D68 was performed. Nasopharyngeal swabs collected from March 2016 through October 2016 in children admitted to the Emergency Room with respiratory diseases and stored in a freezer at −80°C were tested. The sample collection had occurred in the Pediatric Department of the Fondazione IRCCS Policlinico San Matteo, University of Pavia (Italy) and in the Pediatric Highly Intensive Care Unit of the Fondazione IRCCS Ca' Granda Policlinico, University of Milan (Italy). The Ethics Committees of both Institutions approved the study. Written informed consent was obtained from a parent of a legal guardian of each enrolled child. Children aged ≥8 years were asked for their assent. Clinical data were retrieved from archived clinical charts.

Respiratory samples were extracted using a QIAsymphony ® instrument with a QIAsymphony ® DSP Virus/Pathogen Midi Kit (Complex 400 protocol) according to the manufacturer's instructions (QIAGEN, Qiagen, Hilden, Germany). Panels of laboratory-developed real-time RT-PCR or real-time PCR [22, 23] were used to detect and quantify the following viruses: influenza virus A and B (including subtype determination), human Rhinovirus (hRV), human parainfluenza virus (hPIV) 3, respiratory syncytial virus (RSV) types A and B, human coronavirus (hCoV)-OC43, -229E, -NL63, and -HKU1, human metapneumovirus (hMPV) and human adenoviruses (hAdVs). Extracted samples were also amplified using EV-D68-specific real-time RT-PCR as previously described [24] .

Complete VP1 sequences were amplified for all EV-D68-positive samples according to the method of Piralla et al. [4] . Purified PCR products were sequenced using a BigDye Terminator Cycle-Sequencing Kit (Applied Biosystems, Foster City, USA) in an ABI Prism 3130xl Genetic Analyzer (Applied Biosystems, Foster City, USA).

The sequences were assembled using Sequencher software, version 4.6 (Gene Codes Corporation, Ann Arbor, USA). Sequence alignment and phylogenetic analysis were performed using the MEGA software, version 7 with the TN92 + G nucleotide substitution model, which was selected using the hierarchical likelihood ratio test [25] . Statistical support of specific clades in the inferred tree was evaluated by using the bootstrap analysis (1000 replicates).The VP1 secondary structure was investigated using the Jpred4 server [26] .

The nucleotide dataset of the VP1 gene was obtained by including sequences from complete genome (n = 84), partial VP1 (n = 20) from the Netherlands outbreak 2016 [27] as well as from EV-D68 strain sequences that were circulating in Italy during the period 2008-2014 (n = 30). The EV-D68 sequences originated in this study were deposited in the GenBank database (accession numbers: MF073335-MF073354).

Tests for positive selection were conducted using single-likelihood ancestor counting (SLAC), fixed-effects likelihood (FEL), random-effects likelihood (REL), internal branch fixed-effects likelihood (IFEL), the mixed-effects model of evolution (MEME), and fast unconstrained Bayesian approximation (FUBAR) on the Datamonkey server [28] , with the dN/dS ratios being calculated using the SLAC and FEL codon-based maximum likelihood approaches. SLAC counts the number of non-synonymous changes per non-synonymous site (dN) and tests whether it is significantly different from the number of synonymous changes per synonymous site (dS). FEL estimates the ratios of non-synonymous to synonymous changes for each site in an alignment [29] . The IFEL method is similar to FEL but tests for site-by-site selection only along internal branches of the phylogeny. To avoid an excessive false-positive rate, sites with SLAC, FEL, IFEL and MEME p-values of < 0.1 and a FUBAR posterior probability of > 0.90 were accepted as candidates for selection.

Samples obtained from 390 children, 270 from Pavia and 120 from Milan, were tested. Twenty-two patients (13 males; 59.1%), mean age 47 months (range, 1-174 months) were found to be infected by EV-D68 ( Table 1 ). The majority of cases (14/22; 63.6%) occurred in June (n = 9) and July (n = 5). The remaining cases were observed in April (n = 1), May (n = 2), September (n = 2), and October (n = 3). Pneumonia was diagnosed in 12 (54.5%) children, whereas viral wheezing, rhinopharyngitis, asthma and bronchiolitis, were diagnosed in four children, three children, two children, and one child, respectively. The previous clinical history was negative in most of the cases (16/22; 72.7%). However, two (9.1%) children previously suffered from leukemia and had been recently transplanted with hematopoietic stem cells; 2 (9.1%) children had a history of recurrent episodes of viral wheezing; and 2 (9.1%) children had allergic asthma. The disease severity was judged to be mild-to-moderate in all the patients. Lower than 92% oxygen saturation (SaO 2 ) values at the first visit in the emergency room were found only in 5 (22.7%) children. Hospitalization occurred in 19/22 (86.4%) children, with a mean duration of 4 days (range, 4-8 days). None of the patients had to be admitted to the intensive care unit. In 16/22 (72.7%) cases, EV-D68 was the only infectious agent detected in the respiratory samples. Coinfection with human adenovirus (hAdV) and with human rhinovirus (hRV) was detected in 4 (18.2%) and 2 (9.1%) of the children, respectively. In all these cases, coinfecting agents were detected in very low concentrations (almost all below 1.0 × 10 3 copies/mL). The median EV-D68 load observed was 1.3 × 10 5 RNA copies/mL of respiratory sample (range, 9.0 × 10 2 -5.1 × 10 6 RNA copies/mL of respiratory sample). No association was found between disease severity (measured by the SaO 2 values) and EV-D68 load. The 368 samples negative for EV-D68 resulted positive for hRV (60.2%), RSV (26.6%), bocavirus (10.1%), hMPV (8.2%), and hAdV (7.9%). In 48 (13%) cases, particularly when hRV was detected, coinfection between two viruses was evidenced. Ten children with RSV infection, two with hMPV and one with hAdV were hospitalized for severe disease and three of those with RSV required admission to the intensive care unit.

Complete VP1 sequencing was obtained for all the EV-D68-positive patients. Phylogenetic analysis of the VP1 region showed that all the strains identified in this study belonged to clade B3, while EV-D68 strains circulating in Italy in the period 2008-2014, highlighted in red, clustered in the others clades and subclades (Fig. 1) . Within B3 subclade, the Italian EV-D68 strains were most closely related to strains detected in US and the Netherlands in 2016 and Southern China and Japan in 2014-2015.

The average nucleotide identity between Italian EV-D68 subclade B3 strains was 98.9% (range, 96.9-100%), with a maximum of 21 nucleotides and 3 amino acid changes Eight β-sheets (βB to βI) were predicted in the VP1 sequences by the Jpred4 server. Among B3 subclade were also clustered EV-D68 strains identified in Japan and China between 2014 and 2015. In addition, the average nucleotide identity between EV-D68 subclade B3 strains and B1 and B2 subclade strains was 92.3% and 94.7%, respectively. Inside the motifs corresponding to the BC and DE loops, no changes were observed compared with other EV-D68 strains belonging to clade B1-3 (Fig. 2) . The EV-D68 strains belonging to the new circulated B3 clade are characterized by 2 unique substitutions (A6G, S218T; Fig. 2 ). Codon 218 is included in the GH loop that is also involved in the sialic acid receptor binding site [30] . In addition, four strains harbored the I258 V mutation (highlighted with an asterisk in Fig. 1 ) that was only present in the EV-D68 ATCC strains (EVD68/Homo sapiens/XXX/ATCC VR-1197; KT725431). One strain (ITA/MI-ESC128/2016) was characterized by the K256I change, which had never been previously observed. Of note, one strain (ITA/PV-38993/2016) was characterized by the E168K change. This mutation was also present in the Fermon reference strain (AY4263531) but was more present in a few EV-D68 sequences belonging to all the three EV-D68 clades.

To assess the signatures of evolution in the VP1 gene region, the overall dN/dS ratio was estimated as 0.08. Several codons (more than 150) were identified to be under negative selection by at least three of the methods used (SLAC, FEL, REL, FUBAR or MEME). Similarly, no positively selected codons were identified by the IFEL model used to determine the selection pressure acting on the codons along the internal branches of the tree.

This retrospective study did not test all the nasopharyngeal samples collected in children with respiratory symptoms during the period March-October 2016 in Pavia and Milan. Consequently, no definitive data regarding the epidemiology and clinical relevance of this infectious agent in children could be drawn. However, the collected data enables the integration of available information regarding the circulation of and genetic differences among different EV-D68 clades and can help elucidate the molecular evolution of this virus. The four EV-D68 clades have exhibited a differing predominance in recent years. The Italian EV-D68 strains identified in 2008 belonged to clade A and clade C. Clade A was again identified in 2010, but in 2012, clade B became predominant [4] . During the 2014 EV-D68 epidemic in Europe, most of the 205 sequenced viruses collected in 17 countries in the second semester of the year were found to belong to B1 and B2 clades, whereas A clades, although co-circulating, were detected in a smaller number of patients [31] . Since April 2016, EV-D68 infection outbreaks have been reported in several European countries, particularly in Northern Europe [31] [32] [33] . Phylogenetic analyses showed that these viruses belonged to the B3 lineage [32] . The same results were obtained with this study, suggesting that all the 2016 European outbreaks originated from a common source and that the B3 clade has emerged and become preeminent, eclipsing clades A and clades B1 and B2. The worldwide importance of clade B3 appears to be confirmed by the evidence that outbreaks of infections due to EV-D68 infection belonging to clade B3 have been reported in 2015 in China [34] , and almost simultaneously with the European epidemics, in the USA [35] .

In most of the cases, EV-D68 was detected as a single infectious agent. Moreover, when a coinfection was demonstrated, the co-infecting virus was found at a very low concentration, suggesting a bystander role. This result highlights the pathogenetic role of EV-D68, although a strict relationship between viral load and infection site or disease severity is not evidenced. The relatively low number of tested children with EV-D68 might explain this limitation. In this study, EV-D68 infections were diagnosed mainly in June and July and in pre-school children. These findings agree with those reported by other authors, who found a higher incidence of EV-D68 diseases in warm months and in younger subjects [36] [37] [38] [39] . However, in this study, differently from other viral infections, particularly those due to RSV, the EV-D68-related respiratory infections were all of mild or moderate severity, and no patients showed signs and symptoms of nervous system involvement. Even when SaO 2 levels were < 92%, the respiratory signs and symptoms rapidly subsided, and no patients were admitted to the intensive care unit. This result is in contrast to the findings reported by Dyrdak et al. who, among the 74 patients with EV-D68 clade B3 infection they identified, observed ten cases with severe respiratory or neurological symptoms and one death [33] . Knoester et al. also reported a high incidence of severe disease among children infected by EV-D68 [32] . These authors had to admit to the intensive care unit 13/17 (76.4%) of the children seen with respiratory disease. Moreover, in one case, AFP developed. By contrast, a child with EV-D68 B3 clade respiratory infections associated with acute flaccid paralysis has been recently described in Italy, in the same region and in the same period of the cases reported in this study [19] .

In the past, it has been proposed that the development of a severe clinical picture following EV-D68 infection is the consequence of mutations occurring in the virus. In a study in which EV-D68 B1 strains associated with AFP cases diagnosed in the USA in 2014 were examined, 6 polymorphisms were identified and potentially associated with neuropathogenesis [40] . Other viral variants were identified in EV-D68 strains associated with severe disease in Taiwan [41] ; moreover, in this case, it was proposed that the variants were associated with the severity of EV-D68. However, these findings have not been completely confirmed by other studies [20, 42] . By contrast, based on our strains sequences analyses, no evidence of increased clinical severity associated with specific molecular signatures in VP1 sequences has emerged from the present study. Further analysis on full-genome sequencing in a large cohort of positive patients with a well-define clinical syndrome could allow to correlated specific molecular signatures with clinical syndromes.

Based on adaptive evolution analysis, the overall dN/dS ratio (dN/ dS < 1) of the VP1 gene observed in this study was similar to those reported by Lau et al. [43] Although several amino acid changes have been observed in VP1 sequences the meaning of dN/dS ratio reveal a purifying (negative) evolution for the VP1 protein. In addition, in keeping with this report, no codons with a significant signal of positive selective pressure were detected in the VP1 gene. By contrast, the phylogenetic analysis showed how the genetic changes may have driven the emergence of a new lineage of EV-D68 during the last 2 years. In fact, the worldwide circulation of EV-D68 strains belonging to the B2 and B3 lineages might have evolved from a common ancestor originating in 2011-2012 [44] . Recently, the evolutionary analyses performed on complete genome sequences have suggested that subclade B3 was not originated by a mutation occurred on subclade B1 strains. Indeed, it has been suggested a distinct evolutionary origin for the new B3 subclade [45] . Multiple clades of EV-D68 have worldwide circulated simultaneously between 2005 and 2014 [9] , while a predominance of one clade (subclade B3) has recently been observed. Additional both epidemiological and molecular studies are needed to better clarify the reasons of this EV-D68 strain selection and widespread diffusion. It is interesting to note that the emergence of new EV-D68 lineages does not represent the circulation of a new virus. Therefore, diagnostic real-time assays, designed in a well conserved region (5′ UTR), are capable to detect all the new EV-D68 variants.

In conclusion, EV-D68 infections have been confirmed as a cause of lower respiratory illness in pediatric age. The circulation of one EV-D68 lineage has been proven in Italy as well as in the European region during 2016. However, further studies are required to investigate whether some strains or lineages may possess a higher affinity for the lower airway or central nervous system. In addition, the continuous global monitoring of the clinical and molecular epidemiology of EV-D68 should be improved to understand the factors determining possible spatiotemporal differences of EV-D68 infections. 

Centers for disease control and prevention. Enterovirus surveillance-United States

Enterovirus 68 among children with severe acute respiratory infection, the Philippines

Human rhinovirus and human respiratory enterovirus (EV68 and EV104) infections in hospitalized patients in Italy

Phylogenetic characterization of enterovirus 68 strains in patients with respiratory syndromes in Italy

Emergence and epidemic occurrence of enterovirus 68 respiratory infections in The Netherlands in 2010

Amino acid changes in VP1 sequences according to different clades. Amino acid residues characterizing clade B3 EV-D68 strains are indicated with black square. Amino acids included in VP1 motifs corresponding to protein loops are highlighted with a dashed line. The gaps are indicated by a dash (−) and the conserved amino acid residues by a dot (.). Numbering of amino acidic changes is according

Acute respiratory infections due to enterovirus 68 in Yamagata

outbreak of enterovirus D68 in North America

Enterovirus D68: a focused review and clinical highlights from the 2014 United States outbreak

Global reemergence of enterovirus D68 as an important pathogen for acute respiratory infections

Severe respiratory illness associated with enterovirus D68-Missouri and Illinois

Continued seasonal circulation of enterovirus D68 in the Netherlands

Molecular epidemiology of enterovirus D68 from 2013 to 2014 in Philippines

Genome sequence of enterovirus D68 and clinical disease

Genome sequence of enterovirus D68 from St Louis

Enterovirus D68-associated severe pneumonia

Low-level circulation of enterovirus D68-associated acute respiratory infections

A fatal central nervous system enterovirus 68 infection

Enterovirus-D68 in the cerebrospinal fluid of two children with aseptic meningitis

Acute flaccid myelitis associated with enterovirus-D68 infection in an otherwise healthy child

Molecular evolution and intraclade recombination of enterovirus D68 during the 2014 outbreak in the United States

European surveillance for enterovirus D68 during the emerging North-American outbreak in 2014

Phylogenetic patterns of human respiratory picornavirus species, including the newly identified group C rhinoviruses, during a 1-year surveillance of a hospitalized patient population in Italy

FilmArray ® respiratory panel performance in respiratory samples from neonatal care units

A new real-time reverse transcription-PCR assay for detection of human enterovirus 68 in respiratory samples

MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets

JPred4: a protein secondary structure prediction server

Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology

Not so different after all: a comparison of methods for detecting amino acid sites under selection

Sialic acid-dependent cell entry of human enterovirus D68

European surveillance for enterovirus D68 during the emerging North-American outbreak in 2014

Outbreak of enterovirus D68 of the new B3 lineage in

Identification and whole-genome sequencing of four enterovirus D68 strains in Southern China inlate

Complete genome sequences of nine enterovirus D68 strains from patients of the lower Hudson Valley

Enterovirus D68 in hospitalized children: sequence variation, viral loads and clinical outcomes

Clinical course of enterovirus D68 in hospitalized children

Enterovirus D68 infection among children with medically attended acute respiratory illness

Enterovirus D68-associated community-acquired pneumonia in the pediatric age group

A novel outbreak enterovirus D68 strain associated with acute flaccid myelitis cases in the USA (2012-14): a retrospective cohort study

Molecular evolution and the global reemergence of enterovirus D68 by genome-wide analysis

Genetic divergence of enterovirus D68 in China and the United States

Enterovirus D68 infections associated with severe respiratory illness in elderly patients and emergence of a novel clade in Hong Kong

Global distribution and evolutionary history of enterovirus D68, with emphasis on the

Enterovirus D68 subclade B3 strain circulating and causing an outbreak in the United States in 2016

The authors declare that they have no potential conflict of interest.

The Ethics Committees of Fondazione IRCCS Policlinico San Matteo, University of Pavia (Italy) and the Pediatric Highly Intensive Care Unit of the Fondazione IRCCS Ca' Granda Policlinico, University of Milan (Italy) approved the study. Written informed consent was obtained from a parent of a legal guardian of each enrolled child. Children aged ≥8 years were asked for their assent.

