Spread, Circulation, and Evolution of the Middle East Respiratory Syndrome Coronavirus

ABSTRACT The Middle East respiratory syndrome coronavirus (MERS-CoV) was first documented in the Kingdom of Saudi Arabia (KSA) in 2012 and, to date, has been identified in 180 cases with 43% mortality. In this study, we have determined the MERS-CoV evolutionary rate, documented genetic variants of the virus and their distribution throughout the Arabian peninsula, and identified the genome positions under positive selection, important features for monitoring adaptation of MERS-CoV to human transmission and for identifying the source of infections. Respiratory samples from confirmed KSA MERS cases from May to September 2013 were subjected to whole-genome deep sequencing, and 32 complete or partial sequences (20 were ≥99% complete, 7 were 50 to 94% complete, and 5 were 27 to 50% complete) were obtained, bringing the total available MERS-CoV genomic sequences to 65. An evolutionary rate of 1.12 × 10−3 substitutions per site per year (95% credible interval [95% CI], 8.76 × 10−4; 1.37 × 10−3) was estimated, bringing the time to most recent common ancestor to March 2012 (95% CI, December 2011; June 2012). Only one MERS-CoV codon, spike 1020, located in a domain required for cell entry, is under strong positive selection. Four KSA MERS-CoV phylogenetic clades were found, with 3 clades apparently no longer contributing to current cases. The size of the population infected with MERS-CoV showed a gradual increase to June 2013, followed by a decline, possibly due to increased surveillance and infection control measures combined with a basic reproduction number (R0) for the virus that is less than 1. IMPORTANCE MERS-CoV adaptation toward higher rates of sustained human-to-human transmission appears not to have occurred yet. While MERS-CoV transmission currently appears weak, careful monitoring of changes in MERS-CoV genomes and of the MERS epidemic should be maintained. The observation of phylogenetically related MERS-CoV in geographically diverse locations must be taken into account in efforts to identify the animal source and transmission of the virus. MERS-CoV adaptation toward higher rates of sustained human-to-human transmission appears not to have occurred yet. While MERS-CoV transmission currently appears weak, careful monitoring of changes in MERS-CoV genomes and of the MERS epidemic should be maintained. The observation of phylogenetically related MERS-CoV in geographically diverse locations must be taken into account in efforts to identify the animal source and transmission of the virus.

[1]  X. L. Liu,et al.  Isolation and Characterization of Viruses Related to the SARS Coronavirus from Animals in Southern China , 2003, Science.

[2]  M. Suchard,et al.  Bayesian Phylogenetics with BEAUti and the BEAST 1.7 , 2012, Molecular biology and evolution.

[3]  Victor M Corman,et al.  Clinical features and virological analysis of a case of Middle East respiratory syndrome coronavirus infection , 2013, The Lancet Infectious Diseases.

[4]  A. Osterhaus,et al.  Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. , 2012, The New England journal of medicine.

[5]  Julien Riou,et al.  Interhuman transmissibility of Middle East respiratory syndrome coronavirus: estimation of pandemic risk , 2013, The Lancet.

[6]  Andrew Rambaut,et al.  Reconstructing the initial global spread of a human influenza pandemic , 2017 .

[7]  Z. Memish,et al.  Recovery from severe novel coronavirus infection. , 2012, Saudi medical journal.

[8]  Sergei L. Kosakovsky Pond,et al.  FUBAR: a fast, unconstrained bayesian approximation for inferring selection. , 2013, Molecular biology and evolution.

[9]  Linqi Zhang,et al.  Structure of MERS-CoV spike receptor-binding domain complexed with human receptor DPP4 , 2013, Cell Research.

[10]  John P. Huelsenbeck,et al.  MRBAYES: Bayesian inference of phylogenetic trees , 2001, Bioinform..

[11]  Robert C. Edgar,et al.  MUSCLE: a multiple sequence alignment method with reduced time and space complexity , 2004, BMC Bioinformatics.

[12]  Victor M Corman,et al.  Middle East respiratory syndrome coronavirus neutralising serum antibodies in dromedary camels: a comparative serological study , 2013, The Lancet Infectious Diseases.

[13]  R. Baric,et al.  Mechanisms of Zoonotic Severe Acute Respiratory Syndrome Coronavirus Host Range Expansion in Human Airway Epithelium , 2007, Journal of Virology.

[14]  Sergey I. Nikolenko,et al.  SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing , 2012, J. Comput. Biol..

[15]  Ralph S. Baric,et al.  Recombination, Reservoirs, and the Modular Spike: Mechanisms of Coronavirus Cross-Species Transmission , 2009, Journal of Virology.

