key: cord-312307-0hqqheho
authors: Ng, Kim Tien; Oong, Xiang Yong; Pang, Yong Kek; Hanafi, Nik Sherina; Kamarulzaman, Adeeba; Tee, Kok Keng
title: Outbreaks of enterovirus D68 in Malaysia: genetic relatedness to the recent US outbreak strains
date: 2015-08-05
journal: Emerg Microbes Infect
DOI: 10.1038/emi.2015.47
sha: 
doc_id: 312307
cord_uid: 0hqqheho

nan

temperature ( 6 C), relative humidity (%), number of rain days, and the amount of rainfall (mm) collected from the nearest meteorological station (latitude: 3 6 069N, longitude: 101 6 399E) was analyzed using Pearson's correlation. Specimens positive for enteroviruses were further confirmed using standard molecular approaches that involved amplification and sequencing of the human enterovirus VP4/VP2 gene using primers described previously. 10 PCR amplification of the capsid (P1) region (2480 bp) of EV-D68 was performed when VP4/VP2 gene was identified as EV-D68. A combination of previously published primers, namely 9895-forward, 9565-reverse, 10 59fwd1, 59rev1, 59fwd2, 59rev2, 11 and newly designed primers were used to amplify the P1 region by forming several overlapping fragments. These newly designed primers are EV68P1.5F1 (59-TCA AAA TTY ACT GAA CCA GT-39), EV68P1.5R1 (59-GTT GCC ATG AAG CTV CCA CA-39), EV68P1.5R2 (59-GAT ATG TTT CCT ACT ARA GT-39), EV68P1.3F1 (59-CCA GGG CAR GTC CGY AAC ATG-39), EV68P1.3R1 (59-CCA YTT GWA AAA GTT YTT GTC-39), EV68P1.3F2 (59-GTG GAR TCA ATG GAG AT-39), and EV68P1.3R2 (59-GCT GAT TTA TCA CYG TGC GAG-39). A total of 128 VP4/VP2, 376 VP1, and 20 P1 of EV-D68 retrieved from GenBank (accessed on 26 February 2015) were analyzed using neighbor-joining method implemented in MEGA version 6 to deduce the viral phylogenies. 12 The statistical robustness of the branching orders was evaluated by bootstrap analysis of 1000 replicates.

In this molecular epidemiological surveillance, approximately 51% (2009/3935) of patients had at least one viral pathogen detected by the multiplex assay, of whom 0.6% (12/2009) were infected with EV-D68. These EV-D68 cases were detected in the second half of 2012 (June to December) and between December 2013 and January 2014, of which September 2012 and January 2014 were the peak months of infection. No EV-D68 cases were detected in other months, suggesting transient outbreaks of EV-D68. The 12 patients (five males and seven females) from whom EV-D68 was detected ranged in age from 29 to 78 years old. During recruitment, most of the patients experienced mild sneezing and moderate-to-severe cough. Correlation of EV-D68 infections with meteorological factors was not observed (correlation coefficient ,0.3).

Based on previously described EV-D68 classification, 11 the newly sequenced strains from Malaysia were found within clade A (MY-Cluster-1) and clade B (MY-Cluster-2). Phylogenetic analysis of the P1 region indicated that 91.7% (11/12) of the Malaysian EV-D68 formed clusters, suggesting the transient EV-D68 outbreaks were most likely caused by at least two viral lineages ( Figure 1) . A sporadic EV-D68 strain (12MYKL1236) that did not cluster with the major lineages was also observed. Of note, MY-Cluster-1 contains EV-D68 mostly sampled in 2012, while MY-Cluster-2 contains EV-D68 sampled in 2013/2014. Such observation suggests an ongoing ''clade shift'' or lineage replacement of circulating EV-D68 in causing new outbreak, as observed in other enterovirus-associated outbreaks. 13 Sequence homology comparison based on the P1 region indicated that the Malaysian EV-D68 strains shared 86.6%-87.6% (nucleotide) and 93.4%-93.9% (amino acid) sequence identities with the Fermon strain (AY426531), a prototype strain isolated in the USA in 1962. Phylogenetic tree also showed that the representative EV-D68 sequences sampled from recent outbreaks in the USA were grouped into three distinct lineages. These lineages are designated as outbreak lineage 1 (consists of sequences isolated from Missouri, MO), lineage 2 (Illinois, IL and Kentucky, KY), and lineage 3 (KY) (Figure 1) . Interestingly, Malaysian EV-D68 sequences (MY-Cluster-2) was related to outbreak lineage 1 while outbreak lineage 3 (US/KY/14-18953) was clustered with the MY-Cluster-1, supported by strong statistical evidence (bootstrap values of 97% and 100%, respectively) at the internal tree nodes and genetic distance of f0.015. Likewise, outbreak lineage 1 was closely related to the EV-D68 sequences (CU134 and CU171) reported recently in Thailand. All members of MY-Cluster-2 and outbreak lineage 1 exhibited unique amino acid substitution at two residue positions (S647A and T650A), while members of MY-Cluster-1 and outbreak lineage 3 sequences displayed unique amino acid substitutions at four residue positions (N554E, G558E, F655Y, and I739V) relative to the Fermon strain (Figure 1 ). These ''clade-defining'' substitutions differentiate outbreak lineage 3 from other outbreak and non-outbreak lineages. A unique substitution T650A discriminated outbreak lineage 1 from outbreak lineage 2, further suggesting that the former is indeed closely related to the Thai strains. 5 Contribution of these amino acid substitutions on virus pathogenesis however remains unclear. Although the amino acid substitutions are located in the putative immunogenic BC and DE loops, experiments in susceptible cell systems are needed to prove the importance of such substitutions. 14 The divergence dates of the Malaysian outbreak and other major lineages were estimated by the Bayesian coalescent-based relaxed molecular clock model with constant population size, 11 implemented in the BEAST software. 15 . Our data suggest that the recent EV-D68 strains associated with unprecedented severe respiratory outbreaks in the USA in 2014 were probably descended from the recent EV-D68 lineages circulating in Thailand and Malaysia. However such observation remains presumptuous largely because limited EV-D68 sequence data originating from the USA (due to inadequate surveillance) are available for genealogical analysis. Given the close relationship of the Southeast Asian isolates with the US strains from recent outbreak, a more active enterovirus surveillance should be in place to monitor risks of impending outbreak.

Centers for Disease Control and Prevention. Enterovirus surveillance -United States

Clusters of acute respiratory illness associated with human enterovirus 68 -Asia

Severe respiratory illness associated with enterovirus D68 -Missouri and Illinois

Genome sequence of enterovirus D68 from

Genome sequence of enterovirus D68 and clinical disease

Seven strains of enterovirus D68 detected in the United States during the 2014 severe respiratory disease outbreak

Transmission of the common cold to volunteers under controlled conditions. I. The common cold as a clinical entity

Evaluation of NucliSens easyMAG for automated nucleic acid extraction from various clinical specimens

Comparison of the Luminex xTAG respiratory viral panel with xTAG respiratory viral panel fast for diagnosis of respiratory virus infections

Genetic clustering of all 102 human rhinovirus prototype strains: serotype 87 is close to human enterovirus 70

Worldwide emergence of multiple clades of enterovirus 68

MEGA6: molecular evolutionary genetics analysis version 6.0

Evolutionary genetics of human enterovirus 71: origin, population dynamics, natural selection, and seasonal periodicity of the VP1 gene

Can sequence data predict enterovirus D68 infection outcome?

BEAST: Bayesian evolutionary analysis by sampling trees