MERS-CoV recombination: implications about the reservoir and potential for adaptation
Identifieur interne : 001187 ( Pmc/Corpus ); précédent : 001186; suivant : 001188MERS-CoV recombination: implications about the reservoir and potential for adaptation
Auteurs : Gytis Dudas ; Andrew RambautSource :
- Virus Evolution [ 2057-1577 ] ; 2016.
Abstract
Recombination is a process that unlinks neighboring loci allowing for independent evolutionary trajectories within genomes of many organisms. If not properly accounted for, recombination can compromise many evolutionary analyses. In addition, when dealing with organisms that are not obligately sexually reproducing, recombination gives insight into the rate at which distinct genetic lineages come into contact. Since June 2012, Middle East respiratory syndrome coronavirus (MERS-CoV) has caused 1,106 laboratory-confirmed infections, with 421 MERS-CoV-associated deaths as of 16 April 2015. Although bats are considered as the likely ultimate source of zoonotic betacoronaviruses, dromedary camels have been consistently implicated as the source of current human infections in the Middle East. In this article, we use phylogenetic methods and simulations to show that MERS-CoV genome has likely undergone numerous recombinations recently. Recombination in MERS-CoV implies frequent co-infection with distinct lineages of MERS-CoV, probably in camels given the current understanding of MERS-CoV epidemiology.
Url:
DOI: 10.1093/ve/vev023
PubMed: 27774293
PubMed Central: 4989901
Links to Exploration step
PMC:4989901Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">MERS-CoV recombination: implications about the reservoir and potential for adaptation</title><author><name sortKey="Dudas, Gytis" sort="Dudas, Gytis" uniqKey="Dudas G" first="Gytis" last="Dudas">Gytis Dudas</name><affiliation><nlm:aff id="vev023-AFF1">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK,</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="vev023-AFF1">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK,</nlm:aff></affiliation><affiliation><nlm:aff wicri:cut=" and" id="vev023-AFF2">Centre for Immunology, Infection and Evolution at the University of Edinburgh, Edinburgh, UK</nlm:aff></affiliation><affiliation><nlm:aff id="vev023-AFF3">Fogarty International Center, National Institutes of Health, Bethesda, MD, USA</nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">27774293</idno><idno type="pmc">4989901</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4989901</idno><idno type="RBID">PMC:4989901</idno><idno type="doi">10.1093/ve/vev023</idno><date when="2016">2016</date><idno type="wicri:Area/Pmc/Corpus">001187</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001187</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">MERS-CoV recombination: implications about the reservoir and potential for adaptation</title><author><name sortKey="Dudas, Gytis" sort="Dudas, Gytis" uniqKey="Dudas G" first="Gytis" last="Dudas">Gytis Dudas</name><affiliation><nlm:aff id="vev023-AFF1">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK,</nlm:aff></affiliation></author><author><name sortKey="Rambaut, Andrew" sort="Rambaut, Andrew" uniqKey="Rambaut A" first="Andrew" last="Rambaut">Andrew Rambaut</name><affiliation><nlm:aff id="vev023-AFF1">Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK,</nlm:aff></affiliation><affiliation><nlm:aff wicri:cut=" and" id="vev023-AFF2">Centre for Immunology, Infection and Evolution at the University of Edinburgh, Edinburgh, UK</nlm:aff></affiliation><affiliation><nlm:aff id="vev023-AFF3">Fogarty International Center, National Institutes of Health, Bethesda, MD, USA</nlm:aff></affiliation></author></analytic><series><title level="j">Virus Evolution</title><idno type="eISSN">2057-1577</idno><imprint><date when="2016">2016</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Recombination is a process that unlinks neighboring loci allowing for independent evolutionary trajectories within genomes of many organisms. If not properly accounted for, recombination can compromise many evolutionary analyses. In addition, when dealing with organisms that are not obligately sexually reproducing, recombination gives insight into the rate at which distinct genetic lineages come into contact. Since June 2012, Middle East respiratory syndrome coronavirus (MERS-CoV) has caused 1,106 laboratory-confirmed infections, with 421 MERS-CoV-associated deaths as of 16 April 2015. Although bats are considered as the likely ultimate source of zoonotic betacoronaviruses, dromedary camels have been consistently implicated as the source of current human infections in the Middle East. In this article, we use phylogenetic methods and simulations to show that MERS-CoV genome has likely undergone numerous recombinations recently. Recombination in MERS-CoV implies frequent co-infection with distinct lineages of MERS-CoV, probably in camels given the current understanding of MERS-CoV epidemiology.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Aburizaiza, A S" uniqKey="Aburizaiza A">A. S. Aburizaiza</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Adney, D R" uniqKey="Adney D">D. R. Adney</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Anisimova, M" uniqKey="Anisimova M">M. Anisimova</name></author><author><name sortKey="Nielsen, R" uniqKey="Nielsen R">R. Nielsen</name></author><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z. Yang</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="Baele, G" uniqKey="Baele G">G. Baele</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Baric, R S" uniqKey="Baric R">R. S. Baric</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bielejec, F" uniqKey="Bielejec F">F. Bielejec</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boni, M F" uniqKey="Boni M">M. F. Boni</name></author><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author><author><name sortKey="Feldman, M W" uniqKey="Feldman M">M. W. Feldman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boni, M F" uniqKey="Boni M">M. F. Boni</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Boucher, C A" uniqKey="Boucher C">C. A. Boucher</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Bruen, T C" uniqKey="Bruen T">T. C. Bruen</name></author><author><name sortKey="Philippe, H" uniqKey="Philippe H">H. Philippe</name></author><author><name sortKey="Bryant, D" uniqKey="Bryant D">D. Bryant</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="Chang, S F" uniqKey="Chang S">S. F. Chang</name></author><author><name sortKey="Sgro, J Y" uniqKey="Sgro J">J. Y. Sgro</name></author><author><name sortKey="Parrish, C R" uniqKey="Parrish C">C. R. Parrish</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chare, E R" uniqKey="Chare E">E. R. Chare</name></author><author><name sortKey="Gould, E A" uniqKey="Gould E">E. A. Gould</name></author><author><name sortKey="Holmes, E C" uniqKey="Holmes E">E. C. Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chu, D K" uniqKey="Chu D">D. K. Chu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Corman, V M" uniqKey="Corman V">V. M. Corman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Corman, V M" uniqKey="Corman V">V. M. Corman</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="Cotten, M" uniqKey="Cotten M">M. Cotten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Didelot, X" uniqKey="Didelot X">X. Didelot</name></author><author><name sortKey="Falush, D" uniqKey="Falush D">D. Falush</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Dong, G" uniqKey="Dong G">G. Dong</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drosten, C" uniqKey="Drosten C">C. Drosten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drosten, C" uniqKey="Drosten C">C. Drosten</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Excoffier, L" uniqKey="Excoffier L">L. Excoffier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ferreira, M A R" uniqKey="Ferreira M">M. A. R. Ferreira</name></author><author><name sortKey="Suchard, M A" uniqKey="Suchard M">M. A. Suchard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gierer, S" uniqKey="Gierer S">S. Gierer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gill, M S" uniqKey="Gill M">M. S. Gill</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gire, S K" uniqKey="Gire S">S. K. Gire</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Goldman, N" uniqKey="Goldman N">N. Goldman</name></author><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gubareva, L V" uniqKey="Gubareva L">L. V. Gubareva</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S. Guindon</name></author><author><name sortKey="Gascuel, O" uniqKey="Gascuel O">O. Gascuel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hasegawa, M" uniqKey="Hasegawa M">M. Hasegawa</name></author><author><name sortKey="Kishino, H" uniqKey="Kishino H">H. Kishino</name></author><author><name sortKey="Yano, T A" uniqKey="Yano T">T. A. Yano</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hedrick, P W" uniqKey="Hedrick P">P. W. Hedrick</name></author><author><name sortKey="Thomson, G" uniqKey="Thomson G">G. Thomson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Herrewegh, A A" uniqKey="Herrewegh A">A. A. Herrewegh</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hill, W G" uniqKey="Hill W">W. G. Hill</name></author><author><name sortKey="Robertson, A" uniqKey="Robertson A">A. Robertson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Holmes, E C" uniqKey="Holmes E">E. C. Holmes</name></author><author><name sortKey="Urwin, R" uniqKey="Urwin R">R. Urwin</name></author><author><name sortKey="Maiden, M C" uniqKey="Maiden M">M. C. Maiden</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hon, C C" uniqKey="Hon C">C. C. Hon</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Keck, J G" uniqKey="Keck J">J. G. Keck</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kirkegaard, K" uniqKey="Kirkegaard K">K. Kirkegaard</name></author><author><name sortKey="Baltimore, D" uniqKey="Baltimore D">D. Baltimore</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kosakovsky Pond, S L" uniqKey="Kosakovsky Pond S">S. L. Kosakovsky Pond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kottier, S A" uniqKey="Kottier S">S. A. Kottier</name></author><author><name sortKey="Cavanagh, D" uniqKey="Cavanagh D">D. Cavanagh</name></author><author><name sortKey="Britton, P" uniqKey="Britton P">P. Britton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lai, M M" uniqKey="Lai M">M. M. Lai</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lau, S K P" uniqKey="Lau S">S. K. P. Lau</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lemey, P" uniqKey="Lemey P">P. Lemey</name></author><author><name sortKey="Suchard, M" uniqKey="Suchard M">M. Suchard</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lewontin, R C" uniqKey="Lewontin R">R. C. Lewontin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, K S" uniqKey="Li K">K. S. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lu, L" uniqKey="Lu L">L. Lu</name></author><author><name sortKey="Lycett, S J" uniqKey="Lycett S">S. J. Lycett</name></author><author><name sortKey="Brown, A J L" uniqKey="Brown A">A. J. L. Brown</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Makino, S" uniqKey="Makino S">S. Makino</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Maynard Smith, J" uniqKey="Maynard Smith J">J. Maynard Smith</name></author><author><name sortKey="Smith, N H" uniqKey="Smith N">N. H. Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mcvean, G" uniqKey="Mcvean G">G. McVean</name></author><author><name sortKey="Awadalla, P" uniqKey="Awadalla P">P. Awadalla</name></author><author><name sortKey="Fearnhead, P" uniqKey="Fearnhead P">P. Fearnhead</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Memish, Z A" uniqKey="Memish Z">Z. A. Memish</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meunier, J" uniqKey="Meunier J">J. Meunier</name></author><author><name sortKey="Eyre Walker, A" uniqKey="Eyre Walker A">A. Eyre-Walker</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Milkman, R" uniqKey="Milkman R">R. Milkman</name></author><author><name sortKey="Bridges, M M" uniqKey="Bridges M">M. M. Bridges</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mitnaul, L J" uniqKey="Mitnaul L">L. J. Mitnaul</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muller, H J" uniqKey="Muller H">H. J. Muller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muller, M A" uniqKey="Muller M">M. A. Müller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muller, M A" uniqKey="Muller M">M. A. Müller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="O Rien, J D" uniqKey="O Rien J">J. D. O’Brien</name></author><author><name sortKey="Minin, V N" uniqKey="Minin V">V. N. Minin</name></author><author><name sortKey="Suchard, M A" uniqKey="Suchard M">M. A. Suchard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pond, S L K" uniqKey="Pond S">S. L. K. Pond</name></author><author><name sortKey="Frost, S D W" uniqKey="Frost S">S. D. W. Frost</name></author><author><name sortKey="Muse, S V" uniqKey="Muse S">S. V. Muse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author><author><name sortKey="Crandall, K A" uniqKey="Crandall K">K. A. Crandall</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author><author><name sortKey="Crandall, K A" uniqKey="Crandall K">K. A. Crandall</name></author><author><name sortKey="Crandall, K A" uniqKey="Crandall K">K. A. Crandall</name></author><author><name sortKey="Holmes, E C" uniqKey="Holmes E">E. C. Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author><author><name sortKey="Holmes, E" uniqKey="Holmes E">E. Holmes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Reusken, C B" uniqKey="Reusken C">C. B. Reusken</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Reusken, C B" uniqKey="Reusken C">C. B. Reusken</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robertson, D L" uniqKey="Robertson D">D. L. Robertson</name></author><author><name sortKey="Hahn, B H" uniqKey="Hahn B">B. H. Hahn</name></author><author><name sortKey="Sharp, P M" uniqKey="Sharp P">P. M. Sharp</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Robertson, D L" uniqKey="Robertson D">D. L. Robertson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schierup, M H" uniqKey="Schierup M">M. H. Schierup</name></author><author><name sortKey="Hein, J" uniqKey="Hein J">J. Hein</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shriner, D" uniqKey="Shriner D">D. Shriner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Smith, G J D" uniqKey="Smith G">G. J. D. Smith</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Subbarao, E K" uniqKey="Subbarao E">E. K. Subbarao</name></author><author><name sortKey="London, W" uniqKey="London W">W. London</name></author><author><name sortKey="Murphy, B R" uniqKey="Murphy B">B. R. Murphy</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tavare, S" uniqKey="Tavare S">S. Tavaré</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tisdale, M" uniqKey="Tisdale M">M. Tisdale</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, Z" uniqKey="Yang Z">Z. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zaki, A M" uniqKey="Zaki A">A. M. Zaki</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">Virus Evol</journal-id><journal-id journal-id-type="iso-abbrev">Virus Evol</journal-id><journal-id journal-id-type="publisher-id">vevolu</journal-id><journal-id journal-id-type="hwp">vevolu</journal-id><journal-title-group><journal-title>Virus Evolution</journal-title></journal-title-group><issn pub-type="epub">2057-1577</issn><publisher><publisher-name>Oxford University Press</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">27774293</article-id><article-id pub-id-type="pmc">4989901</article-id><article-id pub-id-type="doi">10.1093/ve/vev023</article-id><article-id pub-id-type="publisher-id">vev023</article-id><article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject></subj-group></article-categories><title-group><article-title>MERS-CoV recombination: implications about the reservoir and potential for adaptation</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Dudas</surname><given-names>Gytis</given-names></name><xref ref-type="aff" rid="vev023-AFF1"><sup>1</sup></xref><xref ref-type="corresp" rid="vev023-COR1">*</xref></contrib><contrib contrib-type="author"><name><surname>Rambaut</surname><given-names>Andrew</given-names></name><xref ref-type="aff" rid="vev023-AFF1"><sup>1</sup></xref><xref ref-type="aff" rid="vev023-AFF2"><sup>2</sup></xref><xref ref-type="aff" rid="vev023-AFF3"><sup>3</sup></xref><xref ref-type="corresp" rid="vev023-COR2"><sup>†</sup></xref></contrib><aff id="vev023-AFF1"><sup>1</sup>Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK,</aff><aff id="vev023-AFF2"><sup>2</sup>Centre for Immunology, Infection and Evolution at the University of Edinburgh, Edinburgh, UK and</aff><aff id="vev023-AFF3"><sup>3</sup>Fogarty International Center, National Institutes of Health, Bethesda, MD, USA</aff></contrib-group><author-notes><corresp id="vev023-COR1">*Corresponding author: E-mail: <email>g.dudas@sms.ed.ac.uk</email></corresp><corresp id="vev023-COR2"><sup>†</sup><ext-link ext-link-type="uri" xlink:href="http://orcid.org/0000-0003-4337-3707">http://orcid.org/0000-0003-4337-3707</ext-link></corresp></author-notes><pub-date pub-type="collection"><month>1</month><year>2016</year></pub-date><pub-date pub-type="epub"><day>20</day><month>1</month><year>2016</year></pub-date><volume>2</volume><issue>1</issue><elocation-id>vev023</elocation-id><permissions><copyright-statement>© The Author 2016. Published by Oxford University Press.</copyright-statement><copyright-year>2016</copyright-year><license xlink:href="http://creativecommons.org/licenses/by-nc/4.0/" license-type="creative-commons"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by-nc/4.0/">http://creativecommons.org/licenses/by-nc/4.0/</ext-link>), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com</license-p></license></permissions><abstract><p>Recombination is a process that unlinks neighboring loci allowing for independent evolutionary trajectories within genomes of many organisms. If not properly accounted for, recombination can compromise many evolutionary analyses. In addition, when dealing with organisms that are not obligately sexually reproducing, recombination gives insight into the rate at which distinct genetic lineages come into contact. Since June 2012, Middle East respiratory syndrome coronavirus (MERS-CoV) has caused 1,106 laboratory-confirmed infections, with 421 MERS-CoV-associated deaths as of 16 April 2015. Although bats are considered as the likely ultimate source of zoonotic betacoronaviruses, dromedary camels have been consistently implicated as the source of current human infections in the Middle East. In this article, we use phylogenetic methods and simulations to show that MERS-CoV genome has likely undergone numerous recombinations recently. Recombination in MERS-CoV implies frequent co-infection with distinct lineages of MERS-CoV, probably in camels given the current understanding of MERS-CoV epidemiology.</p></abstract><kwd-group><kwd>MERS-CoV</kwd><kwd>recombination</kwd><kwd>co-infection</kwd><kwd>homoplasy</kwd><kwd>coronavirus</kwd><kwd>MERS</kwd></kwd-group><counts><page-count count="11"></page-count></counts></article-meta></front><body><sec><title>1 Introduction</title><p>Recombination is an important process which expedites selection in many organisms (<xref rid="vev023-B57" ref-type="bibr">Muller 1932</xref>) by unlinking loci. It also leads to different parts of recombining genomes to have different histories which, if not properly accounted for, can interfere with many genetic analyses, of which phylogenetic methods are among the most sensitive. Not accounting for recombination in phylogenetic analyses leads to incorrect (<xref rid="vev023-B70" ref-type="bibr">Schierup and Hein 2000</xref>) and poorly supported genealogies (<xref rid="vev023-B63" ref-type="bibr">Posada and Crandall 2002</xref>) and false inference of selection (<xref rid="vev023-B3" ref-type="bibr">Anisimova, Nielsen, and Yang 2003</xref>; <xref rid="vev023-B71" ref-type="bibr">Shriner et al. 2003</xref>).</p><p>With rising sequence availability during outbreaks of viral infectious disease, phylogenetic methods have been used to supplement our knowledge of epidemics in real time (<xref rid="vev023-B46" ref-type="bibr">Lemey, Suchard, and Rambaut 2009</xref>; <xref rid="vev023-B65" ref-type="bibr">Rambaut and Holmes 2009</xref>; <xref rid="vev023-B72" ref-type="bibr">Smith et al. 2009</xref>; <xref rid="vev023-B23" ref-type="bibr">Drosten et al. 2013</xref>, <xref rid="vev023-B22" ref-type="bibr">2014</xref>; <xref rid="vev023-B18" ref-type="bibr">Cotten et al. 2013</xref>, <xref rid="vev023-B19" ref-type="bibr">2014</xref>; <xref rid="vev023-B30" ref-type="bibr">Gire et al. 2014</xref>). For some outbreaks, there is little reason to suspect recombination, e.g. negative sense single-stranded RNA viruses are thought to recombine over evolutionary, not population-level, time scales (<xref rid="vev023-B14" ref-type="bibr">Chare, Gould, and Holmes 2003</xref>). Observable recombination in RNA viruses requires that two conditions are met that viruses from distinct lineages co-infect a host and that a mechanism for recombination exists. For example, even though influenza A virus co-infection is extremely common in birds based on genome segment reassortment patterns (<xref rid="vev023-B48" ref-type="bibr">Li et al. 2004</xref>; <xref rid="vev023-B21" ref-type="bibr">Dong et al. 2011</xref>; <xref rid="vev023-B49" ref-type="bibr">Lu, Lycett, and Brown 2014</xref>), recombination is extremely rare or absent (<xref rid="vev023-B14" ref-type="bibr">Chare, Gould, and Holmes 2003</xref>; <xref rid="vev023-B9" ref-type="bibr">Boni et al. 2010</xref>). This is thought to be because template switching (<xref rid="vev023-B41" ref-type="bibr">Kirkegaard and Baltimore 1986</xref>; <xref rid="vev023-B6" ref-type="bibr">Baric et al. 1987</xref>), the main mechanism of recombination in RNA viruses, is mechanistically difficult for single-stranded negative sense RNA viruses (<xref rid="vev023-B14" ref-type="bibr">Chare, Gould, and Holmes 2003</xref>), and for influenza A viruses has only been convincingly shown in cell culture under extreme conditions (<xref rid="vev023-B56" ref-type="bibr">Mitnaul et al. 2000</xref>). When the genomic architecture of a virus is permissive to recombination, i.e. template switching occurs and is detectable, the extent of recombination is informative of co-infection and/or duration of infection.</p><p>Here, we focus our attention on the Middle East respiratory syndrome coronavirus (MERS-CoV) (<xref rid="vev023-B77" ref-type="bibr">Zaki et al. 2012</xref>), a recent zoonotic infection with a relatively high case fatality ratio (<xref rid="vev023-B4" ref-type="bibr">Assiri et al. 2013</xref>; <xref rid="vev023-B53" ref-type="bibr">Memish et al. 2013</xref>; <xref rid="vev023-B12" ref-type="bibr">Cauchemez et al. 2014</xref>). Most human infections with MERS-CoV are thought to be the result of contact with <italic>Camelus dromedarius</italic> L., the dromedary camel, which is the presumed host of the virus. MERS-CoV, much like severe acute respiratory syndrome coronavirus (SARS-CoV), is likely ultimately derived from bats (<xref rid="vev023-B16" ref-type="bibr">Corman et al. 2014a</xref>). MERS-CoV, along with Murine hepatitis virus and SARS-CoV, belongs to the Betacoronavirus genus. Betacoronavirus as well as two other genera (Alpha- and Gammacoronavirus) out of four within the subfamily <italic>Coronavirinae</italic> have been shown to recombine in cell culture, <italic>in vivo</italic><italic>,</italic> and in eggs (<xref rid="vev023-B44" ref-type="bibr">Lai et al. 1985</xref>; <xref rid="vev023-B50" ref-type="bibr">Makino et al. 1986</xref>; <xref rid="vev023-B40" ref-type="bibr">Keck et al. 1988</xref>; <xref rid="vev023-B43" ref-type="bibr">Kottier, Cavanagh, and Britton 1995</xref>; <xref rid="vev023-B36" ref-type="bibr">Herrewegh et al. 1998</xref>). Additionally, a coronavirus lineage related to MERS-CoV which was isolated from bats appears to have recombined around the spike (S) protein (<xref rid="vev023-B16" ref-type="bibr">Corman et al. 2014a</xref>). In this article, we show that although the genome of MERS-CoV contains considerable amounts of rate heterogeneity between genomic regions that can interfere with detection of recombination, we do nonetheless find evidence of sustained recombination that cannot be explained by rate heterogeneity alone. This has two important consequences: one is that care has to be taken when constructing phylogenetic trees of MERS-CoV as a single tree cannot accurately describe the complete history of all loci within a recombining genome. Secondly and more importantly, the observed rates of recombination in the MERS-CoV genome is evidence of a large number of MERS-CoV co-infections in some hosts which has implications for understanding the dynamics of the virus in the animal reservoir.</p></sec><sec><title>2 Methods</title><sec sec-type="intro"><title>2.1 Overview</title><p>Recombination leaves several characteristic clues in genomes:
<list list-type="bullet"><list-item><p>Alternative topologies (<xref rid="vev023-B68" ref-type="bibr">Robertson, Hahn, and Sharp 1995</xref>; <xref rid="vev023-B68" ref-type="bibr">Robertson et al. 1995</xref>; <xref rid="vev023-B38" ref-type="bibr">Holmes, Urwin, and Maiden 1999</xref>). In some scenarios, for example, if there has been a single large-scale recombination event, it is possible to clearly identify recombining fragments based on phylogenetic incongruity. Recombination can be inferred by reconstructing two or more phylogenetic trees from a partitioned alignment and looking for topological incongruity between them. Strong support for at least two incompatible phylogenetic trees across well-defined breakpoints is usually the most convincing evidence of recombination.</p></list-item><list-item><p>Excessive homoplasies (<xref rid="vev023-B51" ref-type="bibr">Maynard Smith and Smith 1998</xref>). The transfer of genetic material from one genetic background to another will result in apparent repeat mutations in different parts of a phylogenetic tree. However, it is possible for the same locus to undergo mutation independently, especially if the locus in question is under Darwinian selection. Detecting homoplasies alone is not sufficient to infer recombination but should be demonstrated to occur in excess of expectation.</p></list-item><list-item><p>Linkage disequilibrium (LD) decay (<xref rid="vev023-B54" ref-type="bibr">Meunier and Eyre-Walker 2001</xref>). LD is the non-random association of alleles at different loci. This is a statistic often reported for contemporaneous sequence data. In clonally (i.e. non-recombining) evolving organisms, every allele is linked to every other allele in the genome and requires mutation to break linkage. In recombining organisms, there is an expectation that LD will decay with distance between the loci, i.e. that loci further away from each other are more likely to be unlinked via recombination.</p></list-item></list></p><p>We test for each of these hallmarks of recombination in the MERS coronavirus genome using a combination of phylogenetic and LD metrics. For a more detailed review of recombination detection methods, see <xref rid="vev023-B63" ref-type="bibr">Posada, Crandall, and Holmes (2002)</xref>.</p></sec><sec><title>2.2 Alternative topologies</title><p>We use the Genetic Algorithm for Recombination Detection (GARD) method (<xref rid="vev023-B42" ref-type="bibr">Kosakovsky Pond et al. 2006),</xref> as implemented in the software package HyPhy (<xref rid="vev023-B61" ref-type="bibr">Pond, Frost, and Muse 2005</xref>), to look for alternative tree topologies in sequence data. Briefly, the method compares a model where a single tree is derived from the whole alignment and alternative models where breakpoints are introduced into the alignment and phylogenetic trees are derived independently from the resulting fragments. The presence of recombination, especially if it is recent and concentrated in some parts of the alignment, will result in two or more phylogenetic trees fitting the data better than a single tree model. It is important to note that likelihood estimation also involves other parameters, such as branch lengths, not just topology. We use GARD under a GTR (<xref rid="vev023-B74" ref-type="bibr">Tavaré 1986</xref>) substitution model with <inline-formula><mml:math id="MM1"><mml:mrow><mml:msub><mml:mi>Γ</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>-distributed rate heterogeneity among sites (<xref rid="vev023-B76" ref-type="bibr">Yang 1994</xref>) on a dataset of eighty-five MERS-CoV sequences. GARD was run repeatedly until no more breakpoints could be identified in the resulting fragments.</p><p>In addition to this test, we run BEAST (<xref rid="vev023-B25" ref-type="bibr">Drummond et al. 2012</xref>) on partitioned coding sequences derived from the first well-supported breakpoint inferred by GARD. We extracted the coding sequences from nucleotide positions 1–23722 and 23723–30126 (referred to as Fragments 1 and 2, respectively) of MERS-CoV genomes. Independent HKY+<inline-formula><mml:math id="MM2"><mml:mrow><mml:msub><mml:mi>Γ</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> (<xref rid="vev023-B34" ref-type="bibr">Hasegawa, Kishino, and Yano 1985</xref>; <xref rid="vev023-B76" ref-type="bibr">Yang 1994</xref>) nucleotide substitution models were specified for codon positions 1+2 and 3, and the analyses were run under an uncorrelated relaxed lognormal clock with an uninformative CTMC reference prior (<xref rid="vev023-B27" ref-type="bibr">Ferreira and Suchard 2008</xref>) on the substitution rate for 100 million states, subsampling every 10,000 states. The molecular clocks and trees of each genomic partition were either linked or unlinked, giving a total of four models. We used the multi-locus skygrid (<xref rid="vev023-B29" ref-type="bibr">Gill et al. 2013</xref>) as the demographic model for all analyses. Path-sampling and stepping stone sampling (<xref rid="vev023-B5" ref-type="bibr">Baele et al. 2012</xref>) were used to calculate marginal likelihoods and test the fit of each of the four models, under default parameters. In addition, four similar analyses were set up, but with strict molecular clocks, to contrast the performance of relaxed molecular clocks.</p></sec><sec><title>2.3 Excessive homoplasies</title><p>Testing for recombination by looking for homoplasic mutations in phylogenetic trees requires that two conditions are met. One, that recombination is rare enough, so that there is sufficient phylogenetic signal to reconstruct the ‘correct’ phylogeny otherwise known as the clonal frame (<xref rid="vev023-B55" ref-type="bibr">Milkman and Bridges 1990</xref>). Two, that alternative explanations for homoplasic mutations can be dismissed with some certainty. There is no straightforward way of testing for the former, but the latter is usually dictated by the underlying biology. For example, repeat amino acid substitutions are a well-documented response of influenza viruses and HIV to drug treatment (<xref rid="vev023-B10" ref-type="bibr">Boucher et al. 1993</xref>; <xref rid="vev023-B75" ref-type="bibr">Tisdale et al. 1993</xref>; <xref rid="vev023-B32" ref-type="bibr">Gubareva et al. 2001</xref>).</p><p>We employ two methods to test for excessive homoplasies. First, we use a maximum likelihood phylogeny inferred using PhyML (<xref rid="vev023-B33" ref-type="bibr">Guindon and Gascuel 2003</xref>) under a GTR+<inline-formula><mml:math id="MM3"><mml:mrow><mml:msub><mml:mi>Γ</mml:mi><mml:mn>4</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> (<xref rid="vev023-B74" ref-type="bibr">Tavaré 1986</xref>; <xref rid="vev023-B76" ref-type="bibr">Yang 1994</xref>) nucleotide substitution model to recover a single tree using a MERS-CoV dataset comprising eighty-five sequences. We then reconstruct ancestral sequences at each internal node and identify the mutations that have taken place along each branch using ClonalFrameML (<xref rid="vev023-B20" ref-type="bibr">Didelot and Falush 2007</xref>). Mutations are then classified as either synapomorphies, shared variation derived via common descent, or apparent homoplasies, shared variation derived from convergence, depending on how many times a given mutation has arisen in the phylogeny. The drawback of this method is that it necessarily conditions on a single tree with the highest likelihood.</p><p>We also employ BEAST (<xref rid="vev023-B25" ref-type="bibr">Drummond et al. 2012</xref>) to circumvent the limitation of conditioning the ancestral state reconstruction on a single tree. In addition to sampling various phylogenetic parameters from the posterior distribution, BEAST is also able to map substitutions onto the branches of each MCMC-sampled phylogeny (<xref rid="vev023-B60" ref-type="bibr">O’Brien, Minin, and Suchard 2009</xref>). This method is thus capable of estimating the posterior probability of a given mutation being synapomorphic or homoplasic by integrating over different tree topologies. Homoplasy analyses were performed on the concatenated coding sequences of MERS-CoV after partitioning the alignment into all three codon positions, each with an HKY nucleotide substitution model (<xref rid="vev023-B34" ref-type="bibr">Hasegawa, Kishino, and Yano 1985</xref>) and no Γ-distributed rate heterogeneity among sites. A relaxed uncorrelated molecular clock with lognormally distributed rates (<xref rid="vev023-B24" ref-type="bibr">Drummond et al. 2006</xref>) under a CTMC reference prior (<xref rid="vev023-B27" ref-type="bibr">Ferreira and Suchard 2008</xref>) and the flexible multi-locus skygrid as the demographic model (<xref rid="vev023-B29" ref-type="bibr">Gill et al. 2013</xref>) were used. The MCMC chain was run for 100 million steps, sampling every 10,000 steps.</p><p>Throughout the article, we will refer to the number of branches that have experienced a given mutation as homoplasy degree. We define the homoplasy degree to be the number of times a given mutation has originated independently minus one. For example, a homoplasy degree of 1 indicates that a mutation has occurred on two different branches in the phylogeny. That is, we assume that one of the mutations has arisen through replication error, whereas the other has potential to have been introduced via recombination and thus can be thought of as excessive. Synapomorphies, on the other hand, are states that are shared by two or more taxa through common descent and thus necessarily are those mutations that have occurred exactly once in the phylogeny. They have a homoplasy degree of 0 in all figures.</p><p>Additional tests for recombination were also performed, namely the estimation of the pairwise homoplasy index (<xref rid="vev023-B11" ref-type="bibr">Bruen, Philippe, and Bryant 2006</xref>) and the triplet test implemented in 3Seq (<xref rid="vev023-B8" ref-type="bibr">Boni, Posada, and Feldman 2007</xref>).</p></sec><sec><title>2.4 LD decay</title><p>In the absence of recombination every allele should exhibit a high degree of linkage with other alleles in the genome. Under two extremes—clonal reproduction without recombination and free recombination—there is no correlation between LD and genomic distance and loci should be interchangeable. This is the basis of several non-parametric permutation tests for recombination that are implemented in the software package LDhat (<xref rid="vev023-B52" ref-type="bibr">McVean, Awadalla, and Fearnhead 2002</xref>), which we used in combination with a dataset of 109 MERS-CoV genomes. Other, more complicated tests, such as composite likelihood methods, are also available but in our experience were incompatible with temporal sampling and rate heterogeneity.</p></sec><sec><title>2.5 Sequence simulations</title><p>To test the performance of some of the methods, we simulated two sets of sequences. We use fastsimcoal2 (<xref rid="vev023-B26" ref-type="bibr">Excoffier et al. 2013</xref>) to simulate ten replicate datasets that have the same dates of isolation and similar diversity to the MERS dataset with eighty-five sequences under no recombination.</p><p>Additionally, we use <italic>π</italic>BUSS (<xref rid="vev023-B7" ref-type="bibr">Bielejec et al. 2014</xref>) to simulate sequences down an MCMC-sampled phylogeny drawn at random from a linked-tree unlinked-clocks BEAST analysis described above. We modelled region-specific rate heterogeneity by simulating a 30 kb ‘genome’ and setting the molecular clock rate for the first 20 kb to be 9.5 × 10<sup>–</sup><sup>4</sup> substitutions site<sup>–</sup><sup>1</sup>year<sup>–</sup><sup>1</sup> and the last 10 kb to be 2.85 × 10<sup>–</sup><sup>3</sup>, 1.9 × 10<sup>–</sup><sup>3</sup> or 1.3 × 10<sup>–</sup><sup>3</sup> substitutions site<sup>–</sup><sup>1</sup>year<sup>–</sup><sup>1</sup>, corresponding to roughly 3-, 2-, or 1.3-fold rate heterogeneity between the two parts of the simulated genome. Two replicate datasets were generated for each category of rate heterogeneity. Other than that, all simulations were run under a relaxed lognormal molecular clock (<xref rid="vev023-B24" ref-type="bibr">Drummond et al. 2006</xref>) with standard deviation set to 7.42 × 10<sup>–</sup><sup>7</sup>, HKY substitution model (<xref rid="vev023-B34" ref-type="bibr">Hasegawa, Kishino, and Yano 1985</xref>) with the transition/transversion ratio parameter (<italic>κ</italic>) set to 6.0 and Γ-distributed rate heterogeneity with four categories and shape parameter 0.04 and empirical nucleotide frequencies, all derived from the results of the marginal likelihood analyses described earlier. A MERS-CoV sequence isolated from a camel (NRCE-HKU270) was provided as the starting state at the root. To include the effects of site-specific constraint, we additionally carried out simulations under a Goldman–Yang codon model (Goldman and <xref rid="vev023-B76" ref-type="bibr">Yang 1994</xref>) in <italic>π</italic>BUSS with empirical codon frequencies, <italic>κ</italic> = 6.0 and dN/dS (<italic>ω</italic>) set to 0.1 (i.e. purifying selection) under same levels of rate heterogeneity and on the same phylogeny as the simulations described above. As these sequences were simulated on a tree of MERS-CoV, we refer to these datasets as being empirically simulated. We also reconstructed ancestral states for these sequences using ClonalFrameML, as described above, to arrive at a null expectation for the number of homoplasies we expect to observe under rate heterogeneity but without recombination.</p></sec><sec><title>2.6 Investigating the effects of temporal sampling and rate heterogeneity</title><p>All ten sequence datasets simulated with fastsimcoal2 and twelve sequences empirically simulated in <italic>π</italic>BUSS were analyzed using LDhat (<xref rid="vev023-B52" ref-type="bibr">McVean, Awadalla, and Fearnhead 2002</xref>) to ascertain the effects of temporal sampling, and in the case of <italic>π</italic>BUSS-simulated sequences, the effects of rate heterogeneity in the presence of absence of position-specific constraint. Additionally, empirically simulated sequence datasets were run through GARD (<xref rid="vev023-B42" ref-type="bibr">Kosakovsky Pond et al. 2006),</xref> since the method considers both differences in tree topology and branch lengths when calculating the likelihoods of trees. Stark rate heterogeneity could thus easily be mis-interpreted as evidence for recombination by GARD.</p></sec><sec><title>2.7 Host-association alleles</title><p>To test for the presence of alleles associated with host shifts (presumably camel to human), we adapt the <inline-formula><mml:math id="MM4"><mml:mrow><mml:msubsup><mml:mi>χ</mml:mi><mml:mrow><mml:mtext>df</mml:mtext></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> (<xref rid="vev023-B35" ref-type="bibr">Hedrick and Thomson 1986</xref>) statistic of LD to estimate the association between host (camel or human) and alleles at polymorphic loci. Briefly, we consider the host to act as a polymorphic site (encoded as H or C, for human and camel, respectively) and compare the association between the ‘allele’ or host and alleles at polymorphic sites. A perfect association of 1.0 could mean, for example, that a biallelic site has one allele that is only found in camel viruses and the other allele only in human viruses.</p></sec></sec><sec><title>3 Results</title><sec><title>3.1 MERS-CoV genome shows evidence of alternative tree topologies</title><p>GARD identified a breakpoint at nucleotide position 23722 (corrected ΔAIC = 103.6 between single versus two tree model), roughly in the middle of the coding sequence for the S (spike) protein. The two phylogenies recovered from this breakpoint were incongruent (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S1</ext-link>). Running the resulting Fragment 1 (positions 1–23722) and Fragment 2 (positions 23723–30126) through GARD again yielded a further breakpoint in Fragment 1 at position 12257 (corrected ΔAIC = 33.7), near the boundary between <italic>ORF1a</italic> and <italic>ORF1b</italic> genes. No more breakpoints could be identified by GARD in the resulting Fragments 1.1 (positions 1–12257), 1.2 (positions 12258–23722) and 2 (positions 23723–30126).</p><p>Five out of six empirical simulation alignments under a nucleotide model simulated without recombination were identified by GARD as having breakpoints around position 20000, where the clock rate for the rest of the ‘genome’ was increased to be 1.3, 2, or 3 times higher than the first 20 kb. Similarly, two of six alignments empirically simulated under a codon model were identified as having a breakpoint. Corrected ΔAIC values decreased with decreasing rate heterogeneity, indicating loss of statistical power to detect differences between genomic regions. Analyses in BEAST, where the MERS-CoV genome is partitioned into positions 1–23722 and positions 23723–30126 (corresponding to the first GARD-inferred breakpoint) with each partition having an independent molecular clock rate but the same tree or both independent molecular clock rates and independent trees, showed that rate heterogeneity as expressed by the ratio of second fragment rate to first fragment rate to be on the order of 1.513 (95% highest posterior density 1.275, 1.769) for unlinked clocks and 1.375 (95% HPDs: 1.079–1.707) for unlinked clocks and trees (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S2</ext-link>). As such, empirically simulated sequence data under 2-fold and 3-fold rate heterogeneity should be considered as caricature examples of a MERS-like organism, which we use to test the sensitivity of the methods we employ.</p></sec><sec><title>3.2 MERS-CoV genome exhibits LD decay</title><p>Permutation tests as implemented in LDhat work under the assumption that loci are interchangeable only when there is free recombination or no recombination at all. The tests compare four statistics estimated from the actual data to 1,000 permutations of the data where site numbers for each locus are reshuffled. Correlation coefficient between two measures of LD, <italic>r</italic><sup>2</sup> (<xref rid="vev023-B37" ref-type="bibr">Hill and Robertson 1968</xref>) and <inline-formula><mml:math id="MM5"><mml:mrow><mml:mi>D</mml:mi><mml:mo>′</mml:mo></mml:mrow></mml:math></inline-formula> (<xref rid="vev023-B47" ref-type="bibr">Lewontin 1964</xref>), are expected to show a negative correlation with increasing distance between loci if there is recombination. Permutation of recombining loci will produce a distribution skewed toward more positive values for these two LD statistics and the percentile of the actual observed value can then be used to assess significance.</p><p>G4 is the sum of distances between pairs of loci with four observed haplotypes, which can only occur if there is repeat mutation or recombination at one of the loci. Under recombination, the observed G4 statistic should take a statistically higher value in a distribution of G4 values derived from permuted data. Lkmax is the composite likelihood of pairs of loci under an estimated recombination rate and a given level of sequence diversity. Like the G4 statistic, this statistic is expected to fall in the upper tail of the distribution derived from permuted data in the presence of recombination.</p><p>All four permutation tests show a consistent signal of recombination in the MERS-CoV genome (see <xref ref-type="fig" rid="vev023-F2">Fig. 2</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S3</ext-link>). Data from fastsimcoal2 simulations, which did not have rate heterogeneity, produced values for these statistics which mostly fell inside the range of values generated by permuting the simulated data, as expected (<xref ref-type="fig" rid="vev023-F2">Fig. 2</xref>). On one occasion, this is not the case—Simulation 9 passed the Lkmax test and failed the other three. Empirically simulated data, on the other hand, tended to exhibit extreme values, that is the observed value fell below the 2.5th or above the 97.5th percentile of the permuted data, but in ways which were not consistent with recombination. For example, Replicate 1 of simulation with 3-fold rate heterogeneity under a nucleotide substitution model exhibits extreme values for all four tests, but only one of these – corr(<italic>r</italic><sup>2</sup>,<italic>d</italic>) is consistent with recombination.</p></sec><sec><title>3.3 Composite likelihood methods are susceptible to rate heterogeneity</title><p>The composite likelihood method, which finds the composite likelihood surface of recombination rate, inferred non-zero recombination rates for all simulated datasets (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S4</ext-link>), revealing some degree of susceptibility to both temporal sampling and rate heterogeneity. A window-based approach of this test shows a sharp increase in the recombination rate estimated within 300 nucleotide windows around nucleotide 21000, close to the breakpoint inferred by GARD (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S5</ext-link>). We recovered a qualitatively similar pattern when analyzing empirically simulated sequences. It is important to note, however, that none of the simulated data, even under extreme heterogeneity, reproduced the same scale of the estimated recombination rate. Although in MERS-CoV data, the majority of 300 nucleotide windows after position 23000 have a recombination rate per base consistently higher than 0.005, only data simulated under extreme rate heterogeneity approach values as high as that.</p><p>In addition to the apparently higher recombination rate in regions with higher rates, we expect rate heterogeneity to produce a higher density of polymorphic sites in regions that are evolving faster. This is quite obvious in empirically simulated data with 3-fold rate heterogeneity—the region with higher rate also contains, on average, more polymorphic loci per window in the last third of the ‘genome’ than the first 20 kb (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S6</ext-link>). We only see hints of this in the actual MERS-CoV genome, with an apparent decline in polymorphism density from positions 5000 to 15000 which resembles that of the simulated data with 1.3-fold rate heterogeneity.</p></sec><sec><title>3.4 Homoplasies in MERS-CoV genomes are ubiquitous</title><p>Homoplasy analyses suggest that the MERS-CoV genome is rife with apparent homoplasies. Both maximum likelihood and Bayesian approaches to ancestral sequence reconstruction converge on similar patterns of homoplasy density (<xref ref-type="fig" rid="vev023-F3">Fig. 3</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S7</ext-link>). Both methods identify the region around the <italic>S</italic> (spike) gene as having a high density of synonymous homoplasies.</p><p>Empirically simulated sequences showed that homoplasies are not that unlikely in the absence of recombination. All sequences empirically simulated in <italic>π</italic>BUSS under a nucleotide substitution model had 2-fold homoplasies ranging in frequency from 0.0222 to 0.0550 of all polymorphic sites, with sequences simulated under higher levels of rate heterogeneity having more homoplasies and higher homoplasy degrees (<xref ref-type="fig" rid="vev023-F4">Fig. 4</xref>). However, even under a caricature model of rate heterogeneity, we did not reach the same degree of homoplasy as that observed in MERS-CoV, where homoplasic sites comprise as much as 0.1447 of all polymorphic sites and reach homoplasy degrees as high as 4. We were able to recover similar proportions of homoplasic sites from empirical simulations under a codon model (ranging from 0.0821 to 0.1518 of all polymorphic sites) but only under unrealistic values of rate heterogeneity. Even then, the distributions of homoplasy degrees for all simulated datasets are heavily skewed toward low homoplasy degrees, whereas the homoplasy degree distribution for MERS-CoV has a longer tail, indicating that homoplasies are more repeatable across the phylogeny.</p><p>Additional tests for excessive homoplasies, pairwise homoplasy index, and 3Seq also identify MERS-CoV sequences as being recombinant, albeit both spuriously identify some of the simulated datasets as recombinant (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S10</ext-link>). Similarly to LDhat results, however, there were no cases where both methods falsely inferred the presence of recombination for the same dataset.</p></sec><sec><title>3.5 Model testing supports a model including rate heterogeneity, but not alternative tree topologies</title><p>A model including rate heterogeneity alone across breakpoints inferred by the GARD method (i.e. linked trees, unlinked relaxed clocks) performs best when applied to MERS-CoV data (<xref ref-type="fig" rid="vev023-F5">Fig. 5</xref>, log marginal likelihoods: −48,137.86 and −48,138.91, using path and stepping stone sampling, respectively). The next best-performing model (log Bayes factor <inline-formula><mml:math id="MM6"><mml:mo>≈</mml:mo></mml:math></inline-formula>18) is linked trees and relaxed clocks. Overall, unlinking molecular clock rates between the two genomic partitions appears sufficient to dramatically improve model fit. Additionally, relaxed molecular clocks are preferred over strict molecular clocks (log BF >15, <xref ref-type="fig" rid="vev023-F5">Fig. 5</xref>).</p></sec></sec><sec><title>4 Discussion</title><sec><title>4.1 Recombination tests consistently point to recombination in MERS-CoV</title><p>The majority of methods we used (with the exception of marginal likelihood model testing) point to the combined effects of recombination and rate heterogeneity in the genome of MERS coronavirus. GARD (<xref ref-type="fig" rid="vev023-F1">Fig. 1</xref>) identified two breakpoints in the genome with high support. It is important to note that calculating the likelihood of a phylogenetic tree involves the estimation of multiple parameters, which at the very least must include a particular tree topology and a set of branch lengths. As such, we considered that the inference of a breakpoint in the MERS-CoV genome could be caused by systematic rate heterogeneity, which we address by evaluating rate heterogeneity across the first breakpoint in BEAST. We estimate an empirical rate heterogeneity ratio between MERS-CoV genome positions 23723–30126 and 1–23722 to be on the order between 1.3 and 1.5 (see <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S2</ext-link>). However, the support for this first breakpoint in MERS-CoV is comparable to support for empirically simulated sequences with 2-fold rate heterogeneity, and breakpoints under MERS-like levels of rate heterogeneity are difficult to detect. Overall, this suggests that evidence for differences in likelihoods between MERS Fragments 1 and 2 are beyond what would be expected from rate heterogeneity alone.
<fig id="vev023-F1" orientation="portrait" position="float"><label>Figure 1.</label><caption><p>Summary of GARD results. Colored boxes indicate fragments resulting from GARD-inferred breakpoints with corrected ΔAIC values shown on the right. Dashed line indicates the actual position where the evolutionary model for simulated sequences under three levels of rate heterogeneity is changed. Arrows at the top indicate the positions and names of coding sequences within the MERS-CoV genome.</p></caption><graphic xlink:href="vev023f1p"></graphic></fig><fig id="vev023-F2" orientation="portrait" position="float"><label>Figure 2.</label><caption><p>Summary of non-parametric tests for recombination. The percentile of the observed value for four statistics of LD decay (<italic>y</italic> axis) in the distribution of permuted datasets is indicated by color. Sequence datasets are shown on the <italic>x</italic> axis, starting with MERS-CoV sequences, followed by ten fastsimcoal2-simulated datasets and twelve empirically simulated datasets with different degrees of rate heterogeneity. Expected values for ideal datasets are shown in the last two columns, an ideal positive corresponds to the presence of recombination. Values falling between the 2.5th and 97.5th percentile are shown in green, values falling below the 2.5th percentile are in blue, those that are above the 97.5th percentile in red.</p></caption><graphic xlink:href="vev023f2p"></graphic></fig><fig id="vev023-F3" orientation="portrait" position="float"><label>Figure 3.</label><caption><p>Distribution of apparent homoplasies. Position along the genome is shown on the <italic>x</italic> axis and homoplasy degree, the number of times a particular mutation has occurred in excess in the tree as inferred by maximum likelihood, is shown on the <italic>y</italic> axis (left). Individual mutations are marked by vertical lines, synonymous ones in green and non-synonymous in red. The ratio of apparent homoplasy over synapomorphy kernel density estimates (bandwidth = 0.1) is shown in blue for synonymous (dashed) and non-synonymous (solid) sites separately. Arrows at the top indicate the positions and names of coding sequences within the MERS-CoV genome.</p></caption><graphic xlink:href="vev023f3p"></graphic></fig><fig id="vev023-F4" orientation="portrait" position="float"><label>Figure 4.</label><caption><p>Homoplasy prevalence in MERS-CoV and simulated datasets. Bars show the proportion of all polymorphic sites that are homoplasic, split by homoplasy degree as inferred by maximum likelihood, in MERS-CoV and datasets simulated with different degrees of rate heterogeneity in the presence or absence of site-specific constraint in the form of a codon model. Homoplasy degree indicates how many times a given mutation has occurred in excess in the phylogenetic tree.</p></caption><graphic xlink:href="vev023f4p"></graphic></fig><fig id="vev023-F5" orientation="portrait" position="float"><label>Figure 5.</label><caption><p>Summary of model comparisons. Difference in marginal likelihoods (Bayes factor) estimated by path-sampling between the worst model (linked strict molecular clock, unlinked trees) and all others. Asterisks indicate the best-performing model (unlinked relaxed clocks, linked trees, run 2) for MERS-CoV data. Analyses employing a relaxed molecular clock were run independently 3 times, those with a strict molecular clock 2 times. Marginal likelihoods estimated using stepping stone sampling gave identical results.</p></caption><graphic xlink:href="vev023f5p"></graphic></fig></p><p>Permutation tests show that statistics related to LD decay derived from MERS-CoV sequence data are outliers compared to permuted data (<xref ref-type="fig" rid="vev023-F2">Fig. 2</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S3</ext-link>). Sequences simulated empirically with varying levels of rate heterogeneity, under nucleotide or codon models of substitution, also have a tendency to exhibit extreme values for these statistics. However, only MERS-CoV data have values for all four tests that are in the direction consistent with recombination. For simulated datasets, especially those simulated under extreme levels of rate heterogeneity, values deviating significantly from the permuted data were recovered but often indicated internally contradictory scenarios.</p><p>Homoplasy analyses regardless of inference method show that MERS-CoV sequences contain a large number of homoplasic sites with high homoplasy degrees (<xref ref-type="fig" rid="vev023-F3">Fig. 3</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S7</ext-link>). Through sequence simulation, we also confirmed that both the numbers of homoplasic sites and their homoplasy degrees in MERS-CoV genomes are excessive, even when compared with unrealistic scenarios (e.g. 3-fold rate heterogeneity, <xref ref-type="fig" rid="vev023-F4">Fig. 4</xref>). Homoplasies become increasingly more prevalent when a more realistic codon model is used, due to differences in codon position constraint. Even then, MERS-CoV genomes possess mutations with much higher homoplasy degrees, surpassing simulated datasets with caricature levels of rate heterogeneity. This makes sense under a recombination scenario, as alleles persist in a diverse population and get recombined into novel backgrounds repeatedly, giving an appearance of highly repeatable mutations. Nevertheless, substitution patterns in real genomes are often highly complex and homoplasy-based methods have been shown to be susceptible to rate heterogeneity across sites, especially under higher levels of sequence divergence (<xref rid="vev023-B62" ref-type="bibr">Posada and Crandall 2001</xref>). Although rate heterogeneity certainly exists in MERS-CoV data, the divergence levels are still quite low (<italic>θ</italic>/site = 0.0047), giving us some degree of certainty in our inference of homoplasies. It is also reassuring that both maximum likelihood and Bayesian sequence reconstruction converged on similar patterns of homoplasy and synapomorphy across the genome (<xref ref-type="fig" rid="vev023-F3">Fig. 3</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S7</ext-link>). This is important, since homoplasies inferred using BEAST are integrated over all possible tree topologies, whereas homoplasies inferred by maximum likelihood were conditioned on a single tree. The convergence between these two methods suggests that the data contain enough phylogenetic signal to recover what could be called a ‘true’ tree and that homoplasies, for the most part, can be correctly identified as such.</p><p>One major concern surrounding the inference of homoplasies is host adaptation. There are a number of canonical mutations associated with host shifts, e.g. the glutamic acid to lysine amino acid substitution at position 627 in the PB2 protein of avian influenza A viruses confers the ability of the virus to replicate in mammals (<xref rid="vev023-B73" ref-type="bibr">Subbarao, London, and Murphy 1993</xref>) and a small number of amino acid substitutions in Parvoviruses are associated with adaptations to different hosts (<xref rid="vev023-B13" ref-type="bibr">Chang, Sgro, and Parrish 1992</xref>). If MERS-CoV is repeatedly emerging in humans, convergent mutations would be expected to arise that might allow the virus to adapt to humans.</p><p>However, we expect most host-adaptation mutations to be non-synonymous, whereas we detect both non-synonymous and synonymous homoplasies. This implies the action of recombination, rather than repeated selection for the same host-specific mutations. Furthermore, we do not detect any strong associations between host and particular alleles (<ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S8</ext-link>), although we do not believe that there is a sufficient number of sequences from camels to have much confidence in this result.</p><p>The overall phylogenetic and genomic patterns of homoplasies are consistent with fairly frequent recombination through time (<xref ref-type="fig" rid="vev023-F5">Fig. 5</xref>). Recent recombination should result in long homoplasy tracts shared across branches in the phylogenetic tree. At most we observe two stretches of adjacent homoplasies, one encompassing three homoplasies and another encompassing two homoplasies that are shared between taxa and likely to be caused by recent recombination. The majority of homoplasies that we observe, however, occur on their own. Recombination tracts, rather than single template switches, are not uncommon in other coronaviruses (<xref rid="vev023-B40" ref-type="bibr">Keck et al. 1988</xref>; <xref rid="vev023-B43" ref-type="bibr">Kottier, Cavanagh, and Britton 1995</xref>; <xref rid="vev023-B36" ref-type="bibr">Herrewegh et al. 1998</xref>). Thus, in MERS-CoV, we interpret extremely short homoplasy tracts as evidence of relatively frequent recombination. Alternatively, recombination tracts might be short and thus unable to transfer multiple informative sites across lineages.</p><p>Unlike all other tests we performed, model testing through marginal likelihoods indicates that models including rate heterogeneity explain MERS-CoV data partitioned across a well-supported breakpoint better than models including independent trees. At first, this may seem paradoxical, but we believe this result is due to the combined effects of the way homoplasious sites are distributed across the genome and phylogenetic tree of MERS-CoV (<xref ref-type="fig" rid="vev023-F5">Fig. 5</xref>) and the number of parameters involved. A speckled pattern of homoplasious sites without phylogenetic signal could easily be overwhelmed by the signal coming from the sites that support what could be called ‘the one true tree’, i.e. the clonal frame, in the data. Second, each phylogenetic tree contains at least <italic>n</italic> − 1 free parameters, so it is not surprising then that models attempting to recover two independent trees for both genomic fragments resulting from alternative tree topology analysis of MERS-CoV with highly correlated genealogies are penalized for the extra parameters introduced by a second tree. We would additionally like to point out that the fit of models including relaxed molecular clocks result in dramatic improvements to model fit compared to models with strict molecular clocks (log BF > 15, <xref ref-type="fig" rid="vev023-F6">Fig. 6</xref>). Although this could be interpreted as evidence for a considerable degree of lineage rate heterogeneity, the more parsimonious explanation is the ability of a relaxed molecular clock to accommodate homoplasies of recombinant origin, which do not necessarily accumulate at a relatively constant rate like genuine <italic>de novo</italic> mutations do.
<fig id="vev023-F6" orientation="portrait" position="float"><label>Figure 6.</label><caption><p>Mutations mapped onto a ML phylogeny. A maximum likelihood phylogeny of 85 MERS-CoV sequences with maximum likelihood-mapped mutations. Synapomorphies are shown as coloured ticks (coloured by coding sequence in which they occur) on branches where they occur. Homoplasies are shown as circles connected with coloured lines, colour corresponds with the coding sequence in which the mutation has occured. Mutations are positioned on the branches in proportion to where the mutation occurs in the genome, e.g. mutations shown towards the end of a branch correspond to mutations near the 3' terminus of the genome. Arrows at the top indicate the order, relative length and names of coding sequences within the MERS-CoV genome.</p></caption><graphic xlink:href="vev023f6p"></graphic></fig></p></sec><sec><title>4.2 Implications for future analyses</title><p>Recombination aside, MERS-CoV genomes exhibit a significant degree of rate heterogeneity among sites. Marginal likelihood analyses indicate that estimating independent molecular clocks after partitioning the MERS-CoV genome into two fragments alone substantially increases model fit over a completely linked (trees and clocks) model (log Bayes factor <inline-formula><mml:math id="MM7"><mml:mo>≈</mml:mo></mml:math></inline-formula>18). This highlights the advantage of employing relaxed molecular clocks, as in our case the method is clearly capable of accommodating recombination in an otherwise entirely clonal analysis framework. In addition, previous studies of SARS-like coronaviruses in bats have identified recombination breakpoints in small numbers of isolates falling close to the ‘transition zone’ around site 22000 (<xref rid="vev023-B19" ref-type="bibr">Hon et al. 2008</xref>; <xref rid="vev023-B45" ref-type="bibr">Lau et al. 2010</xref>) which in our analysis of MERS-CoV is where GARD, LDhat, and BEAST identify changes in the underlying model of evolution (<xref ref-type="fig" rid="vev023-F1">Figs. 1</xref> and <xref ref-type="fig" rid="vev023-F3">3</xref> and <ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary Fig. S5</ext-link>). Overall, a more detailed investigation will need to be done to determine if empirical patterns of rate variation in MERS-CoV have the potential to generate apparent recombination signals.</p><p>We also expect that as more sequences of MERS-CoV become available more homoplasies will be detected, some contributing to the homoplasy degree of the homoplasies already reported here, some previously unknown and some turning mutations currently thought of as synapomorphies into homoplasies. Although new sequences are likely to come from human cases, we think that sequencing MERS-CoV circulating in dromedary camels is of extreme importance from both surveillance and epidemiological points of view.</p></sec><sec><title>4.3 Implications about the virus population structure and infection dynamics</title><p>Our results point toward frequent recombination in MERS-CoV in the recent history of the MERS-CoV outbreak. For this to occur, different lineages of the virus must encounter each other often and implies frequent co-infection with MERS-CoV. To date, it is difficult to ascertain whether the human infections with MERS-CoV are a result of substantial asymptomatic transmission among humans or repeated zoonosis of the virus from camels to humans or a combination thereof. Given the severity of MERS, we find it unlikely that humans could be sufficiently frequently co-infected with two or more different lineages of the virus. Previous serological studies have failed to find evidence of prevalent past MERS-CoV infections of humans (<xref rid="vev023-B1" ref-type="bibr">Aburizaiza et al. 2013</xref>; <xref rid="vev023-B28" ref-type="bibr">Gierer et al. 2013</xref>), although a recent nation-wide study in Saudi Arabia has detected non-negligible numbers of individuals with antibodies against MERS-CoV, especially among shepherds and slaughterhouse workers (<xref rid="vev023-B59" ref-type="bibr">Müller et al. 2015</xref>). We thus propose that MERS-CoV mostly infects, and recombines, in camels. A study by <xref rid="vev023-B2" ref-type="bibr">Adney et al. (2014)</xref> has shown that camels only suffer mild symptoms from MERS-CoV infection and numerous other studies indicate an extremely high prevalence of antibodies specific against MERS-CoV in camels (<xref rid="vev023-B15" ref-type="bibr">Chu et al. 2014</xref>; <xref rid="vev023-B17" ref-type="bibr">Corman et al. 2014b</xref>; <xref rid="vev023-B58" ref-type="bibr">Müller et al. 2014</xref>; <xref rid="vev023-B66" ref-type="bibr">Reusken et al. 2013</xref>, <xref rid="vev023-B67" ref-type="bibr">2014</xref>). At the same time, however, sequencing has not indicated the presence of multiple infection in camels or any other animal. We believe that individual MERS-CoV co-infections are rare, but given the size of the epidemic in camels, as inferred from serology, the total number of co-infections is high. In addition, MERS-CoV infection is transient in camels (<xref rid="vev023-B2" ref-type="bibr">Adney et al. 2014</xref>), and thus sequencing efforts, which have been insufficient and very limited in camels, are highly unlikely to capture a co-infection.</p><p>Another point worth considering is that alleles that have arisen through mutation in MERS-CoV can be recombined, increasing the genetic variation of the virus (<xref rid="vev023-B57" ref-type="bibr">Muller 1932</xref>). Whether this is of epidemiological importance for humans depends entirely on what alleles are circulating in the reservoir, although there is no evidence that MERS-CoV is particularly likely to become as transmissible as common human pathogens or even SARS-CoV.</p></sec></sec><sec><title>Data availability</title><p>Python scripts used to process trees and sequences are available at: <ext-link ext-link-type="uri" xlink:href="https://github.com/evogytis/MERS_recombination/tree/master/scripts">https://github.com/evogytis/MERS_recombination/tree/master/scripts</ext-link>. Input and output files for programs used are publicly available at: <ext-link ext-link-type="uri" xlink:href="https://github.com/evogytis/MERS_recombination">https://github.com/evogytis/MERS_recombination</ext-link>.</p></sec><sec sec-type="supplementary-material"><title>Supplementary Material</title><supplementary-material id="PMC_1" content-type="local-data"><caption><title>Supplementary Fig. S1</title></caption><media mimetype="text" mime-subtype="html" xlink:href="ve_vev023_index.html"></media><media xlink:role="associated-file" mimetype="application" mime-subtype="pdf" xlink:href="MERS_rec_supp.pdf"></media></supplementary-material></sec></body><back><ack><title>Acknowledgements</title><p>G.D. was supported by a <funding-source>Natural Environment Research Council</funding-source> studentship <award-id>D76739X</award-id>. The research leading to these results has received funding from the <funding-source>European Research Council under the European Community’s Seventh Framework Programme</funding-source> (<award-id>FP7/2007-2013</award-id>) under Grant Agreement no. <award-id>278433</award-id>-<funding-source>PREDEMICS and ERC</funding-source> Grant agreement no. <award-id>260864</award-id>.</p><sec><title>Supplementary data</title><p><ext-link ext-link-type="uri" xlink:href="http://ve.oxfordjournals.org/lookup/suppl/doi:10.1093/ve/vev023/-/DC1">Supplementary data</ext-link> are available at <italic>Virus Evolution</italic> online.</p><p><bold>Conflict of interest</bold>: None declared.</p></sec></ack><ref-list><title>References</title><ref id="vev023-B1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aburizaiza</surname><given-names>A. S.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Investigation of Anti-MERS-Coronavirus Antibodies in Blood Donors and Abbatoir Workers in Jeddah and Makkah, Kingdom of Saudi Arabia, Fall 2012’</article-title>, <source>Journal of Infectious Diseases</source>, <volume>209</volume>(<issue>2</issue>): <fpage>243</fpage>–<lpage>6</lpage>. <comment>doi: 10.1093/infdis/jit589.</comment><pub-id pub-id-type="pmid">24218504</pub-id></mixed-citation></ref><ref id="vev023-B2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Adney</surname><given-names>D. R.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Replication and Shedding of MERS-CoV in Upper Respiratory Tract of Inoculated Dromedary Camels’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>20</volume>: <fpage>1999</fpage>–<lpage>2005</lpage>.<pub-id pub-id-type="pmid">25418529</pub-id></mixed-citation></ref><ref id="vev023-B3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anisimova</surname><given-names>M.</given-names></name><name><surname>Nielsen</surname><given-names>R.</given-names></name><name><surname>Yang</surname><given-names>Z.</given-names></name></person-group> (<year>2003</year>) <article-title>‘Effect of Recombination on the Accuracy of the Likelihood Method for Detecting Positive Selection at Amino Acid Sites’</article-title>, <source>Genetics</source>, <volume>164</volume>: <fpage>1229</fpage>–<lpage>36</lpage>.<pub-id pub-id-type="pmid">12871927</pub-id></mixed-citation></ref><ref id="vev023-B4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Assiri</surname><given-names>A.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Epidemiological, Demographic, and Clinical Characteristics of 47 Cases of Middle East Respiratory Syndrome Coronavirus Disease from Saudi Arabia: A Descriptive Study’</article-title>, <source>The Lancet Infectious Diseases</source>, <volume>13</volume>: <fpage>752</fpage>–<lpage>61</lpage>.<pub-id pub-id-type="pmid">23891402</pub-id></mixed-citation></ref><ref id="vev023-B5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baele</surname><given-names>G.</given-names></name><etal></etal></person-group> (<year>2012</year>) <article-title>‘Improving the Accuracy of Demographic and Molecular Clock Model Comparison While Accommodating Phylogenetic Uncertainty’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>29</volume>: <fpage>2157</fpage>–<lpage>67</lpage>.<pub-id pub-id-type="pmid">22403239</pub-id></mixed-citation></ref><ref id="vev023-B6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baric</surname><given-names>R. S.</given-names></name><etal></etal></person-group> (<year>1987</year>) <article-title>‘Analysis of Intracellular Small RNAs of Mouse Hepatitis Virus: Evidence for Discontinuous Transcription’</article-title>, <source>Virology</source>, <volume>156</volume>: <fpage>342</fpage>–<lpage>54</lpage>.<pub-id pub-id-type="pmid">3027983</pub-id></mixed-citation></ref><ref id="vev023-B7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bielejec</surname><given-names>F.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘<italic>π</italic>BUSS: A Parallel BEAST/BEAGLE Utility for Sequence Simulation under Complex Evolutionary Scenarios’</article-title>, <source>BMC Bioinformatics</source>, <volume>15</volume>: <fpage>133</fpage>.<pub-id pub-id-type="pmid">24885610</pub-id></mixed-citation></ref><ref id="vev023-B8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boni</surname><given-names>M. F.</given-names></name><name><surname>Posada</surname><given-names>D.</given-names></name><name><surname>Feldman</surname><given-names>M. W.</given-names></name></person-group> (<year>2007</year>) <article-title>‘An Exact Nonparametric Method for Inferring Mosaic Structure in Sequence Triplets’</article-title>, <source>Genetics</source>, <volume>176</volume>: <fpage>1035</fpage>–<lpage>47</lpage>.<pub-id pub-id-type="pmid">17409078</pub-id></mixed-citation></ref><ref id="vev023-B9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boni</surname><given-names>M. F.</given-names></name><etal></etal></person-group> (<year>2010</year>) <article-title>‘Guidelines for Identifying Homologous Recombination Events in Influenza A Virus’</article-title>, <source>PLoS One</source>, <volume>5</volume>: <fpage>e10434</fpage>.<pub-id pub-id-type="pmid">20454662</pub-id></mixed-citation></ref><ref id="vev023-B10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boucher</surname><given-names>C. A.</given-names></name><etal></etal></person-group> (<year>1993</year>) <article-title>‘High-Level Resistance to (-) Enantiomeric 2’-Deoxy-3’-Thiacytidine in vitro is due to One Amino Acid Substitution in the Catalytic Site of Human Immunodeficiency Virus Type 1 reverse Transcriptase’</article-title>, <source>Antimicrobial Agents and Chemotherapy</source>, <volume>37</volume>: <fpage>2231</fpage>–<lpage>4</lpage>.<pub-id pub-id-type="pmid">7504909</pub-id></mixed-citation></ref><ref id="vev023-B11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bruen</surname><given-names>T. C.</given-names></name><name><surname>Philippe</surname><given-names>H.</given-names></name><name><surname>Bryant</surname><given-names>D.</given-names></name></person-group> (<year>2006</year>) <article-title>‘A Simple and Robust Statistical Test for Detecting the Presence of Recombination’</article-title>, <source>Genetics</source>, <volume>172</volume>: <fpage>2665</fpage>–<lpage>81</lpage>.<pub-id pub-id-type="pmid">16489234</pub-id></mixed-citation></ref><ref id="vev023-B12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cauchemez</surname><given-names>S.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Middle East Respiratory Syndrome Coronavirus: Quantification of the Extent of the Epidemic, Surveillance Biases, and Transmissibility’</article-title>, <source>The Lancet Infectious Diseases</source>, <volume>14</volume>: <fpage>50</fpage>–<lpage>6</lpage>.<pub-id pub-id-type="pmid">24239323</pub-id></mixed-citation></ref><ref id="vev023-B13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname><given-names>S. F.</given-names></name><name><surname>Sgro</surname><given-names>J. Y.</given-names></name><name><surname>Parrish</surname><given-names>C. R.</given-names></name></person-group> (<year>1992</year>) <article-title>‘Multiple Amino Acids in the Capsid Structure of Canine Parvovirus Coordinately Determine the Canine Host Range and Specific Antigenic and Hemagglutination Properties’</article-title>, <source>Journal of Virology</source>, <volume>66</volume>: <fpage>6858</fpage>–<lpage>67</lpage>.<pub-id pub-id-type="pmid">1331498</pub-id></mixed-citation></ref><ref id="vev023-B14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chare</surname><given-names>E. R.</given-names></name><name><surname>Gould</surname><given-names>E. A.</given-names></name><name><surname>Holmes</surname><given-names>E. C.</given-names></name></person-group> (<year>2003</year>) <article-title>‘Phylogenetic Analysis Reveals a Low Rate of Homologous Recombination in Negative-Sense RNA Viruses’</article-title>, <source>Journal of General Virology</source>, <volume>84</volume>: <fpage>2691</fpage>–<lpage>703</lpage>.<pub-id pub-id-type="pmid">13679603</pub-id></mixed-citation></ref><ref id="vev023-B15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chu</surname><given-names>D. K.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘MERS Coronaviruses in Dromedary Camels, Egypt’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>20</volume>: <fpage>1049</fpage>–<lpage>53</lpage>.<pub-id pub-id-type="pmid">24856660</pub-id></mixed-citation></ref><ref id="vev023-B16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Corman</surname><given-names>V. M.</given-names></name><etal></etal></person-group> (<year>2014a</year>) <article-title>‘Rooting the Phylogenetic Tree of Middle East Respiratory Syndrome Coronavirus by Characterization of a Conspecific Virus from an African Bat’</article-title>, <source>Journal of Virology</source>, <volume>88</volume>: <fpage>11297</fpage>–<lpage>303</lpage>.<pub-id pub-id-type="pmid">25031349</pub-id></mixed-citation></ref><ref id="vev023-B17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Corman</surname><given-names>V. M.</given-names></name><etal></etal></person-group> (<year>2014b</year>) <article-title>‘Antibodies against MERS Coronavirus in Dromedary Camels, Kenya, 1992–2013’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>20</volume>.</mixed-citation></ref><ref id="vev023-B18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cotten</surname><given-names>M.</given-names></name><etal></etal></person-group> (<year>2013</year>) <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>: <fpage>1993</fpage>–<lpage>2002</lpage>.</mixed-citation></ref><ref id="vev023-B19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cotten</surname><given-names>M.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Spread, Circulation, and Evolution of the Middle East Respiratory Syndrome Coronavirus’</article-title>, <source>MBio</source>, <volume>5</volume>: <fpage>e01062</fpage>–<lpage>13</lpage>.<pub-id pub-id-type="pmid">24549846</pub-id></mixed-citation></ref><ref id="vev023-B20"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Didelot</surname><given-names>X.</given-names></name><name><surname>Falush</surname><given-names>D.</given-names></name></person-group> (<year>2007</year>) <article-title>‘Inference of Bacterial Microevolution Using Multilocus Sequence Data’</article-title>, <source>Genetics</source>, <volume>175</volume>: <fpage>1251</fpage>–<lpage>66</lpage>.<pub-id pub-id-type="pmid">17151252</pub-id></mixed-citation></ref><ref id="vev023-B21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dong</surname><given-names>G.</given-names></name><etal></etal></person-group> (<year>2011</year>) <article-title>‘Phylogenetic Diversity and Genotypical Complexity of H9N2 Influenza a Viruses Revealed by Genomic Sequence Analysis’</article-title>, <source>PLoS One</source>, <volume>6</volume>: <fpage>e17212</fpage>.<pub-id pub-id-type="pmid">21386964</pub-id></mixed-citation></ref><ref id="vev023-B22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drosten</surname><given-names>C.</given-names></name><etal></etal></person-group> (<year>2015</year>) <article-title>‘An Observational, Laboratory-Based Study of Outbreaks of MERS-Coronavirus in Jeddah and Riyadh, Kingdom of Saudi Arabia, 2014’</article-title>, <source>Clinical Infectious Diseases</source>, <volume>60</volume>(<issue>3</issue>): <fpage>369</fpage>–<lpage>77</lpage>. <comment>doi: 10.1093/cid/ciu812.</comment><pub-id pub-id-type="pmid">25323704</pub-id></mixed-citation></ref><ref id="vev023-B23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drosten</surname><given-names>C.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Clinical Features and Virological Analysis of a Case of Middle East Respiratory Syndrome Coronavirus Infection’</article-title>, <source>The Lancet Infectious Diseases</source>, <volume>13</volume>: <fpage>745</fpage>–<lpage>51</lpage>.<pub-id pub-id-type="pmid">23782859</pub-id></mixed-citation></ref><ref id="vev023-B24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname><given-names>A. J.</given-names></name><etal></etal></person-group> (<year>2006</year>) <article-title>‘Relaxed Phylogenetics and Dating with Confidence’</article-title>, <source>PLoS Biology</source>, <volume>4</volume>: <fpage>e88</fpage>.<pub-id pub-id-type="pmid">16683862</pub-id></mixed-citation></ref><ref id="vev023-B25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname><given-names>A. J.</given-names></name><etal></etal></person-group> (<year>2012</year>) <article-title>‘Bayesian Phylogenetics with BEAUti and the BEAST 1.7’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>29</volume>(<issue>8</issue>): <fpage>1969</fpage>–<lpage>73</lpage>.<pub-id pub-id-type="pmid">22367748</pub-id></mixed-citation></ref><ref id="vev023-B26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Excoffier</surname><given-names>L.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Robust Demographic Inference from Genomic and SNP Data’</article-title>, <source>PLoS Genetics</source>, <volume>9</volume>: <fpage>e1003905</fpage>.<pub-id pub-id-type="pmid">24204310</pub-id></mixed-citation></ref><ref id="vev023-B27"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ferreira</surname><given-names>M. A. R.</given-names></name><name><surname>Suchard</surname><given-names>M. A.</given-names></name></person-group> (<year>2008</year>) <article-title>‘Bayesian Analysis of Elapsed Times in Continuous-Time Markov Chains’</article-title>, <source>Canadian Journal of Statistics</source>, <volume>36</volume>: <fpage>355</fpage>–<lpage>68</lpage>.</mixed-citation></ref><ref id="vev023-B28"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gierer</surname><given-names>S.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Lack of MERS Coronavirus Neutralizing Antibodies in Humans, Eastern Province, Saudi Arabia’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>19</volume>: <fpage>2034</fpage>–<lpage>6</lpage>.<pub-id pub-id-type="pmid">24274664</pub-id></mixed-citation></ref><ref id="vev023-B29"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gill</surname><given-names>M. S.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Improving Bayesian Population Dynamics Inference: A Coalescent-Based Model for Multiple Loci’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>30</volume>: <fpage>713</fpage>–<lpage>24</lpage>.<pub-id pub-id-type="pmid">23180580</pub-id></mixed-citation></ref><ref id="vev023-B30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gire</surname><given-names>S. K.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Genomic Surveillance Elucidates Ebola Virus Origin and Transmission during the 2014 Outbreak’</article-title>, <source>Science</source>, <volume>345</volume>: <fpage>1369</fpage>–<lpage>72</lpage>.<pub-id pub-id-type="pmid">25214632</pub-id></mixed-citation></ref><ref id="vev023-B31"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goldman</surname><given-names>N.</given-names></name><name><surname>Yang</surname><given-names>Z.</given-names></name></person-group> (<year>1994</year>) <article-title>‘A Codon-Based Model of Nucleotide Substitution for Protein-Coding DNA Sequences’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>11</volume>: <fpage>725</fpage>–<lpage>36</lpage>.<pub-id pub-id-type="pmid">7968486</pub-id></mixed-citation></ref><ref id="vev023-B32"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gubareva</surname><given-names>L. V.</given-names></name><etal></etal></person-group> (<year>2001</year>) <article-title>‘Selection of Influenza Virus Mutants in Experimentally Infected Volunteers Treated with Oseltamivir’</article-title>, <source>Journal of Infectious Diseases</source>, <volume>183</volume>: <fpage>523</fpage>–<lpage>31</lpage>.<pub-id pub-id-type="pmid">11170976</pub-id></mixed-citation></ref><ref id="vev023-B33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guindon</surname><given-names>S.</given-names></name><name><surname>Gascuel</surname><given-names>O.</given-names></name></person-group> (<year>2003</year>) <article-title>‘A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood’</article-title>, <source>Systematic Biology</source>, <volume>52</volume>: <fpage>696</fpage>–<lpage>704</lpage>.<pub-id pub-id-type="pmid">14530136</pub-id></mixed-citation></ref><ref id="vev023-B34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hasegawa</surname><given-names>M.</given-names></name><name><surname>Kishino</surname><given-names>H.</given-names></name><name><surname>Yano</surname><given-names>T. A.</given-names></name></person-group> (<year>1985</year>) <article-title>‘Dating of the Human-Ape Splitting by a Molecular Clock of Mitochondrial DNA’</article-title>, <source>Journal of Molecular Evolution</source>, <volume>22</volume>: <fpage>160</fpage>–<lpage>74</lpage>.<pub-id pub-id-type="pmid">3934395</pub-id></mixed-citation></ref><ref id="vev023-B35"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hedrick</surname><given-names>P. W.</given-names></name><name><surname>Thomson</surname><given-names>G.</given-names></name></person-group> (<year>1986</year>) <article-title>‘A Two-Locus Neutrality Test: Applications to Humans, <italic>E. coli</italic> and Lodgepole Pine’</article-title>, <source>Genetics</source>, <volume>112</volume>: <fpage>135</fpage>–<lpage>56</lpage>.<pub-id pub-id-type="pmid">3510942</pub-id></mixed-citation></ref><ref id="vev023-B36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Herrewegh</surname><given-names>A. A.</given-names></name><etal></etal></person-group> (<year>1998</year>) <article-title>‘Feline Coronavirus Type II Strains 79-1683 and 79-1146 Originate from a Double Recombination between Feline Coronavirus Type I and Canine Coronavirus’</article-title>, <source>Journal of Virology</source>, <volume>72</volume>: <fpage>4508</fpage>–<lpage>14</lpage>.<pub-id pub-id-type="pmid">9557750</pub-id></mixed-citation></ref><ref id="vev023-B37"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hill</surname><given-names>W. G.</given-names></name><name><surname>Robertson</surname><given-names>A.</given-names></name></person-group> (<year>1968</year>) <article-title>‘Linkage Disequilibrium in Finite Populations’</article-title>, <source>Theoretical and Applied Genetics. Theoretische und Angewandte Genetik</source>, <volume>38</volume>: <fpage>226</fpage>–<lpage>31</lpage>.<pub-id pub-id-type="pmid">24442307</pub-id></mixed-citation></ref><ref id="vev023-B38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Holmes</surname><given-names>E. C.</given-names></name><name><surname>Urwin</surname><given-names>R.</given-names></name><name><surname>Maiden</surname><given-names>M. C.</given-names></name></person-group> (<year>1999</year>) <article-title>‘The Influence of Recombination on the Population Structure and Evolution of the Human Pathogen <italic>Neisseria meningitidis</italic>’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>16</volume>: <fpage>741</fpage>–<lpage>9</lpage>.<pub-id pub-id-type="pmid">10368953</pub-id></mixed-citation></ref><ref id="vev023-B39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hon</surname><given-names>C. C.</given-names></name><etal></etal></person-group> (<year>2008</year>) <article-title>‘Evidence of the Recombinant Origin of a Bat Severe Acute Respiratory Syndrome (SARS)-Like Coronavirus and Its Implications on the Direct Ancestor of SARS Coronavirus’</article-title>, <source>Journal of Virology</source>, <volume>82</volume>: <fpage>1819</fpage>–<lpage>26</lpage>.<pub-id pub-id-type="pmid">18057240</pub-id></mixed-citation></ref><ref id="vev023-B40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Keck</surname><given-names>J. G.</given-names></name><etal></etal></person-group> (<year>1988</year>) <article-title>‘in vivo RNA-RNA Recombination of Coronavirus in Mouse Brain’</article-title>, <source>Journal of Virology</source>, <volume>62</volume>: <fpage>1810</fpage>–<lpage>13</lpage>.<pub-id pub-id-type="pmid">2833625</pub-id></mixed-citation></ref><ref id="vev023-B41"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kirkegaard</surname><given-names>K.</given-names></name><name><surname>Baltimore</surname><given-names>D.</given-names></name></person-group> (<year>1986</year>) <article-title>‘The Mechanism of RNA Recombination in Poliovirus’</article-title>, <source>Cell</source>, <volume>47</volume>: <fpage>433</fpage>–<lpage>43</lpage>.<pub-id pub-id-type="pmid">3021340</pub-id></mixed-citation></ref><ref id="vev023-B42"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kosakovsky Pond</surname><given-names>S. L.</given-names></name><etal></etal></person-group> (<year>2006</year>) <article-title>‘Automated Phylogenetic Detection of Recombination Using a Genetic Algorithm’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>23</volume>: <fpage>1891</fpage>–<lpage>901</lpage>.<pub-id pub-id-type="pmid">16818476</pub-id></mixed-citation></ref><ref id="vev023-B43"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kottier</surname><given-names>S. A.</given-names></name><name><surname>Cavanagh</surname><given-names>D.</given-names></name><name><surname>Britton</surname><given-names>P.</given-names></name></person-group> (<year>1995</year>) <article-title>‘Experimental Evidence of Recombination in Coronavirus Infectious Bronchitis Virus’</article-title>, <source>Virology</source>, <volume>213</volume>: <fpage>569</fpage>–<lpage>80</lpage>.<pub-id pub-id-type="pmid">7491781</pub-id></mixed-citation></ref><ref id="vev023-B44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lai</surname><given-names>M. M.</given-names></name><etal></etal></person-group> (<year>1985</year>) <article-title>‘Recombination between Nonsegmented RNA Genomes of Murine Coronaviruses’</article-title>, <source>Journal of Virology</source>, <volume>56</volume>: <fpage>449</fpage>–<lpage>56</lpage>.<pub-id pub-id-type="pmid">2997467</pub-id></mixed-citation></ref><ref id="vev023-B45"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lau</surname><given-names>S. K. P.</given-names></name><etal></etal></person-group> (<year>2010</year>) <article-title>‘Ecoepidemiology and Complete Genome Comparison of Different Strains of Severe Acute Respiratory Syndrome-Related Rhinolophus Bat Coronavirus in China Reveal Bats as a Reservoir for Acute, Self-Limiting Infection that Allows Recombination Events’</article-title>, <source>Journal of Virology</source>, <volume>84</volume>: <fpage>2808</fpage>–<lpage>19</lpage>.<pub-id pub-id-type="pmid">20071579</pub-id></mixed-citation></ref><ref id="vev023-B46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lemey</surname><given-names>P.</given-names></name><name><surname>Suchard</surname><given-names>M.</given-names></name><name><surname>Rambaut</surname><given-names>A.</given-names></name></person-group> (<year>2009</year>) <article-title>‘Reconstructing the Initial Global Spread of a Human Influenza Pandemic’</article-title>, <source>PLoS Currents</source>, <volume>1</volume><comment>doi: 10.1371/currents.RRN1031</comment>.</mixed-citation></ref><ref id="vev023-B47"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lewontin</surname><given-names>R. C.</given-names></name></person-group> (<year>1964</year>) <article-title>‘The Interaction of Selection and Linkage. I. General Considerations; Heterotic Models’</article-title>, <source>Genetics</source>, <volume>49</volume>: <fpage>49</fpage>–<lpage>67</lpage>.<pub-id pub-id-type="pmid">17248194</pub-id></mixed-citation></ref><ref id="vev023-B48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>K. S.</given-names></name><etal></etal></person-group> (<year>2004</year>) <article-title>‘Genesis of a Highly Pathogenic and Potentially Pandemic H5N1 Influenza Virus in Eastern Asia’</article-title>, <source>Nature</source>, <volume>430</volume>: <fpage>209</fpage>–<lpage>13</lpage>.<pub-id pub-id-type="pmid">15241415</pub-id></mixed-citation></ref><ref id="vev023-B49"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname><given-names>L.</given-names></name><name><surname>Lycett</surname><given-names>S. J.</given-names></name><name><surname>Brown</surname><given-names>A. J. L.</given-names></name></person-group> (<year>2014</year>) <article-title>‘Reassortment Patterns of Avian Influenza Virus Internal Segments among Different Subtypes’</article-title>, <source>BMC Evolutionary Biology</source>, <volume>14</volume>: <fpage>16</fpage>.<pub-id pub-id-type="pmid">24456010</pub-id></mixed-citation></ref><ref id="vev023-B50"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Makino</surname><given-names>S.</given-names></name><etal></etal></person-group> (<year>1986</year>) <article-title>‘High-Frequency RNA Recombination of Murine Coronaviruses’</article-title>, <source>Journal of Virology</source>, <volume>57</volume>: <fpage>729</fpage>–<lpage>37</lpage>.<pub-id pub-id-type="pmid">3005623</pub-id></mixed-citation></ref><ref id="vev023-B51"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Maynard Smith</surname><given-names>J.</given-names></name><name><surname>Smith</surname><given-names>N. H.</given-names></name></person-group> (<year>1998</year>) <article-title>‘Detecting Recombination from Gene Trees’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>15</volume>: <fpage>590</fpage>–<lpage>9</lpage>.<pub-id pub-id-type="pmid">9580989</pub-id></mixed-citation></ref><ref id="vev023-B52"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McVean</surname><given-names>G.</given-names></name><name><surname>Awadalla</surname><given-names>P.</given-names></name><name><surname>Fearnhead</surname><given-names>P.</given-names></name></person-group> (<year>2002</year>) <article-title>‘A Coalescent-Based Method for Detecting and Estimating Recombination from Gene Sequences’</article-title>, <source>Genetics</source>, <volume>160</volume>: <fpage>1231</fpage>–<lpage>41</lpage>.<pub-id pub-id-type="pmid">11901136</pub-id></mixed-citation></ref><ref id="vev023-B53"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Memish</surname><given-names>Z. A.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Family Cluster of Middle East Respiratory Syndrome Coronavirus Infections’</article-title>, <source>New England Journal of Medicine</source>, <volume>368</volume>: <fpage>2487</fpage>–<lpage>94</lpage>.<pub-id pub-id-type="pmid">23718156</pub-id></mixed-citation></ref><ref id="vev023-B54"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meunier</surname><given-names>J.</given-names></name><name><surname>Eyre-Walker</surname><given-names>A.</given-names></name></person-group> (<year>2001</year>) <article-title>‘The Correlation between Linkage Disequilibrium and Distance: Implications for Recombination in Hominid Mitochondria’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>18</volume>: <fpage>2132</fpage>–<lpage>5</lpage>.<pub-id pub-id-type="pmid">11606711</pub-id></mixed-citation></ref><ref id="vev023-B55"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Milkman</surname><given-names>R.</given-names></name><name><surname>Bridges</surname><given-names>M. M.</given-names></name></person-group> (<year>1990</year>) <article-title>‘Molecular Evolution of the <italic>Escherichia coli</italic> Chromosome. III. Clonal Frames’</article-title>, <source>Genetics</source>, <volume>126</volume>: <fpage>505</fpage>–<lpage>17</lpage>.<pub-id pub-id-type="pmid">1979037</pub-id></mixed-citation></ref><ref id="vev023-B56"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mitnaul</surname><given-names>L. J.</given-names></name><etal></etal></person-group> (<year>2000</year>) <article-title>‘Balanced Hemagglutinin and Neuraminidase Activities are Critical for Efficient Replication of Influenza A Virus’</article-title>, <source>Journal of Virology</source>, <volume>74</volume>: <fpage>6015</fpage>–<lpage>20</lpage>.<pub-id pub-id-type="pmid">10846083</pub-id></mixed-citation></ref><ref id="vev023-B57"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Muller</surname><given-names>H. J.</given-names></name></person-group> (<year>1932</year>) <article-title>‘Some Genetic Aspects of Sex’</article-title>, <source>The American Naturalist</source>, <volume>66</volume>: <fpage>118</fpage>–<lpage>38</lpage>.</mixed-citation></ref><ref id="vev023-B58"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Müller</surname><given-names>M. A.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘MERS Coronavirus Neutralizing Antibodies in Camels, Eastern Africa, 1983–1997’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>20</volume>.</mixed-citation></ref><ref id="vev023-B59"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Müller</surname><given-names>M. A.</given-names></name><etal></etal></person-group> (<year>2015</year>) <article-title>‘Presence of Middle East Respiratory Syndrome Coronavirus Antibodies in Saudi Arabia: A Nationwide, Cross-Sectional, Serological Study’</article-title>, <source>The Lancet Infectious Diseases</source>, <volume>15</volume>: <fpage>559</fpage>–<lpage>64</lpage>.<pub-id pub-id-type="pmid">25863564</pub-id></mixed-citation></ref><ref id="vev023-B60"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>O’Brien</surname><given-names>J. D.</given-names></name><name><surname>Minin</surname><given-names>V. N.</given-names></name><name><surname>Suchard</surname><given-names>M. A.</given-names></name></person-group> (<year>2009</year>) <article-title>‘Learning to Count: Robust Estimates for Labeled Distances between Molecular Sequences’</article-title>, <source>Molecular Biology and Evolution</source>, <volume>26</volume>: <fpage>801</fpage>–<lpage>14</lpage>.<pub-id pub-id-type="pmid">19131426</pub-id></mixed-citation></ref><ref id="vev023-B61"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pond</surname><given-names>S. L. K.</given-names></name><name><surname>Frost</surname><given-names>S. D. W.</given-names></name><name><surname>Muse</surname><given-names>S. V.</given-names></name></person-group> (<year>2005</year>) <article-title>‘HyPhy: Hypothesis Testing Using Phylogenies’</article-title>, <source>Bioinformatics</source>, <volume>21</volume>: <fpage>676</fpage>–<lpage>9</lpage>.<pub-id pub-id-type="pmid">15509596</pub-id></mixed-citation></ref><ref id="vev023-B62"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Posada</surname><given-names>D.</given-names></name><name><surname>Crandall</surname><given-names>K. A.</given-names></name></person-group> (<year>2001</year>) <article-title>‘Evaluation of Methods for Detecting Recombination from DNA Sequences: Computer Simulations’</article-title>, <source>Proceedings of the National Academy of Sciences</source>, <volume>98</volume>: <fpage>13757</fpage>–<lpage>62</lpage>.</mixed-citation></ref><ref id="vev023-B63"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Posada</surname><given-names>D.</given-names></name></person-group>(<year>2002</year>) <article-title>‘The Effect of Recombination on the Accuracy of Phylogeny Estimation’</article-title>, <source>Journal of Molecular Evolution</source>, <volume>54</volume>: <fpage>396</fpage>–<lpage>402</lpage>.<pub-id pub-id-type="pmid">11847565</pub-id></mixed-citation></ref><ref id="vev023-B64"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Posada</surname><given-names>D.</given-names></name><name><surname>Crandall</surname><given-names>K. A.</given-names></name><name><surname>Crandall</surname><given-names>K. A.</given-names></name><name><surname>Holmes</surname><given-names>E. C.</given-names></name></person-group> (<year>2002</year>) <article-title>‘Recombination in Evolutionary Genomics’</article-title>, <source>Annual Review of Genetics</source>, <volume>36</volume>: <fpage>75</fpage>–<lpage>97</lpage>.</mixed-citation></ref><ref id="vev023-B65"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rambaut</surname><given-names>A.</given-names></name><name><surname>Holmes</surname><given-names>E.</given-names></name></person-group> (<year>2009</year>) <article-title>‘The Early Molecular Epidemiology of the Swine-Origin A/H1N1 Human Influenza Pandemic’</article-title>, <source>PLoS Currents</source>, <volume>1</volume><comment>doi: 10.1371/currents.RRN1003.</comment></mixed-citation></ref><ref id="vev023-B66"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reusken</surname><given-names>C. B.</given-names></name><etal></etal></person-group> (<year>2013</year>) <article-title>‘Middle East Respiratory Syndrome Coronavirus Neutralising Serum Antibodies in Dromedary Camels: A Comparative Serological Study’</article-title>, <source>The Lancet Infectious Diseases</source>, <volume>13</volume>: <fpage>859</fpage>–<lpage>66</lpage>.<pub-id pub-id-type="pmid">23933067</pub-id></mixed-citation></ref><ref id="vev023-B67"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reusken</surname><given-names>C. B.</given-names></name><etal></etal></person-group> (<year>2014</year>) <article-title>‘Geographic Distribution of MERS Coronavirus among Dromedary Camels, Africa’</article-title>, <source>Emerging Infectious Diseases</source>, <volume>20</volume>: <fpage>1370</fpage>–<lpage>4</lpage>.<pub-id pub-id-type="pmid">25062254</pub-id></mixed-citation></ref><ref id="vev023-B68"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>D. L.</given-names></name><name><surname>Hahn</surname><given-names>B. H.</given-names></name><name><surname>Sharp</surname><given-names>P. M.</given-names></name></person-group> (<year>1995</year>) <article-title>‘Recombination in AIDS Viruses’</article-title>, <source>Journal of Molecular Evolution</source>, <volume>40</volume>: <fpage>249</fpage>–<lpage>59</lpage>.<pub-id pub-id-type="pmid">7723052</pub-id></mixed-citation></ref><ref id="vev023-B69"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Robertson</surname><given-names>D. L.</given-names></name><etal></etal></person-group> (<year>1995</year>) <article-title>‘Recombination in HIV-1’</article-title>, <source>Nature</source>, <volume>374</volume>: <fpage>124</fpage>–<lpage>6</lpage>.<pub-id pub-id-type="pmid">7877682</pub-id></mixed-citation></ref><ref id="vev023-B70"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schierup</surname><given-names>M. H.</given-names></name><name><surname>Hein</surname><given-names>J.</given-names></name></person-group> (<year>2000</year>) <article-title>‘Consequences of Recombination on Traditional Phylogenetic Analysis’</article-title>, <source>Genetics</source>, <volume>156</volume>: <fpage>879</fpage>–<lpage>91</lpage>.<pub-id pub-id-type="pmid">11014833</pub-id></mixed-citation></ref><ref id="vev023-B71"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shriner</surname><given-names>D.</given-names></name><etal></etal></person-group> (<year>2003</year>) <article-title>‘Potential Impact of Recombination on Sitewise Approaches for Detecting Positive Natural Selection’</article-title>, <source>Genetics Research</source>, <volume>81</volume>: <fpage>115</fpage>–<lpage>21</lpage>.</mixed-citation></ref><ref id="vev023-B72"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>G. J. D.</given-names></name><etal></etal></person-group> (<year>2009</year>) <article-title>‘Origins and Evolutionary Genomics of the 2009 Swine-Origin H1N1 Influenza A Epidemic’</article-title>, <source>Nature</source>, <volume>459</volume>: <fpage>1122</fpage>–<lpage>5</lpage>.<pub-id pub-id-type="pmid">19516283</pub-id></mixed-citation></ref><ref id="vev023-B73"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Subbarao</surname><given-names>E. K.</given-names></name><name><surname>London</surname><given-names>W.</given-names></name><name><surname>Murphy</surname><given-names>B. R.</given-names></name></person-group> (<year>1993</year>) <article-title>‘A Single Amino Acid in the PB2 Gene of Influenza A Virus is a Determinant of Host Range’</article-title>, <source>Journal of Virology</source>, <volume>67</volume>: <fpage>1761</fpage>–<lpage>4</lpage>.<pub-id pub-id-type="pmid">8445709</pub-id></mixed-citation></ref><ref id="vev023-B74"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Tavaré</surname><given-names>S.</given-names></name></person-group> (<year>1986</year>) <source><italic>Some Mathematical Questions in Biology: DNA Sequence Analysis. Lectures on Mathematics in the Life Sciences</italic>, Volume 17</source>, pp. <fpage>57</fpage>–<lpage>86</lpage>. <publisher-name>American Mathematical Society</publisher-name>: <publisher-loc>Providence, RI, USA</publisher-loc>.</mixed-citation></ref><ref id="vev023-B75"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tisdale</surname><given-names>M.</given-names></name><etal></etal></person-group> (<year>1993</year>) <article-title>‘Rapid in vitro Selection of Human Immunodeficiency Virus Type 1 Resistant to 3’-Thiacytidine Inhibitors due to a Mutation in the YMDD Region of Reverse Transcriptase’</article-title>, <source>Proceedings of the National Academy of Sciences of the United States of America</source>, <volume>90</volume>: <fpage>5653</fpage>–<lpage>6</lpage>.<pub-id pub-id-type="pmid">7685907</pub-id></mixed-citation></ref><ref id="vev023-B76"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname><given-names>Z.</given-names></name></person-group> (<year>1994</year>) <article-title>‘Maximum Likelihood Phylogenetic Estimation from DNA Sequences with Variable Rates Over Sites: Approximate Methods’</article-title>, <source>Journal of Molecular Evolution</source>, <volume>39</volume>: <fpage>306</fpage>–<lpage>14</lpage>.<pub-id pub-id-type="pmid">7932792</pub-id></mixed-citation></ref><ref id="vev023-B77"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zaki</surname><given-names>A. M.</given-names></name><etal></etal></person-group> (<year>2012</year>) <article-title>‘Isolation of a Novel Coronavirus from a Man with Pneumonia in Saudi Arabia’</article-title>, <source>New England Journal of Medicine</source>, <volume>367</volume>: <fpage>1814</fpage>–<lpage>20</lpage>.<pub-id pub-id-type="pmid">23075143</pub-id></mixed-citation></ref></ref-list></back></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 001187 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 001187 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Sante
|area= MersV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:4989901
|texte= MERS-CoV recombination: implications about the reservoir and potential for adaptation
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:27774293" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a MersV1
| This area was generated with Dilib version V0.6.33. |  |