A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation
Identifieur interne : 000978 ( Pmc/Corpus ); précédent : 000977; suivant : 000979A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation
Auteurs : Tridip Sardar ; Indrajit Ghosh ; Xavier Rod ; Joydev ChattopadhyaySource :
- PLoS Neglected Tropical Diseases [ 1935-2727 ] ; 2020.
Abstract
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe acute respiratory illness with a case fatality rate (CFR) of 35,5%. The highest number of MERS-CoV cases are from Saudi-Arabia, the major worldwide hotspot for this disease. In the absence of neither effective treatment nor a ready-to-use vaccine and with yet an incomplete understanding of its epidemiological cycle, prevention and containment measures can be derived from mathematical models of disease epidemiology. We constructed 2-strain models to predict past outbreaks in the interval 2012–2016 and derive key epidemiological information for Macca, Madina and Riyadh. We approached variability in infection through three different disease incidence functions capturing social behavior in response to an epidemic (e.g. Bilinear, BL; Non-monotone, NM; and Saturated, SAT models). The best model combination successfully anticipated the total number of MERS-CoV clinical cases for the 2015–2016 season and accurately predicted both the number of cases at the peak of seasonal incidence and the overall shape of the epidemic cycle. The evolution in the basic reproduction number (
Url:
DOI: 10.1371/journal.pntd.0008065
PubMed: 32059047
PubMed Central: 7046297
Links to Exploration step
PMC:7046297Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation</title><author><name sortKey="Sardar, Tridip" sort="Sardar, Tridip" uniqKey="Sardar T" first="Tridip" last="Sardar">Tridip Sardar</name><affiliation><nlm:aff id="aff001"><addr-line>Department of Mathematics, Dinabandhu Andrews College, Kolkata, India</addr-line></nlm:aff></affiliation></author><author><name sortKey="Ghosh, Indrajit" sort="Ghosh, Indrajit" uniqKey="Ghosh I" first="Indrajit" last="Ghosh">Indrajit Ghosh</name><affiliation><nlm:aff id="aff002"><addr-line>Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, India</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rod, Xavier" sort="Rod, Xavier" uniqKey="Rod X" first="Xavier" last="Rod">Xavier Rod</name><affiliation><nlm:aff id="aff003"><addr-line>ICREA &CLIMA (Climate and Health Program), ISGlobal, Barcelona Institute for Global Health, Barcelona, Spain</addr-line></nlm:aff></affiliation></author><author><name sortKey="Chattopadhyay, Joydev" sort="Chattopadhyay, Joydev" uniqKey="Chattopadhyay J" first="Joydev" last="Chattopadhyay">Joydev Chattopadhyay</name><affiliation><nlm:aff id="aff002"><addr-line>Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, India</addr-line></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">32059047</idno><idno type="pmc">7046297</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC7046297</idno><idno type="RBID">PMC:7046297</idno><idno type="doi">10.1371/journal.pntd.0008065</idno><date when="2020">2020</date><idno type="wicri:Area/Pmc/Corpus">000978</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000978</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation</title><author><name sortKey="Sardar, Tridip" sort="Sardar, Tridip" uniqKey="Sardar T" first="Tridip" last="Sardar">Tridip Sardar</name><affiliation><nlm:aff id="aff001"><addr-line>Department of Mathematics, Dinabandhu Andrews College, Kolkata, India</addr-line></nlm:aff></affiliation></author><author><name sortKey="Ghosh, Indrajit" sort="Ghosh, Indrajit" uniqKey="Ghosh I" first="Indrajit" last="Ghosh">Indrajit Ghosh</name><affiliation><nlm:aff id="aff002"><addr-line>Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, India</addr-line></nlm:aff></affiliation></author><author><name sortKey="Rod, Xavier" sort="Rod, Xavier" uniqKey="Rod X" first="Xavier" last="Rod">Xavier Rod</name><affiliation><nlm:aff id="aff003"><addr-line>ICREA &CLIMA (Climate and Health Program), ISGlobal, Barcelona Institute for Global Health, Barcelona, Spain</addr-line></nlm:aff></affiliation></author><author><name sortKey="Chattopadhyay, Joydev" sort="Chattopadhyay, Joydev" uniqKey="Chattopadhyay J" first="Joydev" last="Chattopadhyay">Joydev Chattopadhyay</name><affiliation><nlm:aff id="aff002"><addr-line>Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, India</addr-line></nlm:aff></affiliation></author></analytic><series><title level="j">PLoS Neglected Tropical Diseases</title><idno type="ISSN">1935-2727</idno><idno type="eISSN">1935-2735</idno><imprint><date when="2020">2020</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe acute respiratory illness with a case fatality rate (CFR) of 35,5%. The highest number of MERS-CoV cases are from Saudi-Arabia, the major worldwide hotspot for this disease. In the absence of neither effective treatment nor a ready-to-use vaccine and with yet an incomplete understanding of its epidemiological cycle, prevention and containment measures can be derived from mathematical models of disease epidemiology. We constructed 2-strain models to predict past outbreaks in the interval 2012–2016 and derive key epidemiological information for Macca, Madina and Riyadh. We approached variability in infection through three different disease incidence functions capturing social behavior in response to an epidemic (e.g. Bilinear, BL; Non-monotone, NM; and Saturated, SAT models). The best model combination successfully anticipated the total number of MERS-CoV clinical cases for the 2015–2016 season and accurately predicted both the number of cases at the peak of seasonal incidence and the overall shape of the epidemic cycle. The evolution in the basic reproduction number (<italic>R</italic><sub><italic>0</italic></sub>) warns that MERS-CoV may easily take an epidemic form. The best model correctly captures this feature, indicating a high epidemic risk (1≤<italic>R</italic><sub><italic>0</italic></sub>≤2,5) in Riyadh and Macca and confirming the alleged co-circulation of more than one strain. Accurate predictions of the future MERS-CoV peak week, as well as the number of cases at the peak are now possible. These results indicate public health agencies should be aware that measures for strict containment are urgently needed before new epidemics take off in the region.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Sabir, J S M" uniqKey="Sabir J">J.S.M. Sabir</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G. Chowell</name></author><author><name sortKey="Blumberg, S" uniqKey="Blumberg S">S. Blumberg</name></author><author><name sortKey="Simonsen, L" uniqKey="Simonsen L">L. Simonsen</name></author><author><name sortKey="Miller, M A" uniqKey="Miller M">M.A. Miller</name></author><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zaki, A M" uniqKey="Zaki A">A.M. Zaki</name></author><author><name sortKey="Van Boheemen, S" uniqKey="Van Boheemen S">S. Van Boheemen</name></author><author><name sortKey="Bestebroer, T M" uniqKey="Bestebroer T">T.M. Bestebroer</name></author><author><name sortKey="Osterhaus, A D" uniqKey="Osterhaus A">A.D. Osterhaus</name></author><author><name sortKey="Fouchier, R A" uniqKey="Fouchier R">R.A Fouchier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lee, S S" uniqKey="Lee S">S.S. Lee</name></author><author><name sortKey="Wong, N S" uniqKey="Wong N">N.S. Wong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Omrani, A S" uniqKey="Omrani A">A.S. Omrani</name></author><author><name sortKey="Al Tawfiq, J A" uniqKey="Al Tawfiq J">J.A. Al-Tawfiq</name></author><author><name sortKey="Memish, Z A" uniqKey="Memish Z">Z.A Memish</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hemida, M G" uniqKey="Hemida M">M.G. Hemida</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cotten, M" uniqKey="Cotten M">M. Cotten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lingshu, W" uniqKey="Lingshu W">W. Lingshu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Myers, M F" uniqKey="Myers M">M.F. Myers</name></author><author><name sortKey="Rogers, D J" uniqKey="Rogers D">D.J. Rogers</name></author><author><name sortKey="Cox, J" uniqKey="Cox J">J. Cox</name></author><author><name sortKey="Flahault, A" uniqKey="Flahault A">A. Flahault</name></author><author><name sortKey="Hay, S I" uniqKey="Hay S">S.I. Hay</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shaman, J" uniqKey="Shaman J">J. Shaman</name></author><author><name sortKey="Karspeck, A" uniqKey="Karspeck A">A. Karspeck</name></author><author><name sortKey="Yang, W" uniqKey="Yang W">W. Yang</name></author><author><name sortKey="Tamerius, J" uniqKey="Tamerius J">J. Tamerius</name></author><author><name sortKey="Lipsitch, M" uniqKey="Lipsitch M">M Lipsitch</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Biggerstaff, M" uniqKey="Biggerstaff M">M. Biggerstaff</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shaman, J" uniqKey="Shaman J">J. Shaman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yamana, T K" uniqKey="Yamana T">T.K. Yamana</name></author><author><name sortKey="Kandula, S" uniqKey="Kandula S">S. Kandula</name></author><author><name sortKey="Shaman, J" uniqKey="Shaman J">J. Shaman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ray, E L" uniqKey="Ray E">E.L. Ray</name></author><author><name sortKey="Reich, N G" uniqKey="Reich N">N.G. Reich</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Breban, R" uniqKey="Breban R">R. Breban</name></author><author><name sortKey="Riou, J" uniqKey="Riou J">J. Riou</name></author><author><name sortKey="Fontanet, A" uniqKey="Fontanet A">A Fontanet</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cauchemez, S" uniqKey="Cauchemez S">S. Cauchemez</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zumla, A" uniqKey="Zumla A">A. Zumla</name></author><author><name sortKey="Hui, D S" uniqKey="Hui D">D.S. Hui</name></author><author><name sortKey="Perlman, S" uniqKey="Perlman S">S. Perlman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muth, D" uniqKey="Muth D">D. Muth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Capasso, V" uniqKey="Capasso V">V. Capasso</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, W" uniqKey="Wang W">W. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xiao, D" uniqKey="Xiao D">D. Xiao</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Lessler, J" uniqKey="Lessler J">J. Lessler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Assiri, A" uniqKey="Assiri A">A. Assiri</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xia, Z Q" uniqKey="Xia Z">Z.Q. Xia</name></author><author><name sortKey="Zhang, J" uniqKey="Zhang J">J. Zhang</name></author><author><name sortKey="Xue, Y K" uniqKey="Xue Y">Y.K. Xue</name></author><author><name sortKey="Sun, G Q" uniqKey="Sun G">G.Q. Sun</name></author><author><name sortKey="Jin, Z" uniqKey="Jin Z">Z Jin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Majumder, M S" uniqKey="Majumder M">M.S. Majumder</name></author><author><name sortKey="Rivers, C" uniqKey="Rivers C">C. Rivers</name></author><author><name sortKey="Lofgren, E" uniqKey="Lofgren E">E. Lofgren</name></author><author><name sortKey="Fisman, D" uniqKey="Fisman D">D Fisman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Poletto, C" uniqKey="Poletto C">C. Poletto</name></author><author><name sortKey="Colizza, V" uniqKey="Colizza V">V. Colizza</name></author><author><name sortKey="Boelle, P Y" uniqKey="Boelle P">P.Y Boëlle</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H. Nishiura</name></author><author><name sortKey="Miyamatsu, Y" uniqKey="Miyamatsu Y">Y. Miyamatsu</name></author><author><name sortKey="Chowell, G" uniqKey="Chowell G">G. Chowell</name></author><author><name sortKey="Saitoh, M" uniqKey="Saitoh M">M Saitoh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Poletto, C" uniqKey="Poletto C">C. Poletto</name></author><author><name sortKey="Boelle, P Y" uniqKey="Boelle P">P.Y. Boëlle</name></author><author><name sortKey="Colizza, V" uniqKey="Colizza V">V Colizza</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hsieh, Y H" uniqKey="Hsieh Y">Y.H. Hsieh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H. Nishiura</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chan, R W" uniqKey="Chan R">R.W. Chan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kucharski, A J" uniqKey="Kucharski A">A.J. Kucharski</name></author><author><name sortKey="Althaus, C" uniqKey="Althaus C">C Althaus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chowell, G" uniqKey="Chowell G">G. Chowell</name></author><author><name sortKey="Abdirizak, F" uniqKey="Abdirizak F">F. Abdirizak</name></author><author><name sortKey="Lee, S" uniqKey="Lee S">S. Lee</name></author><author><name sortKey="Lee, J" uniqKey="Lee J">J. Lee</name></author><author><name sortKey="Jung, E" uniqKey="Jung E">E. Jung</name></author><author><name sortKey="Nishiura, H" uniqKey="Nishiura H">H. Nishiura</name></author><author><name sortKey="Viboud, C" uniqKey="Viboud C">C Viboud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cotten, M" uniqKey="Cotten M">M. Cotten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Waters, E K" uniqKey="Waters E">E.K. Waters</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lloyd, M" uniqKey="Lloyd M">M. Lloyd</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Mukandavire, Z" uniqKey="Mukandavire Z">Z. Mukandavire</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Johnson, L R" uniqKey="Johnson L">L.R. Johnson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Haario, H" uniqKey="Haario H">H. Haario</name></author><author><name sortKey="Laine, M" uniqKey="Laine M">M. Laine</name></author><author><name sortKey="Mira, A" uniqKey="Mira A">A. Mira</name></author><author><name sortKey="Saksman, E" uniqKey="Saksman E">E Saksman</name></author></analytic></biblStruct><biblStruct></biblStruct><biblStruct></biblStruct><biblStruct><analytic><author><name sortKey="Lin, Q" uniqKey="Lin Q">Q. Lin</name></author><author><name sortKey="Chiu, A P Y" uniqKey="Chiu A">A.P.Y Chiu</name></author><author><name sortKey="Zhao, S" uniqKey="Zhao S">S Zhao</name></author><author><name sortKey="He, D" uniqKey="He D">D He</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS Negl Trop Dis</journal-id><journal-id journal-id-type="iso-abbrev">PLoS Negl Trop Dis</journal-id><journal-id journal-id-type="publisher-id">plos</journal-id><journal-id journal-id-type="pmc">plosntds</journal-id><journal-title-group><journal-title>PLoS Neglected Tropical Diseases</journal-title></journal-title-group><issn pub-type="ppub">1935-2727</issn><issn pub-type="epub">1935-2735</issn><publisher><publisher-name>Public Library of Science</publisher-name><publisher-loc>San Francisco, CA USA</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">32059047</article-id><article-id pub-id-type="pmc">7046297</article-id><article-id pub-id-type="doi">10.1371/journal.pntd.0008065</article-id><article-id pub-id-type="publisher-id">PNTD-D-18-01971</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Viruses</subject><subj-group><subject>RNA viruses</subject><subj-group><subject>Coronaviruses</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Microbiology</subject><subj-group><subject>Medical Microbiology</subject><subj-group><subject>Microbial Pathogens</subject><subj-group><subject>Viral Pathogens</subject><subj-group><subject>Coronaviruses</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Pathology and Laboratory Medicine</subject><subj-group><subject>Pathogens</subject><subj-group><subject>Microbial Pathogens</subject><subj-group><subject>Viral Pathogens</subject><subj-group><subject>Coronaviruses</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Viruses</subject><subj-group><subject>Viral Pathogens</subject><subj-group><subject>Coronaviruses</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Epidemiology</subject><subj-group><subject>Infectious Disease Epidemiology</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Infectious Diseases</subject><subj-group><subject>Infectious Disease Epidemiology</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Epidemiology</subject></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>People and places</subject><subj-group><subject>Geographical locations</subject><subj-group><subject>Asia</subject><subj-group><subject>Saudi Arabia</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and Life Sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Eukaryota</subject><subj-group><subject>Animals</subject><subj-group><subject>Vertebrates</subject><subj-group><subject>Amniotes</subject><subj-group><subject>Mammals</subject><subj-group><subject>Camels</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Pulmonology</subject><subj-group><subject>Respiratory Infections</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Medicine and Health Sciences</subject><subj-group><subject>Infectious Diseases</subject><subj-group><subject>Zoonoses</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Research and Analysis Methods</subject><subj-group><subject>Mathematical and Statistical Techniques</subject><subj-group><subject>Mathematical Models</subject></subj-group></subj-group></subj-group></article-categories><title-group><article-title>A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation</article-title><alt-title alt-title-type="running-head">Seasonal forecasting of MERS-CoV epidemics</alt-title></title-group><contrib-group><contrib contrib-type="author" equal-contrib="yes"><name><surname>Sardar</surname><given-names>Tridip</given-names></name><role content-type="http://credit.casrai.org/">Conceptualization</role><role content-type="http://credit.casrai.org/">Data curation</role><role content-type="http://credit.casrai.org/">Formal analysis</role><role content-type="http://credit.casrai.org/">Investigation</role><role content-type="http://credit.casrai.org/">Methodology</role><role content-type="http://credit.casrai.org/">Validation</role><role content-type="http://credit.casrai.org/">Writing – original draft</role><role content-type="http://credit.casrai.org/">Writing – review & editing</role><xref ref-type="aff" rid="aff001"><sup>1</sup></xref></contrib><contrib contrib-type="author" equal-contrib="yes"><contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0002-0492-2948</contrib-id><name><surname>Ghosh</surname><given-names>Indrajit</given-names></name><role content-type="http://credit.casrai.org/">Data curation</role><role content-type="http://credit.casrai.org/">Formal analysis</role><role content-type="http://credit.casrai.org/">Methodology</role><role content-type="http://credit.casrai.org/">Software</role><role content-type="http://credit.casrai.org/">Validation</role><role content-type="http://credit.casrai.org/">Writing – review & editing</role><xref ref-type="aff" rid="aff002"><sup>2</sup></xref></contrib><contrib contrib-type="author"><contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0003-4843-6180</contrib-id><name><surname>Rodó</surname><given-names>Xavier</given-names></name><role content-type="http://credit.casrai.org/">Conceptualization</role><role content-type="http://credit.casrai.org/">Formal analysis</role><role content-type="http://credit.casrai.org/">Investigation</role><role content-type="http://credit.casrai.org/">Methodology</role><role content-type="http://credit.casrai.org/">Supervision</role><role content-type="http://credit.casrai.org/">Validation</role><role content-type="http://credit.casrai.org/">Writing – original draft</role><role content-type="http://credit.casrai.org/">Writing – review & editing</role><xref ref-type="aff" rid="aff003"><sup>3</sup></xref><xref ref-type="corresp" rid="cor001">*</xref></contrib><contrib contrib-type="author"><name><surname>Chattopadhyay</surname><given-names>Joydev</given-names></name><role content-type="http://credit.casrai.org/">Conceptualization</role><role content-type="http://credit.casrai.org/">Data curation</role><role content-type="http://credit.casrai.org/">Investigation</role><role content-type="http://credit.casrai.org/">Methodology</role><role content-type="http://credit.casrai.org/">Writing – original draft</role><role content-type="http://credit.casrai.org/">Writing – review & editing</role><xref ref-type="aff" rid="aff002"><sup>2</sup></xref></contrib></contrib-group><aff id="aff001"><label>1</label><addr-line>Department of Mathematics, Dinabandhu Andrews College, Kolkata, India</addr-line></aff><aff id="aff002"><label>2</label><addr-line>Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, India</addr-line></aff><aff id="aff003"><label>3</label><addr-line>ICREA &CLIMA (Climate and Health Program), ISGlobal, Barcelona Institute for Global Health, Barcelona, Spain</addr-line></aff><contrib-group><contrib contrib-type="editor"><name><surname>Althouse</surname><given-names>Benjamin</given-names></name><role>Editor</role><xref ref-type="aff" rid="edit1"></xref></contrib></contrib-group><aff id="edit1"><addr-line>Institute for Disease Modeling, UNITED STATES</addr-line></aff><author-notes><fn fn-type="COI-statement" id="coi001"><p>No authors have competing interests.</p></fn><corresp id="cor001">* E-mail: <email>xavier.rodo@isglobal.org</email></corresp></author-notes><pub-date pub-type="epub"><day>14</day><month>2</month><year>2020</year></pub-date><pub-date pub-type="collection"><month>2</month><year>2020</year></pub-date><volume>14</volume><issue>2</issue><elocation-id>e0008065</elocation-id><history><date date-type="received"><day>14</day><month>12</month><year>2018</year></date><date date-type="accepted"><day>15</day><month>1</month><year>2020</year></date></history><permissions><copyright-statement>© 2020 Sardar et al</copyright-statement><copyright-year>2020</copyright-year><copyright-holder>Sardar et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="pntd.0008065.pdf"></self-uri><abstract><p>Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe acute respiratory illness with a case fatality rate (CFR) of 35,5%. The highest number of MERS-CoV cases are from Saudi-Arabia, the major worldwide hotspot for this disease. In the absence of neither effective treatment nor a ready-to-use vaccine and with yet an incomplete understanding of its epidemiological cycle, prevention and containment measures can be derived from mathematical models of disease epidemiology. We constructed 2-strain models to predict past outbreaks in the interval 2012–2016 and derive key epidemiological information for Macca, Madina and Riyadh. We approached variability in infection through three different disease incidence functions capturing social behavior in response to an epidemic (e.g. Bilinear, BL; Non-monotone, NM; and Saturated, SAT models). The best model combination successfully anticipated the total number of MERS-CoV clinical cases for the 2015–2016 season and accurately predicted both the number of cases at the peak of seasonal incidence and the overall shape of the epidemic cycle. The evolution in the basic reproduction number (<italic>R</italic><sub><italic>0</italic></sub>) warns that MERS-CoV may easily take an epidemic form. The best model correctly captures this feature, indicating a high epidemic risk (1≤<italic>R</italic><sub><italic>0</italic></sub>≤2,5) in Riyadh and Macca and confirming the alleged co-circulation of more than one strain. Accurate predictions of the future MERS-CoV peak week, as well as the number of cases at the peak are now possible. These results indicate public health agencies should be aware that measures for strict containment are urgently needed before new epidemics take off in the region.</p></abstract><abstract abstract-type="summary"><title>Author summary</title><p>There is currently no way to anticipate MERS-CoV epidemic outbreaks and strategies for disease prediction and containment are largely undermined by the limited knowledge of its epidemiological cycle. Not an effective treatment nor a vaccine for MERS-CoV exist to date. Instead, using three two-strain mathematical models that incorporate human social behavior as different disease incidence functions (e.g. bilinear, non-monotone and saturated), the best model combinations successfully anticipate the occurrence of the peak week in the season and the incidence at the peak. Our results confirm there are currently 2 strains co-circulating in the most populated regions in Saudi Arabia and highlight the high risk for large epidemic outbreaks, while the role of super-spreaders appears irrelevant for disease spread.</p></abstract><funding-group><award-group id="award001"><funding-source><institution-wrap><institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100002809</institution-id><institution>Generalitat de Catalunya</institution></institution-wrap></funding-source><award-id>SLT002/16/00466</award-id><principal-award-recipient><contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0003-4843-6180</contrib-id><name><surname>Rodó</surname><given-names>Xavier</given-names></name></principal-award-recipient></award-group><award-group id="award002"><funding-source><institution-wrap><institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100010426</institution-id><institution>UGC-DAE Consortium for Scientific Research, University Grants Commission</institution></institution-wrap></funding-source><award-id>UGC-JRF23402016</award-id><principal-award-recipient><contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0002-0492-2948</contrib-id><name><surname>Ghosh</surname><given-names>Indrajit</given-names></name></principal-award-recipient></award-group><funding-statement>I.G. was supported by the research fellowship from University Grants Commission (<ext-link ext-link-type="uri" xlink:href="https://www.ugc.ac.in/">https://www.ugc.ac.in/</ext-link>), Government of India. X.R. acknowledges the support of a fellowship by the PERIS PICAT project SLT002/16/00466 (<email>peris@gencat.cat</email>| <ext-link ext-link-type="uri" xlink:href="http://canalsalut.gencat.cat">http://canalsalut.gencat.cat</ext-link>) from the Catalan Ministry of Health. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement></funding-group><counts><fig-count count="8"></fig-count><table-count count="3"></table-count><page-count count="20"></page-count></counts><custom-meta-group><custom-meta><meta-name>PLOS Publication Stage</meta-name><meta-value>vor-update-to-uncorrected-proof</meta-value></custom-meta><custom-meta><meta-name>Publication Update</meta-name><meta-value>2020-02-27</meta-value></custom-meta><custom-meta id="data-availability"><meta-name>Data Availability</meta-name><meta-value>The website containing all MERS-CoV dataset used is the following: <ext-link ext-link-type="uri" xlink:href="http://empres-i.fao.org/eipws3g/">http://empres-i.fao.org/eipws3g/</ext-link></meta-value></custom-meta></custom-meta-group></article-meta><notes><title>Data Availability</title><p>The website containing all MERS-CoV dataset used is the following: <ext-link ext-link-type="uri" xlink:href="http://empres-i.fao.org/eipws3g/">http://empres-i.fao.org/eipws3g/</ext-link></p></notes></front><body><sec sec-type="intro" id="sec001"><title>Introduction</title><p>The first case of Middle East Respiratory Syndrome Coronavirus (MERS-CoV) infection was identified in Saudi Arabia in 2012[<xref rid="pntd.0008065.ref001" ref-type="bibr">1</xref>–<xref rid="pntd.0008065.ref005" ref-type="bibr">5</xref>]. Since then the country suffers from repeated outbreaks of MERS-CoV in different provinces (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1A and 1B</xref>)[<xref rid="pntd.0008065.ref002" ref-type="bibr">2</xref>]. There are two main routes of MERS-CoV transmission: animal-to-human and human-to-human[<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>–<xref rid="pntd.0008065.ref006" ref-type="bibr">6</xref>]. It is suspected that dromedary camels are the source of human infections[<xref rid="pntd.0008065.ref001" ref-type="bibr">1</xref>] but the transmission route of MERS-CoV to humans is yet not well understood (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1C</xref>). However, transmission from camels to humans is confirmed from the isolation of near-identical strains of MERS-CoV from epidemiologically coexisting camels and humans[<xref rid="pntd.0008065.ref007" ref-type="bibr">7</xref>]. The patients might be exposed to MERS-CoV by consumption of raw camel products[<xref rid="pntd.0008065.ref008" ref-type="bibr">8</xref>] (e.g. milk and dairy products, raw meat, etc.). Meanwhile, human-to-human transmission has been reported in society and hospital settings[<xref rid="pntd.0008065.ref009" ref-type="bibr">9</xref>–<xref rid="pntd.0008065.ref036" ref-type="bibr">36</xref>] with the virus being transmitted between humans during close human-to-human contact through droplets of respiratory secretions. Potential propagation to nearby and more distant regions is also a high-risk possibility as an outbreak of MERS-CoV is likely to emerge in areas such as nearby countries in the Middle East and eastern Africa where the camel trade connects the different regions (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1</xref>). Unfortunately, as for now, MERS-CoV vaccines are only at the preclinical phase[<xref rid="pntd.0008065.ref010" ref-type="bibr">10</xref>], increasing our understanding of its epidemic potential and knowledge on the drivers of MERS-CoV variability might help to achieve better preparedness ahead of forthcoming epidemics.</p><fig id="pntd.0008065.g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g001</object-id><label>Fig 1</label><caption><p>Distribution of MERS-CoV weekly cases in Riyadh, Macca and Madina provinces during July, 2013—June, 2016. (a) The geographic distribution of MERS-CoV cases reported in 2017 in the Middle East. (b) Time series of weekly incidence data of MERS-CoV in three major provinces, Riyadh, Maccaand Madina, respectively. (c) Two strain Models flow diagram considering asymptomatic, symptomatic, hospitalized, and Zoonotic transmission. Open source KML map for the Middle East was kindly obtained and redrawn by Josep-BoyardMicheau from <ext-link ext-link-type="uri" xlink:href="https://community.qlik.com/t5/Qlik-Sense-Enterprise-Documents/GCC-Middle-East-country-boundary-KML-maps-KML-Shapefile/ta-p/1478595">https://community.qlik.com/t5/Qlik-Sense-Enterprise-Documents/GCC-Middle-East-country-boundary-KML-maps-KML-Shapefile/ta-p/1478595</ext-link>.</p></caption><graphic xlink:href="pntd.0008065.g001"></graphic></fig><p>The ability to predict disease outbreaks will provide a mechanism for policymakers and health-care services to respond to epidemics in a timely manner, reducing the impact and maximizing the limited resources available to be deployed[<xref rid="pntd.0008065.ref011" ref-type="bibr">11</xref>]. The timing and severity of infectious disease outbreaks, two matters of considerable public-health relevance, are the main challenges when attempting to predict disease outbreaks[<xref rid="pntd.0008065.ref011" ref-type="bibr">11</xref>–<xref rid="pntd.0008065.ref014" ref-type="bibr">14</xref>]. Attempts to set up prediction frameworks for anticipating epidemics for other diseases such as dengue and influenza were pursued in the recent past with different degree of success, but clear added value[<xref rid="pntd.0008065.ref011" ref-type="bibr">11</xref>–<xref rid="pntd.0008065.ref016" ref-type="bibr">16</xref>]. These systems proved effective to better anticipate future outbreaks and increase our understanding on mechanisms driving disease variability. To this end, a threshold quantity that is most important to anticipate the risk of future outbreaks is the basic reproduction number (<italic>R</italic><sub><italic>0</italic></sub>). There are quite a few studies that estimated <italic>R</italic><sub><italic>0</italic></sub> for the current MERS-CoV outbreak. Majumdar et al.[<xref rid="pntd.0008065.ref028" ref-type="bibr">28</xref>] for instance, estimated <italic>R</italic><sub><italic>0</italic></sub> for Jeddah and Riyadh using a function called Incidence Decay with Exponential Adjustments (IDEA). They found that the estimate of <italic>R</italic><sub><italic>0</italic></sub> in Jeddah and Riyadh are in the ranges (3.5–6.7) and (2.0–2.8), respectively. A stochastic epidemiological model of MERS-CoV with zoonotic and human-to-human transmission was considered by Poletto et al.[<xref rid="pntd.0008065.ref029" ref-type="bibr">29</xref>] to quantify the rates of generation of cases from those two transmission routes. They found that spring 2014 cases led to the increase in transmission rates, which brings <italic>R</italic><sub><italic>0</italic></sub> to values above unity. Heishet al.[<xref rid="pntd.0008065.ref032" ref-type="bibr">32</xref>] used a simple mathematical model to trace the temporal course of South Korea MERS-CoV outbreak. They estimated <italic>R</italic><sub><italic>0</italic></sub> to be in the range of 7.0 to 19.3. Instead, Nishiura <italic>et</italic> al.[<xref rid="pntd.0008065.ref030" ref-type="bibr">30</xref>] estimated the expected number of secondary cases following the importation of an index case (countries other than KSA, Qatar and UAE), using a branching process type modelling approach. They also suggested that even if <italic>R</italic><sub><italic>0</italic></sub> is below unity, a large cluster of cases with multiple generations might occur. The aforementioned studies indicated that the basic reproduction number is above unity, although there are few studies that suggest <italic>R</italic><sub><italic>0</italic></sub> is less than one[<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>,<xref rid="pntd.0008065.ref017" ref-type="bibr">17</xref>,<xref rid="pntd.0008065.ref018" ref-type="bibr">18</xref>]. However, at the same time, intriguing differences in observed and predicted scenarios clearly suggested other important factors might be missing in previous models built for these MERS-CoV case studies. For instance, in a recent survey for Mers-CoV in over 1300 Saudi Arabian camels, Sabiret al.[<xref rid="pntd.0008065.ref001" ref-type="bibr">1</xref>] found that camels share three coronaviruses (CoV) species with humans. Among them, one has been dominant in the region since December 2014 and led to the human MERS-CoV outbreaks occurring in 2015. The wide species range of CoVs and their propensity to cross species boundaries suggest that more MERS-CoV uncontained outbreaks may likely occur in the future. Zumlaet al.[<xref rid="pntd.0008065.ref019" ref-type="bibr">19</xref>] also suggested that MERS-CoV species can mutate to have increased inter-human transmissibility. Cottenet al.[<xref rid="pntd.0008065.ref037" ref-type="bibr">37</xref>] used MERS data from the period May to September 2013. They found four Saudi Arabia MERS-CoV phylogenetic clades, with 3 clades apparently no longer contributing to current cases, therefore not appearing in the saliva of camels. They also show that the ancestors of most of the viral clades originated in Riyadh (See Cotten et al.[<xref rid="pntd.0008065.ref037" ref-type="bibr">37</xref>]). Recently, strain variability in MERS-CoV infection was confirmed in 51 respiratory samples from 32 patients, confirming that more than one strain of human MERS-CoVis present[<xref rid="pntd.0008065.ref020" ref-type="bibr">20</xref>]. The fact that more than one MERS-CoV strain is currently circulating in Saudi Arabia should be properly accounted for as this presence may have substantially contributed to amplify the transmission intensity producing repeated MERS-CoV changes in incidence through their periodic reintroduction into human populations[<xref rid="pntd.0008065.ref019" ref-type="bibr">19</xref>]. Furthermore, Chan et al.[<xref rid="pntd.0008065.ref034" ref-type="bibr">34</xref>] characterized MERS-CoV viruses from dromedaries in Saudi Arabia and Egypt and compared them with a human MERS-CoV reference strain. They isolated three dromedary strains, two from Saudi Arabia and one from Egypt. The human and dromedary MERS-CoV strains had similar viral replication competence in Vero-E6 cells and respiratory tropism in ex-vivo cultures of the human respiratory tract. They also suggested that dromedary viruses from Saudi Arabia and Egypt are probably infectious to human beings[<xref rid="pntd.0008065.ref034" ref-type="bibr">34</xref>].</p></sec><sec sec-type="materials|methods" id="sec002"><title>Methods</title><sec id="sec003"><title>a) Disease incidence functions</title><p>We constructed three different 2-strain models for MERS-CoV that consider community human-human, hospital human-human, and passive zoonotic transmission (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1C</xref> and see also <xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s018">S18</xref> Tables). Models derived also incorporate the effect of strain variability with strain-1 as the dominant strain (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1C</xref> and <xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s018">S18</xref> Tables). Variability among models is based upon three different disease incidence functions (e.g. BL, NM and SAT models [<xref rid="pntd.0008065.ref021" ref-type="bibr">21</xref>–<xref rid="pntd.0008065.ref023" ref-type="bibr">23</xref>] and see formulation in <xref rid="pntd.0008065.t001" ref-type="table">Table 1A</xref> and <xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref>).</p><fig id="pntd.0008065.g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g002</object-id><label>Fig 2</label><caption><p>Representation of the three incidence functions, g(I). Top: g(I) = βI, bilinear incidence function. Bottom-left: g(I) = <inline-formula id="pntd.0008065.e001"><alternatives><graphic id="pntd.0008065.e001g" xlink:href="pntd.0008065.e001"></graphic><mml:math id="M1"><mml:mfrac><mml:mrow><mml:mi>β</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>α</mml:mi><mml:msup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula>, non-monotone incidence function and; bottom-right: g(I) = <inline-formula id="pntd.0008065.e002"><alternatives><graphic id="pntd.0008065.e002g" xlink:href="pntd.0008065.e002"></graphic><mml:math id="M2"><mml:mfrac><mml:mrow><mml:mi>β</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>α</mml:mi><mml:mi>I</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula>, saturated incidence function. Here, N = 1000, α = 1 and β = 1.</p></caption><graphic xlink:href="pntd.0008065.g002"></graphic></fig><table-wrap id="pntd.0008065.t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.t001</object-id><label>Table 1</label><caption><title>(A)Description of the incidence functions for the three models considered. (B)Main parameters of the three 2-strain models.</title></caption><alternatives><graphic id="pntd.0008065.t001g" xlink:href="pntd.0008065.t001"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="center" colspan="6" rowspan="1">A</th></tr></thead><tbody><tr><td align="center" colspan="2" rowspan="1"><underline>Incidence functions</underline></td><td align="center" colspan="2" rowspan="1"><underline>Description</underline></td><td align="center" colspan="2" rowspan="1"><underline>Source</underline></td></tr><tr><td align="center" colspan="2" rowspan="1"><inline-formula id="pntd.0008065.e003"><alternatives><graphic id="pntd.0008065.e003g" xlink:href="pntd.0008065.e003"></graphic><mml:math id="M3"><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>β</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula></td><td align="center" colspan="2" rowspan="1">Bilinear incidence function</td><td align="center" colspan="2" rowspan="1">23</td></tr><tr><td align="center" colspan="2" rowspan="1"><inline-formula id="pntd.0008065.e004"><alternatives><graphic id="pntd.0008065.e004g" xlink:href="pntd.0008065.e004"></graphic><mml:math id="M4"><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>β</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>α</mml:mi><mml:msup><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula></td><td align="center" colspan="2" rowspan="1">Non-monotone incidence function</td><td align="center" colspan="2" rowspan="1">21, 22</td></tr><tr><td align="center" colspan="2" rowspan="1"><inline-formula id="pntd.0008065.e005"><alternatives><graphic id="pntd.0008065.e005g" xlink:href="pntd.0008065.e005"></graphic><mml:math id="M5"><mml:mi>f</mml:mi><mml:mo>(</mml:mo><mml:mrow><mml:mi>I</mml:mi></mml:mrow><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>β</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mi>S</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mo>(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mi>α</mml:mi><mml:mi>I</mml:mi></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula></td><td align="center" colspan="2" rowspan="1">Saturated incidence function</td><td align="center" colspan="2" rowspan="1">21</td></tr><tr><td align="center" colspan="6" rowspan="1">B</td></tr><tr><td align="center" rowspan="1" colspan="1"><underline>Symbol</underline></td><td align="center" colspan="2" rowspan="1"><underline>Description</underline></td><td align="center" colspan="2" rowspan="1"><underline>Range/Fixed</underline></td><td align="center" rowspan="1" colspan="1"><underline>Sources</underline></td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>μ</italic></td><td align="center" colspan="2" rowspan="1">Natural birth/death rate of human</td><td align="center" colspan="2" rowspan="1">2.5924×10<sup>−4</sup>Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">24</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>β</italic><sub>1</sub></td><td align="center" colspan="2" rowspan="1">Transmission coefficients of symptomatic and asymptomatic cases</td><td align="center" colspan="2" rowspan="1"><italic>[0</italic>, <italic>100]</italic> Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>θ</italic></td><td align="center" colspan="2" rowspan="1">Measure of variability of two strains</td><td align="center" colspan="2" rowspan="1"><italic>[0,1]</italic></td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>ρ</italic></td><td align="center" colspan="2" rowspan="1">Reduction in transmission due to less infectiousness of asymptomatic infected individuals</td><td align="center" colspan="2" rowspan="1"><italic>[0,1]</italic></td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>β</italic><sub>2</sub></td><td align="center" colspan="2" rowspan="1">Transmission rate from hospitalized to susceptible individuals</td><td align="center" colspan="2" rowspan="1"><italic>[0</italic>, <italic>100]</italic> Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>β</italic><sub>3</sub></td><td align="center" colspan="2" rowspan="1">Import rate of external zoonotic infection</td><td align="center" colspan="2" rowspan="1"><italic>[0, 10]</italic> Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>p</italic><sub>1</sub></td><td align="center" colspan="2" rowspan="1">Proportion of susceptible that get infected from hospitalized individuals with strain-1</td><td align="center" colspan="2" rowspan="1">[0,1]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>p</italic><sub>2</sub></td><td align="center" colspan="2" rowspan="1">Proportion of susceptible that get infected from zoonotic infection with strain-1</td><td align="center" colspan="2" rowspan="1">[0,1]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>p</italic><sub>3</sub></td><td align="center" colspan="2" rowspan="1">Proportion of exposed individuals that become symptomatic infected</td><td align="center" colspan="2" rowspan="1">0.553</td><td align="center" rowspan="1" colspan="1">4</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>c</italic><sub>1</sub></td><td align="center" colspan="2" rowspan="1">Undetected entry rate of strain-1 asymptomatic class to hospitalized class</td><td align="center" colspan="2" rowspan="1">[0,5]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>c</italic><sub>2</sub></td><td align="center" colspan="2" rowspan="1">Undetected entry rate of strain-2 asymptomatic class to hospitalized class</td><td align="center" colspan="2" rowspan="1">[0,5]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><inline-formula id="pntd.0008065.e006"><alternatives><graphic id="pntd.0008065.e006g" xlink:href="pntd.0008065.e006"></graphic><mml:math id="M6"><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>σ</mml:mi></mml:mrow></mml:mfrac></mml:math></alternatives></inline-formula></td><td align="center" colspan="2" rowspan="1">Incubation period</td><td align="center" colspan="2" rowspan="1">0.7429 Week</td><td align="center" rowspan="1" colspan="1">19,25,26</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>K</italic></td><td align="center" colspan="2" rowspan="1">Recovery rate of asymptomatic individuals</td><td align="center" colspan="2" rowspan="1">1.4 Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">27</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>λ</italic></td><td align="center" colspan="2" rowspan="1">Hospitalization rate of symptomatic individuals</td><td align="center" colspan="2" rowspan="1">4.4824 Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">4</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>γ</italic></td><td align="center" colspan="2" rowspan="1">Recovery rate of hospitalized individuals</td><td align="center" colspan="2" rowspan="1">1 Week<sup>-1</sup></td><td align="center" rowspan="1" colspan="1">4</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>δ</italic><sub>2</sub></td><td align="center" colspan="2" rowspan="1">Disease induced death rate of hospitalized individuals</td><td align="center" colspan="2" rowspan="1">Data</td><td align="center" rowspan="1" colspan="1">2</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>α</italic><sub>1</sub></td><td align="center" colspan="2" rowspan="1">Measure of inhibitory effect of strain-1 (only for non-monotone and saturated incidence models)</td><td align="center" colspan="2" rowspan="1">[0,500]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr><tr><td align="center" rowspan="1" colspan="1"><italic>α</italic><sub>2</sub></td><td align="center" colspan="2" rowspan="1">Measure of inhibitory effect of strain-2 (only for non-monotone and saturated incidence models)</td><td align="center" colspan="2" rowspan="1">[0,500]</td><td align="center" rowspan="1" colspan="1">Estimated</td></tr></tbody></table></alternatives></table-wrap><p>The top panel in <xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref> represents the BL force of infection. Here, as the number of infected individuals increases the disease transmission also increases linearly. The bottom-left panel represents the NM incidence function. Biologically this incidence function can account for “psychological effects” [<xref rid="pntd.0008065.ref038" ref-type="bibr">38</xref>–<xref rid="pntd.0008065.ref039" ref-type="bibr">39</xref>]. In the presence of psychological effects, initially the disease transmission rate increases rapidly when the number of infected individuals is small. However, this rate falls also rapidly in the presence of a large number of infected persons in the community. As in the case of MERS-CoV infections the CFR is about 40%, psychological effects for this infection represents the effect in the community of fear of becoming infected. This fear effect reduces the transmission rate rapidly in the presence of a high number of infections in the community. The bottom-right panel of <xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref> represents the SAT incidence function. Biologically this incidence function represents the crowding effect in disease transmission. This crowding effect explains how the number of new infections becomes constant as the number of infected individual increases. This effect is known to reproduce the awareness effect of the disease during the course of an epidemic.</p></sec><sec id="sec004"><title>b) Super-spreaders and 1-strain vs 2-strain models</title><p>Furthermore, we investigated the potential effects of super-spreading events by incorporating an additional compartment for super-spreaders to our 2-strain model. We also consider the effect of variable zoonotic transmission by incorporating dynamics of the camel population into our 2-strain model. Information about the calculation of the epidemiological parameter values for the newly proposed 2-strain models is provided in <xref rid="pntd.0008065.t001" ref-type="table">Table 1B</xref>. Data used corresponds to weekly cases of MERS-CoV for the three provinces in Saudi Arabia where most clinical cases for MERS-CoV occurred for the time intervals July, 6<sup>th</sup> 2013 till June, 28<sup>th</sup> 2016 (Riyadh), April, 7<sup>th</sup>2014 till June, 26<sup>th</sup>2016 (Macca), and April, 20<sup>th</sup>2014 till June, 26<sup>th</sup>2016 (Madina), respectively (<xref ref-type="fig" rid="pntd.0008065.g001">Fig 1A</xref>). Main parameters of the different models are provided in <xref rid="pntd.0008065.t001" ref-type="table">Table 1B</xref>. Parameters estimates are obtained by fitting the models to new MERS-CoV hospitalized weekly data[<xref rid="pntd.0008065.ref040" ref-type="bibr">40</xref>]. We also estimated from data itself some unknown initial conditions for the model. Delayed Rejection Adaptive Metropolis algorithm is used here to sample the 95% confidence.The goodness of fitsto compare among 2-strain and single strain models were tested by determining the respective AIC and BIC values (See <xref ref-type="supplementary-material" rid="pntd.0008065.s019">S19 Table</xref>). To further compare among predictive capabilities of these models with regard to those of their single-strain versions, predictions of raw clinical cases were attempted. Albeit some discrepancies exist, fairly accurate results were obtained despite the low numbers and highly stochastic nature of the three datasets and the nearly operational capacity of these models (<xref ref-type="fig" rid="pntd.0008065.g003">Fig 3</xref>). In fact, near-operational systems usually require a much longer training period than the scarce three years employed here. Differences arising from the comparison of simulated and observed cases may also be due to the fact that we are fitting the 2-strain and single strain models to cumulative MERS-CoV cases, rather than to new cases.<sup>.</sup>The procedure for model adjustment is as follows: for approaching near-operational predictions in each province, we followed former initiatives (e.g. on influenza[<xref rid="pntd.0008065.ref014" ref-type="bibr">14</xref>] and dengue[<xref rid="pntd.0008065.ref016" ref-type="bibr">16</xref>]) and partitioned the whole clinical cases datasets into two parts. First part was used for calibrating the models and the remaining part (52 weeks) were left aside for out-of-fit prediction. For Riyadh, Macca and Madina prediction periods were, respectively, from weeks 105–156 (7<sup>th</sup> June, 2015 - 11<sup>th</sup> June, 2016), weeks 65–116 (11<sup>th</sup> July, 2015 - 2<sup>nd</sup> July, 2016), and weeks 63–114 (11<sup>th</sup> July, 2015 - 2<sup>nd</sup> July, 2016)(<xref ref-type="fig" rid="pntd.0008065.g003">Fig 3</xref>). We generated predictions for three major characteristics of the epidemiological cycle similar to previous attempts made for cholera, dengue and influenza[<xref rid="pntd.0008065.ref015" ref-type="bibr">15</xref>,<xref rid="pntd.0008065.ref041" ref-type="bibr">41</xref>–<xref rid="pntd.0008065.ref043" ref-type="bibr">43</xref>] namely: (a) peak week: the week during which the maximum number of clinical cases occurred in a season (comprising 52 weeks); (b) peak maximum: the number of cases occurring at the peak week, and (c) season totals: the total number of cases in the entire season. Prediction for each target variable was made every 4 weeks (i.e. week 0, 4, 8, …, 48), with week 0 corresponding to the first week of the prediction interval. In the case of predictions for week <italic>N</italic>, data up to week <italic>N</italic> was used to fit the model and the trajectory predicted for out-of-fit future weeks <italic>N+1</italic> through week 52 (<xref ref-type="fig" rid="pntd.0008065.g003">Fig 3</xref>). Indeed, the best models at predicting MERS-CoV incidence differed among regions (<xref ref-type="fig" rid="pntd.0008065.g004">Fig 4</xref>), highlighting the complexity of the epidemiological situation, particularly in Riyadh. There, where at least there were two strains co-circulating, M2 stands up clearly as the best model function (see bottom panel in <xref ref-type="fig" rid="pntd.0008065.g004">Fig 4</xref>), in comparison with Macca and Madina, with only one strain (<xref ref-type="fig" rid="pntd.0008065.g004">Fig 4</xref> top panels).</p><fig id="pntd.0008065.g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g003</object-id><label>Fig 3</label><caption><title>Model simulations fitted to accumulated MERS-CoV clinical cases in Macca, Madina and Riyadh.</title><p>Observed data points are shown in blue and the solid line depicts the model solutions. The three two-strain models fitted to cumulative MERS-CoV cases are: M1: 2-strain model with bilinear incidence, M2: 2-strain model with non-monotone incidence, and M3: 2-strain model with saturated incidence.</p></caption><graphic xlink:href="pntd.0008065.g003"></graphic></fig><fig id="pntd.0008065.g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g004</object-id><label>Fig 4</label><caption><title>Same as for <xref ref-type="fig" rid="pntd.0008065.g003">Fig 3</xref> but predictions showing the best model for MERS-CoV incidence in each of the three regions.</title><p>Notice that whereas M1 provides the best prediction for both Macca and Madina (top), M2 instead does it for the two-strain situation occurring in Riyadh. Blue line is MERS-CoV cases, solid black line predictions and yellow area denotes 95% confidence interval for predictions.</p></caption><graphic xlink:href="pntd.0008065.g004"></graphic></fig></sec><sec id="sec005"><title>c) Estimation of R<sub>0</sub>, R<sub>C</sub>, and R<sub>H</sub> and prediction framework</title><p>We estimated the basic reproduction number (<italic>R</italic><sub><italic>0</italic></sub>), the community reproduction number (<italic>R</italic><sub><italic>C</italic></sub>), and the hospital reproduction number (<italic>R</italic><sub><italic>H</italic></sub>) for the whole period of data for Riyadh, Macca and Madina (<xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s018">S18</xref> Tables; 44–45). We also estimated weekly values of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub> for each prediction interval in the three provinces (i.e. weeks 105–156 in Riyadh, weeks 65–116 in Macca and weeks 63–114 in Madina, respectively). The predictions of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub> are made for each 4-week time interval in the prediction period. The estimate of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub> during 0<sup>th</sup> week is obtained using the estimated parameter values for the training period of MERS-CoV data (i.e., in Riyadh week 1–104, in Macca week 1–64 and in Madina week 1–62). Afterwards, we keep on adding 4 data points to the previous data and re-estimate the parameters. These estimated parameters were then plugged into obtain values of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub> in 4<sup>th</sup> week to 48<sup>th</sup> week of the prediction period. Thus in intervals of 4 weeks, we obtain an estimate of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub>. Temporal evolution of <italic>R</italic><sub><italic>0</italic></sub>, <italic>R</italic><sub><italic>C</italic></sub>, and <italic>R</italic><sub><italic>H</italic></sub> for the three provinces is depicted in <xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s033">S4</xref> Figs, respectively.</p><p>Using the model that provides the best prediction for the season total cases in the three provinces (i.e., a 2-strain model with BL incidence for Macca and the 2-strain model with SAT incidence for both Madina and Riyadh), we predicted the total number of cases in the following year. Using the parameter estimates in the last prediction point (i.e. 48<sup>th</sup> week) of the current season, we draw 1000 samples. Based on these parameter samples, we determined 1000 estimates of the season total cases in the next year[<xref rid="pntd.0008065.ref043" ref-type="bibr">43</xref>]. We then defined a large outbreak in the forthcoming season in a particular region, as those events when cases exceeded the total number of cases in the previous season by more than one standard deviation (see <xref ref-type="fig" rid="pntd.0008065.g005">Fig 5</xref>). Thus, the probability of a large outbreak (P<sub>l</sub>) in those provinces isdetermined as the ratio of samples that exceed the total number of cases in a season divided by the total number of samples. We also considered the case when two strains of MERS-CoV may be co-circulating in the human population in Saudi Arabia. We also simulated the situation when one strain is assumed to be more active with a higher transmission rate, whilst the other is much less transmissible among individuals in the different provinces of Saudi-Arabia[<xref rid="pntd.0008065.ref037" ref-type="bibr">37</xref>].Although we have considered two strains circulating in the community setting, we do not distinguish among strains in hospital premises. This is due to the lack of strain specific data in hospital settings. However, we are able to distinguish the contribution of infection from Community and Hospital setting by estimating the Community reproduction number (R<sub>C</sub>) and the Hospital reproduction number (R<sub>H</sub>). However, to account for this effect, we propose three different forces of infection (BL, NM and SAT) to model the MERS transmission in the three provinces. Saturated (SAT) functions describe “Crowding effects” in disease transmission, whereas Non-monotone (NM) functions are used to incorporate”psychological effects”. As the CFR in MERS-CoV is about 40%, it is easily expected that awareness in front of a MERS-CoV situation may play an important role during an epidemic[<xref rid="pntd.0008065.ref044" ref-type="bibr">44</xref>–<xref rid="pntd.0008065.ref046" ref-type="bibr">46</xref>]. Thus a model with such a crowding effect might seem a priori more realistic in comparison to the bilinear (BL) transmission (see also <xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref>). The crowding effect in disease transmission[<xref rid="pntd.0008065.ref038" ref-type="bibr">38</xref>,<xref rid="pntd.0008065.ref039" ref-type="bibr">39</xref>] can also be interpreted as the behavioral change of susceptible individuals when the number of infected individuals increases (<xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref>). This fact causes a lower number of new infections when a large number of infected individuals are already present in the community. This behavioral change in the susceptible population may occur due to, for instance, public health awareness campaigns and alerts raised by local or international authorities. The psychological effect is somewhat similar to the crowding effect, but with the effect of fear being added to the awareness. The psychological effect is stronger as the transmission rate decreases more rapidly with an increasing number of infected individuals (<xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref>) [<xref rid="pntd.0008065.ref040" ref-type="bibr">40</xref>–<xref rid="pntd.0008065.ref042" ref-type="bibr">42</xref>]. Recently, Kucharski et al.[<xref rid="pntd.0008065.ref035" ref-type="bibr">35</xref>] suggested that MERS-CoV transmission is over dispersed and hence outbreaks can include super-spreading events. In a similar study, Nishiura et al.[<xref rid="pntd.0008065.ref033" ref-type="bibr">33</xref>] concluded that super-spreaders who visited multiple healthcare facilities drove up the epidemic by generating larger number of secondary cases. Therefore, we also consider the role of potential super-spreaders in the context of a 2-strain model with BL incidence function.</p><fig id="pntd.0008065.g005" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g005</object-id><label>Fig 5</label><caption><title>Best fit obtained from a combination of 2-strain models.</title><p>Each model was aimed at predicting MERS-CoV epidemiological targets (top-peak week (No.), bottom-peak incidence (No.)) in Riyadh. Panels showing peak week and peak incidences are displayed up to the time of occurrence (11th week). Representations of the 2-strain models M1, M2, and M3 are the same as in <xref ref-type="fig" rid="pntd.0008065.g002">Fig 2A</xref>. Dashed lines are observed values.</p></caption><graphic xlink:href="pntd.0008065.g005"></graphic></fig></sec></sec><sec sec-type="conclusions" id="sec006"><title>Results and discussion</title><p>Estimates of the transmission rates for the 2-strain model with BL incidence suggest that MERS-CoV transmission in the three locations is dominated by community and hospital transmission (<xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1 Table</xref> to <xref ref-type="supplementary-material" rid="pntd.0008065.s009">S9 Table</xref>). The former statement is in good agreement with a previous study that suggested that MERS-CoV infections are essentially produced through both hospital and community based human-human transmission[<xref rid="pntd.0008065.ref002" ref-type="bibr">2</xref>]. This fact is well reflected in the estimates of transmission rates obtained from the other two 2-strain models (e.g. NM and SAT incidence functions) for the three locations (<xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s009">S9</xref> Tables). Estimates of the parameter <italic>θ</italic> (a measure of variability within the two strains) in Riyadh and Macca (see <xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s003">S3</xref> Tables) suggest that both strains are currently active in these provinces (albeit strain 2 contributes with only 3% - 50% to community human-human transmission). In Madina instead, a majority of the models (BL and NM) point to only one active strain,with strain 2 contributing less than 1% to community human-human transmission, see <xref ref-type="supplementary-material" rid="pntd.0008065.s001">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s009">S9</xref> Tables.</p><p>For each particular region in this study, two strain models are compared among themselves and also with regard to their single-strain counterparts in order to determine the best model (i.e. single, two-strain or a combination of other two-strain models). Results show that in the three locations, 2-strain models provide a better fitting to cumulative cases in comparison to their single strain counterparts (<xref ref-type="fig" rid="pntd.0008065.g002">Fig 2</xref>). AIC and BIC values for the 2-strain models suggest that the best model is region specific (best model for Riyadh is a 2-strain model with BL incidence, whereas in Macca and Madina the best models are both the model with a NM incidence function and the one incorporating super-spreaders, respectively; see <xref ref-type="supplementary-material" rid="pntd.0008065.s019">S19 Table</xref>). We also compare our 2-strain models with a 2-strain model with variable Zoonotic transmission (SI Eq. A2). Comparing AIC and BIC values, we found that the models with variable Zoonotic transmission (SI Eq. A2) did not improve all the previous results (see <xref ref-type="supplementary-material" rid="pntd.0008065.s019">S19 Table</xref>). Therefore, the best model among all 2-strain models cannot be determined solely from their goodness-of-fit comparison. For this reason, we additionally compared all the two-strain models on the basis of their respective prediction skills. Forecasts for the three models in Riyadh suggest that, a majority of the times, the 2-strain model with SAT incidence can better predict the three targets in comparison to the other two competing 2-strain models (e.g. note that for SAT the average of the mean absolute error, MAE, for peak week is 11.6, whereas for peak maximum is 24.83 and for season totals 55.62) (<xref rid="pntd.0008065.t002" ref-type="table">Table 2A</xref>). For comparison, we also provide predictions in the province of Riyadh of the three targets using the model with super-spreaders. However, prediction with super-spreaders did not improve at all the former results obtained for Riyadh (see <xref rid="pntd.0008065.t002" ref-type="table">Table 2A</xref>). For Madina, both the peak week and season totals are better predicted by the same SAT 2-strain model (i.e. MAE average for peak week being 14.5 and 9.33 for season totals; see also <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>). However, when predicting the peak maximum in Madina, clearly the best model is the 2-strain model with the NM incidence function (i.e. MAE average for peak maximum being 2.62; see also <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>). For Macca, except for the peak week, the other two targets, namely peak maximum and season totals are better predicted with the 2-strain model with BL incidence (i.e. MAE average of 7.08 for peak maximum and of 63.42 for season totals, respectively; see <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>). For peak week in Macca, again a better prediction is achieved by the 2-strain model with SAT incidence (i.e. MAE average for peak week being 23.9, see <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>).</p><table-wrap id="pntd.0008065.t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.t002</object-id><label>Table 2</label><caption><title>(A)Average predictions [Simple average of Mean Absolute Errors (MAE)] over all forecasts of the three 2-strain models and their different combinations for the Riyadh province.</title><p>(B)Estimated values of the Basic reproduction number (R<sub>0</sub>), the Hospital reproduction number (R<sub>H</sub>), and the Community reproduction number (R<sub>C</sub>)for the three provinces of Saudi Arabia using best model (two strain). The best two strain MERS model is with saturated incidence. The data are given as Mean (95% CI). *The MAE for the three models and their combinations during point prediction of peak week and peak incidence were calculated up-to the peak of the prediction season (11th Week) of the Riyadh province.</p></caption><alternatives><graphic id="pntd.0008065.t002g" xlink:href="pntd.0008065.t002"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="center" rowspan="2" style="background-color:#FFFFFF" colspan="1">A<break></break>Observed values (Data)</th><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"><underline>Peak week* (weeks)</underline></th><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"><underline>Peak maximum* (cases)</underline></th><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"><underline>Season totals (cases)</underline></th></tr><tr><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">11</th><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">45</th><th align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">230</th></tr></thead><tbody><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"><underline>Individual model forecasts</underline></td><td align="left" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td><td align="left" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td><td align="left" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M1</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">22.8 [19.5]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">18.4 [32.94]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">478.80 [69.33]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M2</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">28.7 [19.7]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">16.5 [26.25]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">216.60 [75.87]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M3</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">22.1 [11.6]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">19.2 [24.83]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">222.12 [55.62]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M1 with Superspreaders</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">44.1 [34.7]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">68.9 [32.4]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">533.71 [194.5]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"><underline>Average model forecasts</underline></td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1"></td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M1-M2</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">7.3 [12.1]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">42.67 [18.06]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">299.73 [69.92]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M1-M3</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">3.4 [9.5]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">37.63 [22.19]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">296.18 [68.31]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M2-M3</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">2.6 [9.4]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">53.37 [21.74]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">270.47 [42.79]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">M1-M2-M3</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">7.7 [9.6]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">39.37 [17.21]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">306.90 [77.41]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">B<break></break>Provinces</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">R0<break></break>[Mean (95% CI)]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">RH<break></break>[Mean (95% CI)]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">RC<break></break>[Mean (95% CI)]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">Riyadh</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">2.0706<break></break>[2.0629–2.0763]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">2.0508<break></break>[2.0442–2.0543]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">0.5657<break></break>[0.5587–0.5709]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">Macca</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">4.5716<break></break>[2.8781–6.5899]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">2.7239<break></break>[0.4281–4.9401]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">2.0977<break></break>[0.1256–5.7307]</td></tr><tr><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">Madina</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">5.0661<break></break>[2.5264–9.3814]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">4.7808<break></break>[0.7408–9.3132]</td><td align="center" style="background-color:#FFFFFF" rowspan="1" colspan="1">1.6506<break></break>[0.2116–2.9661]</td></tr></tbody></table></alternatives></table-wrap><p>When comparing the three two strain models to their single-strain counterparts, composites of the mean absolute error (MAE) in Riyadh (<xref rid="pntd.0008065.t002" ref-type="table">Table 2A</xref>, and <xref ref-type="supplementary-material" rid="pntd.0008065.s023">S23 Table</xref>) suggest that 2-strain models always outperform the single strain versions in predicting the three targets <italic>(</italic>i.e. peak week, peak maximum, and season totals). Averages of MAE for peak week, peak maximum and for season totals are 11.6, 24.83 and 55.62, respectively. This fact reinforces our earlier inference, confirming that there is more than one strain currently active in Riyadh. Among the single-strain models in Riyadh, the model with SAT incidence is found to have a better predictive capacity for the three targets in comparison to the other two single-strain models (<xref ref-type="supplementary-material" rid="pntd.0008065.s023">S23 Table</xref>). In Macca instead, the single-strain model with SAT incidence provides better predictions of the peak week, the peak maximum and seasontotals in comparison to all the other models (see for instance MAE values in <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>; and values obtained for 1-strain and 2-strain models in <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22 Table</xref>). Conversely, in Madina, both 1-strain and 2-strain models with SAT incidence provide a similar performance in predicting the three targets (e.g. MAE for peak week 16.7 in comparison to 14.5, for peak maximum of 3.33 compared to 3.6 and for season totals of 5.87 against 9.33; see <xref ref-type="supplementary-material" rid="pntd.0008065.s022">S22</xref> and <xref ref-type="supplementary-material" rid="pntd.0008065.s023">S23</xref> Tables). The latter may likely represent the best choice of models to form the basis of a future early-warning system for MERS-CoV prediction in the region.</p><p>The best model, selected on the basis of its capacity to derive skillful predictions, is the model with the SAT incidence function. We used this model to estimate the basic reproduction number, R<sub>0</sub>, the community reproduction number, R<sub>C</sub>, and the hospital reproduction number, R<sub>H</sub>, for the whole period of data available in those three provinces (<xref rid="pntd.0008065.t002" ref-type="table">Table 2B</xref> and <xref ref-type="supplementary-material" rid="pntd.0008065.s020">S20 Table</xref>). In all three provinces, R<sub>0</sub> is estimated to be always greater than unity. In all three provinces, R<sub>H</sub> is found out to be larger than R<sub>C</sub>. This implies, MERS-CoV transmission triggered from the hospital setting (see <xref rid="pntd.0008065.t002" ref-type="table">Table 2B</xref>). These estimates agree with some previous values reported in the literature[<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>,<xref rid="pntd.0008065.ref017" ref-type="bibr">17</xref>–<xref rid="pntd.0008065.ref018" ref-type="bibr">18</xref>], while being a result well supported by other R<sub>0</sub> estimates[<xref rid="pntd.0008065.ref029" ref-type="bibr">29</xref>–<xref rid="pntd.0008065.ref031" ref-type="bibr">31</xref>, <xref rid="pntd.0008065.ref033" ref-type="bibr">33</xref>]. However, previous modeling attempts to model Saudi Arabia MERS-CoV clinical incidence used only a 1-strain model with BL transmission[<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>]. In our case, versions of R<sub>0</sub> were also estimated from our BL 1-strain model in those same three provinces (see <xref ref-type="supplementary-material" rid="pntd.0008065.s021">S21 Table</xref>) and they are in good agreement with the previous values provided by Chowell et al.[<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>]. At the same time, though, 1-strain model show less predictive capacity than their similar two-strain counterparts. To further verify the robustness of these estimates, the temporal evolution of R<sub>0</sub>, R<sub>C</sub>, and R<sub>H</sub> are displayed for different time intervals in the three provinces (<xref ref-type="supplementary-material" rid="pntd.0008065.s030">S1</xref>–<xref ref-type="supplementary-material" rid="pntd.0008065.s033">S4</xref> Figs). Considering the best model configurations, we additionally computed the temporal evolution of R<sub>0</sub>, R<sub>C</sub>, and R<sub>H</sub> in the three provinces. Temporal changes in R<sub>0</sub> (<xref ref-type="supplementary-material" rid="pntd.0008065.s030">S1 Fig</xref>) indicate that in most of the predicted weeks, R<sub>0</sub> stays well above the epidemic threshold (R<sub>0</sub> = 1). This fact is well established from the temporal evolution of R<sub>C</sub>, and R<sub>H</sub> in those three provinces (<xref ref-type="supplementary-material" rid="pntd.0008065.s030">S1</xref> and <xref ref-type="supplementary-material" rid="pntd.0008065.s031">S2</xref> Figs).</p><p>As the basis for an operational EWS for the region, we predicted the total number of cases in the out-of-fit interval covering July 2016 to July 2017. We used the best individual model for the three provinces among all 2-strain and single strain models developed. Figs <xref ref-type="fig" rid="pntd.0008065.g006">6</xref> and <xref ref-type="fig" rid="pntd.0008065.g007">7</xref> and <xref rid="pntd.0008065.t003" ref-type="table">Table 3</xref> clearly indicate the 2016–2017 season might be ripe for a larger outbreak in Macca and Riyadh. Instead, in Madina the likelihood of suffering a larger outbreak is very low (<xref ref-type="fig" rid="pntd.0008065.g006">Fig 6</xref> and <xref rid="pntd.0008065.t003" ref-type="table">Table 3</xref>). <xref ref-type="fig" rid="pntd.0008065.g007">Fig 7</xref> displays the single-strain model fitting to new MERS-CoV cases. Here M1 refers to the single strain model with BL force of infection, M2 to the single strain model with non-monotone incidence, and M3 to the single strain model with saturated incidence. Solid black curve represents model solution and yellow region denotes 95% confidence interval for predictions. Simulations for both Macca and Riyadh reflect quite appropriately the dynamics displayed in observations (<xref ref-type="fig" rid="pntd.0008065.g007">Fig 7</xref>). <xref ref-type="fig" rid="pntd.0008065.g008">Fig 8</xref> shows All-season’s predictions for the three consecutive years from July 6, 2013 to June 28, 2016, with each interval of 52 weeks fitted shown in panels S1-S3.</p><fig id="pntd.0008065.g006" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g006</object-id><label>Fig 6</label><caption><title>Out-of-fit prediction of large outbreaks in Macca, Madina and Riyadh in July 2016 to July 2017.</title><p>Large outbreak size (red circles) are defined as those samples which exceeds previous year (July 2015 to July 2016) total season cases. Blue circles denote those samples that fall below previous year total cases. The black line denotes total cases during July 2015 to July 2016 of the data. M1: two strain model with bilinear incidence, M3: two strain model with saturated incidence.</p></caption><graphic xlink:href="pntd.0008065.g006"></graphic></fig><fig id="pntd.0008065.g007" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g007</object-id><label>Fig 7</label><caption><title>Single-strain model fitting to new MERS-CoV cases.</title><p>Here M1: single strain model with BL force of infection, M2: single strain model with non-monotone incidence, and M3: single strain model with saturated incidence. Solid black curve represents model solution and yellow region denotes 95% confidence interval for predictions.</p></caption><graphic xlink:href="pntd.0008065.g007"></graphic></fig><fig id="pntd.0008065.g008" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.g008</object-id><label>Fig 8</label><caption><title>Season-wise model fitting to cumulative cases of Riyadh from July 6, 2013 to June 28, 2016.</title><p>(Model fitted is a 2-strain model with saturated incidence. One year appears in each panel S1 to S3(S1 = S2 = S3 = 52 weeks). Line and circle line refers to observations and 95%confidence interval for cumulative predictions denoted by the yellow region.</p></caption><graphic xlink:href="pntd.0008065.g008"></graphic></fig><table-wrap id="pntd.0008065.t003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pntd.0008065.t003</object-id><label>Table 3</label><caption><title>Hindcast prediction of season total MERS-CoV cases for July 2016 to July 2017 and probability of large outbreak.</title><p>Prediction based upon the best two-strain model (in terms of season totals forecast) for the three provinces of the Saudi-Arabia.</p></caption><alternatives><graphic id="pntd.0008065.t003g" xlink:href="pntd.0008065.t003"></graphic><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col><col align="left" valign="middle" span="1"></col></colgroup><thead><tr><th align="left" rowspan="1" colspan="1">Two strain models</th><th align="left" rowspan="1" colspan="1">MERS-CoV cases for July 2015 –July 2016</th><th align="left" rowspan="1" colspan="1">Predicted cases of MERS-CoV for the period July 2016 –July 2017 (Mean (95% CI))</th><th align="left" rowspan="1" colspan="1">Probability of larger outbreak ((No. of samples above 2015/2016 cases)/(total samples))</th></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Riyadh (M3)</td><td align="left" rowspan="1" colspan="1">230</td><td align="left" rowspan="1" colspan="1">248.1724 (213.7383–268.26)</td><td align="left" rowspan="1" colspan="1">0.7150</td></tr><tr><td align="left" rowspan="1" colspan="1">Macca (M1)</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">69.0813(59.7604–78.9765)</td><td align="left" rowspan="1" colspan="1">1</td></tr><tr><td align="left" rowspan="1" colspan="1">Madina (M3)</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">7.5268 (7.5087–7.5598)</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives></table-wrap><p>According to WHO reports, 249 MERS-CoV cases including 75 deaths (CFR 30%) were reported from Saudi Arabia between July 2016 and July 2017[<xref rid="pntd.0008065.ref044" ref-type="bibr">44</xref>]. In the aforementioned period, at least 108 new cases[<xref rid="pntd.0008065.ref044" ref-type="bibr">44</xref>] were reported from Riyadh province and at least 5 cases were reported from Madina. These values indicate that a larger epidemic did not occur in the last season (2016–2017), the risk for a larger one in the coming seasons still remains high.</p><sec id="sec007"><title>Conclusions</title><p>Up to July 2018, 2237 new cases were reported, with 1861 only in KSA and 793 deaths[<xref rid="pntd.0008065.ref045" ref-type="bibr">45</xref>](CFR 35.5%). This is an alarming situation as previous predictions on MERS-CoV had instead suggested that MERS-CoV might not sustain as an epidemic in the Arabian Peninsula. The WHO report[<xref rid="pntd.0008065.ref044" ref-type="bibr">44</xref>] suggests that MERS-CoV is still a relatively rare disease about which the medical personnel in health-care facilities have low awareness. Globally, MERS-CoV awareness is limited and because symptoms of MERS-CoV infection are non-specific, initial cases can be sometimes easily missed. With improved compliance in infection prevention and control, namely by stricter adherence to the standard precautions at all times, human-to-human transmission in health-care facilities can be reduced and even possibly eliminated with additional use of transmission-based precautions. In that regard, predictive mathematical models can help strengthen our understanding of both MERS-CoV transmission and control.</p><p>In this study, we addressed the capacity of predictive mathematical models based on two-strain MERS-CoV configurations having different transmission functions. The models differed from each other in their force of infection and in how they cope with heterogeneity in transmission. Estimates of transmission rates suggest that community and hospital transmission are dominant in the case of 2-strain models in Riyadh, Macca and Madina. The majority of the 2-strain models suggest that MERS-CoV transmission is dominated by community and hospital human-human transmissions, a fact that reflects the actual transmission scenario in Saudi Arabia [<xref rid="pntd.0008065.ref002" ref-type="bibr">2</xref>–<xref rid="pntd.0008065.ref004" ref-type="bibr">4</xref>]. Estimates of the parameter that measures transmission diversity between the two strains in the three provinces suggest that two strains are only active in Riyadh. This opposite trend in Riyadh in comparison to the other two provinces may be due to the fact that Riyadh is the capital city with good large health care facilities and a majority of the MERS-CoV patients in Saudi Arabia come to Riyadh for treatment[<xref rid="pntd.0008065.ref002" ref-type="bibr">2</xref>]. These patients may therefore carry different MERS strains, ultimately leading to multiple strains being presently co-circulating in Riyadh. Similarly, Cotten et al.[<xref rid="pntd.0008065.ref037" ref-type="bibr">37</xref>] found that ancestors of most of the viral clades originated from Riyadh.</p><p>We compared among the 2-strain models according to their predictive performance with regard to three targets (i.e. peak week, peak maximum and season totals). Our results suggest that among the three 2-strain models, the model with SAT incidence provides consistently skillful predictions and may be used to date as the best predictive model for MERS-CoV in Riyadh. Riyadh iswhere most of the MERS-CoV cases occur, while for the provinces of Macca and Madina, with lower reported MERS cases, it is difficult to determine the best model among the three 2-strain models. This fact justifies our earlier finding that in Riyadh two different strains are currently active and therefore the performance of the 2-strain models is better there. As per our results for Macca and Madina, only one dominant strain is active in those provinces. Therefore, predictions based on single strain models are there more appropriate. Our results also suggest that among the single strain models, those with SAT incidence always accurately predict the three targets for these two provinces. Thus, a dynamical MERS model considering this crowding effect is the most appropriate configuration to cope with the nature of MERS-CoV transmission.</p><p>We estimated R<sub>0</sub> using the best 2-strain model in Riyadh and estimates are in good agreement with the findings of Majumdar et al. [<xref rid="pntd.0008065.ref028" ref-type="bibr">28</xref>].The finding that R<sub>0</sub> is most of the time above 1 (<xref ref-type="supplementary-material" rid="pntd.0008065.s030">S1 Fig</xref>) is well supported by some previous estimates in literature[<xref rid="pntd.0008065.ref028" ref-type="bibr">28</xref>–<xref rid="pntd.0008065.ref030" ref-type="bibr">30</xref>,<xref rid="pntd.0008065.ref032" ref-type="bibr">32</xref>]. Lower contribution of community transmission in R<sub>0</sub> (See <xref rid="pntd.0008065.t002" ref-type="table">Table 2B</xref>) in Riyadh and Macca suggests that MERS-CoV transmission is triggered from hospital settings in those provinces. Most interestingly, in some forecasted periods, R<sub>0</sub> attains large values, a fact that denotes that a rapid propagation among the susceptible population is indeed possible. More worrisome is the range of values into which R<sub>0</sub> moves, most of the time above 1 and below 2,5, a fact that makes it a dangerous infection in terms of silent and constant potential population spread. Community reproduction number R<sub>C</sub> well above unity (see <xref ref-type="supplementary-material" rid="pntd.0008065.s033">S4 Fig</xref>) in most of the predictive weeks indicate that in a near future a large outbreak may be possible in those provinces. Out-of-fit predictions for the next season totals suggest that there is a high possibility of larger outbreaks in Macca and Riyadh. However, our results instead indicate that there is a very low possibility of larger outbreaks in Madina. The fact that this outbreak did not happen in 2017–2018 does not preclude what may occur in the forthcoming seasons. Under such a scenario, authorities and international health agencies should prepare and actively work towards the prompt implementation of cheap albeit efficient computational platforms ready to assist in the simulation of how a potential outbreak might evolve in the region. More so, given the high probability that another large MERS-CoV outbreak occurs in the Arabian Peninsula or nearby countries. Migration may also play a major role for increased transmission in the provinces of Saudi Arabia and this feature should be properly accounted for in future model configurations. In summary, our findings suggest that in a majority of provinces a single MERS-CoV strain is currently active, conversely to the situation in Riyadh. However, in the near future, it is also possible that more general MERS-CoV transmission occurs from multiple strains in other provinces of Saudi Arabia.</p></sec></sec><sec sec-type="supplementary-material" id="sec008"><title>Supporting information</title><supplementary-material content-type="local-data" id="pntd.0008065.s001"><label>S1 Table</label><caption><title>Estimated parameters for Model-(A) with bilinear incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s001.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s002"><label>S2 Table</label><caption><title>Estimated parameters for Model-(A) with non-monotone incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s002.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s003"><label>S3 Table</label><caption><title>Estimated parameters for Model-(A) with saturated incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s003.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s004"><label>S4 Table</label><caption><title>Estimated parameters for Model-(A) with bilinear incidence for the Mecca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s004.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s005"><label>S5 Table</label><caption><title>Estimated parameters for Model-(A) with non-monotone incidence for the Mecca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s005.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s006"><label>S6 Table</label><caption><title>Estimated parameters for Model-(A) with Saturated incidence for the Mecca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s006.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s007"><label>S7 Table</label><caption><title>Estimated parameters for Model-(A) with bilinear incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s007.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s008"><label>S8 Table</label><caption><title>Estimated parameters for Model-(A) with non-monotone incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s008.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s009"><label>S9 Table</label><caption><title>Estimated parameters for Model-(A) with saturated incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s009.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s010"><label>S10 Table</label><caption><title>Estimated parameters for the Model (B) with bilinear incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s010.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s011"><label>S11 Table</label><caption><title>Estimated parameters for the Model (B) with non-monotone incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s011.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s012"><label>S12 Table</label><caption><title>Estimated parameters for the Model (B) with saturated incidence for the Riyadh province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s012.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s013"><label>S13 Table</label><caption><title>Estimated parameters for the Model (B) with bilinear incidence for the Macca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s013.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s014"><label>S14 Table</label><caption><title>Estimated parameters for the Model (B) with non-monotone incidence for the Macca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s014.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s015"><label>S15 Table</label><caption><title>Estimated parameters for the Model (B) with saturated incidence for the Macca province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s015.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s016"><label>S16 Table</label><caption><title>Estimated parameters for the Model (B) with bilinear incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s016.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s017"><label>S17 Table</label><caption><title>Estimated parameters for the Model (B) with non-monotone incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s017.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s018"><label>S18 Table</label><caption><title>Estimated parameters for the Model (B) with saturated incidence for the Madina province.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s018.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s019"><label>S19 Table</label><caption><title>Multi-model inference quantities (AIC and BIC) for three two-strain and single-strain models.</title><p>Model with smallest AIC and BIC value are given in bold.</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s019.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s020"><label>S20 Table</label><caption><title>Estimated values of the Basic reproduction number (R0), the Hospital reproduction number (RH), and the Community reproduction number (RC), for the three provinces of Saudi Arabia for the two strain Model-1 and Model-2 (Equation (A) with bilinear incidence and non-monotone incidence).</title><p>The data is given in the format (Mean [95%CI]).</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s020.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s021"><label>S21 Table</label><caption><title>Comparison of estimated values of the Basic reproduction number (R0), the Hospital reproduction number (RH), and the Community reproduction number (RC), for the three provinces of Saudi Arabia for the single strain Model-1, Model-2 and Model-3 (Equation (B) with bilinear, non-monotone and saturated incidence).</title><p>The data is given in the format (Mean [95%CI]).</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s021.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s022"><label>S22 Table</label><caption><title>Average predictions [Simple average of Mean Absolute Errors (MAE)] obtained over all the prediction weeks using Model -(A) with different incidence functions.</title><p>Model -1 represents Model-(A) with bilinear incidence function. Model -2 represents Model-(A) with non-monotone incidence and Model -3 represents Model -(A) with saturated incidence.</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s022.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s023"><label>S23 Table</label><caption><title>Aver.age predictions [Simple average of Mean Absolute Errors (MAE)] obtained over all the prediction weeks using a single-strain Model—(B) with different incidence functions.</title><p>Model -1 represents Model-(B) with bilinear incidence function. Model -2 represents Model-(B) with non-monotone incidence and Model-3 represents Model-(B) with saturated incidence.</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s023.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s024"><label>S24 Table</label><caption><title>Estimated parameters for Model-(A1) for Riyadh.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s024.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s025"><label>S25 Table</label><caption><title>Estimated parameters for Model-(A1) for Macca.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s025.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s026"><label>S26 Table</label><caption><title>Estimated parameters for Model-(A1) for Madina.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s026.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s027"><label>S27 Table</label><caption><title>Estimated parameters for the Model (B1) for Riyadh.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s027.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s028"><label>S28 Table</label><caption><title>Estimated parameters for the Model (B1) for Macca.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s028.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s029"><label>S29 Table</label><caption><title>Estimated parameters for the Model (B1) for Madina.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s029.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s030"><label>S1 Fig</label><caption><title>Temporal Evolution of R0 using the best predicted 2-strain model (Saturated incidence) in three provinces; Riyadh, Macca and Madina.</title><p>R0 is estimated over different forecast weeks (0, 4, 8, …, 48).</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s030.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s031"><label>S2 Fig</label><caption><title>Temporal Evolution of Hospital reproduction number (RH) using the best predicted 2-strain model (Saturated incidence) in three provinces; Riyadh, Macca and Madina. RH is estimated over different forecast weeks (0, 4, 8, …, 48).</title><p>Dotted line represents the threshold of the epidemic potential (RH = 1).</p><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s031.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s032"><label>S3 Fig</label><caption><title>Temporal Evolution of Community reproduction number (RC) using the best predicted 2-strain model (Saturated incidence) in three provinces; Riyadh, Macca and Madina. RC is estimated over different forecast weeks (0, 4, 8, …, 48).</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s032.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material><supplementary-material content-type="local-data" id="pntd.0008065.s033"><label>S4 Fig</label><caption><title>Model simulations fitted to accumulated MERS-CoV clinical cases in Macca, Madina and Riyadh.</title><p>(DOCX)</p></caption><media xlink:href="pntd.0008065.s033.docx"><caption><p>Click here for additional data file.</p></caption></media></supplementary-material></sec></body><back><ack><p>The authors thank Joseph Boyard-Micheau for assistance with <xref ref-type="fig" rid="pntd.0008065.g001">Fig 1</xref> drawing.</p></ack><ref-list><title>References</title><ref id="pntd.0008065.ref001"><label>1</label><mixed-citation publication-type="journal"><name><surname>Sabir</surname><given-names>J.S.M.</given-names></name><etal>et al</etal><article-title>Co-circulation of three camel coronavirus species and recombination of MERS-CoVs in Saudi Arabia</article-title>. <source>Science</source><volume>351</volume>, <fpage>81</fpage>–<lpage>84</lpage> (<year>2016</year>). <pub-id pub-id-type="doi">10.1126/science.aac8608</pub-id><pub-id pub-id-type="pmid">26678874</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref002"><label>2</label><mixed-citation publication-type="other">Ministry of Health (Saudi Arabia), Middle East Respiratory Syndrome Coronavirus (MERS-CoV); <ext-link ext-link-type="uri" xlink:href="http://www.moh.gov.sa/en/CCC/PressReleases/Pages/default.aspx?PageIndex=1(2017)">http://www.moh.gov.sa/en/CCC/PressReleases/Pages/default.aspx?PageIndex=1(2017)</ext-link>.</mixed-citation></ref><ref id="pntd.0008065.ref003"><label>3</label><mixed-citation publication-type="other">World Health Organization (WHO), East Respiratory Syndrome Coronavirus (MERS-CoV)—Saudi Arabia; <ext-link ext-link-type="uri" xlink:href="http://www.who.int/csr/don/06-july-2017-mers-saudi-arabia/en/">http://www.who.int/csr/don/06-july-2017-mers-saudi-arabia/en/</ext-link> (2017).</mixed-citation></ref><ref id="pntd.0008065.ref004"><label>4</label><mixed-citation publication-type="journal"><name><surname>Chowell</surname><given-names>G.</given-names></name>, <name><surname>Blumberg</surname><given-names>S.</given-names></name>, <name><surname>Simonsen</surname><given-names>L.</given-names></name>, <name><surname>Miller</surname><given-names>M.A.</given-names></name> &<name><surname>Viboud</surname><given-names>C</given-names></name>. <article-title>Synthesizing data and models for the spread of MERS-CoV, 2013: key role of index cases and hospital transmission</article-title>. <source>Epidemics</source><volume>9</volume>, <fpage>40</fpage>–<lpage>51</lpage> (<year>2014</year>). <pub-id pub-id-type="doi">10.1016/j.epidem.2014.09.011</pub-id><pub-id pub-id-type="pmid">25480133</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref005"><label>5</label><mixed-citation publication-type="journal"><name><surname>Zaki</surname><given-names>A.M.</given-names></name>, <name><surname>Van Boheemen</surname><given-names>S.</given-names></name>, <name><surname>Bestebroer</surname><given-names>T.M.</given-names></name>, <name><surname>Osterhaus</surname><given-names>A.D.</given-names></name> &<name><surname>Fouchier</surname><given-names>R.A</given-names></name>. <article-title>Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia</article-title>. <source>N. Engl. J. Med</source>. <volume>367</volume>, <fpage>1814</fpage>–<lpage>1820</lpage> (<year>2012</year>). <pub-id pub-id-type="doi">10.1056/NEJMoa1211721</pub-id><pub-id pub-id-type="pmid">23075143</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref006"><label>6</label><mixed-citation publication-type="journal"><name><surname>Lee</surname><given-names>S.S.</given-names></name> & <name><surname>Wong</surname><given-names>N.S.</given-names></name><article-title>Probable transmission chains of Middle East respiratory syndrome coronavirus and the multiple generations of secondary infection in South Korea</article-title>. <source>Int. J. Infect. Dis</source>. <volume>38</volume>, <fpage>65</fpage>–<lpage>67</lpage> (<year>2015</year>). <pub-id pub-id-type="doi">10.1016/j.ijid.2015.07.014</pub-id><pub-id pub-id-type="pmid">26216766</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref007"><label>7</label><mixed-citation publication-type="journal"><name><surname>Omrani</surname><given-names>A.S.</given-names></name>, <name><surname>Al-Tawfiq</surname><given-names>J.A.</given-names></name> &<name><surname>Memish</surname><given-names>Z.A</given-names></name>. <article-title>Middle East respiratory syndrome coronavirus (MERS-CoV): animal to human interaction</article-title>. <source>Pathog. Glob. Health</source><volume>109</volume>(<issue>8</issue>), <fpage>354</fpage>–<lpage>362</lpage> (<year>2015</year>). <pub-id pub-id-type="doi">10.1080/20477724.2015.1122852</pub-id><pub-id pub-id-type="pmid">26924345</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref008"><label>8</label><mixed-citation publication-type="journal"><name><surname>Hemida</surname><given-names>M.G.</given-names></name><etal>et al</etal><article-title>Lack of Middle East respiratory syndrome coronavirus transmission from infected camels</article-title>. <source>Emerg. Infect. Diseases</source><volume>21</volume>(<issue>4</issue>), <fpage>699</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25811546</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref009"><label>9</label><mixed-citation publication-type="journal"><name><surname>Cotten</surname><given-names>M.</given-names></name><etal>et al</etal><article-title>Transmission and evolution of the Middle East respiratory syndrome coronavirus in Saudi Arabia: a descriptive genomic study</article-title>. <source>The Lancet</source><volume>382</volume>(<issue>9909</issue>),<fpage>1993</fpage>–<lpage>2002</lpage> (<year>2013</year>).</mixed-citation></ref><ref id="pntd.0008065.ref010"><label>10</label><mixed-citation publication-type="journal"><name><surname>Lingshu</surname><given-names>W.</given-names></name><etal>et al</etal><article-title>Evaluation of candidate vaccine approaches for MERS-CoV</article-title>. <source>Nat. Commun</source>.<volume>6</volume>(<issue>7712</issue>), (<year>2015</year>).</mixed-citation></ref><ref id="pntd.0008065.ref011"><label>11</label><mixed-citation publication-type="journal"><name><surname>Myers</surname><given-names>M.F.</given-names></name>, <name><surname>Rogers</surname><given-names>D.J.</given-names></name>, <name><surname>Cox</surname><given-names>J.</given-names></name>, <name><surname>Flahault</surname><given-names>A.</given-names></name> & <name><surname>Hay</surname><given-names>S.I.</given-names></name><article-title>Forecasting disease risk for increased epidemic preparedness in public health</article-title>. <source>Adv.Parasitol</source>. <volume>47</volume>, <fpage>309</fpage>–<lpage>330</lpage> (<year>2000</year>). <pub-id pub-id-type="doi">10.1016/s0065-308x(00)47013-2</pub-id><pub-id pub-id-type="pmid">10997211</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref012"><label>12</label><mixed-citation publication-type="journal"><name><surname>Shaman</surname><given-names>J.</given-names></name>, <name><surname>Karspeck</surname><given-names>A.</given-names></name>, <name><surname>Yang</surname><given-names>W.</given-names></name>, <name><surname>Tamerius</surname><given-names>J.</given-names></name> &<name><surname>Lipsitch</surname><given-names>M</given-names></name>. <article-title>Real-time influenza forecasts during the 2012–2013 season</article-title>. <source>Nat.Commun</source>. <volume>4</volume>, <fpage>2837</fpage> (<year>2013</year>). <pub-id pub-id-type="doi">10.1038/ncomms3837</pub-id><pub-id pub-id-type="pmid">24302074</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref013"><label>13</label><mixed-citation publication-type="journal"><name><surname>Biggerstaff</surname><given-names>M.</given-names></name><etal>et al</etal><article-title>Results from the centers for disease control and prevention’s predict the 2013–2014 Influenza Season Challenge</article-title>. <source>BMC Infect. Dis</source>. <volume>16</volume>(<issue>1</issue>), <fpage>357</fpage> (<year>2016</year>).<pub-id pub-id-type="pmid">27449080</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref014"><label>14</label><mixed-citation publication-type="journal"><name><surname>Shaman</surname><given-names>J.</given-names></name><article-title>&Karspeck, A. Forecasting seasonal outbreaks of influenza</article-title>. <source>Proc. Natl. Acad. Sci. USA</source><volume>109</volume>(<issue>50</issue>), <fpage>20425</fpage>–<lpage>20430</lpage> (<year>2012</year>). <pub-id pub-id-type="doi">10.1073/pnas.1208772109</pub-id><pub-id pub-id-type="pmid">23184969</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref015"><label>15</label><mixed-citation publication-type="journal"><name><surname>Yamana</surname><given-names>T.K.</given-names></name>, <name><surname>Kandula</surname><given-names>S.</given-names></name> & <name><surname>Shaman</surname><given-names>J.</given-names></name><article-title>Superensemble forecasts of dengue outbreaks</article-title>. <source>J. R. Soc. Interface</source><volume>13</volume>(<issue>123</issue>), <fpage>20160410</fpage> (<year>2016</year>). <pub-id pub-id-type="doi">10.1098/rsif.2016.0410</pub-id><pub-id pub-id-type="pmid">27733698</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref016"><label>16</label><mixed-citation publication-type="journal"><name><surname>Ray</surname><given-names>E.L.</given-names></name> & <name><surname>Reich</surname><given-names>N.G.</given-names></name><article-title>Prediction of infectious disease epidemics via weighted density ensembles</article-title>.<source>PLoSComput. Biol</source>. <volume>14</volume>(<issue>2</issue>), <fpage>e1005910</fpage> (<year>2018</year>).</mixed-citation></ref><ref id="pntd.0008065.ref017"><label>17</label><mixed-citation publication-type="journal"><name><surname>Breban</surname><given-names>R.</given-names></name>, <name><surname>Riou</surname><given-names>J.</given-names></name> &<name><surname>Fontanet</surname><given-names>A</given-names></name>. <article-title>Interhuman transmissibility of Middle East respiratory syndrome coronavirus: estimation of pandemic risk</article-title>. <source>The Lancet</source><volume>382</volume>(<issue>9893</issue>), <fpage>694</fpage>–<lpage>699</lpage> (<year>2013</year>).</mixed-citation></ref><ref id="pntd.0008065.ref018"><label>18</label><mixed-citation publication-type="journal"><name><surname>Cauchemez</surname><given-names>S.</given-names></name><etal>et al</etal><article-title>Unraveling the drivers of MERS-CoV transmission</article-title>. <source>Proc. Natl. Acad. Sci. USA</source><volume>113</volume>(<issue>32</issue>), <fpage>9081</fpage>–<lpage>9086</lpage> (<year>2016</year>). <pub-id pub-id-type="doi">10.1073/pnas.1519235113</pub-id><pub-id pub-id-type="pmid">27457935</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref019"><label>19</label><mixed-citation publication-type="journal"><name><surname>Zumla</surname><given-names>A.</given-names></name>, <name><surname>Hui</surname><given-names>D.S.</given-names></name> & <name><surname>Perlman</surname><given-names>S.</given-names></name><article-title>Middle East respiratory syndrome</article-title>. <source>The Lancet</source><volume>386</volume>(<issue>9997</issue>), <fpage>995</fpage>–<lpage>1007</lpage> (<year>2015</year>).</mixed-citation></ref><ref id="pntd.0008065.ref020"><label>20</label><mixed-citation publication-type="journal"><name><surname>Muth</surname><given-names>D.</given-names></name><etal>et al</etal><article-title>Infectious Middle East respiratory syndrome coronavirus excretion and serotype variability based on live virus isolates from patients in Saudi Arabia</article-title>. <source>J. Clin. Microbiol</source>. <volume>53</volume>(<issue>9</issue>), <fpage>2951</fpage>–<lpage>2955</lpage> (<year>2015</year>). <pub-id pub-id-type="doi">10.1128/JCM.01368-15</pub-id><pub-id pub-id-type="pmid">26157150</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref021"><label>21</label><mixed-citation publication-type="journal"><name><surname>Capasso</surname><given-names>V.</given-names></name><article-title>&Serio, G. A generalization of the Kermack-McKendrick deterministic epidemic model</article-title>. <source>Math.Biosci</source>. <volume>42</volume>(<issue>1–2</issue>), <fpage>43</fpage>–<lpage>61</lpage> (<year>1978</year>).</mixed-citation></ref><ref id="pntd.0008065.ref022"><label>22</label><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>W.</given-names></name><article-title>Epidemic models with nonlinear infection forces. Math.Biosci</article-title>. <source>Eng</source>. <volume>3</volume>(<issue>1</issue>), <fpage>267</fpage>–<lpage>279</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="pntd.0008065.ref023"><label>23</label><mixed-citation publication-type="journal"><name><surname>Xiao</surname><given-names>D.</given-names></name><article-title>&Ruan, S. Global analysis of an epidemic model with nonmonotone incidence rate</article-title>. <source>Math.Biosci</source>. <volume>208</volume>(<issue>2</issue>), <fpage>419</fpage>–<lpage>29</lpage> (<year>2007</year>). <pub-id pub-id-type="doi">10.1016/j.mbs.2006.09.025</pub-id><pub-id pub-id-type="pmid">17303186</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref024"><label>24</label><mixed-citation publication-type="other">Life expectancy, Saudi Arabia; <ext-link ext-link-type="uri" xlink:href="http://countryeconomy.com/demography/">http://countryeconomy.com/demography/</ext-link> (2016).</mixed-citation></ref><ref id="pntd.0008065.ref025"><label>25</label><mixed-citation publication-type="journal"><name><surname>Lessler</surname><given-names>J.</given-names></name><etal>et al</etal><article-title>Incubation periods of acute respiratory viral infections: a systematic review</article-title>. <source>Lancet Infect. Dis</source>. <volume>9</volume>(<issue>5</issue>), <fpage>291</fpage>–<lpage>300</lpage> (<year>2009</year>). <pub-id pub-id-type="doi">10.1016/S1473-3099(09)70069-6</pub-id><pub-id pub-id-type="pmid">19393959</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref026"><label>26</label><mixed-citation publication-type="journal"><name><surname>Assiri</surname><given-names>A.</given-names></name><etal>et al</etal><article-title>Hospital outbreak of Middle East respiratory syndrome coronavirus</article-title>. <source>N. Engl. J. Med</source>. <volume>2013</volume>(<issue>369</issue>), <fpage>407</fpage>–<lpage>16</lpage> (<year>2013</year>).</mixed-citation></ref><ref id="pntd.0008065.ref027"><label>27</label><mixed-citation publication-type="journal"><name><surname>Xia</surname><given-names>Z.Q.</given-names></name>, <name><surname>Zhang</surname><given-names>J.</given-names></name>, <name><surname>Xue</surname><given-names>Y.K.</given-names></name>, <name><surname>Sun</surname><given-names>G.Q.</given-names></name> &<name><surname>Jin</surname><given-names>Z</given-names></name>. <article-title>Modeling the transmission of Middle East respirator syndrome corona virus in the Republic of Korea</article-title>. <source>PloSOne</source><volume>10</volume>(<issue>12</issue>), <fpage>e0144778</fpage> (<year>2015</year>).</mixed-citation></ref><ref id="pntd.0008065.ref028"><label>28</label><mixed-citation publication-type="journal"><name><surname>Majumder</surname><given-names>M.S.</given-names></name>, <name><surname>Rivers</surname><given-names>C.</given-names></name>, <name><surname>Lofgren</surname><given-names>E.</given-names></name> &<name><surname>Fisman</surname><given-names>D</given-names></name>. <article-title>Estimation of MERS-coronavirus reproductive number and case fatality rate for the spring 2014 Saudi Arabia outbreak: insights from publicly available data</article-title>. <source>PLoSCurrents</source><volume>6</volume>, (<year>2014</year>).</mixed-citation></ref><ref id="pntd.0008065.ref029"><label>29</label><mixed-citation publication-type="journal"><name><surname>Poletto</surname><given-names>C.</given-names></name>, <name><surname>Colizza</surname><given-names>V.</given-names></name> &<name><surname>Boëlle</surname><given-names>P.Y</given-names></name>. <article-title>Quantifying spatiotemporal heterogeneity of MERS-CoV transmission in the Middle East region: a combined modelling approach</article-title>. <source>Epidemics</source><volume>15</volume>, <fpage>1</fpage>–<lpage>9</lpage> (<year>2016</year>). <pub-id pub-id-type="doi">10.1016/j.epidem.2015.12.001</pub-id><pub-id pub-id-type="pmid">27266844</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref030"><label>30</label><mixed-citation publication-type="journal"><name><surname>Nishiura</surname><given-names>H.</given-names></name>, <name><surname>Miyamatsu</surname><given-names>Y.</given-names></name>, <name><surname>Chowell</surname><given-names>G.</given-names></name> &<name><surname>Saitoh</surname><given-names>M</given-names></name>. <article-title>Assessing the risk of observing multiple generations of Middle East respiratory syndrome (MERS) cases given an imported case</article-title>. <source>Eurosurveillance</source><volume>20</volume>(<issue>27</issue>), <fpage>6</fpage>–<lpage>11</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26290487</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref031"><label>31</label><mixed-citation publication-type="journal"><name><surname>Poletto</surname><given-names>C.</given-names></name>, <name><surname>Boëlle</surname><given-names>P.Y.</given-names></name> &<name><surname>Colizza</surname><given-names>V</given-names></name>. <article-title>Risk of MERS importation and onward transmission: a systematic review and analysis of cases reported to WHO</article-title>. <source>BMC Infect.Dis</source>. <volume>16</volume>(<issue>1</issue>), <fpage>448</fpage> (<year>2016</year>). <pub-id pub-id-type="doi">10.1186/s12879-016-1787-5</pub-id><pub-id pub-id-type="pmid">27562369</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref032"><label>32</label><mixed-citation publication-type="journal"><name><surname>Hsieh</surname><given-names>Y.H.</given-names></name><article-title>Middle East respiratory syndrome coronavirus (MERS-CoV) nosocomial outbreak in South Korea: insights from modeling</article-title>. <source>PeerJ</source><volume>3</volume>, <fpage>1505</fpage> (<year>2015</year>).</mixed-citation></ref><ref id="pntd.0008065.ref033"><label>33</label><mixed-citation publication-type="journal"><name><surname>Nishiura</surname><given-names>H.</given-names></name><etal>et al</etal><article-title>Identifying determinants of heterogeneous transmission dynamics of the Middle East respiratory syndrome (MERS) outbreak in the Republic of Korea, 2015: a retrospective epidemiological analysis</article-title>. <source>BMJ Open</source><volume>6</volume>(<issue>2</issue>), <fpage>009936</fpage> (<year>2016</year>).</mixed-citation></ref><ref id="pntd.0008065.ref034"><label>34</label><mixed-citation publication-type="journal"><name><surname>Chan</surname><given-names>R.W.</given-names></name><etal>et al</etal><article-title>Tropism and replication of Middle East respiratory syndrome coronavirus from dromedary camels in the human respiratory tract: an in-vitro and ex-vivo study</article-title>. <source>Lancet Respir. Med</source>. <volume>2</volume>(<issue>10</issue>), <fpage>813</fpage>–<lpage>822</lpage> (<year>2014</year>). <pub-id pub-id-type="doi">10.1016/S2213-2600(14)70158-4</pub-id><pub-id pub-id-type="pmid">25174549</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref035"><label>35</label><mixed-citation publication-type="journal"><name><surname>Kucharski</surname><given-names>A.J.</given-names></name> &<name><surname>Althaus</surname><given-names>C</given-names></name>. <article-title>The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission</article-title>. <source>Eurosurveillance</source><volume>20</volume>(<issue>25</issue>), <fpage>21167</fpage> (<year>2015</year>). <pub-id pub-id-type="doi">10.2807/1560-7917.es2015.20.25.21167</pub-id><pub-id pub-id-type="pmid">26132768</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref036"><label>36</label><mixed-citation publication-type="journal"><name><surname>Chowell</surname><given-names>G.</given-names></name>, <name><surname>Abdirizak</surname><given-names>F.</given-names></name>, <name><surname>Lee</surname><given-names>S.</given-names></name>, <name><surname>Lee</surname><given-names>J.</given-names></name>, <name><surname>Jung</surname><given-names>E.</given-names></name>, <name><surname>Nishiura</surname><given-names>H.</given-names></name> &<name><surname>Viboud</surname><given-names>C</given-names></name>. <article-title>Transmission characteristics of MERS and SARS in the healthcare setting: a comparative study</article-title>. <source>BMC Med</source>. <volume>13</volume>(<issue>1</issue>), <fpage>210</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26336062</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref037"><label>37</label><mixed-citation publication-type="journal"><name><surname>Cotten</surname><given-names>M.</given-names></name><etal>et al</etal><article-title>Spread, circulation, and evolution of the Middle East respiratory syndrome coronavirus</article-title>. <source>MBio</source>, <volume>5</volume>(<issue>1</issue>), pii: <fpage>e01062</fpage>–<lpage>13</lpage> (<year>2014</year>). <pub-id pub-id-type="doi">10.1128/mBio.01062-13</pub-id><pub-id pub-id-type="pmid">24549846</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref038"><label>38</label><mixed-citation publication-type="journal"><name><surname>Waters</surname><given-names>E.K.</given-names></name><article-title>Modelling crowding effects in infectious disease transmission</article-title>. <source>B. Aust. Math. Soc</source>. <volume>92</volume>(<issue>3</issue>),<fpage>522</fpage>–<lpage>523</lpage> (<year>2015</year>).</mixed-citation></ref><ref id="pntd.0008065.ref039"><label>39</label><mixed-citation publication-type="journal"><name><surname>Lloyd</surname><given-names>M.</given-names></name><article-title>Mean crowding</article-title>. <source>J. Anim. Ecol</source>. <fpage>1</fpage>–<lpage>30</lpage> (<year>1967</year>).</mixed-citation></ref><ref id="pntd.0008065.ref040"><label>40</label><mixed-citation publication-type="other">World Health Organization (WHO), WHO MERS-CoV Global Summary and Assessment of Risk;<ext-link ext-link-type="uri" xlink:href="http://www.who.int/emergencies/mers-cov/risk-assessment-july-2017.pdf">http://www.who.int/emergencies/mers-cov/risk-assessment-july-2017.pdf</ext-link> (2017).</mixed-citation></ref><ref id="pntd.0008065.ref041"><label>41</label><mixed-citation publication-type="journal"><name><surname>Mukandavire</surname><given-names>Z.</given-names></name><etal>et al.</etal><article-title>Estimating the reproductive numbers for the 2008–2009 cholera outbreaks in Zimbabwe</article-title>. <source>Proc. Natl. Acad. Sci. USA</source><volume>108</volume>(<issue>21</issue>), <fpage>8767</fpage>–<lpage>8772</lpage> (<year>2011</year>). <pub-id pub-id-type="doi">10.1073/pnas.1019712108</pub-id><pub-id pub-id-type="pmid">21518855</pub-id></mixed-citation></ref><ref id="pntd.0008065.ref042"><label>42</label><mixed-citation publication-type="journal"><name><surname>Johnson</surname><given-names>L.R.</given-names></name><etal>et al</etal><article-title>Phenomenological forecasting of disease incidence using heteroskedastic Gaussian processes: a dengue case study</article-title>.<source>Ann. Appl. Stat</source>. <volume>12</volume>(<issue>1</issue>), <fpage>27</fpage>–<lpage>66</lpage> (<year>2018</year>).</mixed-citation></ref><ref id="pntd.0008065.ref043"><label>43</label><mixed-citation publication-type="journal"><name><surname>Haario</surname><given-names>H.</given-names></name>, <name><surname>Laine</surname><given-names>M.</given-names></name>, <name><surname>Mira</surname><given-names>A.</given-names></name> &<name><surname>Saksman</surname><given-names>E</given-names></name>. <article-title>DRAM: efficient adaptive MCMC</article-title>. <source>Stat.Comput</source>. <volume>16</volume>(<issue>4</issue>), <fpage>339</fpage>–<lpage>54</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="pntd.0008065.ref044"><label>44</label><mixed-citation publication-type="other">World Health Organization (WHO), Disease outbreak news MERS-CoV; <ext-link ext-link-type="uri" xlink:href="http://www.who.int/csr/don/archive/disease/coronavirus_infections/en/">http://www.who.int/csr/don/archive/disease/coronavirus_infections/en/</ext-link> (2018).</mixed-citation></ref><ref id="pntd.0008065.ref045"><label>45</label><mixed-citation publication-type="other">World Health Organization (WHO), MERS-CoV situation; <ext-link ext-link-type="uri" xlink:href="http://www.emro.who.int/pandemic-epidemic-diseases/mers-cov/mers-situation-update-july-2018.html">http://www.emro.who.int/pandemic-epidemic-diseases/mers-cov/mers-situation-update-july-2018.html</ext-link> (2018).</mixed-citation></ref><ref id="pntd.0008065.ref046"><label>46</label><mixed-citation publication-type="journal"><name><surname>Lin</surname><given-names>Q.</given-names></name>, <name><surname>Chiu</surname><given-names>A.P.Y</given-names></name>, <name><surname>Zhao</surname><given-names>S</given-names></name>., & <name><surname>He</surname><given-names>D</given-names></name>. <article-title>Modeling the spread of Middle East respiratory syndrome coronavirus in Saudi Arabia</article-title>. <source>Stat Methods Med Res</source>. <volume>27</volume>(<issue>7</issue>), <fpage>1968</fpage>–<lpage>1978</lpage> (<year>2018</year>). <pub-id pub-id-type="doi">10.1177/0962280217746442</pub-id><pub-id pub-id-type="pmid">29846148</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Sante/explor/CovidV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000978 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000978 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Sante
|area= CovidV2
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:7046297
|texte= A realistic two-strain model for MERS-CoV infection uncovers the high risk for epidemic propagation
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:32059047" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a CovidV2
| This area was generated with Dilib version V0.6.33. |  |