Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 22, Number 2—February 2016
Dispatch

Microevolution of Outbreak-Associated Middle East Respiratory Syndrome Coronavirus, South Korea, 2015

Author affiliations: Seoul National University Hospital, Seoul, South Korea (M.-W. Seong, T.S. Kim, S.I. Cho, M.J. Kim, S.J. Lee, J.-S. Lee, S.H. Seo, J.S. Ahn, B.S. Yu, N. Park, M.D. Oh, W.B. Park, E.C. Kim, S.S. Park); National Medical Center, Seoul (S.Y. Kim, J.Y. Lee, G. Kim, J.S. Joh, I. Jeong); University of Bonn Medical Centre, Bonn, Germany (V.M. Corman, C. Drosten); German Centre for Infection Research (DZIF), (V.M. Corman, C. Drosten)

Cite This Article

Abstract

During the 2015 Middle East respiratory syndrome coronavirus outbreak in South Korea, we sequenced full viral genomes of strains isolated from 4 patients early and late during infection. Patients represented at least 4 generations of transmission. We found no evidence of changes in the evolutionary rate and no reason to suspect adaptive changes in viral proteins.

Middle East respiratory syndrome coronavirus (MERS-CoV), first detected in Saudi Arabia in 2012, is a novel human pathogen that causes severe respiratory illness (1). Phylogenetic analyses and transmission studies suggest a zoonotic origin in dromedaries (2,3). Human-to-human transmission among close contacts of patients (e.g., family members and persons in healthcare settings) has been described (4). As of October 2015, the World Health Organization had received reports of 1,593 cases (including at least 568 deaths), most of which were reported from the Arabian Peninsula.

In South Korea, the first imported MERS-CoV case was identified on May 20, 2015, in a 68-year-old man who had traveled to the Middle East 2 weeks earlier (5). Another 185 persons in South Korea were subsequently infected during a 4-week period, mainly through in-hospital transmission. The unusually large number of cases, which occurred during at least 4 sequential generations of human-to-human transmission, raised questions regarding potential adaptations of MERS-CoV to the human host. To determine the possibility of virus adaptation, we repeatedly sequenced complete genomes for MERS-CoV from 4 patients representing different generations of transmission during the South Korea outbreak.

The Study

Figure 1

Thumbnail of Transmission tree timeline for 8 Middle East respiratory syndrome coronavirus strains isolated during an outbreak in South Korea, 2015. Numbers without parentheses indicate patients in this study; numbers inside parentheses indicate patients not included in this study. The index case-patient is represented by (1). Numbers 1 and 2 following patient identification numbers indicate separate samples that were sequenced. The left edge of each shaded box indicates date of symptom onset fo

Figure 1. Transmission tree timeline for 8 Middle East respiratory syndrome coronavirus strains isolated during an outbreak in South Korea, 2015. Numbers without parentheses indicate patients in this study; numbers inside parentheses indicate...

Institutional review boards of the Seoul National University Hospital and the National Medical Center approved this study. We tested samples from 4 patients, designated as patients 14, 35, 163, and 168 (Table 1; Technical Appendix Table 1). Patients 14, 35, and 168 were second-, third-, and fourth-generation case-patients, respectively; each had recorded exposure histories (5) (Figure 1). Patient 163 had a recorded transmission history that traced back to patient 119. However, for 1 interim transmission, the place and approximate time of exposure could be reconstructed, but individual contacts could not be determined. Thus, patient 163 might belong to at least the fourth, but potentially the fifth, generation of transmission (5).To identify virus changes, we obtained 2 samples from each of the 4 patients, 1 at the early and 1 at the late stage of infection. Clinical samples were tested for MERS-CoV RNA by reverse transcription PCR (RT-PCR) targeting the upE (upstream of E) and ORF1 (open reading frame 1) genes (6). Using virus quantity estimates as a basis, we determined full genomes by amplifying overlapping PCR products and sequencing as previously described (7).

Genome sequences with a minimum length of 29,831 bp (99.04% of the genome) were obtained from each specimen and submitted to GenBank (accession nos. KT374050– KT374057). All sequences clustered phylogenetically with MERS-CoV strains identified during the outbreak and with a sequence from a linked case in China (8,9). Our findings confirmed previously described evidence for recombination between MERS-CoV clades (10). MERS-CoV strain Hu/Riyadh KSA_2959_2015, the closest related strain outside the South Korea outbreak, was used as an outgroup to reconstruct the phylogeny of the 8 viral genomes (Technical Appendix Figure); the strain was isolated in Riyadh, Saudi Arabia, during February 2015.

Figure 2

Thumbnail of Location of the 13 variant nucleotide positions identified in genomes of 8 Middle East respiratory syndrome coronavirus strains isolated during an outbreak in South Korea, 2015. Case-patient numbers are indicated on the left, as are GenBank accession numbers for 2 related strains. For each case, the top and bottom horizontal bars represent the genome sequence generated from the first and second samples, respectively. Letters indicate matching nucleotide positions between samples. Sa