[16]  Sergei L. Kosakovsky Pond,et al.  Detecting Individual Sites Subject to Episodic Diversifying Selection , 2012, PLoS genetics.

[17]  Marion P G Koopmans,et al.  Middle East respiratory syndrome coronavirus in dromedary camels: an outbreak investigation , 2013, The Lancet Infectious Diseases.

[18]  M. Suchard,et al.  Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. , 2012, Molecular biology and evolution.

[19]  R. Johnston,et al.  Severe Acute Respiratory Syndrome Coronavirus Papain-Like Protease Ubiquitin-Like Domain and Catalytic Domain Regulate Antagonism of IRF3 and NF-κB Signaling , 2009, Journal of Virology.

[20]  Zhongbin Chen,et al.  Regulation of IRF-3-dependent Innate Immunity by the Papain-like Protease Domain of the Severe Acute Respiratory Syndrome Coronavirus , 2007, Journal of Biological Chemistry.

[21]  R. Bruzzone,et al.  Ezrin Interacts with the SARS Coronavirus Spike Protein and Restrains Infection at the Entry Stage , 2012, PloS one.

[22]  Amit Kapoor,et al.  Middle East Respiratory Syndrome Coronavirus in Bats, Saudi Arabia , 2013, Emerging infectious diseases.

[23]  Z. Memish,et al.  Investigation of Anti–Middle East Respiratory Syndrome Antibodies in Blood Donors and Slaughterhouse Workers in Jeddah and Makkah, Saudi Arabia, Fall 2012 , 2013, The Journal of infectious diseases.

[24]  George F. Gao,et al.  Structure of the Fusion Core and Inhibition of Fusion by a Heptad Repeat Peptide Derived from the S Protein of Middle East Respiratory Syndrome Coronavirus , 2013, Journal of Virology.

[25]  Y. Guan,et al.  Seroepidemiology for MERS coronavirus using microneutralisation and pseudoparticle virus neutralisation assays reveal a high prevalence of antibody in dromedary camels in Egypt, June 2013. , 2013, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[26]  Ziad A Memish,et al.  Family cluster of Middle East respiratory syndrome coronavirus infections. , 2013, The New England journal of medicine.

[27]  B. Bosch,et al.  The Coronavirus Spike Protein Is a Class I Virus Fusion Protein: Structural and Functional Characterization of the Fusion Core Complex , 2003, Journal of Virology.

[28]  D. Cummings,et al.  Hospital outbreak of Middle East respiratory syndrome coronavirus. , 2013, The New England journal of medicine.

[29]  Alexander E. Gorbalenya,et al.  Genomic Characterization of a Newly Discovered Coronavirus Associated with Acute Respiratory Distress Syndrome in Humans , 2012, mBio.

[30]  Y. Guan,et al.  Full-Genome Deep Sequencing and Phylogenetic Analysis of Novel Human Betacoronavirus , 2013, Emerging infectious diseases.

[31]  S. Perlman,et al.  Coronaviruses post-SARS: update on replication and pathogenesis , 2009, Nature Reviews Microbiology.

[32]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[33]  Astrid Gall,et al.  Transmission and evolution of the Middle East respiratory syndrome coronavirus in Saudi Arabia: a descriptive genomic study , 2013, The Lancet.

[34]  Sergei L. Kosakovsky Pond,et al.  HyPhy: hypothesis testing using phylogenies , 2005, Bioinform..

[35]  Andrew Rambaut,et al.  Reconstructing the initial global spread of a human influenza pandemic: A Bayesian spatial-temporal model for the global spread of H1N1pdm. , 2009, PLoS currents.

[36]  Chengsheng Zhang,et al.  Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2 , 2005, The EMBO journal.

[37]  R. Baric,et al.  Amino Acid Substitutions in the S2 Subunit of Mouse Hepatitis Virus Variant V51 Encode Determinants of Host Range Expansion , 2007, Journal of Virology.

[38]  A. Fontanet,et al.  Clinical features and viral diagnosis of two cases of infection with Middle East Respiratory Syndrome coronavirus: a report of nosocomial transmission , 2013, The Lancet.

[39]  P. Kellam,et al.  Viral population analysis and minority-variant detection using short read next-generation sequencing , 2013, Philosophical Transactions of the Royal Society B: Biological Sciences.

[40]  Yi Shi,et al.  Molecular basis of binding between novel human coronavirus MERS-CoV and its receptor CD26 , 2013, Nature.

[41]  Simon Cauchemez,et al.  Edinburgh Research Explorer Middle East respiratory syndrome coronavirus: quantification of the extent of the epidemic, surveillance biases, and transmissibility , 2022 .

[42]  M. Suchard,et al.  Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics. , 2008, Molecular biology and evolution.