Figure 2. Location of the 13 variant nucleotide positions identified in genomes of 8 Middle East respiratory syndrome coronavirus strains isolated during an outbreak in South Korea, 2015. Case-patient numbers are indicated on...

The 8 strains from South Korea shared 99.8%–99.9% nt identity with Hu/Riyadh KSA_2959_2015, deviating by 24–27 positions across the genome. Sequences for the viruses showed 13 variant nucleotide positions: 6 in the ORF1ab gene, 5 in the spike gene, and 1 each in accessory genes ORF4a and ORF5 (Figure 2). Of the 13 variants, 11 caused nonsynonymous substitutions (Technical Appendix Table 2). To analyze substitutions along the transmission tree, we reconstructed the unknown sequence of the index case-patient’s virus based on an isolate from his wife (GenBank accession no. KT029139), a virus from a patient who traveled to China (accession no. KT006149), the co-ancestral strain from Saudi Arabia (accession no. KT026453), and all sequences determined in this study. Considering that sequencing errors were possible in strains KT029139 and KT006149 due to cell culture adaptation and differences in sequencing technique (8,9), and assuming no mutation reversion within short human-to-human passage, an unequivocal ancestral sequence reconstruction was possible (11). The index case-patient showed overt symptoms by May 11, indicating virus exposure occurred around May 1 (5). This date coincides with days (May 1 and 2) that the index case-patient spent in Saudi Arabia, where, at the time, viruses most closely related to the South Korea strains were circulating (10). Thus, that starting sequence, as determined by ancestral sequence reconstruction, was assumed to have existed on May 1.

Substitutions observed along the transmission tree corresponded, on average, to 3.78 × 10−6 substitutions/site/day (Table 2). The average substitution rate within individual patients (comparing first and second samples for each patient) was not significantly different: 3.44 × 10−6 substitutions/site/day. The between-patient evolutionary rate was projected to be 1.3 × 10−3 substitutions/site/year, which is within the CI of an earlier estimate of the evolutionary rate for MERS-CoV: 1.12 × 10−3 substitutions/site/year (95% CI 8.76 × 10−4 to 1.37 × 10−3) (12). This finding suggests that the exclusively human-to-human transmission observed in this outbreak has not led to a change in the apparent evolutionary rate of the virus. No similar mutations occurred in parallel transmission chains leading up to patients 35, 168, and 163, suggesting that quasispecies sampling during transmission events has been stochastic rather than selective. Selective sampling would have been expected in the hypothetical case of emergence of a mutated virus with increased transmissibility or replication level.

The MERS-CoV spike glycoprotein targets the cellular receptor DPP4 (dipeptidyl peptidase 4). The receptor-binding domain (RBD) consists of residues E382 to C585 (13,14). In this study, we identified 2 RBD variants: D510G and I529T. Recent mutagenesis studies indicate that alterations of key residues within this region (D510A and E513A) could substantially reduce the efficiency of binding and virus entry (13,15). Because aspartate and glycine have similar physicochemical properties, the D510G variant found in this study might resemble D510A in its potential to reduce receptor binding.

Patient 168, a 36-year-old man who worked as a radiology technologist at a university hospital in Seoul, South Korea, carried the D510G variant virus. He was exposed to the virus via patient 76 on June 6. Symptoms developed on June 17, and he was immediately admitted to an isolation room in the hospital. RT-PCR was conducted at admission, but MERS-CoV RNA was not detected until 2 days later, suggesting an initially low level of virus replication. After MERS-CoV infection was confirmed by RT-PCR, combined interferon and ribavirin treatment was administered. The patient’s vital signs were stable, and no lung lesions were seen on chest radiographs during hospitalization. He was discharged on June 30 after 2 consecutively negative RT-PCR results (June 28 and 30). The patient’s mild clinical course and indicators of a low level of virus replication warrant further virologic study of a potentially attenuating effect of the D510G variant.

Another RBD variant, I529T, was carried by study patients 14 and 35. Although further studies are needed, physicochemical considerations suggest that this variant is unlikely to affect virus binding. Of note, patient 14 transmitted the virus to >80 tertiary case-patients and was thus considered a superspreader; however, persons infected by this patient, including patient 35, were not superspreaders.

Conclusions

Our study is limited by the absence of full-genome information for the index case-patient’s virus and by the effect of outbreak interventions that may have concealed phenotypic virus changes. However, we found no evidence for an increase in the evolutionary rate or for adaptive changes over at least 4 generations of transmission. Changes that were observed were likely caused by transmission bottleneck effects.

Dr. Seong is a clinical pathologist in the Department of Laboratory Medicine of Seoul National University Hospital in South Korea. His main areas of work are in the molecular diagnostics of emerging infectious diseases.

Top

Acknowledgment

C.D. has been supported by the European Commission research grant PREPARE and by the Deutsche Forschungsgemeinschaft (grant DR772/10-1).

Top

References

  1. Zaki  AM, van Boheemen  S, Bestebroer  TM, Osterhaus  AD, Fouchier  RA. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N Engl J Med. 2012;367:181420. DOIPubMedGoogle Scholar
  2. Reusken  CB, Haagmans  BL, Muller  MA, Gutierrez  C, Godeke  GJ, Meyer  B, Middle East respiratory syndrome coronavirus neutralising serum antibodies in dromedary camels: a comparative serological study. Lancet Infect Dis. 2013;13:85966. DOIPubMedGoogle Scholar
  3. Memish  ZA, Cotten  M, Meyer  B, Watson  SJ, Alsahafi  AJ, Al Rabeeah  AA, Human infection with MERS coronavirus after exposure to infected camels, Saudi Arabia, 2013. Emerg Infect Dis. 2014;20:10125. DOIPubMedGoogle Scholar
  4. Drosten  C, Meyer  B, Muller  MA, Corman  VM, Al-Masri  M, Hossain  R, Transmission of MERS-coronavirus in household contacts. N Engl J Med. 2014;371:82835. DOIPubMedGoogle Scholar
  5. Cowling  BJ, Park  M, Fang  VJ, Wu  P, Leung  GM, Wu  JT. Preliminary epidemiological assessment of MERS-CoV outbreak in South Korea, May to June 2015. Euro Surveill. 2015;20:713. DOIPubMedGoogle Scholar
  6. Corman  VM, Muller  MA, Costabel  U, Timm  J, Binger  T, Meyer  B, Assays for laboratory confirmation of novel human coronavirus (hCoV-EMC) infections. Euro Surveill. 2012;17:20334 .PubMedGoogle Scholar
  7. Corman  VM, Ithete  NL, Richards  LR, Schoeman  MC, Preiser  W, Drosten  C, Rooting the phylogenetic tree of Middle East respiratory syndrome coronavirus by characterization of a conspecific virus from an African bat. J Virol. 2014;88:11297303. DOIPubMedGoogle Scholar
  8. Lu  R, Wang  Y, Wang  W, Nie  K, Zhao  Y, Su  J, Complete genome sequence of Middle East respiratory syndrome coronavirus (MERS-CoV) from the first imported MERS-CoV case in China. Genome Announc. 2015;3:e00818–15.
  9. Kim  YJ, Cho  YJ, Kim  DW, Yang  JS, Kim  H, Park  S, Complete genome sequence of Middle East respiratory syndrome coronavirus KOR/KNIH/002_05_2015, isolated in South Korea. Genome Announc. 2015;3:e00787–15.
  10. Wang  Y, Liu  D, Shi  W, Lu  R, Wang  W, Zhao  Y, Origin and possible genetic recombination of the Middle East respiratory syndrome coronavirus from the first imported case in China: phylogenetics and coalescence analysis. MBio. 2015;6:e0128015. DOIPubMedGoogle Scholar
  11. Cotten  M, Watson  SJ, Kellam  P, Al-Rabeeah  AA, Makhdoom  HQ, Assiri  A, Transmission and evolution of the Middle East respiratory syndrome coronavirus in Saudi Arabia: a descriptive genomic study. Lancet. 2013;382:19932002. DOIPubMedGoogle Scholar
  12. Cotten  M, Watson  SJ, Zumla  AI, Makhdoom  HQ, Palser  AL, Ong  SH, Spread, circulation, and evolution of the Middle East respiratory syndrome coronavirus. MBio. 2014;5:e0106213. DOIPubMedGoogle Scholar
  13. Wang  N, Shi  X, Jiang  L, Zhang  S, Wang  D, Tong  P, Structure of MERS-CoV spike receptor-binding domain complexed with human receptor DPP4. Cell Res. 2013;23:98693. DOIPubMedGoogle Scholar
  14. Lu  L, Liu  Q, Zhu  Y, Chan  KH, Qin  L, Li  Y, Structure-based discovery of Middle East respiratory syndrome coronavirus fusion inhibitor. Nat Commun. 2014;5:3067.
  15. Du  L, Zhao  G, Yang  Y, Qiu  H, Wang  L, Kou  Z, A conformation-dependent neutralizing monoclonal antibody specifically targeting receptor-binding domain in Middle East respiratory syndrome coronavirus spike protein. J Virol. 2014;88:704553. DOIPubMedGoogle Scholar

Top

Figures
Tables

Top

Cite This Article

DOI: 10.3201/eid2202.151700

1These authors contributed equally to this article.

Table of Contents – Volume 22, Number 2—February 2016

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.

Top

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Sung Sup Park, Department of Laboratory Medicine, Seoul National University Hospital, 101, Daehak-ro, Jongno-gu, Seoul 110-744, South Korea; ; Christian Drosten, Institute of Virology, University of Bonn Medical Centre, Sigmund Freud Str 25, 53105 Bonn, Germany

Send To

10000 character(s) remaining.

Top

Page created: January 15, 2016
Page updated: January 15, 2016
Page reviewed: January 15, 2016
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external